Electronic structure calculations of small AI_n (n = 2 - 8) clusters

The electronic states of small AI_n (n = 2 - 8) clusters have been calculated with a relativistic ab-initio MOLCAO Dirac-Fock-Slater method using numerical atomic DFS wave-functions. The excitation energies were obtained from a ground state calculation of neutral clusters, and in addition from negative clusters charged by half an electron in order to account for part of the relaxation. These energies are compared with experimental photoelectron spectra.

Many-electron relativistic calculation and interpretation of atomic processes in time dependent heavy-ion scattering

The time dependence of a heavy-ion-atom collision system is solved via a set of coupled channel equations using energy eigenvalues and matrix elements from a self-consistent field relativistic molecular many-electron Dirac-Fock-Slater calculation. Within this independent particle model we give a full many-particle interpretation by performing a small number of single-particle calculations. First results for the P(b) curves for the Ne K-hole excitation for the systems F{^8+} - Ne and F{^6+} - Ne as examples are discussed.

Study of 3d - 4f and 3d - 5f quartet and doublet transitions of doubly excited three-electron ions

We present the first observation of optical transitions between doubly excited doublet states in the term systems N V, 0 VI and F VII. The spectra were produced by foil excitation of fast ion beams. The assignment of the spectral lines was made by comparison with the results of MCDP calculations along the isoelectronic sequence. The same method also led to the identification of two 3d - 4f quartet transitions in Mg X.

Further investigations on doubly excited four-electron ions: beam-foil experiment and MCDF theory

In continuation of our previous work on the quintet transitions 1s2s2p^2 ^5 P-1s2s2p3d ^5 P^0, ^5 D^0, results on other n = 2 - n' = 3 quintet transitions for elements N, 0 and F are presented. Assignments have been established by comparison with Multi-Configuration Dirac-Fock calculations. High spectral resolution on beam-foil spectroscopy was essential for the identification of most of the lines. For some of the quintet lines decay curves were measured, and the lifetimes extracted were found to be in reasonable agreement with MCDF calculations.

First results of experimental and theoretical investigations on sextet states of doubly-excited five-electron ions

In continuation of our previous work on doubly-excited ions with three and four electrons we present the first results on optical transitions in the term system of doubly-excited ions with five electrons. Transitions between such sextet states were identified in beam-foil spectra of the ions nitrogen, oxygen and fluorine. Assignments were first established by comparison with Multi-Configuration Dirac-Fock calculations. Later assignments were aided by Multi-Configuration Hartree-Fock calculations (see the contribution by G. Miecznik et al. in this issue). Decay curves were recorded for all six candidate lines. The lifetime results are compared to theoretical values which confirm most of the assignments qualitatively.

Relativistic atomic level calculations with nuclei carrying additional fractional charges

One-electron energy levels and wavelengths have been calculated for Na-like ions whose nuclei carry quarks with additional charges ±e/3, ±2e/3. The calculations are based on relativistic self-consistent field procedures. The deviations from experimental values exhibit regularities which allow an extrapolation for the wavelengths of 3s - 3p, 3s - 4p, 3p - 3d, and 3p - 4s transitions for the nuclear charge Z = 11± 1/3, ±2/3. A number of transitions are found in the region of visible light which could be used in an optical search for quark atoms.

Ground state correlation energy of the Be-sequence for Z = 4-20 in MCDF approximation

The ground state (J = 0) electronic correlation energy of the 4-electron Be-sequence is calculated in the Multi-Configuration Dirac-Fock approximation for Z = 4-20. The 4 electrons were distributed over the configurations arising from the 1s, 2s, 2p, 3s, 3p and 3d orbitals. Theoretical values obtained here are in good agreement with experimental correlation energies.

On the total energy of two-electron atoms

Using the Multi-Configuration Dirac-Fock (MCDF) method we calculate with 9 configuration state functions the correlation energy as well as the total energy of the lowest J = 0 ground state of all two-electron systems from H- to Thorium (Z = 90). A comparison with experimental data, which are available only in the low Z region, shows a very good agreement.

