ReSpect: Relativistic spectroscopy DFT program package

Perturbation Theory Treatment of Spin-Orbit Coupling II: A Coupled Perturbed Kohn-Sham Method

A noncanonical coupled perturbed Kohn-Sham density functional theory (KS-DFT)/Hartree-Fock (HF) treatment of spin-orbit coupling (SOC) is provided. We take the scalar-relativistic KS-DFT/HF solution, obtained with a relativistic effective core potential, as the zeroth-order approximation. Explicit expressions are given for the total energy through the 4th order, which satisfy the 2n + 1 rule. Second-order expressions are provided for orbital energies and density variables of spin-current DFT. Test calculations are carried out on the halogen homonuclear diatomic and hydride molecules, including 6p and 7p elements, as well as open-shell negative ions. The computed properties through second or third order match well with those from reference two-component self-consistent field calculations for total and orbital energies as well as spin-current densities. In only one case (At₂^-) did a significant deviation occur for the remaining density variables. Our coupled perturbation theory approach provides an efficient way of adding the effect of SOC to a scalar-relativistic single-reference KS-DFT/HF treatment, in particular because it does not require diagonalization in the two-component spinor basis, leading to saving factors on the number of required floating-point operations that may exceed one order of magnitude.

Crystal and Substituent Effects on Paramagnetic NMR Shifts in Transition-Metal Complexes

Nuclear magnetic resonance (NMR) spectroscopy of paramagnetic molecules provides detailed information about their molecular and electron-spin structure. The paramagnetic NMR spectrum is a very rich source of information about the hyperfine interaction between the atomic nuclei and the unpaired electron density. The Fermi-contact contribution to ligand hyperfine NMR shifts is particularly informative about the nature of the metal–ligand bonding and the structural arrangements of the ligands coordinated to the metal center. In this account, we provide a detailed experimental and theoretical NMR study of compounds of Cr(III) and Cu(II) coordinated with substituted acetylacetonate (acac) ligands in the solid state. For the first time, we report the experimental observation of extremely paramagnetically deshielded ¹³C NMR resonances for these compounds in the range of 900–1200 ppm. We demonstrate an excellent agreement between the experimental NMR shifts and those calculated using relativistic density-functional theory. Crystal packing is shown to significantly influence the NMR shifts in the solid state, as demonstrated by theoretical calculations of various supramolecular clusters. The resonances are assigned to individual atoms in octahedral Cr(acac)₃ and square-planar Cu(acac)₂ compounds and interpreted by different electron configurations and magnetizations at the central metal atoms resulting in different spin delocalizations and polarizations of the ligand atoms. Further, effects of substituents on the ¹³C NMR resonance of the ipso carbon atom reaching almost 700 ppm for Cr(acac)₃ compounds are interpreted based on the analysis of Fermi-contact hyperfine contributions.

The ligand NMR shifts in paramagnetic acetylacetonato Cr(III) and Cu(II) complexes have been predicted and measured in the solid state and interpreted by relativistic DFT calculations. The effects of the metal atom, ligand, and crystal packing on the spin delocalization and polarization reflected in the Fermi-contact contribution to the hyperfine interaction are rationalized.

NMR Indirect Spin-Spin Coupling Constants in a Modern Quasi-Relativistic Density Functional Framework

A quasi-relativistic implementation of NMR indirect spin-spin coupling constants is presented. The exact two-component (X2C) Hamiltonian and its diagonal local approximation to the unitary decoupling transformation (DLU) are utilized together with the (modified) screened nuclear spin-orbit approach. In a restricted kinetic balance, the finite nucleus model is available for both the scalar and vector potentials. The implementation supports density functionals up to the fourth rung of Jacob's ladder, i.e., (range-separated) hybrid and local hybrid functionals based on a seminumerical ansatz. We assess the quality of our quasi-relativistic X2C approach by comparison with "fully" relativistic four-component results for small main-group molecules and alkynyl compounds. The mean absolute error introduced by the DLU scheme is less than 0.05 × 10¹⁹ T J^-2 of the reduced coupling constant for the small main-group molecules and 0.5 Hz for the alkynyl compounds. Thus, the error is significantly smaller than finite nucleus size effects for heavy elements. The basis set convergence and the impact of different density functional approximations are further studied. We propose a simple scheme to develop segmented-contracted relativistic all-electron basis sets for NMR spin-spin couplings. Our implementation allows us to perform calculations of extended molecules with reasonable computational effort, which is illustrated for the ¹J(¹¹⁹Sn, ³¹P) coupling constant of a low-valent tin phosphinidenide complex. The corresponding results are in good agreement with the experimental findings.

