Stress theorem in the determination of static equilibrium by the density functional method

Reversible electric-field-induced phase transition in Ca-modified NaNbO3 perovskites for energy storage applications

Sodium niobate (NaNbO₃) is a potential material for lead-free dielectric ceramic capacitors for energy storage applications because of its antipolar ordering. In principle, a reversible phase transition between antiferroelectric (AFE) and ferroelectric (FE) phases can be induced by an application of electric field (E) and provides a large recoverable energy density. However, an irreversible phase transition from the AFE to the FE phase usually takes place and an AFE-derived polarization feature, a double polarization (P)-E hysteresis loop, does not appear. In this study, we investigate the impact of chemically induced hydrostatic pressure (p_chem) on the phase stability and polarization characteristics of NaNbO₃-based ceramics. We reveal that the cell volume of Ca-modified NaNbO₃ [(Ca_xNa_1-2xV_x)NbO₃], where V is A-site vacancy, decreases with increasing x by a positive p_chem. Structural analysis using micro-X-ray diffraction measurements shows that a reversible AFE-FE phase transition leads to a double P-E hysteresis loop for the sample with x = 0.10. DFT calculations support that a positive p_chem stabilizes the AFE phase even after the electrical poling and provides the reversible phase transition. Our study demonstrates that an application of positive p_chem is effective in delivering the double P-E loop in the NaNbO₃ system for energy storage applications.

First-Principles Calculation of Third-Order Elastic Constants via Numerical Differentiation of the Second Piola-Kirchhoff Stress Tensor

A general method is presented to calculate from first principles the full set of third-order elastic constants of a material of arbitrary symmetry. The method here illustrated relies on a plane-wave density functional theory scheme to calculate the Cauchy stress and the numerical differentiation of the second Piola-Kirchhoff stress tensor to evaluate the elastic constants. It is shown that finite difference formulas lead to a cancellation of the finite basis set errors, whereas simple solutions are proposed to eliminate numerical errors arising from the use of Fourier interpolation techniques. Applications to diamond, silicon, aluminum, magnesium, graphene, and a graphane conformer give results in excellent agreement with both experiments and previous calculations based on fitting energy density curves, demonstrating both the accuracy and generality of our new methodology to investigate nonlinear elastic behaviors of materials.

The phase diagram of high-pressure superionic ice

Superionic ice is a special group of ice phases at high temperature and pressure, which may exist in ice-rich planets and exoplanets. In superionic ice liquid hydrogen coexists with a crystalline oxygen sublattice. At high pressures, the properties of superionic ice are largely unknown. Here we report evidence that from 280 GPa to 1.3 TPa, there are several competing phases within the close-packed oxygen sublattice. At even higher pressure, the close-packed structure of the oxygen sublattice becomes unstable to a new unusual superionic phase in which the oxygen sublattice takes the P2₁/c symmetry. We also discover that higher pressure phases have lower transition temperatures. The diffusive hydrogen in the P2₁/c superionic phase shows strong anisotropic behaviour and forms a quasi-two-dimensional liquid. The ionic conductivity changes abruptly in the solid to close-packed superionic phase transition, but continuously in the solid to P2₁/c superionic phase transition.

At high pressure, water forms superionic ice with an oxygen lattice and melted liquid hydrogens, which could exist on ice-rich planets. Here, Sun et al. predict a new phase of superionic ice where the hydrogens preferentially diffuse in two-dimensions within oxygen superlattice with the P2₁/c symmetry.

Fine-tuning the theoretically predicted structure of MIL-47(V) with the aid of powder X-ray diffraction

The structural characterization of complex crystalline materials can be simplified by closely comparing theoretical and experimental diffraction patterns.

Quasi-1D physics in metal-organic frameworks: MIL-47(V) from first principles

The geometric and electronic structure of the MIL-47(V) metal-organic framework (MOF) is investigated by using ab initio density functional theory (DFT) calculations. Special focus is placed on the relation between the spin configuration and the properties of the MOF. The ground state is found to be antiferromagnetic, with an equilibrium volume of 1554.70 Å(3). The transition pressure of the pressure-induced large-pore-to-narrow-pore phase transition is calculated to be 82 MPa and 124 MPa for systems with ferromagnetic and antiferromagnetic chains, respectively. For a mixed system, the transition pressure is found to be a weighted average of the ferromagnetic and antiferromagnetic transition pressures. Mapping DFT energies onto a simple-spin Hamiltonian shows both the intra- and inter-chain coupling to be antiferromagnetic, with the latter coupling constant being two orders of magnitude smaller than the former, suggesting the MIL-47(V) to present quasi-1D behavior. The electronic structure of the different spin configurations is investigated and it shows that the band gap position varies strongly with the spin configuration. The valence and conduction bands show a clear V d-character. In addition, these bands are flat in directions orthogonal to VO6 chains, while showing dispersion along the the direction of the VO6 chains, similar as for other quasi-1D materials.

Role of Dispersive Interactions in Determining Structural Properties of Organic-Inorganic Halide Perovskites: Insights from First-Principles Calculations

