Ultrafast Dynamics and Coherent Oscillations in Ethylene and Ethylene-d4 Excited at 162 nm

Generalized Semiclassical Ehrenfest Method: A Route to Wave Function-Free Photochemistry and Nonadiabatic Dynamics with Only Potential Energies and Gradients

We reconsider recent methods by which direct dynamics calculations of electronically nonadiabatic processes can be carried out while requiring only adiabatic potential energies and their gradients. We show that these methods can be understood in terms of a new generalization of the well-known semiclassical Ehrenfest method. This is convenient because it eliminates the need to evaluate electronic wave functions and their matrix elements along the mixed quantum-classical trajectories. The new approximations and procedures enabling this advance are the curvature-driven approximation to the time-derivative coupling, the generalized semiclassical Ehrenfest method, and a new gradient correction scheme called the time-derivative matrix (TDM) scheme. When spin-orbit coupling is present, one can carry out dynamics calculations in the fully adiabatic basis using potential energies and gradients calculated without spin-orbit coupling plus the spin-orbit coupling matrix elements. Even when spin-orbit coupling is neglected, the method is useful because it allows calculations by electronic structure methods for which nonadiabatic coupling vectors are unavailable. In order to place the new considerations in context, the article starts out with a review of background material on trajectory surface hopping, the semiclassical Ehrenfest scheme, and methods for incorporating decoherence. We consider both internal conversion and intersystem crossing. We also review several examples from our group of successful applications of the curvature-driven approximation.

Quantum Quality with Classical Cost: Ab Initio Nonadiabatic Dynamics Simulations Using the Mapping Approach to Surface Hopping

Nonadiabatic dynamics methods are an essential tool for investigating photochemical processes. In the context of employing first-principles electronic structure techniques, such simulations can be carried out in a practical manner using semiclassical trajectory-based methods or wave packet approaches. While all approaches applicable to first-principles simulations are necessarily approximate, it is commonly thought that wave packet approaches offer inherent advantages over their semiclassical counterparts in terms of accuracy and that this trait simply comes at a higher computational cost. Here we demonstrate that the mapping approach to surface hopping (MASH), a recently introduced trajectory-based nonadiabatic dynamics method, can be efficiently applied in tandem with ab initio electronic structure. Our results even suggest that MASH may provide more accurate results than on-the-fly wave packet techniques, all at a much lower computational cost.

Benchmarking non-adiabatic quantum dynamics using the molecular Tully models

On-the-fly non-adiabatic dynamics methods are becoming more important as tools to characterise the time evolution of a system after absorbing light. These methods, which calculate quantities such as state energies, gradients and interstate couplings at every time step, circumvent the requirement for pre-computed potential energy surfaces. There are a number of different algorithms used, the most common being Tully Surface Hopping (TSH), but all are approximate solutions to the time-dependent Schrödinger equation and benchmarking is required to understand their accuracy and performance. For this, a common set of systems and observables are required to compare them. In this work, we validate the on-the-fly direct dynamics variational multi-configuration Gaussian (DD-vMCG) method using three molecular systems recently suggested by Ibele and Curchod as molecular versions of the Tully model systems used to test one-dimensional non-adiabatic behaviour [Ibele et al., Phys. Chem. Chem. Phys. 2020, 22, 15183-15196]. Parametrised linear vibronic potential energy surfaces for each of the systems were also tested and compared to on-the-fly results. The molecules, which we term the Ibele-Curchod models, are ethene, DMABN and fulvene and the authors used them to test and compare several versions of the Ab Initio Multiple Spawning (AIMS) method alongside TSH. The three systems present different deactivation pathways after excitation to their ππ* bright states. When comparing DD-vMCG to AIMS and TSH, we obtain crucial differences in some cases, for which an explanation is provided by the classical nature and the chosen initial conditions of the TSH simulations.

Nonadiabatic simulations of photoisomerization and dissociation in ethylene using ab initio classical trajectories

