One-electron spin-orbit contribution by effective nuclear charges

Thio and Seleno-Psoralens as Efficient Triplet Harvesting Photosensitizers for Photodynamic Therapy

The Psoralen (Pso) molecule finds extensive applications in photo-chemotherapy, courtesy of its triplet state forming ability. Sulfur and selenium replacement of exocyclic carbonyl oxygen of organic chromophores foster efficient triplet harvesting with near unity triplet quantum yield. These triplet-forming photosensitizers are useful in Photodynamic Therapy (PDT) applications for selective apoptosis of cancer cells. In this work, we have critically assessed the effect of the sulfur and selenium substitution at the exocyclic carbonyl (TPso and SePso, respectively) and endocyclic oxygen positions of Psoralen. It resulted in a significant redshifted absorption spectrum to access the PDT therapeutic window with increased oscillator strength. The reduction in singlet-triplet energy gap and enhancement in the spin-orbit coupling values increase the number of intersystem crossing (ISC) pathways to the triplet manifold, which shortens the ISC lifetime from 10-5 s for Pso to 10-8 s for TPso and 10-9 s for SePso. The intramolecular photo-induced electron transfer process, a competitive pathway to ISC, is also considerably curbed by exocyclic functionalizations. In addition, a maximum of 115 GM of two-photon absorption (2PA) with IR absorption (660-1050 nm) confirms that the Psoralen skeleton can be effectively tweaked via single chalcogen atom replacement to design a suitable PDT photosensitizer.

Unexpected longer T1 lifetime of 6-sulfur guanine than 6-selenium guanine: the solvent effect of hydrogen bonds to brake the triplet decay

The decay of the T1 state to the ground state is an essential property of photosensitizers because it decides the lifetime of excited states and, thus, the time window for sensitization. The sulfur/selenium substitution of carbonyl groups can red-shift absorption spectra and enhance the triplet yield because of the large spin-orbit coupling, modifying nucleobases to potential photosensitizers for various applications. However, replacing sulfur with selenium will also cause a much shorter T1 lifetime. Experimental studies found that the triplet decay rate of 6-seleno guanine (6SeGua) is 835 times faster than that of 6-thio guanine (6tGua) in aqueous solution. In this work, we reveal the mechanism of the T1 decay difference between 6SeGua and 6tGua by computing the activation energy and spin-orbit coupling for rate calculation. The solvent effect of water is treated with explicit microsolvation and implicit solvent models. We find that the hydrogen bond between the sulfur atom of 6tGua and the water molecule can brake the triplet decay, which is weaker in 6SeGua. This difference is crucial to explain the relatively long T1 lifetime of 6tGua in an aqueous solution. This insight emphasizes the role of solvents in modulating the excited state dynamics and the efficiency of photosensitizers, particularly in aqueous environments.

Role of excited-state hydrogen bonding in CO2 photoreduction catalyzed by sodium magnesium chlorophyll

In this paper, we report a joint experimental and computational study to elaborate the mechanism for the photocatalytic CO2 reduction reaction (CO2RR). Experimental results indicate that the catalyst (sodium magnesium chlorophyll, MgChlNa2), which has a well-defined structure for calculation and understanding, can achieve the photoreduction of CO2 to CO only using water as a dispersant, without adding any photosensitizer or sacrificial agent. Subsequently, a series of structural models of the hydrogen-bonded complexes of the catalyst were constructed and outlined via utilizing density functional theory (DFT) calculations, including photophysical and photochemical processes. The results confirm that the rate-limiting step of the whole CO2RR was the intersystem crossing process. The electron and proton transfers involved in photophysical and photochemical processes are induced by hydrogen bonds in the excited states. The combination of experiments and calculations will provide an important reference for the design of high-efficiency photocatalysts in the photocatalytic CO2RR.

Nonadiabatic Dynamics Simulations for Photoinduced Processes in Molecules and Semiconductors: Methodologies and Applications

