Search

Photochemistry and photophysics processes involve structures far from equilibrium. In these reactions, there is often strong coupling between nuclear and electronic degrees of freedom. For first-row transition metals, K beta X-ray emission spectroscopy (XES) is a sensitive probe of electronic structure due to the direct overlap between the valence orbitals and the 3p hole in the final state. Here the sensitivity of K beta mainline (K beta(1,3)) XES to structural dynamics is analyzed by simulating spectral changes along the excited state dynamics of an iron photosensitizer [Fe-II(bmip)(2)](2+) [bmip = 2,6-bis(3-methyl-imidazole-1-ylidine)-pyridine], using both restricted active space (RAS) multiconfigurational wavefunction theory and a one-electron orbital-energy approach in density-functional theory (1-DFT). Both methods predict a spectral blue-shift with increasing metal-ligand distance, which changes the emission intensity for any given detection energy. These results support the suggestion that the [Fe-II(bmip)(2)](2+) femtosecond K beta XES signal shows oscillations due to coherent wavepacket dynamics. Based on the RAS results, the sensitivity to structural dynamics is twice as high for K beta compared to K alpha, with the drawback of a lower signal-to-noise ratio. K beta sensitivity is favored by a larger spectral blue-shift with increasing metal-ligand distance and larger changes in spectral shape. Comparing the two simulations methods, 1-DFT predicts smaller energy shifts and lower sensitivity, likely due to missing final-state effects. The simulations can be used to design and interpret XES probes of non-equilibrium structures to gain mechanistic insights in photocatalysis.

ESIPT rate, yield and mechanism are revealed in pyrrol pyridine using non-adiabatic dynamics simulations.

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.

Chemiexcitation, the generation of electronic excited states by a thermal reaction initiated on the ground state, is an essential step in chemiluminescence, and it is mediated by the presence of a conical intersection that allows a nonadiabatic transition from ground state to excited state. Conical intersections classified as sloped favor chemiexcitation over ground state relaxation. The chemiexcitation yield of 1,2-dioxetanes is known to increase upon methylation. In this work we explore to which extent this trend can be attributed to changes in the conical intersection topography or accessibility. Since conical intersections are not isolated points, but continuous seams, we locate regions of the conical intersection seams that are close to the configuration space traversed by the molecules as they react on the ground state. We find that conical intersections are energetically and geometrically accessible from the reaction trajectory, and that topographies favorable to chemiexcitation are found in all three molecules studied. Nevertheless, the results suggest that dynamic effects are more important for explaining the different yields than the static features of the potential energy surfaces.

We present a new methodology for measuring few-femtosecond electronic and nuclear dynamics in both atoms and polyatomic molecules using multidimensional high harmonic generation (HHG) spectroscopy measurements, in which the spectra are recorded as a function of the laser intensity to form a two-dimensional data set. The method is applied to xenon atoms and to benzene molecules, the latter exhibiting significant fast nuclear dynamics following ionization. We uncover the signature of the sub-cycle evolution of the returning electron flux in strong-field ionized xenon atoms, implicit in the strong field approximation but not previously observed directly. We furthermore extract the nuclear autocorrelation function in strong field ionized benzene cations, which is determined to have a decay of $$\tau _0 = 4 \pm 1$$ τ 0 = 4 ± 1 fs, in good agreement with the $$ \tau _0 = 3.5$$ τ 0 = 3.5 fs obtained from direct dynamics variational multi-configuration Gaussian calculations. Our method requires minimal assumptions about the system, and is applicable even to un-aligned polyatomic molecules.

The cis → trans photo-isomerisation mechanism of azobenzene, after excitation to the nπ* and ππ* states, is revisited using high-level ab initio surface hopping mixed quantum-classical dynamics in combination with multi-reference CASSCF electronic structure calculations. A reduction of photoisomerisation quantum yield of 0.10 on exciting to the higher energy ππ* state compared to the lower energy nπ* state is obtained, in close agreement with the most recent experimental values [Ladanyi et al., Photochem. Photobiol. Sci., 2017, 16, 1757–1761] which re-examined previous literature values which showed larger changes in quantum yield. By direct comparison of both excitations, we have found that the explanation for the decrease in quantum yield is not the same as for the reduction observed in the trans → cis photoisomerisation. In contrast to the trans → cis scenario, S1 → S0 decay does not occur at ‘earlier’ C–NN–C angles along the central torsional coordinate after ππ* excitation, as in the cis → trans case the rotation about this coordinate occurs too rapidly. The wavelength dependency of the quantum yield is instead found to be due to a potential well on the S2 surface, from which either cis or trans-azobenzene can be formed. While this well is accessible after both excitations, it is more easily accessed after ππ* excitation – an additional 15–17% of photochromes, which under nπ* excitation would have exclusively formed the trans isomer, are trapped in this well after ππ* excitation. The probability of forming the cis isomer when leaving this well is also higher after ππ* excitation, increasing from 9% to 35%. The combination of these two factors results in the reduction of 0.10 of the quantum yield of photoisomerisation on ππ* excitation of cis-azobenzene, compared to nπ* excitation.

