Equation of state for solids with considering cohesive energy and anharmonic effect and its application to MgO

A simple equation of state (EOS) in wide ranges of pressure and temperature is constructed within the Mie-Grüneisen-Debye framework. Instead of the popular Birch-Murnaghan and Vinet EOS, we employ a five-parameter cold energy expression to represent the static EOS term, which can correctly produce cohesive energy without any spurious oscillations in extreme compression and expansion region. We developed a Padé approximation-based analytic Debye quasiharmonic model with high accuracy which improves the performance of EOS in low temperature region. The anharmonic effect is taken into account by using a semi-empirical approach. Its reasonability is verified by the fact that the total thermal pressure tends to the lowest-order anharmonic expansion in the literature at low temperature, and tends to ideal-gas limitation at high temperature, which is physically correct. Besides, based on this approach, the anharmonic thermal pressure can be expressed in the Grüneisen form, which is convenient for applications. The proposed EOS is used to study the thermodynamic properties of MgO including static and shock compression conditions, and the results are very satisfactory as compared with the experimental data.     

High-temperature thermodynamics of silver:Semi-empirical approach

We report high-temperature thermodynamics for fcc silver by combining ab initio phonon dynamics to empirical quadratic temperature-dependent term for anharmonic part of Helmholtz free energy. The electronic free energy is added through an interpolation scheme, which connects ambient condition free electron gas model to Thomas-Fermi results.The present study shows good agreement with experimental and reported findings for several thermal properties, and the discrepancy observed in some caloric properties is addressed. The decreases in the product of volume thermal expansion coefficient and isothermal bulk modulus and in the constant volume anharmonic lattice specific heat at high temperature are the clear evidences of proper account of anharmonicity. The present study also reveals that T~2-dependent anharmonic free energy is sufficient for correct evaluation of thermal pressure and conventional Grüneisen parameter. We observe that the intrinsic phonon anharmonicity starts dominating above characteristic temperature, which is attributed to higher order anharmonicity and can be related to higher order potential derivatives. We conclude that the uncorrelated and largeamplitude lattice vibrations at high temperature raise dominating intrinsic thermal stress mechanism, which surpasses the phonon-anharmonism and requires future consideration.     

Quantum corrected cell model for an anharmonic generalized Lennard-Jones solid

A simple method is proposed to dispose the quantum effect and anharmonic effect at the same time. Considering the quantum effect is remarkable only at low temperature, and tends to zero at high temperature, the potential energy of an atom is expanded harmonically to consider the quantum effect of solids within the harmonic oscillator framework. The anharmonic effect is remarkable only at high temperature, and tends to zero at low temperature, it was disposed by using a classical approximation. The universal formalism is applied to the generalized Lennard-Jones solid. The comparison shows that the results with and without anharmonic effect are in agreement with each other at some low temperature, to which the Einstein model is applicable. The results without anharmonic effect become divergent at slightly higher temperatures; however, the results including anharmonic effect are in good agreement with the experimental data of solid xenon. The method proposed in this paper can be extended to other potentials to develop practical molecular thermodynamic equations of state for solids.     

Equations of State for Regular Solids

Equations of State (EOS) for solids under strong compression and wide ranges in temperature are most commonly represented by "parametric" EOS forms using temperature dependent parameters for the volume V 0 , the bulk modulus K 0 , and its pressure derivative K_{0}^{\prime} for the given ambient (or zero) pressure. Therefore, various common "parametric" EOS forms are compared at first with the recently proposed [1,2] A dapted P olynomical expansion AP2, and in the second part with the Mie-Grüneisen approach, which uses one of the common EOS forms for the pressure of the static lattice or for the zero temperature isotherm, p ZT ( V ), and a detailed modelling of the additional thermal pressure, p th ( V , T ), in the form p(V,T) = p_{\rm ZT}(V) + p_{\rm th}(V,T) . Thereby, it is shown, that "intrinsic" anharmonicity effects have to be taken into account and between the two schemes differences of the order of a few percent in pressure are noticed for regular solids, like Cu, Ag, and Au. These differences are discussed with respect to the present uncertainties in a practical pressure scale for wide ranges in pressure (up to several TPa) and in temperature (up to 1500 K and above).     

