Next Article in Journal
Towards a Software-Defined Industrial IoT-Edge Network for Next-Generation Offshore Wind Farms: State of the Art, Resilience, and Self-X Network and Service Management
Next Article in Special Issue
An Arrhenius-Based Simulation Tool for Predicting Aging of Lithium Manganese Dioxide Primary Batteries in Implantable Medical Devices
Previous Article in Journal
Analyzing the Influence of Energy Consumption and Economic Complexity on Carbon Emissions: Evidence from Malaysia
Previous Article in Special Issue
Intelligent Learning Method for Capacity Estimation of Lithium-Ion Batteries Based on Partial Charging Curves
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Enhanced Ageing Model for Solid-State Batteries

by
Paolo Scaltrito
1,
Amirmasoud Lanjan
2 and
Seshasai Srinivasan
3,*
1
Department of Mechanical and Aerospaae Engineering (DIMEAS), Politecnico di Torino, 10129 Turin, Italy
2
Department of Mechanical Engineering, McMaster University, Hamilton, ON L8S 4L8, Canada
3
W Booth School of Engineering Practice and Technology, McMaster University, Hamilton, ON L8S 4L8, Canada
*
Author to whom correspondence should be addressed.
Energies 2024, 17(12), 2896; https://doi.org/10.3390/en17122896
Submission received: 11 April 2024 / Revised: 5 June 2024 / Accepted: 6 June 2024 / Published: 13 June 2024

Abstract

:
The emphasis in the automotive industry towards sustainable mobility has led to a significant interest in hybrid-electric drive-trains with high energy density batteries. Addressing the needs of this strategy, the battery market is exploring new technologies to improve the safety and lifespan of electric vehicles. To this end, there is a focus on the all-solid-state battery (ASSB) technology for its cycle capabilities. Filling the current void in the literature pertaining to accurate ageing models for ASSBs, in the present work, we propose an enhanced version of the numerical ageing model, originally developed for liquid electrolyte based batteries, to forecast the development of the solid electrolyte interface layer that is the major cause of battery capacity fading. The model has been tested on prototype batteries and reveals an accuracy of 99%. The capacity fade in ASSBs has been investigated under different conditions and the enhanced ageing model has been validated using experimental data from these batteries. The findings suggest that there is potential for solid-state batteries to be commercialized, although significant work is needed to match the manufacturing level of lithium-ion batteries embedded with liquid electrolytes.

1. Introduction