Chemiexcitation is the population of electronic excited states from the electronic ground state via radiationless non-adiabatic transitions upon thermal activation. The subsequent emission of the e ...

International audience; Controlling matter with light has always been a great challenge, leading to the ever-expanding field of photochemistry. In addition, since the first generation of light pulses of attosecond (1 as = 10–18 s) duration, a great deal of effort has been devoted to observing and controlling electrons on their intrinsic time scale. Because of their short duration, attosecond pulses have a large spectral bandwidth populating several electronically excited states in a coherent manner, i.e., an electronic wavepacket. Because of interference, such a wavepacket has a new electronic distribution implying a potentially different and totally new reactivity as compared to traditional photochemistry, leading to the novel concept of “attochemistry”. This nascent field requires the support of theory right from the start. In this Perspective, we discuss the opportunities offered by attochemistry, the related challenges, and the current and future state-of-the-art developments in theoretical chemistry needed to model it accurately.

The possibility of suddenly ionized molecules undergoing extremely fast electron hole (or hole) dynamics prior to significant structural change was first recognized more than 20 years ago and termed charge migration. The accurate probing of ultrafast electron hole dynamics requires measurements that have both sufficient temporal resolution and can detect the localization of a specific hole within the molecule. We report an investigation of the dynamics of inner valence hole states in isopropanol where we use an x-ray pump-x-ray probe experiment, with site and state-specific probing of a transient hole state localized near the oxygen atom in the molecule, together with an ab initio theoretical treatment. We record the signature of transient hole dynamics and make the first tentative observation of dynamics driven by frustrated Auger-Meitner transitions. We verify that the effective hole lifetime is consistent with our theoretical prediction. This state-specific measurement paves the way to widespread application for observations of transient hole dynamics localized in space and time in molecules and thus to charge transfer phenomena that are fundamental in chemical and material physics. (Less)

International audience

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions. ispartof: JOURNAL OF CHEMICAL PHYSICS vol:152 issue:21 ispartof: location:United States status: published

The non-equilibrium dynamics of electrons and nuclei govern the function of photoactive materials. Disentangling these dynamics remains a critical goal for understanding photoactive materials. Here we investigate the photoinduced dynamics of the [Fe(bmip)2]2+ photosensitizer, where bmip = 2,6-bis(3-methyl-imidazole-1-ylidine)-pyridine, with simultaneous femtosecond-resolution Fe Kα and Kβ X-ray emission spectroscopy (XES) and X-ray solution scattering (XSS). This measurement shows temporal oscillations in the XES and XSS difference signals with the same 278 fs period oscillation. These oscillations originate from an Fe-ligand stretching vibrational wavepacket on a triplet metal-centered (3MC) excited state surface. This 3MC state is populated with a 110 fs time constant by 40% of the excited molecules while the rest relax to a 3MLCT excited state. The sensitivity of the Kα XES to molecular structure results from a 0.7% average Fe-ligand bond length shift between the 1 s and 2p core-ionized states surfaces., Photoinduced non-adiabatic intramolecular processes have important applications but their mechanisms are challenging to explore. Here the authors detect and assign vibrational wavepacket dynamics in a Fe carbene complex by ultrafast X-ray emission spectroscopy and X-ray scattering, resolving nuclear and electronic motion.

Post-calculation analyses are often required to extract physical insights from ab initio molecular dynamics simulations. In the present work, we use different machine learning classifiers to take a new perspective on the decomposition reaction of dioxetane. Upon thermally activated decomposition, dioxetane can form products in an electronically excited state and can thus chemiluminesce. Simulated dynamics trajectories exhibit both successful and frustrated dissociations. As an exhaustive and systematic study of the decomposition mechanism “by hand” is beyond feasibility, machine learning models have been employed to study the relevant nuclear distortions governing molecular dissociation. According to all classifiers used in the study, the two sets of geometries differ by the in-phase planarisation of the two formaldehyde moieties. New insights are obtained from this analysis: if both moieties are not planar enough when the dissociation is attempted, it is frustrated and the molecule remains trapped. The postponing of the decomposition reaction by the so-called entropic trap enhances the chemiexcitation efficiency.

