skip to main content


Title: Automatic discovery of photoisomerization mechanisms with nanosecond machine learning photodynamics simulations
Photochemical reactions are widely used by academic and industrial researchers to construct complex molecular architectures via mechanisms that often require harsh reaction conditions. Photodynamics simulations provide time-resolved snapshots of molecular excited-state structures required to understand and predict reactivities and chemoselectivities. Molecular excited-states are often nearly degenerate and require computationally intensive multiconfigurational quantum mechanical methods, especially at conical intersections. Non-adiabatic molecular dynamics require thousands of these computations per trajectory, which limits simulations to ∼1 picosecond for most organic photochemical reactions. Westermayr et al. recently introduced a neural-network-based method to accelerate the predictions of electronic properties and pushed the simulation limit to 1 ns for the model system, methylenimmonium cation (CH 2 NH 2 + ). We have adapted this methodology to develop the Python-based, Python Rapid Artificial Intelligence Ab Initio Molecular Dynamics (PyRAI 2 MD) software for the cis – trans isomerization of trans -hexafluoro-2-butene and the 4π-electrocyclic ring-closing of a norbornyl hexacyclodiene. We performed a 10 ns simulation for trans -hexafluoro-2-butene in just 2 days. The same simulation would take approximately 58 years with traditional multiconfigurational photodynamics simulations. We generated training data by combining Wigner sampling, geometrical interpolations, and short-time quantum chemical trajectories to adaptively sample sparse data regions along reaction coordinates. The final data set of the cis – trans isomerization and the 4π-electrocyclic ring-closing model has 6207 and 6267 data points, respectively. The training errors in energy using feedforward neural networks achieved chemical accuracy (0.023–0.032 eV). The neural network photodynamics simulations of trans -hexafluoro-2-butene agree with the quantum chemical calculations showing the formation of the cis -product and reactive carbene intermediate. The neural network trajectories of the norbornyl cyclohexadiene corroborate the low-yielding syn -product, which was absent in the quantum chemical trajectories, and revealed subsequent thermal reactions in 1 ns.  more » « less
Award ID(s):
1940307
NSF-PAR ID:
10297772
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Chemical Science
Volume:
12
Issue:
14
ISSN:
2041-6520
Page Range / eLocation ID:
5302 to 5314
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Benzene fluorination increases chemoselectivities for Dewar‐benzenes via 4π‐disrotatory electrocyclization. However, the origin of the chemo‐ and regioselectivities of fluorobenzenes remains unexplained because of the experimental limitations in resolving the excited‐state structures on ultrafast timescales. The computational cost of multiconfigurational nonadiabatic molecular dynamics simulations is also currently cost‐prohibitive. We now provide high‐fidelity structural information and reaction outcome predictions with machine‐learning‐accelerated photodynamics simulations of a series of fluorobenzenes, C6F6‐nHn, n=0–3, to study their S1→S0decay in 4 ns. We trained neural networks with XMS‐CASPT2(6,7)/aug‐cc‐pVDZ calculations, which reproduced the S1absorption features with mean absolute errors of 0.04 eV (<2 nm). The predicted nonradiative decay constants for C6F4H2, C6F6, C6F3H3, and C6F5H are 116, 60, 28, and 12 ps, respectively, in broad qualitative agreement with the experiments. Our calculations show that a pseudo Jahn–Teller distortion of fluorinated benzenes leads to an S1local‐minimum region that extends the excited‐state lifetimes of fluorobenzenes. The pseudo Jahn–Teller distortions reduce when fluorination decreases. Our analysis of the S1dynamics shows that the pseudo‐Jahn–Teller distortions promote an excited‐statecis‐transisomerization of a πC‐Cbond. We characterized the surface hopping points from our NAMD simulations and identified instantaneous nuclear momentum as a factor that promotes the electrocyclizations.

     
    more » « less
  2. Bis(bithienyl)-1,2-dicyanoethene (4TCE) is a photoswitch that operates via reversible E / Z photoisomerization following absorption of visible light. cis -to- trans photoisomerization of 4TCE requires excitation below 470 nm, is relatively inefficient (quantum yield < 5%) and occurs via the lowest-lying triplet. We present excitation-wavelength dependent (565–420 nm) transient absorption (TA) studies to probe the photophysics of cis -to- trans isomerization to identify sources of switching inefficiency. TA data reveals contributions from more than one switch conformer and relaxation cascades between multiple states. Fast (∼4 ps) and slow (∼40 ps) components of spectral dynamics observed at low excitation energies (>470 nm) are readily attributed to deactivation of two conformers; this assignment is supported by computed thermal populations and absorption strengths of two molecular geometries (P A and P B ) characterized by roughly parallel dipoles for the thiophenes on opposite sides of the ethene bond. Only the P B conformer is found to contribute to triplet population and the switching of cis -4TCE: high-energy excitation (<470 nm) of P B involves direct excitation to S 2 , relaxation from which prepares an ISC-active S 1 geometry (ISC QY 0.4–0.67, k ISC ∼ 1.6–2.6 × 10 −9 s −1 ) that is the gateway to triplet population and isomerization. We ascribe low cis -to- trans isomerization yield to excitation of the nonreactive P A conformer (75–85% loss) as well as loses along the P B S 2 → S 1 → T 1 cascade (10–20% loss). In contrast, electrocyclization is inhibited by the electronic character of the excited states, as well as a non-existent thermal population of a reactive “antiparallel” ring conformation. 
    more » « less
  3. We report here a study that has revealed two distinct modes of reactivity of azobenzene derivatives (ArNNAr) with benzynes, depending on whether the aryne reacts with a trans - or a cis -azobenzene geometric isomer. Under thermal conditions, trans -azobenzenes engage benzyne via an initial [2 + 2] trapping event, a process analogous to known reactions of benzynes with diarylimines (ArCNAr). This is followed by an electrocyclic ring opening/closing sequence to furnish dihydrophenazine derivatives, subjects of contemporary interest in other fields ( e.g. , electronic and photonic materials). In contrast, when the benzyne is attacked by a cis -azobenzene, formation of aminocarbazole derivatives occurs via an alternative, net (3 + 2) pathway. We have explored these complementary orthogonal processes both experimentally and computationally. 
    more » « less
  4. We demonstrate how two-photon excitation with quantum light can influence elementary photochemical events. The azobenzene trans → cis isomerization following entangled two-photon excitation is simulated using quantum nuclear wave packet dynamics. Photon entanglement modulates the nuclear wave packets by coherently controlling the transition pathways. The photochemical transition state during passage of the reactive conical intersection in azobenzene photoisomerization is strongly affected with a noticeable alteration of the product yield. Quantum entanglement thus provides a novel control knob for photochemical reactions. The distribution of the vibronic coherences during the conical intersection passage strongly depends on the shape of the initial wave packet created upon quantum light excitation. X-ray signals that can experimentally monitor this coherence are simulated. 
    more » « less
  5. A refractive index of dielectrics was modified by several methods and was known to have direct influence on optical forces in nanophotonic structures. The present contribution shows that isomerization of photoswitching molecules can be used to regulate refractive index of dielectrics in-situ. In particular, spectroscopic study of a polydimethylsiloxane–arylazopyrazole (PDMS–AAP) composite revealed that refractive index of the composite shifts from 2.0 to 1.65 in trans and cis states, respectively, of the embedded AAP. Based on this, a proposition is made for a waveguide structure, in which external UV/Vis source reversibly regulates the conformation of the PDMS–AAP core. Computational study is performed using Maxwell’s equations on buried waveguide structure. The simulation, implemented in PYTHON, sequentially utilizes empirical refractive indices of the composite in the isomeric states in lieu of regulation by a source. The simulation revealed highly confined wave propagations for injected signals of 340 and 450 nm wavelengths. It is observed that the cis state suppresses higher order mode when propagating UV wavelength but allows it for visible light. This modal tuning demonstrated that single mode can be selectively excited with appropriate waveguide dimensions. Further impact of the tuning is seen in the optical force between waveguide pair where the forces shift between attractive and repulsive in relation to the isomeric state of the PDMS–AAP core. These effects which stem from the adjustment of refractive index by photoisomerization suggests that in-situ regulation of index is achievable by successful integration of photoswitching molecules in host materials, and the current PDMS–AAP composites investigated in this study can potentially enhance nanophotonic and opto-mechanical platforms. 
    more » « less