Relativistic double-zeta, triple-zeta, and quadruple-zeta basis sets for the 7p elements, with atomic and molecular applications

Relativistic and quantum electrodynamics effects on NMR shielding tensors of TlX (X = H, F, Cl, Br, I, At) molecules

The results of relativistic calculations of nuclear magnetic resonance shielding tensors (σ) for the thallium monocation (Tl+), thallium hydride (TlH), and thallium halides (TlF, TlCl, TlBr, TlI, and TlAt) are presented as obtained within a four-component polarization propagator formalism and a two-component linear response approach within the zeroth-order regular approximation. In addition to a detailed analysis of relativistic effects performed in this work, some quantum electrodynamical (QED) effects on those nuclear magnetic resonance shieldings and other small contributions are estimated. A strong dependence of σ(Tl) on the bonding partner is found, together with a very weak dependence of QED effects with them. In order to explain the trends observed, the excitation patterns associated with relativistic ee (or paramagnetic-like) and pp (or diamagnetic-like) contributions to σ are analyzed. For this purpose, the electronic spin-free and spin-dependent contributions are separated within the two-component zeroth-order regular approximation, and the influence of spin-orbit coupling on involved molecular orbitals is studied, which allows for a thorough understanding of the underlying mechanisms.

Dirac-Coulomb-Breit Molecular Mean-Field Exact-Two-Component Relativistic Equation-of-Motion Coupled-Cluster Theory

We present a relativistic equation-of-motion coupled-cluster with single and double excitation formalism within the exact two-component framework (X2C-EOM-CCSD), where both scalar relativistic effects and spin-orbit coupling are variationally included at the reference level. Three different molecular mean-field treatments of relativistic corrections, including the one-electron, Dirac-Coulomb, and Dirac-Coulomb-Breit Hamiltonian, are considered in this work. Benchmark calculations include atomic excitations and fine-structure splittings arising from spin-orbit coupling. Comparison with experimental values and relativistic time-dependent density functional theory is also carried out. The computation of the oscillator strength using the relativistic X2C-EOM-CCSD approach allows for studies of spin-orbit-driven processes, such as the spontaneous phosphorescence lifetime.

Error of relativistic effective core potentials for closed-shell diatomic molecules of p-block heavy and superheavy elements in DFT and TDDFT calculations

Pseudopotentials (PP) are extensively used in electronic structure calculations, particularly for molecules containing heavy elements. Parameters in PPs are mainly determined from ab initio results, and errors of such PPs in density functional theory (DFT) calculations have been studied previously. However, PP errors on results with spin-orbit coupling and those in time-dependent DFT (TDDFT) calculations have not been reported previously. In this work, we investigate the error of the small-core energy-consistent Stuttgart/Koln pseudopotentials in DFT and TDDFT calculations with and without spin-orbit coupling. Ground state bond lengths, harmonic frequencies, dissociation energies, and vertical excitation energies for a series of closed-shell diatomic heavy and superheavy p-block molecules are calculated using several popular exchange-correlation functionals. PP errors are estimated by comparing with results using the all-electron Dirac-Coulomb (-Gaunt) Hamiltonian. Our results show that the difference between ground state properties and most excitation energies in scalar-relativistic calculations with the PP and those of all-electron calculations is quite small. This difference becomes somewhat larger when spin-orbit coupling (SOC) is present, especially for properties that are affected by SOC to some extent. In addition, the errors of the PPs are insensitive to the employed exchange-correlation functionals in most cases. Our results indicate that reasonable DFT and TDDFT results can be obtained using the small-core energy-consistent Stuttgart/Koln pseudopotentials for heavy and super-heavy p-block molecules.

Novel Computational Chemistry Infrastructure for Simulating Astatide in Water: From Basis Sets to Force Fields Using Particle Swarm Optimization