Improved energy efficiency and advances in renewable energy are essential to achieve zero emissions. With advances in these areas, a significant reduction in CO2 emissions can be accomplished in the transport sector. This sector, as indicated by the 2021 IEA study, is currently contributing nearly 24% of the total CO2 emissions around the globe, and 72% of this comes from passenger cars. This demands a consistent effort towards the decarbonization of the transport sector and energy generation via renewable energy. In response to this, electric vehicles (EVs) are rapidly permeating the market, proposing themselves as an alternative to conventional internal combustion engine (ICE) vehicles. Despite their raising sales numbers, EVs in their current state are not a robust solution for a definite move towards sustainable mobility. That is due to the fact that in terms of life cycle assessment (LCA), commencing from the raw materials and the production to the end of life and recycling, EVs are not a definitive answer compared to ICE-driven vehicles. This assertion is primarily related to the CO2 level adduced by the battery manufacturing, which represents nearly 40% of the emissions emanating from the manufacturing of the EVs [1]. Further, current Li-ion batteries (LIBs) require high production costs and suffer from a limited lifespan, posing a serious environmental problem. Needless to say, extending the battery’s lifespan is of crucial importance for the electrification of the transport sector. To this end, one of the most promising battery designs is the all solid-state lithium-ion battery (ASSB). The replacement of the liquid electrolyte with a solid electrolyte provides a significant improvement in the safety, energy density, and cycle performance of the battery [2,3,4,5]. With this immense potential, from a computational modelling perspective, to design new and novel ASSBs, there is an urgent need to develop advanced ageing models that can be used to analyze the next-generation ASSBs. However, as in the case of LIBs, the precise ageing mechanism in ASSBs is still largely unknown. Since its introduction in 1991, the graphite anode, with a specific charge capacity of 372 mAh g−1 and a low operating potential of about 0.1 V versus Li+/Li, has become a staple in lithium-ion batteries (LIBs) when used with ethylene carbonate (EC)-based liquid carbonate electrolytes and LiPF6 salt [6]. Despite the inherent instability of these electrolytes at low potentials, their decomposition products form a protective layer known as the solid electrolyte interphase (SEI) on the anode surface, preventing further degradation and graphite exfoliation, thereby ensuring the long-term stability of LIBs [7]. Improvements in SEI through advanced electrolyte compositions, additives, and specific formation cycles have been essential in enhancing the performance of modern LIBs [8]. Current commercial LIBs achieve around 270 Wh kg−1 at the cell level [9]. To exceed this and reach specific energies above 400 Wh kg−1, researchers are exploring alternative anode materials [10]. One approach involves replacing graphite with silicon or lithium metal electrodes. Silicon, with a theoretical specific charge capacity of about 3600 mAh g−1, operates at a similarly low potential but suffers from significant volume expansion during lithiation (up to 400%), which damages the SEI and depletes the cell’s lithium inventory [11,12]. Lithium metal, offering a specific charge capacity of 3860 mAh g−1 and an electrode potential of 0 V (−3.04 V vs. the standard hydrogen electrode), presents another promising option, especially in “anode-free” cells where lithium metal forms in situ. However, challenges such as substantial SEI growth and dendritic lithium formation pose safety risks for lithium metal batteries (LMBs) [13]. To mitigate these issues, research is focused on developing new electrolytes and additives, protecting the current collector and lithium metal surface with appropriate coatings, and using structured host materials for lithium deposition. Evaluating interfacial side reactions, including SEI growth, is crucial for assessing the effectiveness of these strategies [14]. All-solid-state batteries (ASSBs) offer a promising avenue for safer lithium metal anode implementation, as they eliminate flammable liquid electrolytes [14]. The high mechanical strength and high lithium transference number (t(Li+) ≈ 1) of solid electrolytes (SEs) help reduce dendrite growth, making them suitable for lithium metal batteries [15]. However, SEs often have a narrow stability window and react with lithium metal at low potentials, forming unstable or stable SEIs. At high current densities, dendrite growth can occur along cracks and grain boundaries in inorganic SEs, exacerbating interfacial reactions and potentially causing short circuits and cell failure [16].
The SEI layer was originally identified and published by Peled in 1979 [17]. This layer is a consequence of specific reduction and oxidation reactions on the electrode surface. More precisely, when the cathode’s potential is lower and/or the anode’s potential is higher than the electrolyte’s highest occupied molecular orbital (HOMO) and/or lowest unoccupied molecular orbital (LUMO), the electrolyte molecules are oxidized and/or reduced on the electrode surface, respectively. The products of these oxidation and reduction reactions will contribute to the formation of the SEI layer [18,19,20]. It must be noted that the operating environment of a battery, such as the loaded voltage, current, and temperature, affect these HOMO/LUMO windows. The SEI production processes deconstruct the electrolyte molecules and deplete active Li-ions, leading to capacity fading [21]. An SEI layer that is appropriately designed inhibits electron tunneling and electrolyte diffusion at the electrode–electrolyte contact. This reduces the rate of SEI generation and reduction processes, and thus the rate of capacity fading [20]. A typical SEI layer consists of the following sections: (I) an inorganic inner layer that is closer to the electrode contact and is only permeable for Li-ions and (II) a highly permeable organic outer layer that is closer to the electrolyte interface and enables the transit of electrolyte molecules as well as Li-ions. In ASSBs, the first layer has similar thermodynamic and kinetic properties as in the presence of a liquid electrolyte and covers a passivating function towards the anode [15,22]. On the other hand, the second layer can be thermodynamically unstable [15], especially when a lithium anode is implemented, and will facilitate the continuous growth of the SEI layer during the operation of the battery. Consequently, characterization and comprehension of the SEI layer’s formation mechanisms and transport properties are crucial to improving ASSBs’ longevity.
Given that some of the SEI formation reactions take place at picosecond (ps) timescales, they cannot be investigated using the traditional experimental methods. Hence, multi-scale multi-physics modelling, which includes a combination of quantum mechanics (QM) calculations, molecular dynamics (MD) simulations, and macro-scale mathematical modelling, is implemented by researchers to study these mechanisms. In solid-state batteries, various aspects pertaining to the growth of the SEI layer are largely unknown because of its complex structure and unstable formation process [23]. Specifically, the main constraint stems from the absence of knowledge on the diffusion coefficient equation in the SEI layer, which is an issue that is also found in LIBs [24,25,26]. In the absence of an accurate mathematical formulation, Deng et al. [27] hypothesized a constant diffusion coefficient in their reduced order model. A more advanced pseudo-two-dimensional model has been proposed for ASSBs that continues to implement a fixed value for the diffusion coefficient over the anode active materials [28].
Clearly, there is a void in the current literature with regards to an accurate description of the diffusion mechanism in the SEI layer. Addressing this gap, this work aims to develop a model for solid-state batteries where a detailed diffusion equation is incorporated into the macro-scale model. The mathematical formulation presented in this work for ASSBs is inspired by our recent proposition for the liquid electrolyte-based LIBs [29]. More precisely, the model proposed in this work is an adaptation of the one developed by Ekström and Lindbergh [30], i.e., a macro-scale continuum mathematical model that quantifies the influence of SEI layer on the ageing of LIBs with a graphite anode material. Their model is a blend of kinetic and transport control systems and employs a constant diffusion coefficient. Additionally, their model requires three lumped fitting parameters, which are substituted in the equations rather than variables such as a diffusion coefficient. These parameters enhance the accuracy of the model with respect to the experimental data over a variety of temperature and concentration values. However, a major drawback of their model is that it needs to be tuned using battery-specific experimental data for investigating a battery design. Further, the fitting parameters will vary with the material of the LIBs, forbidding us from a computational study of new and unique materials in LIBs. Improving upon their model, in our earlier work [21], we combined quantum mechanics (QM) calculations and molecular dynamics (MD) simulations to propose an equation for the diffusion coefficient, dependent on temperature and Li-ion concentration, for each crystal structure in the SEI’s inner section. To accurately represent the physics within the SEI, this macro-scale mathematical model was refined with a single equation for the diffusion coefficient. Specifically, the diffusion coefficient equation from our previous work [21] was integrated into the macro-scale mathematical model (MSMM) developed by Ekström and Lindbergh [30], which is utilized in commercial engineering software such as Comsol Multiphysics, to analyze the capacity fading of LIBs. In implementing an anologous version in this research for the ASSBs, our formulation is verified in relation to the experimental data of ASSB, predicting the battery’s behaviour for a specific temperature under multiple discharging conditions. Adapting the Ekström and Lindbergh MSMM for ASSB, we present a model with just two simplified fitting parameters, omitting the most complex parameter in the primary model. The updated expression compensates for the influence of temperature and concentration on ageing. As part of model validation, it was employed to study the formation of SEI layers and the consequent capacity fading as a function of time and initial SOC, for an extended range of temperature and concentration conditions.
It must be noted that our earlier model for LIBs was designed for batteries that used graphite anodes [29]. However, in ASSBs, there is an acute shortage of experimental data on batteries with these anodes. This is because these anodes have been proven to not work efficiently with solid electrolytes [31,32,33]. Specifically, it has been shown that there is a loss of contact between the solid electrolyte and this anode. Consequently, a large part of battery research in the past years has been focused on the use of the so-called “soft” electrolytes. The study by Kobayashi [34], the one under our scrutiny, describes the adoption of graphite as an anode in conjunction with two classes of solid polymer electrolytes, a “soft” one and a “hard” one. This work focuses on the “soft” polymer, that will be addressed as “SPE1”. The experimental tests from the literature correspond to a 2032 half coin cell with the architecture [Graphite|SPE1|Li], and a prototype pouch cell [Graphite|SPE1|LiFePO4]. In the coin cell case, Li is employed as the counter electrode to permit Li intercalation in the graphite anode. The solid-polymer electrolyte is formed by a high-molecular-weight ether-based polymer (P(EO/MEEGE) = 88/12, DAISO). Ethylene oxide (EO) and 2-(2-methoxyethoxy) ethyl glycidyl ether (MEEGE) are its components. Lithium bis(trifluoromethylsulfonyl) imide (LiTFSI, 3 M) was incorporated in a molar ratio of [O]/[Li] = 16/1. Both P(EO/MEEGE) and LiTFSI were dissolved in acetonitrile (AN) [34]. Although specific details regarding the composition of the SEI layer are not provided, the main components that form the permanent layer of the SEI can be identified within the solid polymer electrolyte (SPE). An additional and significant outcome of this study is the validation of the developed model through a real driving cycle, comparing the model’s ageing data with the LIB experimental one.

