Timing the recoherences of attosecond electronic charge migration by quantum control of femtosecond nuclear dynamics: A case study for HCCI. J Chem Phys 2019;151:244306. [PMID: 31893866 DOI: 10.1063/1.5134665]

Reducing the Cost of TD-CI Simulations of Strong Field Ionization

Strong field ionization of molecules by intense laser pulses can be simulated by time-dependent configuration interaction (TD-CI) with a complex absorbing potential (CAP). Standard molecular basis sets need to be augmented with several sets of diffuse functions for effective interaction with the CAP. This dramatically increases the number of configurations and the cost of the TD-CI simulations as the size of the molecules increases. The cost can be reduced by making use of spin symmetry and by employing an orbital energy cutoff to limit the number of virtual orbitals used to construct the excited configurations. Greater reductions in the number of virtual orbitals can be obtained by examining their interaction with the absorbing potential during simulations and their contributions to the strong field ionization rate. This can be determined from the matrix elements of the absorbing potential and the TD-CI coefficients from test simulations. Compared to a simple 3 hartree cutoff in the orbital energies, these approaches reduce the number of virtual orbitals by 20-35% for neutral molecules and 5-10% for cations. As a result, the cost of simulations is reduced by 35-60% for neutral molecules. The number of virtual orbitals needed can also be estimated by second-order perturbation theory without the need for test simulations. The number of virtual orbitals can be reduced further by adapting orbitals to the laser field using natural orbitals derived from test simulations. This is particularly effective for cations, yielding reductions of more than 20%.

Core-Hole Coherent Spectroscopy in Molecules

We study the ultrafast dynamics initiated by a coherent superposition of core-excited states of nitrous oxide molecule. Using high-level ab initio methods, we show that the decoherence caused by the electronic decay and the nuclear dynamics is substantially slower than the induced ultrafast quantum beatings, allowing the system to undergo several oscillations before it dephases. We propose a proof-of-concept experiment using the harmonic up-conversion scheme available at x-ray free-electron laser facilities to trace the evolution of the created core-excited-state coherence through a time-resolved x-ray photoelectron spectroscopy.

Size Effect in Correlation-Driven Charge Migration in Correlation Bands of Alkyne Chains

Correlation-driven charge migration initiated by inner-valence ionization leading to the population of the correlation bands of alkyne chains containing between 4 and 12 carbon atoms is explored through ab initio simulations. Scaling laws are observed, both for the time scale of the charge migration and for the slope of the density of states of the correlation bands. These can be used for predicting the relaxation time scale in much larger systems from the same molecular family and for finding promising candidates for the development of an attochemistry scheme taking advantages of the specificity of the dynamics in the correlation bands of molecules.

Spin current in the early stage of radical reactions and its mechanisms

We study the electronic spin flux (atomic-scale flow of the spin density in molecules) by a perturbation analysis and ab initio nonadiabatic calculations. We derive a general perturbative expression of the charge and spin fluxes and identify the driving perturbation of the fluxes to be the time derivative of the electron-nucleus interaction term in the Hamiltonian. We then expand the expression in molecular orbitals so as to identify relevant components of the fluxes. Our perturbation theory describes the electronic fluxes in the early stage of reactions in an intuitively clear manner. The perturbation theory is then applied to an analysis of the spin flux obtained in ab initio calculations of the radical reaction of O2 and CH3· starting from three distinct spin configurations; (a) CH3· and triplet O2 with total spin of the system set Stot=1/2 (b) CH3· and singlet O2, Stot=1/2, and (c) CH3· and triplet O2, Stot=3/2. Further analysis of the time-dependent behaviors of the spin flux in these numerical simulations reveals (i) the spin flux induces rearrangement of the local spin structure, such as reduction of the spin polarization arising from the triplet O2 and (ii) the spin flux flows from O2 to CH3· in the reaction starting from spin configuration (a) and from CH3· to O2 in that starting from configuration (b), whereas no major intermolecular spin flux was observed in that starting from configuration (c). Our study thus establishes the mechanism of the spin flux that rearranges the local spin structures associated with chemical bonds.

Charge Migration in HCCI Cations Probed by Strong Field Ionization: Time-Dependent Configuration Interaction and Vibrational Wavepacket Simulations

Strong field ionization of neutral iodoacetylene (HCCI) can produce a coherent superposition of the X and A cations and results in charge migration between the CC π orbital and the iodine π-type lone pair. This charge migration causes oscillations in the rate of strong field ionization of the cation to the dication that can be monitored using intense few-cycle probe pulses. The dynamics and strong field ionization of the coherent superposition the X and A states of HCCI⁺ have been modeled by time-dependent configuration interaction (TDCI) simulations. When the nuclei are allowed to move, the electronic wavefunctions need to be multiplied by vibrational wavefunctions. Nuclear motion has been modeled by vibrational packets moving on quadratic approximations to the potential energy surfaces for the X and A states of the cation. The overlap of the vibrational wavepackets decays in about 10-15 fs. Consequently, the oscillations in the strong field ionization decay on the same time scale. A revival of the vibrational overlap and in the oscillations of the strong field ionization is seen at 60-110 fs. TDCI simulations show that the decay and revival of the charge migration can be monitored by strong field ionization with intense 2- and 4-cycle linearly polarized 800 nm pulses. The revival is also seen with 7-cycle pulses.

Probing Electronic Symmetry Reduction during Charge Migration via Time-Resolved X-Ray Diffraction

