Unified construction of relativistic Hamiltonians

It is shown that the four-component (4C), quasi-four-component (Q4C), and exact two-component (X2C) relativistic Hartree-Fock equations can be implemented in a unified manner by making use of the atomic nature of the small components of molecular 4-spinors. A model density matrix approximation can first be invoked for the small-component charge/current density functions, which gives rise to a static, pre-molecular mean field to be combined with the one-electron term. As a result, only the nonrelativistic-like two-electron term of the 4C/Q4C/X2C Fock matrix needs to be updated during the iterations. A "one-center small-component" approximation can then be invoked in the evaluation of relativistic integrals, that is, all atom-centered small-component basis functions are regarded as extremely localized near the position of the atom to which they belong such that they have vanishing overlaps with all small- or large-component functions centered at other nuclei. Under these approximations, the 4C, Q4C, and X2C mean-field and many-electron Hamiltonians share precisely the same structure and accuracy. Beyond these is the effective quantum electrodynamics Hamiltonian that can be constructed in the same way. Such approximations lead to errors that are orders of magnitude smaller than other sources of errors (e.g., truncation errors in the one- and many-particle bases as well as uncertainties of experimental measurements) and are, hence, safe to use for whatever purposes. The quaternion forms of the 4C, Q4C, and X2C equations are also presented in the most general way, based on which the corresponding Kramers-restricted open-shell variants are formulated for "high-spin" open-shell systems.

Noncollinear and Spin-Flip TDDFT in Multicollinear Approach

Time-dependent density functional theory (TDDFT) is one of the most important tools for investigating the excited states of electrons. The TDDFT calculation for spin-conserving excitation, where collinear functionals are sufficient, has obtained great success and has become routine. However, TDDFT for noncollinear and spin-flip excitations, where noncollinear functionals are needed, is less widespread and still a challenge nowadays. This challenge lies in the severe numerical instabilities that root in the second-order derivatives of commonly used noncollinear functionals. To be free from this problem radically, noncollinear functionals with numerical stable derivatives are desired, and our recently developed approach, called the multicollinear approach, provides an option. In this work, the multicollinear approach is implemented in noncollinear and spin-flip TDDFT, and prototypical tests are given.

Exact Two-Component TDDFT with Simple Two-Electron Picture-Change Corrections: X-ray Absorption Spectra Near L- and M-Edges of Four-Component Quality at Two-Component Cost

X-ray absorption spectroscopy (XAS) has gained popularity in recent years as it probes matter with high spatial and elemental sensitivities. However, the theoretical modeling of XAS is a challenging task since XAS spectra feature a fine structure due to scalar (SC) and spin-orbit (SO) relativistic effects, in particular near L and M absorption edges. While full four-component (4c) calculations of XAS are nowadays feasible, there is still interest in developing approximate relativistic methods that enable XAS calculations at the two-component (2c) level while maintaining the accuracy of the parent 4c approach. In this article we present theoretical and numerical insights into two simple yet accurate 2c approaches based on an (extended) atomic mean-field exact two-component Hamiltonian framework, (e)amfX2C, for the calculation of XAS using linear eigenvalue and damped response time-dependent density functional theory (TDDFT). In contrast to the commonly used one-electron X2C (1eX2C) Hamiltonian, both amfX2C and eamfX2C account for the SC and SO two-electron and exchange-correlation picture-change (PC) effects that arise from the X2C transformation. As we demonstrate on L- and M-edge XAS spectra of transition metal and actinide compounds, the absence of PC corrections in the 1eX2C approximation results in a substantial overestimation of SO splittings, whereas (e)amfX2C Hamiltonians reproduce all essential spectral features such as shape, position, and SO splitting of the 4c references in excellent agreement, while offering significant computational savings. Therefore, the (e)amfX2C PC correction models presented here constitute reliable relativistic 2c quantum-chemical approaches for modeling XAS.

Intramolecular Hydrogen Bond in Pyridine Schiff Bases as Ancillary Ligands of Re(I) Complexes Is a Switcher between Visible and NIR Emissions: A Relativistic Quantum Chemistry Study

