Coupled-cluster method for open-shell heavy-element systems with spin-orbit coupling

Analytic evaluation of energy first derivatives for spin-orbit coupled-cluster singles and doubles augmented with noniterative triples method: General formulation and an implementation for first-order properties

A formulation of analytic energy first derivatives for the coupled-cluster singles and doubles augmented with noniterative triples [CCSD(T)] method with spin-orbit coupling included at the orbital level and an implementation for evaluation of first-order properties are reported. The standard density-matrix formulation for analytic CC gradient theory adapted to complex algebra has been used. The orbital-relaxation contributions from frozen core, occupied, virtual, and frozen virtual orbitals to analytic spin-orbit CCSD(T) gradients are fully taken into account and treated efficiently, which is of importance to calculations of heavy elements. Benchmark calculations of first-order properties including dipole moments and electric-field gradients using the corresponding exact two-component property integrals are presented for heavy-element containing molecules to demonstrate the applicability and usefulness of the present analytic scheme.

Effect of spin-orbit coupling on strong field ionization simulated with time-dependent configuration interaction

Time-dependent configuration interaction with a complex absorbing potential has been used to simulate strong field ionization by intense laser fields. Because spin-orbit coupling changes the energies of the ground and excited states, it can affect the strong field ionization rate for molecules containing heavy atoms. Configuration interaction with single excitations (CIS) has been employed for strong field ionization of closed shell systems. Single and double excitation configuration interaction with ionization (CISD-IP) has been used to treat ionization of degenerate states of cations on an equal footing. The CISD-IP wavefunction consists of ionizing single (one hole) and double (two hole/one particle) excitations from the neutral atom. Spin-orbit coupling has been implemented using an effective one electron spin-orbit coupling operator. The effective nuclear charge in the spin-orbit coupling operator has been optimized for Ar⁺, Kr⁺, Xe⁺, HX⁺ (X = Cl, Br, and I). Spin-orbit effects on angular dependence of the strong field ionization have been studied for HX and HX⁺. The effects of spin-orbit coupling are largest for ionization from the π orbitals of HX⁺. In a static field, oscillations are seen between the ²Π_3/2 and ²Π_1/2 states of HX⁺. For ionization of HX⁺ by a two cycle circularly polarized pulse, a single peak is seen when the maximum in the carrier envelope is perpendicular to the molecular axis and two peaks are seen when it is parallel to the axis. This is the result of the greater ionization rate for the π orbitals than for the σ orbitals.

Relativistic double-ionization equation-of-motion coupled-cluster method: Application to low-lying doubly ionized states

This article deals with the extension of the relativistic double-ionization equation-of-motion coupled-cluster (DI-EOMCC) method [H. Pathak et al. Phys. Rev. A 90, 010501(R) (2014)] for the molecular systems. The Dirac-Coulomb Hamiltonian with four-component spinors is considered to take care of the relativistic effects. The implemented method is employed to compute a few low-lying doubly ionized states of noble gas atoms (Ar, Kr, Xe, and Rn) and Cl₂, Br₂, HBr, and HI. Additionally, we presented results with two intermediate schemes in the four-component relativistic DI-EOMCC framework to understand the role of electron correlation. The computed double ionization spectra for the atomic systems are compared with the values from the non-relativistic DI-EOMCC method with spin-orbit coupling [Z. Wang et al. J. Chem. Phys. 142, 144109 (2015)] and the values from the National Institute of Science and Technology (NIST) database. Our atomic results are found to be in good agreement with the NIST values. Furthermore, the obtained results for the molecular systems agree well with the available experimental values.

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.

Density functional theory investigation on iridium(iii) complexes for efficient blue electrophosphorescence

The geometrical structures, electronic structures, optoelectronic properties and phosphorescence efficiencies of blue-emitting phosphors [Ir(fpmi)2(pyim)], [Ir(pyim)2(fpmi)], [Ir(fpmi)2(fptz)], [Ir(fpmi)2(pypz)] and [Ir(tfmppz)2(pyim)]), were investigated by DFT and TDDFT methods.

Ground-State Gas-Phase Structures of Inorganic Molecules Predicted by Density Functional Theory Methods

We tested a battery of density functional theory (DFT) methods ranging from generalized gradient approximation (GGA) via meta-GGA to hybrid meta-GGA schemes as well as Møller-Plesset perturbation theory of the second order and a single and double excitation coupled-cluster (CCSD) theory for their ability to reproduce accurate gas-phase structures of di- and triatomic molecules derived from microwave spectroscopy. We obtained the most accurate molecular structures using the hybrid and hybrid meta-GGA approximations with B3PW91, APF, TPSSh, mPW1PW91, PBE0, mPW1PBE, B972, and B98 functionals, resulting in lowest errors. We recommend using these methods to predict accurate three-dimensional structures of inorganic molecules when intramolecular dispersion interactions play an insignificant role. The structures that the CCSD method predicts are of similar quality although at considerably larger computational cost. The structures that GGA and meta-GGA schemes predict are less accurate with the largest absolute errors detected with BLYP and M11-L, suggesting that these methods should not be used if accurate three-dimensional molecular structures are required. Because of numerical problems related to the integration of the exchange-correlation part of the functional and large scattering of errors, most of the Minnesota models tested, particularly MN12-L, M11, M06-L, SOGGA11, and VSXC, are also not recommended for geometry optimization. When maintaining a low computational budget is essential, the nonseparable gradient functional N12 might work within an acceptable range of error. As expected, the DFT-D3 dispersion correction had a negligible effect on the internuclear distances when combined with the functionals tested on nonweakly bonded di- and triatomic inorganic molecules. By contrast, the dispersion correction for the APF-D functional has been found to shorten the bonds significantly, up to 0.064 Å (AgI), in Ag halides, BaO, BaS, BaF, BaCl, Cu halides, and Li and Na halides and hydrides. These results do not agree well with very accurate structures derived from microwave spectroscopy; we therefore believe that the dispersion correction in the APF-D method should be reconsidered. Finally, we found that inaccurate structures can easily lead to errors of few kcal/mol in single-point energies.