Quantum Proton Effects from Density Matrix Renormalization Group Calculations

Nuclear-Electronic Orbital Multireference Configuration Interaction for Ground and Excited Vibronic States and Fundamental Insights into Multicomponent Basis Sets

The nuclear-electronic orbital (NEO) approach incorporates nuclear quantum effects into quantum chemistry calculations by treating specified nuclei quantum mechanically, equivalently to the electrons. Within the NEO framework, excited states are vibronic states representing electronic and nuclear vibrational excitations. The NEO multireference configuration interaction (MRCI) method presented herein provides accurate ground and excited vibronic states. The electronic and nuclear orbitals are optimized with a NEO multiconfigurational self-consistent field (NEO-MCSCF) procedure, thereby including both static and dynamic correlation and allowing the description of double and higher excitations. The accuracy of the NEO-MRCI method is illustrated by computing the ground state protonic densities and excitation energies of the vibronic states for four molecular systems with the hydrogen nucleus treated quantum mechanically. In addition, revised conventional electronic basis sets adapted for quantized nuclei are developed and shown to be essential for achieving this level of accuracy. The NEO-MRCI approach, as well as the strategy for revising electronic basis sets, will play a critical role in multicomponent quantum chemistry.

Restoring rotational symmetry of multicomponent wavefunctions with nuclear orbitals

In this work, we present a non-orthogonal configuration interaction (NOCI) approach to address the rotational corrections in multicomponent quantum chemistry calculations where hydrogen nuclei and electrons are described with orbitals under Hartree-Fock (HF) and density functional theory (DFT) frameworks. The rotational corrections are required in systems such as diatomic (HX) and nonlinear triatomic molecules (HXY), where localized broken-symmetry nuclear orbitals have a lower energy than delocalized orbitals with the correct symmetry. By restoring rotational symmetry with the proposed NOCI approach, we demonstrate significant improvements in proton binding energy predictions at the HF level, with average rotational corrections of 0.46 eV for HX and 0.23 eV for HXY molecules. For computing rotational excitation energies, our results indicate that HF kinetic energy corrections are consistently accurate, while discrepancies arise in total energy predictions, primarily from an incomplete treatment of dynamical correlation effects. Rotational energy corrections in multicomponent DFT calculations, using the epc17-2 proton-electron correlation functional, lead to an overestimation of proton binding energies. This is as a result of double-counting of proton-electron correlation effects in the off-diagonal NOCI terms. As a correction, we propose a scaling scheme that effectively adjusts the proton-electron correlation contributions, bringing our results into close agreement with reference CCSD(T) data. The scaled rotational corrections, on average, increase the epc17-2 proton binding energy predictions by 0.055 eV for HX and 0.025 eV for HXY and yield average deviations of 1.0 cm-1 for rotational transitions.

Polaritonic Chemistry Using the Density Matrix Renormalization Group Method

The emerging field of polaritonic chemistry explores the behavior of molecules under strong coupling with cavity modes. Despite recent developments in ab initio polaritonic methods for simulating polaritonic chemistry under electronic strong coupling, their capabilities are limited, especially in cases where the molecule also features strong electronic correlation. To bridge this gap, we have developed a novel method for cavity QED calculations utilizing the Density Matrix Renormalization Group (DMRG) algorithm in conjunction with the Pauli-Fierz Hamiltonian. Our approach is applied to investigate the effect of the cavity on the S0-S1 transition of n-oligoacenes, with n ranging from 2 to 5, encompassing 22 fully correlated π orbitals in the largest pentacene molecule. Our findings indicate that the influence of the cavity intensifies with larger acenes. Additionally, we demonstrate that, unlike the full determinantal representation, DMRG efficiently optimizes and eliminates excess photonic degrees of freedom, resulting in an asymptotically constant computational cost as the photonic basis increases.

SCINE-Software for chemical interaction networks

The software for chemical interaction networks (SCINE) project aims at pushing the frontier of quantum chemical calculations on molecular structures to a new level. While calculations on individual structures as well as on simple relations between them have become routine in chemistry, new developments have pushed the frontier in the field to high-throughput calculations. Chemical relations may be created by a search for specific molecular properties in a molecular design attempt, or they can be defined by a set of elementary reaction steps that form a chemical reaction network. The software modules of SCINE have been designed to facilitate such studies. The features of the modules are (i) general applicability of the applied methodologies ranging from electronic structure (no restriction to specific elements of the periodic table) to microkinetic modeling (with little restrictions on molecularity), full modularity so that SCINE modules can also be applied as stand-alone programs or be exchanged for external software packages that fulfill a similar purpose (to increase options for computational campaigns and to provide alternatives in case of tasks that are hard or impossible to accomplish with certain programs), (ii) high stability and autonomous operations so that control and steering by an operator are as easy as possible, and (iii) easy embedding into complex heterogeneous environments for molecular structures taken individually or in the context of a reaction network. A graphical user interface unites all modules and ensures interoperability. All components of the software have been made available as open source and free of charge.

Nuclear-Electronic Orbital Time-Dependent Configuration Interaction Method