2. The Theoretical Method and Computational Details

In ordinary LIBs, the inorganic inner layer consists of Li2CO3, LiF, and Li2O [15,21,35]. In contrast, the organic outer layer is made of dilithium ethylene glycol dicarbonate (Li2EDC) and ROLi (R depends on the solvent). The inner substrate is mainly composed of fixed materials, which do not depend on the electrolyte composition. This is where the diffusion mechanism is deeply investigated and highlighted in order to integrate it with the single-particle model previously described for SEI formation. Experimental research has established that Li2CO3 is an outcome of the conversion reaction of CoCO3 upon Li-ion insertion, when ethylene carbonate is present in the liquid electrolyte. However, Li2CO3 is unsteady from a thermodynamic point of view and will reduce to Li2C2 and Li2O [36]. Also, Li2C2 will contribute to other processes and create Li+, C2H2, and C [37]. Hence, Li2CO3 cannot be viewed as a permanent element in the inner area of the SEI layer. In the case of an ASSB, the equivalent inner layer composition is confirmed to have LiF and Li2O [22,23,34]. The diffusion equation of Li ions through the crystal structure of the SEI layer can be achieved by implementing a modified version of the Arrhenius law (Equation (2)), where the Li diffusion coefficient is a function of the concentration (C) and temperature (T). The total diffusion coefficient can be computed with the following relation [21]:
D t = δ Li 2 O D Li 2 O + δ LiF D LiF ,
Furthermore, recent advancements in the characterization of SEI layers have revealed that the interplay between inorganic and organic components significantly impacts the stability and performance of LIBs and ASSBs. In particular, understanding the dynamic evolution of the SEI layer during battery cycling is crucial. The continuous formation and repair of the SEI layer contribute to capacity fading over time. Advanced techniques such as in situ electron microscopy and spectroscopy are being employed to observe these changes at the atomic level, providing deeper insights into the mechanisms driving SEI formation and degradation. Additionally, the integration of computational modelling with experimental data is becoming increasingly important. By simulating the electrochemical environment and interactions within the SEI layer, researchers can predict the behaviour of new electrolyte formulations and SEI compositions. This synergy between theory and experiment accelerates the development of more robust SEI layers, ultimately leading to batteries with longer lifespans and higher efficiency. In summary, the ongoing research into the SEI layers of both LIBs and ASSBs is critical for enhancing battery performance. The integration of experimental studies, advanced characterization methods, and computational modelling provides a thorough approach to understanding and enhancing the SEI layer, setting the stage for future energy storage advancements.
In Equation (1), δ i and D i are the fraction of the surface area and the diffusion coefficient of the ith component, respectively. Additionally, for each species the diffusion coefficient can be computed as a function of temperature and concentration, implementing a modified version of the Arrhenius law:
D i ( C , T ) = D 0 e x p A 0 E B ( C ) k b T .
In the above equation, A 0 and D 0 are component dependent constants. Also, EB stands for the energy barrier, which can be expressed as a function of concentration as:
E B ( C L i ) = a 2 C L i + 2 + a 1 C L i + + a 0 ,
where the dependent constants, a 0 , a 1 , and a 2 , are summarized in Table 1.
Since an actual measurement of the fractional surface area of the individual constituents in the SEI layer would be extremely complex [21], the following formulation has been employed in this work:
ω i = m i m Li 2 O + m LiF ,
δ L i 2 O = ω L i 2 O ω L i 2 O + ρ L i 2 O ρ L i F ω L i F and δ L i F = ω L i F ω L i F + ρ L i F ρ L i 2 O ω L i 2 O O ,
In the above equations, m, L, ω , and ρ are the mass, thickness, mass fraction, and density of the material i (in this case, Li2O and LiF), respectively.

A Mathematical Model for Ageing

An upgraded version of the Ekström and Lindbergh model, compatible with a solid-state battery is proposed in this research. The model predicts battery ageing by utilizing data on the SEI growth rate, significantly improving the model’s accuracy. In particular, we provide a zero-dimensional model that employs a constant current control in battery cycling, eliminating the need to define a positive electrode or electrolyte. The accumulated charge (QSEI (C)), which is lost in ASSBs due to the SEI layer formation processes, is computed as follows:
d Q S E I d t = I S E I ,
where I S E I (A) is the current of the parasitic reactions involved in the SEI layer formation. It can be represented as the aggregate of the currents traversing the surfaces that are completely coated by either an intact or a fractured SEI layer. The mathematical derivation of I S E I , presented in the Appendix A, yields the following expression:
I S E I = 1 + H K c r d J I 1 C exp α F η S E I R T + f J Q S E I I 1 C ,
where I S E I and I 1 C denote the currents through the cracked parts and the 1 C-rate charging current, respectively. Furthermore, the original model by Ekström and Lindbergh incorporates the following three lumped fitting parameters:
J = ϵ c o v I 0 I 1 C ,
f = τ c o v V I 1 C 2 ϵ c o v 1 ϵ c o v C D F A 2 ,
H = a c r d ϵ c o v .
The parameters in these equations are described in the nomenclature section. Collectively, these parameters compensate for the lack of information on the growth of the SEI layer and the ageing processes. Further, these fitting parameters must be calibrated for each new battery cell and/or operational condition.
As stated previously, the SEI layer comprises various materials and crystal structures, namely, Li2O and LiF, that allow the diffusion of electrolyte solvent molecules and Li-ions, which are of different sizes and have different charge values. This implies a wide variation in the diffusion coefficients in these constituent structures, which will be governed by the operating temperatures and concentration conditions. Unfortunately, with experimental data for the complex diffusion and reaction processes, Equations (8)–(10) for J , f , and H are adopted for better accuracy and agreement with the experimental data. Since the accumulation of Li-ions increases the gradient of the charge distribution in the SEI layer, it increases the electron leakage, increasing I 0 . Conversely, an increased diffusion coefficient will reduce the gradient of Li-ions within the SEI layer. Specifically, I 0 is inversely related to D T [29]. Thus, in this study, by modifying Equation (8), we utilized the following expression for J, which incorporates a material constant ( J 0 ) and D T :
J = J 0 D T
In this equation, DT is as defined in Equation (1). J 0 , a material constant, is set to a value of 2.2 × 10 15 m 2 s . Thus, with this definition, the lumped fitted parameter (J), defined in Equation (8) is eliminated from the new model. Similarly, updated expressions for f and H parameters were obtained as follows: The parameters f and H were tuned using the experimental data of Kobayashi et al. [34]. Specifically, the values were tuned to the experimental condition of charge/discharge cycling with 1/8 C-rate and 1/16 C-rate load currents at 60 °C. In doing so, the value of J is as prescribed in Equation (11). The tuned values are reported in Table 2.
Although the equation for the diffusion coefficient (Equation (1)) employs the same formulation for all diffusing particles, it appropriately considers the effects of different concentrations, temperatures, and crystal structures on the diffusion coefficient and SEI formation. In the forthcoming section, we establish that our proposed model is accurate with respect to the experimental data, and forecasts the ASSB behaviour under diverse situations, thereby validating the new formulation that requires just two fitting parameters. Furthermore, after considering capacity fading resulting from Li-ion loss during SEI formation, the relative capacity (RC) can be calculated as follows:
RC = Q batt , 0 Q S E l Q batt , 0 ,
where Qbatt,0 is the initial battery capacity. Additionally, the thickness of the SEI layer can be expressed as:
s = Q S E I V 1 ϵ c o v A .
The initial charge accumulation (QSEI) is set to zero at t = 0 , and the relative capacity RC is accordingly set to 100%. The relative capacity for the previously highlighted ageing scenario is evaluated over a period of 200 days, employing the enhanced Ekström and Lindbergh model, and the results are then compared with the experimental data. Further, a parametric analysis has been conducted, supported by evidence and insights from the literature.

