Theoretical Study on the Influence of the Anharmonic Effect on the Ionic Conductivity and Thermal Stability of 8 mol% Yttria-Stabilized Zirconia Solid Electrolyte Material

YSZ is a promising material for resistive memory devices due to its high concentration of oxygen vacancies, which provide the high anion migration rates crucial for the manifestation of resistance switching in metal oxides. Therefore, investigating the ionic conductivity of YSZ is an important issue. The ionic conductivity and thermal stability of 8 mol% YSZ were studied using the theories and methods of solid-state physics and physical chemistry. The impact of anomalous atomic vibrations on the material was also explored, and the variation in the ion vibration frequency, electrical conductivity, and thermal stability coefficient of electrical conductivity with temperature was obtained. The results show that the ion conductivity of an 8 mol% YSZ solid electrolyte increases nonlinearly with temperature, with a smaller increase at lower temperatures and a larger increase at higher temperatures. Considering the anharmonic effect of ion vibrations, the electrolyte conductivity is higher than the result of the harmonic approximation, and the anharmonic effect becomes more significant at higher temperatures. Our research fills the gap in the current literature regarding the theoretical non-harmonic exploration of the ion conductivity and thermal stability factor of YSZ solid electrolytes. These results provide valuable theoretical guidance for the development and application of high-performance YSZ resistive memory devices in high-temperature environments.


Introduction
YSZ (yttria-stabilized zirconia) resistive memory has great potential for use in hightemperature environments, such as the aerospace, automotive, energy, and petrochemical industries [1,2]. In these fields, electronic devices operating in high-temperature environments require high-temperature stability and reliability. However, the thermodynamic properties of solid electrolytes, such as thermal expansion, can significantly impact their conductivity. Therefore, several studies have been conducted to investigate the electrical conductivity and stability of solid electrolyte materials. In terms of experimental research, Kwon [3], Kumar [4], Park [5], Borik [6] et al. prepared doped ZrO 2 solid electrolyte materials and measured the electrical conductivity of the prepared materials by electrical impedance spectroscopy, open-circuit potential measurement, and other methods, and they studied the effects of the defect structure, doping concentration, temperature, and other factors on the electrical conductivity of the material. In terms of theoretical research, Lima [7] studied the thermal stability of yttria-stabilized zirconia coatings obtained by plasma spraying. Glinchuk [8] studied the effect of surface tension on the activation energy of oxygen ion conduction in nanodynamics and calculated the activation energy of oxygen ion diffusion through oxygen vacancies. Zavodinsky [9,10] used the electron density functional method (gradient approximation) and pseudopotential method to study the conduction mechanism of the ionic conductivity of ZrO 2 solid electrolytes and the stability mechanism of cubic zirconia and zirconia nanoparticles with different stoichiometries. Hirschorn [11] introduced a comprehensive model for predicting and explaining the ionic conductivity of YSZ solid electrolytes. Song [12] discussed the effects of temperature and oxygen partial pressure on the electrical conductivity of YSZ electrolytes. Etschmann [13] presented a molecular dynamics study of the ionic conductivity of YSZ. Wang [14], Mahapatra [15], and Yildiz [16] studied the effects of the microstructure, grain size, and grain boundary on YSZ ionic conductivity. These references cover various aspects of YSZ ionic conductivity, including the effects of temperature, oxygen partial pressure, microstructure, and other factors. These studies are of great significance for understanding the electrochemical properties of YSZ and its application in high-temperature environments, such as ion batteries or solid oxide fuel cells.
Currently, research on the conductivity of ZrO 2 solid electrolytes mainly relies on experimental methods. However, these methods cannot provide an analytical formula for the conductivity of solid electrolyte materials, making it difficult to theoretically reveal the mechanism of conductivity changes and other conductivity properties. Although some studies have theoretically investigated the ionic conductivity of YSZ materials, they did not consider the influence of anomalous atomic vibrations, and there are certain gaps between the results obtained and actual measurements. Moreover, the variations in the ionic conductivity and thermal stability of YSZ solid electrolyte materials with temperature have not been extensively studied. In this paper, we aim to apply solid-state physics theory to study the temperature dependence of the conductivity of 8 mol% YSZ solid electrolyte and explore the influence of anharmonic effects to provide valuable theoretical guidance for the development and application of high-performance YSZ resistive memory devices in high-temperature environments. Our research can fill the gap in the current literature regarding the theoretical non-harmonic exploration of the ion conductivity and thermal stability factor of YSZ solid electrolytes. Furthermore, it can provide a theoretical basis for the application research of YSZ in solid oxide fuel cells, oxygen sensors, thermal barrier coatings, and other related fields.