We simulate the nonadiabatic dynamics of photo-induced isomerization and dissociation in ethylene using ab initio classical trajectories in an extended phase space of nuclear and electronic variables. This is achieved by employing the linearized semiclassical initial value representation method for nonadiabatic dynamics, where discrete electronic states are mapped to continuous classical variables using either the Meyer-Miller-Stock-Thoss representation or a more recently introduced spin mapping approach. Trajectory initial conditions are sampled by constraining electronic state variables to a single initial excited state and by drawing nuclear phase space configurations from a Wigner distribution at a finite temperature. An ensemble of classical ab initio trajectories is then generated to compute thermal population correlation functions and analyze the mechanisms of isomerization and dissociation. Our results serve as a demonstration that this parameter-free semiclassical approach is computationally efficient and accurate, identifying mechanistic pathways in agreement with previous theoretical studies and also uncovering dissociation pathways observed experimentally.

Scoring molecular wires subject to an ultrafast laser pulse for molecular electronic devices

A nonionizing ultrafast laser pulse of 20-fs duration with a peak amplitude electric-field ±E = 200 × 10^-4 a.u. was simulated. It was applied to the ethene molecule to consider its effect on the electron dynamics, both during the application of the laser pulse and for up to 100 fs after the pulse was switched off. Four laser pulse frequencies ω = 0.2692, 0.2808, 0.2830, and 0.2900 a.u. were chosen to correspond to excitation energies mid-way between the (S₁ ,S₂ ), (S₂ ,S₃ ), (S₃ ,S₄ ) and (S₄ ,S₅ ) electronic states, respectively. Scalar quantum theory of atoms in molecules (QTAIM) was used to quantify the shifts of the C1C2 bond critical points (BCPs). Depending on the frequencies ω selected, the C1C2 BCP shifts were up to 5.8 times higher after the pulse was switched off compared with a static E-field with the same magnitude. Next generation QTAIM (NG-QTAIM) was used to visualize and quantify the directional chemical character. In particular, polarization effects and bond strengths, in the form of bond-rigidity vs. bond-flexibility, were found, for some laser pulse frequencies, to increase after the laser pulse was switched off. Our analysis demonstrates that NG-QTAIM, in partnership with ultrafast laser irradiation, is useful as a tool in the emerging field of ultrafast electron dynamics, which will be essential for the design, and control of molecular electronic devices.

The vertical excitation energies and a lifetime of the two lowest singlet excited states of the conjugated polyenes from C2 to C22: Ab initio, DFT, and semiclassical MNDO-MD simulations

Electronic excited states in the series of polyene molecules were explored. Optimal ground-state geometry was used for the evaluation of vertical excitation energies. Results of a chosen set of functionals were compared to post-HF methods (EOM-CCSD, NEVPT2, CASPT2, and MRCI). In addition, the semiempirical OM2/MNDO method using MRCISD computational level was confronted with the above-mentioned techniques. Despite the fact that the first excited state has a significant double-excitation character some functionals were able to qualitatively determine the correct state order (where the lowest excited state has a A g - character). The most successful functionals in transition energies predictions were PBE, TPSS and BLYP in Tamm-Dancoff approach (TDA), which had the smallest root-mean-square deviation (RMSD) scoring towards the experimental values. Regarding RMSD scoring, the OM2/MNDO method performed fairly well, too. Besides absorption spectra, lifetimes of the first two excited states were estimated based on a stochastic approach exploring a swarm of OM2/MNDO hopping dynamics using the Tully fewest switch algorithm for each molecule. The longest lifetime of the first excited state (S₁ ) was found for decapentaene (about 5 ps). Further elongation of the conjugated chain caused a mild decrease of this value to ca 1.5 ps for docosaundecaene.

Few-Femtosecond Isotope Effect in Polyatomic Molecules Ionized by Extreme Ultraviolet Attosecond Pulse Trains