Estimation of the ground state correlation energy for isoelectronic series of 2 to 20 electrons

Correlation energies for all isoelectronic sequences of 2 to 20 electrons and Z = 2 to 25 are obtained by taking differences between theoretical total energies of Dirac-Fock calculations and experimental total energies. These are pure relativistic correlation energies because relativistic and QED effects are already taken care of. The theoretical as well as the experimental values are analysed critically in order to get values as accurate as possible. The correlation energies obtained show an essentially consistent behaviour from Z = 2 to 17. For Z > 17 inconsistencies occur indicating errors in the experimental values which become very large for Z > 25.

Transitions between quintet states of doubly excited four-electron ions

Following an earlier observation in F VI we identified the line pair 1s2s2p^2 {^5P}-1s2s2p3d {^5P^0} , {^5D^0} for the elements N, O, Mg, and tentatively for A1 and Si in beam-foil spectra. Assignment was established by comparison with Multi-Configuration Dirac-Fock calculations along the isoelectronic sequence. Using this method we also identified some quartet lines of lithium-like ions with Z > 10.

Ionization potentials and radii of atoms and ions of element 105 (unnilpentium) and ions of tantalum derived from multiconfiguration Dirac-Fock calculations

Multiconfiguration relativistic Dirac-Fock (MCDF) values were calculated for the first five ionization potentials of element 105 (unnilpentium) and of the other group 5b elements (V, Nb, and Ta). Some of these ionization potentials in electron volts (eV) with uncertainties are: 105(0), 7.4±0.4; 105(1 +), 16.3 ±0.2; 105(2 +), 24.3 ± 0.2; 105(3 + ), 34.9 ± 0.5; and 105(4 + ), 44.9 ± 0.1. Ionization potentials for Ta(1+), Ta(2 +), and Ta(3 + ) were also calculated. Accurate experimental values for these ionization potentials are not available. Ionic radii are presented for the 2+, 3+, 4 +, and 5+ ions of element 105 and for the + 2 ions of vanadium and niobium. These radii for vanadium and niobium are not available elsewhere. The ionization potentials and ionic radii obtained are used to determine some standard electrode potentials for element 105. Born-Haber cycles and a form of the Born equation for the Gibbs free energy of hydration of ions were used to calculate the standard electrode potentials.

Relativistic LCAO with Minimax Principle and New Balanced Basis Sets

Relativistic density functional theory is widely applied in molecular calculations with heavy atoms, where relativistic and correlation effects are on the same footing. Variational stability of the Dirac Hamiltonian is a very important field of research from the beginning of relativistic molecular calculations on, among efforts for accuracy, efficiency, and density functional formulation, etc. Approximations of one- or two-component methods and searching for suitable basis sets are two major means for good projection power against the negative continuum. The minimax two-component spinor linear combination of atomic orbitals (LCAO) is applied in the present work for both light and super-heavy one-electron systems, providing good approximations in the whole energy spectrum, being close to the benchmark minimax finite element method (FEM) values and without spurious and contaminated states, in contrast to the presence of these artifacts in the traditional four-component spinor LCAO. The variational stability assures that minimax LCAO is bounded from below. New balanced basis sets, kinetic and potential defect balanced (TVDB), following the minimax idea, are applied with the Dirac Hamiltonian. Its performance in the same super-heavy one-electron quasi-molecules shows also very good projection capability against variational collapse, as the minimax LCAO is taken as the best projection to compare with. The TVDB method has twice as many basis coefficients as four-component spinor LCAO, which becomes now linear and overcomes the disadvantage of great time-consumption in the minimax method. The calculation with both the TVDB method and the traditional LCAO method for the dimers with elements in group 11 of the periodic table investigates their difference. New bigger basis sets are constructed than in previous research, achieving high accuracy within the functionals involved. Their difference in total energy is much smaller than the basis incompleteness error, showing that the traditional four-spinor LCAO keeps enough projection power from the numerical atomic orbitals and is suitable in research on relativistic quantum chemistry. In scattering investigations for the same comparison purpose, the failure of the traditional LCAO method of providing a stable spectrum with increasing size of basis sets is contrasted to the TVDB method, which contains no spurious states already without pre-orthogonalization of basis sets. Keeping the same conditions including the accuracy of matrix elements shows that the variational instability prevails over the linear dependence of the basis sets. The success of the TVDB method manifests its capability not only in relativistic quantum chemistry but also for scattering and under the influence of strong external electronic and magnetic fields. The good accuracy in total energy with large basis sets and the good projection property encourage wider research on different molecules, with better functionals, and on small effects.

