Comment on ‘Molecular structure in non-Born–Oppenheimer quantum mechanics’

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.

MC-QTAIM analysis reveals an exotic bond in coherently quantum superposed malonaldehyde

The proton between the two oxygen atoms of the malonaldehyde molecule experiences an effective double-well potential in which the proton's wavefunction is delocalized between the two wells. Herein we employ a state-of-the-art multi-component quantum theory of atoms in molecules partitioning scheme to obtain the molecular structure, i.e. atoms in molecules and bonding network, from the superposed ab initio wavefunctions of malonaldehyde. In contrast to the familiar clamped-proton portrayal of malonaldehyde, in which the proton forms a hydrogen basin, for the superposed states the hydrogen basin disappears and two novel hybrid oxygen-hydrogen basins appear instead, with an even distribution of the proton population between the two basins. The interaction between the hybrid basins is stabilizing thanks to an unprecedented mechanism. This involves the stabilizing classical Coulomb interaction of the one-proton density in one of the basins with one-electron density in the other basin. This stabilizing mechanism yields a bond foreign to the known bonding modes in chemistry.

Density Functional Computations of Vibrational Circular Dichroism Spectra beyond the Born-Oppenheimer Approximation

Transition-metal complexes provide rich features in vibrational circular dichroism (VCD) spectra, including significant intensity enhancements, and become thus useful in structural and functional studies of molecules. Quite often, however, the vibrational spectral bands are mixed with the electronic ones, and interpretation of such experiments is difficult. In the present study, we elaborate on the theory needed to calculate the VCD intensities beyond the Born-Oppenheimer (BO) approximation. Within a perturbation approach, the coupling between the electronic and vibrational states is estimated using the harmonic approximation and simplified wave functions obtainable from common density functional theory (DFT) computations. Explicit expressions, including Slater determinants and derivatives of molecular orbitals, are given. On a model diamine complex, the implementation is tested and factors affecting spectral intensities and frequencies are investigated. For two larger molecules, the results are in a qualitative agreement with previous experimental data. Typically, the electronic-vibrational interaction Hamiltonian coupling elements are rather small (∼0 to 10 cm^-1), which provides negligible contributions to vibrational frequencies and absorption intensities. However, significant changes in VCD spectra are induced due to the large transition magnetic dipole moment associated with the d-d metal transitions. The possibility to model the spectra beyond the BO limit opens the way to further applications of chiral spectroscopy and transition-metal complexes.

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.

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.

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.

Notes on quantitative structure-properties relationships (QSPR) part 2: the role of the number of atoms as a molecular descriptor

A previous analysis performed in our laboratory about the polynomial dependency of the atomic quantum self-similarity measures on the atomic number, together with recent publications on quantitative structure-properties relationships (QSPR), based on the number of molecular atoms, published by various authors, have driven us to show here that a simplified form of the fundamental quantum QSPR (QQSPR) equation, permits to theoretically demonstrate the important, but obvious, role of the number of atoms in a molecule, as a possible molecular descriptor. A discussion of the practical use of the number of atoms in QSPR is also given at the end, which also contains a discussion on the role of Ockham's razor in descriptor simplification choices.

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.