Path-integral simulation of solids

Nuclear Quantum Effects on the Electronic Structure of Water and Ice

The electronic properties and optical response of ice and water are intricately shaped by their molecular structure, including the quantum mechanical nature of the hydrogen atoms. Despite numerous previous studies, a comprehensive understanding of the nuclear quantum effects (NQEs) on the electronic structure of water and ice at finite temperatures remains elusive. Here, we utilize molecular simulations that harness efficient machine-learning potentials and many-body perturbation theory to assess how NQEs impact the electronic bands of water and hexagonal ice. By comparing path-integral and classical simulations, we find that NQEs lead to a larger renormalization of the fundamental gap of ice, compared to that of water, ultimately yielding similar bandgaps in the two systems, consistent with experimental estimates. Our calculations suggest that the increased quantum mechanical delocalization of protons in ice, relative to water, is a key factor leading to the enhancement of NQEs on the electronic structure of ice.

A neural network-based four-body potential energy surface for parahydrogen

We present an isotropic ab initio (para-H2)4 four-body interaction potential energy surface (PES). The electronic structure calculations are performed at the correlated coupled-cluster theory level, with single, double, and perturbative triple excitations. They use an atom-centered augmented correlation-consistent double zeta basis set, supplemented by a (3s3p2d) midbond function. We use a multilayer perceptron to construct the PES. We apply a rescaling transformation to the output energies during training to improve the prediction of weaker energies in the sample data. At long distances, the interaction energies are adjusted to match the empirically derived four-body dispersion interaction. The four-body interaction energy at short intermolecular separations is net repulsive. The use of this four-body PES, in combination with a first principles pair potential for para-H2 [J. Chem. Phys. 119, 12551 (2015)] and an isotropic ab initio three-body potential for para-H2 [J. Chem. Phys. 156, 044301 (2022)], is expected to provide closer agreement with experimental results.

First-Principles Approach for Coupled Quantum Dynamics of Electrons and Protons in Heterogeneous Systems

The coupled quantum dynamics of electrons and protons is ubiquitous in many dynamical processes involving light-matter interaction, such as solar energy conversion in chemical systems and photosynthesis. A first-principles description of such nuclear-electronic quantum dynamics requires not only the time-dependent treatment of nonequilibrium electron dynamics but also that of quantum protons. Quantum mechanical correlation between electrons and protons adds further complexity to such coupled dynamics. Here we extend real-time nuclear-electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) to periodic systems and perform first-principles simulations of coupled quantum dynamics of electrons and protons in complex heterogeneous systems. The process studied is an electronically excited-state intramolecular proton transfer of o-hydroxybenzaldehyde in water and at a silicon (111) semiconductor-molecule interface. These simulations illustrate how environments such as hydrogen-bonding water molecules and an extended material surface impact the dynamical process on the atomistic level. Depending on how the molecule is chemisorbed on the surface, excited-state electron transfer from the molecule to the semiconductor surface can inhibit ultrafast proton transfer within the molecule. This Letter elucidates how heterogeneous environments influence the balance between the quantum mechanical proton transfer and excited electron dynamics. The periodic RT-NEO-TDDFT approach is applicable to a wide range of other photoinduced heterogeneous processes.

Quantum Vibronic Effects on the Electronic Properties of Molecular Crystals

We present a study of molecular crystals, focused on the effect of nuclear quantum motion and anharmonicity on their electronic properties. We consider a system composed of relatively rigid molecules, a diamondoid crystal, and one composed of floppier molecules, NAI-DMAC, a thermally activated delayed fluorescence compound. We compute fundamental electronic gaps at the density functional theory (DFT) level of theory, with the Perdew-Burke-Erzenhof (PBE) and strongly constrained and approximately normed (SCAN) functionals, by coupling first-principles molecular dynamics with a nuclear quantum thermostat. We find a sizable zero-point renormalization (ZPR) of the band gaps, which is much larger in the case of diamondoids (0.6 eV) than for NAI-DMAC (0.22 eV). We show that the frozen phonon (FP) approximation, which neglects intermolecular anharmonic effects, leads to a large error (∼50%) in the calculation of the band gap ZPR. Instead, when using a stochastic method, we obtain results in good agreement with those of our quantum simulations for the diamondoid crystal. However, the agreement is worse for NAI-DMAC where intramolecular anharmonicities contribute to the ZPR. Our results highlight the importance of accurately including nuclear and anharmonic quantum effects to predict the electronic properties of molecular crystals.

A Glimpse into Quantum Triplet Structures in Supercritical 3He

A methodological study of triplet structures in quantum matter is presented. The focus is on helium-3 under supercritical conditions (4 < T/K < 9; 0.022 < ρN/Å-3 < 0.028), for which strong quantum diffraction effects dominate the behavior. Computational results for the triplet instantaneous structures are reported. Path integral Monte Carlo (PIMC) and several closures are utilized to obtain structure information in the real and the Fourier spaces. PIMC involves the fourth-order propagator and the SAPT2 pair interaction potential. The main triplet closures are: AV3, built as the average of the Kirkwood superposition and the Jackson-Feenberg convolution, and the Barrat-Hansen-Pastore variational approach. The results illustrate the main characteristics of the procedures employed by concentrating on the salient equilateral and isosceles features of the computed structures. Finally, the valuable interpretive role of closures in the triplet context is highlighted.

