scholarly journals Understanding high pressure molecular hydrogen with a hierarchical machine-learned potential

2020 ◽  
Vol 11 (1) ◽  
Author(s):  
Hongxiang Zong ◽  
Heather Wiebe ◽  
Graeme J. Ackland

Abstract The hydrogen phase diagram has several unusual features which are well reproduced by density functional calculations. Unfortunately, these calculations do not provide good physical insights into why those features occur. Here, we present a fast interatomic potential, which reproduces the molecular hydrogen phases: orientationally disordered Phase I; broken-symmetry Phase II and reentrant melt curve. The H2 vibrational frequency drops at high pressure because of increased coupling between neighbouring molecules, not bond weakening. Liquid H2 is denser than coexisting close-packed solid at high pressure because the favored molecular orientation switches from quadrupole-energy-minimizing to steric-repulsion-minimizing. The latter allows molecules to get closer together, without the atoms getting closer, but cannot be achieved within in a close-packed layer due to frustration. A similar effect causes negative thermal expansion. At high pressure, rotation is hindered in Phase I, such that it cannot be regarded as a molecular rotor phase.

Author(s):  
Piotr Rejnhardt ◽  
Marek Drozd ◽  
Marek Daszkiewicz

The phase transition observed in a temperature-dependent experiment at 174 K is unachievable under high-pressure conditions. Negative thermal expansion for phase (II) and negative compressibility for phase (I) were observed. A new salt of 1H-pyrazole-1-carboxamidine, (HPyCA)NO3, for guanylation reaction was obtained in a crystalline form. The compound crystallizes in monoclinic space group P21/c and a phase transition at 174 K to triclinic modification P 1 was found. An unusual increase of the unit-cell volume was observed just after transition. Although the volume decreases upon cooling, it remains higher down to 160 K in comparison to the unit-cell volume of phase (I). The mechanism of the phase transition is connected with a minor movement of the nitrate anions. The triclinic phase was unreachable at room-temperature high-pressure conditions up to 1.27 GPa. On further compression, delamination of the crystal was observed. Phase (I) exhibits negative linear compressibility, whereas abnormal behaviour of the b unit-cell parameter upon cooling was observed, indicating negative thermal linear expansion. The unusual nature of the compound is associated with the two-dimensional hydrogen-bonding network, which is less susceptible to deformation than stacking interactions connecting the layers of hydrogen bonds. Infrared spectroscopy and differential scanning calorimetry measurements were used to investigate the changes of intermolecular interactions during the phase transition.


2020 ◽  
Vol 105 (12) ◽  
pp. 1761-1768
Author(s):  
Gerald V. Gibbs ◽  
David F. Cox ◽  
Nancy L. Ross

Abstract The structures of the silica polymorphs α-quartz and stishovite have been geometry optimized at highly simulated isotropic pressure within the framework of Density Functional Theory. The atoms of the high-pressure polymorph stishovite are virtually incompressible with the bonded radii for Si and O atoms decreasing by only 0.04 and 0.08 Å, respectively, at 100 GPa. In compensating for the increase in the effective interatomic potential associated with the compression of the Si-O bonded interactions, the electron density at the bond critical point between the bonded pair increases from 0.69 to 0.89 e/Å 3. The bonded radii of the Si and O atoms for α-quartz decrease by 0.006 and 0.008 Å, respectively, between 1 bar and 26.4 GPa. The impact of simulated, isotropic pressure on the bonded radii of the atoms for three perovskites YAlO3, LaAlO3, and CaSnO3 was also examined at high pressure. For the YAlO3 perovskite, the bonded radii for Y and Al decrease by 0.06 and 0.05 Å, respectively, at 80 GPa, while the electron density between the bonded atoms increases by 0.12 and 0.15 e/Å3, on average. The calculations also show that the coordination number of the Y atom increases from 9 to 10 while the coordination number of the O1 atom increases concomitantly in the structure from 5 to 6 at 20 GPa. Hence pressure not only promotes an increase in the coordination number of the metal atoms but also a necessary concomitant increase in the coordination number of the O atoms. The bonded radii, determined at a lower pressure between 0.0 and 15 GPa for LaAlO3 and CaSnO3, decrease a smaller amount with the radii for the La and Ca atoms decreasing by 0.03 and 0.04 Å, respectively, while the radii for the smaller Al and Sn atoms decrease by 0.01 and 0.02 Å, respectively. In general, O atoms are more compressible than the metal atoms, but overall the calculations demonstrate that the bonded radii for the atoms in crystals are virtually incompressible when subjected to high pressure. The reason that the bonded radii change little when subjected to high pressure is ascribed to the changes in the effective interatomic potentials that result in increased repulsion when the atoms are squeezed together.


