Multicomponent equation-of-motion coupled cluster singles and doubles: Theory and calculation of excitation energies for positronium hydride

Production of positronium chloride: A study of the charge exchange reaction between Ps and Cl

We present cross sections for the formation of positronium chloride (PsCl) in its ground state from the charge exchange between positronium (Ps) and chloride (Cl-) in the range of 10 meV-100 eV Ps energy. We have used theoretical models based on the first Born approximation in its three-body formulation. We simulated the collisions between Ps and Cl- using ab initio binding energies and positronic wave functions at both the mean-field and correlated levels extrapolated to the complete basis set limit. The accuracy of these ab initio data was benchmarked on the PsF system with the existing highly accurate results, including the very recent quantum Monte Carlo results. We have investigated Ps excited states up to n = 4. The results suggest that the channel Ps(n = 2) is of particular interest for the production of PsCl in the ground state and shows that an accurate treatment of correlation effects (i.e., electron-electron and electron-positron correlations) leads to a significant change in the magnitude of the PsCl production cross section with respect to the mean-field level.

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.

Nuclear-Electronic Orbital Quantum Mechanical/Molecular Mechanical Real-Time Dynamics

Simulating the nuclear-electronic quantum dynamics of large-scale molecular systems in the condensed phase is key for studying biologically and chemically important processes such as proton transfer and proton-coupled electron transfer reactions. Herein, the real-time nuclear-electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) approach is combined with a hybrid quantum mechanical/molecular mechanical (QM/MM) strategy to enable the accurate description of coupled nuclear-electronic quantum dynamics in the presence of heterogeneous environments such as solvent or proteins. The densities of the electrons and quantum protons are propagated in real time, while the other nuclei are propagated classically on the instantaneous electron-proton vibronic surface. This approach is applied to phenol bound to lysozyme, intramolecular proton transfer in malonaldehyde, and nonequilibrium excited-state intramolecular proton transfer in o-hydroxybenzaldehyde. These examples illustrate that the RT-NEO-TDDFT framework, coupled with an atomistic representation of the environment, allows the simulation of condensed-phase systems that exhibit significant nuclear quantum effects.

Nuclear-Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear-electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. These developments and initial examples provide the foundation for future studies of nuclear-electronic quantum dynamics in complex chemical and biological environments.

Time-dependent nuclear-electronic orbital Hartree-Fock theory in a strong uniform magnetic field

In an ultrastrong magnetic field, with field strength B ≈ B₀ = 2.35 × 10⁵ T, molecular structure and dynamics differ strongly from that observed on the Earth. Within the Born-Oppenheimer (BO) approximation, for example, frequent (near) crossings of electronic energy surfaces are induced by the field, suggesting that nonadiabatic phenomena and processes may play a more important role in this mixed-field regime than in the weak-field regime on Earth. To understand the chemistry in the mixed regime, it therefore becomes important to explore non-BO methods. In this work, the nuclear-electronic orbital (NEO) method is employed to study protonic vibrational excitation energies in the presence of a strong magnetic field. The NEO generalized Hartree-Fock theory and time-dependent Hartree-Fock (TDHF) theory are derived and implemented, accounting for all terms that result as a consequence of the nonperturbative treatment of molecular systems in a magnetic field. The NEO results for HCN and FHF^- with clamped heavy nuclei are compared against the quadratic eigenvalue problem. Each molecule has three semi-classical modes owing to the hydrogen-two precession modes that are degenerate in the absence of a field and one stretching mode. The NEO-TDHF model is found to perform well; in particular, it automatically captures the screening effects of the electrons on the nuclei, which are quantified through the difference in energy of the precession modes.

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.

Hydrogen-Atom Electronic Basis Sets for Multicomponent Quantum Chemistry

