Molecular structure in non-Born–Oppenheimer quantum mechanics

Gaussians for Electronic and Rovibrational Quantum Dynamics

The assumptions underpinning the adiabatic Born-Oppenheimer (BO) approximation are broken for molecules interacting with attosecond laser pulses, which generate complicated coupled electronic-nuclear wave packets that generally will have components of electronic and dissociation continua as well as bound-state contributions. The conceptually most straightforward way to overcome this challenge is to treat the electronic and nuclear degrees of freedom on equal quantum-mechanical footing by not invoking the BO approximation at all. Explicitly correlated Gaussian (ECG) basis functions have proved successful for non-BO calculations of stationary molecular states and energies, reproducing rovibrational absorption spectra with very high accuracy. In this Article, we present a proof-of-principle study of the ability of fully flexible ECGs (FFECGs) to capture the intricate electronic and rovibrational dynamics generated by short, high-intensity laser pulses. By fitting linear combinations of FFECGs to accurate wave function histories obtained on a large real-space grid for a regularized 2D model of the hydrogen atom and for the 2D Morse potential, we demonstrate that FFECGs provide a very compact description of laser-driven electronic and rovibrational dynamics.

Quantum Definition of Molecular Structure

Molecular structure, a key concept of chemistry, has remained elusive from the perspective of all-particle quantum mechanics, despite many efforts. Viewing molecular structure as a manifestation of strong statistical correlation between nuclear positions, we propose a practical method based on Markov chain Monte Carlo sampling and unsupervised machine learning. Application to the D3+ molecule unambiguously shows that it possesses an equilateral triangular structure. These results provide a major step forward in our understanding of the molecular structure from fundamental quantum principles.

Representation and conservation of angular momentum in the Born-Oppenheimer theory of polyatomic molecules

This paper concerns the representation of angular momentum operators in the Born-Oppenheimer theory of polyatomic molecules and the various forms of the associated conservation laws. Topics addressed include the question of whether these conservation laws are exactly equivalent or only to some order of the Born-Oppenheimer parameter κ = (m/M)^1/4 and what the correlation is between angular momentum quantum numbers in the various representations. These questions are addressed in both problems involving a single potential energy surface and those with multiple, strongly coupled surfaces and in both the electrostatic model and those for which fine structure and electron spin are important. The analysis leads to an examination of the transformation laws under rotations of the electronic Hamiltonian; of the basis states, both adiabatic and diabatic, along with their phase conventions; of the potential energy matrix; and of the derivative couplings. These transformation laws are placed in the geometrical context of the structures in the nuclear configuration space that are induced by rotations, which include the rotational orbits or fibers, the surfaces upon which the orientation of the molecule changes but not its shape, and the section, an initial value surface that cuts transversally through the fibers. Finally, it is suggested that the usual Born-Oppenheimer approximation can be replaced by a dressing transformation, that is, a sequence of unitary transformations that block-diagonalize the Hamiltonian. When the dressing transformation is carried out, we find that the angular momentum operator does not change. This is a part of a system of exact equivalences among various representations of angular momentum operators in Born-Oppenheimer theory. Our analysis accommodates large-amplitude motions and is not dependent on small-amplitude expansions about an equilibrium position. Our analysis applies to noncollinear configurations of a polyatomic molecule; this covers all but a subset of measure zero (the collinear configurations) in the nuclear configuration space.

On the Probability Density of the Nuclei in a Vibrationally Excited Molecule

For localized and oriented vibrationally excited molecules, the qualitative features of the one-body probability density of the nuclei (one-nucleus density) are investigated. Like the familiar and widely used one-electron density that represents the probability of finding an electron at a given location in space, the one-nucleus density represents the probability of finding a nucleus at a given position in space independent of the location of the other nuclei and independent of their type. In contrast to the electrons, however, the nuclei are comparably localized. Due to this localization of the individual nuclei, the one-nucleus density provides a quantum-mechanical representation of the "chemical picture" of the molecule as an object that can largely be understood in a three-dimensional space, even though its full nuclear probability density is defined on the high-dimensional configuration space of all the nuclei. We study how the nodal structure of the wavefunctions of vibrationally excited states translates to the one-nucleus density. It is found that nodes do not necessarily lead to visible changes in the one-nucleus density: Already for relatively small molecules, only certain vibrational excitations change the one-nucleus density qualitatively compared to the ground state. It turns out that there are simple rules for predicting the shape of the one-nucleus density from the normal mode coordinates. A Python module for the computation of the one-nucleus density is provided at https://gitlab.com/axelschild/mQNMc.

Explicitly correlated Gaussian functions with shifted-center and projection techniques in pre-Born-Oppenheimer calculations

