Molecular dynamics simulation of P-V-T relationship of ZnO with rock-salt structure using pair-wise interactions

Molecular dynamics (MD) method is used to investigate the behavior of the pressure-volume-temperature (P-V-T) relationship, lattice constant and thermal expansivity for ZnO with rock-salt structure at high pressures and temperatures. The interionic potential is taken to be the sum of pair-wise additive Coulomb, van der Waals attraction, and repulsive interactions. The isothermal and isobaric properties are discussed from the corresponding P-V-T relationship, and it is shown that the MD simulation is successful in reproducing the measured volumes of ZnO over a wide range of temperature and pressure. Meanwhile, the equations of state parameters including lattice constant, linear thermal expansion coefficient, and isothermal bulk modulus are calculated and compared with the available experimental data and the latest theoretical results. At an extended pressure and temperature range, P-V-T relationship, lattice constant, and linear thermal expansion coefficient have been predicted. The structural and thermodynamic properties of ZnO with rock-salt structure are summarized in the pressure 0-100 GPa ranges and the temperature up to 3100 K.     

Effects of pressure and temperature on the isothermal bulk modulus of CaO

Molecular dynamics (MD) simulations have been performed to investigate the effects of pressure and temperature on the isothermal bulk modulus of CaO using pair-wise interactions that include polarization effects through the shell model (SM). The dependence of isothermal bulk modulus B_T of CaO on the compression ratio V/V₀ and pressure P have been obtained from MD runs at T=300 K, and compared with the available experimental data and other theoretical results. A good agreement between theory and experiment is obtained. Meanwhile, B_T dependence on temperature T at zero pressure is investigated. At extended pressure and temperature ranges, SM-MD method has also been carried out for predicting the P-V-T equation of state and isothermal bulk modulus at different temperatures along the isotherms 0, 1000, 2000, 3000, and 4000 K, and at different pressures along the isobars 5, 15, 30, 40, and 50 GPa for CaO, respectively.     

Analysis of P-V-T relationships and thermodynamic properties for some alkali halides

In this paper, Pandey approximation for the volume dependence of Anderson-Grüneisen parameter at fixed pressure, Anderson approximation for the temperature dependence of thermal expansivity, the equations of thermal expansivity along isobars derived by Shanker et al., and the presented approximation for the volume dependence of Anderson-Grüneisen parameter along isobars, have been used to study and predict the pressure-volume-temperature (P-V-T) data and the variations of the volume expansion coefficient and isothermal bulk modulus with temperature and pressure for NaCl, CsCl, LiF, NaF crystals, up to 30 GPa and in the temperature range 298-1073 K. The calculated values are compared with each other. It is found that these equations-of-state are valid and present good agreement with the available experimental data.     

Thermal equation-of-state of osmium: a synchrotron X-ray diffraction study

The pressure-volume-temperature behavior of osmium was studied at pressures and temperatures up to 15 GPa and 1273 K. In situ measurements were conducted using energy-dispersive synchrotron X-ray diffraction in a T-cup 6-8 high pressure apparatus. A fit of room-temperature data by the third-order Birch-Murnaghan equation-of-state yielded isothermal bulk modulus K₀=435(19) GPa and its pressure derivative K′₀=3.5(0.8) GPa. High-temperature data were analyzed using Birch-Murnaghan equation of state and thermal pressure approach. The temperature derivative of bulk modulus was found to be −0.061(9) GPa K⁻¹. Significant anisotropy of osmium compressibility was observed.     

Elastic, thermal and structural properties of platinum

The thermal equation of state (EOS) for platinum has been calculated to 300 GPa and 3000 K using ab initio molecular dynamics employing the local density approximation (LDA) and the projector augmented-wave methods (PAW). Direct ab initio molecular dynamics avoids the simplifying assumptions inherent in empirical treatments of thermoelasticity. A third-order Birch-Murnaghan equation EOS fitted to the 300 K data yielded an isothermal bulk modulus of B_T0=290.8 GPa and a pressure derivative of B_T′=5.11, which are in better agreement with the measured values than those obtained by previous calculations. The high-temperature data were fitted to a thermal pressure EOS and a Mie-Grüneisen-Debye EOS. The resulting calculated thermal expansion coefficient, α₀, temperature derivative of the isothermal bulk modulus, (∂B_T/∂T)_V, and second temperature derivative of the pressure, (∂²P/∂T²)_V, were 1.94×10⁻⁵ K⁻¹, −0.0038 GPa K⁻¹, and 1.7×10⁻⁷ GPa² K⁻², respectively. A fit to the Mie-Grüneisen-Debye EOS yielded values for the Grüneisen parameter, γ₀, and its volume dependence parameter, q, of 2.18 and 1.75, respectively. An analysis of our data revealed a strong volume dependence of the thermal pressure of platinum. We also present a qualitative analysis of the effects of intrinsic anharmonicity from the calculated Grüneisen parameter at high temperatures.     