3. Results and Discussion

Model Validation: Due to the lack of sufficient experimental data, we investigated the capacity fade as a function of time for one temperature. The outcome of the investigation is presented in Figure 1. The results demonstrate the ability of the enhanced MSMM in forecasting a capacity fade, consistent with the experimental results for both full-cell and half-cell configurations. This is because of the two fitting parameters considered in this model, which can be adjusted for different cell types and chemistries. Experimental data are from a custom 2032 coin cell whose specific capacity is retrieved from the reference battery.
Furthermore, to demonstrate the ability of the model to adapt to different geometries and cathode compositions, the MSMM is calibrated with an upgraded set of f and H values reported in Table 2, and tested for the ASSB pouch cell made by Kobayashi et al. [34] (see Figure 1). As seen in this figure, there is an excellent agreement with the experimental data with an R 2 value of 0.99 in both cases.
Additionally, we also applied our enhanced model to a LIB 63 Ah Nickel/Manganese/Cobalt–Lithium Nickel Manganese Cobalt Oxide(NMC-LMO) battery examined through a Worldwide Harmonized Light Vehicle Test Procedure (WLTP) cycle in Ref. [38]. The testing current profile is retrieved by applying a backward modelling approach to the vehicle data. The predictions from our model along with the experimental data of Micari et al. [39] are shown in Figure 2. As seen in this figure, there is an excellent agreement between our model and the experimental data.
To further establish the accuracy and thereby the validity of our model, we studied the open-circuit storage conditions and a driving cycle test using our model. The results from these evaluations are presented in Figure 3 and Figure 4. To match the experimental data, the values of f and H were retuned for the two operating temperatures of 25 °C and 45 °C (Table 3).
As seen in Figure 4, the WLTP results are extremely satisfactory. A certain amount of divergence is registered in the first 2000 cycles because of the oscillatory nature of the experimental data [38]. To summarize, the determination of coefficient is reported in Table 4.
Parametric Analysis: The validated model has been used to evaluate the ageing process of the ASSB coin cell at different operating conditions. In absence of data for additional temperatures, the following temperature-dependent linear profile was assumed for the parameters f and H:
f ( T ) = f 0 T + f 1 ,
H ( T ) = H 0 T + H 1 ,
where f 0 and H 0 are summarized in Table 5. The choice of this linear profile is based on our earlier work [29] where we successfully applied an analogous model to study the effect of temperature on LFP commercial battery, assuming a linear temperature-dependent function for the parameters f and H. It must be noted that the values highlighted in our prior work [29] were for liquid electrolyte-based batteries and as such would not be applicable to ASSBs. In fact, with those constants, the produced results were not feasible and in the 40–50 °C range they did not demonstrate sufficient sensitiveness to temperature variation. Hence, with the intent of forecasting a reliable trend for the ASSB capacity fade, the linear profile coefficients, f 0 , f 1 , and H 1 were re-tuned (see Table 5). The value of H 0 was the same as in our earlier model [29].
First, the ASSB cycling capabilities were analyzed for various temperature levels. The battery performance was evaluated within a range limited by the available experimental data. Specifically, the effect of temperature in the range of [40 °C, 60 °C] on the ageing of a battery cycling at 1 / 8 C rate has been evaluated, and the results are summarized in Figure 5.
As seen in this figure, lowering the operating temperature significantly retards the ageing of a battery. Though we expect a higher lithiation temperature to not only enhance the ionic conductivity of the SPE but also mitigate the mechanical stress due to the softening effect, there is evidence in the literature that a raise in the temperature has a deleterious effect on solid polymer electrolyte, accelerating the ageing process [40]. Additionally, the electrochemical performance of solid-state lithium batteries significantly declines at extreme temperatures, leading to considerable energy and power losses, reduced cycling lifespan, and increased safety risks [41]. Our trends are in excellent agreement with these propositions. A similar result is obtained when the model is applied to study the effect on the battery cycling at 1 C-rate, the theoretical upper limit value at which the enhanced model has been properly tested and validated (see Figure 6). Comparing the ageing results from the different C-rates shown in Figure 6, it is clear that the raise of the operating c-rate has an adverse effect on the ageing of the solid-state battery. Once again, this is in agreement with the experimental work of Fang et al. [42]. The findings indicate that capacity loss in all-solid-state batteries (ASSBs) is initially driven by a reduction in the contact area between the electrolyte and electrode interfaces. Over time, this degradation is mainly influenced by the mobility of lithium ions within the ageing electrolyte. Furthermore, enhancing the discharge rate can magnify this deterioration [42]. More precisely, their results indicated a 10% increase in capacity loss when they raised the continuous discharging rate of 1 C.
In Figure 7a, the relative capacity is computed for an extended range of initial SOC. The slowest ageing rate was recorded near 50% SOC. However, the initial SOC has a minor influence on the rate of capacity loss with respect to the LIB case. This is due to the higher operating temperature, which mitigates the initial SOC leverage, a tendency consistent with our previous observations [29].
We also investigated the growth of an SEI layer that grows rapidly in the battery. As shown in Figure 7b, in the ASSB coin cell, the maximum layer thickness was 79 nm at the extremes of the SOC window. At 50% SOC, the SEI layer thickness is about 75 nm. This lack of variation in the thickness is consistent with the findings of Kobayashi et al. [34], who stated that unlike the liquid electrolyte systems, the SEI formation occurs over a wider voltage region in the SPE systems. Note that the SEI formation process is more evident in an ASSB because of the high difference in reactivity between the SPE and the graphite anode [34]. In Figure 7a, the relative capacity has been calculated as a function of time and initial state of charge (SOC). Due to the limited experimental data, it is not possible to compare different initial SOCs directly or to confirm the impact of initial SOC variations on battery ageing. However, for an initial SOC of 50%, the trend aligns well with the experimental data, achieving an accuracy of over 98%. Conversely, Figure 7b illustrates the SEI thickness over time. Unfortunately, experimental tests on battery systems do not report SEI thickness, making direct comparisons difficult. In the model, SEI thickness is directly correlated with capacity fade and is a primary factor in capacity retention. Given the high accuracy of capacity retention results relative to the experimental data, the SEI thickness results generated by this model can be considered reliable.

