Universal Gaussian basis set for accurate ab initio /P relat ivistic Dirac-Fock calculations

Atomic Electronic Structure Calculations with Hermite Interpolating Polynomials

We have recently described the implementation of atomic electronic structure calculations within the finite element method with numerical radial basis functions of the form χμ(r) = r-1Bμ(r), where high-order Lagrange interpolating polynomials (LIPs) were used as the shape functions Bμ(r). In this work, we discuss how χμ(r) can be evaluated in a stable manner at small r and also revisit the choice of the shape functions Bμ(r). Three kinds of shape functions are considered: in addition to the C 0 continuous LIPs, we consider the analytical implementation of first-order Hermite interpolating polynomials (HIPs) that are C 1 continuous, as well as numerical implementations of n-th order (C n continuous) HIPs that are expressed in terms of an underlying high-order LIP basis. Furnished with the new implementation, we demonstrate that the first-order HIPs are reliable even with large numbers of nodes and that they also work with nonuniform element grids, affording even better results in atomic electronic structure calculations than LIPs with the same total number of basis functions. We demonstrate that discontinuities can be observed in the spin-σ local kinetic energy τσ in small LIP basis sets, while HIP basis sets do not suffer from such issues; however, either set can be used to reach the complete basis set limit with smooth τσ. Moreover, we discuss the implications of HIPs on calculations with meta-GGA functionals with a number of recent meta-GGA functionals, and we find most Minnesota functionals to be ill-behaved. We also examine the potential usefulness of the explicit control over the derivative in HIPs for forming numerical atomic orbital basis sets, but we find that confining potentials are still likely a better option.

Douglas-Kroll and infinite order two-component transformations of Dirac-Fock operator

We extended the conventional Douglas-Kroll (DK) and infinite order two-component (IOTC) methods to a technique applicable to Fock matrices, called extended DK (EDK) and extended IOTC (EIOTC), respectively. First, we defined a strategy to divide the Dirac-Fock operator into zero- and first-order terms. We then demonstrated that the first-order extended DK transformation, which is the Foldy-Wouthuysen transformation for the zero-order term, as well as the second- and third-order EDK and EIOTC, could be well defined. The EDK- and EIOTC-transformed Fock matrix, kinetic energy operator, nuclear attraction operator, and density matrix were derived. These equations were numerically evaluated, and it was found that these methods were accurate. In particular, EIOTC was consistent with the four-component approach. Four-component and extended two-component calculations are more expensive than non-relativistic calculations due to small-component-type two-electron integrals. We developed a new approximation formula, RIS-V, for small-component-type two-electron integrals, including the spin-orbit interaction between electrons. These results suggest that the RIS-V formula effectively accelerates the four-component and extended two-component methods.

Relativistic Fock Space Coupled Cluster Method for Many-Electron Systems: Non-Perturbative Account for Connected Triple Excitations

The Fock space relativistic coupled cluster method (FS-RCC) is one of the most promising tools of electronic structure modeling for atomic and molecular systems containing heavy nuclei. Until recently, capabilities of the FS-RCC method were severely restricted by the fact that only single and double excitations in the exponential parametrization of the wave operator were considered. We report the design and the first computer implementation of FS-RCC schemes with full and simplified non-perturbative account for triple excitations in the cluster operator. Numerical stability of the new computational scheme and thus its applicability to a wide variety of molecular electronic states is ensured using the dynamic shift technique combined with the extrapolation to zero-shift limit. Pilot applications to atomic (Tl, Pb) and molecular (TlH) systems reported in the paper indicate that the breakthrough in accuracy and predictive power of the electronic structure calculations for heavy-element compounds can be achieved. Moreover, the described approach can provide a firm basis for high-precision modeling of heavy molecular systems with several open shells, including actinide compounds.

Detection of the 5p - 4f orbital crossing and its optical clock transition in Pr9

Recent theoretical works have proposed atomic clocks based on narrow optical transitions in highly charged ions. The most interesting candidates for searches of physics beyond the Standard Model are those which occur at rare orbital crossings where the shell structure of the periodic table is reordered. There are only three such crossings expected to be accessible in highly charged ions, and hitherto none have been observed as both experiment and theory have proven difficult. In this work we observe an orbital crossing in a system chosen to be tractable from both sides: Pr\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}^{9+}$$\end{document}9+. We present electron beam ion trap measurements of its spectra, including the inter-configuration lines that reveal the sought-after crossing. With state-of-the-art calculations we show that the proposed nHz-wide clock line has a very high sensitivity to variation of the fine-structure constant, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}α, and violation of local Lorentz invariance; and has extremely low sensitivity to external perturbations.