Physical Model
From a mesoscopic perspective, the solid electrolyte can be viewed as a cylindrical thin film made up of YSZ particles with a larger particle size and a thickness of L, as shown in Figure 1, in contact with the electrodes. YSZ has a cubic crystal structure [3]. For computational simplicity, it is assumed that in 8 mol% YSZ, the ZrO 2 lattice consists of a complex lattice of cubic cells with a lattice constant of a, where the O atoms are located at the centers of 8 vertices and 6 faces, and the Zr atoms are located at the midpoints of 12 edges. The doped Y atoms replace the Zr atoms located at the body center. The crystal structures of undoped ZrO 2 and Y-doped ZrO 2 cells are shown in Figure 2a,b, respectively. the conduction mechanism of the ionic conductivity of ZrO2 solid electrolytes and the s bility mechanism of cubic zirconia and zirconia nanoparticles with different stoic ometries. Hirschorn [11] introduced a comprehensive model for predicting and explain the ionic conductivity of YSZ solid electrolytes. Song [12] discussed the effects of temp ature and oxygen partial pressure on the electrical conductivity of YSZ electrolytes. Ets mann [13] presented a molecular dynamics study of the ionic conductivity of YSZ. Wa [14], Mahapatra [15], and Yildiz [16] studied the effects of the microstructure, grain si and grain boundary on YSZ ionic conductivity. These references cover various aspects YSZ ionic conductivity, including the effects of temperature, oxygen partial pressure, m crostructure, and other factors. These studies are of great significance for understand the electrochemical properties of YSZ and its application in high-temperature envir ments, such as ion batteries or solid oxide fuel cells.
Currently, research on the conductivity of ZrO2 solid electrolytes mainly relies experimental methods. However, these methods cannot provide an analytical formula the conductivity of solid electrolyte materials, making it difficult to theoretically rev the mechanism of conductivity changes and other conductivity properties. Althou some studies have theoretically investigated the ionic conductivity of YSZ materials, th did not consider the influence of anomalous atomic vibrations, and there are certain ga between the results obtained and actual measurements. Moreover, the variations in ionic conductivity and thermal stability of YSZ solid electrolyte materials with tempe ture have not been extensively studied. In this paper, we aim to apply solid-state phys theory to study the temperature dependence of the conductivity of 8 mol% YSZ solid el trolyte and explore the influence of anharmonic effects to provide valuable theoreti guidance for the development and application of high-performance YSZ resistive memo devices in high-temperature environments. Our research can fill the gap in the curr literature regarding the theoretical non-harmonic exploration of the ion conductivity a thermal stability factor of YSZ solid electrolytes. Furthermore, it can provide a theoreti basis for the application research of YSZ in solid oxide fuel cells, oxygen sensors, therm barrier coatings, and other related fields.

