Inclusion of explicit electron-proton correlation in the nuclear-electronic orbital approach using Gaussian-type geminal functions

Nuclear-electronic all-particle density matrix renormalization group

We introduce the Nuclear-Electronic All-Particle Density Matrix Renormalization Group (NEAP-DMRG) method for solving the time-independent Schrödinger equation simultaneously for electrons and other quantum species. In contrast to the already existing multicomponent approaches, in this work, we construct from the outset a multi-reference trial wave function with stochastically optimized non-orthogonal Gaussian orbitals. By iterative refining of the Gaussians' positions and widths, we obtain a compact multi-reference expansion for the multicomponent wave function. We extend the DMRG algorithm to multicomponent wave functions to take into account inter- and intra-species correlation effects. The efficient parameterization of the total wave function as a matrix product state allows NEAP-DMRG to accurately approximate the full configuration interaction energies of molecular systems with more than three nuclei and 12 particles in total, which is currently a major challenge for other multicomponent approaches. We present the NEAP-DMRG results for two few-body systems, i.e., H₂ and H₃ ⁺, and one larger system, namely, BH₃.

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.

Developing effective electronic-only coupled-cluster and Møller-Plesset perturbation theories for the muonic molecules

Recently we have proposed an effective Hartree-Fock (EHF) theory for the electrons of the muonic molecules that is formally equivalent to the HF theory within the context of the nuclear-electronic orbital theory [Phys. Chem. Chem. Phys., 2018, 20, 4466]. In the present report we extend the muon-specific effective electronic structure theory beyond the EHF level by introducing the effective second order Møller-Plesset perturbation theory (EMP2) and the effective coupled-cluster theory at single and double excitation levels (ECCSD) as well as an improved version including perturbative triple excitations (ECCSD(T)). These theories incorporate electron-electron correlation into the effective paradigm and through their computational implementation, a diverse set of small muonic species is considered as a benchmark at these post-EHF levels. A comparative computational study on this set demonstrates that the muonic bond length is in general non-negligibly longer than corresponding hydrogenic analogs. Next, the developed post-EHF theories are applied for the muoniated N-heterocyclic carbene/silylene/germylene and the muoniated triazolium cation revealing the relative stability of the sticking sites of the muon in each species. The computational results, in line with previously reported experimental data demonstrate that the muon generally prefers to attach to the divalent atom with carbeneic nature. A detailed comparison of these muonic adducts with the corresponding hydrogenic adducts reveals subtle differences that have already been overlooked.

Toward a muon-specific electronic structure theory: effective electronic Hartree-Fock equations for muonic molecules

An effective set of Hartree-Fock (HF) equations are derived for electrons of muonic systems, i.e., molecules containing a positively charged muon, conceiving the muon as a quantum oscillator, which are completely equivalent to the usual two-component HF equations used to derive stationary states of the muonic molecules. In these effective equations, a non-Coulombic potential is added to the orthodox coulomb and exchange potential energy terms, which describes the interaction of the muon and the electrons effectively and is optimized during the self-consistent field cycles. While in the two-component HF equations a muon is treated as a quantum particle, in the effective HF equations it is absorbed into the effective potential and practically transformed into an effective potential field experienced by electrons. The explicit form of the effective potential depends on the nature of muon's vibrations and is derivable from the basis set used to expand the muonic spatial orbital. The resulting effective Hartree-Fock equations are implemented computationally and used successfully, as a proof of concept, in a series of muonic molecules containing all atoms from the second and third rows of the Periodic Table. To solve the algebraic version of the equations muon-specific Gaussian basis sets are designed for both muon and surrounding electrons and it is demonstrated that the optimized exponents are quite distinct from those derived for the hydrogen isotopes. The developed effective HF theory is quite general and in principle can be used for any muonic system while it is the starting point for a general effective electronic structure theory that incorporates various types of quantum correlations into the muonic systems beyond the HF equations.

Multicomponent Density Functional Theory: Impact of Nuclear Quantum Effects on Proton Affinities and Geometries

Nuclear quantum effects such as zero point energy play a critical role in computational chemistry and often are included as energetic corrections following geometry optimizations. The nuclear-electronic orbital (NEO) multicomponent density functional theory (DFT) method treats select nuclei, typically protons, quantum mechanically on the same level as the electrons. Electron-proton correlation is highly significant, and inadequate treatments lead to highly overlocalized nuclear densities. A recently developed electron-proton correlation functional, epc17, has been shown to provide accurate nuclear densities for molecular systems. Herein, the NEO-DFT/epc17 method is used to compute the proton affinities for a set of molecules and to examine the role of nuclear quantum effects on the equilibrium geometry of FHF^-. The agreement of the computed results with experimental and benchmark values demonstrates the promise of this approach for including nuclear quantum effects in calculations of proton affinities, pK_a's, optimized geometries, and reaction paths.

On the physical interpretation of the nuclear molecular orbital energy

Recently, several groups have extended and implemented molecular orbital (MO) schemes to simultaneously obtain wave functions for electrons and selected nuclei. Many of these schemes employ an extended Hartree-Fock approach as a first step to find approximate electron-nuclear wave functions and energies. Numerous studies conducted with these extended MO methodologies have explored various effects of quantum nuclei on physical and chemical properties. However, to the best of our knowledge no physical interpretation has been assigned to the nuclear molecular orbital energy (NMOE) resulting after solving extended Hartree-Fock equations. This study confirms that the NMOE is directly related to the molecular electrostatic potential at the position of the nucleus.