Nonadiabatic dynamics (NAMD) simulations have become powerful tools for elucidating complicated photoinduced processes in various systems from molecules to semiconductor materials. In this review, we present an overview of our recent research on photophysics of molecular systems and periodic semiconductor materials with the aid of ab initio NAMD simulation methods implemented in the generalized trajectory surface-hopping (GTSH) package. Both theoretical backgrounds and applications of the developed NAMD methods are presented in detail. For molecular systems, the linear-response time-dependent density functional theory (LR-TDDFT) method is primarily used to model electronic structures in NAMD simulations owing to its balanced efficiency and accuracy. Moreover, the efficient algorithms for calculating nonadiabatic coupling terms (NACTs) and spin-orbit couplings (SOCs) have been coded into the package to increase the simulation efficiency. In combination with various analysis techniques, we can explore the mechanistic details of the photoinduced dynamics of a range of molecular systems, including charge separation and energy transfer processes in organic donor-acceptor structures, ultrafast intersystem crossing (ISC) processes in transition metal complexes (TMCs), and exciton dynamics in molecular aggregates. For semiconductor materials, we developed the NAMD methods for simulating the photoinduced carrier dynamics within the framework of the Kohn-Sham density functional theory (KS-DFT), in which SOC effects are explicitly accounted for using the two-component, noncollinear DFT method. Using this method, we have investigated the photoinduced carrier dynamics at the interface of a variety of van der Waals (vdW) heterojunctions, such as two-dimensional transition metal dichalcogenides (TMDs), carbon nanotubes (CNTs), and perovskites-related systems. Recently, we extended the LR-TDDFT-based NAMD method for semiconductor materials, allowing us to study the excitonic effects in the photoinduced energy transfer process. These results demonstrate that the NAMD simulations are powerful tools for exploring the photodynamics of molecular systems and semiconductor materials. In future studies, the NAMD simulation methods can be employed to elucidate experimental phenomena and reveal microscopic details as well as rationally design novel photofunctional materials with desired properties.

State Interaction Linear Response Time-Dependent Density Functional Theory with Perturbative Spin-Orbit Coupling: Benchmark and Perspectives

Spin-orbit coupling (SOC) is an important driving force in photochemistry. In this work, we develop a perturbative spin-orbit coupling method within the linear response time-dependent density function theory framework (TDDFT-SO). A full state interaction scheme, including singlet-triplet and triplet-triplet coupling, is introduced to describe not only the coupling between the ground and excited states, but also between excited states with all couplings between spin microstates. In addition, expressions to compute spectral oscillator strengths are presented. Scalar relativity is included variationally using the second-order Douglas-Kroll-Hess Hamiltonian, and the TDDFT-SO method is validated against variational SOC relativistic methods for atomic, diatomic, and transition metal complexes to determine the range of applicability and potential limitations. To demonstrate the robustness of TDDFT-SO for large-scale chemical systems, the UV-Vis spectrum of Au₂₅(SR)₁₈ ^- is computed and compared to experiment. Perspectives on the limitation, accuracy, and capability of perturbative TDDFT-SO are presented via analyses of benchmark calculations. Additionally, an open-source Python software package (PyTDDFT-SO) is developed and released to interface with the Gaussian 16 quantum chemistry software package to perform this calculation.

Generalized Ab Initio Nonadiabatic Dynamics Simulation Methods from Molecular to Extended Systems

Nonadiabatic dynamics simulation has become a powerful tool to describe nonadiabatic effects involved in photophysical processes and photochemical reactions. In the past decade, our group has developed generalized trajectory-based ab initio surface-hopping (GTSH) dynamics simulation methods, which can be used to describe a series of nonadiabatic processes, such as internal conversion, intersystem crossing, excitation energy transfer and charge transfer of molecular systems, and photoinduced nonadiabatic carrier dynamics of extended systems with and without spin-orbit couplings. In this contribution, we will first give a brief introduction to our recently developed methods and related numerical implementations at different computational levels. Later, we will present some of our latest applications in realistic systems, which cover organic molecules, biological proteins, organometallic compounds, periodic organic and inorganic materials, etc. Final discussion is given to challenges and outlooks of ab initio nonadiabatic dynamics simulations.

Enhanced Reverse Intersystem Crossing Promoted by Triplet Exciton-Photon Coupling

Polaritons are hybrid light-matter states formed via strong coupling between excitons and photons inside a microcavity, leading to upper and lower polariton (LP) bands splitting from the exciton. The LP has been applied to reduce the energy barrier of the reverse intersystem crossing (rISC) process from T₁, harvesting triplet energy for fluorescence through thermally activated delayed fluorescence. The spin-orbit coupling between T₁ and the excitonic part of the LP was considered as the origin for such an rISC transition. Here we propose a mechanism, namely, rISC promoted by the light-matter coupling (LMC) between T₁ and the photonic part of LP, which is originated from the ISC-induced transition dipole moment of T₁. This mechanism was excluded in previous studies. Our calculations demonstrate that the experimentally observed enhancement to the rISC process of the erythrosine B molecule can be effectively promoted by the LMC between T₁ and a photon. The proposed mechanism would substantially broaden the scope of the molecular design toward highly efficient cavity-promoted light-emitting materials and immediately benefit the illumination of related experimental phenomena.

