Dariusz Kȩdziera

Two-component relativistic methods for the heaviest elements

Different generalized Douglas-Kroll transformed Hamiltonians (DKn, n=1, 2,...,5) proposed recently by Hess et al. are investigated with respect to their performance in calculations of the spin-orbit splittings. The results are compared with those obtained in the exact infinite-order two-component (IOTC) formalism which is fully equivalent to the four-component Dirac approach. This is a comprehensive investigation of the ability of approximate DKn methods to correctly predict the spin-orbit splittings. On comparing the DKn results with the IOTC (Dirac) data one finds that the calculated spin-orbit splittings are systematically improved with the increasing order of the DK approximation. However, even the highest-order approximate two-component DK5 scheme shows certain deficiencies with respect to the treatment of the spin-orbit coupling terms in very heavy systems. The meaning of the removal of the spin-dependent terms in the so-called spin-free (scalar) relativistic methods for many-electron systems is discussed and a computational investigation of the performance of the spin-free DKn and IOTC methods for many-electron Hamiltonians is carried out. It is argued that the spin-free IOTC rather than the Dirac-Coulomb results give the appropriate reference for other spin-free schemes which are based on approximate two-component Hamiltonians. This is illustrated by calculations of spin-free DKn and IOTC total energies, r(-1) expectation values, ionization potentials, and electron affinities of heavy atomic systems.

Orbit-orbit relativistic corrections to the pure vibrational non-Born-Oppenheimer energies of H2

We report the derivation of the orbit-orbit relativistic correction for calculating pure vibrational states of diatomic molecular systems with sigma electrons within the framework that does not assume the Born-Oppenheimer (BO) approximation. The correction is calculated as the expectation value of the orbit-orbit interaction operator with the non-BO wave function expressed in terms of explicitly correlated Gaussian functions multiplied by even powers of the internuclear distance. With that we can now calculate the complete relativistic correction of the order of alpha(2) (where alpha=1/c). The new algorithm is applied to determine the full set of the rotationless vibrational levels and the corresponding transition frequencies of the H(2) molecule. The results are compared with the previous calculations, as well as with the frequencies obtained from the experimental spectra. The comparison shows the need to include corrections higher than second order in alpha to further improve the agreement between the theory and the experiment.

Electron affinity of Li7 calculated with the inclusion of nuclear motion and relativistic corrections

Explicitly correlated Gaussian functions have been used to perform very accurate variational calculations for the ground states of (7)Li and (7)Li(-). The nuclear motion has been explicitly included in the calculations (i.e., they have been done without assuming the Born-Oppenheimer (BO) approximation). An approach based on the analytical energy gradient calculated with respect to the Gaussian exponential parameters was employed. This led to a noticeable improvement of the previously determined variational upper bound to the nonrelativistic energy of Li(-). The Li energy obtained in the calculations matches those of the most accurate results obtained with Hylleraas functions. The finite-mass (non-BO) wave functions were used to calculate the alpha(2) relativistic corrections (alpha=1c). With those corrections and the alpha(3) and alpha(4) corrections taken from Pachucki and Komasa [J. Chem. Phys. 125, 204304 (2006)], the electron affinity (EA) of (7)Li was determined. It agrees very well with the most recent experimental EA.

Convergence of approximate two-component Hamiltonians: How far is the Dirac limit

A systematic elimination of the off-diagonal parts of the Dirac Hamiltonian is carried out in the spirit of the Douglas-Kroll [Ann. Phys. 82, 87 1974] approach and the recently proposed infinite-order two-component method. The present approach leads to a series of approximate two-component Hamiltonians which are exact through a certain order in the external potential. These Hamiltonians are used to study the convergence pattern of approximate two-component theories. It is shown that to achieve an acceptably high accuracy for low-lying one-electron levels in heavy and superheavy systems one needs to use approximate Hamiltonians of prohibitively high order in the external potential. One can conclude that the finite-order two-component Hamiltonians are of limited usefulness in accurate relativistic calculations for heavy and superheavy systems.

Relativistic corrections to the non-Born-Oppenheimer energies of the lowest singlet Rydberg states of He3 and He4

In this work the authors present an approach to calculate the leading-order relativistic corrections for ground and excited states of helium isotopomers. In the calculations they used variational wave functions expanded in terms of explicitly correlated Gaussians obtained without assuming the Born-Oppenheimer approximation.

Electronic spectrum of the confined auride ion

Several models of the confinement of the auride ion are discussed in the context of the colour of the tetramethylammonium auride crystals. These models correspond to either neutral or charged nearest neighbourhood of the auride ion as represented by a cluster of helium and He-like atoms, respectively. The models involving neutral helium clusters account for what is referred to as the confinement effect and simulate the Pauli repulsion between the auride ion and its surrounding. The charged clusters supply additional attractive potential which is shown to be indispensable for the interpretation of the large blue shift of the first absorption band of the auride ion. It is also found that the Pauli repulsion and the effect of the increase of the attractive potential due to the charged neighbourhood of the auride ion are not additive. The confinement models discussed in this paper are used to calculate the energy of the first electronic transition in the confined ion. The calculations have been carried out in the TDHF (RPA) and CASSCF/CASPT2 approximations with relativistic effects accounted for by using the Douglas-Kroll approximation. The charged confinement models are found to increase the importance of both the electron correlation and relativistic contributions to the energy of the first electronic transition in the encapsulated auride ion. The computed data are used to discuss the observed spectral transparency of the tetramethylammonium auride crystals.

Orbit-orbit relativistic correction calculated with all-electron molecular explicitly correlated Gaussians

An algorithm for calculating the first-order electronic orbit-orbit magnetic interaction correction for an electronic wave function expanded in terms of all-electron explicitly correlated molecular Gaussian (ECG) functions with shifted centers is derived and implemented. The algorithm is tested in calculations concerning the H2 molecule. It is also applied in calculations for LiH and H3+ molecular systems. The implementation completes our work on the leading relativistic correction for ECGs and paves the way for very accurate ECG calculations of ground and excited potential energy surfaces (PESs) of small molecules with two and more nuclei and two and more electrons, such as HeH-, H3+, HeH2+, and LiH2+. The PESs will be used to determine rovibrational spectra of the systems.