Density functional theory for hyperfine coupling constants with the restricted-unrestricted approach

Mechanisms of Ligand Hyperfine Coupling in Transition-Metal Complexes: σ and π Transmission Pathways

Theoretical interpretation of hyperfine interactions was pioneered in the 1950s-1960s by the seminal works of McConnell, Karplus, and others for organic radicals and by Watson and Freeman for transition-metal (TM) complexes. In this work, we investigate a series of octahedral Ru(III) complexes with aromatic ligands to understand the mechanism of transmission of the spin density from the d-orbital of the metal to the s-orbitals of the ligand atoms. Spin densities and spin populations underlying ligand hyperfine couplings are analyzed in terms of π-conjugative or σ-hyperconjugative delocalization vs spin polarization based on symmetry considerations and restricted open-shell vs unrestricted wave function analysis. The transmission of spin density is shown to be most efficient in the case of symmetry-allowed π-conjugative delocalization, but when the π-conjugation is partially or fully symmetry-forbidden, it can be surpassed by σ-hyperconjugative delocalization. Despite a lower spin population of the ligand in σ-hyperconjugative transmission, the hyperfine couplings can be larger because of the direct involvement of the ligand s-orbitals in this delocalization pathway. We demonstrate a quantitative correlation between the hyperfine couplings of aromatic ligand atoms and the characteristics of the metal-ligand bond modulated by the trans substituent, a hyperfine trans effect.

Bonding in Nitrile Photo-dissociating Ruthenium Drug Candidates --A Local Vibrational Mode Study