An efficient implementation of spin-orbit coupling within the framework of semiempirical orthogonalization-corrected methods for ultrafast intersystem crossing dynamics

We implement spin-orbit coupling (SOC) within the framework of semiempirical orthogonalization-corrected methods (OMx). The excited-state wavefunction is generated from configuration interaction with single excitations (CIS). The SOC Hamiltonian in terms of the one-electron Breit-Pauli operator with effective nuclear charges is adopted in this work. Benchmark calculations show that SOCs evaluated using the OMx/CIS method agree very well with those obtained from time-dependent density functional theory. As a particularly attractive application, we incorporate SOCs between singlet and triplet states into Tully's fewest switches surface hopping algorithm to enable excited-state nonadiabatic dynamics simulations, treating internal conversion and intersystem crossing on an equal footing. This semiempirical dynamics simulation approach is applied to investigate ultrafast intersystem crossing processes in core-substituted naphthalenediimides.

Substitution effect on the nonradiative decay and trans → cis photoisomerization route: a guideline to develop efficient cinnamate-based sunscreens

Cinnamate derivatives are very useful as UV protectors in nature and as sunscreen reagents in daily life. They convert harmful UV energy to thermal energy through effective nonradiative decay (NRD) including trans → cis photoisomerization. However, the mechanism is not simple because different photoisomeirzation routes have been observed for different substituted cinnamates. Here, we theoretically examined the substitution effects at the phenyl ring of methylcinnamate (MC), a non-substituted cinnamate, on the electronic structure and the NRD route involving trans → cis isomerization based on time-dependent density functional theory. A systematic reaction pathway search using the single-component artificial force-induced reaction method shows that the very efficient photoisomerization route of MC can be essentially described as "1ππ* (trans) → 1nπ* → T1 (3ππ*) → S0 (trans or cis)". We found that for efficient 1ππ* (trans) → 1nπ* internal conversion (IC), MC should have the substituent at the appropriate position of the phenyl ring to stabilize the highest occupied π orbital. Substitution at the para position of MC slightly lowers the 1ππ* state energy and photoisomerization occurs via a slightly less efficient "1ππ* (trans) → 3nπ* → T1 (3ππ*) → S0 (trans or cis)" pathway. Substitution at the meta or ortho positions of MC significantly lowers the 1ππ* state energy so that the energy barrier of IC (1ππ* → 1nπ*) becomes very high. This substitution leads to a much longer 1ππ* state lifetime than that of MC and para-substituted MC, and a change in the dominant photoisomerization route to "1ππ* (trans) → C[double bond, length as m-dash]C bond twisting on 1ππ* → S0 (trans or cis)". As a whole, the "1ππ* → 1nπ*" IC observed in MC is the most important initial step for the rapid change of UV energy to thermal energy. We also found that the stabilization of the π orbital (i) minimizes the energy gap between 1ππ* and 1nπ* at the 1ππ* minimum and (ii) makes the 0-0 level of 1ππ* higher than 1nπ* as observed in MC. These MC-like relationships between the 1ππ* and 1nπ* energies should be ideal to maximize the "1ππ* → 1nπ*" IC rate constant according to Marcus theory.

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.

Theoretical survey of the photochemical deracemization mechanism of chiral allene 3-(3,3-dimethyl-1-buten-1-ylidene)-2-piperidinone

The photochemical deracemization mechanism of the chiral allene 3-(3,3-dimethyl-1-buten-1-ylidene)-2-piperidinone in the presence of photosensitizer thioxanthone is theoretically investigated by using a density functional theory (DFT) method.

Fine Modulation of the Higher-Order Excitonic States toward More Efficient Conversion from Upper-Level Triplet to Singlet

Hot exciton luminogens capable of harvesting nonemissive triplet excitons via reverse intersystem crossing from high-order triplet (hRISC) to singlet have great potential in high-efficiency fluorescent organic light-emitting diodes (OLEDs). Although spin-orbit coupling (SOC) is regarded as a key factor affecting the RISC process, its effects on hot exciton materials are poorly understood. Herein, we design and synthesize two blue-emitting hot exciton luminogens, PABP and PAIDO, to study this issue by modulating the excited-state properties. Theoretical and experimental research contributions demonstrate that a stronger SOC between energetically close S₁ (π-π*) and T_n (T₃, n-π*) of PAIDO gives rise to faster and more efficient hRISC in comparison to that of PABP, leading to a higher external quantum efficiency and a higher exciton utilization efficiency. Crucially, the experimentally measured hRISC rate (k_hRISC) of hot exciton materials is on the order of 10⁷ s^-1, which is much faster than that of the thermally activated delayed fluorescence materials.

