Mixed-sector intermediate Hamiltonian Fock-space coupled cluster approach

Triple Electron Attachments with a New Intermediate-Hamiltonian Fock-Space Coupled-Cluster Method

The implementation of a new intermediate-Hamiltonian Fock-space coupled-cluster (IHFSCC) scheme for the (3,0) sector of the Fock space is reported. In this IHFSCC approach, the three-body contributions in the cluster operator S(3,0) corresponding to the (3,0) sector of the Fock space are considered, while S(1,0) and S(2,0) at the (1,0) and (2,0) level only include one- and two-body operators. By introducing a suitable partition of the wave operator, the intermediate Hamiltonian, which only depends on the two-body operator of S(1,0), is obtained. S(2,0) and S(3,0) are not required within this new IHFSCC scheme, and a large reference space can be possibly employed. The performance of this (3,0) IHFSCC method in calculating triple ionization potentials and excitation energies for atoms and cations with a ground p3 configuration as well as adiabatic excitation energies for some molecules is investigated. The effect of the number of active virtual orbitals and three different types of orbitals, i.e., reference orbitals, restricted open-shell Hartree-Fock orbitals (ROHF) of the target systems, and canonicalized ROHF orbitals, on IHFSCC results, is also studied. Our calculations indicate that reasonable results can be achieved with this (3,0) IHFSCC method when a minimal reference space is employed. Further increasing the number of active orbitals does not necessarily improve the results. In addition, the IHFSCC results using canonicalized ROHF orbitals generally agree well with reference values, and they are not very sensitive to the number of active orbitals compared with results using the reference orbitals. The new (3,0) IHFSCC method can be applied to open-shell systems with three unpaired electrons with reasonable accuracy at a relatively low computational cost.

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.

Tailored coupled cluster theory in varying correlation regimes

The tailored coupled cluster (TCC) approach is a promising ansatz that preserves the simplicity of single-reference coupled cluster theory while incorporating a multi-reference wave function through amplitudes obtained from a preceding multi-configurational calculation. Here, we present a detailed analysis of the TCC wave function based on model systems, which require an accurate description of both static and dynamic correlation. We investigate the reliability of the TCC approach with respect to the exact wave function. In addition to the error in the electronic energy and standard coupled cluster diagnostics, we exploit the overlap of TCC and full configuration interaction wave functions as a quality measure. We critically review issues, such as the required size of the active space, size-consistency, symmetry breaking in the wave function, and the dependence of TCC on the reference wave function. We observe that possible errors caused by symmetry breaking can be mitigated by employing the determinant with the largest weight in the active space as reference for the TCC calculation. We find the TCC model to be promising in calculations with active orbital spaces which include all orbitals with a large single-orbital entropy, even if the active spaces become very large and then may require modern active-space approaches that are not restricted to comparatively small numbers of orbitals. Furthermore, utilizing large active spaces can improve on the TCC wave function approximation and reduce the size-consistency error because the presence of highly excited determinants affects the accuracy of the coefficients of low-excited determinants in the active space.

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.

QED effects on individual atomic orbital energies

Several issues, concerning QED corrections, that are important in precise atomic calculations are presented. The leading QED corrections, self-energy and vacuum polarization, to the orbital energy for selected atoms with 30 ≤ Z ≤ 118 have been calculated. The sum of QED and Breit contributions to the orbital energy is analyzed. It has been found that for ns subshells the Breit and QED contributions are of comparative size, but for np and nd subshells the Breit contribution takes a major part of the QED+Breit sum. It has also, been found that the Breit to leading QED contributions ratio for ns subshells is almost independent of Z. The Z-dependence of QED and Breit+QED contributions per subshell is shown. The fitting coefficients may be used to estimate QED effects on inner molecular orbitals. We present results of our calculations for QED contributions to orbital energy of valence ns-subshell for group 1 and 11 atoms and discuss about the reliability of these numbers by comparing them with experimental first ionization potential data.

Spin–orbit effects in optical spectra of gold–silver trimers

Photodissociation spectra of cationic gold–silver trimers are analysed using relativistic electronic structure theories paying special attention to the importance of spin–orbit and charge transfer effects.