Multicomponent methods are a conceptually simple way to include nuclear quantum effects into quantum chemistry calculations. In multicomponent methods, the electronic molecular orbitals are described using the linear combination of atomic orbitals approximation. This requires the selection of a one-particle electronic basis set which, in practice, is commonly a correlation-consistent basis set. In multicomponent method studies, it has been demonstrated that large electronic basis sets are required for quantum hydrogen nuclei to accurately describe electron-nuclear correlation. However, as we show in this study, much of the need for large electronic basis sets is due to the correlation-consistent electronic basis sets not being optimized to describe nuclear properties and electron-nuclear correlation. Herein, we introduce a series of correlation-consistent electronic basis sets for hydrogen atoms called cc-pVnZ-mc with additional basis functions optimized to reproduce multicomponent density functional theory protonic densities. These new electronic basis sets are shown to yield better protonic densities with fewer electronic basis functions than the standard correlation-consistent basis sets and reproduce other protonic properties such as proton affinities and protonic excitation energies, even though they were not optimized for these purposes. The cc-pVnZ-mc basis sets should enable multicomponent many-body calculations on larger systems due to the improved computational efficiency they provide for a given level of accuracy.

(T) Correction for Multicomponent Coupled-Cluster Theory for a Single Quantum Proton

(T) and [T] perturbative corrections are derived for multicomponent coupled-cluster theory with single and double excitations (CCSD). Benchmarking for systems with a single quantum proton shows that multicomponent CCSD methods that include perturbative corrections are more accurate than multicomponent CCSD for the calculation of proton affinities and absolute energies. An approximation is introduced that includes only (T) or [T] contributions from mixed electron-nuclear excitations.

Triple electron-electron-proton excitations and second-order approximations in nuclear-electronic orbital coupled cluster methods

The accurate description of nuclear quantum effects, such as zero-point energy, is important for modeling a wide range of chemical and biological processes. Within the nuclear-electronic orbital (NEO) approach, such effects are incorporated in a computationally efficient way by treating electrons and select nuclei, typically protons, quantum mechanically with molecular orbital techniques. Herein, we implement and test a NEO coupled cluster method that explicitly includes the triple electron-proton excitations, where two electrons and one proton are excited simultaneously. Our calculations show that this NEO-CCSD(eep) method provides highly accurate proton densities and proton affinities, outperforming any previously studied NEO method. These examples highlight the importance of the triple electron-electron-proton excitations for an accurate description of nuclear quantum effects. Additionally, we also implement and test the second-order approximate coupled cluster with singles and doubles (NEO-CC2) method, as well as its scaled-opposite-spin (SOS) versions. The NEO-SOS$'$-CC2 method, which scales the electron-proton correlation energy as well as the opposite-spin and same-spin components of the electron-electron correlation energy, achieves nearly the same accuracy as the NEO-CCSD(eep) method for the properties studied. Because of its low computational cost, this method will enable a wide range of chemical and photochemical applications for large molecular systems. This work sets the stage for a wide range of developments and applications within the NEO framework.

Theoretical perspectives on non-Born-Oppenheimer effects in chemistry

The Born-Oppenheimer approximation, which assumes that the electrons respond instantaneously to the motion of the nuclei, breaks down for a wide range of chemical and biological processes. The rate constants of such nonadiabatic processes can be calculated using analytical theories, and the real-time nonequilibrium dynamics can be described using numerical atomistic simulations. The selection of an approach depends on the desired balance between accuracy and efficiency. The computational expense of generating potential energy surfaces on-the-fly often favours the use of approximate, robust and efficient methods such as trajectory surface hopping for large, complex systems. The development of formally exact non-Born-Oppenheimer methods and the exploration of well-defined approximations to such methods are critical for providing benchmarks and preparing for the next generation of faster computers. Thus, the parallel development of rigorous but computationally expensive methods and more approximate but computationally efficient methods is optimal. This Perspective briefly summarizes the available theoretical and computational non-Born-Oppenheimer methods and presents examples illustrating how analytical theories and nonadiabatic dynamics simulations can elucidate the fundamental principles of chemical and biological processes. These examples also highlight how theoretical calculations are able to guide the interpretation of experimental data and provide experimentally testable predictions for nonadiabatic processes. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.

Multicomponent heat-bath configuration interaction with the perturbative correction for the calculation of protonic excited states