Machine learning techniques, specifically Gradient-Enhanced Kriging (GEK), has been implemented for molecular geometry optimization.GEK has many advantages as compared to conventional -- step-restricted second-order truncated -- molecular optimization methods.In particular, the surrogate model associated with GEK can have multiple stationary points, will smoothly converge to theexact model as the size of the data set increases, and contains an explicit expression for the expected average error of the model functionat an arbitrary point in space.In this respect GEK can be of interest for methods used in molecular geometry optimizations.GEK is usually, however, associated with abundance of data, contrary to the situation desired forefficient geometry optimizations.In the paper we will demonstrate how the GEK procedure can be utilized in a fashion such that in the presence of few data points, thesurrogate surface will in a robust way guide the optimization to a minimum of a molecular structure.In this respect the GEK procedure will be used to mimic the behavior of a conventional second-order scheme, but retaining theflexibility of the superior machine learning approach -- GEK is an exact interpolator.Moreover, the expected variance will be used in the optimization to facilitate restricted-variance rational function optimizations (RV-RFO).A procedure which relates the eigenvalues of the Hessian-model-function Hessian with the individual characteristiclengths, used in the GEK, reduces the number of empirical parameters to two -- the value of the trend function and themaximum allowed variance. These parameters are determined using the extended Baker (e-Baker) test suite, at the Hartree-Fock level of approximation,and a single reaction of the Baker transition-state (Baker-TS) test suite as a training set. The so-created optimizationprocedure -- RV-RFO-GEK -- is tested using the e-Baker, the full Baker-TS, and the S22 test suites, at the density-functional-theory level for the two Baker test suitesand at the second order Møller-Plesset level of approximation for the S22 test suite, respectively.The tests show that the new method is generally on par with a state-of-the-art conventional method, while for difficult cases it exhibits a definite advantage.

International audience; Recent work, particularly by Cederbaum and co-workers, has identified the phenomenon of charge migration, whereby charge flow occurs over a static molecular framework after the creation of an electronic wavepacket. In a real molecule, this charge migration competes with charge transfer, whereby the nuclear motion also results in the re-distribution of charge. To study this competition, quantum dynamics simulations need to be performed. To break the exponential scaling of standard grid-based algorithms, approximate methods need to be developed that are efficient yet able to follow the coupled electronic-nuclear motion of these systems. Using a simple model Hamiltonian based on the ionisation of the allene molecule, the performance of different methods based on Gaussian Wavepackets is demonstrated.

In this article we describe the OpenMolcas environment and invite the computational chemistry community to collaborate. The open-source project already includes a large number of new developments realized during the transition from the commercial MOLCAS product to the open-source platform. The paper initially describes the technical details of the new software development platform. This is followed by brief presentations of many new methods, implementations, and features of the OpenMolcas program suite. These developments include novel wave function methods such as stochastic complete active space self-consistent field, density matrix renormalization group (DMRG) methods, and hybrid multiconfigurational wave function and density functional theory models. Some of these implementations include an array of additional options and functionalities. The paper proceeds and describes developments related to explorations of potential energy surfaces. Here we present methods for the optimization of conical intersections, the simulation of adiabatic and nonadiabatic molecular dynamics and interfaces to tools for semiclassical and quantum mechanical nuclear dynamics. Furthermore, the article describes features unique to simulations of spectroscopic and magnetic phenomena such as the exact semiclassical description of the interaction between light and matter, various X-ray processes, magnetic circular dichroism and properties. Finally, the paper describes a number of built-in and add-on features to support the OpenMolcas platform with post calculation analysis and visualization, and new electronic and muonic basis sets.

