Deciphering the Internal Conversion Processes Involved in the Photochemical Ring-Opening of 1,3-Cyclohexadiene by Symmetric sp2-Carbon Substitutions

Chemical substituents hold the potential to markedly influence the photochemical behavior in molecular systems and assist in gaining a comprehensive understanding of nonadiabatic phenomena. In this study, we have conducted a comparative analysis of the influence of chemical substituents on the photochemical ring-opening of 1,3-cyclohexadiene (CHD), considering four systems: CHD, 2,3-dimethylcyclohexadiene (CHD-Me2-1), 1,4-dimethylcyclohexadiene (CHD-Me2-2), and 1,2,3,4-tetramethylcyclohexadiene (CHD-Me4), using electronic structure theory calculations and nonadiabatic molecular dynamics simulations. Employing extended multistate complete active space second-order perturbation (XMS-CASPT2) theory, we optimized reactants, S1 states, conical intersections (CIns), and products, revealing structural and energetic variations consistent with prior research. Nonadiabatic molecular dynamics simulation was used to gain insights into photochemical dynamics at state-averaged complete active space self-consistent field (SA-CASSCF) theory. CHD-Me4 exhibited reduced carbon-carbon single bond rupture rates, responsible for ring-opening, due to substituent proximity. Further, CHD-Me2-2 and CHD-Me4 displayed prolonged excited-state relaxation times, highlighting notable substituents' impact. Analysis of kinetic energy profiles of specific carbon atoms also revealed restrained atomic displacements, particularly in CHD-Me2-2 and CHD-Me4. These findings advance our understanding of how substituents modulate photochemical reactions in cyclohexadiene derivatives, guiding new molecular design and future research.

Classically forbidden nonadiabatic transitions in multidimensional chemical dynamics

An accurate method is proposed to deal with such nonadiabatic transitions as those energetically inaccessible, namely, classically forbidden transitions. This is formulated by using the corresponding Zhu-Nakamura formulas and finding the optimal paths in the classically forbidden tunneling regions that maximize the overall transition probabilities. This can be done for both the nonadiabatic tunneling type (so-called normal case in electron transfer) in which two diabatic potentials have opposite signs of slopes and the Landau-Zener type (inverted case) in which two diabatic potentials have the same sign of slopes. The method is numerically demonstrated to be useful for clarifying chemical and biological dynamics.

Dense and Acidic Organelle-Targeted Visualization in Living Cells: Application of Viscosity-Responsive Fluorescence Utilizing Restricted Access to Minimum Energy Conical Intersection

Cell-imaging methods with functional fluorescent probes are an indispensable technique to evaluate physical parameters in cellular microenvironments. In particular, molecular rotors, which take advantage of the twisted intramolecular charge transfer (TICT) process, have helped evaluate microviscosity. However, the involvement of charge-separated species in the fluorescence process potentially limits the quantitative evaluation of viscosity. Herein, we developed viscosity-responsive fluorescent probes for cell imaging that are not dependent on the TICT process. We synthesized AnP₂-H and AnP₂-OEG, both of which contain 9,10-di(piperazinyl)anthracene, based on 9,10-bis(N,N-dialkylamino)anthracene that adopts a nonflat geometry at minimum energy conical intersection. AnP₂-H and AnP₂-OEG exhibited enhanced fluorescence as the viscosity increased, with sensitivities comparable to those of conventional molecular rotors. In living cell systems, AnP₂-OEG showed low cytotoxicity and, reflecting its viscosity-responsive property, allowed specific visualization of dense and acidic organelles such as lysosomes, secretory granules, and melanosomes under washout-free conditions. These results provide a new direction for developing functional fluorescent probes targeting dense organelles.

Oxidation and Reduction Pathways in the Knowles Hydroamination via a Photoredox-Catalyzed Radical Reaction

Systematic reaction path exploration revealed the entire mechanism of Knowles's light-promoted catalytic intramolecular hydroamination. Bond formation/cleavage competes with single electron transfer (SET) between the catalyst and substrate. These processes are described by adiabatic processes through transition states in an electronic state and non-radiative transitions through the seam of crossings (SX) between different electronic states. This study determined the energetically favorable SET path by introducing a practical computational model representing SET as non-adiabatic transitions via SXs between substrate's potential energy surfaces for different charge states adjusted based on the catalyst's redox potential. Calculations showed that the reduction and proton shuttle process proceeded concertedly. Also, the relative importance of SET paths (giving the product and leading back to the reactant) varies depending on the catalyst's redox potential, affecting the yield.

