The relativistic Fock-space coupled-cluster method for molecules: CdH and its ions

Relativistic Coupled Cluster with Completely Renormalized and Perturbative Triples Corrections

We have implemented noniterative triples corrections to the energy from coupled-cluster with single and double excitations (CCSD) within the 1-electron exact two-component (1eX2C) relativistic framework. The effectiveness of both the CCSD(T) and the completely renormalized (CR) CC(2,3) approaches are demonstrated by performing all-electron computations of the potential energy curves and spectroscopic constants of copper, silver, and gold dimers in their ground electronic states. Spin-orbit coupling effects captured via the 1eX2C framework are shown to be crucial for recovering the correct shape of the potential energy curves, and the correlation effects due to triples in these systems change the dissociation energies by about 0.1-0.2 eV or about 4-7%. We also demonstrate that relativistic effects and basis set size and contraction scheme are significantly more important in Au2 than in Ag2 or Cu2.

Relativistic Exact Two-Component Coupled-Cluster Study of Molecular Sensitivity Factors for Nuclear Schiff Moments

Relativistic exact two-component coupled-cluster calculations of molecular sensitivity factors for nuclear Schiff moments (NSMs) are reported. We focus on molecules containing heavy nuclei, especially octupole-deformed nuclei. Analytic relativistic coupled-cluster gradient techniques are used and serve as useful tools for identifying candidate molecules that sensitively probe for physics beyond the Standard Model in the hadronic sector. Notably, these tools enable straightforward "black-box" calculations. Two competing chemical mechanisms that contribute to the NSM are analyzed, illuminating the physics of ligand effects on NSM sensitivity factors.

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.

Intermediate Hamiltonian Fock-space coupled-cluster theory for excitation energies, double ionization potentials, and double electron attachments with spin–orbit coupling

The intermediate Hamiltonian Fock-space coupled-cluster methods at the singles and doubles level (IHFSCCSD) for excitation energies in the (1p, 1h) sector, double ionization potentials in the (0p, 2h) sector, and double electron attachments in the (2p, 0h) sector of the Fock space are implemented based on the CCSD method with spin–orbit coupling (SOC) included in the post-Hartree–Fock treatment using a closed-shell reference in this work. The active space is chosen to contain those orbitals that have the largest contribution to principal ionized or electron-attached states obtained from the equation-of-motion coupled-cluster calculations. Both time-reversal symmetry and spatial symmetry are exploited in the implementation. Our results show that the accuracy of IHFSCCSD results is closely related to the active space, and the sufficiency of the active space can be assessed from the percentage of transitions within the active space. In addition, unreasonable results may be encountered when the ionized or electron-attached states with a somewhat larger contribution from double excitations are included to determine the active space and cluster operators in the (0p, 1h) or (1p, 0h) sector of the Fock space. A larger active space may be required to describe SO splitting reliably than that in the scalar-relativistic calculations in some cases. The IHFSCCSD method with SOC developed in this work can provide reliable results for heavy-element systems when a sufficient active space built upon the principal ionization potential/electron affinity states is adopted.

Ab initio calculations in atoms, molecules, and solids, treating spin-orbit coupling and electron interaction on an equal footing

We incorporate explicit, non-perturbative treatment of spin-orbit coupling into ab initio auxiliary-field quantum Monte Carlo (AFQMC) calculations. The approach allows a general computational framework for molecular and bulk systems in which material specificity, electron correlation, and spin-orbit coupling effects can be captured accurately and on an equal footing, with favorable computational scaling vs system size. We adopt relativistic effective-core potentials that have been obtained by fitting to fully relativistic data and that have demonstrated a high degree of reliability and transferability in molecular systems. This results in a two-component spin-coupled Hamiltonian, which is then treated by generalizing the ab initio AFQMC approach. We demonstrate the method by computing the electron affinity in Pb, the bond dissociation energy in Br₂ and I₂, and solid Bi.

Equation-of-motion coupled-cluster theory for double electron attachment with spin-orbit coupling

We report implementation of the equation-of-motion coupled-cluster (EOM-CC) method for double electron-attachment (DEA) with spin-orbit coupling (SOC) at the CC singles and doubles (CCSD) level using a closed-shell reference in this work. The DEA operator employed in this work contains two-particle and three-particle one-hole excitations, and SOC is included in post-Hartree-Fock treatment. Time-reversal symmetry and spatial symmetry are exploited to reduce computational cost. The EOM-DEA-CCSD method with SOC allows us to investigate SOC effects of systems with two-unpaired electrons. According to our results on atoms, double ionization potentials (DIPs), excitation energies (EEs), and SO splittings of low-lying states are calculated reliably using the EOM-DEA-CCSD method with SOC. Its accuracy is usually higher than that of EOM-CCSD for EEs or DIPs if the same target can be reached from single excitations by choosing a proper closed-shell reference. However, performance of the EOM-DEA-CCSD method with SOC on molecules is not as good as that for atoms. Bond lengths for the ground and the several lowest excited states of GaH, InH, and TlH are underestimated pronouncedly, although reasonable EEs are obtained, and splittings of the ³Σ^- state from the π² configuration are calculated to be too small with EOM-DEA-CCSD.

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.

Splittings of d8 configurations of late-transition metals with EOM-DIP-CCSD and FSCCSD methods