4. Further Investigations

We have reported the efficacy of an enhanced model applied to LIBs evaluated with a fast-varying discharge rate and to solid-state batteries. In the latter instance, implementing only two fitting parameters allows for estimating the ASSB’s electrochemical performance for diverse temperatures and concentrations. Nevertheless, it impels us to acquire new experimental evidence to absolutely validate our model. Furthermore, in the macro-scale model we have assumed the same diffusion profile through the whole SEI layer for all the particles [29], which could explain for both battery technologies the slight divergence from experimental data. On the ASSB’s side, an additional factor to be addressed concerns the composition of the SEI layer. Specifically, the outer SEI layer, whose species are significantly reliant on the solid electrolyte composition, has to be examined as it could experience an uncontrolled expansion. Hence, additional research is needed to gain a better understanding of the SEI layer and to provide a diffusion equation that gathers the effects of particle charge and size, and contact resistance for ASSBs, on the diffusion coefficient and SEI species.

5. Conclusions

We have presented an adaptation to solid-state batteries of the macro-scale model previously developed [29]. The calibration of the two modified lumped fitting parameters in conjunction with the upgraded equation of DT (Equation (1)) is shown to adequately predict the ASSB capacity decay as a function of time. This model is validated with regard to experimental data at a single temperature under charge/discharge cycling with 1/8 C-rate load current ageing condition. Additionally, acting only on the two lumped parameters allows for adjusting the model to a different ASSB geometry and cathode composition. The model estimates a capacity fade trend in solid agreement with actual data in the instance of an ASSB pouch cell cycled with a 1/16 C-rate charge/discharge load current. Further, the upgraded model can also determine the SEI formation rate as a function of temperature and Li-ion concentration. The LIB analysis gave further insights; in this case, the model is verified at a specific temperature for multiple ageing conditions: (I) open circuit state at 80% and 100% SOC and (II) WLTP cycling with a fast-varying load current. In conclusion, the macro-scale model, whose fitting parameters need to be assessed only for each distinct operating temperature, poses itself as a solid computational tool for accurately predicting the LIBs’ electrochemical performance and conducting a primary study on ASSBs; for the latter, this could be beneficial in the selection of optimal materials.

Author Contributions

Conceptualization, S.S.; methodology, P.S. and A.L.; software, P.S. and A.L.; validation, P.S. and A.L.; formal analysis, P.S. and A.L.; investigation, P.S. and A.L.; resources, S.S.; writing—original draft preparation, P.S. and A.L.; writing—review and editing, S.S.; supervision, S.S.; project administration, S.S.; funding acquisition, S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by NSERC-Canada through the Discovery grants program (grant number RGPIN-2022-04988).

Data Availability Statement

There are no new data in this research. All results are from simulations.

Acknowledgments

The authors are grateful to the reviewers for their constructive criticism to improve this manuscript. The research has been partially funded by NSERC-Canada through the Discovery grants program.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

IndexDescription
covSurfaces covered by an microporous SEI layer
crdSurfaces where the SEI layer has cracked
icalIntercalating reaction
negNegative electrode
SEISEI layer or SEI layer forming reaction
sSolid (Electrode) phase
lLiquid (Electrolyte) phase
SymbolUnitDescription
a c r d 1Proportionality factor
A m 2 Electrode surface area
C m o l m 3 Reactant concentration of SEI formation
C b a t t 1Relative capacity
D i m 2 s Diffusion coefficient
E e q , i VEquilibrium potential
f1Lumped fitting parameter
F C m o l Faraday’s constant, 96485
I 0 AExchange current
I 1 C A1C charge/discharge current
I k i n AKinetic current
I l i m ALimiting current
I l o a d AApplied current on the battery
J1Lumped fitting parameter
H1Lumped fitting parameter
Q b a t t . 0 CInitial battery capacity
Q n e g CCharge accumulated in the negative electrode
Q S E I CCharge depleted to SEI-forming reactions
R J m o l K Molar gas constant, 8.3145
smSEI layer thickness
TkTemperature
V m 3 C Coulombic volume for forming the SEI
x1Stoichiometric coefficient in LixC6
α 1Transfer coefficient
ϵ i 1Porosity
η i VOverpotential
Φ i VPotential of i phase
τ i 1Tortuosity

Appendix A