In this study, we extend the multicomponent heat-bath configuration interaction (HCI) method to excited states. Previous multicomponent HCI studies have been performed using only the variational stage of the HCI algorithm as they have largely focused on the calculation of protonic densities. Because this study focuses on energetic quantities, a second-order perturbative correction after the variational stage is essential. Therefore, this study implements the second-order Epstein-Nesbet correction to the variational stage of multicomponent HCI for the first time. Additionally, this study introduces a new procedure for calculating reference excitation energies for multicomponent methods using the Fourier-grid Hamiltonian (FGH) method, which should allow the one-particle electronic basis set errors to be better isolated from errors arising from an incomplete description of electron-proton correlation. The excited-state multicomponent HCI method is benchmarked by computing protonic excitations of the HCN and FHF^- molecules and is shown to be of similar accuracy to previous excited-state multicomponent methods such as the multicomponent time-dependent density-functional theory and equation-of-motion coupled-cluster theory relative to the new FGH reference values.

Direct Dynamics with Nuclear-Electronic Orbital Density Functional Theory

Direct dynamics simulations of chemical reactions typically require the selection of a method for generating the potential energy surfaces and a method for the dynamical propagation of the nuclei on these surfaces. The nuclear-electronic orbital (NEO) framework avoids this Born-Oppenheimer separation by treating specified nuclei on the same level as the electrons with wave function methods or density functional theory (DFT). The NEO approach is particularly applicable to proton, hydride, and proton-coupled electron transfer reactions, where the transferring proton(s) and all electrons are treated quantum mechanically. In this manner, the zero-point energy, density delocalization, and anharmonicity of the transferring protons are inherently and efficiently included in the energies, optimized geometries, and dynamics.This Account describes how various NEO methods can be used for direct dynamics simulations on electron-proton vibronic surfaces. The strengths and limitations of these approaches are discussed, and illustrative examples are presented. The NEO-DFT method can be used to simulate chemical reactions on the ground state vibronic surface, as illustrated by the application to hydride transfer in C₄H₉⁺. The NEO multistate DFT (NEO-MSDFT) method is useful for simulating ground state reactions in which the proton density becomes bilobal during the dynamics, a characteristic of hydrogen tunneling, as illustrated by proton transfer in malonaldehyde. The NEO time-dependent DFT (NEO-TDDFT) method produces excited electronic, vibrational, and vibronic surfaces. The application of linear-response NEO-TDDFT to H₂ and H₃⁺, as well as the partially and fully deuterated counterparts, shows that this approach produces accurate fundamental vibrational excitation energies when all nuclei and all electrons are treated quantum mechanically. Moreover, when only specified nuclei are treated quantum mechanically, this approach can be used to optimize geometries on excited state vibronic surfaces, as illustrated by photoinduced single and double proton transfer systems, and to conduct adiabatic dynamics on these surfaces. The real-time NEO-TDDFT method provides an alternative approach for simulating nonequilibrium nuclear-electronic dynamics of such systems. These various NEO methods can be combined with nonadiabatic dynamics methods such as Ehrenfest and surface hopping dynamics to include the nonadiabatic effects between the quantum and classical subsystems. The real-time NEO-TDDFT Ehrenfest dynamics simulation of excited state intramolecular proton transfer in o-hydroxybenzaldehyde illustrates the power of this type of combined approach. The field of multicomponent quantum chemistry is in the early stages, and the methods discussed herein provide the foundation for a wide range of promising future directions to be explored. An appealing future direction is the expansion of the real-time NEO-TDDFT method to describe the dynamics of all nuclei and electrons on the same level. Direct dynamics simulations using NEO wave function methods such as equation-of-motion coupled cluster or multiconfigurational approaches are also attractive but computationally expensive options. The further development of NEO direct dynamics methods will enable the simulation of the nuclear-electronic dynamics for a vast array of chemical and biological processes that extend beyond the Born-Oppenheimer approximation.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Nuclear-electronic orbital methods: Foundations and prospects