Incorporating nuclear vibrational energies into the "atom in molecules" analysis: An analytical study

The quantum theory of atoms in molecules (QTAIM) is based on the clamped nucleus paradigm and solely working with the electronic wavefunctions, so does not include nuclear vibrations in the AIM analysis. On the other hand, the recently extended version of the QTAIM, called the multi-component QTAIM (MC-QTAIM), incorporates both electrons and quantum nuclei, i.e., those nuclei treated as quantum waves instead of clamped point charges, into the AIM analysis using non-adiabatic wavefunctions. Thus, the MC-QTAIM is the natural framework to incorporate the role of nuclear vibrations into the AIM analysis. In this study, within the context of the MC-QTAIM, the formalism of including nuclear vibrational energy in the atomic basin energy is developed in detail and its contribution is derived analytically using the recently proposed non-adiabatic Hartree product nuclear wavefunction. It is demonstrated that within the context of this wavefunction, the quantum nuclei may be conceived pseudo-adiabatically as quantum oscillators and both isotropic harmonic and anisotropic anharmonic oscillator models are used to compute the zero-point nuclear vibrational energy contribution to the basin energies explicitly. Inspired by the results gained within the context of the MC-QTAIM analysis, a heuristic approach is proposed within the context of the QTAIM to include nuclear vibrational energy in the basin energy from the vibrational wavefunction derived adiabatically. The explicit calculation of the basin contribution of the zero-point vibrational energy using the uncoupled harmonic oscillator model leads to results consistent with those derived from the MC-QTAIM.

Extending the Domain-Averaged Exchange-Correlation Energies Within the Context of the MC-QTAIM: Tracing Subtle Variations Induced by Isotope Substitution

Recently, it has been demonstrated that the domain-averaged exchange-correlation energies, V^xc , are capable of tracing the covalent character of atom-atom interactions unequivocally and thus pave the way for detailed bonding analysis within the context of the quantum theory of atoms in molecules (QTAIM) [M. García-Revilla, E. Francisco, P. L. Popelier, A. Martín Pendás, ChemPhysChem 2013, 14, 1211-1218]. Herein, the concept of V^xc is extended within the context of the newly developed multicomponent QTAIM (MC-QTAIM). The extended version, Veexc , is capable of analyzing nonadiabatic wavefunctions and thus is sensitive to the mass of nuclei and can trace "locally" the subtle electronic variations induced by isotope substitution. To demonstrate this capability in practice, ab initio nonadiabatic wavefunctions for three isotopically substituted hydrogen cyanide molecules, in which the hydrogen nucleus was assumed to be a proton, deuterium, or tritium, were derived. The resulting wavefunctions were then used to compute Veexc and it emerged that for the hydrogen-carbon bond, the Veexc was distinct for each isotopic composition and varied in line with chemical expectations. Indeed, the introduction of Veexc paves the way for the investigation of vast numbers of structural and kinetic isotope effects within the context of the MC-QTAIM.

The divide-and-conquer second-order proton propagator method based on nuclear orbital plus molecular orbital theory for the efficient computation of proton binding energies

An efficient computational method to evaluate the binding energies of many protons in large systems was developed.

Electric transition dipole moment in pre-Born-Oppenheimer molecular structure theory

This paper presents the calculation of the electric transition dipole moment in a pre-Born-Oppenheimer framework. Electrons and nuclei are treated equally in terms of the parametrization of the non-relativistic total wave function, which is written as a linear combination of basis functions constructed from explicitly correlated Gaussian functions and the global vector representation. The integrals of the electric transition dipole moment are derived corresponding to these basis functions in both the length and the velocity representation. The calculations are performed in laboratory-fixed Cartesian coordinates without relying on coordinates which separate the center of mass from the translationally invariant degrees of freedom. The effect of the overall motion is eliminated through translationally invariant integral expressions. The electric transition dipole moment is calculated between two rovibronic levels of the H2 molecule assignable to the lowest rovibrational states of the X (1)Σ(g)(+) and B (1)Σ(u)(+) electronic states in the clamped-nuclei framework. This is the first evaluation of this quantity in a full quantum mechanical treatment without relying on the Born-Oppenheimer approximation.

Calculation of electron-hole recombination probability using explicitly correlated Hartree-Fock method

The electron-hole explicitly correlated Hartree-Fock method (eh-XCHF) is presented as a general strategy for investigation of electron-hole correlation and computation of electron-hole recombination probability. The eh-XCHF method is a variational method which uses explicitly correlated wavefunction that depends on the electron-hole inter-particle distances. It is shown that the explicitly correlated ansatz provides a systematic route to variationally minimize the total energy. The parabolic quantum dot is used as the benchmark system and the eh-XCHF method is used for computation of the ground state energy and electron-hole recombination probability. The results are compared to Hartree-Fock and explicitly correlated full configuration interaction (R12-FCI) calculations. The results indicate that an accurate description of the electron-hole wavefunction at short electron-hole inter-particle distances is crucial for qualitative description of the electron-hole recombination probability. The eh-XCHF method successfully addresses this issue and comparison of eh-XCHF calculations with R12-FCI shows good agreement. The quality of the mean field approximation for electron-hole system is also investigated by comparing HF and R12-FCI energies for electron-electron and electron-hole systems. It was found that performance of the mean field approximation is worse for the electron-hole system as compared to the corresponding electron-electron system.