International audience; Electronically excited states play important roles in many chemical reactions and spectroscopic techniques. In quantum chemistry, a common technique to solve excited states is the multiroot Davidson algorithm, but it is not designed for processes like X-ray spectroscopy that involves hundreds of highly excited states. We show how the use of a restricted active space wavefunction together with a projection operator to remove low-lying electronic states offers an efficient way to reach single and double-core-hole states. Additionally, several improvements to the stability and efficiency of the configuration interaction (CI) algorithm for a large number of states are suggested. When applied to a series of transition metal complexes the new CI algorithm does not only resolve divergence issues but also leads to typical reduction in computational time by 70%, with the largest savings for small molecules and large active spaces. Together, the projection operator and the improved CI algorithm now make it possible to simulate a wide range of single- and two-photon spectroscopies. © 2019 Wiley Periodicals, Inc.

In this article we describe the OpenMolcas environment and invite the computational chemistry community to collaborate. The open-source project already includes a large number of new developments realized during the transition from the commercial MOLCAS product to the open-source platform. The paper initially describes the technical details of the new software development platform. This is followed by brief presentations of many new methods, implementations, and features of the OpenMolcas program suite. These developments include novel wave function methods such as stochastic complete active space self-consistent field, density matrix renormalization group (DMRG) methods, and hybrid multiconfigurational wave function and density functional theory models. Some of these implementations include an array of additional options and functionalities. The paper proceeds and describes developments related to explorations of potential energy surfaces. Here we present methods for the optimization of conical intersections, the simulation of adiabatic and nonadiabatic molecular dynamics, and interfaces to tools for semiclassical and quantum mechanical nuclear dynamics. Furthermore, the article describes features unique to simulations of spectroscopic and magnetic phenomena such as the exact semiclassical description of the interaction between light and matter, various X-ray processes, magnetic circular dichroism, and properties. Finally, the paper describes a number of built-in and add-on features to support the OpenMolcas platform with postcalculation analysis and visualization, a multiscale simulation option using frozen-density embedding theory, and new electronic and muonic basis sets.

International audience; Molecular dynamics simulations are often key to the understanding of the mechanism, rate and yield of chemical reactions. One current challenge is the in-depth analysis of the large amount of data produced by the simulations, in order to produce valuable insight and general trends. In the present study, we propose to employ recent machine learning analysis tools to extract relevant information from simulation data without a priori knowledge on chemical reactions. This is demonstrated by training machine learning models to predict directly a specific outcome quantity of ab initio molecular dynamics simulations - the timescale of the decomposition of 1,2-dioxetane. The machine learning models accurately reproduce the dissociation time of the compound. Keeping the aim of gaining physical insight, it is demonstrated that, in order to make accurate predictions, the models evidence empirical rules that are, today, part of the common chemical knowledge. This opens the way for conceptual breakthroughs in chemistry where machine analysis would provide a source of inspiration to humans.

International audience

International audience; The influence of ring-puckering on the light-induced ring-opening dynamics of heterocyclic compounds was studied on the sample 5-membered ring molecules γ-valerolactone and 5H-furan-2-one using time-resolved photoelectron spectroscopy and ab initio molecular dynamics simulations. In γ-valerolactone, ring-puckering is not a viable relaxation channel and the only available reaction pathway is ring-opening, which occurs within one vibrational period along the C-O bond. In 5H-furan-2-one, the C=C double bond in the ring allows for ring-puckering which slows down the ring-opening process by about 150 fs while only marginally reducing its quantum yield. This demonstrates that ring-puckering is an ultrafast process, which is directly accessible upon excitation and which spreads the excited state wave packet quickly enough to influence even the outcome of an otherwise expectedly direct ring-opening reaction.

Machine learning models, trained to reproduce molecular dynamics results, help interpreting simulations and extracting new understanding of chemistry., Molecular dynamics simulations are often key to the understanding of the mechanism, rate and yield of chemical reactions. One current challenge is the in-depth analysis of the large amount of data produced by the simulations, in order to produce valuable insight and general trends. In the present study, we propose to employ recent machine learning analysis tools to extract relevant information from simulation data without a priori knowledge on chemical reactions. This is demonstrated by training machine learning models to predict directly a specific outcome quantity of ab initio molecular dynamics simulations – the timescale of the decomposition of 1,2-dioxetane. The machine learning models accurately reproduce the dissociation time of the compound. Keeping the aim of gaining physical insight, it is demonstrated that, in order to make accurate predictions, the models evidence empirical rules that are, today, part of the common chemical knowledge. This opens the way for conceptual breakthroughs in chemistry where machine analysis would provide a source of inspiration to humans.