Incorporating Nuclear Quantum Effects in Molecular Dynamics with a Constrained Minimized Energy Surface

The accurate incorporation of nuclear quantum effects in large-scale molecular dynamics (MD) simulations remains a significant challenge. Recently, we combined constrained nuclear-electronic orbital (CNEO) theory with classical MD and obtained a new approach (CNEO-MD) that can accurately and efficiently incorporate nuclear quantum effects into classical simulations. In this Letter, we provide the theoretical foundation for CNEO-MD by developing an alternative formulation of the equations of motion for MD. In this new formulation, the expectation values of quantum nuclear positions evolve classically on an effective energy surface that is obtained from a constrained energy minimization procedure when solving for the quantum nuclear wave function, thus enabling the incorporation of nuclear quantum effects in classical MD simulations. For comparison with other existing approaches, we examined a series of model systems and found that this new MD approach is significantly more accurate than the conventional way of performing classical MD and generally outperforms centroid MD and ring-polymer MD in describing vibrations in these model systems.

Influence of nuclear quantum effects on the electronic properties of amorphous carbon

In crystalline solids, atoms are arranged in periodic patterns on regular lattices. Amorphous solids, instead, lack long-range order; namely, a regular array of atoms beyond first or second nearest neighbors is absent. The lack of periodicity influences many properties of amorphous materials, including the coupling of electronic and nuclear motion. Here we study amorphous carbon, a system composed of a relatively light atom. We show that to understand its electronic properties, a quantum mechanical treatment of electron–nuclear coupling is essential, and we illustrate a simulation framework based on first principles to do so. We also discuss the role of specific defect states in the disordered network in determining the physical properties of amorphous carbon.

We carry out quantum simulations to study the physical properties of diamond-like amorphous carbon by coupling first-principles molecular dynamics with a quantum thermostat, and we analyze multiple samples representative of different defective sites present in the disordered network. We show that quantum vibronic coupling is critical in determining the electronic properties of the system, in particular its electronic and mobility gaps, while it has a moderate influence on the structural properties. We find that despite localized electronic states near the Fermi level, the quantum nature of the nuclear motion leads to a renormalization of the electronic gap surprisingly similar to that found in crystalline diamond. We also discuss the notable influence of nuclear quantum effects on band-like and variable-hopping mechanisms contributing to electrical conduction. Our calculations indicate that methods often used to evaluate electron–phonon coupling in ordered solids are inaccurate to study the electronic and transport properties of amorphous semiconductors composed of light atoms.

Molecular Dynamics with Constrained Nuclear Electronic Orbital Density Functional Theory: Accurate Vibrational Spectra from Efficient Incorporation of Nuclear Quantum Effects

Nuclear quantum effects play a crucial role in many chemical and biological systems involving hydrogen atoms yet are difficult to include in practical molecular simulations. In this paper, we combine our recently developed methods of constrained nuclear-electronic orbital density functional theory (cNEO-DFT) and constrained minimized energy surface molecular dynamics (CMES-MD) to create a new method for accurately and efficiently describing nuclear quantum effects in molecular simulations. By use of this new method, dubbed cNEO-MD, the vibrational spectra of a set of small molecules are calculated and compared with those from conventional ab initio molecular dynamics (AIMD) as well as from experiments. With the same formal scaling, cNEO-MD greatly outperforms AIMD in describing the vibrational modes with significant hydrogen motion characters, demonstrating the promise of cNEO-MD for simulating chemical and biological systems with significant nuclear quantum effects.

Renormalization group analysis of weakly interacting van der Waals Fermi system

Weak-coupling phenomena of the two-dimensional Hubbard model is gaining momentum as a new interesting research field due to its extraordinarily rich behavior as a function of the carrier density and model parameters. Salmhofer (1998Commun.Math.Phys.194249; 2001Phys.Rev.Lett.87187004) developed a new renormalization-group method for interacting Fermi systems and Metzner (2000Phys.Rev. B617364; 2000Phys.Rev.Lett.855162) implemented this renormalization group analysis of the two-dimensional Hubbard model. In this work, we demonstrate the spin-wave dependent electronic structure and susceptibility behavior of model graphene-phosphorene van der Waals heterostructure in the framework of renormalization group approach. We implement singlet vertex response function for the weakly interacting van der Waals Fermi system with nearest-neighbor hopping amplitudes. This analytical approach is further extended for spin-wave dependent susceptibility behavior. We present the resulting compressibility and phase diagram in the vicinity of half-filling, and also results for the density dependence of the critical energy scale.

Real Space Triplets in Quantum Condensed Matter: Numerical Experiments Using Path Integrals, Closures, and Hard Spheres