The present work focuses on probing ultrafast charge migration after symmetry-breaking excitation using ultrashort laser pulses. LiCN is chosen as prototypical system because it can be oriented in the laboratory frame and it possesses optically-accessible charge transfer states at low energies. The charge migration is simulated within the hybrid time-dependent density functional theory/configuration interaction framework. Time-resolved electronic current densities and simulated time-resolved x-ray diffraction signals are used to unravel the mechanism of charge migration. Our simulations demonstrate that specific choices of laser polarization lead to a control over the symmetry of the induced charge migration. Moreover, time-resolved x-ray diffraction signals are shown to encode transient symmetry reduction at intermediate times.

Time-resolved X-ray and XUV based spectroscopic methods for nonadiabatic processes in photochemistry

The photochemistry of numerous molecular systems is influenced by conical intersections (CIs). These omnipresent nonadiabatic phenomena provide ultra-fast radiationless relaxation channels by creating degeneracies between electronic states and decide over the final photoproducts. In their presence, the Born-Oppenheimer approximation breaks down, and the timescales of the electron and nuclear dynamics become comparable. Due to the ultra-fast dynamics and the complex interplay between nuclear and electronic degrees of freedom, the direct experimental observation of nonadiabatic processes close to CIs remains challenging. In this article, we give a theoretical perspective on novel spectroscopic techniques capable of observing clear signatures of CIs. We discuss methods that are based on ultra-short laser pulses in the extreme ultraviolet and X-ray regime, as their spectral and temporal resolution allow for resolving the ultra-fast dynamics near CIs.

Ionization of HCCI Neutral and Cations by Strong Laser Fields Simulated With Time Dependent Configuration Interaction

Strong field ionization of neutral iodoacetylene (HCCI) can produce a coherent superposition of the X and A cations. This superposition results in charge migration between the CC π orbital and the iodine π-type lone pair which can be monitored by strong field ionization with short, intense probe pulses. Strong field ionization of the X and A states of HCCI cation was simulated with time-dependent configuration interaction using singly ionized configurations and singly excited, singly ionized configurations (TD-CISD-IP) and an absorbing boundary. Studies with static fields were used to obtain the 3-dimensional angular dependence of instantaneous ionization rates by strong fields and the orbitals involved in producing the cations and dications. The frequency of charge oscillation is determined by the energy separation of the X and A states; this separation can change depending on the direction and strength of the field. Furthermore, fields along the molecular axis can cause extensive mixing between the field-free X and A configurations. For coherent superpositions of the X and A states, the charge oscillations are characterized by two frequencies-the driving frequency of the laser field of the probe pulse and the intrinsic frequency due to the energy separation between the X and A states. For linear and circularly polarized pulses, the ionization rates show marked differences that depend on the polarization direction of the pulse, the carrier envelope phase and initial phase of the superposition. Varying the initial phase of the superposition at the beginning of the probe pulse is analogous to changing the delay between the pump and probe pulses. The charge oscillation in the coherent superposition of the X and A states results in maxima and minima in the ionization yield as a function of the superposition phase.

Systematic Investigation of the Reliability of the Frozen Nuclei Approximation for Short-Pulse Excitation: The Example of HCCI+

In this work we quantitatively study the reliability of the frozen nuclei approximation for ultrafast dynamics. Specifically we study laser excitation of HCCI⁺ from its ground state to the first electronically excited state. The population of the first excited state is obtained by both the frozen nuclei approximation and by multidimensional nuclear dynamics. Detailed comparison of the results by the two methods are performed to provide quantitative criteria for the reliability of the frozen nuclei approximation for this system.

Transient Valence Charge Localization in Strong-Field Dissociative Ionization of HCl Molecules

Probing transient charge localization in the innershell orbital of atoms and molecules has been made possible by the recent progress of advanced light sources. Here, we demonstrate that the ultrafast electron tunneling ionization by an intense femtosecond laser pulse could induce an asymmetric transient charge localization in the valence shell of the HCl molecule during the dissociative ionization process. The transient charge localization is encoded in the laser impulse acquired by the outgoing ionic fragments, and the asymmetry is revealed by carefully examining the electron tunneling-site distinguished momentum angular distribution of the ejected H^{+} fragments. Our work proposes a way to visualize the transient valence charge motion and will stimulate further investigations of the tunneling-site-sensitive ultrafast dynamics of molecules in strong laser fields.

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.

Theoretical Study of Dynamic Stark-Induced π-Electron Rotations in Low-Symmetry Aromatic Ring Molecules beyond the Frozen Nuclear Approximation

The effects of vibrational motions on dynamic Stark-induced π-electron rotations in a low-symmetry aromatic ring molecule are theoretically studied in the adiabatic approximation. We adopt a simplified three-electronic state model with a few vibronic states. A pair of the lowest vibronic states in two electronic excited states is set degenerate by irradiation of two linearly polarized UV lasers. The resultant degenerate state is named the dynamic Stark-induced degenerate vibronic state (DSIDVS). The laser parameters (intensities and central frequencies) are determined under the conditions of DSIDVS formation. The aromatic ring molecules of interest are supposed to belong to the weak coupling case. The analytical expressions for the DSIDVS and coherent angular momentum L_Z(t) are derived in the displaced harmonic oscillator (DHO) model. Two horizontal potential displacements (δ_α, δ_β) between the two electronic excited states (α and β) and the ground state are the parameters in the DHO model. The L_Z(t) calculated with δ_α = δ_β is characterized by a regular sequence of the angular momentum pulses with a positive (or negative) constant. For a more general case with δ_α ≠ δ_β, the regular sequence is broken down because of the contribution of the first excited vibronic state in each electronic state to L_Z(t).