International audience; An algorithm is described for quantum dynamics where an Ehrenfest potential is combined with fully quantum nuclear motion (Quantum-Ehrenfest, Qu-Eh). The method is related to the single-set variational multi-configuration Gaussian approach (vMCG) but has the advantage that only a single quantum chemistry computation is required at each time step since there is only a single time-dependent potential surface. Also shown is the close relationship to the "exact factorization method." The quantum Ehrenfest method is compared with vMCG for study of electron dynamics in a modified bismethylene-adamantane cation system. Illustrative examples of electron-nuclear dynamics are presented for a distorted allene system and for HCCI+ where one has a degenerate Π system.

Knowledge about the electron dynamics in molecules is essential for our understanding of chemical and biological processes. Because of their light mass, electrons are expected to move on the attosecond (1 as = 10− 18 s) timescale. The first synthesis of attosecond pulses in 2001 has opened up the possibility of probing electronic motion on its intrinsic timescale. Excitation or ionisation of a molecule with such a short pulse leads to the coherent population of several electronic states, called an electronic wavepacket. The interference between electronic states in such a superposition, alternating between constructive and destructive, leads to oscillating motion of the electron cloud. This purely quantum process relies on the coherence of the electronic wavepacket. A fundamental challenge is to understand to what extent the electronic wavepacket retains its coherence, i.e., how long the oscillations in the electron cloud survive, in the presence of interactions with the nuclei of the molecule. To address this question, we have developed semi-classical and quantum mechanical methods to simulate the dynamics upon ionisation of polyatomic molecules. The chapter contains a review of the theoretical methods we have developed and some applications illustrating new important physical insights about the predicted decoherence process.

We demonstrate that charge migration can be ‘engineered’ in arbitrary molecular systems if a single localised orbital – that diabatically follows nuclear displacements – is ionised. Specifically, we describe the use of natural bonding orbitals in Complete Active Space Configuration Interaction (CASCI) calculations to form cationic states with localised charge, providing consistently well-defined initial conditions across a zero point energy vibrational ensemble of molecular geometries. In Ehrenfest dynamics simulations following localised ionisation of -electrons in model polyenes (hexatriene and decapentaene) and -electrons in glycine, oscillatory charge migration can be observed for several femtoseconds before dephasing. Including nuclear motion leads to slower dephasing compared to fixed-geometry electron-only dynamics results. For future work, we discuss the possibility of designing laser pulses that would lead to charge migration that is experimentally observable, based on the proposed diabatic orbital approach. KEYWORDS: Ehrenfest method, coupled electron-nuclear dynamics, charge migration, localised orbital

Bioluminescence is a phenomenon that has fascinated mankind for centuries. Today the phenomenon and its sibling, chemiluminescence, have impacted society with a number of useful applications in fields like analytical chemistry and medicine, just to mention two. In this review, a molecular-orbital perspective is adopted to explain the chemistry behind chemiexcitation in both chemi- and bioluminescence. First, the uncatalyzed thermal dissociation of 1,2-dioxetane is presented and analyzed to explain, for example, the preference for triplet excited product states and increased yield with larger nonreactive substituents. The catalyzed fragmentation reaction and related details are then exemplified with substituted 1,2-dioxetanone species. In particular, the preference for singlet excited product states in that case is explained. The review also examines the diversity of specific solutions both in Nature and in artificial systems and the difficulties in identifying the emitting species and unraveling the color modulation process. The related subject of excited-state chemistry without light absorption is finally discussed. The content of this review should be an inspiration to human design of new molecular systems expressing unique light-emitting properties. An appendix describing the state-of-the-art experimental and theoretical methods used to study the phenomena serves as a complement.

International audience; Chemiluminescence is the emission of light as a result of a nonadiabatic chemical reaction. The present work is concerned with understanding the yield of chemiluminescence, in particular how it dramatically increases upon methylation of 1,2-dioxetane. Both ground-state and nonadiabatic dynamics (including singlet excited states) of the decomposition reaction of various methyl-substituted dioxetanes have been simulated. Methyl-substitution leads to a significant increase in the dissociation time scale. The rotation around the O-C-C-O dihedral angle is slowed; thus, the molecular system stays longer in the "entropic trap" region. A simple kinetic model is proposed to explain how this leads to a higher chemiluminescence yield. These results have important implications for the design of efficient chemiluminescent systems in medical, environmental, and industrial applications.