Physical Model
From a mesoscopic perspective, the solid electrolyte can be viewed as a cylindri thin film made up of YSZ particles with a larger particle size and a thickness of L, as sho in Figure 1, in contact with the electrodes. YSZ has a cubic crystal structure [3]. For co putational simplicity, it is assumed that in 8 mol% YSZ, the ZrO2 lattice consists of a co plex lattice of cubic cells with a lattice constant of a, where the O atoms are located at centers of 8 vertices and 6 faces, and the Zr atoms are located at the midpoints of 12 edg The doped Y atoms replace the Zr atoms located at the body center. The crystal structu of undoped ZrO2 and Y-doped ZrO2 cells are shown in Figure 2a,b, respectively.   In the model, atoms are treated as particles that interact with each other and anharmonic vibrations near their equilibrium positions. It is assumed that the between the two nearest atoms is r, and the average interaction energy is u(r), w pands near the equilibrium position r0. When deviating from the equilibrium p represented by ξ = r − r0, its value is small:

Pt particles holes
where ε0 is the simple harmonic coefficient, and ε1 and ε2 are the first and secon monic coefficients, respectively. For transition metals such as zirconium, the interaction potential between Z where λD = 4(πFC/6)a is the average Debye wavelength, a is the lattice constant, represents the structure factor for the face-centered cubic crystal, g is related to b n is the bond strength parameter, and r0 is the equilibrium distance between the tw est atoms, and the parameter values of Zr and Y [18,19] are seen in Table 1.  represents the repulsive interaction range parameter, which can be set at 5 times th of the Zr ion. λ represents the repulsive energy parameter, with a value of λZ0 = 2 erg [18]. The specific values of these parameters are provided in Table 2. In the model, atoms are treated as particles that interact with each other and undergo anharmonic vibrations near their equilibrium positions. It is assumed that the distance between the two nearest atoms is r, and the average interaction energy is u(r), which expands near the equilibrium position r 0 . When deviating from the equilibrium position, represented by ξ = r − r 0 , its value is small: where ε 0 is the simple harmonic coefficient, and ε 1 and ε 2 are the first and second anharmonic coefficients, respectively. For transition metals such as zirconium, the interaction potential between Zr atoms [17] is where λ D = 4(πF C /6)a is the average Debye wavelength, a is the lattice constant, F C = 1/4 represents the structure factor for the face-centered cubic crystal, g is related to bonding, n is the bond strength parameter, and r 0 is the equilibrium distance between the two nearest atoms, and the parameter values of Zr and Y [18,19] are seen in Table 1. For Zr-O ionic bonding, the Zr-O atomic interaction potential [18] is: where q represents the charge of the zirconium ion, and α = 1.6381 is the Madelung constant. z 0 represents the number of ions in the nearest neighbor of the Zr ion, while ρ represents the repulsive interaction range parameter, which can be set at 5 times the radius of the Zr ion. λ represents the repulsive energy parameter, with a value of λ Z0 = 2.1 × 10 −8 erg [18].
The specific values of these parameters are provided in Table 2. Table 2. Born-Mayer-Buckingham potential parameters [18][19][20]. The interaction potential for O-O can be written in the form of Van der Waals potential [21]:

Ion Pairs
where D is the well depth, and r 0 is the distance between two atoms when u(r) = 0. a, b, r 0 , and D can be determined by the specific potential energy curve.
Based on the structure shown in Figure 2 and Equations (2)- (4), the values of ε 0 , ε 1 , and ε 2 were obtained by considering the nearest neighbor, second nearest neighbor, and next nearest neighbor, as shown in Table 3.

Conduction Mechanism of YSZ Electrolyte
ZrO 2 is an ionic conductor crystal that conducts electricity through oxygen ion vacancies in its crystal lattice. Oxygen molecules are adsorbed on the electrode surface and transform into oxygen atoms. These atoms then diffuse to the three-phase interface on the electrode surface, where they become oxygen ions due to electrochemical reactions. These ions enter the YSZ electrolyte, migrate, and diffuse in the YSZ electrolyte along the path shown in Figure 3, before moving to the other electrode through the YSZ electrolyte to complete the circuit, thereby forming an electrical current. The addition of oxides generates a large number of oxygen ion vacancies in ZrO 2 . For Y 2 O 3 , every time two Y ions are added, one oxygen ion vacancy is generated. The concentration of oxygen ion vacancies mainly depends on the composition and quantity of the dopant [22]. tial [21]: where D is the well depth, and r0 is the distance between two atoms when and D can be determined by the specific potential energy curve. Based on the structure shown in Figure 2 and Equations (2)-(4), the and ε2 were obtained by considering the nearest neighbor, second neares next nearest neighbor, as shown in Table 3.