Using the example of astatine, the heaviest naturally occurring halogen whose isotope At-211 has promising medical applications, we propose a new infrastructure for large-scale computational models of heavy elements with strong relativistic effects. In particular, we focus on developing an accurate force field for At- in water based on reliable relativistic density functional theory (DFT) calculations. To ensure the reliability of such calculations, we design novel basis sets for relativistic DFT, via the particle swarm optimization algorithm to optimize the coefficients of the new basis sets and the polarization-consistent basis set idea's extension to heavy elements to eliminate the basis set error from DFT calculations. The resulting basis sets enable the well-grounded evaluation of relativistic DFT against "gold-standard" CCSD(T) results. Accounting for strong relativistic effects, including spin-orbit interaction, via our redesigned infrastructure, we elucidate a noticeable dissimilarity between At- and I- in halide-water force field parameters, radial distribution functions, diffusion coefficients, and hydration energies. This work establishes the framework for the systematic development of polarization-consistent basis sets for relativistic DFT and accurate force fields for molecular dynamics simulations to be used in large-scale models of complex molecular systems with elements from the bottom of the periodic table, including actinides and even superheavy elements.

Relativistic resolution-of-the-identity with Cholesky integral decomposition

In this study, we present an efficient integral decomposition approach called the restricted-kinetic-balance resolution-of-the-identity (RKB-RI) algorithm, which utilizes a tunable RI method based on the Cholesky integral decomposition for in-core relativistic quantum chemistry calculations. The RKB-RI algorithm incorporates the restricted-kinetic-balance condition and offers a versatile framework for accurate computations. Notably, the Cholesky integral decomposition is employed not only to approximate symmetric large-component electron repulsion integrals but also those involving small-component basis functions. In addition to comprehensive error analysis, we investigate crucial conditions, such as the kinetic balance condition and variational stability, which underlie the applicability of Dirac relativistic electronic structure theory. We compare the computational cost of the RKB-RI approach with the full in-core method to assess its efficiency. To evaluate the accuracy and reliability of the RKB-RI method proposed in this work, we employ actinyl oxides as benchmark systems, leveraging their properties for validation purposes. This investigation provides valuable insights into the capabilities and performance of the RKB-RI algorithm and establishes its potential as a powerful tool in the field of relativistic quantum chemistry.

Exact two-component Hamiltonians for relativistic quantum chemistry: Two-electron picture-change corrections made simple

Based on self-consistent field (SCF) atomic mean-field (amf) quantities, we present two simple yet computationally efficient and numerically accurate matrix-algebraic approaches to correct both scalar-relativistic and spin-orbit two-electron picture-change effects (PCEs) arising within an exact two-component (X2C) Hamiltonian framework. Both approaches, dubbed amfX2C and e(xtended)amfX2C, allow us to uniquely tailor PCE corrections to mean-field models, viz. Hartree-Fock or Kohn-Sham DFT, in the latter case also avoiding the need for a point-wise calculation of exchange-correlation PCE corrections. We assess the numerical performance of these PCE correction models on spinor energies of group 18 (closed-shell) and group 16 (open-shell) diatomic molecules, achieving a consistent ≈10^-5 Hartree accuracy compared to reference four-component data. Additional tests include SCF calculations of molecular properties such as absolute contact density and contact density shifts in copernicium fluoride compounds (CnF_n, n = 2,4,6), as well as equation-of-motion coupled-cluster calculations of x-ray core-ionization energies of 5d- and 6d-containing molecules, where we observe an excellent agreement with reference data. To conclude, we are confident that our (e)amfX2C PCE correction models constitute a fundamental milestone toward a universal and reliable relativistic two-component quantum-chemical approach, maintaining the accuracy of the parent four-component one at a fraction of its computational cost.

From the gas phase to the solid state: The chemical bonding in the superheavy element flerovium

As early as 1975, Pitzer suggested that copernicium, flerovium and oganesson are volatile substances behaving noble-gas like because of their closed-shell configurations and accompanying relativistic effects. It is, however, precarious to predict the chemical bonding and physical behavior of a solid by knowledge of the atomic or molecular properties only. Copernicium and oganesson have been analyzed very recently by our group. Both are predicted to be semi-conductors and volatile substances with rather low melting and boiling points, which may justify a comparison with the noble gas elements. Here we study closed-shell flerovium in detail to predict solid-state properties including the melting point from a decomposition of the total energy into many-body forces derived from relativistic coupled-cluster and from density functional theory. The convergence of such a decomposition for flerovium is critically analyzed, and the problem of using density functional theory is highlighted. We predict that flerovium is in many ways not behaving like a typical noble gas element despite its closed-shell 7$p_{1/2}^2$ configuration and resulting weak interactions. Unlike for the noble gases, the many-body expansion in terms of the interaction energy is not converging smoothly. This makes the accurate prediction of phase transitions very difficult. Nevertheless, a first prediction by Monte-Carlo simulation estimates the melting point at $284\pm 50$ K. Furthermore, calculations for the electronic band gap suggests that flerovium is a semi-conductor similar to copernicium

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree-Fock virtual orbitals by natural orbitals (NOs). In this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller-Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties-electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F-Ts), and parity-violating energy differences for H₂Z₂ (Z = O-Se)-we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