International audience; We review various aspects of photoemission dynamics in the case of two-photon ionisation. We first recall the definition of a transition phase specific to two-photon transitions. Numerical experiments on model atoms are used to show how the group delay associated with the transition phase is actually representative of the early dynamics of the detected photoelectron wave packets. Then we address the question of measuring these transition delays using a standard interferometric technique of experimental attosecond physics, so-called RABBIT. Finally, we outline different reinterpretations of RABBITgiving access to the more fundamental scattering dynamics affecting any photoemission processes.

International audience

International audience; We theoretically study high-harmonic generation in toluene, ortho-xylene and fluorobenzene driven by a 1.8 μm ultrashort pulse. We find that the chemical substitutions have a strong influence on the amplitude and phase of the emission from the highest occupied molecular orbital, despite having a small influence on the orbital itself. We show that this influence is due to the tunnel ionization step, which depends critically on the sign and amplitude of the asymptotic part of the wave function. We discuss how these effects would manifest in phase-sensitive high-harmonic generation spectroscopy experiments.

International audience; We have simulated the coupled electron and nuclear dynamics using the Ehrenfest method upon valence ionisation of modified bismethylene-adamantane (BMA) molecules where there is an electron transfer between the two π bonds. We have shown that the nuclear motion significantly affects the electron dynamics after a few fs when the electronic states involved are close in energy. We have also demonstrated how the non-stationary electronic wave packet determines the nuclear motion, more precisely the asymmetric stretching of the two π bonds, illustrating "charge-directed reactivity". Taking into account the nuclear wave packet width results in the dephasing of electron dynamics with a half-life of 8 fs; this eventually leads to the equal delocalisation of the hole density over the two methylene groups and thus symmetric bond lengths.

International audience; In this article, we outline the current state-of-theart “on-the-fly” methods for non-adiabatic dynamics, highlighting the similarities and differences between them. We derive the equations of motion for both the Ehrenfest and variational multi-configuration Gaussian (vMCG) methods from the Dirac–Frenkel variational principle. We explore the connections between these two methods by presenting an alternative derivation of the vMCG method, which gives the Ehrenfest equations of motion when taking the appropriate limits.

International audience; We simulate electron dynamics following ionization in 2-phenyl-ethyl-amine and 2-phenylethyl-N,N-dimethylamine as examples of systems where 3 coupled cationic states are involved. We study two nuclear effects on electron dynamics: (i) coupled electron-nuclear motion and (ii) nuclear spatial delocalization as a result of the zero-point energy in the neutral molecule. Within the Ehrenfest approximation, our calculations show that the coherent electron dynamics in these molecules is not lost as a result of coupled electron-nuclear motion. In contrast, as a result of nuclear spatial delocalization, dephasing of the oscillations occurs on a time scale of only a few fs, long before any significant nuclear motion can occur. The results have been rationalized using a semi-quantitative model based upon the gradients of the potential energy surfaces.

The dynamics of photoemission is fully encoded in the continuum wave functions selected by the transitions. Using numerical simulations on simple benchmark models, we show how scattering phase shifts and photoemission delays can be retrieved from this unambiguously defined class of wave functions. In contrast with standard scattering waves inherited from collision theory, they are real-valued for one-photon transitions and provide a clear-cut interpretation of the delays recently discussed in the framework of attosecond science.

International audience; The advent of attosecond techniques opens up the possibility to observe experimentally electron dynamics following ionization of molecules. Theoretical studies of pure electron dynamics at single fixed nuclear geometries in molecules have demonstrated oscillatory charge migration at a well-defined frequency but often neglecting the natural width of the nuclear wave packet. The effect on electron dynamics of the spatial delocalization of the nuclei is an outstanding question. Here, we show how the inherent distribution of nuclear geometries leads to dephasing. Using a simple analytical model, we demonstrate that the conditions for a long-lived electronic coherence are a narrow nuclear wave packet and almost parallel potential-energy surfaces of the states involved. We demonstrate with numerical simulations the decoherence of electron dynamics for two real molecular systems (paraxylene and polycyclic norbornadiene), which exhibit different decoherence time scales. To represent the quantum distribution of geometries of the nuclear wave packet, the Wigner distribution function is used. The electron dynamics decoherence result has significant implications for the interpretation of attosecond spectroscopy experiments since one no longer expects long-lived oscillations.