Multireference methods are usually required for transition metal systems due to the partially filled d electrons. In this work, the single-reference equation-of-motion coupled-cluster method at the singles and doubles level for double ionization potentials (EOM-DIP-CCSD) is employed to calculate energies of states from the d⁸ configuration of late-transition metal atoms starting from a closed-shell reference. Its results are compared with those from the multireference Fock-space coupled-cluster method at the CCSD level (FSCCSD) for DIP from the same closed-shell reference. Both scalar-relativistic effects and spin-orbit coupling are considered in these calculations. Compared with all-electron FSCCSD results with four-component Dirac-Coulomb Hamiltonian, FSCCSD with relativistic effective core potentials can provide reasonable results, except for atoms with unstable reference. Excitation energies for states in the (n - 1)d⁸ns² configuration are overestimated pronouncedly with these two methods, and this overestimation is more severe than those in the (n - 1)d⁹ns¹ configuration. Error of EOM-CCSD on these excitation energies is generally larger than that of FSCCSD. On the other hand, relative energies of most of the states in the d⁸ configuration with respect to the lowest state in the same configuration are predicted reliably with EOM-DIP-CCSD, except for the ³P₀ state of Hg²⁺ and states in Ir⁺. FSCCSD can provide reasonable relative energies for the several lowest states, while its error tends to be larger for higher states.

Analytical energy gradients for ionized states using equation-of-motion coupled-cluster theory with spin-orbit coupling

Spin-orbit coupling (SOC) may have a significant effect on the structure and harmonic frequencies of particularly heavy p-block element compounds. However, reports on analytical energy gradients with SOC are scarce, especially for excited states. In this work, we implemented analytical energy gradients for ionized states using the equation-of-motion coupled-cluster (CC) theory at the CC singles and doubles level (EOM-IP-CCSD) with SOC. Effects of SOC on structure and harmonic frequencies as well as properties for both the ground and some excited states of open-shell compounds with one unpaired electron can be investigated efficiently with the present implementation. A closed-shell reference is required in the calculations, and SOC is included in post-Hartree-Fock treatment. Relativistic effective core potentials are employed in dealing with both scalar relativistic effects and SOC, and we treat perturbations that are even under time reversal in this work. Both time-reversal symmetry and double point group symmetry for D_2h ^* and its subgroups are exploited in the implementation. The method is applicable to states which can be reached by removing one electron from a closed-shell reference state. The results of some open-shell cations indicate the importance of SOC on structures and harmonic frequencies of heavy element compounds.

Combining the spin-separated exact two-component relativistic Hamiltonian with the equation-of-motion coupled-cluster method for the treatment of spin–orbit splittings of light and heavy elements

An accurate and efficient treatment of spin–orbit splittings has been achieved by combining the sf-X2C+soc-DKH1 Hamiltonian with the equation-of-motion coupled-cluster method.

Equation-of-Motion Coupled-Cluster Theory for Excitation Energies of Closed-Shell Systems with Spin-Orbit Coupling

Excitation energies of closed-shell systems based on the equation-of-motion (EOM) coupled-cluster theory at the singles and doubles (CCSD) level with spin-orbit coupling (SOC) included in the post-Hartree-Fock treatment are implemented in the present work. SOC can be included in both the CC and EOM steps (EOM-SOC-CCSD) or only in the EOM part (SOC-EOM-CCSD). The latter approach is an economical way to account for SOC effects, but excitation energies with this approach are not size-intensive. When the unlinked term in the latter approach is neglected (cSOC-EOM-CCSD), size-intensive excitation energies can be obtained. Time-reversal symmetry and spatial symmetry are exploited to reduce the computational effort. Imposing time-reversal symmetry results in a real matrix representation for the similarity-transformed Hamiltonian, which facilitates the requirement of time-reversal symmetry for new trial vectors in Davidson's algorithm. Results on some closed-shell atoms and molecules containing heavy elements show that EOM-SOC-CCSD can provide excitation energies and spin-orbit splittings with reasonable accuracy. On the other hand, the SOC-EOM-CCSD approach is able to afford accurate estimates of SOC effects for valence electrons of systems containing elements up to the fifth row, while cSOC-EOM-CCSD is less accurate for spin-orbit splittings of transitions involving p1/2 spinors, even for Kr.

Chemically Accurate Thermochemistry of Cadmium: An ab Initio Study of Cd + XY (X = H, O, Cl, Br; Y = Cl, Br)

Using a composite coupled cluster method employing sequences of correlation consistent basis sets for complete basis set (CBS) extrapolations and with explicit treatment of core-valence correlation and scalar and spin-orbit relativistic effects, the 0 K enthalpies of a wide range of cadmium-halide reactions, namely, Cd + (HCl, HBr, ClO, BrO, Cl2, BrCl, Br2) have been determined to an estimated accuracy of +/-1 kcal/mol. In addition, accurate equilibrium geometries, harmonic frequencies, and dissociation energies have been calculated at the same level of theory for all the diatomic (e.g., CdH, CdO, CdCl, CdBr) and triatomic (CdHCl, CdHBr, CdClO, CdBrO, CdCl2, CdBrCl, CdBr2) species involved in these reactions, some for the very first time. Like their mercury analogues, all of the abstraction reactions are predicted to be endothermic, while the insertion reactions are strongly exothermic with the formation of stable linear, Cd-centric complexes. With the exception of CdH and the reactions involving this species, the present results for the remaining Cd-containing systems are believed to be the most accurate to date.