Path integral Monte Carlo and closure computations are utilized to study real space triplet correlations in the quantum hard-sphere system. The conditions cover from the normal fluid phase to the solid phases face-centered cubic (FCC) and cI16 (de Broglie wavelengths , densities ). The focus is on the equilateral and isosceles features of the path-integral centroid and instantaneous structures. Complementary calculations of the associated pair structures are also carried out to strengthen structural identifications and facilitate closure evaluations. The three closures employed are Kirkwood superposition, Jackson-Feenberg convolution, and their average (AV3). A large quantity of new data are reported, and conclusions are drawn regarding (i) the remarkable performance of AV3 for the centroid and instantaneous correlations, (ii) the correspondences between the fluid and FCC salient features on the coexistence line, and (iii) the most conspicuous differences between FCC and cI16 at the pair and the triplet levels at moderately high densities (. This research is expected to provide low-temperature insights useful for the future related studies of properties of real systems (e.g., helium, alkali metals, and general colloidal systems).

Phase behavior of the quantum Lennard-Jones solid

The Lennard-Jones (LJ) potential is perhaps one of the most widely used models for the interaction of uncharged particles, such as noble gas solids. The phase diagram of the classical LJ solid is known to exhibit transitions between hcp and fcc phases. However, the phase behavior of the quantum LJ solid remains unknown. Thermodynamic integration based on path integral molecular dynamics (PIMD) and lattice dynamics calculations are used to study the phase stability of the hcp and fcc LJ solids. The hcp phase is shown to be stabilized by quantum effects in PIMD, while fcc is shown to be favored by lattice dynamics, which suggests a possible re-entrant low pressure fcc phase for highly quantum systems. Implications for the phase stability of noble gas solids are discussed. For parameters equating to helium, the expansion due to zero-point vibrations is associated with quantum melting: neither crystal structure is stable at zero pressure.

Nuclear quantum effect with pure anharmonicity and the anomalous thermal expansion of silicon

Despite the widespread use of silicon in modern technology, its peculiar thermal expansion is not well understood. Adapting harmonic phonons to the specific volume at temperature, the quasiharmonic approximation, has become accepted for simulating the thermal expansion, but has given ambiguous interpretations for microscopic mechanisms. To test atomistic mechanisms, we performed inelastic neutron scattering experiments from 100 K to 1,500 K on a single crystal of silicon to measure the changes in phonon frequencies. Our state-of-the-art ab initio calculations, which fully account for phonon anharmonicity and nuclear quantum effects, reproduced the measured shifts of individual phonons with temperature, whereas quasiharmonic shifts were mostly of the wrong sign. Surprisingly, the accepted quasiharmonic model was found to predict the thermal expansion owing to a large cancellation of contributions from individual phonons.

Inverse Temperature Dependence of Nuclear Quantum Effects in DNA Base Pairs

Despite the inherently quantum mechanical nature of hydrogen bonding, it is unclear how nuclear quantum effects (NQEs) alter the strengths of hydrogen bonds. With this in mind, we use ab initio path integral molecular dynamics to determine the absolute contribution of NQEs to the binding in DNA base pair complexes, arguably the most important hydrogen-bonded systems of all. We find that depending on the temperature, NQEs can either strengthen or weaken the binding within the hydrogen-bonded complexes. As a somewhat counterintuitive consequence, NQEs can have a smaller impact on hydrogen bond strengths at cryogenic temperatures than at room temperature. We rationalize this in terms of a competition of NQEs between low-frequency and high-frequency vibrational modes. Extending this idea, we also propose a simple model to predict the temperature dependence of NQEs on hydrogen bond strengths in general.

Very fast averaging of thermal properties of crystals by molecular simulation

Knowledge of approximate harmonic behavior of crystals is introduced into a new "mapped averaging" framework to yield alternative expressions for the thermodynamic properties of crystalline systems. The expressions separate the known harmonic behavior from residual averages, which thus encapsulate anharmonic contributions to the properties. With harmonic contributions removed, direct measurement of these anharmonic contributions by molecular simulation can be accomplished without contamination by noise produced by the already-known harmonic behavior. We show with application to the Lennard-Jones model that first-derivative properties (pressure, energy) can be obtained to a given precision via this harmonically mapped averaging at least 10 times faster than by using conventional averaging, and second-derivative properties (e.g., heat capacity) are obtained at least 100 times faster; in more favorable cases, the speedup exceeds a millionfold. Free-energy calculations are accelerated by 50 to 1000 times. Data obtained using these formulations are rigorous and not subject to any added approximation, and in fact are less sensitive to inaccuracies relating to finite-size effects, potential truncation, equilibration, and similar considerations. Moreover, the approach does not require any alteration in how sampling is performed during the simulation, so it may be used with standard Monte Carlo or molecular dynamics methods. However, the mapped averages do require evaluation of first and second derivatives of the intermolecular potential, for evaluation of first and second thermodynamic-derivative properties, respectively. Apart from its usefulness to simulation, the formalism developed here may constitute a basis for new theoretical treatments of crystals.