Role of Triplet States in the Photodynamics of Aniline

Ultrafast Intersystem Crossing in Benzanthrone: Effect of Hydrogen Bonding and Viscosity

Understanding the intricate factors governing intersystem crossing (ISC) in aromatic carbonyl compounds remains a long-standing interest among researchers. This study unveils the crucial roles of vibration in influencing the ISC of a typical aromatic carbonyl chromophore, benzanthrone, and how hydrogen bonding and solvent viscosity affect these vibrations and, thus, the associated ISC kinetics. We demonstrate that for benzanthrone, the ISC is exceedingly facile in an aprotic solvent, while in protic solvents, the ISC is significantly suppressed through the formation of the hydrogen-bonded state. Moreover, in a high-viscosity medium, ISC is further retarded due to restrictions of volume-changing motions, which may assist ISC. Theoretical calculations revealed that the C═O bond vibration and specific out-of-plane vibrations accompanying a volume change could be the probable coordinates for ISC. These findings provide valuable insights for tailoring the excited-state behavior of carbonyl-functionalized materials for diverse applications in photocatalysis, organic electronics, and biomedicine.

Room Temperature Phosphorescence in Crystalline Iodinated Eumelanin Monomer

We report the room temperature phosphorescence upon iodination on a crystalline eumelanin monomer with shielded hydroxyl moieties, ethyl 5,6-dimethoxyindole-2-carboxylate (DMICE). Ultrafast intersystem crossing (ISC) is observed in the iodinated (IDMICE) as well as brominated (BDMICE) analogues of the eumelanin monomer derivative in solution. The triplet quantum yields (φT) and intersystem crossing rates (kISC) of the halogenated eumelanin derivatives areφ T B D M I C E ${{\phi{} }_{T}^{BDMICE}}$ =25.4±1.1 %;k I S C B D M I C E ${{k}_{ISC}^{BDMICE}}$ =1.95×109 s-1 andφ T I D M I C E ${{\phi{} }_{T}^{IDMICE}}$ =59.1±1.6 %;k I S C I D M I C E = ${{k}_{ISC}^{IDMICE}=}$ 1.36×1010 s-1, as monitored using transient absorption spectroscopy. Theoretical calculations based on nuclear ensemble method reveal that computed kISC and spin-orbit coupling matrix elements for eumelanin derivatives are larger for IDMICE relative to BDMICE. The halogen and π-π interactions, with distinct excitonic coupling and higher ISC rate promote phosphorescence in IDMICE molecular crystals. Accessing triplet excited states and resultant photoluminescence through structural modification of eumelanin scaffolds paves way for exploring the versatility of eumelanin-inspired molecules as bio-functional materials.

Generalized Born-Huang expansion under macroscopic quantum electrodynamics framework

Born-Huang expansion is the cornerstone for studying potential energy surfaces and non-adiabatic couplings (NACs) in molecular systems. However, the traditional approach is insufficient to describe the molecular system, which strongly interacts with quantum light. Inspired by the work by Schäfer et al., we develop the generalized Born-Huang expansion theory within a macroscopic quantum electrodynamics (QED) framework. The theory we present allows us to describe electromagnetic vacuum fluctuations in dielectric media and incorporate the effects of dressed photons (or polaritons) into NACs. With the help of the generalized Born-Huang expansion, we clearly classify electronic nuclear NACs, polaritonic nuclear NACs, and polaritonic electronic NACs. Furthermore, to demonstrate the advantage of the macroscopic QED framework, we estimate polaritonic electronic NACs without any free parameter, such as the effective mode volume, and demonstrate the distance dependence of the polaritonic electronic NACs in a silver planar system. In addition, we take a hydrogen atom in free space as an example and derive spontaneous emission rates from photonic electronic NACs (polaritonic electronic NACs are reduced to photonic electronic NACs). We believe that this work not only provides an avenue for the theoretical exploration of NACs in a nucleus-electron-polariton coupled system but also offers a more comprehensive understanding for molecules coupled with quantum light.

Heavy-Atom Quantum Tunnelling in Spin Crossovers of Nitrenes

