Can we predict the structure and stability of molecular crystals under increased pressure? First‐principles study of glycine phase transitions

Abstract

The aim of this study was to determine whether the periodic density functional theory (DFT) calculations can be used for accurate prediction of the influence of the increased pressure on crystal structure and stability of molecular solids. To achieve this goal a series of geometry optimization and thermodynamic parameters calculations were performed for γ-glycine and δ-glycine structures at different pressure values using CASTEP program. In order to perform most accurate geometry optimization various exchange-correlation functionals defined within generalized gradient approximation (GGA): PBE, PW91, RPBE, WC, PBESOL as well as defined within local density approximation (LDA), i.e. CAPZ, were tested. Geometry optimization was carried out using different dispersion correction methods (i.e. Grimme, TS, OBS) or without them. The linear response density functional perturbation theory (DFPT) was used to obtain the phonon dispersion curves and phonon density of states from which thermodynamic parameters, such as: free energy (ΔG), enthalpy (ΔH) and entropy (ΔS) were evaluated. The results of the geometry optimization depend strongly on the choice of the DFT functional. Calculated differences between the free energy of the studied polymorphic forms at the studied pressure values were consistent with experimental observations on their stability. The computations of thermodynamic properties not only confirmed the order of stability of two studied forms, but also enabled to predict the pressure at which this order is reversed. The results obtained in this study have proven that the plane-wave basis set first principles calculations under periodic conditions is suitable for accurate prediction of crystal structure and stability. © 2018 Wiley Periodicals, Inc.

How does the NMR thermometer work? Application of combined quantum molecular dynamics and GIPAW calculations into the study of lead nitrate

“…Finally, the PBESOL was especially created for the densely packed solids, like lead nitrate. This functional was recently proven to accurately predict the unit cell dimensions under increased pressure as well …”

Structural, elastic, mechanical and thermodynamic properties of HfB 4 under high pressure

“…Pseudo-atom calculations are performed for B: 2s 2 2p 1 and Hf: 5p 6 5d 2 6s 2 . To compare the performance of different approximations of exchange–correlation interaction, here we used the generalized gradient approximation (GGA) with the Perdew–Burke–Ernzerhof (PBE) functional [ 41 ] and the modified PBE functional (PBESOL) [ 42 ], which was proved to provide good results for solids of high density [ 43 ], as well as the local density approximation (LDA) with the form of Ceperley–Adler parametrized by Perdew & Zunger (CAPZ) [ 44 ], used to describe the exchange–correlation potentials. The structures were relaxed using the Broyden, Fletcher, Goldfarb and Shannon (BFGS) minimization method algorithm [ 45 ].…”

Can We Predict the Pressure Induced Phase Transition of Urea? Application of Quantum Molecular Dynamics

“…Fortunately, calculations on molecular crystals using density functional theory (DFT) based programs that enable to include the periodic boundary conditions of a studied system and a planewave basis set, such as CASTEP [9], have proven to be very accurate. Still, in most of the reported studies on the relative stability of the polymorphic forms solely the lattice energies [10][11][12] or, in less number of cases, free energy differences [13,14] of the structures are being calculated and compared. While those computational studies are indeed very interesting and appreciated, since their results enable the insight into the structure and stability of polymorphic forms, accurate phase transition modeling-defined here as a possibility to predict the changes in the crystal structure when exposed to the high pressure-is a much more complicated and challenging task.…”