Geometric phase effects in the coherent control of the branching ratio of photodissociation products of phenol

Probing avoided crossings and conical intersections by two-pulse femtosecond stimulated Raman spectroscopy: Theoretical study

This study leverages two-pulse femtosecond stimulated Raman spectroscopy (2FSRS) to characterize molecular systems with avoided crossings (ACs) and conical intersections (CIs) in their low-lying excited electronic states. By simulating 2FSRS spectra of microscopically inspired ACs and CIs models, we demonstrate that 2FSRS not only delivers valuable information on the molecular parameters characterizing ACs and CIs but also helps distinguish between these two systems.

On the multiphoton ionisation photoelectron spectra of phenol

The phenol molecule is a prototype for non-adiabatic dynamics and the excited-state photochemistry of biomolecules. In this article, we report a joint theoretical and experimental investigation on the resonance enhanced multiphoton ionisation photoelectron (REMPI) spectra of the two lowest ionisation bands of phenol. The focus is on the theoretical interpretation of the measured spectra using quantum dynamics simulations. These were performed by numerically solving the time-dependent Schrödinger equation using the multi-layer variant of the multiconfiguration time-dependent Hartree algorithm together with a vibronic coupling Hamiltonian model. The ionising laser pulse is modelled explicitly within the ionisation continuum model to simulate experimental femtosecond 1+1 REMPI photoelectron spectra. These measured spectra are sensitive to very short lived electronically excited states, providing a rigorous benchmark for our theoretical methods. The match between experiment and theory allows for an interpretation of the features of the spectra at different wavelengths and shows that there are features due to both 'direct' and 'indirect' ionisation, resulting from non-resonant and resonant excitation by the pump pulse.

Influence of Mode-Specific Excitation on the Nonadiabatic Dynamics of Methyl Nitrate (CH3ONO2)

The impact of mode-specific vibrational excitations on initial-preparation conditions was studied by examining the excited-state population decay rates in the nonadiabatic dynamics of methyl nitrate (CH₃ONO₂). In particular, exciting a few specific modes by adding a single quantum of energy clearly decelerated the nonadiabatic dynamics population decay rates. The underlying reason for this slower population decay was explained by analyzing the profiles of the excited-state potential energy surfaces in the Franck-Condon regions and the topology of the S₁/S₀ conical intersection. This study not only provides physical insights into the key mechanisms controlling nonadiabatic dynamics but also shows the possibility of controlling nonadiabatic dynamics via mode-specific vibrational excitations.

Mode-dependent H atom tunneling dynamics of the S1 phenol is resolved by the simple topographic view of the potential energy surfaces along the conical intersection seam

Mode-dependent H atom tunneling dynamics of the O-H bond predissociation of the S₁ phenol has been theoretically analyzed. As the tunneling is governed by the complicated multi-dimensional potential energy surfaces that are dynamically shaped by the upper-lying S₁(ππ*)/S₂(πσ*) conical intersection, the mode-specific tunneling dynamics of phenol (S₁) has been quite formidable to be understood. Herein, we have examined the topography of the potential energy surface along the particular S₁ vibrational mode of interest at the nuclear configurations of the S₁ minimum and S₁/S₂ conical intersection. The effective adiabatic tunneling barrier experienced by the reactive flux at the particular S₁ vibrational mode excitation is then uniquely determined by the topographic shape of the potential energy surface extended along the conical intersection seam coordinate associated with the particular vibrational mode. The resultant multi-dimensional coupling of the specific vibrational mode to the tunneling coordinate is then reflected in the mode-dependent tunneling rate as well as nonadiabatic transition probability. Remarkably, the mode-specific experimental result of the S₁ phenol tunneling reaction [K. C. Woo and S. K. Kim, J. Phys. Chem. A 123, 1529-1537 (2019)] (in terms of the qualitative and relative mode-dependent dynamic behavior) could be well rationalized by semi-classical calculations based on the mode-specific topography of the effective tunneling barrier, providing the clear conceptual insight that the skewed potential energy surfaces along the conical intersection seam (strongly or weakly coupled to the tunneling reaction coordinate) may dictate the tunneling dynamics in the proximity of the conical intersection.

Optimal control of photodissociation of phenol using genetic algorithm

