Few-electron Molecules Research Articles

Precision Measurements in Few-Electron Molecules: The Ionization Energy of Metastable ^{4}He_{2} and the First Rotational Interval of ^{4}He_{2}^{+}.

Molecular helium represents a benchmark system for testing abinitio calculations on few-electron molecules. We report on the determination of the adiabatic ionization energy of the a ^{3}Σ_{u}^{+} state of He_{2}, corresponding to the energy interval between the a ^{3}Σ_{u}^{+} (v^{''}=0, N^{''}=1) state of He_{2} and the X^{+} ^{2}Σ_{u}^{+} (v^{+}=0, N^{+}=1) state of He_{2}^{+}, and of the lowest rotational interval of He_{2}^{+}. These measurements rely on the excitation of metastable He_{2} molecules to high Rydberg states using frequency-comb-calibrated continuous-wave UV radiation in a counterpropagating laser-beam setup. The observed Rydberg states were extrapolated to their series limit using multichannel quantum-defect theory. The ionization energy of He_{2} (a ^{3}Σ_{u}^{+}) and the lowest rotational interval of He_{2}^{+} (X^{+} ^{2}Σ_{u}^{+}) are 34 301.207 002(23)±0.000 037_{syst} cm^{-1} and 70.937 589(23)±0.000 060_{syst} cm^{-1}, respectively.

Explicitly-correlated Gaussian geminals in electronic structure calculations

Explicitly correlated functions have been used since 1929, but initially only for two-electron systems. In 1960, Boys and Singer showed that if the correlating factor is of Gaussian form, many-electron integrals can be computed for general molecules. The capability of explicitly correlated Gaussian (ECG) functions to accurately describe many-electron atoms and molecules was demonstrated only in the early 1980s when Monkhorst, Zabolitzky and the present authors cast the many-body perturbation theory (MBPT) and coupled cluster (CC) equations as a system of integro-differential equations and developed techniques of solving these equations with two-electron ECG functions (Gaussian-type geminals, GTG). This work brought a new accuracy standard to MBPT/CC calculations. In 1985, Kutzelnigg suggested that the linear r 12 correlating factor can also be employed if n-electron integrals, n > 2, are factorised with the resolution of identity. Later, this factor was replaced by more general functions f (r 12), most often by , usually represented as linear combinations of Gaussian functions which makes the resulting approach (called F12) a special case of the original GTG expansion. The current state-of-art is that, for few-electron molecules, ECGs provide more accurate results than any other basis available, but for larger systems the F12 approach is the method of choice, giving significant improvements over orbital calculations.

An improved density matrix functional by physically motivated repulsive corrections

An improved density matrix functional [correction to Buijse and Baerends functional (BBC)] is proposed, in which a hierarchy of physically motivated repulsive corrections is employed to the strongly overbinding functional of Buijse and Baerends (BB). The first correction C1 restores the repulsive exchange-correlation (xc) interaction between electrons in weakly occupied natural orbitals (NOs) as it appears in the exact electron pair density rho(2) for the limiting two-electron case. The second correction C2 reduces the xc interaction of the BB functional between electrons in strongly occupied NOs to an exchange-type interaction. The third correction C3 employs a similar reduction for the interaction of the antibonding orbital of a dissociating molecular bond. In addition, C3 applies a selective cancellation of diagonal terms in the Coulomb and xc energies (not for the frontier orbitals). With these corrections, BBC still retains a correct description of strong nondynamical correlation for the dissociating electron pair bond. BBC greatly improves the quality of the BB potential energy curves for the prototype few-electron molecules and in several cases BBC reproduces very well the benchmark ab initio potential curves. The average error of the self-consistent correlation energies obtained with BBC3 for prototype atomic systems and molecular systems at the equilibrium geometry is only ca. 6%.