The mapping of the local contributions of Fermi and Coulomb correlation into intracule and extracule density distributions

The contributions of the correlated and uncorrelated components of the electron-pair density to atomic and molecular intracule I(r) and extracule E(R) densities and its Laplacian functions ∇2I(r) and ∇2E(R) are analyzed at the Hartree-Fock (HF) and configuration interaction (CI) levels of theory. The topologies of the uncorrelated components of these functions can be rationalized in terms of the corresponding one-electron densities. In contrast, by analyzing the correlated components of I(r) and E(R), namely, IC(r) and EC(R), the effect of electron Fermi and Coulomb correlation can be assessed at the HF and CI levels of theory. Moreover, the contribution of Coulomb correlation can be isolated by means of difference maps between IC(r) and EC(R) distributions calculated at the two levels of theory. As application examples, the He, Ne, and Ar atomic series, the C2-2, N2, O2+2 molecular series, and the C2H4 molecule have been investigated. For these atoms and molecules, it is found that Fermi correlation accounts for the main characteristics of IC(r) and EC(R), with Coulomb correlation increasing slightly the locality of these functions at the CI level of theory. Furthermore, IC(r), EC(R), and the associated Laplacian functions, reveal the short-ranged nature and high isotropy of Fermi and Coulomb correlation in atoms and molecules

A chemical Hamiltonian approach study of the basis set superposition error changes on electron densities and one- and two-center energy components

A comparision of the local effects of the basis set superposition error (BSSE) on the electron densities and energy components of three representative H-bonded complexes was carried out. The electron densities were obtained with Hartee-Fock and density functional theory versions of the chemical Hamiltonian approach (CHA) methodology. It was shown that the effects of the BSSE were common for all complexes studied. The electron density difference maps and the chemical energy component analysis (CECA) analysis confirmed that the local effects of the BSSE were different when diffuse functions were present in the calculations

Estimates of electronic coupling for excess electron transfer in DNA

Electronic coupling Vda is one of the key parameters that determine the rate of charge transfer through DNA. While there have been several computational studies of Vda for hole transfer, estimates of electronic couplings for excess electron transfer (ET) in DNA remain unavailable. In the paper, an efficient strategy is established for calculating the ET matrix elements between base pairs in a π stack. Two approaches are considered. First, we employ the diabatic-state (DS) method in which donor and acceptor are represented with radical anions of the canonical base pairs adenine-thymine (AT) and guanine-cytosine (GC). In this approach, similar values of Vda are obtained with the standard 6-31 G* and extended 6-31++ G* basis sets. Second, the electronic couplings are derived from lowest unoccupied molecular orbitals (LUMOs) of neutral systems by using the generalized Mulliken-Hush or fragment charge methods. Because the radical-anion states of AT and GC are well reproduced by LUMOs of the neutral base pairs calculated without diffuse functions, the estimated values of Vda are in good agreement with the couplings obtained for radical-anion states using the DS method. However, when the calculation of a neutral stack is carried out with diffuse functions, LUMOs of the system exhibit the dipole-bound character and cannot be used for estimating electronic couplings. Our calculations suggest that the ET matrix elements Vda for models containing intrastrand thymine and cytosine bases are essentially larger than the couplings in complexes with interstrand pyrimidine bases. The matrix elements for excess electron transfer are found to be considerably smaller than the corresponding values for hole transfer and to be very responsive to structural changes in a DNA stack