Toward the Understanding of the Photophysics and Photochemistry of 1-Nitronaphthalene under Solar Radiation: The First Theoretical Evidence of a Photodegradation Intramolecular Rearrangement Mechanism Involving the Triplet States

XMS-CASPT2//XMS-CASPT2 and XMS-CASPT2//CASSCF at comparison: The impact of dynamic correlation in the excited state optimization of nitronaphthalene

Here, analytical extended multi-state complete active space second-order perturbation method (XMS-CASPT2) gradients are used to rationalize the decreasing triplet quantum yield trend in 2-nitronaphthalene, 1-nitronaphthalene, and 2-methyl-1-nitronaphthalene, a series of nitro-substituted aromatic compounds. Comparison with the XMS-CASPT2//CASSCF (where CASSCF stands for complete active space self-consistent field method) results highlights the importance of dynamic correlation in geometry optimization and challenges the validity of an XMS-CASPT2//CASSCF approach: XMS-CASPT2 S1 minima leads to planar structures, while CASSCF optimizations trigger a pyramidalization of the nitro group. The XMS-CASPT2 results correlate the reported decreasing triplet quantum yield trend in these species to a decrease in S1 to T2 population transfer and an increase in S1-S0 decay, while no such correlation is observed when using XMS-CASPT2//CASSCF data.

A First Proposal on the Nitrobenzene Photorelease Mechanism of NO2 and Its Relation to NO Formation through a Roaming Mechanism

Despite the fact that NO2 is considered to be the main photoproduct of nitrobenzene photochemistry, no mechanism has ever been proposed to rationalize its formation. NO photorelease is instead a more studied process, probably due to its application in the drug delivery sector and the study of roaming mechanisms. In this contribution, a photoinduced mechanism accounting for the formation of NO2 in nitrobenzene is theorized based on CASPT2, CASSCF, and DFT electronic structure calculations and CASSCF classical dynamics. A triplet nπ* state is shown to evolve toward C-NO2 dissociation, being, in fact, the only low-lying excited state favoring such a deformation. Along the triplet dissociation path, the possibility to decay to the singlet ground state results in the frustration of the dissociation and in the recombination of the fragments, either back to the nitro or the nitrite isomer. The thermal decomposition of the latter to NO constitutes globally a roaming mechanism of NO formation.

Solvent-polarity dependence of ultrafast excited-state dynamics of trans-4-nitrostilbene

Ultrafast excited-state dynamics of the simplest nitrostilbenes, namely trans-4-nitrostilbene (t-NSB), was studied in solvents of various polarities with ultrafast broadband time-resolved fluorescence and transient absorption spectroscopies, and by quantum-chemical computations. The results revealed that the initially excited S1(ππ*) state deactivation dynamics is strongly influenced by the solvent polarity. Specifically, the t-NSB S1-state lifetime decreases by three orders of magnitude from ∼60 ps in high-polarity solvents to ∼60 fs in nonpolar solvents. The strong solvent-polarity dependence arises from the differences in dipole moments among the S1 and relevant states, including the major intersystem crossing (ISC) receiver triplet states, and therefore, the solvent polarity can modulate their relative energies and ISC rates. In nonpolar solvents, the sub-100 fs lifetime is due to a combination of efficient ISC and internal conversion. In medium-polarity solvents, the S1-state population decays via a competing ISC relaxation mechanism in a biphasic manner, and the ISC rates are found to obey the inverse energy gap law of the strong coupling case. In high-polarity solvents, the S1 state is stabilized to a much lower energy such that ISC becomes energetically infeasible, and the S1 state decays via barrier crossing along the torsion angle of the central ethylenic bond to the nonfluorescent perpendicular configuration. Regardless of the initial S1-state deactivation pathways in various solvents, the excited-state population is ultimately trapped in the metastable T1-state perpendicular configuration, at which a slower ISC occurs to bring the system to the ground state and bifurcate into either trans or cis form of NSB.

Overlooked Transformation of Nitrated Polycyclic Aromatic Hydrocarbons in Natural Waters: Role of Self-Photosensitization