Rhenium(I) tricarbonyl complexes have been described as suitable fluorophores, particularly for biological applications. fac-[Re(CO)₃(N,N)L]^(0 or 1+) complexes, where N,N is a substituted dinitrogenated ligand (bipyridine or derivatives with relatively small substituents) and L the ancillary ligand [a pyridine Schiff base (PSB) harboring an intramolecular hydrogen bond (IHB)], have presented promissory results concerning their use as fluorophores, especially for walled cells (i.e., bacteria and fungi). In this work, we present a relativistic theoretical analysis of two series of fac-[Re(CO)₃(N,N)PSB]¹⁺ complexes to predict the role of the IHB in the ancillary ligand concerning their photophysical behavior. N,N corresponds to 2,2'-bipyridine (bpy) (series A) or 4,4'-bis(ethoxycarbonyl)-2,2'-bipyridine (deeb) (series B). We found that all the complexes present absorption in the visible light range. In addition, complexes presenting a PSB with an IHB exhibit luminescent emission suitable for biological purposes: large Stokes shift, emission in the range of 600-700 nm, and τ in the order of 10^-2 to 10^-3 s. By contrast, complexes with PSB lacking the IHB show a predicted emission with the lowest triplet excited-state energy entering the NIR region. These results suggest a role of the IHB as an important switcher between visible and NIR emissions in this kind of complexes. Since the PSB can be substituted to modulate the properties of the whole Re(I) complex, it will be interesting to explore whether other substitutions can also affect the photophysical properties, mainly the emission range.

Relativistic Kramers-Unrestricted Exact-Two-Component Density Matrix Renormalization Group

In this work we develop a variational relativistic density matrix renormalization group (DMRG) approach within the exact-two-component (X2C) framework (X2C-DMRG), using spinor orbitals optimized with the two-component relativistic complete active space self-consistent field. We investigate fine-structure splittings of p- (Ga, In, Tl) and d-block (Sc, Y, La) atoms and excitation energies of monohydride molecules (GeH, SnH, and TlH) with X2C-DMRG calculations using an all-electron relativistic Hamiltonian in a Kramers-unrestricted basis. We find that X2C-DMRG yields accurate ²P and ²D splittings compared to multireference configuration interaction with singles and doubles (MRCISD). We also investigated the degree of symmetry breaking in the atomic multiplets and convergence of electron correlation in the total energies. Symmetry breaking can be large in some cases (∼30 meV); however, increasing the number of renormalized block states m for the DMRG optimization recovers the symmetry breaking by several orders of magnitude. Encouragingly, we find the convergence of electron correlation to be close to MRCISDTQ5 quality. Relativistic X2C-DMRG approaches are important for cases where spin-orbit coupling is significant and the underlying reference wave function requires a large determinantal space. We are able to obtain quantitatively correct fine-structure splittings for systems up to 10¹⁹ number of determinants with traditional CI approaches, which are currently unfeasible to converge for the field.

Self-Adaptive Real-Time Time-Dependent Density Functional Theory for X-ray Absorptions

Real-time time-dependent density functional theory (RT-TDDFT) can in principle access the whole absorption spectrum of a many-electron system exposed to a narrow pulse. However, this requires an accurate and efficient propagator for the numerical integration of the time-dependent Kohn-Sham equation. While a low-order time propagator is already sufficient for the low-lying valence absorption spectra, it is no longer the case for the X-ray absorption spectra (XAS) of systems composed even only of light elements, for which the use of a high-order propagator is indispensable. It is then crucial to choose a largest possible time step and a shortest possible simulation time, so as to minimize the computational cost. To this end, we propose here a robust AutoPST approach to determine automatically (Auto) the propagator (P), step (S), and time (T) for relativistic RT-TDDFT simulations of XAS.

Evaluation of molecular photophysical and photochemical properties using linear response time-dependent density functional theory with classical embedding: Successes and challenges