Atomic clocks are based on the frequency of optical transitions and offer high precision. Here the authors demonstrate a configuration crossing in the highly charged ion praseodymium (Pr\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}^{9+}$$\end{document}9+) and determine the frequency of a potential reference transition for a highly charged ion clock.

Identification of the predicted 5s-4f level crossing optical lines with applications to metrology and searches for the variation of fundamental constants

We measure optical spectra of Nd-like W, Re, Os, Ir, and Pt ions of particular interest for studies of a possibly varying fine-structure constant. Exploiting characteristic energy scalings we identify the strongest lines, confirm the predicted 5s-4f level crossing, and benchmark advanced calculations. We infer two possible values for optical M2/E3 and E1 transitions in Ir^{17+} that have the highest predicted sensitivity to a variation of the fine-structure constant among stable atomic systems. Furthermore, we determine the energies of proposed frequency standards in Hf^{12+} and W^{14+}.

Reappraisal of nuclear quadrupole moments of atomic halogens via relativistic coupled cluster linear response theory for the ionization process

The coupled cluster based linear response theory (CCLRT) with four-component relativistic spinors is employed to compute the electric field gradients (EFG) of (35)Cl, (79)Br, and (127)I nuclei. The EFGs resulting from these calculations are combined with experimental nuclear quadrupole coupling constants (NQCC) to determine the nuclear quadrupole moments (NQM), Q of the halide nuclei. Our estimated NQMs [(35)Cl = -81.12 mb, (79)Br = 307.98 mb, and (127)I = -688.22 mb] agree well with the new atomic values [(35)Cl = -81.1(1.2), (79)Br = 302(5), and (127)I = -680(10) mb] obtained via Fock space multireference coupled cluster method with the Dirac-Coulomb-Breit Hamiltonian. Although our estimated Q((79)Br) value deviates from the accepted reference value of 313(3) mb, it agrees well with the recently recommended value, Q((79)Br) = 308.7(20) mb. Good agreement with current reference data indicates the accuracy of the proposed value for these halogen nuclei and lends credence to the results obtained via CCLRT approach. The electron affinities yielded by this method with no extra cost are also in good agreement with experimental values, which bolster our belief that the NQMs values for halogen nuclei derived here are reliable.

Dynamic phase transition from localized to spatiotemporal chaos in coupled circle map with feedback

We investigate coupled circle maps in the presence of feedback and explore various dynamical phases observed in this system of coupled high dimensional maps. We observe an interesting transition from localized chaos to spatiotemporal chaos. We study this transition as a dynamic phase transition. We observe that persistence acts as an excellent quantifier to describe this transition. Taking the location of the fixed point of circle map (which does not change with feedback) as a reference point, we compute a number of sites which have been greater than (less than) the fixed point until time t. Though local dynamics is high dimensional in this case, this definition of persistence which tracks a single variable is an excellent quantifier for this transition. In most cases, we also obtain a well defined persistence exponent at the critical point and observe conventional scaling as seen in second order phase transitions. This indicates that persistence could work as a good order parameter for transitions from fully or partially arrested phase. We also give an explanation of gaps in eigenvalue spectrum of the Jacobian of localized state.

Relativistic double-zeta, triple-zeta, and quadruple-zeta basis sets for the 4s, 5s, 6s, and 7s elements

Relativistic basis sets of double-zeta, triple-zeta, and quadruple-zeta quality have been optimized in Dirac-Hartree-Fock calculations for the 4s, 5s, 6s, and 7s elements: K, Ca, Rb, Sr, Cs, Ba, Fr, and Ra. The basis sets include SCF exponents for the occupied spinors and for the np shell, exponents of correlating and polarizing functions for the (n - 1) shell and correlating functions for the (n - 2) shell. For the group 2 elements, correlating functions are given for the ns and np shells, whereas for the group 1 elements, functions for polarization of the ns shell are provided. A finite nuclear size was used in all optimizations. Prescriptions are given for constructing contracted basis sets by addition of primitives to the SCF occupied functions.

A Fock space coupled cluster study on the electronic structure of the UO2, UO2+, U4+, and U5+ species

The ground and excited states of the UO(2) molecule have been studied using a Dirac-Coulomb intermediate Hamiltonian Fock-space coupled cluster approach (DC-IHFSCC). This method is unique in describing dynamic and nondynamic correlation energies at relatively low computational cost. Spin-orbit coupling effects have been fully included by utilizing the four-component Dirac-Coulomb Hamiltonian from the outset. Complementary calculations on the ionized systems UO(2) (+) and UO(2) (2+) as well as on the ions U(4+) and U(5+) were performed to assess the accuracy of this method. The latter calculations improve upon previously published theoretical work. Our calculations confirm the assignment of the ground state of the UO(2) molecule as a (3)Phi(2u) state that arises from the 5f(1)7s(1) configuration. The first state from the 5f(2) configuration is found above 10,000 cm(-1), whereas the first state from the 5f(1)6d(1) configuration is found at 5,047 cm(-1).