Photodissociation dynamics of the OH bond of phenol is studied with an optimally shaped laser pulse. The theoretical model consists of three electronic states (the ground electronic state, ππ* state, and πσ* state) in two nuclear coordinates (the OH stretching coordinate as a reaction coordinate, r, and the CCOH dihedral angle as a coupling coordinate, θ). The optimal UV laser pulse is designed using the genetic algorithm, which optimizes the total dissociative flux of the wave packet. The latter is calculated in the adiabatic asymptotes of the S0 and S1 electronic states of phenol. The initial state corresponds to the vibrational levels of the electronic ground state and is defined as | n r, n θ⟩, where n r and n θ represent the number of nodes along r and θ, respectively. The optimal UV field excites the system to the optically dark πσ* state predominantly over the optically bright ππ* state with the intensity borrowing effect for the |0, 0⟩ and |0, 1⟩ initial states. For the |0, 0⟩ initial condition, the photodissociation to the S1 asymptotic channel is favored slightly over the S0 asymptotic channel. Addition of one quantum of energy along the coupling coordinate increases the dissociation probability in the S1 channel. This is because the wave packet spreads along the coupling coordinate on the πσ* state and follows the adiabatic path. Hence, the S1 asymptotic channel gets more ([Formula: see text]11%) dissociative flux as compared to the S0 asymptotic channel for the |0, 1⟩ initial condition. The |1, 0⟩ and |1, 1⟩ states are initially excited to both the ππ* and πσ* states in the presence of the optimal UV pulse. For these initial conditions, the S1 channel gets more dissociative flux as compared to the S0 channel. This is because the high energy components of the wave packet readily reach the S1 channel. The central frequency of the optimal UV pulse for the |0, 0⟩ and |0, 1⟩ initial states has a higher value as compared to the |1, 0⟩ and |1, 1⟩ initial states. This is explained with the help of an excitation mechanism of a given initial state in relation to its energy.

Deprotonation of Phenol linked to a silicon dioxide surface using Adaptive Feedback Laser Control with a Heterodyne Detected Sum Frequency Generation Signal

The development of laser-controlled surface reactions has been limited by the lack of decisive methods for detecting evolving changes in the surface chemistry. In this work, we demonstrate successful laser...

Conical intersection seams in spin-orbit coupled systems with an even number of electrons: A numerical study based on neural network fit surfaces

In this work, we consider the existence and topography of seams of conical intersections (CIs) for two key singlet-triplet systems, including a uniformly scaled spin-orbit interaction. The basic one triplet and one singlet state system denoted as (S₀,T₁) and the two singlets and one triplet system denoted as (S₀,S₁,T₁) are treated. Essential to this analysis are realistic electronic structure data taken from a recently reported neural network fit for the 1,2¹A and 1³A states of NH₃, including H_sf (spin-free) and H_so (spin-orbit) surfaces derived from high quality ab initio wavefunctions. Three types of seams for the (S₀,S₁,T₁) system are reported, which depend on the choice of the electronic Hamiltonian, H^e. The nonrelativistic CI seam [H^e = H_sf, (S₀,S₁)], the energy minimized nonrelativistic singlet-triplet intersection seam [H^e = H_sf, (S₀,T₁)], and the fully relativistic seam in the spin-diabatic representation (H^e = H_tot = H_sf + H_so) are reported as functions of R(N-H). The derivative couplings are computed using H^e = H_tot and H_sf from the fit data. The line integral of the derivative coupling is employed to juxtapose the geometric phase in the relativistic, H^e = H_tot, and nonrelativistic, H^e = H_sf, cases. It is found for the (S₀,T₁) system that there is no CI in the spin-adiabatic representation, while for the (S₀,S₁,T₁) system, CI can only be formed for two pairs of spin-adiabatic electronic states. The geometric phase effect thus needs to be handled with care when it comes to spin-nonconserving dynamics simulations.

Femtosecond Wavepacket Dynamics Reveals the Molecular Structures in the Excited (S1) and Cationic (D0) States

Molecular structures in the electronically excited (S₁) and cationic (D₀) states of 2-fluorothioanisole (2-FTA) have been precisely refined from the real-time dynamics of the femtosecond (fs) wavepacket prepared by the coherent excitation of the Franck-Condon active out-of-plane torsional modes in the S₁ ← S₀ transition at 285 nm. The simulation to reproduce the experiment in terms of the beating frequencies gives the nonplanar geometry of 2-FTA in S₁, where the out-of-plane dihedral angle (φ) of the S-CH₃ moiety is 51° with respect to the molecular plane. The behavior of the fs wavepacket in terms of the amplitudes and phases with the change of the probe (ionization) wavelength (λ_probe = 300-330 nm) provides the otherwise veiled structure of the cationic D₀ state. While the 2-FTA cation adopts the planar geometry (φ = 0°) at the global minimum, it is found to have a vertical minimum at φ ≈ 135° from the perspective of the D₀ ← S₁ vertical transition. Ab initio calculations support the experiment quite well although the simulation using the model potentials could improve the match with the experiment, giving the new interpretation for the previously disputed photoelectron spectroscopic results.