Multi-state Energy Landscape for Photoreaction of Stilbene and Dimethyl-stilbene

We have recently developed the reaction space projector (ReSPer) method, which constructs a reduced-dimensionality reaction space uniquely determined from reference reaction paths for a polyatomic molecular system and projects classical trajectories into the same reaction space. In this paper, we extend ReSPer to the analysis of photoreaction dynamics and relaxation processes of stilbene and present the concept of a "multi-state energy landscape," incorporating the ground- and excited-state reaction subspaces. The multi-state energy landscape successfully explains the previously established photoreaction processes of cis-stilbene, such as the cis-trans photoisomerization and photocyclization. In addition, we discuss the difference in the excited-state reaction dynamics between stilbene and 1,1'-dimethyl stilbene based on a common reaction subspace determined from the framework part of reference structures with different number of atoms. This approach allows us to target any molecule with a common framework, greatly expanding the applicability of the ReSPer analysis. The multi-state energy landscape provides fruitful insight into photochemical reactions, exploring the excited- and ground-state potential energy surfaces, as well as comprehensive reaction processes with nonradiative transitions between adiabatic states, within the stage of a reduced-dimensionality reaction space.

Nonadiabatic Dynamics Algorithms with Only Potential Energies and Gradients: Curvature-Driven Coherent Switching with Decay of Mixing and Curvature-Driven Trajectory Surface Hopping

Direct dynamics by mixed quantum-classical nonadiabatic methods is an important tool for understanding processes involving multiple electronic states. Very often, the computational bottleneck of such direct simulation comes from electronic structure theory. For example, at every time step of a trajectory, nonadiabatic dynamics requires potential energy surfaces, their gradients, and the matrix elements coupling the surfaces. The need for the couplings can be alleviated by employing the time derivatives of the wave functions, which can be evaluated from overlaps of electronic wave functions at successive time steps. However, evaluation of overlap integrals is still expensive for large systems. In addition, for electronic structure methods for which the wave functions or the coupling matrix elements are not available, nonadiabatic dynamics algorithms become inapplicable. In this work, building on recent work by Baeck and An, we propose new nonadiabatic dynamics algorithms that only require adiabatic potential energies and their gradients. The new methods are named curvature-driven coherent switching with decay of mixing (κCSDM) and curvature-driven trajectory surface hopping (κTSH). We show how powerful these new methods are in terms of computation time and accuracy as compared to previous mixed quantum-classical nonadiabatic dynamics algorithms. The lowering of the computational cost will allow longer nonadiabatic trajectories and greater ensemble averaging to be affordable, and the ability to calculate the dynamics without electronic structure coupling matrix elements extends the dynamics capability to new classes of electronic structure methods.

A systematic model study quantifying how conical intersection topography modulates photochemical reactions

Despite their important role in photochemistry and expected presence in most polyatomic molecules, conical intersections have been thoroughly characterized in a comparatively small number of systems. Conical intersections can confer molecular photoreactivity or photostability, often with remarkable efficacy, due to their unique structure: at a conical intersection, the adiabatic potential energy surfaces of two or more electronic states are degenerate, enabling ultrafast decay from an excited state without radiative emission, known as nonadiabatic transfer. Furthermore, the precise conical intersection topography determines fundamental properties of photochemical processes, including excited-state decay rate, efficacy, and molecular products that are formed. However, these relationships have yet to be defined comprehensively. In this article, we use an adaptable computational model to investigate a variety of conical intersection topographies, simulate resulting nonadiabatic dynamics, and calculate key photochemical observables. We varied the vibrational mode frequencies to modify conical intersection topography systematically in four primary classes of conical intersections and quantified the resulting rate, total yield, and product yield of nonadiabatic decay. The results reveal that higher vibrational mode frequencies reduce nonadiabatic transfer, but increase the transfer rate and resulting photoproduct formation. These trends can inform progress toward experimental control of photochemical reactions or tuning of molecules' photochemical properties based on conical intersections and their topography.

Anomalous Linear Dichroism in Bent Chromophores of π-conjugated Polymers: Departure from the Franck-Condon Principle

