Non-Born–Oppenheimer nuclear and electronic densities for a three-particle Hooke–Coulomb model
Graphical abstract
Introduction
The Born–Oppenheimer, BO, approximation has been, without any doubt, a fundamental concept underlying the development of quantum-mechanical approaches to the electronic structure of molecular and crystalline systems. However, the BO approximation breaks down for many systems and processes of chemical and physical interest (see, for example, Refs. [1], [2], [3], [4], [5], [6], [7], [8]).
Thus, in sight of this breakdown of the BO approximation and bearing in mind that some important chemical and physical phenomena occurring at the atomic level strongly depend on nuclear quantum effects, recent efforts in the field of quantum theory have been directed toward the development of methods based on a non-Born–Oppenheimer, nBO, regime [2], [9], [10], [11], [12], [13], [14], [15], [16], [17], [18], [19], [20], [21], [22], [23], [24]. In the particular case of chemistry, the attention placed on nBO treatments has reopened the debate of whether molecular structure, a basic concept in chemistry [25] can be derived from wave functions that include the behavior of non-clamped nuclei [22], [26], [27], [28], [29], [30], [31].
The problem is that the average distances and angles obtained from nBO wave functions generally are not sufficient to extract from them molecular structure [32], [33]. However, let us mention that as shown in recent work of Mátyus et al. [31] the triangular structure of the ion can be retrieved from the angular density corresponding to a highly accurate nBO variational wave function.
In the context of Bader’s Quantum Theory of Atoms in Molecules, QTAIM [34] a very basic relation has been shown to exist between the concept of molecular structure and the topology of the BO one-electron density, obtained either through theoretical or experimental means (in the latter case for solids) [34], [35], [36], [37]. In QTAIM, the one-electron density comes from a clamped-nuclei wave function. Recently, however, this restriction has been removed in the TC-QTAIM where topological atoms are constructed from a joint density made up from electronic and nuclear ones calculated in a nBO context [38], [39], [40], [41], [42].
Since nBO treatments are becoming available, the open question is whether from the topological features of nBO one-particle densities (i.e., electronic and nuclear) we may obtain molecular structure. This question, however, does not have an easy answer since the inclusion of the nuclear motion in the Schrödinger equation makes its solution much more difficult to obtain. As a result, analytical and variational wave functions, and associated one-particle densities, have been computed only for a small group of simple systems consisting of few particles: [22], [30], [43], [44] [45] [46] [22] [47] [48] [49] and [50]. Moreover, in these calculations there has been an implicit use of a given reference point from which the density is defined. As different reference points have been employed in the above calculations, it is difficult to compare the topological features of the ensuing densities [51], [52].
An alternative way for examining this question is provided by the use of exactly solvable models for few-particle systems as they yield analytic expressions for the one-particle densities. In particular, the Hooke–Calogero three-particle model has been previously employed in order to examine whether molecular structure arises in a nBO regime [53], [51]. However, as pointed out by Ludeña et al. [51] the topology of the nBO one-particle density depends upon the selection of the reference point. This non-uniqueness of the one-particle densities, observed for the Hooke–Calogero model showed that there does not exist a one-to-one relationship between density and molecular structure.
In the present work, we go beyond the Hooke–Calogero model [51] and introduce a more realistic Hooke–Coulomb model for a three-particle system where the harmonic interaction between identical particles are replaced by the true Coulomb potential. As in the present case, we can also obtain analytic solutions for the nBO Schrödinger equation, we calculate and analyze close-form expressions for electronic and nuclear one-particle densities for this three-particle system. However, since the analytic solutions are restricted to selected values of the coupling constant, we also obtain general variational solutions to this problem. We examine for both types of wave functions the effect that the selection of the reference points has on the topologies of the nuclear and electronic one-particle densities. In Section 2, we describe the Hooke–Coulomb model employed in the present work, In Section 3, we discuss the calculation of one-particle densities for different reference points, and consider the specific cases of the , , and systems. Finally in Section 4, we present some conclusions.
Section snippets
Non-Born–Oppenheimer Hooke–Coulomb three-particle model
The model for the present theoretical study consists of a system of three interacting particles, where particles 1 and 2 are identical charged particles with masses , and where particle 3 has a different mass as well as a different charge. The position of each particle is defined by a set of vectors arbitrarily chosen in the laboratory frame. The Hamiltonian associated with the model is the following:where the explicit form of the
One-particle densities for the three-particle Hooke–Coulomb system with different reference points
The one-particle density for identical particles 1 and 2 is computed by means of the following expression:where and the corresponding one-particle density operator isThe vectors denote the position of particle i from the standpoint of a given reference point : where the position of the reference point may, in general, be given by a linear combination of the vectors , satisfying the
Conclusions
In the present article, we have extended the analysis carried out previously for the Hooke-Calogero model for a three-particle system to a more realistic Hooke–Coulomb model. In contrast with the former model, for which we can obtain analytic solutions for all values of the masses and coupling constants, in the present case, we have to complement the analytic solutions with variational wave functions. However, in order to facilitate the treatment of this system, we have constructed these
Acknowledgements
E.V.L. would like to express his gratitude to SENESCYT of Ecuador for having given him the opportunity to participate in the Prometheus Program.
References (61)
The decoupling of electronic and nuclear motions in the isolated molecule Schrödinger Hamiltonian
Adv. Quantum Chem.
(2005)- et al.
Molecular structure in non-Born–Oppenheimer quantum mechanics
Chem. Phys. Lett.
(2004) - et al.
Comment on ‘Molecular structure in non-Born–Oppenheimer quantum mechanics’
Chem. Phys. Lett.
(2005) The density in density functional theory
J. Mol. Struct. THEOCHEM
(2010)- et al.
Nucleus-nucleus correlation function in non-Born–Oppenheimer molecular calculations: vibrationally excited states of HD+
Chem. Phys. Lett.
(2005) - et al.
The topological analysis of electronic charge densities: a reassessment of foundations
J. Mol. Struct. THEOCHEM
(2008) Chemical reaction dynamics beyond the born–oppenheimer approximation
Annu. Rev. Phys. Chem.
(1998)- et al.
Non-Born–Oppenheimer molecular dynamics
Acc. Chem.
(2006) - et al.
Breakdown of the adiabatic Born–Oppenheimer in graphene
Nat. Mat.
(2007) - et al.
Breakdown of the Born–Oppenheimer approximation in the reaction
Science
(2007)
Direct observation of Born–Oppenheimer approximation breakdown in carbon nanotubes
Nano Lett.
Excited electronic states and nonadiabatic effects in contemporary chemical dynamics
Acc. Chem. Res.
Time-resolved photoelectron spectroscopy of coupled electron-nuclear motion
J. Chem. Phys.
Probing and controlling non-Born–Oppenheimer dynamics in highly excited molecular ions
Nat. Phys.
The decoupling of electronic and nuclear motions in the isolated molecule Schrödinger Hamiltonian
Adv. Chem. Phys.
Nuclear dynamics in quantum clusters
Phys. Chem. Chem. Phys.
Convergence of experiment and theory on the pure vibrational spectrum of HeH+
Phys. Rev. Lett.
A non-Born–Oppenheimer self-consistent field method
J. Math. Chem.
Can chemistry be derived from quantum mechanics? Chemical dynamics and structure
J. Math. Chem.
Nuclear orbital plus molecular orbital theory: simultaneous determination of nuclear and electronic wave functions without Born–Oppenheimer approximation
Int. J. Quantum Chem.
Generalized electronic diabatic approach to structural similarity and the Hammond postulate
Int. J. Quantum Chem.
Nuclear and electronic motion in molecules can be computed separately
Phys. Today
Born–Oppenheimer approximation and beyond for time-dependent electronic processes
J. Chem. Phys.
Extension of density functional theory to nuclear orbital plus molecular orbital theory: self-consistent field calculations with the Colle-Salvetti electron-nucleus correlation functional
Bull. Chem. Soc. Jpn.
Nuclear quantum effects on the He2H+ complex with the nuclear molecular orbital approach
Int. J. Quantum Chem.
Molecular structure calculations: a unified quantum mechanical description of electrons and nuclei using explicitly correlated Gaussian functions and the global vector representation
J. Chem. Phys.
Analysis of electron-positron wavefunctions in the nuclear-electronic orbital framework
J. Chem. Phys.
Multicomponent density functional theory study of the interplay between electron-electron and electron-proton correlation
J. Chem. Phys.
Chemical bonding: from Lewis to atoms in molecules
J. Comput. Chem.
Cited by (8)
Non-deformed singular and non-singular exponential-type potentials
2017, Journal of Molecular ModelingHidden aspects of the Structural theory of chemistry: MC-QTAIM analysis reveals "alchemical" transformation from a triatomic to a diatomic structure
2015, Physical Chemistry Chemical Physics