Time-dependent density functional theory (TDDFT) based approaches have been developed in recent years to model the excited-state properties and transition processes of the molecules in the gas-phase and in a condensed medium, such as in a solution and protein microenvironment or near semiconductor and metal surfaces. In the latter case, usually, classical embedding models have been adopted to account for the molecular environmental effects, leading to the multi-scale approaches of TDDFT/polarizable continuum model (PCM) and TDDFT/molecular mechanics (MM), where a molecular system of interest is designated as the quantum mechanical region and treated with TDDFT, while the environment is usually described using either a PCM or (non-polarizable or polarizable) MM force fields. In this Perspective, we briefly review these TDDFT-related multi-scale models with a specific emphasis on the implementation of analytical energy derivatives, such as the energy gradient and Hessian, the nonadiabatic coupling, the spin-orbit coupling, and the transition dipole moment as well as their nuclear derivatives for various radiative and radiativeless transition processes among electronic states. Three variations of the TDDFT method, the Tamm-Dancoff approximation to TDDFT, spin-flip DFT, and spin-adiabatic TDDFT, are discussed. Moreover, using a model system (pyridine-Ag20 complex), we emphasize that caution is needed to properly account for system-environment interactions within the TDDFT/MM models. Specifically, one should appropriately damp the electrostatic embedding potential from MM atoms and carefully tune the van der Waals interaction potential between the system and the environment. We also highlight the lack of proper treatment of charge transfer between the quantum mechanics and MM regions as well as the need for accelerated TDDFT modelings and interpretability, which calls for new method developments.

Polarizable Embedding Complex Polarization Propagator in Four- and Two-Component Frameworks

Explicit embedding methods combined with the complex polarization propagator (CPP) enable the modeling of spectroscopy for increasingly complex systems with a high density of states. We present the first derivation and implementation of the CPP in four- and exact-two-component (X2C) polarizable embedding (PE) frameworks. We denote the developed methods PE-4c-CPP and PE-X2C-CPP, respectively. We illustrate the methods by estimating the solvent effect on ultraviolet-visible (UV-vis) and X-ray atomic absorption (XAS) spectra of [Rh(H₂O)₆]³⁺ and [Ir(H₂O)₆]³⁺ immersed in aqueous solution. We moreover estimate solvent effects on UV-vis spectra of a platinum complex that can be photochemically activated (in water) to kill cancer cells. Our results clearly show that the inclusion of the environment is required: UV-vis and (to a lesser degree) XAS spectra can become qualitatively different from vacuum calculations. Comparison of PE-4c-CPP and PE-X2C-CPP methods shows that X2C essentially reproduces the solvent effect obtained with the 4c methods.

Relativistic nonorthogonal configuration interaction: application to L2,3-edge X-ray spectroscopy

In this article, we develop a relativistic exact-two-component nonorthogonal configuration interaction (X2C-NOCI) for computing L-edge X-ray spectra. This article to our knowledge is the first time NOCI has been used for relativistic wave functions. A set of molecular complexes, including SF₆, SiCl₄ and [FeCl₆]^3-, are used to demonstrate the accuracy and computational scaling of the X2C-NOCI method. Our results suggest that X2C-NOCI is able to satisfactorily capture the main features of the L_2,3-edge X-ray absorption spectra. Excitations from the core require a large amount of orbital relaxation to yield reasonable energies and X2C-NOCI allows us to treat orbital optimization explicitly. However, the cost of computing the nonorthogonal coupling is higher than in conventional CI. Here, we propose an improved integral screening using overlap-scaled density combined with a continuous measure of the generalized Slater-Condon rules that allows us to estimate if an element is zero before attempting a two-electron integral contraction.

Exact-Two-Component Multiconfiguration Pair-Density Functional Theory

Molecules containing late-row elements exhibit large relativistic effects. To account for both relativistic effects and electron correlation in a computationally inexpensive way, we derived a formulation of multiconfiguration pair-density functional theory with the relativistic exact-two-component Hamiltonian (X2C-MC-PDFT). In this new method, relativistic effects are included during variational optimization of a reference wave function by exact-two-component complete active-space self-consistent-field (X2C-CASSCF) theory, followed by an energy evaluation using pair-density functional theory. Benchmark studies of excited-state and ground-state fine-structure splitting of atomic species show that X2C-MC-PDFT can significantly improve the X2C-CASSCF results by introducing additional state-specific electron correlation.