Nuclear quadrupole moment of Au197 from high-accuracy atomic calculations

The electric field gradient (EFG) at the gold nucleus is calculated using a finite field approach, to make the extraction of the nuclear quadrupole moment Q from experimental nuclear quadrupole coupling constants possible. The four-component Dirac-Coulomb Hamiltonian serves as the framework, 51 of the 79 electrons are correlated by the relativistic Fock-space coupled cluster method with single and double excitations, and the contribution of the Gaunt term, the main part of the Breit interaction, is evaluated. Large basis sets (up to 26s22p18d12f8g5h uncontracted Gaussians) are employed. Energy splittings of the 2D5/2 and 2D3/2 levels, rather than level shifts, are used to extract the EFG constants, as the former remain linear with Q up to 10(-5) a.u., whereas the latter display significant nonlinearity even at Q=10(-8) a.u. Larger Q values lead to larger energy changes and better precision. Excellent agreement (0.1%) is obtained between Q values derived from 2D5/2 and 2D3/2 data. Systematic errors connected with neglecting triple and higher excitations, truncating the basis and orbital active space, and approximating the Gaunt contribution are evaluated. The final value of Q(197Au) is 521(7) mb. It is lower than the muonic 547(16) mb and agrees within error bounds with the recent value of 510(15) mb obtained from molecular calculations.

Relativistic correlating basis sets for actinide atoms from90Th to103Lr

For 14 actinide atoms from (90)Th to (103)Lr, contracted Gaussian-type function sets are developed for the description of correlations of the 5f, 6d, and 7s electrons. Basis sets for the 6d orbitals are also prepared, since the orbitals are important in molecular environments despite their vacancy in the ground state of some actinides. A segmented contraction scheme is employed for the compactness and efficiency. Contraction coefficients and exponents are so determined as to minimize the deviation from accurate natural orbitals of the lowest term arising from the 5f(n-1)6d(1)7s(2) configuration. The spin-free relativistic effects are considered through the third-order Douglas-Kroll approximation. To test the present correlating sets, all-electron calculations are performed on the ground state of (90)ThO molecule. The calculated spectroscopic constants are in excellent agreement with experimental values.

Exact decoupling of the Dirac Hamiltonian. II. The generalized Douglas-Kroll-Hess transformation up to arbitrary order

In order to achieve exact decoupling of the Dirac Hamiltonian within a unitary transformation scheme, we have discussed in part I of this series that either a purely numerical iterative technique (the Barysz-Sadlej-Snijders method) or a stepwise analytic approach (the Douglas-Kroll-Hess method) are possible. For the evaluation of Douglas-Kroll-Hess Hamiltonians up to a pre-defined order it was shown that a symbolic scheme has to be employed. In this work, an algorithm for this analytic derivation of Douglas-Kroll-Hess Hamiltonians up to any arbitrary order in the external potential is presented. We discuss how an estimate for the necessary order for exact decoupling (within machine precision) for a given system can be determined from the convergence behavior of the Douglas-Kroll-Hess expansion prior to a quantum chemical calculation. Once this maximum order has been accomplished, the spectrum of the positive-energy part of the decoupled Hamiltonian, e.g., for electronic bound states, cannot be distinguished from the corresponding part of the spectrum of the Dirac operator. An efficient scalar-relativistic implementation of the symbolic operations for the evaluation of the positive-energy part of the block-diagonal Hamiltonian is presented, and its accuracy is tested for ground-state energies of one-electron ions over the whole periodic table. Furthermore, the first many-electron calculations employing sixth up to fourteenth order DKH Hamiltonians are presented.

Retarded boson-fermion interaction in atomic systems

The retarded interaction between an electron and a spin-0 nucleus is derived from electrodynamical perturbation theory. The contribution of retardation at order v(2)c(2) mimics the Breit interaction [Phys. Rev. 34, 553 (1929); 36, 388 (1930); 39, 616 (1932)] with the Dirac matrix alpha(2) being replaced by p(2)m(2)c where p(2) is the linear momentum operator for the nucleus. An effective one-electron retardation operator is obtained in relative coordinates, and this can be used through all orders in perturbation theory without any problem of infinite degeneracy. A few steps of unitary transformation lead to the nonrelativistic limit. The leading terms in retardation corrections to energy are of order (m(e)m(n))alpha(2)Z(4)(alpha(2)m(e)c(2)). The implications for atomic systems are discussed.

Dissociation energy of ekaplutonium fluoride E126F: The first diatomic with molecular spinors consisting of g atomic spinors