Conduction Mechanism of YSZ Electrolyte
ZrO2 is an ionic conductor crystal that conducts electricity through cancies in its crystal lattice. Oxygen molecules are adsorbed on the electr transform into oxygen atoms. These atoms then diffuse to the three-phase electrode surface, where they become oxygen ions due to electrochemical ions enter the YSZ electrolyte, migrate, and diffuse in the YSZ electrolyte shown in Figure 3, before moving to the other electrode through the YS complete the circuit, thereby forming an electrical current. The addition o ates a large number of oxygen ion vacancies in ZrO2. For Y2O3, every tim added, one oxygen ion vacancy is generated. The concentration of oxyge mainly depends on the composition and quantity of the dopant [22]. The chemical reaction that occurs in the ZrO2 solid electrolyte with Y2 can be expressed as [23]: The chemical reaction that occurs in the ZrO 2 solid electrolyte with Y 2 O 3 as a stabilizer can be expressed as [23]: The formation of oxygen ion vacancies can be described by the Kroger-Yink mechanism: where O 2− represents the oxygen ion in the lattice, V O represents the oxygen ion vacancy, Zr 4+ represents the zirconium ion, and Y 3+ represents the yttrium ion.
In addition to O 2− , Zr 4+ , and Y 3+ , the charged particles in the electrolyte, electrons and holes h, also contribute to its conductivity. Assuming that the conductivities are σO 2− , σ e , and σ h , respectively, the total electrical conductivity of the electrolyte can be expressed as σ = σO 2− + σ e + σ h . However, at temperatures below 1200 K, only oxygen ion conduction plays a major role. Therefore, the total electrical conductivity can be simplified to σ = σO 2− .

Ionic Conductivity in YSZ
The conduction of solid electrolytes is primarily ionic conduction, which results from the migration and diffusion of oxygen ions. As mentioned earlier, the addition of oxides generates a significant number of oxygen ion vacancies in ZrO 2 . Due to thermal excitation, when the ion energy exceeds the potential barrier between adjacent ions (referred to as the diffusion activation energy), the ion transitions from its original equilibrium position to the adjacent vacancy, creating vacancies at the original position. When the ions migrate and diffuse to the right, the vacancies migrate in the opposite direction.
Using the theory of solid-state physics and the vacancy diffusion ion conduction mechanism, the migration speed of oxygen ions along the direction of the electric field can be calculated [23]: For general sensors, the electric field strength E involved is not large, Eqd k B T: where µ is the mobility, and its relationship with temperature is obtained from Equation (8) as where γ represents the geometric factor, ν 0 represents the vibrational frequency of the vacancy (i.e., the vibrational frequency of the ion), d represents the transition distance, which is approximately equal to the lattice constant a, q represents the charge of the oxygen ion, and ε represents the change in the Gibbs function caused by the vacancy transition, i.e., the diffusion activation energy. Considering the ion's anomalous vibration, the relationship between the vibration frequency of the ion and the temperature [23] is Assuming that the number of oxygen ions per unit volume is C i , the current density is obtained from the speed of the oxygen ion migration in the direction of the electric field using Equation (8): By substituting Equation (10) into Equation (9), and then into Equations (8) and (11), and considering only the oxygen ion conductivity, the relationship between the electrolyte conductivity and temperature can be expressed as The thermal expansion coefficient is a physical quantity that describes the proportional change in the length, area, or volume of a material with temperature. It is a measure of the deformation of the material during thermal expansion, indicating the ratio between the change in length, area, or volume and the original length, area, or volume. The thermal expansion coefficient is an important property in materials, especially at high temperatures, as the coefficient increases with temperature and may cause the material to deform or crack under heat. Therefore, for materials and devices that operate at high temperatures, controlling and optimizing the thermal expansion coefficient is critical. The thermal stability of the material is often described by the temperature stability coefficient, and the temperature stability coefficient α σ of the conductivity of the electrolyte material is defined as [24]