We examine the influence of bending of π-conjugated chromophores on photoluminescence (PL) by spectrally resolving the depolarization of fluorescence on the single-molecule level. The effect of excited-state mixing mediated by molecular vibrations is manifested in the departure from the usual achromatic linear dichroism of fluorescence, with the polarization anisotropy decreasing in the vibronic progression. Bent chromophores reveal an overall increase in vibronic PL intensity with polarization orthogonal to the molecular long axis. This manifestation of the Renner-Herzberg-Teller (RHT) effect illustrates the breakdown of the Franck-Condon principle in macromolecules used in organic electronics, providing information on the orientation of transition-dipole moments and the origin of spectral broadening. While some of the spectral signatures of the RHT effect appear similar to those of H aggregation in molecular dimers, discrimination between the two phenomena is straightforward since H aggregation does not induce anomalous linear dichroism.

An excited-state Wolff rearrangement reaction of 5-diazo Meldrum's acid: an ab initio on-the-fly nonadiabatic dynamics simulation

A global switching on-the-fly trajectory surface hopping dynamics simulation at the 3SA-CASSCF(12,11)/6-31G* quantum level has been employed to probe the photo-induced Wolff rearrangement (WR) reaction of 5-diazo Meldrum's acid (DMA) within three low-lying electronic excited states. The present simulation predicted that the branching ratios for relaxing back to the ground state, isomerizing to diazirine, and reaction to ketene I via carbene I are 69% ± 0.1, 3% ± 0.4, and 28% ± 0.1, which are in excellent agreement with those obtained by the femtosecond spectroscopy experiment, 67%, 3% and 30%, respectively. In particular, the present simulation revealed that the major WR reaction to ketene I pathway is stepwise via the excited-state to carbene I (17.8% ± 0.2) and via the ground-state to carbene I (8.7% ± 0.2), and the minor pathway is concerted synchronous (1.5% ± 0.6). The photo-induced WR reaction of DMA has been quantitatively interpreted in terms of the distribution of extended seam surfaces as a function of CN dissociation bonds for two important conical intersections within three low-lying electronic excited states. Ultrafast dynamic time constants have been estimated to be about 500 fs ± 120 fs and 180 fs ± 80 fs for the stepwise and the concerted WR reaction to ketene I which are also in good agreement with those determined by the experiment. Therefore, the photo-induced excited-state WR reaction mechanism has been quantitatively revealed by the present real-time dynamics simulation.

The photoinduced isomerization mechanism of the 2-(1-(methylimino)methyl)-6-chlorophenol (SMAC): Nonadiabatic surface hopping dynamics simulations

The photophysical properties of the Schiff base family are crucial for their applications such as molecular switches and molecular memories. However, it was found that the photophysical behavior is not uniform for all Schiff base molecules, which shows a significant substituent dependent property. In this article, we studied the photoisomerization mechanism of one Schiff base chlorosubstituted derivative 2-(1-(methylimino)methyl)-6-chlorophenol by employing geometrical optimization, energy profiles scanning, and on-the-fly dynamical simulations. Three types of minimum energy conical intersections were located on the S₁/S₀ crossing seam, with two characterized by twisting motion of the C=N bond and one featured with the excited state intramolecular proton transfer process and then twisting motion around the C=C bond [excited-state intramolecular proton transfer process (ESIPT)-then-twisting]. By a combination of the dynamics simulation results with the energy profiles scanned along with the ESIPT coordinate, it was found that the photophysical property of the targeted molecule is different from that of most Schiff base members, which prefer to decay by a twisting motion around the C=N bridge bond rather than the ESIPT-then-twisting channel. The minor ESIPT channel is probably governed by a tunneling mechanism. The proposed deactivation mechanism can provide a reasonable explanation for the observations in the experiment and would provide fundamental indications for further design of new and efficient photochromic products.

Exploring radiative and nonradiative decay paths in indole, isoindole, quinoline, and isoquinoline

Radiative and nonradiative decay paths from the first excited singlet electronic state (S₁) in four heteroaromatics, indole, isoindole, quinoline, and isoquinoline, were systematically explored.

Ultrafast intersystem crossing for nitrophenols: ab initio nonadiabatic molecular dynamics simulation

Ultrafast intersystem crossing mechanisms for two p- and m-nitrophenol groups (PNP and MNP) have been investigated using ab initio nonadiabatic molecular dynamics simulations at the 6SA-CASSCF level of theory.

