Search

Relativistic basis sets of double-zeta, triple-zeta, and quadruple-zeta quality have been optimized for the lanthanide elements La-Lu. The basis sets include SCF exponents for the occupied spinors and for the 6p shell, exponents of correlating functions for the valence shells (4f, 5d and 6s) and the outer core shells (4d, 5s and 5p), and diffuse functions, including functions for dipole polarization of the 4f shell. A finite nuclear size was used in all optimizations. The basis sets are illustrated by calculations on YbF. Prescriptions are given for constructing contracted basis sets. The basis sets are available as an internet archive and from the Dirac program web site, . [ABSTRACT FROM AUTHOR]

In this work a scheme for constructing systematic sequences of relativistic SCF basis sets at a reasonable computational cost is presented and applied to atoms of the s- and p-block up to Xe. This scheme, which couples simplex optimization and the use of geometric series given by four-term polynomial expressions for the logarithm of the exponents, allows for the construction of basis sets that exhibit very regular patterns of convergence to the numerical reference values of atomic total energies, spinor energies and radial expectation values. This regularity, together with the broad range of basis set sizes presented, enables these sets to be used as building blocks for basis sets applicable in both routine and benchmark relativistic calculations on atomic and molecular systems. [ABSTRACT FROM AUTHOR]

The electrostatic contribution to the Mössbauer isomer shift of mercury for the series HgF ( n = 1, 2, 4) with respect to the neutral atom has been investigated in the framework of four- and two-component relativistic theory. Replacing the integration of the electron density over the nuclear volume by the contact density (that is, the electron density at the nucleus) leads to a 10% overestimation of the isomer shift. The systematic nature of this error suggests that it can be incorporated into a correction factor, thus justifying the use of the contact density for the calculation of the Mössbauer isomer shift. The performance of a large selection of density functionals for the calculation of contact densities has been assessed by comparing with finite-field four-component relativistic coupled-cluster with single and double and perturbative triple excitations [CCSD(T)] calculations. For the absolute contact density of the mercury atom, the Density Functional Theory (DFT) calculations are in error by about 0.5%, a result that must be judged against the observation that the change in contact density along the series HgF ( n = 1, 2, 4), relevant for the isomer shift, is on the order of 50 ppm with respect to absolute densities. Contrary to previous studies of the Fe isomer shift (F Neese, Inorg Chim Acta 332:181, 2002), for mercury, DFT is not able to reproduce the trends in the isomer shift provided by reference data, in our case CCSD(T) calculations, notably the non-monotonous decrease in the contact density along the series HgF ( n = 1, 2, 4). Projection analysis shows the expected reduction of the 6 s population at the mercury center with an increasing number of ligands, but also brings into light an opposing effect, namely the increasing polarization of the 6 s orbital due to increasing effective charge of the mercury atom, which explains the non-monotonous behavior of the contact density along the series. The same analysis shows increasing covalent contributions to bonding along the series with the effective charge of the mercury atom reaching a maximum of around +2 for HgF at the DFT level, far from the formal charge +4 suggested by the oxidation state of this recently observed species. Whereas the geometries for the linear HgF and square-planar HgF molecules were taken from previous computational studies, we optimized the equilibrium distance of HgF at the four-component Fock-space CCSD/aug-cc-pVQZ level, giving spectroscopic constants r = 2.007 Å and ω = 513.5 cm. [ABSTRACT FROM AUTHOR]