Results and Discussion
In ref. [25], the masses of O, Zr, and Y atoms are reported as M O = 0.265686 × 10 −25 kg, M Zr = 1.514811 × 10 −25 kg, and M Y = 1.476357 × 10 −25 kg, respectively. Furthermore, the vibration frequencies at absolute zero, ω 0O = 9.375247 × 10 13 s −1 , ω 0zr = 7.853824 × 10 13 s −1 , and ω 0Y = 6.198625 × 10 13 s −1 , of O, Zr, and Y, respectively, are obtained from ω = (ε 0 /M) 1/2 . By substituting the above data into Equation (10), the change in the oxygen ion vibration frequency ν = ω/2π with temperature is obtained, as shown in Figure 4. Curves 0, 1, and 2 represent the results of simple harmonic approximations considering only the first anharmonic term and considering both the first and second anharmonic terms, respectively. The same notation is used below.  It can be seen that the vibration frequency of oxygen ions remains constant simple harmonic approximation. However, after considering the anharmonic eff ions, the vibration frequency of oxygen ions increases linearly with increasing tem It can be seen that the vibration frequency of oxygen ions remains constant in the simple harmonic approximation. However, after considering the anharmonic effect of ions, the vibration frequency of oxygen ions increases linearly with increasing temperature. The vibration frequency of only the first anharmonic term is greater than that of considering both the first and second anharmonic terms, and it is also greater than that of the simple harmonic approximation. Furthermore, the higher the temperature, the greater the difference between the vibration frequencies, indicating a more significant anharmonic effect. The influence of the anharmonic effect on the vibrational frequency of oxygen ions increases with increasing temperature. For example, when T is 700 K, the influence of the anharmonic effect on the vibrational frequency of the oxygen ion is 2.14%, while when T is 1200 K, it is 3.62%.
The anharmonic effect arises due to the nonlinearity of the interatomic potential energy function. In the simple harmonic approximation, the potential energy function is approximated as a symmetric parabolic well, and the motion of ions is assumed to be harmonic. This approximation works well at low temperatures when the amplitude of motion is small. However, at higher temperatures, the amplitude of motion is larger, and the potential energy function becomes asymmetric, leading to anharmonic behavior. The anharmonicity causes the interatomic potential energy function to deviate from the parabolic well, resulting in the variation in the vibrational frequency with increasing temperature. The anharmonic effect is more significant at higher temperatures, leading to a larger deviation from the simple harmonic approximation. Therefore, at high temperatures, the anharmonic effect becomes more pronounced, leading to an increase in the vibrational frequency of oxygen ions.
Ref. [6] provides the activation energy of oxygen ions doped with 8 mol% YSZ. In the temperature range of 300~1500 K, the lattice constant of the YSZ electrolyte with a doping concentration of 8 mol% varies with temperature approximately as a = (5.1218 + 4.2056 × 10 −5 T) Å, as reported in ref. [26]. The quantity of oxygen ions is q = 2e, and the number of oxygen ions per unit volume is C i = 2.208 × 10 27 . By substituting the above data into Equation (12), the change in the electrolyte's conductivity with temperature and the conductivity curve of the YSZ crystal with a doping concentration of 8 mol% given in ref. [26] are shown in Figure 5a. For comparison, the research results of ref. [27] are presented in Figure 5b. The consistency of the results in this paper with the changes in YSZ electrolyte conductivity with temperature reported in other documents, as indicated by refs. [26,27], further supports the validity of the calculations in this paper.
The conductivity of 8 mol% YSZ increases nonlinearly with increasing temperature. The conductivity is low when T is low, but it increases rapidly with temperature when T > 1000 K. For instance, when 800 K < T < 1000 K, the conductivity changes from 0.0024 S·cm −1 to 0.025 S·cm −1 , which is an increase of 0.0226 S·cm −1 . When 1000 K < T < 1200 K, it changes from 0.025 S·cm −1 to 0.115 S·cm −1 , which is an increase of 0.090 S·cm −1 . After taking into account the anharmonic effect of ion vibrations, the conductivity of the electrolyte is greater than that of the harmonic approximation. Moreover, the higher the temperature, the greater the difference in conductivity between the anharmonic and harmonic approximations. This indicates that the anharmonic effect becomes more significant at higher temperatures. When T is 1173 K, the conductivity of 8 mol% YSZ electrolyte is 0.096 S·cm −1 in this paper, while it is 0.063 S·cm −1 in ref. [6].
The conductivity of a solid electrolyte is related to the mobility of ions, which increases with temperature due to increased thermal energy. However, the conductivity does not increase linearly with temperature, and it may exhibit nonlinear behavior due to various factors, such as the concentration of defects, the crystal structure, and the anharmonic effect of ion vibrations. In the case of 8 mol% YSZ, the conductivity is low at low temperatures, but it increases rapidly with temperature when T > 1000 K. This nonlinear behavior can be attributed to the increased mobility of ions at higher temperatures, which leads to an increase in conductivity. However, the anharmonic effect of ion vibrations can also contribute to this behavior. As the temperature increases, the anharmonic effect becomes more significant, causing the ion vibrations to deviate from the harmonic approximation. This leads to an increase in conductivity compared to the harmonic approximation. The difference in conductivity between the anharmonic and harmonic approximations becomes greater at higher temperatures, indicating that the anharmonic effect becomes more significant at higher temperatures. Therefore, taking into account the anharmonic effect of ion vibrations is necessary for accurate predictions of conductivity at high temperatures. The conductivity of 8 mol% YSZ increases nonlinearly with increasing tempe The conductivity is low when T is low, but it increases rapidly with temperature > 1000 K. For instance, when 800 K < T < 1000 K, the conductivity changes from S•cm −1 to 0.025 S•cm −1 , which is an increase of 0.0226 S•cm −1 . When 1000 K < T < 12 changes from 0.025 S•cm −1 to 0.115 S•cm −1 , which is an increase of 0.090 S•cm −1 . Afte into account the anharmonic effect of ion vibrations, the conductivity of the electr greater than that of the harmonic approximation. Moreover, the higher the tempe Substituting Equation (12) into Equation (13), the variation curve of the conductivity's temperature stability coefficient of the YSZ electrolyte with temperature can be obtained, as shown in Figure 6.