Photochemical transformation is an important process that involves trace organic contaminants (TrOCs) in sunlit surface waters. However, the environmental implications of their self-photosensitization pathway have been largely overlooked. Here, we selected 1-nitronaphthalene (1NN), a representative nitrated polycyclic aromatic hydrocarbon, to study the self-photosensitization process. We investigated the excited-state properties and relaxation kinetics of 1NN after sunlight absorption. The intrinsic decay rate constants of triplet (³1NN*) and singlet (¹1NN*) excited states were estimated to be 1.5 × 10⁶ and 2.5 × 10⁸ s^-1, respectively. Our results provided quantitative evidence for the environmental relevance of ³1NN* in waters. Possible reactions of ³1NN* with various water components were evaluated. With the reduction and oxidation potentials of -0.37 and 1.95 V, ³1NN* can be either oxidized or reduced by dissolved organic matter isolates and surrogates. We also showed that hydroxyl (^•OH) and sulfate (SO₄^•-) radicals can be generated via the ³1NN*-induced oxidation of inorganic ions (OH^- and SO₄^2-, respectively). We further investigated the reaction kinetics of ³1NN* and OH^- forming ^•OH, an important photoinduced reactive intermediate, through complementary experimental and theoretical approaches. The rate constants for the reactions of ³1NN* with OH^- and 1NN with ^•OH were determined to be 4.22 × 10⁷ and 3.95 ± 0.01 × 10⁹ M^-1 s^-1, respectively. These findings yield new insights into self-photosensitization as a pathway for TrOC attenuation and provide more mechanistic details into their environmental fate.

Effect of electron-donating substitution on the triplet state reactivities of 1-nitronaphthalene

Nitro-polycyclic aromatic hydrocarbons (nitro-PAHs), often found in polluted air, are carcinogenic and mutagenic. The nitro group increases the spin-orbit coupling and results in the lowest excited triplet (T₁) on the picosecond time scale with a high yield. The electron-donating substituents have a significant influence on the photophysics and photochemistry of nitro-PAHs. We used transient absorption spectroscopy and kinetic analysis to investigate the reactivities of the T₁ state 1-methoxy-4-nitronaphthalene (³MeO-NN) and 1-methyl-4-nitronaphthalene (³Me-NN). The results show that the methoxy and methyl substitutions have a minor effect on their hydrogen abstraction and electron accepting abilities. The main distinction is their reaction rates towards protons. The second order rate constant of ³MeO-NN towards protons is three orders of magnitude greater than that of ³Me-NN, indicating that ³MeO-NN has a stronger hydrogen bond accepting ability. The kinetic analysis reveals that the dimer of 2,2,2-trifluoroethanol participates in the reaction with ³MeO-NN. These results suggest that the formation of the hydrogen-bonded complex is responsible for the unusually short lifetime of ³MeO-NN in methanol solution and the lack of hydrogen abstraction radicals during the decay of ³MeO-NN in methanol.

On the photorelease of nitric oxide by nitrobenzene derivatives: A CASPT2//CASSCF model

Nitroaromatic compounds can photorelease nitric oxide after UV absorption. The efficiency of the photoreaction depends on the molecular structure, and two features have been pointed out as particularly important for the yield of the process: the presence of methyl groups at the ortho position with respect to the nitro group and the degree of conjugation of the molecule. In this paper, we provide a theoretical characterization at the CASPT2//CASSCF (complete active space second-order perturbation theory//complete active space self-consistent field) level of theory of the photorelease of NO for four molecules derived from nitrobenzene through the addition of ortho methyl groups and/or the elongation of the conjugation. Our previously described mechanism obtained for the photorelease of NO in nitrobenzene has been adopted as a model for the process. According to this model, the process proceeds through a reactive singlet-triplet crossing (STC) region that the system can reach from the triplet 3(πOπ*) minimum. The energy barrier that must be surmounted in order to populate the reactive STC can be associated with the efficiency of the photoreaction. Here, the obtained results display clear differences in the efficiency of the photoreaction in the studied systems and can be correlated with experimental results. Thus, the model proves its ability to highlight the differences in the photoreaction efficiency for the nitroaromatic compounds studied here.

Assessing the Performances of CASPT2 and NEVPT2 for Vertical Excitation Energies