The material presented in this chapter is based on important advances realized in “attophysics” which make feasible to follow the motion of electrons in atoms and molecules with attosecond-level time resolution. In this context, time-delays have been recently determined in the process of photoionization by extreme-ultra-violet (XUV) pulses and the question of the significance of these measured delays arises. As we shall outline here, numerical experiments show that they are intimately related to the structure of the ionized species’ continuous spectrum. Another point addressed here is that, in experiments, the measurements have the common characteristic to be performed in the presence of an auxiliary infra-red (IR) field, used to “clock” the timing of the process. This implies to adapt the theory treatment to handle such “two-color” photoionization processes. We review a systematic analysis of these features that are characteristic of this class of electronic transitions, when viewed in the time domain.

International audience; Nonadiabatic dynamics in the vicinity of conical intersections is of essential importance in photochemistry. It is well known that if the branching space is represented in polar coordinates, then for a geometry represented by angle θ, the corresponding adiabatic states are obtained from the diabatic states with the mixing angle θ/2. In an equivalent way, one can study the relation between the real rotation of diabatic states and the resulting nuclear gradient. In this work, we extend the concept to allow a complex rotation of diabatic states to form a nonstationary superposition of electronic states. Our main result is that this leads to an elliptical transformation of the effective potential energy surfaces; i.e., the magnitude of the initial nuclear gradient changes as well as its direction. We fully explore gradient changes that result from varying both θ and ϕ (the complex rotation angle) as a way of electronically controlling nuclear motion, through Ehrenfest dynamics simulations for benzene cation.

International audience; Photoionization can create a nonstationary electronic state and therefore initiates coupled electron-nuclear dynamics in molecules. Using a CASSCF implementation of the Ehrenfest method, we study the nuclear dynamics following vertical ionization of toluene, starting close to the conical intersection between ground and first excited states of its cation. The results show how the initial nuclear dynamics is controlled by the nonstationary electronic state character. In particular, ionization of this system leading to an equal superposition of the two lowest energy states can initiate nuclear dynamics in an orthogonal direction in the branching space to dynamics on the ground or first excited state potential energy surfaces alone.

This article describes the Ehrenfest method and our second-order implementation (with approximate gradient and Hessian) within a CASSCF formalism. We demonstrate that the second-order implementation with the predictor–corrector integration method improves the accuracy of the simulation significantly in terms of energy conservation. Although the method is general and can be used to study any coupled electron–nuclear dynamics, we apply it to investigate charge migration upon ionization of small organic molecules, focusing on benzene cation. Using this approach, we can study the evolution of a nonstationary electronic wavefunction for fixed atomic nuclei, and where the nuclei are allowed to move, to investigate the interplay between them for the first time. Analysis methods for the interpretation of the electronic and nuclear dynamics are suggested: we monitor the electronic dynamics by calculating the spin density of the system as a function of time.

Coupled electron-nuclear dynamics, implemented using the Ehrenfest method, has been used to study charge migration with fixed nuclei, together with charge transfer when nuclei are allowed to move. Simulations were initiated at reference geometries of neutral benzene and 2-phenylethylamine (PEA), and at geometries close to potential energy surface crossings in the cations. Cationic eigenstates, and the so-called sudden approximation, involving removal of an electron from a correlated ground-state wavefunction for the neutral species, were used as initial conditions. Charge migration without coupled nuclear motion could be observed if the Ehrenfest simulation, using the sudden approximation, was started near a conical intersection where the states were both strongly coupled and quasi-degenerate. Further, the main features associated with charge migration were still recognizable when the nuclear motion was allowed to couple. In the benzene radical cation, starting from the reference neutral geometry with the sudden approximation, one could observe sub-femtosecond charge migration with a small amplitude, which results from weak interaction with higher electronic states. However, we were able to engineer large amplitude charge migration, with a period between 10 and 100 fs, corresponding to oscillation of the electronic structure between the quinoid and anti-quinoid cationic electronic configurations, by distorting the geometry along the derivative coupling vector from the D6h Jahn-Teller crossing to lower symmetry where the states are not degenerate. When the nuclear motion becomes coupled, the period changes only slightly. In PEA, in an Ehrenfest trajectory starting from the D2 eigenstate and reference geometry, a partial charge transfer occurs after about 12 fs near the first crossing between D1, D2 (N(+)-Phenyl, N-Phenyl(+)). If the Ehrenfest propagation is started near this point, using the sudden approximation without coupled nuclear motion, one observes an oscillation of the spin density--charge migration--between the N atom and the phenyl ring with a period of 4 fs. When the nuclear motion becomes coupled, this oscillation persists in a damped form, followed by an effective charge transfer after 30 fs.