Intersystem crossing-branched excited-state intramolecular proton transfer for o-nitrophenol: An ab initio on-the-fly nonadiabatic molecular dynamic simulation

The 6SA-CASSCF(10, 10)/6-31G (d, p) quantum chemistry method has been applied to perform on-the-fly trajectory surface hopping simulation with global switching algorithm and to explore excited-state intramolecular proton transfer reactions for the o-nitrophenol molecule within low-lying electronic singlet states (S0 and S1) and triplet states (T1 and T2). The decisive photoisomerization mechanisms of o-nitrophenol upon S1 excitation are found by three intersystem crossings and one conical intersection between two triplet states, in which T1 state plays an essential role. The present simulation shows branch ratios and timescales of three key processes via T1 state, non-hydrogen transfer with ratio 48% and timescale 300 fs, the tunneling hydrogen transfer with ratios 36% and timescale 10 ps, and the direct hydrogen transfer with ratios 13% and timescale 40 fs. The present simulated timescales might be close to low limit of the recent experiment results.

Control of Nonadiabatic Passage through a Conical Intersection by a Dynamic Resonance

Nonadiabatic processes, dominated by dynamic passage of reactive fluxes through conical intersections (CIs), are considered to be appealing means for manipulating reaction paths, particularly via initial vibrational preparation. Nevertheless, obtaining direct experimental evidence of whether specific-mode excitation affects the passage at the CI is challenging, requiring well-resolved time- or frequency-domain experiments. Here promotion of methylamine-d2 (CH3ND2) molecules to spectral-resolved rovibronic states on the excited S1 potential energy surface, coupled to sensitive D photofragment probing, allowed us to follow the N-D bond fission dynamics. The branching ratios between slow and fast D photofragments and the internal energies of the CH3ND(X̃) photofragments confirm correlated anomalies for predissociation initiated from specific rovibronic states. These anomalies reflect the existence of a dynamic resonance that strongly depends on the energy of the initially excited rovibronic states, the evolving vibrational mode on the repulsive S1 part during N-D bond elongation, and the manipulated passage through the CI that leads to CH3ND radicals excited with C-N-D bending. This resonance plays an important role in the bifurcation dynamics at the CI and can be foreseen to exist in other photoinitiated processes and to control their outcome.

Development of semiclassical molecular dynamics simulation method

Various quantum mechanical effects such as nonadiabatic transitions, quantum mechanical tunneling and coherence play crucial roles in a variety of chemical and biological systems.

Exploration of Quenching Pathways of Multiluminescent Acenes Using the GRRM Method with the SF-TDDFT Method

The quenching pathways were investigated for three types of multiluminescent acene derivatives, which show environment-dependent fluorescence. Spin-flip time dependent density functional theory (SF-TDDFT) combined with the Global Reaction Route mapping (GRRM) strategy is employed to locate minimum-energy conical intersections (MECIs). The energies and geometries of the MECIs relative to the Franck-Condon (FC) state control the difference in fluorescence behavior among the three derivatives. For the molecule with a phenyamide moiety, a MECI with energy lower than the FC state with large geometrical change from V-type to flat structure provides an efficient internal conversion (quenching) pathway in solution. For the same molecule, in a solid, this large geometrical change is inhibited, and the second MECI, with an energy lower than FC but higher than the first MECI requiring only a small geometry change of CH out-of-plane bending, contributes to the quenching. The molecule with the napthaleneimide moiety has only one low-energy MECI that requires large geometrical change from the V-type to flat structure. Although this MECI provides the quenching pathway in solution, in the solid, this large motion is inhibited, and the molecule will stay in the excited state and emit. The molecule with an anthraceneimide moiety has no conical intersection lower than the FC state, and no quenching pathway is available in solution or solid. In addition, in this molecule, at the local minimum of the excited state, the dipole transition to the ground state is allowed, and this molecule prefers emission rather than internal conversion.

Tropospheric photolysis rates of the acetaldehyde isotopologues CD₃CHO and CD₃CDO relative to CH₃CHO measured at the European Photoreactor Facility