一个等温状态方程（Ⅴ）——不同温度的高压状态方程

 本文以NaCl为例，证实了在我们的状态方程体系中，不同温度的状态方程可以认为是在p_x坐标系中不同起点开始的同一状态方程。于是，可由室温状态方程和常压热膨胀的测量数据估计高温状态方程。     

First principles studies of elastic and thermodynamic properties of fcc-VH2 with pressure and temperature

In this work, we study elastic and thermodynamic properties of VH₂ at different pressures and temperatures. Elastic constants and bulk modulus of VH₂ decrease with increase in temperature, and hence increase with pressure. Thermal expansion of the crystal lattice will be suppressed by high pressure. When the temperature is 1500 K, 15.99 GPa of pressure can completely restrain the volume expansion caused by temperature. At a given pressure, the lower the temperature, the easier the cell compression. At low temperatures, C_v is proportional to T³, and C_v tends to the Dulong-Petit limit at higher temperatures. The Debye temperature increases with pressure, but decreases with temperature. At low temperature and low pressure, thermal expansion coefficient increases sharply with temperature. At high temperature and high pressure, the increasing trend slows down.     

Thermal equation of state, and melting and thermoelastic properties of bcc tantalum from molecular dynamics

We performed molecular dynamics simulations with the extended Finnis-Sinclair (EFS) potential to investigate thermal equation of state (EOS), and melting and thermoelastic properties of tantalum. The agreement of the obtained thermal EOS with experiments at ambient conditions is reasonably good. The EFS potential with the two-phase method also reproduced very satisfyingly the high-pressure melting curve, excellently consistent with both the experiments of melting temperature at ambient pressure and shock melting at high pressure. From molecular dynamics simulations, we also obtained the thermoelastic properties of Ta for temperatures up to 3000 K at ambient pressure. Fully including anharmonic effects in molecular dynamics, our calculated elastic constants are in excellent agreement with experimental data. Shear modulus G decreases quickly with increasing temperature.     

硅纳米晶体薄膜热膨胀性质的分子动力学研究

基于Stillinger-Weber势对硅纳米晶体薄膜的热膨胀性质进行了分子动力学模拟。研究表明,硅纳米晶体薄膜表面层原子的二聚现象引起薄膜收缩,而原子之间的非和谐势能引起薄膜膨胀;在约400K以下的低温段,由于硅纳米晶体薄膜表面层原子发生二聚的原子数目随温度的升高而明显增多,而原子间非和谐势能很小,故此时二聚主导热膨胀性质,热膨胀系数为负;在高温段（约400K以上）,由于发生二聚的原子数目随温度升高不再显著地增加并渐趋于稳定,而原子间非和谐势能逐渐显著并主导热膨胀性质,故热膨胀系数为正。     

Thermodynamic properties of cubic boron nitride based on an analytic mean field approach

We present a detailed derivation of the analytic expressions for the equation of state (EOS) and internal energy of Morse model solids based on an analytic mean field potential (AMFP) method. The formalism is applied to cubic boron nitride (c-BN). One set of potential parameters are determined by fitting the experimental P-V-T data of c-BN up to 160 GPa at 295 K and 80 GPa in the range 500–900 K. Various physical quantities including the isothermals, thermal expansion, isochoric heat capacity, Helmholtz free energy and internal energy are calculated and analyzed. The theoretical results are consistent with the available experimental data and those calculated by others. These results presented in this paper verify that the AMFP method is a useful approach to consider the anharmonic effect at high temperature. Numerous reasonable predictions and the change trend of the properties for c-BN at extreme conditions have been given.     

Five-parameter equation of solids considering thermal effect which correctly incorporates cohesive energy

A five-parameter equation of state (EOS) is proposed, which correctly incorporates the cohesive energy data without physically incorrect oscillations in both extreme high pressure and expansion regions. Based on a modified Einstein model, the thermal effect is included in the proposed EOS complying with the zero-pressure condition. With this thermal EOS applied to five solids (Ar, Al, Au, Cu and Li), some important thermodynamic properties as isotherm, isochore, thermal expansion coefficient, volume modulus, heat capacity and Hugoniot are calculated for each selected solid with good agreement with experimental data, which confirms the validity of the present work     