Design, Synthesis, and Evaluation of Near-Infrared Fluorescent Molecules Based on 4H-1-Benzopyran Core

A series of novel fluorescent 4H-1-benzopyrans was designed and developed as near-infrared fluorescent molecules with a compact donor–acceptor-donor architecture. Spectral intensity of the fluorescent molecules M-1, M-2, M-3 varied significantly with the increasing polarities of solvents, where M-3 showed high viscosity sensitivity in glycerol-ethanol system with a 3-fold increase in emission intensity. Increasing concentrations of compound M-3 to 5% BSA in PBS elicited a 4-fold increase in fluorescence intensity, exhibiting a superior environmental sensitivity. Furthermore, the in vitro cellular uptake behavior and CLSM assay of cancer cell lines demonstrated that M-3 could easily enter the cell nucleus and bind to proteins with low toxicity. Therefore, the synthesized near-infrared fluorescent molecules could provide a new direction for the development of optical imaging probes and potential further drugs.

NAC-TDDFT: Time-Dependent Density Functional Theory for Nonadiabatic Couplings

First-order nonadiabatic coupling (NAC) matrix elements (fo-NACMEs) are the basic quantities in theoretical descriptions of electronically nonadiabatic processes that are ubiquitous in molecular physics and chemistry. Given the large size of systems of chemical interests, time-dependent density functional theory (TDDFT) is usually the first choice of methods. However, the lack of many-electron wave functions in TDDFT renders the formulation of NAC-TDDFT for fo-NACMEs conceptually difficult. Because of this, various variants of NAC-TDDFT have been proposed in the literature from different standing points, including the Hellmann-Feynman-like expression and auxiliary/pseudo-wave function (AWF)-, equation-of-motion (EOM)-, and time-dependent perturbation theory (TDPT)-based formulations. Based on critical analyses, the following conclusions are made here: (1) The Hellmann-Feynman-like expression, which is rooted in exact wave function theory, is hardly useful due to huge demand on basis sets. (2) Although most popular, the AWF variants of NAC-TDDFT are not theoretically founded and become ambiguous particularly for the fo-NACMEs between two excited states, although they do agree with the EOM and TDPT variants under the Tamm-Dancoff approximation. (3) The TDPT variant of NAC-TDDFT is theoretically most rigorous but suffers from numerical instabilities on the one hand and does not differ to a significant extent from the EOM variant on the other hand. (4) As such, the EOM variant of NAC-TDDFT for the fo-NACMEs between the ground and excited states and between two excited states is solely the right choice in practice. These formal analyses are fully supported by numerical experimentations, taking azulene as a showcase. The proper implementation of the EOM variant of NAC-TDDFT is also highlighted, showing that the fo-NACMEs between the ground and excited states and between two excited states are computationally very much the same as the analytic energy gradients of DFT and TDDFT, respectively. Possible future developments of the EOM variant of NAC-TDDFT are also highlighted. Its extensions to spin-adapted open-shell TDDFT and proper treatment of spin-orbit couplings (which are another source of force for electronically nonadiabatic processes) are particularly warranted in the near future.

Exact-two-component block-localized wave function: A simple scheme for the automatic computation of relativistic ΔSCF

Block-localized wave function is a useful method for optimizing constrained determinants. In this article, we extend the generalized block-localized wave function technique to a relativistic two-component framework. Optimization of excited state determinants for two-component wave functions presents a unique challenge because the excited state manifold is often quite dense with degenerate states. Furthermore, we test the degree to which certain symmetries result naturally from the ΔSCF optimization such as time-reversal symmetry and symmetry with respect to the total angular momentum operator on a series of atomic systems. Variational optimizations may often break the symmetry in order to lower the overall energy, just as unrestricted Hartree-Fock breaks spin symmetry. Overall, we demonstrate that time-reversal symmetry is roughly maintained when using Hartree-Fock, but less so when using Kohn-Sham density functional theory. Additionally, maintaining total angular momentum symmetry appears to be system dependent and not guaranteed. Finally, we were able to trace the breaking of total angular momentum symmetry to the relaxation of core electrons.

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.