Methods able to simultaneously account for both static and dynamic electron correlations have often been employed, not only to model photochemical events but also to provide reference values for vertical transition energies, hence allowing benchmarking of lower-order models. In this category, both the complete-active-space second-order perturbation theory (CASPT2) and the N-electron valence state second-order perturbation theory (NEVPT2) are certainly popular, the latter presenting the advantage of not requiring the application of the empirical ionization-potential-electron-affinity (IPEA) and level shifts. However, the actual accuracy of these multiconfigurational approaches is not settled yet. In this context, to assess the performances of these approaches, the present work relies on highly accurate (±0.03 eV) aug-cc-pVTZ vertical transition energies for 284 excited states of diverse character (174 singlet, 110 triplet, 206 valence, 78 Rydberg, 78 n → π*, 119 π → π*, and 9 double excitations) determined in 35 small- to medium-sized organic molecules containing from three to six non-hydrogen atoms. The CASPT2 calculations are performed with and without IPEA shift and compared to the partially contracted (PC) and strongly contracted (SC) variants of NEVPT2. We find that both CASPT2 with IPEA shift and PC-NEVPT2 provide fairly reliable vertical transition energy estimates, with slight overestimations and mean absolute errors of 0.11 and 0.13 eV, respectively. These values are found to be rather uniform for the various subgroups of transitions. The present work completes our previous benchmarks focused on single-reference wave function methods ( J. Chem. Theory Comput. 2018, 14, 4360; J. Chem. Theory Comput. 2020, 16, 1711), hence allowing for a fair comparison between various families of electronic structure methods. In particular, we show that ADC(2), CCSD, and CASPT2 deliver similar accuracies for excited states with a dominant single-excitation character.

Solvent dielectric delimited nitro–nitrito photorearrangement in a perylenediimide derivative

The discovery of vibrant excited-state dynamics and distinctive photochemistry has established nitrated polycyclic aromatic hydrocarbons as an exhilarating class of organic compounds. Herein, we report the atypical photorearrangement of nitro-perylenediimide (NO₂-PDI) to nitrito-perylenediimide (ONO-PDI), triggered by visible-light excitation and giving rise to linkage isomers in the polar aprotic solvent acetonitrile. ONO-PDI has been isolated and unambiguously characterized using standard spectroscopic, spectrometric, and elemental composition techniques. Although nitritoaromatic compounds are conventionally considered to be crucial intermediates in the photodissociation of nitroaromatics, experimental evidence for this has not been observed heretofore. Ultrafast transient absorption spectroscopy combined with computational investigations revealed the prominence of a conformationally relaxed singlet excited-state (S^CR₁) of NO₂-PDI in the photoisomerization pathway. Theoretical transition state (TS) analysis indicated the presence of a six-membered cyclic TS, which is pivotal in connecting the S^CR₁ state to the photoproduct state. This article addresses prevailing knowledge gaps in the field of organic linkage isomers and provides a comprehensive understanding of the unprecedented photoisomerization mechanism operating in the case of NO₂-PDI.

The unprecedented photorearrangement of nitro-perylenediimide (NO₂-PDI) to nitrito-perylenediimide (ONO-PDI) is shown to occur through a cyclic six-membered transition state triggered by visible-light excitation.

Excited state dynamics and photochemistry of nitroaromatic compounds

Nitrated aromatic molecules have unique photoinduced channels. Due to the presence of oxygen-centered non-bonding orbitals, they can undergo sub-picosecond intersystem crossing showing one of the strongest couplings between the singlet and triplet manifolds among organic molecules. Several nitroaromatic compounds also have a distinctive nitric oxide photodissociation channel which occurs through a complex sequence of atom rearrangements and state changes. These remarkable processes have stimulated the attention of several research groups over the last few years who have applied modern femtosecond spectroscopies and new computational methods to these topics. Nitroaromatic molecules also have demonstrated their value as case-studies, where they can serve to understand the influence of torsional motions between the nitro substituent and the aromatic system in the conversions between states. In this contribution we highlight several of the recent results in this area. Due to the importance of the atmospheric photochemistry of nitrated compounds and their accumulating applications as nitric oxide release agents, continued research about the effects of the different state orderings, substitution patterns, and solvent effects is central to the development of future applications and for a better understanding of their environmental pathways. From this analysis, several pending issues are highlighted, which include the nature of the dominant singlet state involved in intersystem crossing, the role of the formation of charge-transfer states, the yield of the internal conversion channel to the electronic ground state, and a more generalized understanding of the sequence of steps which lead to nitric oxide dissociation.