Analysis of energy distribution in a loaded quantum anharmonic oscillator in a wide temperature range

Analysis of the energy distribution in an ensemble of quantum anharmonic oscillators loaded by an external force in a wide temperature range (from T = 0) is carried out using a general approach based on the virial theorem. At T = 0, anharmonic effects are observed: a linear variation of zero-point energy of an oscillator under loading (energy decrease during extension and increase under compression) and a linear variation of the average kinetic and potential energy components. At high temperatures, at which the dynamics of the oscillators becomes classical, the anharmonic effects are manifested in a linear variation in the vibrational energy and a linear variation in the average kinetic and potential energy components upon an increase in force. Mutually compensating variation in the average kinetic and potential energy components of the internal dynamic energy of an oscillator (energy redistribution upon loading) takes place both at low and high temperatures.     

极端条件下锆的动力学稳定性研究

过渡金属Zr具有优良的物理、化学及力学性能, 具有广泛的应用价值. 主要通过新近发展的自洽晶格动力学方法, 充分考虑声子间的相互作用, 成功获得了β-Zr的高温高压声子色散曲线, 预测了β-Zr在相图中能够稳定存在的区域, 进一步比较α-Zr, ω-Zr和β-Zr的自由能, 获得了α-β 及ω-β 相变的相边界, 构建了Zr的参考相图. 同时, 也获得了β-Zr的高温状态方程及热膨胀系数, 能够为构建Zr的全区物态方程提供有益的参考.     

石墨烯低温热膨胀和声子弛豫时间随温度的变化规律

考虑到原子非简谐振动和电子-声子相互作用,用固体物理理论和方法研究了石墨烯格林艾森参量和低温热膨胀系数以及声子弛豫时间随温度的变化规律,探讨了原子非简谐振动项对它们的影响.结果表明:1)在低于室温的温度范围内,石墨烯的热膨胀系数为负值,随着温度的升高,其热膨胀系数的绝对值单调增加,室温热膨胀系数为-3.64×10~(-6)K~(-1);2)简谐近似下的格林艾森参量为零.考虑到非简谐项后,格林艾森参量在1.40-1.42之间并随温度升高而缓慢增大,几乎成线性关系,第二非简谐项对格林艾森参量的影响小于第一非简谐项;3)石墨烯声子弛豫时间随着温度的升高而减小,其中,温度很低(T10 K)时变化很快,此后变化很慢,当温度不太低(T300 K)时,声子弛豫时间与温度几乎成反比关系.     

Lattice dynamics,third order elastic constants and thermal expansion of rhenium

A systematic investigation of the lattice dynamics, second and third order elastic constants and thermal expansion has been carried out in rhenium using Keating's approach. The ten third-order elastic constants of rhenium are calculated using four anharmonic parameters. The present model reproduces the measured pressure derivatives of the second order elastic constants well. The low and high temperature limits γ̄L and γ̄H of the lattice thermal expansion are calculated. The agreement between the calculated γ̄H and that obtained by Gschneidner from the thermal expansion and specific heat data of rhenium is satisfactory.     

Dispersion relations and lattice thermal expansion of dysprosium

The lattice dynamics, second and third order elastic constants and the lattice thermal expansion of dysprosium have been calculated using Keating's approach. The ten third order elastic constants are calculated using four anharmonic parameters. The present model reproduces the measured pressure derivatives of the second order elastic constants of dysprosium well. The low and high temperature limits γ?_I and γ?_H of the lattice thermal expansion are evaluated and the agreement between the calculated γ?_H and that obtained from the thermal expansion and specific heat data is satisfactory.     

Thermal Expansion and Deformation of Graphene