Coupled nuclear and electron dynamics in the vicinity of a conical intersection

Ultrafast optical techniques allow us to study ultrafast molecular dynamics involving both nuclear and electronic motion. To support interpretation, theoretical approaches are needed that can describe both the nuclear and electron dynamics. Hence, we revisit and expand our ansatz for the coupled description of the nuclear and electron dynamics in molecular systems (NEMol). In this purely quantum mechanical ansatz, the quantum-dynamical description of the nuclear motion is combined with the calculation of the electron dynamics in the eigenfunction basis. The NEMol ansatz is applied to simulate the coupled dynamics of the molecule NO₂ in the vicinity of a conical intersection (CoIn) with a special focus on the coherent electron dynamics induced by the non-adiabatic coupling. Furthermore, we aim to control the dynamics of the system when passing the CoIn. The control scheme relies on the carrier envelope phase of a few-cycle IR pulse. The laser pulse influences both the movement of the nuclei and the electrons during the population transfer through the CoIn.

Geometric phase effects on photodissociation dynamics of diatomics

We investigated the effect of the geometric phase (GP) on photodissociation dynamics at a light-induced conical intersection (LICI) through exact quantum dynamical calculations. By taking the one-photon photodissociation of H 2 + ionic molecules as an example, we explored the conditions wherein the LICI associated GP affects dissociation dynamics. We found that GP leads to a phase shift between the angular distributions of GP included and GP excluded photofragments. This effect is more pronounced when the energy of the initial vibrational level is above the energy of the LICI point.

Excited-state dissociation dynamics of phenol studied by a new time-resolved technique

Phenol is an important model molecule for the theoretical and experimental investigation of dissociation in the multistate potential energy surfaces. Recent theoretical calculations [X. Xu et al., J. Am. Chem. Soc. 136, 16378 (2014)] suggest that the phenoxyl radical produced in both the X and A states from the O-H bond fission in phenol can contribute substantially to the slow component of photofragment translational energy distribution. However, current experimental techniques struggle to separate the contributions from different dissociation pathways. A new type of time-resolved pump-probe experiment is described that enables the selection of the products generated from a specific time window after molecules are excited by a pump laser pulse and can quantitatively characterize the translational energy distribution and branching ratio of each dissociation pathway. This method modifies conventional photofragment translational spectroscopy by reducing the acceptance angles of the detection region and changing the interaction region of the pump laser beam and the molecular beam along the molecular beam axis. The translational energy distributions and branching ratios of the phenoxyl radicals produced in the X, A, and B states from the photodissociation of phenol at 213 and 193 nm are reported. Unlike other techniques, this method has no interference from the undissociated hot molecules. It can ultimately become a standard pump-probe technique for the study of large molecule photodissociation in multistates.

Unambiguous Signature of the Berry Phase in Intense Laser Dissociation of Diatomic Molecules

We report strong evidence of Berry phase effects in intense laser dissociation of D₂⁺ molecules, manifested as Aharonov-Bohm-like oscillations in the photofragment angular distribution (PAD). Our calculations show that this interference pattern strongly depends on the parity of the diatom initial rotational state, (-1) ^j. Indeed, the PAD local maxima (minima) observed in one case ( j odd) correspond to local minima (maxima) in the other case ( j even). Using simple topological arguments, we clearly show that such interference conversion is a direct signature of the Berry phase. The sole effect of the latter on the rovibrational wave function is a sign change of the relative phase between two interfering components, which wind in opposite senses around a light-induced conical intersection (LICI). Therefore, encirclement of the LICI leads to constructive ( j odd) or destructive ( j even) self-interference of the initial nuclear wavepacket in the dissociative limit. To corroborate our theoretical findings, we suggest an experiment of strong-field indirect dissociation of D₂⁺ molecules, comparing the PAD of the ortho and para molecular species in directions nearly perpendicular to the laser polarization axis.

