Zentropy theory for accurate prediction of free energy, volume, and thermal expansion without fitting parameters
Based on statistical mechanics, a macroscopically homogeneous system, i.e., a single phase in the present context, is composed of many independent configurations that the system embraces. The macroscopical properties of the system are determined by the properties and statistical probabilities of those configurations with respect to external conditions. The volume of a single phase is thus the weighted sum of the volumes of all configurations. Consequently, the derivative of the volume to temperature of a single phase depends on both the derivatives of the volumes of every configuration to temperature and the derivatives of their statistical probabilities to temperature, with the latter introducing nonlinear emergent behaviors. It is shown that the derivative of the volume to the temperature of the single phase can be negative, i.e., negative thermal expansion, due to the symmetry-breaking non-ground-state configurations with smaller volumes than that of the ground-state configuration and the rapid increase of the statistical probabilities of the former, and negative thermal expansion can be predicted without fitting parameters from the zentropy theory that combines quantum mechanics and statistical mechanics with the free energy of each configuration predicted from quantum mechanics and the partition function of each configuration calculated from its free energy.
The stability criteria of a macroscopically homogeneous system require that the derivatives between conjugate variables in the combined law of thermodynamics must be positive except at the limit of stability where they become zero. However, there is no such requirement for derivatives between non-conjugate variables. The derivative of volume (V) to temperature (T) under constant pressure, i.e., thermal expansion, is thus not required to be positive in a stable system because V and T are not conjugate variables. Although the derivative of volume to temperature of a macroscopically homogeneous material, i.e., a single phase, is commonly thought to be positive, i.e., positive thermal expansion (PTE), it is known that the volume of water does decrease with the increase of temperature at temperature below 4 ℃, i.e., negative thermal expansion (NTE). The first manmade NTE material was an Fe65Ni35 alloy in 1897 by Guillaume with near zero thermal expansion (ZTE) at room temperature, commonly referred to as INVAR.
While there have been many theories that aim to describe NTE mechanisms in various materials, they are mostly phenomenological and focused on interpreting experimental observations[2–9]. Nevertheless, a fundamental understanding and a predictive theory without experimental inputs applicable to all materials are still lacking. Based on the Maxwell relation, the derivative of V to T under constant pressure equals the negative derivative of entropy to pressure under constant temperature. Therefore, the fundamental understanding of NTE is related to the entropy of the system. While the total entropy of a system can be accurately obtained as a function of temperature from integration of experimentally measured heat capacity, its pressure dependence is more difficult to comprehend, and its theoretical prediction remains elusive. In the present paper, the efforts by the author's group to accurately predict entropy of a single phase as a function of temperature and pressure are discussed in terms of the recently termed zentropy theory, along with its predictive capability of PTE, ZTE, and NTE without experimental inputs.
STATISTICAL MECHANICS AND ENTROPY
Based on statistical mechanics, a single phase at finite temperature is composed of various independent configurations that are in statistical equilibrium with each other and its surroundings under given external constraints. The probability of each configuration is related to its own and the system's partition functions as follows
For a closed system under hydrostatic pressure (
For a canonical ensemble, the combined law is written in terms of Helmholtz energy as follows
The partition functions of the system and its configurations are written as
When there is only one configuration in the system, Eq. 6 becomes
ZENTROPY THEORY AND DFT-BASED CALCULATIONS
For configurations with non-zero entropy, it is necessary to add their contributions to the total entropy of the system. Our zentropy theory is schematically shown in Figure 1 with the following equation for the total entropy of the system[10,14,15]
The first summation in Eq. 8 reflects the bottom-up view of the system by considering the contribution from individual configurations, and the second summation represents the top-down view of the system, seeing the statistical fluctuations of configurations, as discussed in Section 2 above. This nested formula of the zentropy theory can be extended to consider more complex systems such as forests, planets, and black holes with more degrees of freedom, with Eq. 8 representing one of the subsystems of the system[14,15]. This nested formula can also be extended in another direction to configurations with fewer degrees of freedom within the configuration
The Helmholtz energy of the system can thus be obtained as
Eq. 8 to Eq. 11 reduce to standard statistical mechanics when
Formulated as an exact theory of many-body systems, DFT articulates that for an interacting electron gas, there exists a universal function of electron density such that the energy is at its minimum value, i.e., the ground-state energy with a unique ground-state electron density. The numerical solution is formulated by explicitly separating the independent-electron kinetic energy and long-range Coulomb interaction energy and replacing the many-body electron problem using independent valence electrons with an exchange-correlation functional of the electron density and an associated exchange-correlation energy and potential, i.e., coarse-graining of electrons. Kohn and Sham used the finite temperature generalization of ground-state energy of an interacting inhomogeneous electron gas by Mermin and formulated the entropy of thermal electrons at finite temperatures. Wang et al. added the vibrational contribution and presented the Helmholtz energy as follows
As the electron and phonon degrees of freedom are included in each configuration from the DFT-based calculations, all configurations are ergodic and symmetry-breaking in terms of magnetic spin, electric polarization, atomic short-range ordering, and defects such as vacancy, dislocation, and stacking faults. For their
DERIVATIVE OF VOLUME TO TEMPERATURE
The pressure of a system can be calculated from the derivative of Helmholtz energy to volume as follows
As mentioned in our previous publications[10,19], we pointed out that one may consider the Gibbs ensemble, i.e., constant
From the Maxwell relation, one has
The statistical mechanics in terms of a Gibbs ensemble is shown below
The first term, represented by Eq. 21, is a linear combination of contributions from each configuration, and the second term gives nonlinear behavior. The condition for
As all variables in Eq. 22 are positive, and there are no contributions from
Both cases are shown in Figure 2 for Ce and Fe3Pt, respectively, in the temperature-volume phase diagrams. For Ce, the volumes of antiferromagnetic (AFM) and ferromagnetic (FM) symmetry-breaking non-ground-state configurations are larger than that of the non-magnetic (NM) ground-state configuration, while for Fe3Pt, the FM ground-state configuration has the largest volume. There is a critical point in both systems where the stable high-temperature single phase becomes unstable and separates into two phases with the same crystal structure and different molar quantities such as volume, entropy, and magnetic spin configurations. The two phases at lower temperature and the stable single phase at high temperature are all composed of the same configurations, and the only difference among them is the probabilities of various configurations. As each configuration in all phases is stable, the instability of the macroscopically homogeneous single phase originates from the competition of various configurations when viewed from high temperature rather than the conventional interpretation of phonon softening that considers phonon being stationary. While viewed from low temperature, the macroscopically homogeneous single-phase results from the mixture of two macroscopically homogeneous phases, and there is no instability involved.
Figure 2. Temperature-volume phase diagrams with isobaric volumes at various pressures of (A) Ce and (B) Fe3Pt. The volume (V) is normalized to their respective equilibrium volume (
ANHARMONICITY AND EMERGENT BEHAVIORS
Anharmonicity is usually represented by the deviation of entropy or heat capacity away from quasiharmonic behavior. It is noted in Figure 2 that the volume change for Fe3Pt, at a given pressure from 0 K to 600 K, is rather small, supported by the experimental data, as shown by the symbols on the isobaric volume curve under the ambient pressure. This indicates that the quasiharmonic approximation can give an accurate prediction of the entropy of each configuration. From Eq. 8, it can be seen that the first summation is the linear combination of entropies of individual configurations, and the emergent behaviors, i.e., the behaviors that none of the individual configurations possess, originate from the second summation in the equation. This is the same for the derivative of volume to temperature shown by Eq. 21 where the emergent behavior of NTE is due to the rapid increase of the symmetry-breaking non-ground-state configurations, and their volumes are smaller than that of the ground-state configuration, and none of the individual configurations possess NTE.
For a stable system, the derivatives between a potential and its conjugate molar quantity in the combined law are positive, i.e.,
When these derivatives become zero, the macroscopic system reaches its limit of stability and the extreme of anharmonicity, and the inverses of these derivatives diverge positively, i.e.,
Eq. 26 represents the heat capacity under constant pressure and can be derived from Eq. 8 as follows
Again, the first summation is the linear combination of the heat capacity of each configuration, and the nonlinear emergent behavior comes from the second summation.
However, derivatives between non-conjugate variables are not required to be positive, such as the derivative between volume and temperature. As they will also diverge at the limit of stability, they could be either positive or negative, i.e.,
As discussed above, the negative divergency occurs when the volume of the ground-state configuration is larger than those of non-ground-state configurations. The NTE spreads to single-phase regions far away from the critical point, as shown in Figure 2B, with significant anharmonic behaviors.
More recently, the zentropy theory was applied to predict the ferroelectric-paraelectric (FE-PE) transitions in PbTiO3 with three configurations considered, i.e., tetragonal polarized configurations without domain wall (FEG), with 90o domain wall (90 DW), and with 180o domain wall (180 DW). With two sets of domain wall energies predicted by DFT-based calculations at 0 K in the literature, the predicted FE-PE transition temperatures are 776 and 653 K, respectively, in comparison with experimental 763 K in the literature. The present author's group is currently computing the free energies of the three configurations, aiming for a more accurate prediction. Our preliminary results on the energy-volume curves at 0 K are plotted in Figure 3, showing that equilibrium volumes at 0 K are 603.79, 600.50, and 597.62 Å3> for FEG, 90 DW, and 180 DW, respectively, in agreement with 603.42 and 599.88 Å3 for FEG and 180 DW configurations reported in the literature by the DFT-based calculations. Based on the zentropy theory, the NTE in PbTiO3 originates from the fact that the volume of the FEG ground-state configuration is larger than those of 90 DW and 180 DW symmetry-breaking non-ground-state configurations.
Figure 3. Predicted energy-volume curves of PbTiO3 at 0 K for FEG, 90 DW, and 180 DW configurations, respectively. The filled circles are from DFT-based calculations using the LDA exchange-correlation functional, while the curves were fitted using the third-order Birch Murnaghan equation of states (EOS). The open circles represent the energy minimum from EOS fitting.
The zentropy theory postulates that the entropy of a system contains contributions from entropies of the ground-state and symmetry-breaking non-ground-state configurations of a system and the statistical configurational entropy among these configurations. With the free energies of individual configurations predicted from DFT, the zentropy theory integrates quantum and statistical thermodynamics with the partition function of each configuration calculated from its free energy instead of the total energy commonly used in the literature. With the accurate free energy landscape of the system predicted by the zentropy theory, the properties of the system can be predicted by the first and higher-order derivatives of free energy with respect to its natural variables such as volume as the first derivative and thermal expansion as the second derivative. While the derivatives between conjugate variables are always positive for a stable system, the derivatives between non-conjugate variables, such as thermal expansion, can be either positive or negative, with the latter due to the larger volume of the ground-state configuration than those of non-ground-state configurations in the system. It is articulated that the emergent behaviors and anharmonicity originate from the competition among the configurations, and their accurate predictions can be realized by the zentropy theory.
Writting: Liu ZK
Revised; figure 3: Hew NLE, Shang SL
Availability of data and materials
Financial support and sponsorship
The present review article covers research outcomes supported by multiple funding agencies over multiple years, with the most recent ones including the Endowed Dorothy Pate Enright Professorship at the Pennsylvania State University, U.S. Department of Energy (DOE) Grant No. DE-SC0023185, No. AR0001435, No. DE-NE0008945, and No. DE-NE0009288, along with U.S. National Science Foundation (NSF) Grant No. NSF-2229690
Conflicts of interest
All authors declared that there are no conflicts of interest.
Ethical approval and consent to participate
Consent for publication
© The Author(s) 2024.
2. Wittenauer J. The invar effect: a centennial symposium. Mine Met Mater Soc 1997.
3. van Schilfgaarde M, Abrikosov IA, Johansson B. Origin of the invar effect in iron-nickel alloys. Nature 1999;400:46-9.
4. Miller W, Smith CW, Mackenzie DS, Evans KE. Negative thermal expansion: a review. J Mater Sci 2009;44:5441-51.
6. Mittal R, Gupta M, Chaplot S. Phonons and anomalous thermal expansion behaviour in crystalline solids. Prog Mater Sci 2018;92:360-445.
7. Liang E, Sun Q, Yuan H, Wang J, Zeng G, Gao Q. Negative thermal expansion: Mechanisms and materials. Front Phys 2021;16:53302.
8. Wendt D, Bozin E, Neuefeind J, et al. Entropic elasticity and negative thermal expansion in a simple cubic crystal. Sci Adv 2019;5:eaay2748.
9. Lohaus SH, Heine M, Guzman P, et al. A thermodynamic explanation of the Invar effect. Nat Phys 2023;19:1642-8.
10. Liu Z, Wang Y, Shang S. Zentropy theory for positive and negative thermal expansion. J Phase Equilib Diffus 2022;43:598-605.
11. Landau LD, Lifshitz EM. Statistical physics Oxford: Pergamon Press Ltd.; 1970.
13. Kohn W, Sham LJ. Self-consistent equations including exchange and correlation effects. Phys Rev 1965;140:A1133-8.
14. Liu Z. Theory of cross phenomena and their coefficients beyond Onsager theorem. Mater Res Lett 2022;10:393-439.
15. Liu Z. Thermodynamics and its prediction and CALPHAD modeling: review, state of the art, and perspectives. Calphad 2023;82:102580.
16. Liu Z, Li B, Lin H. Multiscale entropy and its implications to critical phenomena, emergent behaviors, and information. J Phase Equilib Diffus 2019;40:508-21.
18. Wang Y, Liu Z, Chen L. Thermodynamic properties of Al, Ni, NiAl, and Ni3Al from first-principles calculations. Acta Mater 2004;52:2665-71.
20. Liu Z, Shang S, Du J, Wang Y. Parameter-free prediction of phase transition in PbTiO3 through combination of quantum mechanics and statistical mechanics. Scr Mater 2023;232:115480.
Cite This Article
Liu ZK, Hew NLE, Shang SL. Zentropy theory for accurate prediction of free energy, volume, and thermal expansion without fitting parameters. Microstructures 2024;4:2024009. http://dx.doi.org/10.20517/microstructures.2023.56
Liu ZK, Hew NLE, Shang SL. Zentropy theory for accurate prediction of free energy, volume, and thermal expansion without fitting parameters. Microstructures. 2024; 4(1): 2024009. http://dx.doi.org/10.20517/microstructures.2023.56
Liu, Zi-Kui, Nigel L. E. Hew, Shun-Li Shang. 2024. "Zentropy theory for accurate prediction of free energy, volume, and thermal expansion without fitting parameters" Microstructures. 4, no.1: 2024009. http://dx.doi.org/10.20517/microstructures.2023.56
Liu, Z.K.; Hew NLE.; Shang S.L. Zentropy theory for accurate prediction of free energy, volume, and thermal expansion without fitting parameters. Microstructures. 2024, 4, 2024009. http://dx.doi.org/10.20517/microstructures.2023.56
About This Article
Data & Comments
Comments must be written in English. Spam, offensive content, impersonation, and private information will not be permitted. If any comment is reported and identified as inappropriate content by OAE staff, the comment will be removed without notice. If you have any queries or need any help, please contact us at email@example.com.