Modeling the Formation, Degradation, and Spatiotemporal Distribution of 2-Nitrofluoranthene and 2-Nitropyrene in the Global Atmosphere

Polycyclic aromatic hydrocarbons (PAHs) are common atmospheric pollutants and known to cause adverse health effects. Nitrated PAHs (NPAHs) are formed in combustion activities and by nitration of PAHs in the atmosphere and may be equally or more toxic, but their spatial and temporal distribution in the atmosphere is not well characterized. Using the global EMAC model with atmospheric chemistry and surface compartments coupled, we investigate the formation, abundance, and fate of two secondarily formed NPAHs, 2-nitrofluoranthene (2-NFLT) and 2-nitropyrene (2-NPYR). The default reactivity scenario, the model with the simplest interpretation of parameters from the literature, tends to overestimate both absolute concentrations and NPAH/PAH ratios at observational sites. Sensitivity scenarios indicate that NO₂-dependent NPAH formation leads to better agreement between measured and predicted NPAH concentrations and that photodegradation is the most important loss process of 2-NFLT and 2-NPYR. The highest concentrations of 2-NFLT and 2-NPYR are found in regions with strong PAH emissions, but because of continued secondary formation from the PAH precursors, these two NPAHs are predicted to be spread across the globe.

Theoretical studies on the photochemistry of 2-nitrofluorene in the gas phase and acetonitrile solution

The photophysical and photochemical mechanisms of 2-nitrofluorene (2-NF) in the gas phase and acetonitrile solution have been studied theoretically. Upon ∼330 nm irradiation to the first bright state (¹ππ*), the 2-NF system can decay to triplet excited states via rapid intersystem crossing (ISC) processes through different surface crossing points or to the ground state via an ultrafast internal conversion (IC) process through the S₁/S₀ conical intersection. The ¹nπ* dark state will serve as a bridge when the system leaves the Franck-Condon (FC) region and approaches to the S₁ minimum. The molecule maintains a planar geometry during the excited-state relaxation processes. The differences on excitation properties such as electronic configurations and spin-orbit coupling (SOC) interactions between those in the gas phase and acetonitrile solution cannot be neglected, indicating possible changes on the efficiency of the related ISC processes for the 2-NF system in solution. Once arrived at the T₁ state, it would further decay to the S₀ state or photodegrade into the Ar-O˙ and NO˙ free radicals. During the intramolecular rearrangement process, the twisting of the nitro group out of the aromatic-ring plane is regarded as a critical structural variation for the photodegradation of the 2-NF system. The free radicals finally form through oxaziridine-type intermediate and transition state structures. The present work provides important mechanistic insights to the photochemistry of nitro-substituted polyaromatic compounds.

How important is roaming in the photodegradation of nitrobenzene?

Three pathways have been found for the formation of NO from nitrobenzene photodegradation that lead to either low or high translational energy, with a roaming mechanism involved at high excitation energies.

Nonadiabatic Dynamics Simulation Predict Intersystem Crossing in Nitroaromatic Molecules on a Picosecond Time Scale

Previous time-resolved spectroscopic experiments and static quantum-chemical calculations attributed nitronaphthalene derivatives one of the fastest time scales for intersystem crossing within organic molecules, reaching the 100 fs mark. Nonadiabatic dynamics simulations on three nitronaphthalene derivatives challenge this view, showing that the experimentally observed ∼100 fs process corresponds to internal conversion in the singlet manifolds. Intersystem crossing, instead, takes place on a longer time scale of ∼1 ps. The dynamics simulations further reveal that the spin transitions occur via two distinct pathways with different contribution for the three systems, which are determined by electronic factors and the torsion of the nitro group. This study, therefore, indicates that the existence of sub-picosecond intersystem crossing in other nitroaromatic molecules should be questioned.

The Kinetics and Mechanisms for Photodegradation of Nitrated Polycyclic Aromatic Hydrocarbons on Lettuce Leaf Surfaces: An In Vivo Study

