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 non-linear emergent behaviors. It is shown that the derivative of the volume to temperature of the single phase can be negative, i.e., negative thermal expansion (NTE), 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 NTE 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.


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.Though 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 [1] with near zero thermal expansion (ZTE) at room temperature, commonly referred as INVAR.
While there have been many theories that aim to describe NTE mechanisms in various materials, they are mostly phenomenological aiming for interpretation of experimental observations [2][3][4][5][6][7][8][9], and a fundamental understanding and a predictive theory without experimental inputs applicable to all materials are still lacking.Based on Maxwell relation, the derivative of V to T under constant pressure equals to 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 and the system's partition functions as follows where  ! and  ! are the probability and partition function of configuration , and is the partition function of the system or the phase with  being the number of independent configurations.The configurational entropy among the configuration is obtained as where  % is the Boltzmann constant.
For a closed system under hydrostatic pressure (), the combined law of thermodynamics in terms of internal/total energy, , is written as 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  ! is the total energy of the configuration .The Helmholtz energy of the phase can be obtained from above equation as follows When there is only one configuration in the system, Eq. 6 becomes  =  !Eq. 7 Since  =  !−  ! by definition, Eq. 7 gives  != 0 at finite temperature, indicating that the configurations are all pure quantum states without any unspecified internal degrees of freedom as implicitly implied by Gibbs as the quantum mechanics was not invented yet at that time and envisioned by Landau and Lifshitz [11].For systems of practical interest, the number of pure quantum states is very large, and their complete sampling is in general intractable.The current 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] 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 forest, planet, 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  until it reaches the ground-state configuration with its properties predicted by the density functional theory (DFT) [16].The latter may provide some insights into superconducting and other interesting ground-state configurations as postulated by the present author [14,15].
All with the same NVT Nested formula è Down to quantum scale

Zentropy Top Down
Figure 1: Schematic top-down and bottom-up integration of the zentropy theory [15].
The Helmholtz energy of the system can thus be obtained as Eq. 9 where  !=  !−  ! is the Helmholtz energy of configuration .Re-arranging Eq. 9 in the form of partition function, one obtains Eq. 10 Eq. 11 Eq. 8 to Eq. 11 reduce to standard statistical mechanics when  != 0, i.e., pure quantum configurations with  !=  ! as discussed in Section 2.
Formulated as an exact theory of many-body systems, DFT [12] articulates that for an interacting electron gas there exists a universal functional 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 [18] added the vibrational contribution and presented the Helmholtz energy as follows  !=  !,+ +  !,,-+  !,./0 =  !−  !Eq. 12  !=  !,+ +  !,,-+  !,./0 Eq. 13 !=  !,,-+  !,./0 Eq. 14 where  !,,-,  !,,-, and  !,,-are the contributions of thermal electron to Helmholtz energy, total energy, and entropy of configuration  based on the Fermi-Dirac statistics for electrons, and  !,./0 ,  !,./0 , and  !,./0 are the vibrational contributions to Helmholtz energy, total energy, and entropy of configuration  based on the Bose-Einstein statistics for phonons, respectively.
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-rang ordering, and defects such as vacancy, dislocation, and stacking faults.For their  ! to be predicted from DFT as a function of external stimuli, the non-ground-state configurations must be metastable.With their  !calculated from partition functions using  !, the zentropy theory enables the integration of the quantum and statistical mechanics through Eq. 8 to Eq. 14 and is capable to predicting how a system responds macroscopically to external stimuli.

Derivative of volume to temperature
The pressure of a system can be calculated from the derivative of Helmholtz energy to volume as follows where  ! is the pressure of configuration  evaluated at the system volume.In our previous publications as reviewed in refs.[14,15], the volume for 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 singlephase 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 the free energy barrier 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 , , and , in order to evaluate the derivative of volume to temperature under constant pressure.The combined law of thermodynamics is thus written in terms of Gibbs energy as follows,  = − +  Eq. 16 From the Maxwell relation, one has   = −   Eq. 17 The statistical mechanics in terms of Gibbs ensemble is shown below Eq. 18 where  ! is the Gibbs energy of configuration .The volume and the derivative of volume to temperature is obtained as follows Eq. 20 where  2 is the volume of the ground-state configuration.From Eq. 20, it can be seen that if With increase of temperature, the probability of the groundstate configuration decreases, i.e., 34 # 3( < 0, while the probabilities of non-ground-state configurations increases, i.e., The first term Eq.21 is a linear combination of contributions from each configuration, and the second term gives nonlinear behavior.The condition for 35 3( = 0, i.e., ZTE, can be obtained as follows As all variables in Eq. 22 are positive, and no contributions from Eq. 23 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 separate 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.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 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 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 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.,   > 0 Eq.24 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 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 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.,

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 total energy commonly used in the literature.With 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.

34 # 3 (
, the necessary condition to have a solution is for  !<  2 , and the sufficient condition is that the volume decrease due to the second summation surpasses the weighted sum of the volume increase of individual configurations as shown by the first summation in the equation, i.e.,

Figure 2 :
Figure 2: Temperature-volume phase diagrams with isobaric volumes at various pressures of (a) 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 Figure2(b) with significant anharmonic behaviors.More recently, the zentropy theory was applied to predict the ferroelectric-paraelectric (FE-PE) transitions in PbTiO3[20] with three configurations considered, i.e., tetragonal polarized configurations without domain wall (FEG), with 90° domain wall (90DW), and with 180° domain wall (180DW).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 more accurate prediction.Our preliminary results on the energy-volume curves at 0 K are plotted in Figure3, showing that equilibrium volumes at 0K are 603.79,600.50, and 597.62 Å 3 for FEG, 90DW, and 180DW, respectively, in agreement with 603.42 and 599.88 Å 3 for FEG and 180DW configurations reported in the literature by the DFT-based calculations[21].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 90DW and 180DW symmetry-breaking nonground-state configurations.