We simulate two recent matrix-isolation experiments at cryogenic temperatures, in which a nitrene undergoes spin crossover from its triplet state to a singlet state via quantum tunnelling. We detail the failure of the commonly applied weak-coupling method (based on a linear approximation of the potentials) in describing these deep-tunnelling reactions. The more rigorous approach of semiclassical golden-rule instanton theory in conjunction with double-hybrid density-functional theory and multireference perturbation theory does, however, provide rate constants and kinetic isotope effects in good agreement with experiment. In addition, these calculations locate the optimal tunnelling pathways, which provide a molecular picture of the reaction mechanism. The reactions involve substantial heavy-atom quantum tunnelling of carbon, nitrogen and oxygen atoms, which unexpectedly even continues to play a role at room temperature.

IR absorption spectra of aniline cation, anilino radical, and phenylnitrene isolated in solid argon

Electron bombardment of aniline (PhNH₂) in an Ar matrix mainly generated the aniline cation (PhNH₂⁺), anilino (PhNH) and phenyl (Ph) radicals, and phenylnitrene (PhN). Further irradiation of the electron-bombarded matrix sample at 365 nm depleted PhNH₂⁺ and PhN, and resulted in the formation of PhNH₂, PhNH, and Ph. In separate experiments, irradiation of the PhNH₂/Ar matrix samples at 265 or 160 nm mainly generated PhNH and Ph radicals, but without the formation of PhNH₂⁺ and PhN. According to the observed photochemical behaviors, quantum-chemically predicted harmonic vibrational wavenumbers of each species, and the information reported in previous photodissociation studies, we unambiguously characterized the IR features of the aromatic species. The information of the vibrational fundamentals of PhNH is new and the formation mechanism is discussed.

Internal conversion and intersystem crossing dynamics based on coupled potential energy surfaces with full geometry-dependent spin-orbit and derivative couplings. Nonadiabatic photodissociation dynamics of NH3(A) leading to the NH(X3Σ-, a1Δ) + H2 channel

We simulate the photodissociation of NH₃ originating from its first excited singlet state S₁ into the NH₂ + H (radical) and NH + H₂ (molecular) channels. The states considered are the ground singlet state S₀, the first excited singlet state S₁ and the lowest-lying triplet state T₁, which permit for the first time a uniform treatment of the internal conversion and intersystem crossing. The simulations are based on a diabatic potential energy matrix (DPEM) of S₀, S₁ coupled by a conical intersection seam, as well as a potential energy surface (PES) for T₁ coupled by spin-orbit coupling (SOC) to the two singlet states. The DPEM and PES are fitted to ab initio electronic structure data (ESD) including energies, energy gradients, and derivative couplings. The DPEM also defines an adiabatic to diabatic state (AtD) transformation, which is used to transform the singular adiabatic SOC into a smooth function of the nuclear coordinates in the diabatic representation, allowing the diabatic SOC to be fit to an analytical functional form. ESD and SOC data obtained from these surfaces can serve as input for either quantum or semi-classical characterization of the nonadiabatic dynamics. Using the SHARC suite of programs, nonadiabatic simulations based on over 40 000 semi-classical trajectories assess the convergence of our results. The production of NH + H₂ is not direct, but is only achieved through a quasi-statistical dissociation mechanism after internal conversion to the ground electronic state. This leads to a much lower yield comparing with the main NH₂ + H channel. The NH(X³Σ^_) radical produced through the intersystem crossing from S₀ to T₁ is rare (∼0.2%) compared to NH(a¹Δ) due to the process being spin forbidden.

NAST: Nonadiabatic Statistical Theory Package for Predicting Kinetics of Spin-Dependent Processes

We present a nonadiabatic statistical theory (NAST) package for predicting kinetics of spin-dependent processes, such as intersystem crossings, spin-forbidden unimolecular reactions, and spin crossovers. The NAST package can calculate the probabilities and rates of transitions between the electronic states of different spin multiplicities. Both the microcanonical (energy-dependent) and canonical (temperature-dependent) rate constants can be obtained. Quantum effects, including tunneling, zero-point vibrational energy, and reaction path interference, can be accounted for. In the limit of an adiabatic unimolecular reaction proceeding on a single electronic state, NAST reduces to the traditional transition state theory. Because NAST requires molecular properties at only a few points on potential energy surfaces, it can be applied to large molecular systems, used with accurate high-level electronic structure methods, and employed to study slow nonadiabatic processes. The essential NAST input data include the nuclear Hessian at the reactant minimum, as well as the nuclear Hessians, energy gradients, and spin-orbit coupling at the minimum energy crossing point (MECP) between two states. The additional computational tools included in the NAST package can be used to extract the required input data from the output files of electronic structure packages, calculate the effective Hessian at the MECP, and fit the reaction coordinate for more advanced NAST calculations. We describe the theory, its implementation, and three examples of application to different molecular systems.