Spin multiplicity effects in doublet versus singlet emission: the photophysical consequences of a single electron

Ambient-stable fluorescent radicals have recently emerged as promising luminescent materials; however, tailoring their properties has been difficult due to the limited photophysical understanding of open-shell organic systems. Here we report the experimental and computational analysis of a redox pair of π-conjugated fluorescent molecules that differ by one electron. A π-dication (DC) and π-radical cation (RC) demonstrate different absorption spectra, but similar red emission (λ_emiss,max = ∼630 nm), excitation maxima (λ_exc,max = ∼530 nm), fluorescence lifetimes (1–10 ns), and even excited-state (non-emissive) lifetimes when measured by transient absorption spectroscopy. Despite their experimental similarities, time-dependent density functional theory (TDDFT) studies reveal that DC and RC emission mechanisms are distinct and rely on different electronic transitions. Excited-state reorganization occurs by hole relaxation in singlet DC, while doublet RC undergoes a Jahn-Teller distortion by bending its π-backbone in order to facilitate spin-pairing between singly occupied molecular orbitals. This relationship between the excited-state dynamics of RC and its π-backbone geometry illustrates a potential strategy for developing π-conjugated radicals with new emission properties. Additionally, by comparing TDDFT and CIS (configuration interaction singles) excitations, we show that unrestricted TDDFT accurately reproduces experimental absorption spectra and provides an opportunity to examine the relaxed excited-state properties of large open-shell molecules like RC.

Experimental and computational studies reveal mechanistic differences in the photophysics of an open- versus closed-shell π-conjugated redox pair.

Chemical Bonding: The Journey from Miniature Hooks to Density Functional Theory

Our modern understanding of chemistry is predicated upon bonding interactions between atoms and ions resulting in the assembly of all of the forms of matter that we encounter in our daily life. It was not always so. This review article traces the development of our understanding of bonding from prehistory, through the debates in the 19th century C.E. bearing on valence, to modern quantum chemical models and beyond.

Spin-orbit coupling and vibronic transitions of Ce(C3H4) and Ce(C3H6) formed by the Ce reaction with propene: Mass-analyzed threshold ionization and relativistic quantum computation

A Ce atom reaction with propene is carried out in a pulsed laser vaporization molecule beam source. Several Ce-hydrocarbon species formed by the C-H and C-C bond activation of propene are observed by time-of-flight mass spectrometry, and Ce(C₃H_n) (n = 4 and 6) are characterized by mass-analyzed threshold ionization (MATI) spectroscopy and density functional theory, multiconfiguration, and relativistic quantum chemical calculations. The MATI spectrum of each species consists of two vibronic band systems, each with several vibronic bands. Ce(C₃H₆) is identified as an inserted species with Ce inserting into an allylic C-H bond of propene and Ce(C₃H₄) as a metallocycle through 1,2-vinylic dehydrogenation. Both species have a C_s structure with the Ce 4f¹6s¹ ground valence electron configuration in the neutral molecule and the Ce 4f¹ configuration in the singly charged ion. The two vibronic band systems observed for each species are attributed to the ionization of two pairs of the lowest spin-orbit coupled states with each pair being nearly degenerate.

Natural transition orbitals for complex two-component excited state calculations

While the natural transition orbital (NTO) method has allowed electronic excitations from time-dependent Hartree-Fock and density functional theory to be viewed in a traditional orbital picture, the extension to multicomponent molecular orbitals such as those used in relativistic two-component methods or generalized Hartree-Fock (GHF) or generalized Kohn-Sham (GKS) is less straightforward due to mixing of spin-components and the inherent inclusion of spin-flip transitions in time-dependent GHF/GKS. An extension of single-component NTOs to the two-component framework is presented, in addition to a brief discussion of the practical aspects of visualizing two-component complex orbitals. Unlike the single-component analog, the method explicitly describes the spin and frequently obtains solutions with several significant orbital pairs. The method is presented using calculations on a mercury atom and a CrO₂ Cl₂ complex.