A microscopic picture of structure and bonding in organic-inorganic perovskites is imperative to understanding their remarkable semiconducting and photovoltaic properties. On the basis of a density functional theory treatment that includes both spin-orbit coupling and dispersive interactions, we provide detailed insight into the crystal binding of lead-halide perovskites and quantify the effect of different types of interactions on the structural properties. Our analysis reveals that cohesion in these materials is characterized by a variety of interactions that includes important contributions from both van der Waals interactions among the halide atoms and hydrogen bonding. We also assess the role of spin-orbit coupling and show that it causes slight changes in lead-halide bonding that do not significantly affect the lattice parameters. Our results establish that consideration of dispersive effects is essential for understanding the structure and bonding in organic-inorganic perovskites in general and for providing reliable theoretical predictions of structural parameters in particular.

Melting of cubic boron nitride at extreme pressures

Due to its large pressure range of stability and inert nature, cubic boron nitride has been proposed as a potential pressure standard for high pressure experiments. It is extremely refractive upon compression, although its melting temperature is not known beyond 10 GPa. We apply first-principles molecular dynamics to evaluate the thermodynamics of zincblende structured (cubic) and liquid boron nitride at extreme temperatures and pressures, and compute the melting curve up to 1 TPa by integration of the Clapeyron equation. The resulting equations of state reveal that liquid boron nitride becomes denser than the solid phase at pressures of around 0.5 TPa. This is expressed as a turnover in the melting curve, which reaches a maximum at 510 GPa and 6550 ± 700 K. The origin of this density crossover is explained in terms of the underlying liquid structure, which diverges from that of the zincblende structured solid as the phases are compressed.

First Principles Calculations of Surface Stress

Only recently have there been fully quantum-mechanical calculations of two-dimensional surface stress tensors. We have calculated total energies and stresses of semiconductor surfaces within the Local Density Approximation, using norm-conserving pseudopotentials. In order to hasten convergence of the stress with respect to basis size, it is useful to remove a fictitious tensile stress. We have calculated surface stress for the relaxed Si (111) 1×1 and 2×2-adatom surfaces, as well as for the relaxed Ge (111) 1×1 and 2×2-adatom surfaces. We have also calculated the surface stress for several chemisorbed systems, including Ga, Ge and As chemisorbed onto Si. We find a dramatic correlation between the electronic structure and chemistry of the surface, and its elastic properties.

Solid-State Phase Transition Induced by Pressure in LiOH·H2O

When the free energy surface of the lithium hydroxide monohydrate crystal was explored, the high-pressure solid-state phase transition was determined. The high-pressure phase has been obtained through ab initio Car-Parrinello molecular dynamics simulation in the isothermic-isobaric ensemble. The recent metadynamics method has been applied to overcome the high activation energy barriers typical of rare events, like solid-state phase transition at high pressures. In the LiOH x H2O system, there are two kinds of H bonds: water-water and hydroxyl-water. The effect of the pressure has been investigated, to give further insight into the high-pressure phase. The strengthening of the H bonds of the system produces modifications in the water and the hydroxyl ion dipole electronic environment. The infrared spectra of both phases have been calculated and compared with experiments, and the assignment of the external modes has been discussed.

Isobaric-Isothermal Monte Carlo Simulations from First Principles: Application to Liquid Water at Ambient Conditions

A series of first-principles Monte Carlo simulations in the isobaric-isothermal ensemble were carried out for liquid water at ambient conditions (T=298 K and p=1 atm). The Becke-Lee-Yang-Parr (BLYP) exchange and correlation energy functionals and norm-conserving Goedecker-Teter-Hutter (GTH) pseudopotentials were employed with the CP2 K simulation package to examine systems consisting of 64 water molecules. The fluctuations in the system volume encountered in simulations in the isobaric-isothermal ensemble require a reconsideration of the suitability of the typical charge-density cutoff and the regular grid-generation method previously used for the computation of the electrostatic energy in first-principles simulations in the microcanonical or canonical ensembles. In particular, it is noted that a much higher cutoff is needed and that the most computationally efficient method of creating grids can result in poor simulations. Analysis of the simulation trajectories using a very large charge-density cutoff at 1200 Ry and four different grid-generation methods point to a significantly underestimated liquid density of about 0.8 g cm-3 resulting in a somewhat understructured liquid (with a value of about 2.7 for the height of the first peak in the oxygen-oxygen radial distribution function) for BLYP-GTH water at ambient conditions. In addition, a simulation using a charge-density cutoff at 280 Ry yields a higher density of 0.9 g cm-3, showing the sensitivity of the simulation outcome to this parameter.

Liquid-liquid phase transitions of phosphorus via constant-pressure first-principles molecular dynamics simulations

Pressure-induced phase transitions in liquid phosphorus have been studied by constant-pressure first-principles molecular dynamics simulations. By compressing a low-pressure liquid which consists of the tetrahedral P4 molecules, a structural phase transition from the molecular to polymeric liquid (a high-pressure phase) observed in the recent experiment by Katayama et al. [Nature (London) 403, 170 (2000)] was successfully realized. It is found that this transition is caused by a breakup of the tetrahedral molecules with large volume contraction. The same transition is also realized by heating. This indicates that only the polymeric liquid can stably exist at high temperature.