The incorporation of nuclear quantum effects and non-Born-Oppenheimer behavior into quantum chemistry calculations and molecular dynamics simulations is a longstanding challenge. The nuclear-electronic orbital (NEO) approach treats specified nuclei, typically protons, quantum mechanically on the same level as the electrons with wave function and density functional theory methods. This approach inherently includes nuclear delocalization and zero-point energy in molecular energy calculations, geometry optimizations, reaction paths, and dynamics. It can also provide accurate descriptions of excited electronic, vibrational, and vibronic states as well as nuclear tunneling and nonadiabatic dynamics. Nonequilibrium nuclear-electronic dynamics simulations beyond the Born-Oppenheimer approximation can be used to investigate a wide range of excited state processes. This Perspective provides an overview of the foundational NEO methods and enumerates the prospects for using these methods as building blocks for future developments. The conceptual simplicity and computational efficiency of the NEO approach will enhance its accessibility and applicability to diverse chemical and biological systems.

Analytical Gradients for Nuclear-Electronic Orbital Time-Dependent Density Functional Theory: Excited-State Geometry Optimizations and Adiabatic Excitation Energies

The computational investigation of photochemical processes often entails the calculation of excited-state geometries, energies, and energy gradients. The nuclear-electronic orbital (NEO) approach treats specified nuclei, typically protons, quantum mechanically on the same level as the electrons, thereby including the associated nuclear quantum effects and non-Born-Oppenheimer behavior into quantum chemistry calculations. The multicomponent density functional theory (NEO-DFT) and time-dependent DFT (NEO-TDDFT) methods allow efficient calculations of ground and excited states, respectively. Herein, the analytical gradients are derived and implemented for the NEO-TDDFT method and the associated Tamm-Dancoff approximation (NEO-TDA). The programmable equations for these analytical gradients as well as the NEO-DFT analytical Hessian are provided. The NEO approach includes the anharmonic zero-point energy (ZPE) and density delocalization associated with the quantum protons as well as vibronic mixing in geometry optimizations and energy calculations of ground and excited states. The harmonic ZPE associated with the other nuclei can be computed via the NEO Hessian. This approach is used to compute the 0-0 adiabatic excitation energies for a set of nine small molecules with all protons quantized, exhibiting slight improvement over the conventional electronic approach. Geometry optimizations of two excited-state intramolecular proton-transfer systems, [2,2'-bipyridyl]-3-ol and [2,2'-bipyridyl]-3,3'-diol, are performed with one and two quantized protons, respectively. The NEO calculations for these systems produce electronically excited-state geometries with stronger intramolecular hydrogen bonds and similar relative stabilities compared to conventional electronic methods. This work provides the foundation for nonadiabatic dynamics simulations of fundamental processes such as photoinduced proton transfer and proton-coupled electron transfer.

Multicomponent Unitary Coupled Cluster and Equation-of-Motion for Quantum Computation

The variational quantum eigensolver (VQE) algorithm combined with the unitary coupled cluster (UCC) ansatz has been developed for the quantum computation of molecular energies and wave functions within the Born-Oppenheimer approximation. Herein, this approach is extended to multicomponent systems to enable the quantum mechanical treatment of more than one type of particle, such as electrons and positrons or electrons and nuclei, without invoking the Born-Oppenheimer approximation. Specifically, we introduce the multicomponent unitary coupled cluster (mcUCC) method combined with the VQE algorithm for the calculation of ground-state energies and wave functions as well as the multicomponent equation-of-motion (mcEOM) method for the calculation of excitation energies. These methods are developed within the nuclear-electronic orbital (NEO) framework and are formulated in the qubit basis to enable implementations on quantum computers. Moreover, these methods are used to calculate the ground-state energy and excitation energies of positronium hydride, where both electrons and the positron are treated quantum mechanically, as well as the H₂ molecule, where both electrons and one proton are treated quantum mechanically. These applications validate the implementation and provide benchmark data for future calculations. The errors due to Trotterization of the mcUCC ansatz are also analyzed. This formalism, as well as the accompanying computer code, will serve as the basis for applications to more complex multicomponent systems, such as simulations of photoinduced nonadiabatic molecular processes, on both classical and quantum computers.

Multicomponent CASSCF Revisited: Large Active Spaces Are Needed for Qualitatively Accurate Protonic Densities