Minerals ◽  
2021 ◽  
Vol 11 (6) ◽  
pp. 607
Author(s):  
Raquel Chuliá-Jordán ◽  
David Santamaria-Perez ◽  
Javier Ruiz-Fuertes ◽  
Alberto Otero-de-la-Roza ◽  
Catalin Popescu

The structure of the naturally occurring, iron-rich mineral Ca1.08(6)Mg0.24(2)Fe0.64(4)Mn0.04(1)(CO3)2 ankerite was studied in a joint experimental and computational study. Synchrotron X-ray powder diffraction measurements up to 20 GPa were complemented by density functional theory calculations. The rhombohedral ankerite structure is stable under compression up to 12 GPa. A third-order Birch–Murnaghan equation of state yields V0 = 328.2(3) Å3, bulk modulus B0 = 89(4) GPa, and its first-pressure derivative B’0 = 5.3(8)—values which are in good agreement with those obtained in our calculations for an ideal CaFe(CO3)2 ankerite composition. At 12 GPa, the iron-rich ankerite structure undergoes a reversible phase transition that could be a consequence of increasingly non-hydrostatic conditions above 10 GPa. The high-pressure phase could not be characterized. DFT calculations were used to explore the relative stability of several potential high-pressure phases (dolomite-II-, dolomite-III- and dolomite-V-type structures), and suggest that the dolomite-V phase is the thermodynamically stable phase above 5 GPa. A novel high-pressure polymorph more stable than the dolomite-III-type phase for ideal CaFe(CO3)2 ankerite was also proposed. This high-pressure phase consists of Fe and Ca atoms in sevenfold and ninefold coordination, respectively, while carbonate groups remain in a trigonal planar configuration. This phase could be a candidate structure for dense carbonates in other compositional systems.


1985 ◽  
Vol 77-78 ◽  
pp. 201-212 ◽  
Author(s):  
Y.J Chabal ◽  
C.K.N Patel

Author(s):  
David R. Allan ◽  
Daniel Bailey ◽  
Nigel Bird ◽  
Alexander J. Blake ◽  
Neil R. Champness ◽  
...  

The mononuclear macrocyclic PdIIcomplexcis-[PdCl2([9]aneS3)] ([9]aneS3= 1,4,7-trithiacyclo-nonane) converts at 44 kbar into an intensely coloured chain polymer exhibiting distorted octahedral coordination at the metal centre and an unprecedented [1233] conformation for the thioether ligand. The evolution of an intramolecular axial sulfur–metal interaction and an intermolecular equatorial sulfur–metal interaction is central to these changes. High-pressure crystallographic experiments have also been undertaken on the related complexes [PtCl2([9]aneS3)], [PdBr2([9]aneS3)], [PtBr2([9]aneS3)], [PdI2([9]aneS3)] and [PtI2([9]aneS3)] in order to establish the effects of changing the halide ligands and the metal centre on the behaviour of these complexes under pressure. While all complexes undergo contraction of the various interaction distances with increasing pressure, only [PdCl2([9]aneS3)] undergoes a phase change. Pressure-induced I...I interactions were observed for [PdI2([9]aneS3)] and [PtI2([9]aneS3)] at 19 kbar, but the corresponding Br...Br interactions in [PdBr2([9]aneS3)] and [PtBr2([9]aneS3)] only become significant at much higher pressure (58 kbar). Accompanying density functional theory (DFT) calculations have yielded interaction energies and bond orders for the sulfur–metal interactions.


Molecules ◽  
2021 ◽  
Vol 26 (22) ◽  
pp. 6831
Author(s):  
Tianming Li ◽  
Junyu Fan ◽  
Zhuoran Wang ◽  
Hanhan Qi ◽  
Yan Su ◽  
...  

