Thermodynamic properties of the alkali metals at high temperatures and high pressures using mean-field potential model

We have carried out theoretical calculations for the thermodynamic properties of alkali metals at high temperatures and pressures. In the present study, we have used recently proposed mean-field potential (MFP) model of Wang and Li for evaluating the vibrational contribution to the total free energy. The recently proposed local pseudopotential of Filohais et al. is used to construct mean-field potential in terms of the total energy-volume relation. We have calculated static and shock compressions, bulk moduli (B_T and B_S), thermal expansion, specific heats (C_P and C_V), anharmonic contribution to the specific heat , enthalpy (E_H) and temperature along shock Hugoniot. We find that our computed results (except for C_P and C_V) are in good agreement with experiments as well as comparable with those generated by using first principles methods and other theoretical models. The present study indicates that in comparison with other theoretical models the present model is computationally simple, physically transparent and reliable to study the thermodynamic properties in the high pressures and high temperatures environment.     

High temperature and pressure thermodynamics of strontium: A macroscopic approach

Close proximity of d-bands (above) to the Fermi level (E _F) makes the heavy alkaline earth metals (Ca, Sr and Ba) fairly sensitive to external influences like temperature and pressure. Softening of some of the phonon modes at high temperatures and/or pressures implies that anharmonic effects can play an important role in determining lattice dynamics and related properties. In the conventional approach, phonon density of states (p-dos) have to be calculated at each volume to compute free energy and thereby the other thermodynamic properties, which is computationally quite demanding. Using an alternative technique, the mean-field potential (MFP) approach was combined with the relatively soft local pseudopotential to obtain the free energy at different temperatures and pressures. The results for phonon frequency shifts at finite temperatures using the MFP approach and those calculated from p-dos within the quasiharmonic approximation are very similar. This validates the use of the MFP approach coupled with the local pseudopotential to estimate vibrational response of the system at high-temperature and high-pressure environments. The present scheme was used to study various thermophysical properties for elemental strontium at elevated temperatures and pressures, including the high-pressure melting curve and temperature along the shock Hugoniot. Computed results are affirmatively compared and analyzed with other reported data. The present scheme completely bypasses traditional cumbersome calculations, and it is computationally convenient yet accurate.     

Mean-field potential calculations of high-pressure equation of state for BeO

A systematic study of the Hugoniot equation of state, phase transition, and the other thermodynamic properties including the Hugoniot temperature, the electronic and ionic heat capacities, and the Gruneisen parameter for shockcompressed BeO, has been carried out by calculating the total free energy. The method of calculations combines first-principles treatment for 0 K and finite-T electronic contribution and the mean-field-potential approach for the vibrational contribution of the lattice ion to the total energy. Our calculated Hugoniot is in good agreement with the experimental data.     

First-Principles Calculations of Elastic and Thermal Properties of Lanthanum Hexaboride

The plane-wave pseudopotentiai method using the generaiized gradient approximation within the framework of density functional theory is applied to anaylse the bulk modulus, thermal expansion coefficient and heat capacity of LaB6. The quasi-harmonic Debye model, using a set of total energy versus volume obtained with the plane-wave pseudopotential method, is applied to the study of the thermal properties and vibrationai effects. We analyse the bulk modulus of LaB6 up to 1500 K. The elastic properties calculations show that our system is mechanically stable. For the heat capacity and the thermal expansion, significant differences in properties are observed above 300K. The calculated zero pressure bulk modulus is in good agreement with the experimental data. Moreover, the Debye temperatures are determined from the non-equilibrium Gibbs functions and compared to available data.     

General construction of mean-field potential and its application to the multiphase equations of state of tin

The mean-field potential (MFP) approach is an efficient way to evaluate the free energy contribution of ion motions for both solid and liquid states. In this paper the MFP is generally constructed with a volume-dependent term and a shape function. The former is derived in accordance with quasi-harmonic approximation. The latter is given semi-empirically. Application to multiphase equations of state for β-, γ- and liquid-tin has been examined. The theoretical phase diagram and thermodynamic properties of isotherm, thermal expansion, heat capacity, Hugoniot states as well as phase transitions are all in excellent agreement with experiments.     