Multicomponent methods seek to treat select nuclei, typically protons, fully quantum mechanically and equivalent to the electrons of a chemical system. In such methods, it is well-known that due to the neglect of electron-proton correlation, a Hartree-Fock (HF) description of the electron-proton interaction catastrophically fails leading to qualitatively incorrect protonic properties. In single-component quantum chemistry, the qualitative failure of HF is normally indicative of the need for multireference methods such as complete active space self-consistent field (CASSCF). While a multicomponent CASSCF method was implemented nearly 20 years ago, it is only able to perform calculations with very small active spaces (∼10⁵ multicomponent configurations). Therefore, in order to extend the realm of applicability of the multicomponent CASSCF method, this study derives and implements a new two-step multicomponent CASSCF method that uses multicomponent heat-bath configuration interaction for the configuration interaction step, enabling calculations with very large active spaces (up to 16 electrons in 48 orbitals). We find that large electronic active spaces are needed to obtain qualitatively accurate protonic densities for the HCN and FHF^- molecules. Additionally, the multicomponent CASSCF method implemented here should have further applications for double-well protonic potentials and systems that are inherently electronically multireference.

Frequency and Time Domain Nuclear-Electronic Orbital Equation-of-Motion Coupled Cluster Methods: Combination Bands and Electronic-Protonic Double Excitations

The accurate description of excited vibronic states is important for modeling a wide range of photoinduced processes. The nuclear-electronic orbital (NEO) approach, which treats specified protons on the same level as the electrons, can describe excited electronic-protonic states. Herein the multicomponent equation-of-motion coupled cluster with singles and doubles (NEO-EOM-CCSD) method and its time-domain counterpart, TD-NEO-EOM-CCSD, are developed and implemented. The application of these methods to the HCN molecule highlights their capabilities. These methods predict qualitatively reasonable energies and intensities for a combination band corresponding to simultaneous excitation of two vibrational modes, as well as an overtone. These methods also describe states with double excitation character, such as excited electronic-protonic states corresponding to the simultaneous excitation of an electron and a proton. The ability of the NEO-EOM-CCSD method and its time-dependent counterpart to describe combination bands, overtones, and double excitations will enable a wide range of photochemical applications.

Development of nuclear basis sets for multicomponent quantum chemistry methods

The nuclear-electronic orbital (NEO) framework provides a practical approach for directly incorporating nuclear quantum effects and non-Born-Oppenheimer effects of specified nuclei, typically protons, into quantum chemistry calculations. Multicomponent wave function based methods, such as NEO coupled cluster singles and doubles, and multicomponent density functional theory (DFT), such as NEO-DFT, require the appropriate selection of electronic and nuclear basis sets. Although a wide array of electronic basis sets are available, systematically developed nuclear basis sets that balance accuracy and efficiency have been lacking. Herein, a series of nuclear basis sets are developed and shown to be accurate and efficient for describing both ground and excited state properties of multicomponent systems in which electrons and specified protons are treated quantum mechanically. Three series of Gaussian-type nuclear basis sets, denoted PB4, PB5, and PB6, are developed with varying levels of angular momentum. A machine-learning optimization procedure relying on the Gaussian process regression method is utilized to accelerate the optimization process. The basis sets are validated in terms of predictions of ground state energies, proton densities, proton affinities, and proton vibrational excitation energies, allowing the user to select the desired balance between accuracy and efficiency for the properties of interest. These nuclear basis sets will enhance the tractability of NEO methods for applications to a wide range of chemical systems.

Real-Time Time-Dependent Nuclear-Electronic Orbital Approach: Dynamics beyond the Born-Oppenheimer Approximation

The quantum mechanical treatment of both electrons and nuclei is crucial in nonadiabatic dynamical processes such as proton-coupled electron transfer. The nuclear-electronic orbital (NEO) method provides an elegant framework for including nuclear quantum effects beyond the Born-Oppenheimer approximation. To enable the study of nonequilibrium properties, we derive and implement a real-time NEO (RT-NEO) approach based on time-dependent Hatree-Fock or density functional theory, in which the electronic and nuclear degrees of freedom are propagated in a time-dependent variational framework. Nuclear and electronic spectral features can be resolved from the time-dependent dipole moment computed using the RT-NEO method. The test cases show the dynamical interplay between the quantum nuclei and the electrons through vibronic coupling. Moreover, vibrational excitation in the RT-NEO approach is demonstrated by applying a resonant driving field, and electronic excitation is demonstrated by simulating excited state intramolecular proton transfer. This work shows that the RT-NEO approach is a promising tool to study nonadiabatic quantum dynamical processes within a time-dependent variational description for the coupled electronic and nuclear degrees of freedom.