Following ionization by an extreme ultraviolet (XUV) attosecond pulse train, a polyatomic molecule can be promoted to more-than-one excited states of the residual ion. The ensuing relaxation dynamics is often facilitated by several reaction coordinates, making them difficult to disentangle by the usual spectroscopic means. Here, we show that in atto-chemistry isotope labeling can be an efficient tool for unraveling the relaxation pathways in highly excited photoionized molecules. Employing an XUV pump pulse and a near-infrared probe pulse, we found the nuclear as well as coupled electron-nuclear dynamics in ethylene to be almost 40% faster compared to that of its deuterated counterpart. The findings, which are supported by advanced nonadiabatic dynamics calculations, led to the identification of the relevant nuclear coordinates controlling the relaxation. Our experiment highlights the relevance of ultrashort XUV pulses to capture the isotopic effect in few-femtosecond molecular photodynamics.

Ultrafast Photoisomerization of Ethylene Studied Using Time-Resolved Extreme Ultraviolet Photoelectron Spectroscopy

The photoisomerization of isolated ethylene (ethene) was observed in real time from the Franck-Condon region in the ¹ππ* state to ground-state products using time-resolved photoelectron spectroscopy with extreme ultraviolet (EUV, 21.7 eV) probe pulses. A combination of filamentation four-wave mixing and high-order harmonic generation was employed to obtain a temporal resolution of 31 ± 2 fs. The nuclear wave packet created by a 160 nm pump pulse accesses C═C twisted geometries within 10 fs, and the population transfer from the excited to the ground state occurs within the next 20-30 fs. Formation of vibrationally highly excited ground-state molecules was observed in less than 45 fs, and they decayed with two time constants of 0.87 and >5 ps. The interpretation of the photoelectron spectra is supported by vertical ionization energies calculated using XMS-CASPT2 along geodesically interpolated reaction paths from the Franck-Condon region to the products.

Fewest switches surface hopping with Baeck-An couplings

In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.

Fewest switches surface hopping with Baeck-An couplings

In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.

Velocity Adjustment in Surface Hopping: Ethylene as a Case Study of the Maximum Error Caused by Direction Choice

The most common surface hopping dynamics algorithms require velocity adjustment after hopping to ensure total-energy conservation. Based on the semiclassical analysis, this adjustment must be made parallel to the nonadiabatic coupling vector's direction. Nevertheless, this direction is not always known, and the common practice has been to adjust the velocity in either the linear momentum or velocity directions. This paper benchmarks surface hopping dynamics of photoexcited ethylene with velocity adjustment in several directions, including those of the nonadiabatic coupling vector, the momentum, and the energy gradient difference. It is shown that differences in time constants and structural evolution fall within the statistical uncertainty of the method considering up to 500 trajectories in each dynamics set, rendering the three approaches statistically equivalent. For larger ensembles beyond 1000 trajectories, significant differences between the results arise, limiting the validity of adjustment in alternative directions. Other possible adjustment directions (velocity, single-state gradients, angular momentum) are evaluated as well. Given the small size of ethylene, the results reported in this paper should be considered an upper limit for the error caused by the choice of the velocity-adjustment direction on surface hopping dynamics.

Accessing the Activation Mechanisms of Ethylene Photo-Polymerization under Pressure by Transient Infrared Absorption Spectroscopy

The ambient temperature photoinduced polymerization of compressed (P < 1 GPa) fluid ethylene was characterized by transient infrared absorption spectroscopy with a resolution of few nanoseconds, 3 orders of magnitude higher than previously reported. The reaction has been studied under both one- and two-photon excitation evidencing in the latter case its occurrence only in the presence of different transition metal oxides. Their photocatalytic activity is ascribed to the stabilization of the excited biradicals through electron density exchange between the d orbitals of the metal and the π antibonding orbitals of ethylene which lengthens the lifetime of the biradicals. In both one- and two-photon activation cases the polymerization is characterized by an initial step distinguished by a molecularity of 0.15 ± 0.02 identified as the activation step of the reaction lasting, in the one-photon excitation case, a few hundreds of nanoseconds. Using pulsed excitation the reaction evolves toward a free radical polymerization only under one-photon excitation whereas the critical concentration of radicals required to propagate the reaction is never achieved in the two-photon excitation case. Comparison with continuous wave excitation unambiguously identifies in the average power released to the sample the key factor to drive quantitatively and qualitatively the polymerization.