Thermophysical properties of iridium at finite temperature

The bulk properties of materials in an extreme environment such as high temperature and high pressure can be understood by studying anharmonic effects due to the vibration of lattice ions and thermally excited electrons.In this spirit,in the present paper,anharmonic effects are studied by using the recently proposed mean-field potential(MFP) approach and Mermin functional which arise due to the vibration of lattice ions and thermally excited electrons,respectively.The MFP experienced by a wanderer atom in the presence of surrounding atoms is constructed in terms of cold energy using the local form of the pseudopotential.We have calculated the temperature variation of several thermophysical properties in an extreme environment up to melting temperature.The results of our calculations are in excellent agreement with the experimental findings as well as the theoretical results obtained by using first principle methods.We conclude that presently used conjunction scheme(MFP+pseudo potential) is simple computationally,transparent physically,and accurate in the sense that the results generated are comparable and sometimes better than the results obtained by first principle methods.Local pseudopotential used is transferable to extreme environment without adjusting its parameters.     

Thermodynamic properties of rhodium at high temperature and pressure by using mean field potential approach

The thermophysical properties of rhodium are studied up to melting temperature by incorporating anharmonic effects due to lattice ions and thermally excited electrons. In order to account anharmonic effects due to lattice vibrations, we have employed mean field potential (MFP) approach and for thermally excited electrons Mermin functional. The local form of the pseudopotential with only one effective adjustable parameter r_c is used to construct MFP and hence vibrational free energy due to ions – F_ion. We have studied equation of state at 300 K and further, to access the applicability of present conjunction scheme, we have also estimated shock-Hugoniot and temperature along principle Hugoniot. We have carried out the study of temperature variation of several thermophysical properties like thermal expansion (β_P), enthalpy (E_H), specific heats at constant pressure and volume (C_P and C_V), specific heats due to lattice ions and thermally excited electrons (C_V^ion and C_V^el, isothermal and adiabatic bulk moduli (B_T and B_s) and thermodynamic Gruneisen parameter (γ_th) in order to examine the inclusion of anharmonic effects in the present study. The computed results are compared with available experimental results measured by using different methods and previously obtained theoretical results using different theoretical philosophy. Our computed results are in good agreement with experimental findings and for some physical quantities better or comparable with other theoretical results. We conclude that local form of the pseudopotential used accounts s-p-d hybridization properly and found to be transferable at extreme environment without changing the values of the parameter. Thus, even the behavior of transition metals having complexity in electronic structure can be well understood with local pseudopotential without any modification in the potential at extreme environment. Looking to the success of present scheme (MFP + pseudopotential) we would like to extend it further for the study of liquid state properties as well as thermophysical properties of d and f block metals.     

First-principles study of lattice dynamics and thermodynamics ofosmium under pressure

We have performed the first-principles linear response calculations of the lattice dynamics, thermal equation of state and thermodynamical properties of hcp Os metal by using the plane-wave pseudopotential method. The thermodynamical properties are deduced from the calculated Helmholtz free energy by taking into account the electronic contribution and lattice vibrational contribution. The phonon frequencies at Gamma point are consistent with experimental values and the dispersion curves at various pressures have been determined. The calculated volume, bulk modulus and their pressure derivatives as a function of temperature are in excellent agreement with the experimental results. The calculated specific heat indicates that the electronic contribution is important not only at very low temperatures but also at high temperatures due to the electronic thermal excitation. The calculated Debye temperature at a very low temperature is in good agreement with experimental values and drops to a constant until 100~K.     

金属铍热力学性质的理论研究

使用第一性原理方法结合平均场模型研究了压力从0到150GPa、温度从0到1500K，金属铍六角密排结构(hcp)的热力学性质，包括铍的常态性质，等温高压物态方程，以及常压下平衡体积、体弹模量随温度的变化，Hugoniot曲线等.0K物态方程由广义梯度近似下的密度泛函理论计算，粒子热运动的贡献由平均场模型计算.由于铍的Debye温度比较高，计算自由能时考虑了零点振动能修正.计算结果与已有的静力学和冲击波实验数据符合得非常好. 关键词： 热力学性质 物态方程 第一原理计算     

Structural and thermodynamic properties of the cubic perovskite BiAlO3

The structural and elastic properties of the cubic perovskite-type BiAlO₃ are studied using the pseudopotential plane wave method within the local density approximation. The calculated structural parameters are in good agreement with previous calculations. The elastic constants are calculated using the static finite strain technique. Thermal effects on some macroscopic properties of BiAlO₃ are predicted using the quasi-harmonic Debye model, in which the lattice vibrations are taken in account. We have obtained successfully the variations of the lattice constant, volume expansion coefficient, heat capacities and Debye temperature with pressure and temperature in the ranges of 0-30 GPa and 0-1000 K.     

First-Principles Calculations of Elastic and Thermal Properties of Molybdenum Disilicide

The first-principles plane-wave pseudopotential method using the generalized gradient approximation within the framework of density functional theory is applied to anaylse the equilibrium lattice parameters, six independent elastic constants, bulk moduli, thermal expansions and heat capacities of MoSi2. The quasi-harmonic Debye model, using a set of total energy versus cell volume obtained with the plane-wave pseudopotential method, is applied to the study of the elastic properties, thermodynamic properties and vibrational effects. The calculated zero pressure elastic constants are in overall good agreement with the experimental data. The calculated heat capacities and the thermal expansions agree well with the observed values under ambient conditions and those calculated by others. The results show that the temperature has hardly any effect under high pressure.     

Theoretical Study of Elastic Properties of Tungsten Disilicide

The plane-wave pseudopotential method using the generalized gradient approximation within the framework of density functional theory is applied to analyse the lattice parameters, elastic constants, bulk moduli, shear moduli and Young's moduli of WSi₂. The quasi-harmonic Debye model, using a set of total energy versus cell volume obtained with the plane-wave pseudopotential method, is applied to the study of the elastic properties and vibrational effects. The athermal elastic constants of WSi₂ are calculated as a function of pressure up to 35GPa. The relationship between bulk modulus and temperature up to 1200K is also obtained. Moreover, the Debye temperature is determined from the non-equilibrium Gibbs function. The calculated results are in good agreement with the experimental data.     

Phase Transition and Phonon Spectrum of Zinc-Blende Structure ZnX （X = S,Se, Te）

Calculations have been performed to investigate the pressure-induced solid-solid phase transitions and the mechanical stability for three zinc-blende II-VI semiconductor compounds： ZnS, ZnSe, ZnTe by ab initio plane-wave pseudopotential density functional theory （DFT）. Using the generalized gradient approximation （GGA） for exchange and correlation in the scheme of Perdew-Wang 1991 （P Wgl ）, the ground state properties and equation of state are obtained, which are well consistent with the experimental data available and other calculations. On the basis of the forth-order Birch-Murnaghan equation of states, the transition pressures Pt are determined through the analysis of enthalpy variation with pressure. A linear-response approach is used to calculate the frequencies of the phonon dispersion. Finally, by the calculations of phonon frequencies, some thermodynamic properties such as the vibrational contribution to the Helmholtz free energy （F）, enthedpy （H）, entropy （S）, and the heat capacity （Cv ） are also successfully obtained.     

First-principles study of the anisotropic thermal expansion of hcp metals Be and Y

A first-principles study of the anisotropic thermal expansion of hcp metals Be and Y is reported. According to quasiharmonic approximation, the phonon spectra were computed at a set of lattice parameters using the pseudopotential plane wave method with the local density approximation in the framework of the density functional perturbation theory. The free energies were obtained according to the calculated phonon spectra and thermal properties such as specific heat at constant volume (pressure) were calculated. The electronic contribution to specific heat was found important to metal Y not only at very low temperature but also over room temperature. The calculated results are in good agreement with available experimental data in a wide range of temperature.     

On the stability of the polaron in one dimensional disordered systems

A one-dimensional disordered system of electrons described by a tight binding model interacting with vibrational degrees of freedom (in harmonic approximation) is considered. A stable configuration is determined by a numerical minimization of the total energy which is based on the adiabatic approximation. The behaviour of the electron density (charge density wave) and the density of states is analysed. The localization properties are investigated as well. In contrast to the corresponding disordered system with vanishing electron-phonon coupling the present model has an energy gap. The formation of the gap and the polaron band is shown to be quite different for both onsite and intersite types of coupling terms. For large disorder, the lattice distortion and the gap disappear if only the vibrational contribution to the intersite coupling is important. They increase, however, if only the vibrational contribution to the site energies is taken into account. In both cases the localization length decreases upon increasing the electron-phonon coupling energy. The results are discussed with respect to low dimensional organic materials and amorphous semiconductors.     

Discontinuities in the specific heat of magnesium and associated latent heat at pressure-induced structural phase transitions using a local first principles pseudopotential

After a successful application of a local first principles pseudopotential to the study of the structural pressure-induced phase transition for magnesium, I now use the same local first principles type of pseudopotential, to study the specific heat as function of the pressure, at 300 K, in the region around the predicted phase transitions. I found that the specific presents two discontinuities, one for each transition of phase. These discontinuities are associated to the existence of latent heat at each transition, which has not been yet reported experimentally.     

Pseudo-potential investigations of structural, elastic and thermal properties of tungsten disilicide

The plane-wave pseudopotential method using the generalized gradient approximation within the density functional theory is used to investigate the structure and bulk modulus of WSi₂. The quasi-harmonic Debye model, using a set of total energy versus cell volume obtained with the plane-wave pseudopotential method, is applied to the study of the elastic properties and vibrational effects. We have analysed the bulk modulus of WSi₂ up to 1600~K. The major trend shows that the WSi₂ crystal becomes more compressible when the temperature rises and the increase of compressibility leads to the decrease of Debye temperature. The predicted temperature and pressure effects on the thermal expansion, heat capacity and Debye temperatures are determined from the non-equilibrium Gibbs functions and compared with the data available.     

Theoretical investigations on the structural, lattice dynamical and thermodynamic properties of XC (X=Si, Ge, Sn)

First-principles calculations, which is based on the plane-wave pseudopotential approach to the density functional perturbation theory within the local density approximation, have been performed to investigate the structural, lattice dynamical, and thermodynamic properties of SiC, GeC, and SnC. The results of ground state parameters, phase transition pressure and phonon dispersion are compared and agree well with the experimental and theoretical data in the previous literature. The obtained phonon frequencies at the zone-center are analyzed. We also used the phonon density of states and quasiharmonic approximation to calculate and predict some thermodynamic properties such as entropy, heat capacity, internal energy, and phonon free energy of SiC, GeC, and SnC in B3 phase.     

Pressure induced transitions in calcium: a tight-binding approach

A tight-binding (TB) hamiltonian for calcium is built with a high precision parametrization technique based on density functional calculations of the energy bands and the total energy at various lattice volumes. The new set of TB parameters is appropriate to study phenomena under pressures as high as 20 GPa. Specifically, both the metal to nonmetal transition at 4 GPa and the structural transition fcc to bcc at 19 GPa are well reproduced. These transitions and several static properties are in excellent agreement with experiments. Phonon frequencies, plasmon energy, melting temperature and the coefficient of thermal expansion were calculated with a molecular dynamics scheme of this TB hamiltonian.     

Phase transition and thermodynamic properties of ThO2：Quasi-harmonic approximation calculations and anharmonic effects

The thermodynamic properties and the phase transition of ThO2 from the cubic structure to the orthorhombic structure are investigated using the first-principles projector-augmented wave method. The vibrational contribution to Helmholtz free energy is evaluated from the first-principles phonon calculations. The anharmonic contribution to quasi-harmonic free energy is accounted for by using an effective method(2010 Phys. Rev. B 81 172301). The results reveal that at ambient temperature, the phase transition from the cubic phase to the orthorhombic phase occurs at 26.45 GPa, which is consistent with the experimental and theoretical data. With increasing temperature, the transition pressure decreases almost linearly. By comparing the experimental results with the calculation results, it is shown that the thermodynamic properties of ThO2 at high temperature improve substantially after including the anharmonic correction to quasi-harmonic free energy.