The Born-Oppenheimer Approximation Revisited

We present an improved Born-Oppenheimer method for the treatment of molecules. Our development is based on taking into account a portion of the kinetic energy which was inadvertently omitted by Born and Oppenheimer. Our theory replaces the set of standard atom-atom potentials by an evolution operator. In particular, in contrast to the traditional Born-Oppenheimer method our amended basis describes the so-called Wannier phenomenon which deforms potential surfaces, and traps one or more threshold electrons into an unstable equilibrium. Within our framework that trapping manifests itself as a Fresnel distribution. Finally we discover a π/2 phase jump in the evolution operator caused by the diffraction of a two-electron wave from a potential ridge.


Introduction
Since a century molecules are usually treated within the Born-Oppenheimer (hereafter shorty BO) approximation [1], that method has been also used to solve other nonseparable problems like many-electron atoms in terms of hypersherical coordinates [2] [3]. Usually the comunity believes that the success of the BO method rests on the existence of a small parameter like the electron to nucleus mass ratio. That is, however, not quite true. Important for the validity of the BO approximation is a smooth variation of the electronic Hamiltonian as function of the internuclear separation. That property rather than a small numerical parameter delivers a week non-adiabatic channel coupling. In atomic and molecular applications is that mostly the case due to smooth Coulomb interactions.
There are, however, situations where BO breaks really down. This is, for instance, the case at ionization/dissociation thresholds. Threshold ionization of H by electron impact was investigated long ago by Wannier [4] in terms of classical mechanics. Its quantum analog is entirely foreign to BO. Further, electronic transitions are not described by the BO standard basis. It is therefore the aim of this paper to develop a better basis which removes these shortcomings. At the time when the BO paper was written there was due to the experimental situation no need to treat such critical situations. However, today there are much more sophisticated experimental possibilities available. Therefore, there is a need for an improvement of the BO technique. Before we go into the new development we summarize shortly the standard BO procedure.
To this end we consider a non-rotating diatomic molecule with nuclear masses A m and B m , charges A Z and B Z , respectively. Its center-of-mass Hamiltonian reads ( ) where M is the reduced mass of the two nuclei; R is the internuclear separation; r summarizes all elctron coordinates; and ( ) ; h r R is the electronic part of the Hamiltonian. The latter contains the kinetic energy of all elctrons, all electron-nuclei interactions, and all electron-electron interactions. We consider here only light nuclei in nonrelativistic description, i.e. with pure Coulomb interactions.
The standard BO strategy is as follows. Because of the big mass of the nuclei they may be regarded to be at rest during the electrons motion. Therefore one constructs in a first step electronic eigenfunctions with nuclei at rest, i.e.
where the eigenvalue U(R) plays the role of an atom-atom potential.
BO use then a product form for the whole wave function, Mathematically, the above U plays the role of a separation parameter. Substitution of (2) into (1) delivers easily the nuclear wave equation given by The last step would be correct provided we were allowed to neglect the derivative terms R ∂Φ ∂ and BO were aware of that point. However they claim that the terms (5) are small because of the small parameter 1/M. That statement is indeed correct for the lower well separated channels; actually the BO channels are coupled by the derivarive terms (5). It is important to realize that the electron's motion occurs in a H. Klar Journal of Applied Mathematics and Physics moving frame even if that motion is slow. The motion occurs along the nuclear separation coordinate R. As novelty Section 2 takes the nuclear motion into account.

New Development
The above sketched BO derivation suffers from one weak point as may be seen below. It was assumed that the kinetic energy of the whole complex is equal to the sum of nuclear kinetic energy and electronic kinetic energy. That conclusion, however, is not true within the frame of a product wavefunction (3). The following trivial calculation shows The standard version of the BO procedure disregards the last two terms in (6).
Born and Oppenheimmer were aware of these terms which describe non-adiabatic couplings. The authors claim that these terms are small because of the mass ratio. It is, however, clear from (6) that all three portions carry the same mass ratio factor. Further it is not evident that the channel functions Φ depend only weakly on R. Below we will see that particularly the last term, a momentum-momentum coupling causes electronic transitions. Moreover, that coupling produces the Wannier phenomenon, i.e. an unusual phenomenon caused by equilibrium electronic configurations. Therefore it is the aim of this paper to take that momentum- To this end we modify the l.h.s. of the channel Equation (2) Note, however, that the radial function F is not yet known in (7). Its derivative we can, however, express by its logarithmic derivative corresponding to an and employ, following Sommerfeld [5], an Abstrahlungsbedingung 1 , i.e. we adress that condition to a specific collision channel labelled below by the index µ, and use above the µ-threshold is the wavenumber in the collision channel labelled by μ.
Along these lines we arrive at ( ) ( ) 1 Sommerfeld [5] used (9) in the context of emitted radiation. We do it here for particles. Journal of Applied Mathematics and Physics Because the operator in curly bracketts (10) depends on all coordinates there exists nolonger a R-dependent potential but (10) must be equal to a constant  .
As part of the Hamiltonian this constant has the dimension of an energy. For a finite radial function F we complete expression (10) to the equation where we employ  to be any energy level of the whole (N-1)-electron system.
The l.h.s. of (11) has the structure of a time-dependent wave equation where time t has been replaced by the separation R. Therefore (11) describes an evolution, for instance from the united atom limit to the limit of separated atoms. We rewrite (11) for convenience as Thus we intend to solve (12) with the Ansatz  is an eigenstate of the system in the united atom limit with energy  . The evolution operator satisfies then the same equation as (11), i.e.