Acetaldehyde is a hazardous pollutant found in indoor and ambient air. Acetaldehyde photolysis is pressure- and wavelength-dependent with three distinct product channels. In this study, the photolysis rates of CH3CHO, CD3CDO, and CD3CHO are studied in natural tropospheric conditions using long path FTIR spectroscopy, at the European Photoreactor Facility (EUPHORE) in Valencia, Spain. The average relative photolysis rate as an average of four experiments for the fully deuterated isotopologue is j(CH3CHO)/j(CD3CDO) = 1.75 ± 0.04, and as a result of a single experiment j(CH3CHO)/j(CD3CHO) = 1.10 ± 0.10. These results, combined with our previous determination of j(CH3CHO)/j(CH3CDO) = 1.26 ± 0.03, provide mechanistic insight into the photodissociation dynamics of the photoexcited species. Despite the extensive isotopic scrambling in photoexcited acetaldehyde that has recently been reported, the position of the substitution has a clear effect on the relative photolysis rates.

Evidence for quantum effects in the predissociation of methylamine isotopologues

The H product distributions obtained from the predissociation of methylamine isotopologues are extremely sensitive to the energy difference between the initially prepared vibrational states and the conical intersections and not only to the nature of the pre-excited nuclear motions.

Exploration of minimum energy conical intersection structures of small polycyclic aromatic hydrocarbons: toward an understanding of the size dependence of fluorescence quantum yields

A correlation between the fluorescence quantum yields and the energy barrier to the conical intersection region was discovered for five small PAHs.

Exploring Multiple Potential Energy Surfaces: Photochemistry of Small Carbonyl Compounds

In theoretical studies of chemical reactions involving multiple potential energy surfaces (PESs) such as photochemical reactions, seams of intersection among the PESs often complicate the analysis. In this paper, we review our recipe for exploring multiple PESs by using an automated reaction path search method which has previously been applied to single PESs. Although any such methods for single PESs can be employed in the recipe, the global reaction route mapping (GRRM) method was employed in this study. By combining GRRM with the proposed recipe, all critical regions, that is, transition states, conical intersections, intersection seams, and local minima, associated with multiple PESs, can be explored automatically. As illustrative examples, applications to photochemistry of formaldehyde and acetone are described. In these examples as well as in recent applications to other systems, the present approach led to discovery of many unexpected nonadiabatic pathways, by which some complicated experimental data have been explained very clearly.

Slice imaging and wave packet study of the photodissociation of CH3I in the blue edge of the A-band: evidence of reverse 3Q0 ← 1Q1 non-adiabatic dynamics

The photodissociation of CH(3)I in the blue edge (217-230 nm) of the A-band has been studied using a combination of slice imaging and resonance enhanced multiphoton ionization (REMPI) detection of the CH(3) fragment in the vibrational ground state (ν = 0). The profiles of the CH(3) (ν = 0) kinetic energy distributions and the photofragment anisotropies are interpreted in terms of the contribution of the excited surfaces involved in the photodissociation process, as well as the probability of non-adiabatic curve crossing between the (3)Q(0) and (1)Q(1) states. In the studied region, unlike in the central part of the A-band where absorption to the (3)Q(0) state dominates, the I((2)P(J)), with J = 1/2, 3/2, in correlation with CH(3) (ν = 0) kinetic energy distributions show clearly two contributions of different anisotropy, signature of the competing adiabatic and non-adiabatic dynamics, whose ratio strongly depends on the photolysis wavelength. The experimental results are compared with multisurface wave packet calculations carried out using the available ab initio potential energy surfaces, transition moments, and non-adiabatic couplings, employing a reduced dimensionality model. A good qualitative agreement is found between experiment and theory and both show evidence of reverse (3)Q(0)←(1)Q(1) non-adiabatic dynamics at the bluest excitation wavelengths both in the fragment kinetic energy and angular distributions.

Finding Minimum Structures on the Seam of Crossing in Reactions of Type A + B → X: Exploration of Nonadiabatic Ignition Pathways of Unsaturated Hydrocarbons

A new theoretical approach is proposed for finding automatically minimum structures on the seam of crossing (MSX) in reactions of type A + B → X, where the artificial-force-induced reaction (AFIR) method is combined with the seam model function (SMF) approach. Its application to reactions between triplet dioxygen and unsaturated hydrocarbons provided many MSX structures. In addition to known ignition pathways, we discovered a pathway through a new type of MSX in the reaction of dioxygen with aromatic hydrocarbons; for benzene, this new pathway requires a lower energy than those of three known ignition pathways and is likely to be the most important. This demonstrates that the AFIR-SMF approach has the ability to discover unknown/unexpected MSX structures without prejudice for presumed pathways or mechanisms.