Taking into consideration short-atomic-range interactions and anharmonic effects,we calculate the thermal expansion coefficients,Gruneisen parameters,the elastic modulus of graphene varying with temperature and the phonon frequency.The anharmonic effects associated with the graphene deformation are also discussed.The results show that the value of thermal expansion coefficient is negative in the moderate temperature range,and it becomes positive when the temperature grows to be higher than a certain value.The change rate of elastic modulus with respect to temperature and pressure are calculated,and phonon frequencies are estimated.In the process of graphene thermal expansion,it is accompanied with the change of bond length and the rotation around the axis normal to the plane.Our results indicate that the effects due to the bond change are more significant than that of the rotation.We also show that if anharmonic effects are ignored,the thermal expansion coefficient and the Gruneisen parameters are zero,and the elastic modulus and the phonon frequency are constant.If anharmonic effects are considered up to the second term,these values will vary with temperature,and become closer to the experimental value.The higher the temperature is,the more significant the anharmonic effects become.     

Analytic Debye-Grüneisen equation of state for a generalized Lennard-Jones solids

The approximate method to treat the practical quantum anharmonic solids proposed by Hardy,Lacks and Shukla is reformulated with explicit physical meanings.It is shown that the quantum effect is important at low temperature,it can be treated in the harmonic framework.and the anharmonic effect is important at high temperature and tends to zero at low temperature,it can be treated by using a classical approximation.The alternative formulation is easier for various applications,and is applied to a Debye-Grueneisen solid with the generalized Lennard-Jones intermolecular interaction.The expressions for the Debys temperature and Grueneisen parameter as a function of volume are analytically derived.The analytic equation of state is applied to predict the thermodynamic properties of solid xenon at normal-pressure with the nearest-neighbour Lennard-Jones interaction,and is further applied to research the properties of solid xenon and krypton at high pressure by using an all-neighbour Lennard-Jones interaction.The theoretical results are in agreement with the experiments.     

Effects of temperature on the structure of neutron stars at high temperature

In Newtonian physics, higher temperature leads to higher thermal pressure, which provides stronger support against the gravitational contraction of stars. However, in the temperature range of tens of MeV involved in the evolution of a proto-neutron star or a higher massive neutron star, the effects of temperature are richer. We showed that, for a high temperature neutron star (HTNS) constructed with a realistic equation of state (EOS), the HTNS may expand or contract during cooling, the central density may increase or decrease, the quasi-normal mode oscillation frequencies may increase or decrease, and in particular, (i) independent of the EOS, for a HTNS of a given mass, there exists a maximum temperature \(T_{max}\) that it could ever attend at birth (with the value of \(T_{max}\) different for different EOS), and (ii) for the Hempel EOS and the Shen EOS, there is a range of mass that the HTNS may gravitationally collapse after a period of radiative cooling; however, for the Lattimer–Swesty EOS and Banik EOS, no delayed collapse is possible. Our study, which describes the cooling of HTNSs with simple quasi-stationary TOV sequences, provides an understanding of the effects of the thermal energy/pressure at high temperature, and a demonstration that different EOSs can lead to qualitatively different evolution paths.     

Lattice dynamics,third order elastic constants and thermal expansion of gadolinium

A systematic investigation of the lattice dynamics, second and third order elastic constants and the temperature variation of the effective Grüneisen functions has been carried out in gadolinium using Keating's approach. The ten third order elastic constants are calculated using five anharmonic parameters. The present model reproduces the measured pressure derivatives of the second order elastic constants of gadolinium well. The low and high temperature limits gg_L and gg_H of the lattice thermal expansion are evaluated. The agreement between the calculated gg_H and that obtained from the thermal expansion and specific heat data of gadolinium is good.     

沿等压路径求解疏松材料Hugoniot关系的微分方程组及其求解

 根据固体材料的三项式物态方程和Grüneisen物态方程,导出了沿等压路径求解疏松材料冲击温度和压缩体积随初始密度变化的微分方程组。从体积的微分方程出发,在假定Wu-Jing参量为常数的前提下,导出了冲击压缩体积和体积-焓物态方程的Wu-Jing表达式。采用数值差分方法求解微分方程组,计算了疏松铜的冲击压缩特性,并与文献中部分实验数据进行了比较,特别强调了热电子对冲击压缩体积、冲击温度和Wu-Jing参数的贡献。还讨论了Grüneisen物态方程与Wu-Jing物态方程的内在联系及后者的适用范围。