Role of nonlocal exchange in molecular crystals: The case of two proton-ordered phases of ice

Range-Separated Double-Hybrid Functional from Nonempirical Constraints

On the basis of our previous developments in the field of nonempirical double hybrids, we present here a new exchange-correlation functional based on a range-separated model for the exchange part and integrating a nonlocal perturbative correction to the electron correlation contribution. Named RSX-QIDH, the functional is free from any kind of empirical parametrization. Its range-separation parameter is set to recover the total energy of the hydrogen atom, thus eliminating the self-interaction error for this one-electron system. Subsequent tests on some relevant benchmark data sets confirm that the self-interaction error is particularly low for RSX-QIDH. This new functional provides also correct dissociation profiles for charged rare-gas dimers and very accurate ionization potentials directly from Kohn-Sham orbital energies. Above all, these good results are not obtained at the expense of other properties. Indeed, further tests on standard benchmarks show that RSX-QIDH is competitive with the more empirical ωB97X-2 double hybrid and outperforms the parent LC-PBE long-range corrected hybrid, thus underlining the important role of the nonlocal perturbative correlation.

Conceptual design of tetraazaporphyrin- and subtetraazaporphyrin-based functional nanocarbon materials: electronic structures, topologies, optical properties, and methane storage capacities

A large variety of conceptual three- and fourfold tetraazaporphyrin- and subtetraazaporphyrin-based functional 3D nanocage and nanobarrel structures have been proposed on the basis of in silico design. The designed structures differ in their sizes, topology, porosity, and conjugation properties. The stability of nanocages of Oh symmetry and nanobarrels of D4h symmetry was revealed on the basis of DFT and MD calculations, whereas their optical properties were assessed using a TDDFT approach and a long-range corrected LC-wPBE exchange-correlation functional. It was shown that the electronic structures and vertical excitation energies of the functional nanocage and nanobarrel structures could be easily tuned via their size, topology, and the presence of bridging sp(3) carbon atoms. TDDFT calculations suggest significantly lower excitation energies in fully conjugated nanocages and nanobarrels compared with systems with bridging sp(3) carbon fragments. Based on DFT and TDDFT calculations, the optical properties of the new materials can rival those of known quantum dots and are superior to those of monomeric phthalocyanines and their analogues. The methane gas adsorption properties of the new nanostructures and nanotubes generated by conversion from nanobarrels were studied using an MD simulation approach. The ability to store large quantities of methane (106-216 cm(3) (STP) cm(-3)) was observed in all cases with several compounds being close to or exceeding the DOE target of 180 cm(3) (STP) cm(-3) for material-based methane storage at a pressure of 3.5 MPa and room temperature.

Probing the structural and dynamical properties of liquid water with models including non-local electron correlation

Water is a ubiquitous liquid that displays a wide range of anomalous properties and has a delicate structure that challenges experiment and simulation alike. The various intermolecular interactions that play an important role, such as repulsion, polarization, hydrogen bonding, and van der Waals interactions, are often difficult to reproduce faithfully in atomistic models. Here, electronic structure theories including all these interactions at equal footing, which requires the inclusion of non-local electron correlation, are used to describe structure and dynamics of bulk liquid water. Isobaric-isothermal (NpT) ensemble simulations based on the Random Phase Approximation (RPA) yield excellent density (0.994 g/ml) and fair radial distribution functions, while various other density functional approximations produce scattered results (0.8-1.2 g/ml). Molecular dynamics simulation in the microcanonical (NVE) ensemble based on Møller-Plesset perturbation theory (MP2) yields dynamical properties in the condensed phase, namely, the infrared spectrum and diffusion constant. At the MP2 and RPA levels of theory, ice is correctly predicted to float on water, resolving one of the anomalies as resulting from a delicate balance between van der Waals and hydrogen bonding interactions. For several properties, obtaining quantitative agreement with experiment requires correction for nuclear quantum effects (NQEs), highlighting their importance, for structure, dynamics, and electronic properties. A computed NQE shift of 0.6 eV for the band gap and absorption spectrum illustrates the latter. Giving access to both structure and dynamics of condensed phase systems, non-local electron correlation will increasingly be used to study systems where weak interactions are of paramount importance.

Communication: On the stability of ice 0, ice i, and I(h)