Insights into the environmental fates of nitrated polycyclic aromatic hydrocarbons (NPAHs) in edible vegetables are of great significance for better evaluating human exposure to NPAHs through the dietary pathway. In this work, a fluorescence quenching method using graphene quantum dots as a fluorescent probe was first applied for the in vivo determination of 9-nitroanthracene (9-NAnt) and 1-nitropyrene (1-NPyr) adsorbed on the leaf surfaces of living lettuce (Lactuca sativa L.) seedlings. Moreover, the photolysis kinetics and mechanisms of the two adsorbed NPAHs were discussed. The photodegradation kinetics followed the pseudo-first-order equation, and the photodegradation half-life of 1-NPyr (7.4 ± 0.2 h) was greater than that of 9-NAnt (2.3 ± 0.1 h). Anthraquinone and pyrenediones were identified to be the main photolytic products of 9-NAnt and 1-NPyr, respectively. Intramolecular rearrangement was the most reasonable mechanism for the NPAH photolysis. The photolysis-driven degradation exhibited a key role in scavenging NPAHs from the vegetable leaf, indicating the reduction of NPAH transportation in the food chain.

Similar chemical structures, dissimilar triplet quantum yields: a CASPT2 model rationalizing the trend of triplet quantum yields in nitroaromatic systems

S₁/S₀ accessibility strongly influences the triplet quantum yields of nitronaphthalenes.

Mechanism of Ultrafast Intersystem Crossing in 2-Nitronaphthalene

Nitronaphthalene derivatives efficiently populate their electronically excited triplet states upon photoexcitation through ultrafast intersystem crossing (ISC). Despite having been studied extensively by time-resolved spectroscopy, the reasons behind their ultrafast ISC remain unknown. Herein, we present the first ab initio nonadiabatic molecular dynamics study of a nitronaphthalene derivative, 2-nitronaphthalene, including singlet and triplet states. We find that there are two distinct ISC reaction pathways involving different electronic states at distinct nuclear configurations. The high ISC efficiency is explained by the very small electronic and nuclear alterations that the chromophore needs to undergo during the singlet-triplet transition in the dominating ISC pathway after initial dynamics in the singlet manifold. The insights gained in this work are expected to shed new light on the photochemistry of other nitro polycyclic aromatic hydrocarbons that exhibit ultrafast intersystem crossing.

Ab Initio Study of Decay Dynamics of 1-Nitronaphthalene Initiated from the S2(ππ* + nNOπ*) State

Irradiation of nitro-PAHs in solution at ambient conditions leads to formation of its lowest excited triplet, dissociation intermediates nitrogen oxide (NO^•) and aryloxy radical (Ar-O^•). Experimental and theoretical studies demonstrated that Franck-Condon excited singlet state S_FC(ππ*) to a receiver, higher-energy triplet state T_n(nπ*) controlled the ultrafast population of the triplet state and, hence, the slight fluorescence yield of nitronaphthalenes. However, the detailed information about the curve-crossings of potential energy surfaces and the major channels for forming T₁ species and Ar-O^• radical were unclear. Here, by using the CASSCF//CASPT2 method, an efficient decay channel is revealed: S_2-FC-1NN → S_2-MIN-1NN or S₂T_3-MIN-1NN → T_3-MIN-1NN or T₃T_2-MIN-1NN→ T_2-MIN-1NN or T₂T_1-MIN-1NN → T_1-MIN-1NN. This explains the high yield of T_1-1NN species and minor yield of Ar-O^• and NO^• radicals. The calculation results suggest the bifurcation processes take place predominantly after the internal conversion to the T_1-1NN state via T₂T_1-MIN-1NN, one leads to T_1-MIN-1NN, while the other to T_1-MIN-ISO to produce Ar-O^• and NO^• radicals.

Intermolecular Hydrogen Abstraction from Hydroxy Group and Alkyl by T1(ππ*) of 1-Chloro-4-nitronaphthalene

Nanosecond transient absorption and theoretical calculations have been used to investigate the intermolecular hydrogen abstractions from alcohols and 1-naphthol by the lowest excited triplet (T₁) of 1-chloro-4-nitronaphthalene upon excitation of the compound in organic solvents. The hydrogen abstraction of T₁ from hydroxy group of 1-naphthol takes place through an electron transfer followed by a proton transfer through hydrogen bonding interaction with rate constants of ∼10⁹ M^-1 s^-1. Hydrogen-bonding is crucial in this process, indicated by the observation of a half reduction for T₁ yield when increasing the concentration of 1-naphthol. The hydrogen abstraction in this way can be decelerated by increasing solvent polarity and hydrogen-bonding donor ability. The T₁ of 1-chloro-4-nitronaphthalene can undergo one-step H atom abstraction from alkyl hydrogen in alcoholic solvents, with rate constants of ∼10⁴ M^-1 s^-1, and produce radical intermediates with the absorption maximum at 368 nm. DFT calculation results indicate both oxygens of the nitro group are active sites for hydrogen abstraction, and the difference of activation barriers for formation of two radical isomers is only 1.0 kcal/mol.