temperatures.
Substituting Equation (12) into Equation (13), the variation curve of the co ity's temperature stability coefficient of the YSZ electrolyte with temperature c tained, as shown in Figure 6. The temperature stability coefficient ασ of YSZ electrolyte conductivity d nonlinearly with increasing temperature. As the temperature continues to rise, α be constant, indicating that the higher the temperature, the smaller the temper bility coefficient of the conductivity. This suggests that the thermal stability of the tivity of the YSZ electrolyte material is better. For example, when the temperatu K < T < 1000 K and 1000 K < T < 1200 K, the temperature stability coefficient of e conductivity decreases by 52.6% and 31.5%, respectively. The harmonic approxi consistent with the temperature stability coefficient of the electrolyte conductiv the anharmonic effect is considered, indicating that the anharmonic effect has li on the temperature stability of the conductivity.
The temperature stability coefficient of YSZ electrolyte conductivity is a m how much the conductivity changes with temperature. A small temperature sta efficient indicates that the conductivity is less sensitive to temperature changes, desirable for practical applications. In the case of YSZ electrolyte conductivity, perature stability coefficient decreases nonlinearly with increasing temperature havior can be attributed to the fact that at higher temperatures, the mobility o creases, which leads to a larger change in conductivity with temperature. Howev temperature continues to rise, the effect of increased ion mobility on conductivity less significant, resulting in a more stable conductivity with respect to tempera The temperature stability coefficient α σ of YSZ electrolyte conductivity decreases nonlinearly with increasing temperature. As the temperature continues to rise, α σ tends to be constant, indicating that the higher the temperature, the smaller the temperature stability coefficient of the conductivity. This suggests that the thermal stability of the conductivity of the YSZ electrolyte material is better. For example, when the temperature is 800 K < T < 1000 K and 1000 K < T < 1200 K, the temperature stability coefficient of electrolyte conductivity decreases by 52.6% and 31.5%, respectively. The harmonic approximation is consistent with the temperature stability coefficient of the electrolyte conductivity when the anharmonic effect is considered, indicating that the anharmonic effect has little effect on the temperature stability of the conductivity.
The temperature stability coefficient of YSZ electrolyte conductivity is a measure of how much the conductivity changes with temperature. A small temperature stability coefficient indicates that the conductivity is less sensitive to temperature changes, which is desirable for practical applications. In the case of YSZ electrolyte conductivity, the temperature stability coefficient decreases nonlinearly with increasing temperature. This behavior can be attributed to the fact that at higher temperatures, the mobility of ions increases, which leads to a larger change in conductivity with temperature. However, as the temperature continues to rise, the effect of increased ion mobility on conductivity becomes less significant, resulting in a more stable conductivity with respect to temperature. The decrease in the temperature stability coefficient of electrolyte conductivity with increasing temperature suggests that the thermal stability of the YSZ electrolyte material is better at higher temperatures. This is because the conductivity of the material changes less with temperature, making it more reliable for practical applications at high temperatures. The fact that the harmonic approximation is consistent with the temperature stability coefficient of the electrolyte conductivity when the anharmonic effect is considered indicates that the anharmonic effect has little effect on the temperature stability of the conductivity. This suggests that the anharmonic effect mainly affects the absolute value of the conductivity, rather than its temperature dependence.