X-ray diffraction study of molybdenum diselenide to 35.9 GPa

In situ high-pressure angle dispersive synchrotron X-ray diffraction studies of molybdenum diselenide (MoSe₂) were carried out in a diamond-anvil cell to 35.9 GPa. No evidence of a phase transformation was observed in the pressure range. By fitting the pressure-volume data to the third-order Birch-Murnaghan equation of state, the bulk modulus, K_0T, was determined to be 45.7±0.3 GPa with its pressure derivative, K′_0T, being 11.6±0.1. It was found that the c-axis decreased linearly with pressure at a slope of −0.1593 when pressures were lower than 10 GPa. It showed different linear decrease with the slope of a −0.0236 at pressures higher than 10 GPa.     

High pressure X-ray diffraction study of ReS2

The high-pressure behavior of rhenium disulfide (ReS₂) has been investigated to 51.0 GPa by in situ synchrotron X-ray diffraction in a diamond anvil cell at room temperature. The results demonstrate that the ReS₂ triclinic phase is stable up to 11.3 GPa, at which pressure the ReS₂ transforms to a new high-pressure phase, which is tentatively identified with a hexagonal lattice in space group P6?m2. The high-pressure phase is stable up to the highest pressure in this study (51.0 GPa) and not quenchable upon decompression to ambient pressure. The compressibility of the triclinic phase exhibits anisotropy, meaning that it is more compressive along interlayer directions than intralayer directions, which demonstrates the properties of the weak interlayer van der Waals interactions and the strong intralayer covalent bonds. The largest change in the unit cell angles with increasing pressures is the increase of β, which indicates a rotation of the sulfur atoms around the rhenium atoms during the compression. Fitting the experimental data of the triclinic phase to the third-order Birch-Murnaghan EOS yields a bulk modulus of K_OT=23±4 GPa with its pressure derivative K_OT′= 29±8, and the second-order yields K_OT=49±3 GPa.     

Volume and structural study of Fe64Mn36 anti-ferromagnetic Invar alloy under high pressure

We have investigated the pressure variation of the volume and structure of an FCC Fe₆₄Mn₃₆ anti-ferromagnetic Invar alloy. The inclination of the pressure-volume (P-V) curve of the FCC structure becomes discontinuous at a pressure of 4 GPa. According to the bulk modulus at zero pressure estimated by the Birch-Murnaghan equation of state, the pressure between 4 and 10 GPa is 33 GPa larger than that at a pressure below 4 GPa. Considering previous experiments on magnetism at high pressure the Neel temperature at 4 GPa almost decreases to room temperature. These results suggest that the increase in the bulk modulus by 33 GPa can be attributed to the pressure-induced magnetic phase transition from anti-ferromagnetism to paramagnetism. Volume at zero pressure was estimated using the Birch-Murnaghan equation of state. The volume of FCC structure in the anti-ferromagnetic state was 1.17% larger than the volume in the paramagnetic state, namely, the spontaneous magnetostriction was 1.17%. Pressure-induced structural transition from FCC to HCP occurs with an increase in the pressure, especially at up to 5 GPa. The value of c/a is 1.62; this value almost corresponds to that of an ideal HCP structure. The bulk modulus of the HCP structure estimated by the Birch-Murnaghan equation of state is larger than that of the FCC structure, and the volume/atom ratio is smaller than that of the FCC structure.     

X-ray diffraction and infrared spectroscopy of monazite-structured CaSO4 at high pressures: Implications for shocked anhydrite

X-ray diffraction and infrared spectroscopy of CaSO₄ are conducted to pressures of 28 and 25 GPa, respectively. A reversible phase transition to the monoclinic monazite-structure occurs gradually between 2 and ∼5 GPa with a highly pressure-dependent volume change of ∼6-8%. A second-order fit of the X-ray data to the Birch-Murnaghan equation of state yields a bulk modulus (K) of 151.2 (±21.4) GPa for the high-pressure monoclinic phase. In the high-pressure infrared spectrum, the infrared-active asymmetric stretching and bending vibrations of the sulfate tetrahedra split at the phase transition, in accord with the results of factor group analysis. Additionally, the tetrahedral symmetric stretching vibration, which is weak in the anhydrite phase, becomes strongly resolved at the transition to the monazite structure. The infrared results indicate that the sulfate tetrahedra are more distorted in the monazite-structured phase than in anhydrite. Kinetic calculations indicate that the anhydrite to monazite transformation may generate the phase transition observed near 30 GPa under shock loading in CaSO₄. Our results indicate that the anhydrite- and monazite-structured phases may be the only phases that occur under shock loading of CaSO₄ to pressures in excess of 100 GPa.     