Electrochemical Potential of the Metal Organic Framework MIL-101(Fe) as Cathode Material in Li-Ion Batteries

We discuss the characteristic factors that determine the electrochemical potentials in a metal-organic framework used as cathode for Li-ion batteries via density functional theory-based simulations. Our focus is on MIL-101(Fe) cathode material. Our study gives insight into the role of local atomic environment and structural deformations in generating electrochemical potential.

Relativistic EOM-CCSD for Core-Excited and Core-Ionized State Energies Based on the Four-Component Dirac-Coulomb(-Gaunt) Hamiltonian

We report an implementation of the core-valence separation approach to the four-component relativistic Hamiltonian-based equation-of-motion coupled-cluster with singles and doubles theory (CVS-EOM-CCSD) for the calculation of relativistic core-ionization potentials and core-excitation energies. With this implementation, which is capable of exploiting double group symmetry, we investigate the effects of the different CVS-EOM-CCSD variants and the use of different Hamiltonians based on the exact two-component (X2C) framework on the energies of different core-ionized and -excited states in halogen- (CH₃I, HX, and X^-, X = Cl-At) and xenon-containing (Xe, XeF₂) species. Our results show that the X2C molecular mean-field approach [Sikkema, J.; J. Chem. Phys. 2009, 131, 124116], based on four-component Dirac-Coulomb mean-field calculations (²DC^M), is capable of providing core excitations and ionization energies that are nearly indistinguishable from the reference four-component energies for up to and including fifth-row elements. We observe that two-electron integrals over the small-component basis sets lead to non-negligible contributions to core binding energies for the K and L edges for atoms such as iodine or astatine and that the approach based on Dirac-Coulomb-Gaunt mean-field calculations (²DCG^M) are significantly more accurate than X2C calculations for which screened two-electron spin-orbit interactions are included via atomic mean-field integrals.

Efficient Four-Component Dirac-Coulomb-Gaunt Hartree-Fock in the Pauli Spinor Representation

Four-component Dirac-Hartree-Fock is an accurate mean-field method for treating molecular systems where relativistic effects are important. However, the computational cost and complexity of the two-electron interaction make this method less common, even though we can consider the Dirac-Hartree-Fock Hamiltonian the "ground truth" of the electronic structure, barring explicit quantum electrodynamical effects. Being able to calculate these effects is then vital to the design of lower scaling methods for accurate predictions in computational spectroscopy and properties of heavy element complexes that must include relativistic effects for even qualitative accuracy. In this work, we present a Pauli quaternion formalism of maximal component and spin separation for computing the Dirac-Coulomb-Gaunt Hartree-Fock ground state, with a minimal floating point operation count algorithm. This approach also allows one to explicitly separate different spin physics from the two-body interactions, such as spin-free, spin-orbit, and spin-spin contributions. Additionally, we use this formalism to examine relativistic trends in the periodic table and analyze the basis set dependence of atomic gold and gold dimer systems.

Quantum chemical calculations of 77 Se and 125 Te nuclear magnetic resonance spectral parameters and their structural applications

An accurate quantum chemical (QC) modeling of ⁷⁷ Se and ¹²⁵ Te nuclear magnetic resonance (NMR) spectra is deeply involved in the NMR structural assignment for selenium and tellurium compounds that are of utmost importance both in organic and inorganic chemistry nowadays due to their huge application potential in many fields, like biology, medicine, and metallurgy. The main interest of this review is focused on the progress in QC computations of ⁷⁷ Se and ¹²⁵ Te NMR chemical shifts and indirect spin-spin coupling constants involving these nuclei. Different computational methodologies that have been used to simulate the NMR spectra of selenium and tellurium compounds since the middle of the 1990s are discussed with a strong emphasis on their accuracy. A special accent is placed on the calculations resorting to the relativistic methodologies, because taking into account the relativistic effects appreciably influences the precision of NMR calculations of selenium and, especially, tellurium compounds. Stereochemical applications of quantum chemical calculations of ⁷⁷ Se and ¹²⁵ Te NMR parameters are discussed so as to exemplify the importance of integrated approach of experimental and computational NMR techniques.

Generalization of Intrinsic Orbitals to Kramers-Paired Quaternion Spinors, Molecular Fragments, and Valence Virtual Spinors