Numerical projection methods are elaborated for the calculation of eigenstates of the non-relativistic many-particle Coulomb Hamiltonian with selected rotational and parity quantum numbers employing shifted explicitly correlated Gaussian functions, which are, in general, not eigenfunctions of the total angular momentum and parity operators. The increased computational cost of numerically projecting the basis functions onto the irreducible representations of the three dimensional rotation-inversion group is the price to pay for the increased flexibility of the basis functions. This increased flexibility allowed us to achieve a substantial improvement for the variational upper bound to the Pauli-allowed ground-state energy of the H 3 + = { p + , p + , p + , e - , e - } molecular ion treated as an explicit five-particle system. We compare our pre-Born-Oppenheimer result obtained for this molecular ion with rotational-vibrational calculations carried out on a potential energy surface.

Hidden aspects of the Structural theory of chemistry: MC-QTAIM analysis reveals "alchemical" transformation from a triatomic to a diatomic structure

The Structural theory of chemistry introduces chemical/molecular structure as a combination of relative arrangement and bonding patterns of atoms in a molecule. Nowadays, the structure of atoms in molecules is derived from the topological analysis of the quantum theory of atoms in molecules (QTAIM). In this context, a molecular structure is varied by large geometrical variations and concomitant reorganization of electronic structure that usually take place in chemical reactions or under extreme hydrostatic pressure. In this report, a new mode of structural variation is introduced within the context of the newly proposed multi-component QTAIM (MC-QTAIM) that originates from the mass variation of nuclei. Accordingly, XCN and CNX series of species are introduced where X stands for a quantum particle with a unit of positive charge and a variable mass that is varied in discrete steps between the mass of a proton and a positron. Ab initio non-Born-Oppenheimer (non-BO) calculations are done on both series of species and the resulting non-BO wavefunctions are used for the MC-QTAIM analysis, revealing a triatomic structure for the proton mass and a diatomic structure for the positron mass. In both series of species, a critical mass between that of proton and positron mass is discovered where the transition from triatomic to diatomic structure takes place. This abrupt structural transformation has a topological nature resembling the usual phase transitions in thermodynamics. The discovered mass-induced structural transformation is a hidden aspect of the Structural theory which is revealed only beyond the BO paradigm, when nuclei are treated as quantum waves instead of clamped point charges.

Non-adiabatic dynamics close to conical intersections and the surface hopping perspective

Conical intersections play a major role in the current understanding of electronic de-excitation in polyatomic molecules, and thus in the description of photochemistry and photophysics of molecular systems. This article reviews aspects of the basic theory underlying the description of non-adiabatic transitions at conical intersections, with particular emphasis on the important case when the dynamics of the nuclei are treated classically. Within this classical nuclear motion framework, the main aspects of the surface hopping methodology in the conical intersection context are presented. The emerging picture from this treatment is that of electronic transitions around conical intersections dominated by the interplay of the nuclear velocity and the derivative non-adiabatic coupling vector field.

High angular momentum states of lithium atom, studied with symmetry-projected explicitly correlated Gaussian lobe functions

Method of construction of wave functions approximating eigenfunctions of the [Formula: see text] operator is proposed for high angular momentum states of few-electron atoms. Basis functions are explicitly correlated Gaussian lobes, projected onto irreducible representations of finite point groups. Variational calculations have been carried out for the lowest states of lithium atom, with quantum number L in the range from 1 to 8. Nonrelativistic energies accurate to several dozens of nanohartree have been obtained. For 2(2)P, 3(2)D, and 4(2)F states they agree well with the reference results. Transition frequencies have been computed and compared with available experimental data.

Non-Born-Oppenheimer electronic and nuclear densities for a Hooke-Calogero three-particle model: non-uniqueness of density-derived molecular structure

We consider the calculation of non-Born-Oppenheimer, nBO, one-particle densities for both electrons and nuclei. We show that the nBO one-particle densities evaluated in terms of translationally invariant coordinates are independent of the wavefunction describing the motion of center of mass of the whole system. We show that they depend, however, on an arbitrary reference point from which the positions of the vectors labeling the particles are determined. We examine the effect that this arbitrary choice has on the topology of the one-particle density by selecting the Hooke-Calogero model of a three-body system for which expressions for the one-particle densities can be readily obtained in analytic form. We extend this analysis to the one-particle densities obtained from full Coulomb interaction wavefunctions for three-body systems. We conclude, in view of the fact that there is a close link between the choice of the reference point and the topology of one-particle densities that the molecular structure inferred from the topology of these densities is not unique. We analyze the behavior of one-particle densities for the Hooke-Calogero Born-Oppenheimer, BO, wavefunction and show that topological transitions are also present in this case for a particular mass value of the light particles even though in the BO regime the nuclear masses are infinite. In this vein, we argue that the change in topology caused by variation of the mass ratio between light and heavy particles does not constitute a true indication in the nBO regime of the emergence of molecular structure.

Matrix elements of N-particle explicitly correlated Gaussian basis functions with complex exponential parameters