Equation of state of aluminum carbide Al4C3

Quasi-hydrostatic compression of aluminum carbide, Al₄C₃ has been studied to 6 GPa at room temperature using energy-dispersive X-ray powder diffraction with synchrotron radiation. A fit of the experimental p-V data to the Birch equation of state yields the values of the bulk modulus, B₀, of 130(5) GPa and the first pressure derivative of the bulk modulus, B′₀, of 4.6(9). The compression is found to be anisotropic, with the a-axis being more compressible than the c-axis.     

Compression behavior of nanocrystalline anatase TiO2

We present a synchrotron X-ray diffraction study of pressure-induced changes in nanocrystalline anatase (with a crystallite size of 30-40 nm) to 35 GPa. The nanoanatase was observed to a pressure above 20 GPa. Direct transformation to the baddeleyite-TiO₂ polymorph was seen at 18 GPa. A fit of the pressure versus volume data to a Birch-Murnaghan equation yielded the following parameters: zero-pressure volume, V₀=136.15 Å³, bulk modulus, K_T=243(3) GPa, and the pressure derivative of bulk modulus, K′=4 (fixed). The bulk modulus value obtained for the nanocrystalline anatase is about 35% larger than that of the macrocrystalline counterpart.     

Phase diagram of kaolinite from Molecular Dynamics calculations

Structural and thermal behaviors of uncharged 1:1 phyllosilicates kaolinite were investigated from molecular dynamics simulations based on the CLAYFF force field. The focus is on the variation of structural properties including bulk modulus with pressure from 0 to 20 GPa under various range of temperature. The largest bulk modulus between the pressures of 200 and 800 MPa varies from 80 GPa at 298 K to 50 GPa at 1473 K. The obtained value of C_p varies between 7.8 and 13.6 Kcal mol⁻¹ K⁻¹ in the pressure range of 0.1 MPa–20 GPa. Besides, a huge difference was noticed regarding the computed properties at the superheating point. Finally, we show the relationship between superheating point temperature and pressure leading to a phase diagram of kaolinite.     

Equation of state and thermal stability of Al3BC

The lattice parameters of Al₃BC have been measured up to 5 GPa at ambient temperature using energy-dispersive X-ray powder diffraction with synchrotron radiation. A fit to the experimental p-V data using Birch-Murnaghan equation of state gives values of the Al₃BC bulk modulus 116(4) GPa and its first pressure derivative 9(2). In the 1.6-4.8 GPa range at temperatures above 1700 K Al₃BC undergoes incongruent melting that results in the formation of Al₃BC₃, AlB₂ and liquid aluminum.     

Thermal expansivity and bulk modulus of ZnO with NaCl-type cubic structure at high pressures and temperatures

The thermal expansivity and bulk modulus of ZnO with NaCl-type cubic structure were estimated by using the constant temperature and pressure molecular dynamics technique with effective pair potentials which consist of the Coulomb, dispersion, and repulsion interaction at high pressures and temperatures. It is shown that the calculated thermodynamic parameters including linear thermal expansion coefficient, isothermal bulk modulus and its pressure derivative are in good agreement with the available experimental data and the latest theoretical results. At an extended pressure and temperature ranges, linear thermal expansion coefficient and isothermal bulk modus have also been predicted. The thermodynamic properties of ZnO with NaCl-type cubic structure are summarized in the pressure 0–150 GPa ranges and the temperature up to 3000 K.     

Structural phase transitions in brookite-type TiO2 under high pressure

We have studied polycrystalline brookite TiO₂ using energy-dispersive X-ray diffraction at pressures up to 27.8 GPa and derived an ambient-pressure bulk modulus of 255 GPa using Birch-Murnaghan's equations of state with a fixed value of 4 as its first derivative. The transition from brookite-type to baddeleyite-type was observed to start at 15.8 GPa and finished at 22.8 GPa. Upon decompression, the α-PbO₂ structure appeared at 3.5 GPa and the baddeleyite-type structure remained down to 1.6 GPa, the lowest pressure in the present work.     

High-pressure powder diffraction study of TaO2F