Using ab initio methods, we examine the stability of ice 0, a recently proposed tetragonal form of ice implicated in the homogeneous freezing of water [J. Russo, F. Romano, and H. Tanaka, Nat. Mater. 13, 670 (2014)]. Vibrational frequencies are computed across the complete Brillouin Zone using Density Functional Theory (DFT), to confirm mechanical stability and quantify the free energy of ice 0 relative to ice I(h). The robustness of this result is tested via dispersion corrected semi-local and hybrid DFT, and Quantum Monte-Carlo calculation of lattice energies. Results indicate that popular molecular models only slightly overestimate the stability of ice zero. In addition, we study all possible realisations of proton disorder within the ice zero unit cell, and identify the ground state as ferroelectric. Comparisons are made to other low density metastable forms of ice, suggesting that the ice i structure [C. J. Fennel and J. D. Gezelter, J. Chem. Theory Comput. 1, 662 (2005)] may be equally relevant to ice formation.

The individual and collective effects of exact exchange and dispersion interactions on the ab initio structure of liquid water

In this work, we report the results of a series of density functional theory (DFT) based ab initio molecular dynamics (AIMD) simulations of ambient liquid water using a hierarchy of exchange-correlation (XC) functionals to investigate the individual and collective effects of exact exchange (Exx), via the PBE0 hybrid functional, non-local van der Waals/dispersion (vdW) interactions, via a fully self-consistent density-dependent dispersion correction, and an approximate treatment of nuclear quantum effects, via a 30 K increase in the simulation temperature, on the microscopic structure of liquid water. Based on these AIMD simulations, we found that the collective inclusion of Exx and vdW as resulting from a large-scale AIMD simulation of (H2O)128 significantly softens the structure of ambient liquid water and yields an oxygen-oxygen structure factor, SOO(Q), and corresponding oxygen-oxygen radial distribution function, gOO(r), that are now in quantitative agreement with the best available experimental data. This level of agreement between simulation and experiment demonstrated herein originates from an increase in the relative population of water molecules in the interstitial region between the first and second coordination shells, a collective reorganization in the liquid phase which is facilitated by a weakening of the hydrogen bond strength by the use of a hybrid XC functional, coupled with a relative stabilization of the resultant disordered liquid water configurations by the inclusion of non-local vdW/dispersion interactions. This increasingly more accurate description of the underlying hydrogen bond network in liquid water also yields higher-order correlation functions, such as the oxygen-oxygen-oxygen triplet angular distribution, POOO(θ), and therefore the degree of local tetrahedrality, as well as electrostatic properties, such as the effective molecular dipole moment, that are in much better agreement with experiment.

Periodic MP2, RPA, and Boundary Condition Assessment of Hydrogen Ordering in Ice XV

Ice XV is the hydrogen-ordered form of the ice VI phase whose structure was predicted to be Cc and ferroelectric using periodic DFT approaches. However, neutron diffraction and Raman spectroscopy data show the structure to have P1̅ symmetry and to be antiferroelectric. Recent work1 using fragment-based MP2 and CCSD(T) approaches predicts the experimental structure as the ground state. We have reconsidered this problem using fully periodic MP2 and RPA approaches and find that the ferroelectric Cc structure is the lowest energy configuration. However, ubiquitously employed tinfoil boundary conditions stabilize polar structures. We suggest that ferroelectric Cc crystals can grow within a paraelectric ice VI matrix but may become unstable once a fraction of the matrix has become hydrogen-ordered. The reduction in dielectric constant causes P1̅ and other structures with small polarization to become favored, providing a possible resolution between observation and theoretical predictions.

Hydrogen bonds and van der waals forces in ice at ambient and high pressures

The first principles methods, density-functional theory and quantum Monte Carlo, have been used to examine the balance between van der Waals (vdW) forces and hydrogen bonding in ambient and high-pressure phases of ice. At higher pressure, the contribution to the lattice energy from vdW increases and that from hydrogen bonding decreases, leading vdW to have a substantial effect on the transition pressures between the crystalline ice phases. An important consequence, likely to be of relevance to molecular crystals in general, is that transition pressures obtained from density-functional theory exchange-correlation functionals which neglect vdW forces are greatly overestimated.