Relativistic Coupled Cluster Calculations with Variational Quantum Electrodynamics Resolve the Discrepancy between Experiment and Theory Concerning the Electron Affinity and Ionization Potential of Gold

The first ionization potential (IP) and electron affinity (EA) of the gold atom have been determined to an unprecedented accuracy using relativistic coupled cluster calculations up to the pentuple excitation level including the Breit and QED contributions. We reach meV accuracy (with respect to the experimental values) by carefully accounting for all individual contributions beyond the standard relativistic coupled cluster approach. Thus, we are able to resolve the long-standing discrepancy between experimental and theoretical IP and EA of gold.

Effect of Triples to Dipole Moments in Fock-Space Multireference Coupled Cluster Method

In this paper, we present the new implementation of partial triples for the dipole moment of doublet radicals in Lagrangian formulation of Fock-space multireference coupled cluster (Λ-FSMRCC) response method. We have implemented a specific scheme of noniterative triples, in addition to singles and doubles schemes, which accounts for the effects appearing at least at the third order in dipole moments. The method is applied to the ground states of OH, OOH, HCOO, CN, CH, and PO radicals.

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).

On the performance of the intermediate Hamiltonian Fock-space coupled-cluster method on linear triatomic molecules: The electronic spectra of NpO2+, NpO22+, and PuO22+

In this paper we explore the use of the novel relativistic intermediate Hamiltonian Fock-space coupled-cluster method in the calculation of the electronic spectrum for small actinyl ions (NpO2+, NpO2(2+), and PuO2(2+)). It is established that the method, in combination with uncontracted double-zeta quality basis sets, yields excitation energies in good agreement with experimental values, and better than those obtained previously with other theoretical methods. We propose the reassignment of some of the peaks that were observed experimentally, and confirm other assignments.

Extrapolated intermediate Hamiltonian coupled-cluster approach: Theory and pilot application to electron affinities of alkali atoms

The intermediate Hamiltonian (IH) coupled-cluster method makes possible the use of very large model spaces in coupled-cluster calculations without running into intruder states. This is achieved at the cost of approximating some of the IH matrix elements, which are not taken at their rigorous effective Hamiltonian (EH) value. The extrapolated intermediate Hamiltonian (XIH) approach proposed here uses a parametrized IH and extrapolates it to the full EH, with model spaces larger by several orders of magnitude than those possible in EH coupled-cluster methods. The flexibility and resistance to intruders of the IH approach are thus combined with the accuracy of full EH. Various extrapolation schemes are described. A pilot application to the electron affinities (EAs) of alkali atoms is presented, where converged EH results are obtained by XIH for model spaces of approximately 20,000 determinants; direct EH calculations converge only for a one-dimensional model space. Including quantum electrodynamic effects, the average XIH error for the EAs is 0.6 meV and the largest error is 1.6 meV. A new reference estimate for the EA of Fr is proposed at 486+/-2 meV.

A Case Study of State-Specific and State-Averaged Brueckner Equation-of-Motion Coupled-Cluster Theory: The Ionic-Covalent Avoided Crossing in Lithium Fluoride

State-specific Brueckner equation-of-motion coupled-cluster theory (SS-B-EOMCC) is summarized, which can be considered an internally contracted version of a state-selective multireference coupled-cluster theory, which, however, is not entirely size-consistent. The method is applicable to general multireference problems, adheres to the space and spin symmetries of the molecular system, is straightforwardly extended to a state-averaged version, and has an associated perturbative variant which yields results close to the full coupled-cluster treatment. A key strength is that Brueckner orbitals are used, such that orbitals are optimized in the presence of dynamic correlation. A number of variations on the theme of SS-EOMCC is applied to study the ionic-covalent avoided crossing in LiF in a 6-311++G(3df,3pd) basis set. While reasonable results are obtained at the state-averaged level, the iterative solution process does not consistently converge for SS-EOMCC, due to the non-Hermiticity of the transformed Hamiltonian which may yield complex eigenvalues upon truncated diagonalization. This leads to an irrevocable breakdown of the state-specific EOMCC approach. We indicate some future directions that can resolve some of the problems with the SS-EOMCC methodology, as revealed by the demanding test case of the LiF potential energy curves.