The entire current through an SEI layer can be expressed as:
I S E I = I c o v + I c r d
where I c o v and I c r d are the currents through the covered and cracked parts, respectively.
In addition, I cov ( or crd ) can be written relying on a first-order mass transfer limiting current through a Nernst boundary layer as:
I c o v ( o r c r d ) = I k i n , c o v ( or crd ) 1 + I k i n , c o v ( o r c r d ) I lim . cov ( o r c r d ) .
In the reported equation, I kin , cov ( or crd ) (A) is the kinetic current and I lim is the highest current of the SEI formation, restricted by the mass transport. These currents can be computed as:
I kin , cov ( o r c r d ) = ϵ cov ( o r c r d d ) I 0 exp α η S E I F R T
where I 0 denotes the exchange current and η S E I signifies the overpotential of the SEI formation reaction. The definitions of all other symbols and notations are provided in the nomenclature area. Additionally, ϵ c r d and ϵ c o v are computed as follows:
ϵ c r d = a c r d K c r d ,
ϵ c o v = 1 ϵ c r d ,
where K c r d and a c r d are the expansion factor and a dimensionless proportionality factor, respectively. Moreover, for K c r d , we have:
K cra = 2 I ical I 1 C I ical < 0 and x < 0.3 0 I ical < 0 and 0.3 x 0.7 I ical I 1 C I kcal < 0 and 0.7 x 0 I kcal 0
The overpotential of the SEI formation reaction can be obtained by imposing the SEI reaction equilibrium potential ( E e q , S E l ) to zero and calculating the difference between the electrode potential Φ s and electrolyte potential Φ l as follows:
η S E I = Φ s Φ l E e q , S E I
It should be noted that the liquid phase potential is explicitly associated with the equilibrium potential of the intercalation reaction, E e q , i c a l , as well as the corresponding overpotential, η i c a l , and can be formulated as follows:
Φ l = E eq , ical + η ical .
By inverting the Butler–Volmer equation, the relationship between the overpotential and the intercalation current can be reported as:
η ical = R T 0.5 F arcsinh I ical 2 k ical I 1 C ( ( 1 x ) x ) 0.5 ,
where x and I ical are the SOC and intercalation current of the anode, respectively. I 1 C is the battery’s nominal 1 C charge, and k ical is the proportionality factor that provides the overpotential. In the Nernst boundary layer, the limiting current density and the accumulated SEI layer thickness are inversely connected:
I lim . cov ( or crd ) = ϵ cov ( or crd ) C D cov ( or crd ) F A s
where C , F and A are the concentration, Faraday’s constant, and electrode surface area, respectively. D cov ( or crd ) is inversely related to the tortuosity of the layer as:
D cov ( or crd ) = D T τ c o v .
In conclusion, to mitigate discrepancies and refine the precision of the model, Ekström and Lindbergh [30] introduced the following three consolidated fitting parameters:
J = ϵ c o v I 0 I 1 C ,
f = τ c o v V I 1 C 2 ϵ c o v 1 ϵ c o v C D F A 2 ,
H = a c r d ϵ c o v .