The 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105) is a newly energetic material with an excellent performance and low sensitivity and has attracted considerable attention. On the basis of the dispersion-corrected density functional theory (DFT-D), the high-pressure responses of vibrational properties, in conjunction with structural properties, are used to understand its intermolecular interactions and anisotropic properties under hydrostatic and uniaxial compressions. At ambient and pressure conditions, the DFT-D scheme could reasonably describe the structural parameters of LLM-105. The hydrogen bond network, resembling a parallelogram shape, links two adjacent molecules and contributes to the structure stability under hydrostatic compression. The anisotropy of LLM-105 is pronounced, especially for Raman spectra under uniaxial compression. Specifically, the red-shifts of modes are obtained for [100] and [010] compressions, which are caused by the pressure-induced enhance of the strength of the hydrogen bonds. Importantly, coupling modes and discontinuous Raman shifts are observed along [010] and [001] compressions, which are related to the intramolecular vibrational redistribution and possible structural transformations under uniaxial compressions. Overall, the detailed knowledge of the high-pressure responses of LLM-105 is established from the atomistic level. Uniaxial compression responses provide useful insights for realistic shock conditions.


Author(s):  
Keivan Esfarjani ◽  
Gang Chen ◽  
Asegun Henry

Based on first-principles density-functional calculations, we have developed and tested a force-field for silicon, which can be used for molecular dynamics simulations and the calculation of its thermal properties. This force field uses the exact Taylor expansion of the total energy about the equilibrium positions up to 4th order. In this sense, it becomes systematically exact for small enough displacements, and can reproduce the thermodynamic properties of Si with high fidelity. Having the harmonic force constants, one can easily calculate the phonon spectrum of this system. The cubic force constants, on the other hand, will allow us to compute phonon lifetimes and scattering rates. Results on equilibrium Green-Kubo molecular dynamics simulations of thermal conductivity as well as an alternative calculation of the latter based on the relaxation-time approximation will be reported. The accuracy and ease of computation of the lattice thermal conductivity using these methods will be compared. This approach paves the way for the construction of accurate bulk interatomic potentials database, from which lattice dynamics and thermal properties can be calculated and used in larger scale simulation methods such as Monte Carlo.


Author(s):  
Andrea N. Bootsma ◽  
Steven Wheeler

<div>Density functional theory (DFT) has emerged as a powerful tool for analyzing organic and organometallic systems and proved remarkably accurate in computing the small free energy differences that underpin many chemical phenomena (e.g. regio- and stereoselective reactions). We show that the lack of rotational invariance of popular DFT integration grids reveals large uncertainties in computed free energies for isomerizations, torsional barriers, and regio- and stereoselective reactions. The result is that predictions based on DFT-computed free energies for many systems can change qualitatively depending on molecular orientation. For example, for a metal-free propargylation of benzaldehyde, predicted enantioselectivities based on B97-D/def2-TZVP free energies using the popular (75,302) integration grid can vary from 62:38 to 99:1 by simply rotating the transition state structures. Relative free energies for the regiocontrolling transition state structures for an Ir-catalyzed C–H functionalization reaction computed using M06/6-31G(d,p)/LANL2DZ and the same grid can vary by more than 5 kcal mol–1, resulting in predicted regioselectivities that range anywhere from 14:86 to >99:1. Errors of these magnitudes occur for different functionals and basis sets, are widespread among modern applications of DFT, and can be reduced by using much denser integration grids than commonly employed.</div>


2014 ◽  
Vol 92 (7) ◽  
pp. 616-624 ◽  
Author(s):  
Zhichao Liu ◽  
Qiong Wu ◽  
Weihua Zhu ◽  
Heming Xiao

Density functional theory was used to study the structural, electronic, and optical properties of crystalline 4,10-dinitro-2,6,8,12-tetraoxa-4,10-diazatetracyclo-[5.5.0.05,903,11]-dodecane (TEX) under hydrostatic pressure. The results indicate that there is a displacive transition in TEX under compression that has never been found in experiments. As the pressure increases, the band gap gradually decreases but presents an abnormal increase at 61 GPa, called the structural transition; moreover, the gap reduction is more pronounced in the low-pressure range compared with the high-pressure range. An analysis of density of states shows that the electronic delocalization in TEX is enhanced gradually under the influence of pressure. The peaks of the imaginary parts of the dielectric functions, energy-loss function, and reflectivity may come mainly from the electron transitions between the oxygen 2p and nitrogen 2p states. The electron energy-loss function presents a blue shift under compression. TEX has relatively higher optical activity at high pressure than at ambient conditions.


Sign in / Sign up

Export Citation Format

Share Document