Conformer specific nonadiabatic reaction dynamics in the photodissociation of partially deuterated thioanisoles (C6H5S-CH2D and C6H5S-CHD2)

In this work, we have investigated nonadiabatic dynamics in the vicinity of conical intersections for predissociation reactions of partially deuterated thioanisole molecules: C₆H₅S-CH₂D and C₆H₅S-CHD₂. Each isotopomer has two distinct rotational conformers according to the geometrical position of D or H of the methyl moiety with respect to the molecular plane for C₆H₅S-CH₂D or C₆H₅S-CHD₂, respectively, as spectroscopically characterized in our earlier report [J. Lee, S.-Y. Kim and S. K. Kim, J. Phys. Chem. A, 2014, 118, 1850]. Since identification and separation of two different rotational conformers of each isotopomer have been unambiguously done, we could interrogate nonadiabatic dynamics of thioanisole in terms of both H/D substitutional and conformational structural effects. Nonadiabatic transition probability, estimated by the experimentally measured branching ratio of the nonadiabatically produced ground-state channel giving C₆H₅S·(X[combining tilde]) versus the adiabatic excited-state channel leading to the C₆H₅S·(Ã) radical, shows resonance-like increases at symmetric (ν_s) or asymmetric (7a) S-CH₂D (or S-CHD₂) stretching mode excitation in S₁ for all conformational isomers of two isotopomers. However, absolute probabilistic value of the nonadiabatic transition is found to vary quite drastically depending on different conformers and isotopomers. The experimental finding that nonadiabatic transition dynamics are very sensitive to subtle changes in the nuclear configuration within the Franck-Condon region induced by the H/D substitution indicates that the S₁/S₂ conical intersection seam is quite narrowly defined in the multi-dimensional nuclear configurational space as far as the S-methyl predissociation reaction is concerned. In order to understand the relation between molecular structure and nonadiabaticity of reaction, potential energy surfaces near S₁/S₂ conical intersections have been theoretically calculated along ν_s and 7a normal mode coordinates for all conformational isomers. Slow-electron velocity map imaging (SEVI) spectroscopy is employed to unravel the extent of intramolecular vibrational redistribution (IVR) for particular mode excitations of S₁, providing insights into the dynamic interplay between IVR and nonadiabatic transition probability near the conical intersection seam.

An alternative laser driven photodissociation mechanism of pyrrole via πσ*1∕S0 conical intersection

A first principles quantum dynamics study of N-H photodissociation of pyrrole on the S₀-¹πσ^*(A21) coupled electronic states is carried out with the aid of an optimally designed UV-laser pulse. A new photodissociation path, as compared to the conventional barrier crossing on the πσ*1 state, opens up upon electronic transitions under the influence of pump-dump laser pulses, which efficiently populate both the dissociation channels. The interplay of electronic transitions due both to vibronic coupling and the laser pulse is observed in the control mechanism and discussed in detail. The proposed control mechanism seems to be robust, and not discussed in the literature so far, and is expected to trigger future experiments on the πσ*1 photochemistry of molecules of chemical and biological importance. The design of the optimal pulses and their application to enhance the overall dissociation probability is carried out within the framework of optimal control theory. The quantum dynamics of the system in the presence of pulse is treated by solving the time-dependent Schrödinger equation in the semi-classical dipole approximation.

Crossing Electronic States in the Franck-Condon Zone of Carbon Dioxide: A Five-Fold Closed Seam of Conical and Glancing Intersections

Intersections of electronic states are recognized as natural molecular hubs routing initial photoexcitation toward specific products. In this Letter, we investigate a five-fold intersection between the 2,3(1)A' and 1,2,3(1)A″ electronic states of CO2 responsible for the ultraviolet absorption of the molecule between 120 and 160 nm. It is demonstrated that the intersection is not an isolated point but a seam with the following properties: (1) The seam comprises both conical and glancing intersections. (2) The seam is a closed loop traced by the concerted displacements of two CO bonds. (3) The seam is constrained to linear geometries and confined to the CO bond plane. (4) The seam passes directly through the Franck-Condon zone. Topography and properties of the intersection seam are validated using assignment of the electronic wave functions and illustrated using circulation of the derivative couplings along closed contours. The possibility of observation of strongly curved seams is briefly discussed.

Direct Observation of Hydrogen Tunneling Dynamics in Photoexcited Phenol