Photochemistry of Amylene Double Bond in Clusters on Free Argon Nanoparticles

We have investigated reactivity of double bond in 2-methyl-2-butene (also trimethylethylene or amylene) in the excited and ionized states. In a combined experimental and theoretical study, we focused on both the intermolecular and intramolecular reactions. In a molecular beam experiment, we have sequentially picked up several amylene molecules on the surface of argon nanoparticles Ar_M, M̅ ≈ 90, acting as a cold support. Ionization with 70 eV electrons yields mass spectra strongly dominated by amylene cluster ions Am(Am)_n⁺. Interestingly, upon multiphoton ionization with 193 nm (6.4 eV) photons, a new strong fragment series appears in the spectra, nominally corresponding to an addition of two carbon atoms, i.e., (Am)_nC₂⁺. This difference between electron and photoionization suggests a reaction in an excited state of amylene with a neighboring amylene molecule. We used techniques of nonadiabatic molecular dynamics to study the reactivity of amylene molecules both in the excited and in ionized states. Possible reaction pathways are proposed, substantiating the observed differences between the electron ionization and photoionization mass spectra.

The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

Quasi-Diabatic Scheme for Nonadiabatic On-the-Fly Simulations

We use the quasi-diabatic (QD) propagation scheme to perform on-the-fly nonadiabatic simulations of the photodynamics of ethylene. The QD scheme enables a seamless interface between accurate diabatic-based quantum dynamics approaches and adiabatic electronic structure calculations, explicitly avoiding any efforts to construct global diabatic states or reformulate the diabatic dynamics approach to the adiabatic representation. Using the partial linearized path-integral approach and the symmetrical quasi-classical approach as the diabatic dynamics methods, the QD propagation scheme enables direct nonadiabatic simulation with complete active space self-consistent field on-the-fly electronic structure calculations. The population dynamics obtained from both approaches are in a close agreement with the quantum wavepacket-based method and outperform the widely used trajectory surface-hopping approach. Further analysis of the ethylene photodeactivation pathways demonstrates the correct predictions of competing processes of nonradiative relaxation mechanism through various conical intersections. This work provides the foundation of using accurate diabatic dynamics approaches and on-the-fly adiabatic electronic structure information to perform ab initio nonadiabatic simulation.

The 3s Rydberg state as a doorway state in the ultrafast dynamics of 1,1-difluoroethylene

The deactivation dynamics of 1,1-difluoroethylene after light excitation is studied within the surface hopping formalism in the presence of 3s and 3p Rydberg states using multi-state second order perturbation theory (MS-CASPT2).

Photoinduced C–H bond fission in prototypical organic molecules and radicals

We survey and assess current knowledge regarding the primary photochemistry of hydrocarbon molecules and radicals.

Time-resolved nuclear dynamics in bound and dissociating acetylene

We have investigated nuclear dynamics in bound and dissociating acetylene molecular ions in a time-resolved reaction microscopy experiment with a pair of few-cycle pulses. Vibrating bound acetylene cations or dissociating dications are produced by the first pulse. The second pulse probes the nuclear dynamics by ionization to higher charge states and Coulomb explosion of the molecule. For the bound cations, we observed vibrations in acetylene (HCCH) and its isomer vinylidene (CCHH) along the CC-bond with a periodicity of around 26 fs. For dissociating dication molecules, a clear indication of enhanced ionization is found to occur along the CH- and CC-bonds after 10 fs to 40 fs. The time-dependent ionization processes are simulated using semi-classical on-the-fly dynamics revealing the underling mechanisms.

Time-dependent view of an isotope effect in electron-nuclear nonequilibrium dynamics with applications to N2

