Search

In the previous paper of this series [JCTC 2020, 16, 3757], we presented a theoretical and algorithmic framework based on a localized representation of the occupied space that exploits the inherent sparsity in the real-space evaluation of the EXX interaction in finite-gap systems. This was accompanied by a detailed description of exx, a massively parallel hybrid MPI/OpenMP implementation of this approach in Quantum ESPRESSO that enables linear-scaling hybrid DFT based AIMD in the NVE/NVT ensembles of condensed-phase systems containing 500--1000 atoms (in fixed orthorhombic cells) with a wall time cost comparable to semi-local DFT. In this work, we extend exx to enable hybrid DFT based AIMD of large-scale condensed-phase systems with general and fluctuating cells in the NpH/NpT ensembles. Our theoretical extension includes an analytical derivation of the EXX contribution to the stress tensor for systems in general cells with a computational complexity that scales linearly with system size. The corresponding algorithmic extensions to exx include optimized routines that: (i) handle static/fluctuating cells with non-orthogonal lattice symmetries, (ii) solve Poisson's equation in general cells via an automated selection of the auxiliary grid directions in the Natan-Kronik representation of the discrete Laplacian operator, and (iii) evaluate the EXX contribution to the stress tensor. We also critically assess the computational performance of the extended exx module across several different HPC architectures via case studies on ice Ih, II, and III as well as ambient liquid water. We find that the extended exx can evaluate the EXX contribution to the stress tensor with negligible cost (< 1%) and remains highly scalable, thereby bringing us another step closer to routinely performing hybrid DFT based AIMD for large-scale condensed-phase systems across a wide range of thermodynamic conditions., 30 pages and 5 figures

We carry out X-ray absorption spectroscopy experiment at oxygen K-edge in croconic acid (C5H2O5) crystal as a prototype of ferroelectric organic molecular solid, whose electric polarization is generated by proton transfer. The experimental spectrum is well reproduced by the electron-hole excitation theory simulations from configuration generated by ab initio molecular dynamics simulation. When inversion symmetry is broken in ferroelectric state, the hydrogen bonding environment on the two bonded molecules become inequivalent. Such a difference is sensitively probed by the bound excitation in the pre-edge, which are strongly localized on the excited molecules. Our analysis shows that a satellite peak in the pre-edge will emerge at higher excitation energy which serves as a clear signature of ferroelectricity in the material., Comment: 6 pages, 3 figures

We explore how anharmonicity, nuclear quantum effects (NQE), many-body dispersion interactions, and Pauli repulsion influence thermal properties of dispersion-bound molecular crystals. Accounting for anharmonicity with ab initio molecular dynamics yields cell parameters accurate to within $2%$ of experiment for a set of pyridinelike molecular crystals at finite temperatures and pressures. From the experimental thermal expansion curve, we find that pyridine-I has a Debye temperature just above its melting point, indicating sizable NQE across the entire crystalline range of stability. We find that NQE lead to a substantial volume increase in pyridine-I $(\ensuremath{\approx}40$% more than classical thermal expansion at 153 K) and attribute this to intermolecular Pauli repulsion promoted by intramolecular quantum fluctuations. When predicting delicate properties such as the thermal expansivity, we show that many-body dispersion interactions and more sophisticated density functional approximations improve the accuracy of the theoretical model.

We introduce a local order metric (LOM) that measures the degree of order in the neighborhood of an atomic or molecular site in a condensed medium. The LOM maximizes the overlap between the spatial distribution of sites belonging to that neighborhood and the corresponding distribution in a suitable reference system. The LOM takes a value tending to zero for completely disordered environments and tending to one for environments that perfectly match the reference. The site-averaged LOM and its standard deviation define two scalar order parameters, $S$ and $\ensuremath{\delta}S$, that characterize with excellent resolution crystals, liquids, and amorphous materials. We show with molecular dynamics simulations that $S$, $\ensuremath{\delta}S$, and the LOM provide very insightful information in the study of structural transformations, such as those occurring when ice spontaneously nucleates from supercooled water or when a supercooled water sample becomes amorphous upon progressive cooling.

We perform ab initio molecular dynamics (AIMD) simulation of liquid water in the canonical ensemble at ambient conditions using the strongly constrained and appropriately normed (SCAN) meta-generalized-gradient approximation (GGA) functional approximation and carry out systematic comparisons with the results obtained from the GGA-level Perdew-Burke-Ernzerhof (PBE) functional and Tkatchenko-Scheffler van der Waals (vdW) dispersion correction inclusive PBE functional. We analyze various properties of liquid water including radial distribution functions, oxygen-oxygen-oxygen triplet angular distribution, tetrahedrality, hydrogen bonds, diffusion coefficients, ring statistics, density of states, band gaps, and dipole moments. We find that the SCAN functional is generally more accurate than the other two functionals for liquid water by not only capturing the intermediate-range vdW interactions but also mitigating the overly strong hydrogen bonds prescribed in PBE simulations. We also compare the results of SCAN-based AIMD simulations in the canonical and isothermal-isobaric ensembles. Our results suggest that SCAN provides a reliable description for most structural, electronic, and dynamical properties in liquid water.

Quantum ESPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudo-potential and projector-augmented-wave approaches. Quantum ESPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement theirs ideas. In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software., Comment: Psi-k highlight September 2017: psi-k.net/dowlnload/highlights/Highlight_137.pdf; J. Phys.: Condens. Matter, accepted

Significance Water is vital to our everyday life, but its structure at a molecular level is still not fully understood from either experiment or theory. The latter is hampered by our inability to construct a purely predictive, first principles model. The difficulty in modeling water lies in capturing the delicate interplay among the many strong and weak forces that govern its behavior and phase diagram. Herein, molecular simulations with a recently proposed nonempirical quantum mechanical approach (the SCAN density functional) yield an excellent description of the structural, electronic, and dynamic properties of liquid water. SCAN (strongly constrained and appropriately normed)-based approaches, which describe diverse types of bonds in materials on an equal, accurate footing, will likely enable efficient and reliable modeling of aqueous phase chemistry.

Metal string complexes or extended metal atom chains (EMACs) belong to a family of molecules that consist of a linear chain of directly bonded metal atoms embraced helically by four multidentate organic ligands. These four organic ligands are usually made up of repeating pyridyl units, single-nitrogen-substituted heterocyclic annulenes, bridged by independent amido groups. Here, in this paper, we show that these heterocyclic annulenes are actually nanoscale molecular split-ring resonators (SRRs) that can exhibit simultaneous negative electric permittivity and magnetic permeability in the UV-Vis region. Moreover, a monolayer of self-assembled EMACs is a periodic array of molecular SRRs which can be considered as a negative refractive index material. In the molecular scale, where the quantum-size effect is significant, we apply the tight-binding method to obtain the frequency-dependent permittivity and permeability of these molecular SRRs with their tensorial properties carefully considered., 8 pages, 8 figures