The excited-state dynamics of phenol following ultraviolet (UV) irradiation have received considerable interest in recent years, most notably because they can provide a model for understanding the UV-induced dynamics of the aromatic amino acid tyrosine. Despite this, there has been some debate as to whether hydrogen tunneling dynamics play a significant role in phenol's excited-state O-H bond fission when UV excitation occurs below the (1)ππ*/(1)πσ* conical intersection (CI). In this Letter, we present direct evidence that (1)πσ*-mediated O-H bond fission below the (1)ππ*/(1)πσ* CI proceeds exclusively through hydrogen tunneling dynamics. The observation of hydrogen tunneling may have some parallels with proton tunneling dynamics from tyrosine residues (along the O-H bond of the phenol moiety) in a wide range of natural enzymes, potentially adding further justification for utilizing phenols as model systems for investigating tyrosine-based dynamics.

Experimental probing of conical intersection dynamics in the photodissociation of thioanisole

Chemical reactions that occur in the ground electronic state are described well by invoking the Born-Oppenheimer approximation, which allows their development to be rationalized by nuclear rearrangements that smoothly traverse an adiabatic potential energy surface. The situation is different, however, for reactions in electronically excited states, where non-adiabatic transitions occur between adiabatic surfaces. The conical intersection, in which two adiabatic surfaces touch, is accepted widely as the dynamic funnel for efficient non-adiabatic transitions, but its direct experimental probing is rare. Here, we investigate the photodissociation of thioanisole and observe a striking dependence of the relative yields of two reaction channels on the photoexcitation energy as indicated by a dynamic resonance in the product branching ratio. This results from the interference of two different adiabatic states that are in close proximity in the region of a conical intersection. The location of the observed resonance on the multidimensional potential energy surface thus reveals the nuclear configuration of the conical intersection and its dynamic role in the non-adiabatic transition.

Photodissociation dynamics of thiophenol-d1: the nature of excited electronic states along the S-D bond dissociation coordinate

The S-D bond dissociation dynamics of thiophenol-d1 (C6H5SD) pumped at 266, 243, and 224 nm are examined using the velocity map ion imaging technique. At both 266 and 243 nm, distinct peaks associated with X and A states of the phenylthiyl radical (C6H5S*) are observed in the D+ image at high and low kinetic energy regions, respectively. The partitioning of the available energy into the vibrational energy of the phenylthiyl radical is found to be enhanced much more strongly at 266 nm compared to that at 243 nm. This indicates that the pipi* electronic excitation at 266 nm is accompanied by significant vibrational excitation. Given the relatively large anisotropy parameter of -0.6, the S-D dissociation at 266 nm is prompt and should involve the efficient coupling to the upper-lying n(pi)sigma* repulsive potential energy surface. The optical excitation of thiophenol at 224 nm is tentatively assigned to the pisigma* transition, which leads to the fast dissociation on the repulsive potential energy surface along the S-D coordinate. The nature of the electronic transitions associated with UV absorption bands is investigated with high-level ab initio calculations. Excitations to different electronic states of thiophenol result in unique branching ratios and vibrational excitations for the fragment of the phenylthiyl radical in the two lowest electronic states.

A computational study on the mechanism of intramolecular oxo-hydroxy phototautomerism driven by repulsive pi sigma* state

Potential energy (PE) surfaces of the lowest excited states of the 4(3H)-pyrimidinone/4-hydroxypyrimidine system were investigated with the aid of the CC2 and CASSCF methods of the electronic structure theory. These studies resulted in identification of a low-lying pi sigma* state, which is dissociative with respect to the stretching of the N-H or O-H bonds in the oxo and hydroxy structures of the compound, respectively. After initial excitation to the lowest local n pi* and/or pi pi* singlet states, the system can access the PE surface of the pi sigma* state by crossing a low barrier. It was computationally demonstrated that the system should evolve on the PE surface of the repulsive pi sigma* state toward a broad seam of intersection with the PE surface of the ground state. At the intersection, the nonadiabatic transition to the ground electronic state takes place and the system can either evolve to a minimum of the initially excited tautomer or to the ground-state minimum of the other tautomer. The steps listed above provide a mechanism of photoinduced dissociation-association (PIDA) phototautomerism, experimentally observed for a number of monomeric molecules, structurally similar to 4(3H)-pyrimidinone/4-hydroxypyrimidine. This mechanism describes a new class of intramolecular phototautomeric reactions driven by a repulsive pi sigma* state.