Isotopic fractionation in the photodissociation of N₂ could explain the considerable variation in the ¹⁴N/¹⁵N ratio in different regions of our galaxy. We previously proposed that such an isotope effect is due to coupling of photoexcited bound valence and Rydberg electronic states in the frequency range where there is strong state mixing. We here identify features of the role of the mass in the dynamics through a time-dependent quantum-mechanical simulation. The photoexcitation of N₂ is by an ultrashort pulse so that the process has a sharply defined origin in time and so that we can monitor the isolated molecule dynamics in time. An ultrafast pulse is necessarily broad in frequency and spans several excited electronic states. Each excited molecule is therefore not in a given electronic state but in a superposition state. A short time after excitation, there is a fairly sharp onset of a mass-dependent large population transfer when wave packets on two different electronic states in the same molecule overlap. This coherent overlap of the wave packets on different electronic states in the region of strong coupling allows an effective transfer of population that is very mass dependent. The extent of the transfer depends on the product of the populations on the two different electronic states and on their relative phase. It is as if two molecules collide but the process occurs within one molecule, a molecule that is simultaneously in both states. An analytical toy model recovers the (strong) mass and energy dependence.

Ultrafast Relaxation Dynamics of the Ethylene Cation C(2)H(4)+

We present a combined experimental and computational study of the relaxation dynamics of the ethylene cation. In the experiment, we apply an extreme-ultraviolet-pump/infrared-probe scheme that permits us to resolve time scales on the order of 10 fs. The photoionization of ethylene followed by an infrared (IR) probe pulse leads to a rich structure in the fragment ion yields reflecting the fast response of the molecule and its nuclei. The temporal resolution of our setup enables us to pinpoint an upper bound of the previously defined ethylene-ethylidene isomerization time to 30 ± 3 fs. Time-dependent density functional based trajectory surface hopping simulations show that internal relaxation between the first excited states and the ground state occurs via three different conical intersections. This relaxation unfolds on femtosecond time scales and can be probed by ultrashort IR pulses. Through this probe mechanism, we demonstrate a route to optical control of the important dissociation pathways leading to separation of H or H2.

How does tetraphenylethylene relax from its excited states?

Photocyclization play a key role in the deactivation mechanism of tetraphenylethylene.

Probing the dynamics of highly excited toluene on the fs timescale

Investigation of the dynamics of toluene-h8 (C6H5CH3), toluene-d8 (C6D5CD3) and toluene-α,α,α-d3 (C6H5CD3) has been performed utilizing the VUV pump-IR probe technique on the fs timescale. Using the 5th harmonic (∼160 nm) of a Ti:sapphire laser as the pump beam, two superimposed electronic states, the valence S3 and the Rydberg 4p, were excited by one-photon absorption, followed by ionization and dissociation induced by the probe beam (800 nm). Analysis of the transient signal of the parent (P(+)) and fragment ions ([P-H](+) or [P-D](+)) implies the existence of two different relaxation processes: (i) from the Rydberg and (ii) from the S3 valence state. Using a rate equation model, the decay times have been determined and comparison between the different isotopologues has been made. Conclusions on the relaxation path, the relative displacements of the potential energy surfaces and the activation energies needed have been drawn from the decay times. The signals corresponding to the fragment ions present a small in amplitude, but nonetheless, unambiguous periodical modulation, which is attributed to out-of-plane bending oscillation, involving also the methyl group. The dynamics of the H- and D-loss channels has been investigated. Especially for the case of toluene-α,α,α-d3, where both channels are in operation, it was found that the ratio of the abundance of H/D-loss dissociation reactions decreases as the pump-probe delay time increases.

Efficient generation of below-threshold harmonics for high-fidelity multi-photon physics in the VUV spectral range

We demonstrate the generation of microjoule level, sub-20-fs, Ti:Sa fifth-harmonic pulses utilizing a loose-focusing geometry in a long Ar gas cell. The VUV pulses centered at 161.8 nm reach pulse energies of 1.1 μJ per pulse, while the corresponding pulse duration is measured with a second-order, fringe-resolved autocorrelation scheme to be 18±1  fs. Nonresonant, two-photon ionization of Kr and three-photon ionization of Ne verify the fifth-harmonic pulse high-intensity content and indicate the feasibility of multi-photon VUV pump-VUV probe studies of ultrafast atomic and molecular dynamics.