Innovative Organic Electroluminescent Materials with a Doublet Ground State: A Theoretical Investigation

The displaced and distorted harmonic oscillator model, which has been proven to be appropriate in calculating vibronic spectra, is employed to treat the emission spectrum of title molecules in combination with a thermal vibration correlation function. The calculated results indicate that the main peak of the emission spectrum is visibly impacted by the normal modes with lower frequencies and that the shoulder peak is originated from the middle-frequency modes. On the level of time-dependent density functional theory (TDDFT), the calculated fluorescence lifetimes of TTM-3NCz and TTM-3PCz are 22.1 and 26.0 ns, respectively, which happen to coincide with the observed values of TTM-3NCz (17.2 ns) and TTM-3PCz (21.2 ns). The above data indicate that both the calculated radiative decay rates are reasonable at room temperature. Furthermore, we investigate the influence of the Duschinsky effect on the fluorescence quantum efficiency (FQE). When it is considered, the predicted FQE of the TTM-3NCz molecule is only 0.11%, and the observed value (49% in toluene) deviates significantly. If we ignore the Duschinsky effect, the FQE of TTM-3NCz increases dramatically to 41.8%. For the TTM-3PCz molecule (the FQE is 46% in toluene), the calculated FQE is 0.042% with the Duschinsky effect and increases to 45.2% without the Duschinsky effect. This phenomenon might be related to external factors and the nature of the TDDFT only considering a single configuration. In addition, the fluorescent properties of the fluorinated TTM-3NCz molecules are studied predictably. The obtained results show that the perfluorinated TTM-3NCz shows better luminous performance due to larger oscillator strength. Finally, the dimers, which are composed of both single title molecules, are explored theoretically to determine how they impact the fluorescent property; however, the effect can be nearly eliminated because of the small binding energies.

Toward the evaluation of intersystem crossing rates with variational relativistic methods

The change in electronic state from one spin multiplicity to another, known as intersystem crossing, occurs in molecules via the relativistic phenomenon of spin-orbit coupling. Current means of estimating intersystem crossing rates rely on the perturbative evaluation of spin-orbit coupling effects. This perturbative approach, valid in lighter atoms where spin-orbit coupling is weaker, is expected to break down for heavier elements where relativistic effects become dominant. Methods which incorporate spin-orbit effects variationally, such as the exact-two-component (X2C) method, will be necessary to treat this strong-coupling regime. We present a novel procedure which produces a diabatic basis of spin-pure electronic states coupled by spin-orbit terms, generated from fully variational relativistic calculations. This method is implemented within X2C using time-dependent density-functional theory and is compared to results from a perturbative relativistic study in the weak spin-orbit coupling regime. Additional calculations on a more strongly spin-orbit-coupled [UO₂Cl₄]^2- complex further illustrate the strengths of this method. This procedure will be valuable in the estimation of intersystem crossing rates within strongly spin-coupled species.

Spin trapping and flipping in FeCO through relativistic electron dynamics

Transition metal compounds are very versatile, and their characteristics can differ profoundly depending on their electronic structure. Compounds in which a spin transition from a low-spin to a high-spin state can be achieved through means of an optical excitation are particularly intriguing, as a controlled spin-flip opens promising avenues in areas such as sensing, information technology, molecular switches and energy technology. The fundamental mechanisms in spin crossover and spin transitions remain unanswered, due to the complexity of electronic structure and interplay of relativistic effects. Presented here is a new approach that allows the first direct study of spin flip dynamics through a mapping of spin-mixed to spin-pure states. The method is applied to FeCO and addresses the spin-flip dynamics during a spin transition. Wave packets that combine different spin states are generated through optical excitation and relevant mechanisms in optically triggered spin transitions are discussed.