In this work, we investigated bonding features 15 ruthenium complexes of the type [Ru(tpy)(L)-(CH3CN)]n+, containing the tridentate tpy ligand (tpy = 2,2':6',2'--terpyridine) and various bidentate ancillary ligands, 12 compounds originally synthesized by Loftus et al. (J. Phys. Chem. C 123, 10291-10299 (2019)) complemented with three additional complexes. The main focus of our work was to relate these local features to the experimental data of Loftus et al. which assess the efficiency of nitrile release in an indirect way via observed quantum yields for ruthenium water association after nitrile release. As a tool to quantitatively assess Ru-NC and Ru-L bonding we utilized the local vibrational mode analysis complemented by the topological analysis of the electron density and the natural bond orbital analysis. Interestingly, the stronger Ru-NC bonds have the greater observed quantum yields, leading to the conclusion that the observed quantum yields are a result of a complex interplay of several processes excluding a direct relationship between QY and Ru-NC or Ru-L bond strengths. We identified the ST splitting as one of the key players and not the Ru-NC bond strength, as one may have thought. In summary, this work has presented a modern computational tool set for the investigation of bonding features applied to nitrile photo-dissociating ruthenium drug candidates forming a valuable basis for future design and fine tuning of nitrile releasing ruthenium compounds, as well as for the understanding of how local properties affect overall experimental outcomes.

Dalton Project: A Python platform for molecular- and electronic-structure simulations of complex systems

The Dalton Project provides a uniform platform access to the underlying full-fledged quantum chemistry codes Dalton and LSDalton as well as the PyFraME package for automatized fragmentation and parameterization of complex molecular environments. The platform is written in Python and defines a means for library communication and interaction. Intermediate data such as integrals are exposed to the platform and made accessible to the user in the form of NumPy arrays, and the resulting data are extracted, analyzed, and visualized. Complex computational protocols that may, for instance, arise due to a need for environment fragmentation and configuration-space sampling of biochemical systems are readily assisted by the platform. The platform is designed to host additional software libraries and will serve as a hub for future modular software development efforts in the distributed Dalton community.

Selected-Nuclei Method for the Computation of Hyperfine Coupling Constants within Second-Order Møller-Plesset Perturbation Theory

We introduce a new ansatz to compute hyperfine coupling constants of selected nuclei at the level of second-order Møller-Plesset perturbation (MP2) and double-hybrid density functional theory with reduced computational effort, opening the route to the analyis of hyperfine coupling constants of large molecular structures. Our approach is based on a reformulation of the canonical MP2 term in atomic orbitals, thus exploiting the locality of electron correlation. We show that a perturbation-including integral screening reduces the scaling behavior of the number of significant two-electron integrals to sublinear. This selected-nuclei approach allows for an efficient computation within scaled-opposite spin (SOS) RI-MP2 on massively parallelized architectures such as graphical processor units (GPUs), thus enabling studies on the influence of the environment on hyperfine coupling constants.

Ro-vibrational averaging of the isotropic hyperfine coupling constant for the methyl radical

We present the first variational calculation of the isotropic hyperfine coupling constant of the carbon-13 atom in the CH3 radical for temperatures T = 0, 96, and 300 K. It is based on a newly calculated high level ab initio potential energy surface and hyperfine coupling constant surface of CH3 in the ground electronic state. The ro-vibrational energy levels, expectation values for the coupling constant, and its temperature dependence were calculated variationally by using the methods implemented in the computer program TROVE. Vibrational energies and vibrational and temperature effects for coupling constant are found to be in very good agreement with the available experimental data. We found, in agreement with previous studies, that the vibrational effects constitute about 44% of the constant's equilibrium value, originating mainly from the large amplitude out-of-plane bending motion and that the temperature effects play a minor role.

The Dalton quantum chemistry program system

Dalton is a powerful general-purpose program system for the study of molecular electronic structure at the Hartree-Fock, Kohn-Sham, multiconfigurational self-consistent-field, Møller-Plesset, configuration-interaction, and coupled-cluster levels of theory. Apart from the total energy, a wide variety of molecular properties may be calculated using these electronic-structure models. Molecular gradients and Hessians are available for geometry optimizations, molecular dynamics, and vibrational studies, whereas magnetic resonance and optical activity can be studied in a gauge-origin-invariant manner. Frequency-dependent molecular properties can be calculated using linear, quadratic, and cubic response theory. A large number of singlet and triplet perturbation operators are available for the study of one-, two-, and three-photon processes. Environmental effects may be included using various dielectric-medium and quantum-mechanics/molecular-mechanics models. Large molecules may be studied using linear-scaling and massively parallel algorithms. Dalton is distributed at no cost from http://www.daltonprogram.org for a number of UNIX platforms.

π-Stacking effects on the EPR parameters of a prototypical DNA spin label

The character and value of spin labels for probing environments like double-stranded DNA depend on the degree of change in the spin Hamiltonian parameters of the spin label induced by the environment. Herein we provide a systematic theoretical investigation of this issue, based on a density functional theory method applied to a spin labeled DNA model system, focusing on the dependence of the EPR properties of the spin label on the π stacking and hydrogen bonding that occur upon incorporating the spin label into the selected base pair inside DNA. It is found that the EPR spin Hamiltonian parameters of the spin label are only negligibly affected by its incorporation into DNA, when compared to its free form. This result gives a theoretical ground for the common empirical assumption regarding the behaviour of spin Hamiltonian parameters made in EPR based measurements of the distance between spin labels incorporated into DNA.

Role of zero-point vibrational corrections to carbon hyperfine coupling constants in organic π radicals

By analyzing a set of organic π radicals, we demonstrate that zero-point vibrational corrections give significant contributions to carbon hyperfine coupling constants, in one case even inducing a sign reversal for the coupling constant. We discuss the implications of these findings for the computational analysis of electron paramagnetic spectra based on hyperfine coupling constants evaluated at the equilibrium geometry of radicals. In particular, we note that a dynamical description that involves the nuclear motion is in many cases necessary in order to achieve a semi-quantitatively predictive theory for carbon hyperfine coupling constants. In addition, we discuss the implications of the strong dependence of the carbon hyperfine coupling constants on the zero-point vibrational corrections for the selection of exchange-correlation functionals in density functional theory studies of these constants.

Spin-adapted open-shell random phase approximation and time-dependent density functional theory. I. Theory

The spin-adaptation of single-reference quantum chemical methods for excited states of open-shell systems has been nontrivial. The primary reason is that the configuration space, generated by a truncated rank of excitations from only one component of a reference multiplet, is spin-incomplete. Those "missing" configurations are of higher ranks and can, in principle, be recaptured by a particular class of excitation operators. However, the resulting formalisms are then quite involved and there are situations [e.g., time-dependent density functional theory (TD-DFT) under the adiabatic approximation] that prevent one from doing so. To solve this issue, we propose here a tensor-coupling scheme that invokes all the components of a reference multiplet (i.e., a tensor reference) rather than increases the excitation ranks. A minimal spin-adapted n-tuply excited configuration space can readily be constructed by tensor products between the n-tuple tensor excitation operators and the chosen tensor reference. Further combined with the tensor equation-of-motion formalism, very compact expressions for excitation energies can be obtained. As a first application of this general idea, a spin-adapted open-shell random phase approximation is first developed. The so-called "translation rule" is then adopted to formulate a spin-adapted, restricted open-shell Kohn-Sham (ROKS)-based TD-DFT (ROKS-TD-DFT). Here, a particular symmetry structure has to be imposed on the exchange-correlation kernel. While the standard ROKS-TD-DFT can access only excited states due to singlet-coupled single excitations, i.e., only some of the singly excited states of the same spin (S(i)) as the reference, the new scheme can capture all the excited states of spin S(i)-1, S(i), or S(i)+1 due to both singlet- and triplet-coupled single excitations. The actual implementation and computation are very much like the (spin-contaminated) unrestricted Kohn-Sham-based TD-DFT. It is also shown that spin-contaminated spin-flip configuration interaction approaches can easily be spin-adapted via the tensor-coupling scheme.

Zero-point vibrational corrections to isotropic hyperfine coupling constants in polyatomic molecules

The present work addresses isotropic hyperfine coupling constants in polyatomic systems with a particular emphasis on a largely neglected, but a posteriori significant, effect, namely zero-point vibrational corrections. Using the density functional restricted-unrestricted approach, the zero-point vibrational corrections are evaluated for the allyl radical and four of its derivatives. In addition for establishing the numerical size of the zero-point vibrational corrections to the isotropic hyperfine coupling constants, we present simple guidelines useful for identifying hydrogens for which such corrections are significant. Based on our findings, we critically re-examine the computational procedures used for the determination of hyperfine coupling constants in general as well as the practice of using experimental hyperfine coupling constants as reference data when benchmarking and optimizing exchange-correlation functionals and basis sets for such calculations.

Density functional restricted-unrestricted approach for nonlinear properties: application to electron paramagnetic resonance parameters of square planar copper complexes

The density functional restricted-unrestricted approach for treatments of spin polarization effects in molecular properties using spin restricted Kohn-Sham theory has been extended from linear to nonlinear properties. It is shown that the spin polarization contribution to a nonlinear property has the form of a quadratic response function that includes the zero-order Kohn-Sham operator, in analogy to the lower order case where the spin polarization correction to an expectation value has the form of a linear response function. The developed approach is used to formulate new schemes for computation of electronic g-tensors and hyperfine coupling constants, which include spin polarization effects within the framework of spin restricted Kohn-Sham theory. The proposed computational schemes are in the present work employed to study the spin polarization effects on electron paramagnetic resonance spin Hamiltonian parameters of square planar copper complexes. The obtained results indicate that spin polarization gives rise to sizable contributions to the hyperfine coupling tensor of copper in all investigated complexes, while the electronic g-tensors of these complexes are only marginally affected by spin polarization and other factors, such as choice of exchange-correlation functional or molecular structures, will have more pronounced impact on the accuracy of the results.

Time-dependent density functional theory for nonlinear properties of open-shell systems

This paper presents response theory based on a spin-restricted Kohn-Sham formalism for computation of time-dependent and time-independent nonlinear properties of molecules with a high spin ground state. The developed approach is capable to handle arbitrary perturbations and constitutes an efficient procedure for evaluation of electric, magnetic, and mixed properties. Apart from presenting the derivation of the proposed approach, we show results from illustrating calculations of static and dynamic hyperpolarizabilities of small Si(3n+1)H(6n+3) (n=0,1,2) clusters which mimic Si(111) surfaces with dangling bond defects. The results indicate that the first hyperpolarizability tensor components of Si(3n+1)H(6n+3) have an ordering compatible with the measurements of second harmonic generation in SiO2/Si(111) interfaces and, therefore, support the hypothesis that silicon surface defects with dangling bonds are responsible for this phenomenon. The results exhibit a strong dependence on the quality of basis set and exchange-correlation functional, showing that an appropriate set of diffuse functions is required for reliable predictions of the first hyperpolarizability of open-shell compounds.

DFT Calculations of EPR Parameters in an Ionic Lattice of [M(CN)4]3-(M = Ni, Pd, Fe, Ru, Os) Complexes

The electronic g-tensor and hyperfine coupling constants were calculated for cyanide coordination complexes [M(CN)4]3- (M = Ni, Pd, Fe, Ru, Os) in KCl or NaCl host lattices through an embedded calculation approach using the Density Functional Theory and compared with previous experiments. For all tested complexes, the B3LYP functional is in good agreement with the experiments for the hyperfine coupling constants. For the electronic g-tensor calculations, performed using the coupled perturbed SCF theory, some discrepancies were found, and the best agreements with the experimental values were achieved by the B3LYP functional.

Theoretical investigation of the EPR hyperfine coupling constants in amino derivatives

The HFCCs of the radical cations of a series of amines have been determined at different levels of approximation including the CISD, QCISD, and CCSD ab initio correlated methods and density functional theory approaches employing the B3LYP, PBE0, BHandHLYP, TPSS, and BLYP exchange-correlation functionals. Although quantitative differences with respect to experimental data have been noticed, these are mostly systematic within a given class of N and H atoms. As a consequence, these different levels of theory are reliable in most cases to account for the substituent and structure effects on the HFCCs of amines. Linear regression fits have then been performed to reach quantitative agreement between the theoretical and experimental values. This has finally been substantiated by considering the EPR signal of the recently synthesized radical cations of two derivatives of [10-(4-aminophenyl)-9-anthryl]aniline as well as in confirming a recent assignment of the EPR signal of n-propylamine.

Time-dependent density functional theory with the generalized restricted-unrestricted approach

This work presents the derivation, implementation, and first applications of the generalized restricted-unrestricted method based on the density functional Kohn-Sham formalism. By using a spin-restricted Kohn-Sham representation for the reference state the well-known spin contamination problem is avoided, while the unrestricted representation of the perturbation response retains a proper description of spin polarization. The formulation is a generalization of our previous implementation of the restricted-unrestricted method [Z. Rinkevicius et al., J. Chem. Phys. 121, 7614 (2004)], as it accounts for the full unrestricted response instead of describing the spin polarization in terms of triplet operators only. The purpose of this paper is to investigate the role of the generalization employed and demonstrate its numerical performance. For this purpose we focus on isotropic hyperfine coupling constants of a set of organic radicals and transition metal compounds. For both classes of molecules we observe that the effect of neglecting singlet excitation operators in the response part of restricted-unrestricted formalism changes the calculated hyperfine coupling parameters by a few MHz. The obtained results confirm the validity of the approximation used in the simplified restricted-unrestricted approach, i.e. spin polarization can, in most cases, be adequately described by an account of only triplet operators in the response term.

Time-dependent density-functional theory calculations of triplet-triplet absorption

We present density-functional theory calculations of triplet-triplet absorption by three different approaches based on time-dependent density-functional theory (DFT): unrestricted DFT linear response, open-shell restricted DFT linear response applied to the triplet state, and quadratic response with triplet excitations applied to the ground state. Comparison is also made with corresponding results obtained by Hartree-Fock and multiconfiguration self-consistent-field response theory. Two main conclusions concerning triplet-triplet transitions are drawn in this study: First, the very good agreement between unrestricted and restricted DFT results indicates that spin contamination of the triplet state is not a serious problem when computing triplet-triplet spectra of common organic molecules. Second, DFT response calculations of triplet-triplet transitions can be affected by triplet instability problems, especially for the combination of DFT quadratic response with functionals containing fractional exact Hartree-Fock exchange.