Recent advances in computational 31 P NMR: Part 1. Chemical shifts

This is the first part of two closely related reviews dealing with the computation of phosphorus-31 nuclear magnetic resonance chemical shifts in a wide series of organophosphorus compounds including complexes, clusters, and bioorganic phosphorus compounds. In particular, the analysis of the accuracy factors, such as substitution effects, solvent effects, vibrational corrections, and relativistic effects, is presented. This review is dedicated to the Full Member of the Russian Academy of Sciences Professor Boris A. Trofimov in view of his invaluable contribution to the field of synthesis, nuclear magnetic resonance, and computation studies of organophosphorus compounds.

The DIRAC code for relativistic molecular calculations

DIRAC is a freely distributed general-purpose program system for one-, two-, and four-component relativistic molecular calculations at the level of Hartree-Fock, Kohn-Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, electron propagator, and various flavors of coupled cluster theory. At the self-consistent-field level, a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows for the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding model, and the frozen density embedding model.

Ionization potentials of the superheavy element livermorium (Z = 116)

Ionization potentials (IPs) of the superheavy element (SHE) livermorium (Lv) and its ions Lvⁿ⁺ (n = 1, …, 6) are obtained using the multiconfiguration Dirac-Hartree-Fock method. The effects of electron correlation in the subshells {6s, 6p, 6d, 7s, 7p} are taken into account, together with the Breit interaction and quantum electrodynamic (QED) effects. In Lv, the strong relativistic effect causes a large splitting between the energies of the 7p_1/2 and 7p_3/2 orbitals, which results in a large difference between IP₃ and IP₂. As a consequence, the behavior of the IPs of Lv differs from that of the lighter oxygen group elements: among the IP_n ^Z (Z = Se, Te, Po, Lv; n = 1, …, 6), IP_1,2 ^Lv are the smallest, whereas IP_3,4,5,6 ^Lv are the second largest among the IP_3,4,5,6 ^Z. This jump in IP can be taken to be a natural characteristic of SHEs because the calculations of the difference between IP₃ and IP₂ are only weakly affected by electron correlation, the Breit interaction, and QED effects, with only the relativistic effect being significant. We also show that the energies and IPs of neutral Lv and Lv⁺ are clearly influenced by the electron correlation effect in the subshells {6s, 6p, 6d}. The Breit interaction and QED have an effect on the energies that has an exponential dependence on the atomic number, although they have only a weak influence on the IPs. The analysis of the stabilities of the 2+, 4+, and 6 + states of Lv shows good agreement with predictions from other studies.

Relativistic effects on the aromaticity of E3M3H3 (E = C–Pb; M = N–Bi) benzene analogues

The relativistic effects on the aromaticity of a set of benzene analogues, E₃M₃H₃ (E = C–Pb; M = N–Bi) heterocycles, using magnetically induced current density (MICD) and the NICS_zz component of the conventional nucleus independent chemical shift (NICS), is hereby examined.

Mössbauer isomer shifts and effective contact densities obtained by the exact two-component (X2C) relativistic method and its local variants

A standalone program to calculate scalar relativistic effective contact densities.

Computational 1 H NMR: Part 1. Theoretical background

This is the first one of the three closely interrelated reviews to be published in Magnetic Resonance in Chemistry dealing with accordingly theoretical background, chemical applications, and biochemical studies of and by means of computational ¹ H NMR. Presented in the first part of the review is a general outline of the modern theoretical methods and accuracy factors of computational ¹ H NMR involving locally dense basis set schemes, solvent effects, vibrational corrections, and relativistic effects performed at the density functional theory and/or nonempirical levels. This review is dedicated to Prof. Stephan Sauer in view of his invaluable contribution to the field of computational nuclear magnetic resonance.