Our ab initio all-electron fully relativistic Dirac-Fock (DF) and nonrelativistic (NR) Hartree-Fock (HF) self-consistent field (SCF) calculations predict the superheavy diatomic ekaplutonium fluoride E126F to be bound with the calculated dissociation energy of 7.44 and 10.46 eV at the predicted E126-F bond lengths of 2.03 and 2.18 Angstroms, respectively. The antibinding effects of relativity to the dissociation energy of E126F are approximately 3 eV. The predicted dissociation energy with both our NR HF and relativistic DF SCF wave functions is fairly large and is comparable to that for very stable diatomics. This is the first case, where in a diatomic, an atom has g orbital (l = 4) occupied in its ground state electronic configuration and such superheavy diatomics would have occupied molecular spinors (orbitals) consisting of g atomic spinors (orbitals). This opens up a whole new field of chemistry where g atomic spinors (orbitals) may be involved in electronic structure and chemical bonding of systems of superheavy elements with Z> or =122.

Electronic structure and prediction of atomization energy of naked homoleptic uranium hexacarbonyl U(CO)6

Our ab initio all-electron fully relativistic Dirac-Fock and nonrelativistic Hartree-Fock self-consistent field (SCF) calculations predict the octahedral (Oh) uranium hexacarbonyl U(CO)6 to be bound with the calculated atomization energy of 49.84 and 48.76 eV at the predicted U-C bond lengths (assuming the C-O bond distance fixed at 1.17 A) of 2.53 and 2.63 A, respectively. Moreover, our all-electron fully relativistic Dirac-Fock SCF calculations predict U(CO)6 to be lower in energy by 3.90 eV with respect to dissociation into U plus six CO ligands. We predict U(CO)6 (Oh) to be very stable in view of our predicted large atomization energy (approximately 49 eV) and stability (approximately 4 eV) with respect to dissociation into U plus six CO molecules. Innovative techniques should be devised for the synthesis of uranium hexacarbonyl since the usual synthetic methods have failed so far for this naked actinide hexacarbonyl.

Relativistic correlating basis sets for lanthanide atoms from Ce to Lu

Contracted Gaussian-type function (CGTF) sets for the description of the 4f subshell correlation and of the 6s and 5d subshell correlation are developed for lanthanide atoms from Ce to Yb. Also prepared are basis sets for the 5d orbitals, which are vacant in the ground states of most lanthanide atoms but are essential in molecular environments. In addition, correlating CGTF sets for the 4f subshell correlation are supplemented for the Lu atom. A segmented contraction scheme is employed for their compactness and efficiency. Contraction coefficients and exponents are determined by minimizing the deviation from accurate natural orbitals generated from configuration interaction calculations that include relativistic effects through the third-order Douglas-Kroll approximation. All-electron and model core potential calculations with the present correlating sets are performed on the ground state of the diatomic CeO molecule. The calculated spectroscopic constants are in good agreement with experimental values.

Relativistic calculation of indirect NMR spin-spin couplings using the Douglas-Kroll-Hess approximation

We have employed the Douglas-Kroll-Hess approximation to derive the perturbative Hamiltonians involved in the calculation of NMR spin-spin couplings in molecules containing heavy elements. We have applied this two-component quasirelativistic approach using finite perturbation theory in combination with a generalized Kohn-Sham code that includes the spin-orbit interaction self-consistently and works with Hartree-Fock and both pure and hybrid density functionals. We present numerical results for one-bond spin-spin couplings in the series of tetrahydrides CH(4), SiH(4), GeH(4), and SnH(4). Our two-component Hartree-Fock results are in good agreement with four-component Dirac-Hartree-Fock calculations, although a density-functional treatment better reproduces the available experimental data.

Energy band gaps and lattice parameters evaluated with the Heyd-Scuseria-Ernzerhof screened hybrid functional

This work assesses the Heyd-Scuseria-Ernzerhof (HSE) screened Coulomb hybrid density functional for the prediction of lattice constants and band gaps using a set of 40 simple and binary semiconductors. An extensive analysis of both basis set and relativistic effects is given. Results are compared with established pure density functionals. For lattice constants, HSE outperforms local spin-density approximation (LSDA) with a mean absolute error (MAE) of 0.037 A for HSE vs 0.047 A for LSDA. For this specific test set, all pure functionals tested produce MAEs for band gaps of 1.0-1.3 eV, consistent with the very well-known fact that pure functionals severely underestimate this property. On the other hand, HSE yields a MAE smaller than 0.3 eV. Importantly, HSE correctly predicts semiconducting behavior in systems where pure functionals erroneously predict a metal, such as, for instance, Ge. The short-range nature of the exchange integrals involved in HSE calculations makes their computation notably faster than regular hybrid functionals. The current results, paired with earlier work, suggest that HSE is a fast and accurate alternative to established density functionals, especially for solid state calculations.