References

  1. Shafique, M.; Azam, A.; Rafiq, M.; Luo, X. Life cycle assessment of electric vehicles and internal combustion engine vehicles: A case study of Hong Kong. Res. Transp. Econ. 2022, 91, 101112. [Google Scholar] [CrossRef]
  2. Janek, J.; Zeier, W.G. A solid future for battery development. Nat. Energy 2016, 1, 1–4. [Google Scholar] [CrossRef]
  3. Kato, Y.; Hori, S.; Saito, T.; Suzuki, K.; Hirayama, M.; Mitsui, A.; Yonemura, M.; Iba, H.; Kanno, R. High-power all-solid-state batteries using sulfide superionic conductors. Nat. Energy 2016, 1, 16030. [Google Scholar] [CrossRef]
  4. Kim, J.G.; Son, B.; Mukherjee, S.; Schuppert, N.; Bates, A.; Kwon, O.; Choi, M.J.; Chung, H.Y.; Park, S. A review of lithium and non-lithium based solid state batteries. J. Power Sources 2015, 282, 299–322. [Google Scholar] [CrossRef]
  5. Park, K.H.; Kaup, K.; Assoud, A.; Zhang, Q.; Wu, X.; Nazar, L.F. High-Voltage Superionic Halide Solid Electrolytes for All-Solid-State Li-Ion Batteries. Acs Energy Lett. 2020, 5, 533–539. [Google Scholar] [CrossRef]
  6. Winter, M.; Barnett, B.; Xu, K. Before Li Ion Batteries. Chem. Rev. 2018, 118, 11433–11456. [Google Scholar] [CrossRef] [PubMed]
  7. Diana, G.; Steve, G.; Kristina, E.; Stefano, P. Interfacial phenomena in lithium batteries and beyond. J. Power Sources 2023, 562, 232673. [Google Scholar]
  8. Heiskanen, S.K.; Kim, J.; Lucht, B.L. Generation and Evolution of the Solid Electrolyte Interphase of Lithium-Ion Batteries; Elsevier: Amsterdam, The Netherlands, 2019; Volume 3, pp. 2322–2333. [Google Scholar]
  9. Duffner, F.; Kronemeyer, N.; Tübke, J.; Leker, J.; Winter, M.; Schmuch, R. Post-lithium-ion battery cell production and its compatibility with lithium-ion cell production infrastructure. Nat. Energy 2021, 6, 123–134. [Google Scholar] [CrossRef]
  10. Randau, S.; Weber, D.A.; Kötz, O.; Koerver, R.; Braun, P.; Weber, A.; Ivers-Tiffée, E.; Adermann, T.; Kulisch, J.; Zeier, W.G.; et al. Benchmarking the performance of all-solid-state lithium batteries. Nat. Energy 2020, 5, 259–270. [Google Scholar] [CrossRef]
  11. Borodin, O.; Self, J.; Persson, K.; Wang, C.; Xu, K. Uncharted Waters: Super-Concentrated Electrolytes; Elsevier: Amsterdam, The Netherlands, 2020; Volume 4, pp. 69–100. [Google Scholar]
  12. Louli, A.J.; Eldesoky, A.; Weber, R.; Genovese, M.; Coon, M.; deGooyer, J.; Deng, Z.; White, R.T.; Lee, J.; Rodgers, T.; et al. Diagnosing and correcting anode-free cell failure via electrolyte and morphological analysis. Nat. Energy 2020, 5, 693–702. [Google Scholar] [CrossRef]
  13. Wu, Z.; Wang, C.; Hui, Z.; Liu, H.; Wang, S.; Yu, S.; Xing, X.; Holoubek, J.; Miao, Q.; Xin, H.L.; et al. Growing single-crystalline seeds on lithiophobic substrates to enable fast-charging lithium-metal batteries. Nat. Energy 2023, 8, 340–350. [Google Scholar] [CrossRef]
  14. Aktekin, B.; Riegger, L.M.; Otto, S.; Fuchs, T.; Henss, A.; Janek, J. SEI growth on Lithium metal anodes in solid-state batteries quantified with coulometric titration time analysis. Nat. Commun. 2023, 14, 6946. [Google Scholar] [CrossRef] [PubMed]
  15. Banerjee, A.; Wang, X.; Fang, C.; Wu, E.A.; Meng, Y.S. Interfaces and interphases in all-solid-state batteries with inorganic solid electrolytes. Chem. Rev. 2020, 120, 6878–6933. [Google Scholar] [CrossRef] [PubMed]
  16. Xiao, Y.; Wang, Y.; Bo, S.; Kim, J.C.; Miara, L.J.; Ceder, G. Understanding interface stability in solid-state batteries. Nat. Rev. Mater. 2020, 5, 105–126. [Google Scholar] [CrossRef]
  17. Peled, E. The Electrochemical Behavior of Alkali and Alkaline Earth Metals in Nonaqueous Battery Systems: The Solid Electrolyte Interphase Model. J. Electrochem. Soc. 1979, 126, 2047. [Google Scholar] [CrossRef]
  18. Xu, K. Nonaqueous Liquid Electrolytes for Lithium-Based Rechargeable Batteries. Chem. Rev. 2004, 104, 4303–4418. [Google Scholar] [CrossRef] [PubMed]
  19. Plett, G.L. Battery Management Systems. Volume II, Equivalent-Circuit Methods; Artech: Morristown, NJ, USA, 2015; p. 329. [Google Scholar]
  20. Wang, A.; Kadam, S.; Li, H.; Shi, S.; Qi, Y. Review on modeling of the anode solid electrolyte interphase (SEI) for lithium-ion batteries. Npj Comput. Mater. 2018, 4, 1–26. [Google Scholar]
  21. Lanjan, A.; Moradi, Z.; Srinivasan, S. Multiscale Investigation of the Diffusion Mechanism within the Solid-Electrolyte Interface Layer: Coupling Quantum Mechanics, Molecular Dynamics, and Macroscale Mathematical Modeling. ACS Appl. Mater. Interfaces 2021, 13, 42220–42229. [Google Scholar] [CrossRef]
  22. Leung, K.; Soto, F.; Hankins, K.; Balbuena, P.B.; Harrison, K.L. Stability of Solid Electrolyte Interphase Components on Lithium Metal and Reactive Anode Material Surfaces. J. Phys. Chem. C 2016, 120, 6302–6313. [Google Scholar] [CrossRef]
  23. Multi-scale Imaging of Solid-State Battery Interfaces: From Atomic Scale to Macroscopic Scale. Chem 2020, 6, 2199–2218. [CrossRef]
  24. Christensen, J.; Newman, J. A mathematical model for the lithium-ion negative electrode solid electrolyte interphase. J. Electrochem. Soc. 2004, 151, A1977–A1988. [Google Scholar] [CrossRef]
  25. Deng, J.; Wagner, G.J.; Muller, R.P. Phase Field Modeling of Solid Electrolyte Interface Formation in Lithium Ion Batteries. J. Electrochem. Soc. 2013, 160, A487–A496. [Google Scholar] [CrossRef]
  26. Liu, L.; Park, J.; Lin, X.; Sastry, A.M.; Lu, W. A thermal-electrochemical model that gives spatial-dependent growth of solid electrolyte interphase in a Li-ion battery. J. Power Sources 2014, 268, 482–490. [Google Scholar] [CrossRef]
  27. Deng, Z.; Hu, X.; Lin, X.; Xu, L.; Li, J.; Guo, W. A Reduced-Order Electrochemical Model for All-Solid-State Batteries. IEEE Trans. Transp. Electrif. 2021, 7, 464–473. [Google Scholar] [CrossRef]
  28. Toghyani, S.; Baakes, F.; Zhang, N.; Kühnelt, H.; Cistjakov, W.; Krewer, U. Model-Based Design of High Energy All-Solid-State Li Batteries with Hybrid Electrolytes. J. Electrochem. Soc. 2022, 169, 040550. [Google Scholar] [CrossRef]
  29. Lanjan, A.; Srinivasan, S. An Enhanced Battery Aging Model Based on a Detailed Diffusing Mechanism in the SEI Layer. ECS Adv. 2022, 1, 030504. [Google Scholar] [CrossRef]
  30. Ekström, H.; Lindbergh, G. A Model for Predicting Capacity Fade due to SEI Formation in a Commercial Graphite/LiFePO 4 Cell. J. Electrochem. Soc. 2015, 162, A1003–A1007. [Google Scholar] [CrossRef]
  31. Huu, N.; Phuc, H.; Hikima, K.; Muto, H.; Matsuda, A. Recent developments in materials design for all-solid-state Li-S batteries. Crit. Rev. Solid State Mater. Sci. 2021, 47, 283–308. [Google Scholar] [CrossRef]
  32. Otoyama, M.; Sakuda, A.; Hayashi, A.; Tatsumisago, M. Optical microscopic observation of graphite composite negative electrodes in all-solid-state lithium batteries. Solid State Ionics 2018, 323, 123–129. [Google Scholar] [CrossRef]
  33. Otoyama, M.; Kowada, H.; Sakuda, A.; Tatsumisago, M.; Hayashi, A. Operando Confocal Microscopy for Dynamic Changes of Li + Ion Conduction Path in Graphite Electrode Layers of All-Solid-State Batteries. J. Phys. Chem. Lett. 2020, 11, 900–904. [Google Scholar] [CrossRef]
  34. Kobayashi, Y.; Seki, S.; Mita, Y.; Ohno, Y.; Miyashiro, H.; Charest, P.; Guerfi, A.; Zaghib, K. High reversible capacities of graphite and SiO/graphite with solvent-free solid polymer electrolyte for lithium-ion batteries. J. Power Sources 2008, 185, 542–548. [Google Scholar] [CrossRef]
  35. Cresce, A.V.; Russell, S.M.; Baker, D.R.; Gaskell, K.J.; Xu, K. In situ and quantitative characterization of solid electrolyte interphases. Nano Lett. 2014, 14, 1405–1412. [Google Scholar]
  36. McShane, E.J.; McShane, E.J.; Colclasure, A.M.; Brown, D.E.; Brown, D.E.; Konz, Z.M.; Konz, Z.M.; Smith, K.; McCloskey, B.D.; McCloskey, B.D. Quantification of Inactive Lithium and Solid-Electrolyte Interphase Species on Graphite Electrodes after Fast Charging. ACS Energy Lett. 2020, 5, 2045–2051. [Google Scholar]
  37. Cabras, L.; Oancea, V.; Salvadori, A. A novel two-mechanism model for all solid state Li-ion batteries: Review and comparisons. arXiv 2021, arXiv:2104.05424. [Google Scholar]
  38. Micari, S.; Foti, S.; Testa, A.; De Caro, S.; Sergi, F.; Andaloro, L.; Aloisio, D.; Leonardi, S.G.; Napoli, G. Effect of WLTP CLASS 3B Driving Cycle on Lithium-Ion Battery for Electric Vehicles. Energies 2022, 15, 6703. [Google Scholar] [CrossRef]
  39. Micari, S.; Foti, S.; Testa, A.; De Caro, S.; Sergi, F.; Andaloro, L.; Aloisio, D.; Napoli, G. Reliability assessment and lifetime prediction of Li-ion batteries for electric vehicles. Electr. Eng. 2022, 104, 165–177. [Google Scholar] [CrossRef]
  40. He, Y.; Li, S.; Zhou, S.; Hu, H. Mechanical Integrity Degradation and Control of All-Solid-State Lithium Battery with Physical Aging Poly (Vinyl Alcohol)-Based Electrolyte. Polymers 2020, 12, 1886. [Google Scholar] [CrossRef] [PubMed]
  41. Li, Z.; Ren, Y.; Guo, X. Polymer-based electrolytes for solid-state lithium batteries with a wide operating temperature range. Mater. Chem. Front. 2023, 7, 6305–6317. [Google Scholar] [CrossRef]
  42. Fang, X.; He, Y.; Fan, X.; Zhang, D.; Hu, H. Modeling and Simulation in Capacity Degradation and Control of All-Solid-State Lithium Battery Based on Time-Aging Polymer Electrolyte. Polymers 2021, 13, 1206. [Google Scholar] [CrossRef]