Insights into the hydrated electron generation from UV/aniline: Mechanism and quantum efficiency

The insight into the neglected reduction process accompanied by UV-based oxidation process may provide new ideas for the development of advanced oxidation and reduction technologies. In this study, aniline was comprehensively investigated as an unexpected indicator of hydrated electron (e_aq^-) under UV irradiation. Monochloroacetic acid (MCAA) was selected as the probe of e_aq^- and the balance of chloride ions indicated the reduction of MCAA. Further, laser flash photolysis experiments demonstrated the generation of e_aq^- in the UV/aniline process and the half-life period of formed e_aq^- was demonstrated to be 0.13 μs. The photolysis of aniline along with the decay of the excited state of aniline was responsible for the e_aq^- generation. Besides, the hydrogen atom (H) generated from the photolysis can subsequently reacted with OH^- to generate e_aq^-. The photolysis pathways of aniline were proposed by the results of GC-MS. Aniline was abstracted of H in solution to the formation of aniline radical (PhNH) or form aminophenol in three different isomers (orto-, meta- and para-aminophenol). Moreover, UV/aniline showed a higher reducing capacity of MCAA compared with other organic electron donors and sustained a highly reducing ability in a wide pH. And the calculation results of quantum efficiency (Φ) showed that excessive aniline was not conducive to the elevation of Φ. This study introduced a novel pathway of e_aq^- generation during the photolysis of aniline and provided a new perspective for e_aq^--based advanced reduction processes.

Electronic structure of strongly correlated systems: recent developments in multiconfiguration pair-density functional theory and multiconfiguration nonclassical-energy functional theory

Strong electron correlation plays an important role in transition-metal and heavy-metal chemistry, magnetic molecules, bond breaking, biradicals, excited states, and many functional materials, but it provides a significant challenge for modern electronic structure theory. The treatment of strongly correlated systems usually requires a multireference method to adequately describe spin densities and near-degeneracy correlation. However, quantitative computation of dynamic correlation with multireference wave functions is often difficult or impractical. Multiconfiguration pair-density functional theory (MC-PDFT) provides a way to blend multiconfiguration wave function theory and density functional theory to quantitatively treat both near-degeneracy correlation and dynamic correlation in strongly correlated systems; it is more affordable than multireference perturbation theory, multireference configuration interaction, or multireference coupled cluster theory and more accurate for many properties than Kohn–Sham density functional theory. This perspective article provides a brief introduction to strongly correlated systems and previously reviewed progress on MC-PDFT followed by a discussion of several recent developments and applications of MC-PDFT and related methods, including localized-active-space MC-PDFT, generalized active-space MC-PDFT, density-matrix-renormalization-group MC-PDFT, hybrid MC-PDFT, multistate MC-PDFT, spin–orbit coupling, analytic gradients, and dipole moments. We also review the more recently introduced multiconfiguration nonclassical-energy functional theory (MC-NEFT), which is like MC-PDFT but allows for other ingredients in the nonclassical-energy functional. We discuss two new kinds of MC-NEFT methods, namely multiconfiguration density coherence functional theory and machine-learned functionals.

This feature article overviews recent work on active spaces, matrix product reference states, treatment of quasidegeneracy, hybrid theory, density-coherence functionals, machine-learned functionals, spin–orbit coupling, gradients, and dipole moments.

Quantum-classical hybrid algorithm for the simulation of all-electron correlation

While chemical systems containing hundreds to thousands of electrons remain beyond the reach of quantum devices, hybrid quantum-classical algorithms present a promising pathway toward a quantum advantage. Hybrid algorithms treat the exponentially scaling part of the calculation-the static correlation-on the quantum computer and the non-exponentially scaling part-the dynamic correlation-on the classical computer. While a variety of algorithms have been proposed, the dependence of many methods on the total wave function limits the development of easy-to-use classical post-processing implementations. Here, we present a novel combination of quantum and classical algorithms, which computes the all-electron energy of a strongly correlated molecular system on the classical computer from the 2-electron reduced density matrix (2-RDM) evaluated on the quantum device. Significantly, we circumvent the wave function in the all-electron calculations by using density matrix methods that only require input of the statically correlated 2-RDM. Although the algorithm is completely general, we test it with two classical density matrix methods, the anti-Hermitian contracted Schrödinger equation (ACSE) and multiconfiguration pair-density functional theories, using the recently developed quantum ACSE method for simulating the statically correlated 2-RDM. We obtain experimental accuracy for the relative energies of all three benzyne isomers and thereby demonstrate the ability of the developed algorithm to achieve chemically relevant and accurate results on noisy intermediate-scale quantum devices.