TaO₂F, with a ReO₃-type structure, has been studied at up to 12.8 GPa using monochromatic synchrotron powder diffraction and diamond anvil cells. Two-phase transitions at ∼0.7 and 4 GPa were observed on compression. Below ∼0.7 GPa the cubic material was found to have a bulk modulus (K₀) of 36(3) GPa (K_p fixed at 4.0), similar to that reported for NbO₂F but much smaller than that of ReO₃. Immediately above 0.7 GPa on compression, the diffraction data were not fully consistent with a VF₃-type structure as previously proposed for NbO₂F. On decompression, the data between 8 and 4 GPa could be satisfactorily attributed to a single R-3c phase with a VF₃-type structure and an average bulk modulus of 60(2) GPa.     

Analysis of thermoelastic behaviour of MgO at high temperatures and high pressures

The thermoelastic behaviour of MgO has been studied for the temperature range (300-3000 K) under different compressions down to V/V₀=0.3. It has been shown that a comprehensive study of the thermoelastic properties of MgO can be made with the help of the Anderson-Isaak equation for thermal expansivity and the Vinet equation of state taken together. We have estimated the values of thermal expansivity α, isothermal bulk modulus K_T, their variations with pressure and temperature, the Anderson-Gruneisen parameter and the change in entropy with compression for MgO along isotherms at different temperatures. The results have been discussed and compared with the corresponding values reported in the recent literature.     

Structure and high pressure behaviour of organic crystals based on aromatically substituted 1,3,4-oxadiazoles

The crystalline structure of a new compound containing the 1,3,4-oxadiazole moiety, 4-(5-methyl-1,3,4-oxadiazole-2yl-)-N,N′-dimethyl-phenylamine (MODPA) was determined. It shows a monoclinic structure with space group P2₁/c and lattice parameters: a=1.02997(6), b=0.64840(4), c=1.58117(10) nm and β=99.4820(10)°. To study the intermolecular interactions in oxadiazole containing organic crystals, X-ray studies on MODPA and 2,5-diphenyl-1,3,4-oxadiazole (DPO) were performed up to 5 GPa at room temperature. The Murnaghan equation of state is used to describe the compression behaviour of both substances. From these results, the bulk modulus and its pressure derivative were determined. The values obtained are: K₀=6.3 GPa and K′₀=6.8 for MODPA and K₀=7.3 GPa and K′₀=6.7 for DPO. Additionally, measurements under increasing temperature at ambient pressure were carried out to evaluate the thermal expansion coefficient: α=1.8×10⁻⁴ K⁻¹ for MODPA and α=1.9×10⁻⁴ K⁻¹ for DPO.     

Electronic and magnetic properties and their response to pressure in Ni2MnB

The electronic structure and magnetic properties of Ni₂MnB upon pressure up to 20 GPa have been studied by using the density functional theory (DFT) method. The results indicate that ferromagnetic ordered Ni₂MnB in L2₁ structure is more stable than the nonmagnetic one. The magnetic moments of Ni and Mn atoms as well as the total magnetic moment of Ni₂MnB are found to decrease weakly with increasing pressure. The pressure derivative of the total magnetic moment is −3.07×10⁻³ GPa⁻¹. The equilibrium bulk modulus and its derivative from the Murnaghan equation of state (EOS) are B₀=247.7 GPa, B′=4.98.     

In situ observations of temperature- and pressure-induced phase transitions in TiH2: Angle-dispersive and synchrotron energy-dispersive X-ray diffraction studies

We investigated the behavior of the structure of titanium hydride (TiH₂), an important compound in hydrogen storage research, at elevated temperatures (0-120 °C) and high pressures (1 bar-34 GPa). Temperature-induced changes of TiH₂ as indicated in the alteration of the ambient X-ray demonstrated a cubic to tetragonal phase transition occurring at about 17 °C. The main focus of this study was to identify any pressure-induced structural transformations, including possible phase transitions, in TiH₂. Synchrotron X-ray diffraction studies were carried out in situ (diamond anvil cell) in a compression sequence up to 34 GPa and in subsequent decompression to ambient pressure. The pressure evolution of the diffraction patterns revealed a cubic (Fm-3m) to tetragonal (I4/mmm) phase transition at 2.2 GPa. The high-pressure phase persisted up to 34 GPa. After decompression to ambient conditions the observed phase transition was completely reversible. A Birch-Murnaghan fit of the unit cell volume as a function of pressure yielded a zero-pressure bulk modulus K₀=146(14) GPa, and its pressure derivative K′₀=6(1) for the high-pressure tetragonal phase of TiH₂.