Resolution-of-identity accelerated relativistic two- and four-component electron dynamics approach to chiroptical spectroscopies

We present an implementation and application of electron dynamics based on real-time time-dependent density functional theory (RT-TDDFT) and relativistic 2-component X2C and 4-component Dirac-Coulomb (4c) Hamiltonians to the calculation of electron circular dichroism and optical rotatory dispersion spectra. In addition, the resolution-of-identity approximation for the Coulomb term (RI-J) is introduced into RT-TDDFT and formulated entirely in terms of complex quaternion algebra. The proposed methodology was assessed on the dimethylchalcogenirane series, C₄H₈X (X = O, S, Se, Te, Po, Lv), and the spectra obtained by non-relativistic and relativistic methods start to disagree for Se and Te, while dramatic differences are observed for Po and Lv. The X2C approach, even in its simplest one-particle form, reproduces the reference 4c results surprisingly well across the entire series while offering an 8-fold speed-up of the simulations. An overall acceleration of RT-TDDFT by means of X2C and RI-J increases with system size and approaches a factor of almost 25 when compared to the full 4c treatment, without compromising the accuracy of the final spectra. These results suggest that one-particle X2C electron dynamics with RI-J acceleration is an attractive method for the calculation of chiroptical spectra in the valence region.

Normal halogen dependence of 13 C NMR chemical shifts of halogenomethanes revisited at the four-component relativistic level

The 'Normal Halogen Dependence' of ¹³ C NMR chemical shifts in the series of halogenomethanes is revisited at the four-component relativistic level. Calculations of ¹³ C NMR chemical shifts of 70 halogenomethanes have been carried out at the density functional theory (DFT) and MP2 levels with taking into account relativistic effects using the four-component relativistic theory of Dirac-Coulomb within the different computational methods (4RPA, 4OPW91) and hybrid computational schemes (MP2 + 4RPA, MP2 + 4OPW91). The most efficient computational protocols are derived for practical purposes. Relativistic shielding effect reaches as much as several hundreds of ppm for heavy halogenomethanes, and to account for this effect in comparison with experiment at the qualitative level, relativistic Dyall's basis sets of triple-zeta quality or higher are to be used within the framework of the four-component relativistic theory taking into account solvent effects. Relativistic geometrical optimization (as compared with the non-relativistic level) is essential for the molecules containing at least two iodines at one carbon atom. Copyright © 2016 John Wiley & Sons, Ltd.

MP2 calculation of (77) Se NMR chemical shifts taking into account relativistic corrections

The main factors affecting the accuracy and computational cost of the Second-order Möller-Plesset perturbation theory (MP2) calculation of (77) Se NMR chemical shifts (methods and basis sets, relativistic corrections, and solvent effects) are addressed with a special emphasis on relativistic effects. For the latter, paramagnetic contribution (390-466 ppm) dominates over diamagnetic term (192-198 ppm) resulting in a total shielding relativistic correction of about 230-260 ppm (some 15% of the total values of selenium absolute shielding constants). Diamagnetic term is practically constant, while paramagnetic contribution spans over 70-80 ppm. In the (77) Se NMR chemical shifts scale, relativistic corrections are about 20-30 ppm (some 5% of the total values of selenium chemical shifts). Solvent effects evaluated within the polarizable continuum solvation model are of the same order of magnitude as relativistic corrections (about 5%). For the practical calculations of (77) Se NMR chemical shifts of the medium-sized organoselenium compounds, the most efficient computational protocols employing relativistic Dyall's basis sets and taking into account relativistic and solvent corrections are suggested. The best result is characterized by a mean absolute error of 17 ppm for the span of (77) Se NMR chemical shifts reaching 2500 ppm resulting in a mean absolute percentage error of 0.7%.

A theoretical benchmark study of the spectroscopic constants of the very heavy rare gas dimers

The binding energy of the superheavy dimer Uuo₂ is considerably larger than that of its lighter homologues, despite a 40% reduction due to spin-other orbit interaction.