Conclusions
This paper investigates the variation in ionic conductivity in 8 mol% YSZ solid electrolyte with temperature and the effect of the anharmonic vibration of ions by applying the theory and methods of solid-state physics. The results indicate the following: 1. The ionic conductivity of 8 mol% YSZ solid electrolyte increases nonlinearly with the increase in temperature. 2. The conductivity is low at low temperatures but increases rapidly with temperature when the temperature is higher than 1000 K. 3. After considering the anharmonic effect of ionic vibrations, the ionic conductivity is larger than the result of simple harmonic approximation, and the higher the temperature, the greater the difference between the anharmonic and simple harmonic approximations. This suggests that the anharmonic effect becomes more significant at higher temperatures. 4. These findings are consistent with other results in the literature on the variation in the conductivity of 8 mol% YSZ electrolyte with temperature, which demonstrates the validity of the calculation. By understanding the electrical conductivity of YSZ materials at different temperatures, researchers can better optimize the design and performance of YSZ-based resistive memory devices for specific applications. However, this paper also has some shortcomings: the studied YSZ solid electrolyte has a concentration of 8 mol%, and different concentrations of YSZ solid electrolytes may have different performance characteristics. So the results of this study may not be applicable to other concentrations of YSZ solid electrolytes. Furthermore, this study mainly focused on basic properties, such as ionic conductivity and thermal stability, and did not deeply investigate other properties relevant to practical application scenarios, such as reliability and durability. Therefore, a more comprehensive consideration is needed when applying these results to the design of practical devices.