In this work we present analytical expressions for Hamiltonian matrix elements with spherically symmetric, explicitly correlated Gaussian basis functions with complex exponential parameters for an arbitrary number of particles. The expressions are derived using the formalism of matrix differential calculus. In addition, we present expressions for the energy gradient that includes derivatives of the Hamiltonian integrals with respect to the exponential parameters. The gradient is used in the variational optimization of the parameters. All the expressions are presented in the matrix form suitable for both numerical implementation and theoretical analysis. The energy and gradient formulas have been programmed and used to calculate ground and excited states of the He atom using an approach that does not involve the Born-Oppenheimer approximation.

On equilibrium structures of the water molecule

Equilibrium structures are fundamental entities in molecular sciences. They can be inferred from experimental data by complicated inverse procedures which often rely on several assumptions, including the Born-Oppenheimer approximation. Theory provides a direct route to equilibrium geometries. A recent high-quality ab initio semiglobal adiabatic potential-energy surface (PES) of the electronic ground state of water, reported by Polyansky et al. [ ibid. 299, 539 (2003)] and called CVRQD here, is analyzed in this respect. The equilibrium geometries resulting from this direct route are deemed to be of higher accuracy than those that can be determined by analyzing experimental data. Detailed investigation of the effect of the breakdown of the Born-Oppenheimer approximation suggests that the concept of an isotope-independent equilibrium structure holds to about 3 x 10(-5) A and 0.02 degrees for water. The mass-independent [Born-Oppenheimer (BO)] equilibrium bond length and bond angle on the ground electronic state PES of water is r(e) (BO)=0.957 82 A and theta e (BO)=104.48(5) degrees , respectively. The related mass-dependent (adiabatic) equilibrium bond length and bond angle of H2 (16)O is r(e) (ad)=0.957 85 A and theta e (ad)=104.50(0) degrees , respectively, while those of D2 (16)O are r(e) (ad)=0.957 83 A and theta e (ad)=104.49(0) degrees . Pure ab initio prediction of J=1 and 2 rotational levels on the vibrational ground state by the CVRQD PESs is accurate to better than 0.002 cm(-1) for all isotopologs of water considered. Elaborate adjustment of the CVRQD PESs to reproduce all observed rovibrational transitions to better than 0.05 cm(-1) (or the lower ones to better than 0.0035 cm(-1)) does not result in noticeable changes in the adiabatic equilibrium structure parameters. The expectation values of the ground vibrational state rotational constants of the water isotopologs, computed in the Eckart frame using the CVRQD PESs and atomic masses, deviate from the experimentally measured ones only marginally, especially for A0 and B0. The small residual deviations in the effective rotational constants are due to centrifugal distortion, electronic, and non-Born-Oppenheimer effects. The spectroscopic (nonadiabatic) equilibrium structural parameters of H2 16O, obtained from experimentally determined A'0 and B'0 rotational constants corrected empirically to obtain equilibrium rotational constants, are r(e) (sp)=0.957 77 A and theta e (sp)=104.48 degrees .

Non-Born–Oppenheimer molecular structure and one-particle densities for H2D+

We show that the nonadiabatic (non-Born-Oppenheimer) ground state of a three-nuclei system can be effectively calculated with the use of an explicitly correlated Gaussian basis set with floating centers. Sample calculations performed for the H2D+ system with various basis set sizes show good convergence with respect to both the total energy and the expectation values of the internuclear distances (molecular geometry), the distances between the nuclei and the electrons, and between the electrons. We also provide a derivation of the formulas for one-particle density calculations and some density plots showing the spatial distribution of the H2D+ nuclear and electronic densities.

Non-Born–Oppenheimer variational calculations of HT+ bound states with zero angular momentum

We report fully nonadiabatic calculations of all rotationless bound states of HT+ molecular ion (t+p+e-) carried out in the framework of the variational method. We show that, in all the states, except the two highest ones, the bond in the system can be described as covalent. In the highest two states the bond becomes essentially ionic and HT+ can be described as a T+H+ complex. The wave function of the system was expanded in terms of spherically symmetric, explicitly correlated Gaussian functions with preexponential multipliers consisting of powers of the internuclear distance. Apart from the total energies of the states, we have calculated the expectation values of the t-p, t-e, and p-e interparticle distances, their squares, and the nucleus-nucleus correlation functions.

Molecular structure calculations without clamping the nuclei

A number of recent papers have considered ways in which molecular structure may be calculated when both the electrons and the nuclei are treated from the outset as quantum particles. This is in contrast to the conventional approach in which the nuclei initially have their positions fixed and so merely provide a potential for electronic motion. The usual approach is generally assumed to be justified by the 1927 work of Born and Oppenheimer. In this paper we discuss what precisely might be anticipated in the way of molecular structure from a mathematical consideration of the spectral properties of the full Coulomb Hamiltonian, to what extent the very idea of molecular structure might be dependent upon treating the nuclei simply as providing a potential and the extent to which the work of Born and Oppenheimer can be used to support this position.