Spin Crossover of Thiophosgene via Multidimensional Heavy-Atom Quantum Tunneling

The spin-crossover reaction of thiophosgene has drawn broad attention from both experimenters and theoreticians as a prime example of radiationless intramolecular decay via intersystem crossing. Despite multiple attempts over 20 years, theoretical predictions have typically been orders of magnitude in error relative to the experimentally measured triplet lifetime. We address the T₁ → S₀ transition by the first application of semiclassical golden-rule instanton theory in conjunction with on-the-fly electronic-structure calculations based on multireference perturbation theory. Our first-principles approach provides excellent agreement with the experimental rates. This was only possible because instanton theory goes beyond previous methods by locating the optimal tunneling pathway in full dimensionality and thus captures "corner cutting" effects. Since the reaction is situated in the Marcus inverted regime, the tunneling mechanism can be interpreted in terms of two classical trajectories, one traveling forward and one backward in imaginary time, which are connected by particle-antiparticle creation and annihilation events. The calculated mechanism indicates that the spin crossover is sped up by many orders of magnitude due to multidimensional quantum tunneling of the carbon atom even at room temperature.

Dipole Moment Calculations Using Multiconfiguration Pair-Density Functional Theory and Hybrid Multiconfiguration Pair-Density Functional Theory

The dipole moment is the molecular property that most directly indicates molecular polarity. The accuracy of computed dipole moments depends strongly on the quality of the calculated electron density, and the breakdown of single-reference methods for strongly correlated systems can lead to poor predictions of the dipole moments in those cases. Here, we derive the analytical expression for obtaining the electric dipole moment by multiconfiguration pair-density functional theory (MC-PDFT), and we assess the accuracy of MC-PDFT for predicting dipole moments at equilibrium and nonequilibrium geometries. We show that MC-PDFT dipole moment curves have reasonable behavior even for stretched geometries, and they significantly improve upon the CASSCF results by capturing more electron correlation. The analysis of a dataset consisting of 18 first-row transition-metal diatomics and 6 main-group polyatomic molecules with a multireference character suggests that MC-PDFT and its hybrid extension (HMC-PDFT) perform comparably to CASPT2 and MRCISD+Q methods and have a mean unsigned deviation of 0.2-0.3 D with respect to the best available dipole moment reference values. We explored the dependence of the predicted dipole moments upon the choice of the on-top density functional and active space, and we recommend the tPBE and hybrid tPBE0 on-top choices for the functionals combined with the moderate correlated-participating-orbitals scheme for selecting the active space. With these choices, the mean unsigned deviations (in debyes) of the calculated equilibrium dipole moments from the best estimates are 0.77 for CASSCF, 0.29 for MC-PDFT, 0.24 for HMC-PDFT, 0.28 for CASPT2, and 0.25 for MRCISD+Q. These results are encouraging because the computational cost of MC-PDFT or HMC-PDFT is largely reduced compared to the CASPT2 and MRCISD+Q methods.

Benzene, Toluene, and Monosubstituted Derivatives: Diabatic Nature of the Oscillator Strengths of S1 ← S0 Transitions

For benzene, toluene, aniline, fluorobenzene, and phenol, even sophisticated treatments of electron correlation, such as MRCI and XMS-CASPT2 calculations, show oscillator strengths typically lower than experiment. Inclusion of a simple pseudo-diabatization approach to perturb the S1 state with approximate vibronic coupling to the S2 state for each molecule results in more accurate oscillator strengths. Their absolute values agree better with experiment for all molecules except aniline. When the coupling between the S1 and S2 states is strong at the S0 geometry, the simple diabatization scheme performs less well with respect to the oscillator strengths relative to the adiabatic values. However, we expect the scheme to be useful in many cases where the coupling is weak to moderate (where the maximum component of the coupling has a magnitude less than 1.5 au). Such calculations give an insight into the effects of vibronic coupling of excited states on UV/vis spectra.