[Paper Review] Electronic non-adiabatic states

Nikitas I. Gidopoulos, E. K. U. Gross|arXiv (Cornell University)|Feb 17, 2005

Mechanical and Optical Resonators3 references3 citations

TL;DR

This paper proposes a formally exact framework for electronic non-adiabatic states by expressing the full molecular wavefunction as an exact product of a conditional electronic wavefunction Φ_R(r) and a nuclear wavefunction X(R). It derives self-consistent equations for both, introducing non-adiabatic couplings via an optimized effective potential (OEP) that captures electron-nuclear correlation beyond the Born-Oppenheimer approximation, with successful numerical validation on H₂⁺.

ABSTRACT

A novel treatment of non-adiabatic couplings is proposed. The derivation starts from the long-known, but not well-known, fact that the wave function of the complete system of elctrons and nuclei can be written, without approximation, as a Born-Oppenheimer-type product of a nuclear wavefunction, X(R), and an electronic one, Phi_R(r), which depends parametrically on the nuclear configuration R. From the variational principle we deduce formally exact equations for Phi_R(r) and X(R). The algebraic structure of the exact nuclear equation coincides with the corresponding one in the adiabatic approximation. The electronic equation, however, contains terms not appearing in the adiabatic case, which couple the electronic and the nuclear wavefunctions and account for the electron-nuclear correlation beyond the Born-Oppenheimer level. It is proposed that these terms can be incorporated using an optimized local effective potential.

Motivation & Objective

• To develop a rigorous, non-adiabatic extension of the Born-Oppenheimer approximation that preserves its conceptual simplicity while improving accuracy.
• To address the breakdown of the adiabatic approximation in systems with strong electron-nuclear correlation or near-degeneracy.
• To derive formally exact equations for electronic and nuclear wavefunctions that include non-adiabatic couplings through a variational principle.
• To formulate a practical method for incorporating non-adiabatic effects using an optimized effective potential (OEP) within existing electronic structure frameworks.

Proposed method

• Starts from the exact factorization of the molecular wavefunction as Ψ(r,R) = Φ_R(r) X(R), where Φ_R(r) is the conditional electronic wavefunction and X(R) the nuclear marginal wavefunction.
• Derives exact variational equations for Φ_R(r) and X(R) using the molecular Hamiltonian, with the nuclear equation retaining the same algebraic structure as in the adiabatic case.
• The electronic equation includes additional terms coupling Φ_R and X(R), arising from non-adiabatic effects, which are expressed via a Green's function operator Ĝ_R.
• Introduces an optimized effective potential (OEP) formalism to represent the non-adiabatic coupling terms as an effective potential V_R(r), enabling practical implementation.
• Uses the Krieger-Li-Iafrate (KLI) approximation to simplify the OEP equation for practical computation, especially in Kohn-Sham DFT frameworks.
• Validates the method numerically on H₂⁺ by comparing the non-adiabatic OEP solution to the standard Born-Oppenheimer result, showing convergence to the correct wavefunction.

Experimental results

Research questions

• RQ1 Can the non-adiabatic coupling between electrons and nuclei be systematically incorporated into the electronic structure framework while preserving the simplicity of the adiabatic approximation?
• RQ2 What is the exact form of the electronic and nuclear equations that include non-adiabatic effects beyond the Born-Oppenheimer level?
• RQ3 Can the non-adiabatic coupling terms be represented as an effective potential within a self-consistent framework?
• RQ4 How can the optimized effective potential (OEP) method be adapted to include electron-nuclear correlation effects in a way compatible with existing electronic structure codes?
• RQ5 Does the proposed method recover the correct molecular wavefunction in a system like H₂⁺ where the adiabatic approximation is known to be accurate?

Key findings

• The exact factorization Ψ(r,R) = Φ_R(r) X(R) provides a rigorous foundation for non-adiabatic theory, with Φ_R(r) and X(R) satisfying formally exact variational equations.
• The electronic equation contains additional terms involving the nuclear wavefunction gradient and the Green's function Ĝ_R, which account for non-adiabatic coupling.
• The non-adiabatic coupling terms can be mapped exactly onto an effective potential V_R(r) via the OEP formalism, enabling practical implementation.
• The KLI approximation allows the OEP equation to be solved within standard Kohn-Sham DFT self-consistency loops, making the method computationally feasible.
• Numerical tests on H₂⁺ show that the non-adiabatic OEP solution converges to the correct wavefunction obtained via standard Born-Oppenheimer separation, validating the method.
• The method successfully captures electron-nuclear correlation beyond the adiabatic level without introducing approximations to the factorization.