Promoting Intersystem Crossing of a Fluorescent Molecule via Single Functional Group Modification

Pure light-atom organic phosphorescent molecules have been under scientific scrutiny because they are inexpensive, flexible, and environment friendly. The development of such materials, however, faces a bottleneck problem of intrinsically small spin-orbit couplings, which can be addressed by seeking a proper balance between intersystem crossing (ISC) and fluorescence rates. Using N-substituted naphthalimides as the prototype molecule, we applied chemical modifications with several electrophilic and nucleophilic functional groups, to approach the goal. The selected electron donating groups actively restrain the fluorescence, enabling an efficient ISC to the triplet manifold. Electron withdrawing groups do not change the luminescent properties of the parent species. The changes in ISC and fluorescence rates are related to the nature of the lowest singlet state, which changes from localized excitation into charge-transfer excitation upon configuration change of excited molecules. This finding opens an alternative strategy for designing pure light-atom organic phosphorescent molecules for emerging luminescent materials applications.

Different photoisomerization routes found in the structural isomers of hydroxy methylcinnamate

An experimental and theoretical study has been carried out to elucidate the nonradiative decay (NRD) and trans(E) → cis(Z) isomerization from the S1 (1ππ*) state of structural isomers of hydroxy methylcinnamate (HMC); ortho-, meta- and para-HMC (o-, m- and p-HMC). A low temperature matrix-isolation Fourier Transform Infrared (FTIR) spectroscopic study revealed that all the HMCs are cis-isomerized upon UV irradiation. A variety of laser spectroscopic methods have been utilized for jet-cooled gas phase molecules to investigate the vibronic structure and lifetimes of the S1 state, and to detect the transient state appearing in the NRD process. In p-HMC, the zero-point level of the S1 state decays as quickly as 9 ps. A transient electronic state reported by Tan et al. (Faraday Discuss. 2013, 163, 321-340) was reinvestigated by nanosecond UV-tunable deep UV pump-probe spectroscopy and was assigned to the T1 state. For m- and o-HMC, the lifetime at the zero-point energy level of S1 is 10 ns and 6 ns, respectively, but it becomes substantially shorter at an excess energy higher than 1000 cm-1 and 600 cm-1, respectively, indicating the onset of NRD. Different from p-HMC, no transient state (T1) was observed in m- nor o-HMC. These experimental results are interpreted with the aid of TDDFT calculations by considering the excited-state reaction pathways and the radiative/nonradiative rate constants. It is concluded that in p-HMC, the trans → cis isomerization proceeds via a [trans-S1 → 1nπ* → T1 → cis-S0] scheme. On the other hand, in o- and m-HMC, the isomerization proceeds via a [trans-S1 → twisting along the C[double bond, length as m-dash]C double bond by 90° on S1 → cis-S0] scheme. The calculated barrier height along the twisting coordinate agrees well with the observed onset of the NRD channel for both o- and m-HMC.

Mechanism of enhanced triplet decay of thionucleobase by glycosylation and rate-modulating strategies

Functionalization of the sugar group can be used to control the triplet decay rate of thionucleosides.

Incorporating spin–orbit effects into surface hopping dynamics using the diagonal representation: a linear-response time-dependent density functional theory implementation with applications to 2-thiouracil

Evaluation of SOC values employs Casida's wave functions and the Breit–Pauli spin–orbit Hamiltonian with effective charge approximation.

Combining the spin-separated exact two-component relativistic Hamiltonian with the equation-of-motion coupled-cluster method for the treatment of spin–orbit splittings of light and heavy elements

An accurate and efficient treatment of spin–orbit splittings has been achieved by combining the sf-X2C+soc-DKH1 Hamiltonian with the equation-of-motion coupled-cluster method.

A quantitative structure–property relationship (QSPR) study of singlet oxygen generation by pteridines

Singlet oxygen production quantum yields of pteridine photosensitizers were analyzed with the QSPR method. The ability of pterins and flavins to generate¹O₂in D₂O correlated withE_HOMOand electronegativity, as well as with the dipole moment and some other parameters.