Figure 1. The relative capacity over a duration of 200 days at 60 °C for the coin cell (cycling at 1/8 C-rate) and the pouch cell (1/16 C-rate). The experimental data are from Kobayashi et al. [34].
Figure 1. The relative capacity over a duration of 200 days at 60 °C for the coin cell (cycling at 1/8 C-rate) and the pouch cell (1/16 C-rate). The experimental data are from Kobayashi et al. [34].
Energies 17 02896 g001
Figure 2. A comparison of the predicted SOC values from the vehicle model and the experimental data of Micari et al. [39], for the WLTP cycles.
Figure 2. A comparison of the predicted SOC values from the vehicle model and the experimental data of Micari et al. [39], for the WLTP cycles.
Energies 17 02896 g002
Figure 3. Relative capacity over a duration of 700 days at 25 °C and 45 °C for a LIB-NMC-LMO-Prismatic battery. The experimental data are from Micari et al. [38].
Figure 3. Relative capacity over a duration of 700 days at 25 °C and 45 °C for a LIB-NMC-LMO-Prismatic battery. The experimental data are from Micari et al. [38].
Energies 17 02896 g003
Figure 4. The relative capacity over a duration of 4184 WLTP cycles at 25 °C for a LIB-NMC-LMO-Prismatic battery. The experimental data are from Micari et al. [38].
Figure 4. The relative capacity over a duration of 4184 WLTP cycles at 25 °C for a LIB-NMC-LMO-Prismatic battery. The experimental data are from Micari et al. [38].
Energies 17 02896 g004
Figure 5. Temperature sensitivity analysis: relative capacity over a duration of 200 days.
Figure 5. Temperature sensitivity analysis: relative capacity over a duration of 200 days.
Energies 17 02896 g005
Figure 6. Current sensitivity analysis: relative capacity over a duration of 200 days at 60 °C.
Figure 6. Current sensitivity analysis: relative capacity over a duration of 200 days at 60 °C.
Energies 17 02896 g006
Figure 7. The relative capacity and SEI thickness for different operating conditions in an ASSB.
Figure 7. The relative capacity and SEI thickness for different operating conditions in an ASSB.
Energies 17 02896 g007
Table 1. Coefficients of the second-degree polynomial described in Equation (3) for Li2O and LiF [21].
Table 1. Coefficients of the second-degree polynomial described in Equation (3) for Li2O and LiF [21].
a 0 [eV] a 1 [eV] a 2 [eV]
Li2O3.9488 8.9294 12.0460
LiF1.9886 2.5607 3.5237
Table 2. Parameters for the modified MSMM ageing model for ASSBs.
Table 2. Parameters for the modified MSMM ageing model for ASSBs.
ASSB FormatTemperaturefH
Coin cell60 °C1.5 × 1053.5
Pouch cell60 °C3.4 × 1042.5
Table 3. Parameters for the modified MSMM ageing model for LIBs.
Table 3. Parameters for the modified MSMM ageing model for LIBs.
TemperaturefH
25 °C3.1 × 10674.8
45 °C4.9 × 10542.1
Table 4. R 2 values of the computed relative capacity derived from the modified MSMM models in relation to the experimental results.
Table 4. R 2 values of the computed relative capacity derived from the modified MSMM models in relation to the experimental results.
TemperatureOperating Condition R 2 Value Using MSMM
25 °COpen Circuit-80% SOC 0.99
45 °COpen Circuit-100% SOC 0.98
25 °CWLTP 0.96
Table 5. The constants in Equations (14) and (15).
Table 5. The constants in Equations (14) and (15).
f 0 [ 1 K ] f 1 H 0 [ 1 K ] H 1
4.5 × 10 3 1.7 × 106 0.14 50.1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Scaltrito, P.; Lanjan, A.; Srinivasan, S. An Enhanced Ageing Model for Solid-State Batteries. Energies 2024, 17, 2896. https://doi.org/10.3390/en17122896

AMA Style

Scaltrito P, Lanjan A, Srinivasan S. An Enhanced Ageing Model for Solid-State Batteries. Energies. 2024; 17(12):2896. https://doi.org/10.3390/en17122896

Chicago/Turabian Style

Scaltrito, Paolo, Amirmasoud Lanjan, and Seshasai Srinivasan. 2024. "An Enhanced Ageing Model for Solid-State Batteries" Energies 17, no. 12: 2896. https://doi.org/10.3390/en17122896

APA Style

Scaltrito, P., Lanjan, A., & Srinivasan, S. (2024). An Enhanced Ageing Model for Solid-State Batteries. Energies, 17(12), 2896. https://doi.org/10.3390/en17122896

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop