materials
Article Fundamentals of Thermal Expansion and Thermal Contraction
Zi-Kui Liu *, Shun-Li Shang and Yi Wang Department of Materials Science and Engineering, The Pennsylvania State University, University Park, PA 16802, USA; [email protected] (S.-L.S.); [email protected] (Y.W.) * Correspondence: [email protected]
Academic Editor: Thomas Fiedler Received: 29 March 2017; Accepted: 11 April 2017; Published: 14 April 2017
Abstract: Thermal expansion is an important property of substances. Its theoretical prediction has been challenging, particularly in cases the volume decreases with temperature, i.e., thermal contraction or negative thermal expansion at high temperatures. In this paper, a new theory recently developed by the authors has been reviewed and further examined in the framework of fundamental thermodynamics and statistical mechanics. Its applications to cerium with colossal thermal expansion and Fe3Pt with thermal contraction in certain temperature ranges are discussed. It is anticipated that this theory is not limited to volume only and can be used to predict a wide range of properties at finite temperatures.
Keywords: negative thermal expansion (NTE); thermal contraction; theory; thermodynamics; statistical mechanics
1. Introduction It is well known that thermal expansion originates from the effect of anharmonic terms in the potential energy on the mean separation of atoms at a temperature [1], as schematically shown in Figure1a in the vicinity of equilibrium separation distance. When the temperature is increased, the kinetic energy of atoms increases, and the atoms vibrate and move, resulting in a greater average separation of atoms and thus thermal expansion, i.e., the vibrational origin of thermal expansion. It is therefore self-evident that for a substance with a constant separation of atoms at a temperature, i.e., zero thermal expansion, its potential energy must be symmetric with respect to the separation (see Figure1b). Furthermore, if the separation of the substance at a temperature decreases with the increase in temperature, i.e., thermal contraction or negative thermal expansion, its potential energy must be again asymmetric as shown in Figure1c, but in the opposite direction of thermal expansion. The state with zero thermal expansion is thus the boundary between the thermal expansion and thermal contraction. Thermal contraction is counterintuitive as the increase in kinetic energy with temperature should push atoms apart rather then pull them together, but in reality many substances experience thermal contraction in certain temperature and pressure ranges [2]. There are a number of excellent reviews on this topic in the literature [3–7]. Many theoretical interpretations for thermal contraction have been developed for individual or groups of substances, but there is a lack of fundamental understanding how the potential energy of a substance can change its asymmetry into the opposite direction with increasing temperature from Figure1a, to Figure1b, to Figure1c, and then back to Figure1a,b. Consequently, thermal contraction could not be predicted with the existing theories and models, particularly at high temperatures. In our previous works [2,8–10], we proposed that a single-phase substance at high temperatures may consist of many states, each with its own potential energy similar to Figure1a. All these states, except one and its degeneracies, are metastable in terms of total energy. When the statistical probabilities of the metastable states
Materials 2017, 10, 410; doi:10.3390/ma10040410 www.mdpi.com/journal/materials Materials 2017, 10, 410 2 of 11 Materials 2017, 10, 410 2 of 10 increaseincrease withwith temperature, and and their their physical physical proper properties,ties, such such as as volume, volume, may may have have values values larger larger or orsmaller smaller than than those those of ofthe the stable stable state, state, the the behaviors behaviors of of the the substance substance can can be significantlysignificantly altered, resultingresulting inin anomaliesanomalies suchsuch asas thermalthermal contraction.contraction. Furthermore,Furthermore, ourour computationalcomputational approachapproach basedbased onon first-principlesfirst-principles calculationscalculations andand statisticalstatistical mechanicsmechanics was able toto predictpredict thethe criticalcritical pointspoints ofof ceriumcerium andand Fe Fe3Pt3Pt and and the the associated associated maximum maximum anomalies anomalies of thermal of thermal expansion expansion in cerium in andcerium thermal and contractionthermal contraction in Fe3Pt, in i.e., Fe positive3Pt, i.e., andpositive negative and negative infinite at infinite their critical at their points, critical respectively, points, respectively, in addition in toaddition thermal to expansion thermal andexpansion contraction and temperaturecontraction te rangesmperature in their ranges single-phase in their regions single-phase under variousregions pressuresunder various [10]. pressures [10].
FigureFigure 1.1. Schematic potential energy of aa substancesubstance withwith ((aa)) thermalthermal expansion;expansion; ((bb)) zerozero thermalthermal expansion;expansion; ((cc)) thermalthermal contractioncontraction oror negativenegative thermalthermal expansion.expansion.
In the present paper, the fundamentals of thermal expansion and thermal contraction are In the present paper, the fundamentals of thermal expansion and thermal contraction are discussed discussed in terms of the asymmetry of potential energy. Our previously proposed theory is in terms of the asymmetry of potential energy. Our previously proposed theory is reviewed and further reviewed and further analyzed based on these fundamentals, and thus becomes more general and analyzed based on these fundamentals, and thus becomes more general and comprehensive. comprehensive. 2. Thermodynamics of a Substance 2. Thermodynamics of a Substance For the sake of simplicity, let us consider a closed system at equilibrium with its combined law of thermodynamicsFor the sake writtenof simplicity, as let us consider a closed system at equilibrium with its combined law of thermodynamics written as dG = −SdT − Vd(−P) (1) = − − − (1) where G, S, and V are the Gibbs energy, entropy, and volume of the system, and T and P the temperature andwhere pressure, G, S, respectively.and V are the It isGibbs noted energy, that S and entropy,T are conjugate and volume variables, of the so aresystem,−P and andV .T The and stability P the oftemperature the system and requires pressure, that therespectively. conjugate It variables is noted changethat S and in the T are same conjugate direction variables, as follows so are [11, 12−P]: and V. The stability of the system requires that the conjugate variables change in the same direction as ∂T follows [11,12]: > 0 (2) ∂S V or P 0 (2) ∂(− P) > 0 (3) ∂ V− S or T 0 (3)
Materials 2017, 10, 410 3 of 11
When both derivatives become zero, the system reaches the critical point and becomes unstable, resulting in the maximum anomalies of entropy and volume changes as shown below:
∂S ∂V = = +∞ (4) ∂T V or P ∂(−P) S or T
A positive sign is put in front of infinity because the derivatives are positive before they approach infinity as shown by Equations (2) and (3). For an unstable system, one or more of its vibrational frequencies become imaginary, commonly referred as soft vibrational modes [13]. It is evident that the derivative of volume to temperature under constant pressure, i.e., isobaric thermal expansion, at the critical point is also infinite, but there are no fundamental relations indicating whether the derivative is positive or negative because volume and temperature are not conjugate variables. What is known is the Maxwell relation derived from the second derivatives of Gibbs energy from Equation (1) as follows [11,12]:
∂V ∂2G ∂2G ∂S = = = (5) ∂T P ∂T∂(−P) ∂(−P)∂T ∂(−P) T
It can be seen from Equation (5) that the thermal expansion is related to entropy change with respect to pressure. Even though entropy in a stable system increases with its conjugate variable, temperature, as shown by Equation (2), its relation with its non-conjugate variable, pressure, is not defined by any fundamental laws either, the same as the thermal expansion. In Figure1, the negative pressure, −P, is represented by the slope of the potential energy curve. In Figure1a, the higher entropy at higher temperature corresponds to positive values of ( −P) with larger separations of atoms, and S and (−P) thus vary in the same direction, resulting in positive values for Equation (5), i.e., thermal expansion. While in Figure1c, the higher entropy at higher temperature corresponds to negative (−P) with smaller separations of atoms, and S and (−P) thus vary in the opposite directions, resulting in negative values for Equation (5), i.e., thermal contraction. In the next section, the transitions between potential energies shown in Figure1 are discussed.
3. Potential Energy on Mean Separation of Atoms/Molecules For any given states of a substance, their potential energies at 0 K can be calculated from first-principles calculations based on the density functional theory (DFT) [14,15]. All exhibit shapes similar to the one shown in Figure1a with the ground state having the lowest minimum potential energy for a classical oscillator [1], and other metastable states having higher minimum potential energies and larger or smaller equilibrium atomic/molecular distances and volumes. It is evident that the potential energy curve at the vicinities of the minimum potential energy in Figure1b,c may be mathematically represented by the weighted average of potential energies of the ground state and those metastable states as schematically shown in Figure2. The combined potential energy, Ec, can be written as follows: Ec = Eg + ∑ pi Ei − Eg (6) where Ec and Eg are the potential energies of the combined and ground/stable states, respectively, and pi and Ei the weight and potential energy of metastable state i. In Figure2a, the equilibrium separation distance of the metastable state is smaller than that of the ground state, and their potential energy curves cross each other at rig. The combination of the metastable state with the ground state thus increases the energy of the ground state at the separation distances larger than rig, and decreases the energy at the separation distances smaller than rig. The equilibrium separation distance of the combined state is thus smaller than that of the ground state though with a higher potential energy. In case the potential energy curves of two states do not cross each other, the addition of the metastable state increases the energy of the ground state more at the larger separation distance than at the smaller separation distance, and still decreases the Materials 2017, 10, 410 4 of 11
Materials 2017, 10, 410 4 of 10 equilibrium separation distance of the combined state. On the other hand, Figure2b shows the case thatequilibrium the equilibrium separation separation distance distance of the metastable of the metastable state is state larger is largerthan that than of that the of ground the ground state. state. The Thecombined combined state state thus thus has hasa larger a larger equilibrium equilibrium separation separation distance distance than than the theground ground state. state.
FigureFigure 2.2. SchematicSchematic diagrams diagrams depicting depicting the the combination combination of potential energies ofof variousvarious configurations/statesconfigurations/states of of a substance to show (a) thermalthermal contractioncontraction oror negativenegative thermalthermal expansion;expansion; andand ((bb)) thermalthermal expansion.expansion.
Mathematically,Mathematically, it it is is then then possible possible to findto find a set a of se metastablet of metastable states sostates that theso combinedthat the combined potential energypotential represented energy represented by Equation by 6 becomes Equation quadratic 6 become (Figures quadratic1b) first and(Figure then 1b) tilts first to smaller and then separation tilts to distancesmaller separation (Figure1c) distance at the vicinity (Figure of 1c) its at new the equilibriumvicinity of its separation new equilibrium distance. separation The potential distance. energy The thuspotential becomes energy symmetric thus becomes with zero symmetric thermal expansionwith zero orthermal exhibits expansion anomaly or with exhibits thermal anomaly contraction. with Inthermal the next contraction. section, energetics In the next of groundsection, andenergetics metastable of ground states and theirmetastable weights states at finite and temperaturestheir weights willat finite be discussed temperatures and evaluated will be discussed in terms ofand DFT-based evaluated first-principles in terms of calculationsDFT-based first-principles and statistical mechanics.calculations Instead and statistical of the separation mechanics. distance Instead used of the in this separation section fordistance the sake used of simplicity,in this section the volume for the issake used of insimplicity, the rest of the the volume paper. is used in the rest of the paper. ItIt shouldshould bebe mentionedmentioned thatthat BarreraBarrera etet al.al. [[4]4] discusseddiscussed thethe conceptconcept of “apparent” and “true” separationseparation distancesdistances for for the the distance distance between between the meanthe mean positions positions of two of atoms two andatoms the and mean the distance mean betweendistance twobetween atoms, two respectively. atoms, respectively. The “apparent” The separation “apparent” distance separation is obtained distance from is X-ray obtained or neutron from diffractionX-ray or neutron measurements diffraction with measurements low temporal with resolution, low temporal and theresolution, “true” separation and the “true” distances separation from extendeddistances X-rayfrom extended absorption X-ray fine structureabsorption (XAFS) fine stru withcture high (XAFS) temporal with resolution. high temporal One mayresolution. attempt One to connectmay attempt the “apparent” to connect and the “true” “apparent” separation and distances “true” with separation the separation distances distances with ofthe the separation combined statedistances and theof the individual combined states state in an thed presentthe individual work, respectively.states in the Thispresent was work, demonstrated respectively. in our This recent was demonstrated in our recent ab initio molecular dynamic simulations of PbTiO3 in comparison with ab initio molecular dynamic simulations of PbTiO3 in comparison with the experimental X-ray and XAFSthe experimental data in the literatureX-ray and showing XAFS data the apparentin the literature cubic structure showing by the X-ray apparent and neutron cubic structure diffraction by measurementsX-ray and neutron and thediffraction significant measurements amount of “true” and the local significant tetragonal amount structures of “true” by XAFS local [16 tetragonal]. structures by XAFS [16]. 4. Statistical Mechanics of Ground and Metastable States and Their Interactions 4. Statistical Mechanics of Ground and Metastable States and Their Interactions At a given temperature and pressure (zero external pressure used herein for simplicity), the equilibriumAt a given volume temperature of a state i ,andVi, of pressure a substance (zero is determinedexternal pressure when itsused Helmholtz herein energy,for simplicity), defined the equilibrium volume of a state i, , of a substance is determined when its Helmholtz energy, defined as = − , is minimized, where the entropy of the state, , is typically obtained from quasi-harmonic first-principles phonon calculations or Debye models including both vibrational and
Materials 2017, 10, 410 5 of 11
asMaterialsFi = 2017Ei, −10, TS410i , is minimized, where the entropy of the state, Si, is typically obtained from5 of 10 quasi-harmonic first-principles phonon calculations or Debye models including both vibrational andthermal thermal electronic electronic contributions contributions [17–20]. [17 –It20 should]. It should be noted be notedthat the that phonon the phonon calculations calculations are usually are usuallybased on based the potential on the potentialenergy at energy0 K. In atthe 0 case K. Ins where the cases the vibrational where the energy vibrational of the energy ground of state the groundreaches statethe potential reaches energy the potential of metastable energy states, of metastable the calculated states, vibrationa the calculatedl energy vibrational of the substance energy ofcontains the substance thus vibrational contains thus contributions vibrational contributionsfrom both the from ground both thestate ground and statethe metastable and the metastable state at statetemperatures at temperatures close to close 0 K, to representing 0 K, representing the thequan quantumtum origin origin of ofthe the vibrational vibrational energy energy for thermalthermal expansion.expansion. WithWith thethe metastablemetastable statestate havinghaving smallersmaller volumesvolumes thanthan thethe groundground state,state, asas schematicallyschematically shownshown inin FigureFigure3 ,3, the the substance substance can can show show thermal thermal contraction contraction as demonstratedas demonstrated for for ice ice and and Si close Si close to 0to K 0 in K ourin our previous previous publication publication [10 ][10] and and in the in the literature literature [21], [21], representing representing the quantumthe quantum origin origin of the of vibrationalthe vibrational energy energy for thermal for thermal contraction. contraction. The rest The of the rest paper of the addresses paper addresses cases when cases the vibrational when the energyvibrational of the energy ground of the state ground is smaller state thanis smaller the energy than the difference energy difference between thebetween ground the state ground and state the metastableand the metastable states, i.e., states, typical i.e., phonon typical calculationsphonon calculations do not show do not property show property anomalies. anomalies.
FigureFigure 3.3. SchematicSchematic diagramdiagram showingshowing thethe groundground statestate andand aa metastablemetastable statestate withwith thethe vibrationalvibrational energyenergy ofof thethe groundground statestate reachingreaching the the potential potential energy energy of of the the metastable metastable state. state.
InIn thethe canonicalcanonical ensembleensemble withwith NN,, VV,, andand TT asas thermodynamicthermodynamic naturalnatural variables,variables, thethe canonicalcanonical partitionpartition functionfunction ofof aa statestate isis writtenwritten asas = − F i = − Eik (7) kBT e kBT Zi = e = ∑ k (7) where are the energy eigenvalues of microstate k in the state i, including the contributions from wherephononsEik andare thermal the energy electrons eigenvalues [12]. Fo ofr a microstate system withk inmultiple the state statesi, including schematically the contributions shown in Figure from 2, phononsthe canonical and thermal partition electrons function [12 of]. the For system a system is as with follows multiple [8,9,22]: states schematically shown in Figure2, the canonical partition function of the system is as follows [8,9,22]: = = E (8) Fc cj − k T e − k T Zc = e B = ∑ j B (8) where is the Helmholtz energy of the system with multiple states, and the energy whereeigenvaluesFc is the of Helmholtz microstate energy j in the of thecombined system withstate, multiple c, including states, contributions and Ecj the energy from eigenvaluesphonons and of microstatethermal electronsj in the combined [9,12]. The state, contric, includingbution of contributions each state to from the phononscombined and state, thermal i.e., electronsthe weights [9,12 in]. TheEquation contribution 6, can be of defined each state by toth thee partition combined functions state, i.e., as follows: the weights in Equation 6, can be defined by the partition functions as follows: Z == i . (9) pi . (9) Zc CombiningCombining EquationsEquations (8)(8) andand (9)(9) andand rearrangingrearranging them,them, oneone obtainsobtains [[8–10,22]8–10,22] =− =− Zi − Zi Zi = − (10) Fc = −kBTlnZc = −kBT ∑ lnZi − ∑ ln = ∑ pi Fi − TSSCE (10) Zc Zc Zc with =− ∑ defined as the state configurational entropy [10]. The total entropy of the system of the combined state can then be written as