( ) ( )
(14) may be regarded as a generalization of the Born-Oppenheimer equation (2) as may be seen as follows. Actually, if we confine ourselves to standing wave solutions we rediscover the BO eigenvalue problem as may be seen as follows. If the operator ( ) ; h r R depends only weakly on R we can solve our equation (11) with help of a phase, i.e. we put Substitution into (11) or (12) delivers immediately after identification of Entirely different is the situation for an ionizing collision with a charged collision partner at low energy, i.e. near threshold. We are then confronted with three charged particles in the continuum at nearly zero energy. We know from inelastic electron atom scattering that we are confronted in that situation with H. Klar the so-called Wannier [4] phenomenon. This is foreign to the standard BO approximation, but the next section will show that our evolution operator describes that critical situation as well.

Example: H(n ≥ 1) + Bare Nucleus
We illustrate below the power of our method treating as example an ionizing nucleus-atom collision. We use hydrogen as simplest target H(nl) with principal quantum number 1 n ≥ and energy levels To this end we give at first the location of the equilibrium point. The electron-proton separation in that point is whereas the electron-nucleus separation reads where Z is the charge of the colliding nucleus.
Since we intend to explore the threshold region, we determine the evolution operator near the nuclear charge center. The stable Coulomb attraction forces the elctron onto the internuclear axis, i.e. a collinear model of the whole system is justified. To describe the electron motion in the vicinity of that point we introduce at fixed separation R a coordinate x to describe departures from that point. To that end we write for the distances ( ) ( ) with 1 x ≤ . The electronic Hamiltonian reads therefore ( ) ( ) , with the total Coulomb potential given by Journal of Applied Mathematics and Physics ( ) ( ) For small values of x , i.e. near the nuclear charge-center this simplifies to ( ) with the net charge on the top of the (x,R)-potential ridge given by ( ) and the curvature across the ridge given by ( ) In (14) we write thus ( ) and put the ridge charge into the wavenumber now given by At threshold corresponding to n → ∞ we put The equation for the evolution reads therefore ( ) with the initial condition ( ) It is easy to find the solution in the Coulomb zone which extends at threshold to infinity. To this end we put where c is a constant to be determined below. Substitution of (31) into (30) confirms our Ansatz (31) with the condition The second solution combines incoming waves along R and x which represents a dominantly correlated shrinking solution which corresponds to Wannier's converging trajectory whose c-value is given by The algebraic structure of the evolution (31) resambles a Fresnel distribution [6] denoted in the mathematical literature by… The expression (35) converges to a Dirac delta distribution according to valid for x ∈  and real 0 a > . Therefore, an evolution starting from the united atom attracts the electron to its equilibrium point. Precisely, we get where we have used 1 R η = [6]. Moreover, we remark that the threshold evolution (31) delivers on the r.h.s. of (30) the quantity ( ) The Hermiticity of the evolution equation forbids an imaginary eigenvalue like (38). We read, therefore, our (38) in the form That interpretation tells us that the correlated three-body complex has picked up during the converging mode of motion a phase-jump of π/2 in consequence of the diffraction from the ridge. An effect similar to ours was reported by Fano [7] when he derived a phase jump due to the embedding of a Feshbach resonance into a continuum. Our quantity 3 2 2c R − may now be regarded as attractive potential curve which crosses all adiabatic BO curves. At these crossings we expect to occur electronic transitions. The evolution reported here plays therefore a katalytic role in electronic excitation processes.
The latter result is entirely foreign to the standard Born-Oppenheimer approximation. In the BO picture this unusual attraction is hidden in the infinity of strongly coupled adiabatic channels. It is remarkable that an electron is attracted to a statically unstable location. The result achieved in §3 represents a novel effect stemming from the diffraction of an electron wave from a many-body potential ridge. Here, we have demonstrated the effect in the simplest case of Journal of Applied Mathematics and Physics three bodies. It is straight forward to generalize the above effect to more electrons and multiply charged ions because; for instance, two or even more than two electrons can be trapped onto multidimensional potential ridges [8]. Finally, we remark that the present development did not use the smallness of any parameter like a mass ratio. Therefore it applies to electron/ion-atom scattering as well where the mass ratio 1/M is close to one.

Conflicts of Interest
The author declares no conflicts of interest regarding the publication of this paper.