Resolving multi-exciton generation by attosecond spectroscopy

We propose an experimentally viable attosecond transient absorption spectroscopy scheme to resolve controversies regarding multiexciton (ME) generation in nanoscale systems. Absence of oscillations indicates that light excites single excitons, and MEs are created by incoherent impact ionization. An oscillation indicates the coherent mechanism, involving excitation of superpositions of single and MEs. The oscillation decay, ranging from 5 fs at ambient temperature to 20 fs at 100 K, gives the elastic exciton-phonon scattering time. The signal is best observed with multiple-cycle pump pulses.

Photochemistry and photophysics at extended seams of conical intersection

The role of extended seams of conical intersection in excited-state mechanisms is reviewed. Seams are crossings of the potential energy surface in many dimensions where the decay from the excited to the ground state can occur, and the extended seam is composed of different segments lying along a reaction coordinate. Every segment is associated with a different primary photoproduct, which gives rise to competing pathways. This idea is first illustrated for fulvene and ethylene, and then it is used to explain more complex cases such as the dependence of the isomerisation of retinal chromophore isomers on the protein environment, the dependence of the efficiency of the azobenzene photochemical switch on the wavelength of irradiation and the direction of the isomerisation, and the coexistence of different mechanisms in the photo-induced Wolff rearrangement of diazonaphthoquinone. The role of extended seams in the photophysics of the DNA nucleobases and the relationship between two-state seams and three-state crossings is also discussed. As an outlook, the design of optical control strategies based on the passage of the excited molecule through the seam is considered, and it is shown how the excited-state lifetime of fulvene can be modulated by shaping the energy of the seam.

Cyclobutadiene automerization and rotation of ethylene: energetics of the barriers by using spin-polarized wave functions

Spin-projected spin polarized Møller-Plesset and spin polarized coupled clusters calculations have been made to estimate the cyclobutadiene automerization, the ethylene torsion barriers in their ground state, and the gap between the singlet and triplet states of ethylene. The results have been obtained optimizing the geometries at MP4 and/or CCSD levels, by an extensive Gaussian basis set. A comparative analysis with more complex calculations, up to MP5 and CCSDTQP, together with others from the literature, have also been made, showing the efficacy of using spin-polarized wave functions as a reference wave function for Møller-Plesset and coupled clusters calculations, in such problems.

Hexamethylcyclopentadiene: time-resolved photoelectron spectroscopy and ab initio multiple spawning simulations

Time-resolved photoelectron spectroscopy and ab initio multiple spawning dynamical simulations of hexamethylcyclopentadiene reveal wavepacket evolution in a distinct degree of freedom.

Nonlinear dimensionality reduction for nonadiabatic dynamics: the influence of conical intersection topography on population transfer rates

Conical intersections play a critical role in the nonadiabatic relaxation of excited electronic states. However, there are an infinite number of these intersections and it is difficult to predict which are actually relevant. Furthermore, traditional descriptors such as intrinsic reaction coordinates and steepest descent paths often fail to adequately characterize excited state reactions due to their highly nonequilibrium nature. To address these deficiencies in the characterization of excited state mechanisms, we apply a nonlinear dimensionality reduction scheme (diffusion mapping) to generate reaction coordinates directly from ab initio multiple spawning dynamics calculations. As illustrated with various examples of photoisomerization dynamics, excited state reaction pathways can be derived directly from simulation data without any a priori specification of relevant coordinates. Furthermore, diffusion maps also reveal the influence of intersection topography on the efficiency of electronic population transfer, providing further evidence that peaked intersections promote nonadiabatic transitions more effectively than sloped intersections. Our results demonstrate the usefulness of nonlinear dimensionality reduction techniques as powerful tools for elucidating reaction mechanisms beyond the statistical description of processes on ground state potential energy surfaces.