# Zentropy theory for accurate prediction of free energy, volume, and thermal expansion without fitting parameters

*Microstructures*2024;4:2024009.

## Abstract

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.

## Keywords

*,*entropy

*,*density functional theory

*,*statistical mechanics

*,*anharmonicity

*,*emergent behaviors

*,*negative thermal expansion

## INTRODUCTION

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 Fe_{65}Ni_{35} alloy in 1897 by Guillaume^{[1]} 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^{[10]}, 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

where

where

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

where

When there is only one configuration in the system, Eq. 6 becomes

Since ^{[11]}. For systems of practical interest, the number of pure quantum states is very large, and their complete sampling is, in general, intractable. The currently available solution is their coarse-graining through density functional theory (DFT)^{[12,13]}, resulting in a non-zero entropy for each configuration at finite temperature and, thus, the necessity to modify the formula of entropy and partition function in terms of the zentropy theory as discussed below.

## 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]}

Figure 1. Schematic top-down and bottom-up integration of the zentropy theory^{[15]}. Reproduced with the permission of Ref.^{[15]} Copyright © 2023, Elsevier.

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 ^{[16]}. The latter may provide some insights into superconducting and other interesting ground-state configurations, as postulated by the present author^{[14,15]}.

The Helmholtz energy of the system can thus be obtained as

where

Eq. 8 to Eq. 11 reduce to standard statistical mechanics when

Formulated as an exact theory of many-body systems, DFT^{[12]} 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^{[13]}, i.e., *coarse-graining of electrons*. Kohn and Sham^{[13]} used the finite temperature generalization of ground-state energy of an interacting inhomogeneous electron gas by Mermin^{[17]} and formulated the entropy of thermal electrons at finite temperatures. Wang *et al*. added the vibrational contribution and presented the Helmholtz energy as follows^{[18]}

where

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

where ^{[14,15]}, the volume for the given temperature and pressure was numerically evaluated from Eq. 15, along with the Helmholtz energy of the system. When the Helmholtz energy of the system is at its lowest value with one well, the system is in a single-phase region. When Helmholtz energy of the system with double or more wells can be lowered by separating into two phases with different volumes, the system is in a two or more-phase region. This indicates that the zentropy theory can predict the Helmholtz energy of the system under metastable and unstable states; thus, it can also estimate the free energy barrier existing between stable and metastable states. This is significant because the common wisdom is that the free energy of an unstable state could not be defined due to the imaginary vibrational modes that prevent the evaluation of its entropy. However, this view assumes that the atoms are static when evaluating the phonon properties, while the atoms are constantly moving at finite temperature. As all configurations used in the zentropy theory are stable, it does not depend on this unrealistic assumption in evaluating the free energy of unstable states.

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

where

where

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 Fe_{3}Pt, 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 Fe_{3}Pt, 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) Fe_{3}Pt. The volume (V) is normalized to their respective equilibrium volume (_{3}Pt. Below the critical point marked by the green open circle, the single phase is no longer stable and decomposes into a two-phase mixture in the region of a miscibility gap. Symbols are from various experimental measurements in the literature^{[19]}. Reproduced with the permission of Ref.^{[19]} Copyright © 2014, The Authors.

## ANHARMONICITY AND EMERGENT BEHAVIORS

Anharmonicity is usually represented by the deviation of entropy or heat capacity away from quasiharmonic behavior^{[5]}. It is noted in Figure 2 that the volume change for Fe_{3}Pt, 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 PbTiO_{3}^{[20]} with three configurations considered, i.e., tetragonal polarized configurations without domain wall (FEG), with 90^{o} domain wall (90 DW), and with 180^{o} 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^{[20]}. 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^{[21]}. Based on the zentropy theory, the NTE in PbTiO_{3} 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 PbTiO_{3} 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.

## SUMMARY

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.

## DECLARATIONS

### Availability of data and materials

Not applicable.

### 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

Not applicable.

### Consent for publication

Not applicable.

### Copyright

© The Author(s) 2024.

## REFERENCES

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 Ni_{3}Al 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 PbTiO_{3} through combination of quantum mechanics and statistical mechanics. *Scr Mater* 2023;232:115480.

## Cite This Article

## How to Cite

Liu, Z. K.; Hew N. L. E.; 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

## Download Citation

## Export Citation File:

## Type of Import

### Tips on Downloading Citation

### Citation Manager File Format

### Type of Import

**Direct Import:**When the Direct Import option is selected (the default state), a dialogue box will give you the option to Save or Open the downloaded citation data. Choosing Open will either launch your citation manager or give you a choice of applications with which to use the metadata. The Save option saves the file locally for later use.

**Indirect Import:**When the Indirect Import option is selected, the metadata is displayed and may be copied and pasted as needed.

## About This Article

### Special Issue

### Copyright

**Open Access**This article is licensed under a Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, sharing, adaptation, distribution and reproduction in any medium or format, for any purpose, even commercially, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

## Data & Comments

### 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 support@oaepublish.com.

^{0}