Combining real-time electronic structure with the nuclear-electronic orbital (NEO) method has enabled the simulation of complex nonadiabatic chemical processes. However, accurate descriptions of hydrogen tunneling and double excitations require multiconfigurational treatments. Herein, we develop and implement the real-time NEO time-dependent configuration interaction (NEO-TDCI) approach. Comparison to NEO-full CI calculations of absorption spectra for a molecular system shows that the NEO-TDCI approach can accurately capture the tunneling splitting associated with the electronic ground state as well as vibronic progressions corresponding to double electron-proton excitations associated with excited electronic states. Both of these features are absent from spectra obtained with single reference real-time NEO methods. Our simulations of hydrogen tunneling dynamics illustrate the oscillation of the proton density from one side to the other via a delocalized, bilobal proton wave function. These results indicate that the NEO-TDCI approach is highly suitable for studying hydrogen tunneling and other inherently multiconfigurational systems.

Toward Accurate Post-Born-Oppenheimer Molecular Simulations on Quantum Computers: An Adaptive Variational Eigensolver with Nuclear-Electronic Frozen Natural Orbitals

Nuclear quantum effects such as zero-point energy and hydrogen tunneling play a central role in many biological and chemical processes. The nuclear-electronic orbital (NEO) approach captures these effects by treating selected nuclei quantum mechanically on the same footing as electrons. On classical computers, the resources required for an exact solution of NEO-based models grow exponentially with system size. By contrast, quantum computers offer a means of solving this problem with polynomial scaling. However, due to the limitations of current quantum devices, NEO simulations are confined to the smallest systems described by minimal basis sets, whereas realistic simulations beyond the Born-Oppenheimer approximation require more sophisticated basis sets. For this purpose, we herein extend a hardware-efficient ADAPT-VQE method to the NEO framework in the frozen natural orbital (FNO) basis. We demonstrate on H2 and D2 molecules that the NEO-FNO-ADAPT-VQE method reduces the CNOT count by several orders of magnitude relative to the NEO unitary coupled cluster method with singles and doubles while maintaining the desired accuracy. This extreme reduction in the CNOT gate count is sufficient to permit practical computations employing the NEO method─an important step toward accurate simulations involving nonclassical nuclei and non-Born-Oppenheimer effects on near-term quantum devices. We further show that the method can capture isotope effects, and we demonstrate that inclusion of correlation energy systematically improves the prediction of difference in the zero-point energy (ΔZPE) between isotopes.

Symmetry-Projected Nuclear-Electronic Hartree-Fock: Eliminating Rotational Energy Contamination

We present a symmetry projection technique for enforcing rotational and parity symmetries in nuclear-electronic Hartree-Fock wave functions, which treat electrons and nuclei on equal footing. The molecular Hamiltonian obeys rotational and parity inversion symmetries, which are, however, broken by expanding in Gaussian basis sets that are fixed in space. We generate a trial wave function with the correct symmetry properties by projecting the wave function onto representations of the three-dimensional rotation group, i.e., the special orthogonal group in three dimensions SO(3). As a consequence, the wave function becomes an eigenfunction of the angular momentum operator which (i) eliminates the contamination of the ground-state wave function by highly excited rotational states arising from the broken rotational symmetry and (ii) enables the targeting of specific rotational states of the molecule. We demonstrate the efficiency of the symmetry projection technique by calculating the energies of the low-lying rotational states of the H2 and H3+ molecules.

Second-Order Self-Consistent Field Algorithms: From Classical to Quantum Nuclei

This work presents a general framework for deriving exact and approximate Newton self-consistent field (SCF) orbital optimization algorithms by leveraging concepts borrowed from differential geometry. Within this framework, we extend the augmented Roothaan-Hall (ARH) algorithm to unrestricted electronic and nuclear-electronic calculations. We demonstrate that ARH yields an excellent compromise between stability and computational cost for SCF problems that are hard to converge with conventional first-order optimization strategies. In the electronic case, we show that ARH overcomes the slow convergence of orbitals in strongly correlated molecules with the example of several iron-sulfur clusters. For nuclear-electronic calculations, ARH significantly enhances the convergence already for small molecules, as demonstrated for a series of protonated water clusters.

Divide-and-Conquer Linear-Scaling Quantum Chemical Computations

Fragmentation and embedding schemes are of great importance when applying quantum-chemical calculations to more complex and attractive targets. The divide-and-conquer (DC)-based quantum-chemical model is a fragmentation scheme that can be connected to embedding schemes. This feature article explains several DC-based schemes developed by the authors over the last two decades, which was inspired by the pioneering study of DC self-consistent field (SCF) method by Yang and Lee (J. Chem. Phys. 1995, 103, 5674-5678). First, the theoretical aspects of the DC-based SCF, electron correlation, excited-state, and nuclear orbital methods are described, followed by the two-component relativistic theory, quantum-mechanical molecular dynamics simulation, and the introduction of three programs, including DC-based schemes. Illustrative applications confirmed the accuracy and feasibility of the DC-based schemes.