Localization of molecular orbitals finds its importance in the representation of chemical bonding (and antibonding) and in the local correlation treatments beyond mean-field approximation. In this paper, we generalize the intrinsic atomic and bonding orbitals [G. Knizia, J. Chem. Theory Comput. 2013, 9, 11, 4834-4843] to relativistic applications using complex and quaternion spinors, as well as to molecular fragments instead of atomic fragments only. By performing a singular value decomposition, we show how localized valence virtual orbitals can be expressed on this intrinsic minimal basis. We demonstrate our method on systems of increasing complexity, starting from simple cases such as benzene, acrylic acid, and ferrocene molecules, and then demonstrate the use of molecular fragments and inclusion of relativistic effects for complexes containing heavy elements such as tellurium, iridium, and astatine. The aforementioned scheme is implemented into a standalone program interfaced with several different quantum chemistry packages.

Investigating the influence of relativistic effects on absorption spectra for platinum complexes with light-activated activity against cancer cells

We report the first systematic investigation of relativistic effects on the UV-vis spectra of two prototype complexes for so-called photo-activated chemotherapy (PACT), trans-trans-trans-[Pt(N3)2(OH)2(NH3)2] and cis-trans-cis-[Pt(N3)2(OH)2(NH3)2]. In PACT, design of new drugs requires in-depth understanding of the photo-activation mechanisms. A first step is usually to rationalize their UV-vis spectra for which time-dependent density functional theory (TD-DFT) is an indispensable tool. We carried out TD-DFT calculations with a systematic series of non-relativistic (NR), scalar-relativistic (SR), and four-component (4c) Hamiltonians. As expected, large differences are found between spectra calculated within 4c and NR frameworks, while the most intense features (found at higher energies below 300 nm) can be reasonably well reproduced within a SR framework. It is also shown that effective core potentials (ECPs) yield essentially similar results as all-electron SR calculations. Yet the underlying transitions can be strongly influenced by spin-orbit coupling, which is only present in the 4c framework: while this can affect both intense and less intense transitions in the spectra, the effect is most pronounced for weaker transitions at lower energies, above 300 nm. Since the investigated complexes are activated with light of wavelengths above 300 nm, employing a method with explicit inclusion of spin-orbit coupling may be crucial to rationalize the activation mechanism.

Behind the scenes of spin-forbidden decay pathways in transition metal complexes

The interpretation of the ultrafast photophysics of transition metal complexes following photo-absorption is quite involved as the heavy metal center leads to a complicated and entangled singlet-triplet manifold. This opens up multiple pathways for deactivation, often with competitive rates. As a result, intersystem crossing (ISC) and phosphorescence are commonly observed in transition metal complexes. A detailed understanding of such an excited-state structure and dynamics calls for state-of-the-art experimental and theoretical methodologies. In this review, we delve into the inability of non-relativistic quantum theory to describe spin-forbidden transitions, which can be overcome by taking into account spin-orbit coupling, whose importance grows with increasing atomic number. We present the quantum chemical theory of phosphorescence and ISC together with illustrative examples. Finally, a few applications are highlighted, bridging the gap between theoretical studies and experimental applications, such as photofunctional materials.

First-Principles Calculation of 1H NMR Chemical Shifts of Complex Metal Polyhydrides: The Essential Inclusion of Relativity and Dynamics

1H NMR spectroscopy has become an important technique for the characterization of transition-metal hydride complexes, whose metal-bound hydrides are often difficult to locate by X-ray diffraction. In this regard, the accurate prediction of 1H NMR chemical shifts provides a useful, but challenging, strategy to help in the interpretation of the experimental spectra. In this work, we establish a density-functional-theory protocol that includes relativistic, solvent, and dynamic effects at a high level of theory, allowing us to report an accurate and reliable interpretation of 1H NMR hydride chemical shifts of iridium polyhydride complexes. In particular, we have studied in detail the hydride chemical shifts of the [Ir6(IMe)8(CO)2H14]2+ complex in order to validate previous assignments. The computed 1H NMR chemical shifts are strongly dependent on the relativistic treatment, the choice of the DFT exchange-correlation functional, and the conformational dynamics. By combining a fully relativistic four-component electronic-structure treatment with ab initio molecular dynamics, we were able to reliably model both the terminal and bridging hydride chemical shifts and to show that two NMR hydride signals were inversely assigned in the experiment.

Carbonyl compounds of Rh, Ir, and Mt: electronic structure, bonding and volatility

First bond dissociation energies and other properties have been predicted for carbonyl compounds of group-9 elements including those of element 109, Mt, from relativistic DFT and CC calculations. A remarkable Λ-shape of the trends is observed, caused by strong relativistic effects on the valence AOs of Mt.

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.