The diagonal Born–Oppenheimer correction to molecular dynamical properties

Nucleus-electron correlation revising molecular bonding fingerprints from the exact wavefunction factorization

We present a novel theory and implementation for computing coupled electronic and quantal nuclear subsystems on a single potential energy surface, moving beyond the standard Born-Oppenheimer (BO) separation of nuclei and electrons. We formulate an exact self-consistent nucleus-electron embedding potential from the single product molecular wavefunction and demonstrate that the fundamental behavior of the correlated nucleus-electron can be computed for mean-field electrons that are responsive to a quantal anharmonic vibration of selected nuclei in a discrete variable representation. Geometric gauge choices are discussed and necessary for formulating energy invariant biorthogonal electronic equations. Our method is further applied to characterize vibrationally averaged molecular bonding properties of molecular energetics, bond lengths, and protonic and electron densities. Moreover, post-Hartree-Fock electron correlation can be conveniently computed on the basis of nucleus-electron coupled molecular orbitals, as demonstrated for correlated models of second-order Møllet-Plesset perturbation and full configuration interaction theories. Our approach not only accurately quantifies non-classical nucleus-electron couplings for revising molecular bonding properties but also provides an alternative time-independent approach for deploying non-BO molecular quantum chemistry.

Heavy Element Effects in the Diagonal Born-Oppenheimer Correction within a Relativistic Spin-Free Hamiltonian

Methodologies beyond the Born-Oppenheimer (BO) approximation are nowadays important to explain high precision spectroscopic measurements. Most previous evaluations of the BO correction are, however, focused on light-element molecules and based on a nonrelativistic Hamiltonian, so no information about the BO approximation (BOA) breakdown in heavy-element molecules is available. The present work is the first to investigate the BOA breakdown for the entire periodic table, by considering scalar relativistic effects in the Diagonal BO correction (DBOC). In closed shell atoms, the relativistic EDBOC scales as Z(1.25) and the nonrelativistic EDBOC scales as Z(1.17), where Z is the atomic number. Hence, we found that EDBOC becomes larger in heavy element atoms and molecules, and the relativistic EDBOC increases faster than nonrelativistic EDBOC. We have further investigated the DBOC effects on properties such as potential energy curves, spectroscopic parameters, and various energetic properties. The DBOC effects for these properties are mostly affected by the lightest atom in the molecule. Hence, in X2 or XAt molecule (X = H, Li, Na, K, Rb, and Cs) the effect of DBOC systematically decreases when X becomes heavier but in HX molecules, the effect of DBOC seems relatively similar among all the molecules.

Benchmark calculations of the complete configuration-interaction limit of Born–Oppenheimer diagonal corrections to the saddle points of isotopomers of the H+H2 reaction

We present a detailed ab initio study of the effect that the Born-Oppenheimer diagonal correction (BODC) has on the saddle-point properties of the H3 system and its isotopomers. Benchmark values are presented that are estimated to be within 0.1 cm(-1) of the complete configuration-interaction limit. We consider the basis set and correlation treatment requirements for accurate BODC calculations, and both are observed to be more favorable than for the Born-Oppenheimer energies. The BODC raises the H + H2 barrier height by 0.1532 kcal/mol and slightly narrows the barrier--with the imaginary frequency increasing by approximately 2%.