Molecular Vibrational Frequencies with Multiple Quantum Protons within the Nuclear-Electronic Orbital Framework

The nuclear-electronic orbital (NEO) approach treats all electrons and specified nuclei, typically protons, on the same quantum mechanical level. Proton vibrational excitations can be calculated using multicomponent time-dependent density functional theory (NEO-TDDFT) for fixed classical nuclei. Recently the NEO-DFT(V) approach was developed to enable the calculation of molecular vibrational frequencies for modes composed of both classical and quantum nuclei. This approach uses input from NEO-TDDFT to construct an extended NEO Hessian that depends on the expectation values of the quantum protons as well as the classical nuclear coordinates. Herein strategies are devised for extending these approaches to molecules with multiple quantum protons in a self-contained, effective, and computationally practical manner. The NEO-TDDFT method is shown to describe vibrational excitations corresponding to collective nuclear motions, such as linear combinations of proton vibrational excitations. The NEO-DFT(V) approach is shown to incorporate the most significant anharmonic effects in the molecular vibrations, particularly for the hydrogen stretching modes. These theoretical strategies pave the way for a wide range of multicomponent quantum chemistry applications.

Covalent bonds in positron dihalides

We report a computational study on homo- and heteronuclear e⁺[X^–Y^–] compounds formed by two halide anions (X^–, Y^– = F^–, Cl^–, Br^–) and one positron.

We report a computational study on homo- and heteronuclear e⁺[X^–Y^–] compounds formed by two halide anions (X^–, Y^– = F^–, Cl^–, Br^–) and one positron. Our results indicate the formation of energetically stable positronic molecules in all cases. Analysis of the electron and positron densities points out that the formation of positron covalent bonds underlies the stabilization of the otherwise repelling dihalides, revealing that positronic bonding can reach far beyond the previously addressed e⁺[H^–H^–] molecule [J. Charry, M. T. do N. Varella and A. Reyes, Angew. Chem. Int. Ed., 2018, 57, 8859–8864.]. To a significant extent, the properties of the positron dihalides are similar to those of the purely electronic analogs, e^–[A⁺B⁺], molecular cations with isoelectronic atomic cores (A⁺, B⁺ = Na⁺, K⁺, Rb⁺) bound by one electron. The positron bonds in the e⁺[X^–Y^–] complexes are however stronger than those in the isoelectronic e^–[A⁺B⁺] counterparts, as the former have shorter bond lengths and higher bond energies. While an energy decomposition analysis points out that both electronic and positronic bonds essentially arise from electrostatic interactions, the more stable positron bonds are partly due to the higher polarizabilities of the dihalide anions, and partly to more significant contributions from correlation and relaxation effects.

Multicomponent coupled cluster singles and doubles and Brueckner doubles methods: Proton densities and energies

The nuclear-electronic orbital (NEO) framework enables computationally practical coupled cluster calculations of multicomponent molecular systems, in which all electrons and specified nuclei, typically protons, are treated quantum mechanically. In addition to energies, computing accurate proton densities is essential for the calculation of reliable molecular properties, including vibrationally averaged geometries and vibrational frequencies. Herein, the Lagrangian formalism for the multicomponent coupled cluster with single and double excitations (NEO-CCSD) method is derived and implemented. The multicomponent coupled cluster with double excitations method using optimized Brueckner orbitals, denoted as NEO-BCCD, is also developed. Both of these methods are used to compute the proton densities for two molecular systems. The results illustrate that orbital relaxation effects, which can be included either indirectly with the NEO-CCSD method or directly with the NEO-BCCD method, are critical for computing even qualitatively accurate proton densities. Both methods are also able to provide accurate proton affinities and vibrationally averaged optimized geometries. This Lagrangian formalism will enable the calculation of other properties such as analytical nuclear gradients and Hessians with NEO coupled cluster methods. Moreover, the accuracy of these methods may be improved systematically by the inclusion of higher-order excitations. Thus, this work provides the foundation for a wide range of future methodological developments and applications within the NEO framework.