Photoinduced formation mechanism of the thymine–thymine (6–4) adduct in DNA; a QM(CASPT2//CASSCF):MM(AMBER) study

Out-of-plane puckering of the C_4′ atom forming a new CC bond seems essential for oxetane formation.

The IPEA dilemma in CASPT2

Multi-configurational second order perturbation theory (CASPT2) has become a very popular method for describing excited-state properties since its development in 1990. To account for systematic errors found in the calculation of dissociation energies, an empirical correction applied to the zeroth-order Hamiltonian, called the IPEA shift, was introduced in 2004. The errors were attributed to an unbalanced description of open-shell versus closed-shell electronic states and is believed to also lead to an underestimation of excitation energies. Here we show that the use of the IPEA shift is not justified and the IPEA should not be used to calculate excited states, at least for organic chromophores. This conclusion is the result of three extensive analyses. Firstly, we survey the literature for excitation energies of organic molecules that have been calculated with the unmodified CASPT2 method. We find that the excitation energies of 356 reference values are negligibly underestimated by 0.02 eV. This value is an order of magnitude smaller than the expected error based on the calculation of dissociation energies. Secondly, we perform benchmark full configuration interaction calculations on 137 states of 13 di- and triatomic molecules and compare the results with CASPT2. Also in this case, the excited states are underestimated by only 0.05 eV. Finally, we perform CASPT2 calculations with different IPEA shift values on 309 excited states of 28 organic small and medium-sized organic chromophores. We demonstrate that the size of the IPEA correction scales with the amount of dynamical correlation energy (and thus with the size of the system), and gets immoderate already for the molecules considered here, leading to an overestimation of the excitation energies. It is also found that the IPEA correction strongly depends on the size of the basis set. The dependency on both the size of the system and of the basis set, contradicts the idea of a universal IPEA shift which is able to compensate for systematic CASPT2 errors in the calculation of excited states.

Modeling the high-energy electronic state manifold of adenine: Calibration for nonlinear electronic spectroscopy

Pump-probe electronic spectroscopy using femtosecond laser pulses has evolved into a standard tool for tracking ultrafast excited state dynamics. Its two-dimensional (2D) counterpart is becoming an increasingly available and promising technique for resolving many of the limitations of pump-probe caused by spectral congestion. The ability to simulate pump-probe and 2D spectra from ab initio computations would allow one to link mechanistic observables like molecular motions and the making/breaking of chemical bonds to experimental observables like excited state lifetimes and quantum yields. From a theoretical standpoint, the characterization of the electronic transitions in the visible (Vis)/ultraviolet (UV), which are excited via the interaction of a molecular system with the incoming pump/probe pulses, translates into the determination of a computationally challenging number of excited states (going over 100) even for small/medium sized systems. A protocol is therefore required to evaluate the fluctuations of spectral properties like transition energies and dipole moments as a function of the computational parameters and to estimate the effect of these fluctuations on the transient spectral appearance. In the present contribution such a protocol is presented within the framework of complete and restricted active space self-consistent field theory and its second-order perturbation theory extensions. The electronic excited states of adenine have been carefully characterized through a previously presented computational recipe [Nenov et al., Comput. Theor. Chem. 1040-1041, 295-303 (2014)]. A wise reduction of the level of theory has then been performed in order to obtain a computationally less demanding approach that is still able to reproduce the characteristic features of the reference data. Foreseeing the potentiality of 2D electronic spectroscopy to track polynucleotide ground and excited state dynamics, and in particular its expected ability to provide conformational dependent fingerprints in dimeric systems, the performances of the selected reduced level of calculations have been tested in the construction of 2D electronic spectra for the in vacuo adenine monomer and the unstacked adenine homodimer, thereby exciting the Lb/La transitions with the pump pulse pair and probing in the Vis to near ultraviolet spectral window.