Full Ab Initio Many-Electron Simulation of Attosecond Molecular Pump–Probe Spectroscopy

Multidimensional high-harmonic echo spectroscopy: Resolving coherent electron dynamics in the EUV regime

We theoretically propose a multidimensional high-harmonic echo spectroscopy technique which utilizes strong optical fields to resolve coherent electron dynamics spanning an energy range of multiple electronvolts. Using our recently developed semi-perturbative approach, we can describe the coherent valence electron dynamics driven by a sequence of phase-matched and well-separated short few-cycle strong infrared laser pulses. The recombination of tunnel-ionized electrons by each pulse coherently populates the valence states of a molecule, which allows for a direct observation of its dynamics via the high harmonic echo signal. The broad bandwidth of the effective dipole between valence states originated from the strong-field excitation results in nontrivial ultra-delayed partial rephasing echo, which is not observed in standard two-dimensional optical spectroscopic techniques in a two-level molecular systems. We demonstrate the results of simulations for the anionic molecular system and show that the ultrafast valence electron dynamics can be well captured with femtosecond resolution.

Nonadiabatic dynamics with classical trajectories: The problem of an initial coherent superposition of electronic states

Advances in coherent light sources and development of pump-probe techniques in recent decades have opened the way to study electronic motion in its natural time scale. When an ultrashort laser pulse interacts with a molecular target, a coherent superposition of electronic states is created and the triggered electron dynamics is coupled to the nuclear motion. A natural and computationally efficient choice to simulate this correlated dynamics is a trajectory-based method where the quantum-mechanical electronic evolution is coupled to a classical-like nuclear dynamics. These methods must approximate the initial correlated electron-nuclear state by associating an initial electronic wavefunction to each classical trajectory in the ensemble. Different possibilities exist that reproduce the initial populations of the exact molecular wavefunction when represented in a basis. We show that different choices yield different dynamics and explore the effect of this choice in Ehrenfest, surface hopping, and exact-factorization-based coupled-trajectory schemes in a one-dimensional two-electronic-state model system that can be solved numerically exactly. This work aims to clarify the problems that standard trajectory-based techniques might have when a coherent superposition of electronic states is created to initialize the dynamics, to discuss what properties and observables are affected by different choices of electronic initial conditions and to point out the importance of quantum-momentum-induced electronic transitions in coupled-trajectory schemes.

Irradiation of Plutonium Tributyl Phosphate Complexes by Ionizing Alpha Particles: A Computational Study

The PUREX solvent extraction process, widely used for recovering uranium and plutonium from spent nuclear fuel, utilizes an organic solvent composed of tributyl phosphate (TBP). The emission of ionizing particles such as alpha particles, resulting from the decay of plutonium, makes the organic solvent vulnerable to degradation. Here, we study the ultrashort time alpha irradiation of tributylphosphate (TBP) and Pu(NO3)4(TBP)2 complex formed in the PUREX process. Electron dynamics is propagated by Real-Time-Dependent Auxiliary Density Functional Theory (RT-TD-ADFT). We investigate the use of previously proposed absorption boundary conditions (ABC) in the molecular orbital space to treat secondary electron emission. Basis set and exchange correlation functional effects with ABC are reported as well as a detailed analysis of the ABC parametrization. Preliminary results on the water molecule and then on TBP show that the phenomenological nature of the ABC parameters necessitates selecting appropriate values for each system under study. Irradiation of free and complexed TBP shows an influence of the ligands on the variation of atomic charges on the femtosecond time scale. An accumulation of atomic charges in the alkyl chains of TBP is observed in the case where the nitrate groups are predominantly irradiated. In addition, we find that the Pu atom regains its electric charge very rapidly after being hit by the projectile, with the coordination sphere serving as an electron reservoir to preserve its formal redox state. This study paves the road toward a full understanding of the degradation of organic extracants employed in the nuclear industry.

Quantum coherence in molecular photoionization

The study of onset and decay, as well as control of ultrafast quantum coherence in many-electron systems is in the focus of interest of attosecond physics. Interpretation of attosecond experiments detecting the ultrafast quantum coherence requires application of advanced theoretical and computational tools combining many-electron theory, description of the electronic continuum, including in the strong laser field scenario, as well as nuclear dynamics theory. This perspective reviews the recent theoretical advances in understanding the attosecond dynamics of quantum coherence in photoionized molecular systems and outlines possible future directions of theoretical and experimental study of coherence and entanglement in the attosecond regime.

Dialogue on analytical and ab initio methods in attoscience

The perceived dichotomy between analytical and ab initio approaches to theory in attosecond science is often seen as a source of tension and misconceptions. This Topical Review compiles the discussions held during a round-table panel at the 'Quantum Battles in Attoscience' cecam virtual workshop, to explore the sources of tension and attempt to dispel them. We survey the main theoretical tools of attoscience-covering both analytical and numerical methods-and we examine common misconceptions, including the relationship between ab initio approaches and the broader numerical methods, as well as the role of numerical methods in 'analytical' techniques. We also evaluate the relative advantages and disadvantages of analytical as well as numerical and ab initio methods, together with their role in scientific discovery, told through the case studies of two representative attosecond processes: non-sequential double ionisation and resonant high-harmonic generation. We present the discussion in the form of a dialogue between two hypothetical theoreticians, a numericist and an analytician, who introduce and challenge the broader opinions expressed in the attoscience community.

Quantum electronic coherences by attosecond transient absorption spectroscopy: ab initio B-spline RCS-ADC study

Here I present a fully ab initio time-resolved study of X-ray attosecond transient absorption spectroscopy (ATAS) in a prototypical polyatomic molecule, pyrazine, and demonstrate the possibility of retrieving the many-electron quantum ionic coherences arising in attosecond molecular photoionisation and pre-determining the subsequent charge-directed photochemical reactivity. Advanced first-principles many-electron simulations are performed, within a hybrid XUV pump/X-ray probe setup, to describe the interaction of pyrazine with both XUV pump and X-ray probe pulses, and study the triggered correlated many-electron dynamics. The calculations are carried out by means of the recently-developed ab initio method for many-electron dynamics in polyatomic molecules, the time-dependent (TD) B-spline Restricted Correlation Space-Algebraic Diagrammatic Construction (RCS-ADC). RCS-ADC simulates molecular ionisation from first principles, combining the accurate description of electron correlation of quantum chemistry with the full account of the continuum dynamics of the photoelectron. Complete theoretical characterisation of the atto-ionised many-electron state and photo-induced attosecond charge dynamics is achieved by calculating the reduced ionic density matrix (R-IDM) of the bipartite ion-photoelectron system, with full inclusion of the correlated shakeup states. Deviations from the sudden approximation picture of photoionisation, in the low-photon-energy limit, are presented. The effect of the multi-channel interaction between the parent-ion and the emitted photoelectron on the onset of the quantum electronic coherences is analysed. Moreover, I show how the Schmidt decomposition of the R-IDM unravels the many-electron dynamics triggered by the pump, allowing for the identification of the key channels involved. Finally, I calculate the X-ray attosecond transient absorption spectra of XUV-ionised pyrazine. The results unveil the mapping of the ATAS measurement onto the quantum electronic coherences, and related non-zero R-IDM matrix elements, produced upon ionisation by the XUV pump laser pulse.

Control of Attosecond Entanglement and Coherence

Calculations are presented of vibrational wave packet dynamics in H_{2}^{+} ions formed by ionization of neutral H_{2} by a pair of attosecond extreme ultraviolet laser pulses, using time-delayed dissociation of the cation by an ultraviolet probe pulse. The strength of experimentally observable two-level quantum beats as a function of the attosecond two-pulse delay can be related to ion+photoelectron entanglement resulting from the ionization process. This conclusion is supported by an evaluation of the purity of the reduced ion and photoelectron density matrices.

A systematic construction of Gaussian basis sets for the description of laser field ionization and high-harmonic generation

A precise understanding of mechanisms governing the dynamics of electrons in atoms and molecules subjected to intense laser fields has a key importance for the description of attosecond processes such as the high-harmonic generation and ionization. From the theoretical point of view, this is still a challenging task, as new approaches to solve the time-dependent Schrödinger equation with both good accuracy and efficiency are still emerging. Until recently, the purely numerical methods of real-time propagation of the wavefunction using finite grids have been frequently and successfully used to capture the electron dynamics in small one- or two-electron systems. However, as the main focus of attoscience shifts toward many-electron systems, such techniques are no longer effective and need to be replaced by more approximate but computationally efficient ones. In this paper, we explore the increasingly popular method of expanding the wavefunction of the examined system into a linear combination of atomic orbitals and present a novel systematic scheme for constructing an optimal Gaussian basis set suitable for the description of excited and continuum atomic or molecular states. We analyze the performance of the proposed basis sets by carrying out a series of time-dependent configuration interaction calculations for the hydrogen atom in fields of intensity varying from 5 × 10¹³ W/cm² to 5 × 10¹⁴ W/cm². We also compare the results with the data obtained using Gaussian basis sets proposed previously by other authors.

Multi-wave mixing in the high harmonic regime: monitoring electronic dynamics

It has been demonstrated that electronic coherences across many eV can be detected in pump-probe experiments involving high harmonic sources. An additional degree of control over the phase matching can be employed by investigating a more general class of multi-wave mixing. Non-collinear multi-wave mixing of high harmonics with energy (q₁ω₁ + q₂ω₂) can be selectively detected along the direction of (q₁k₁ + q₂k₂). Simulations based on a recently developed semi-perturbative approach show that only the specific harmonic signals with q₁ω₁ close to the energy difference between ground state and excited states are observable when the two input pulses are well separated in time. The coherent dynamics between different states can be selectively tracked by detecting the time-delay dependent signals with different q₁k₁, which can overcome the potential spectral congestion in real experiments. Additionally, such non-collinear geometry can be used to separate the dephasing induced decay and collision induced recovery behaviors of pump-probe high harmonic signal typically observed in the time-resolved high harmonic pump-probe signals.

Detecting electronic coherences by time-domain high-harmonic spectroscopy

Ultrafast spectroscopy is capable of monitoring electronic and vibrational states. For electronic states a few eV apart, an X-ray laser source is required. We propose an alternative method based on the time-domain high-order harmonic spectroscopy where a coherent superposition of the electronic states is first prepared by the strong optical laser pulse. The coherent dynamics can then be probed by the higher-order harmonics generated by the delayed probe pulse. The high nonlinearity typically modeled by the three-step mechanism introduced by Lewenstein and Corkum can serve as a recipe for generation of the coherent excitation with broad bandwidth. The main advantage of the method is that only optical (non-X-ray) lasers are needed. A semiperturbative model based on the Liouville space superoperator approach is developed for the bookkeeping of the different orders of the nonlinear response for the high-order harmonic generation using multiple pulses. Coherence between bound electronic states is monitored in the harmonic spectra from both first- and second-order responses.

Competition between charge migration and charge transfer induced by nuclear motion following core ionization: Model systems and application to Li2. J Chem Phys 2019;151:124108. [PMID: 31575213 DOI: 10.1063/1.5117246]

Attosecond and femtosecond spectroscopies present opportunities for the control of chemical reaction dynamics and products, as well as for quantum information processing; we address the somewhat unique situation of core-ionization spectroscopy which, for dimeric chromophores, leads to strong valence charge localization and hence tightly paired potential-energy surfaces of very similar shape. Application is made to the quantum dynamics of core-ionized Li₂ ⁺. This system is chosen as Li₂ is the simplest stable molecule facilitating both core ionization and valence ionization. First, the quantum dynamics of some model surfaces are considered, with the surprising result that subtle differences in shape between core-ionization paired surfaces can lead to dramatic differences in the interplay between electronic charge migration and charge transfer induced by nuclear motion. Then, equation-of-motion coupled-cluster calculations are applied to determine potential-energy surfaces for 8 core-excited state pairs, calculations believed to be the first of their type for other than the lowest-energy core-ionized molecular pair. While known results for the lowest-energy pair suggest that Li₂ ⁺ is unsuitable for studying charge migration, higher-energy pairs are predicted to yield results showing competition between charge migration and charge transfer. Central is a focus on the application of Hush's 1975 theory for core-ionized X-ray photoelectron spectroscopy to understand the shapes of the potential-energy surfaces and hence predict key features of charge migration.

Onset of ionic coherence and ultrafast charge dynamics in attosecond molecular ionisation

Here is presented a fully ab initio theoretical framework for simulating the correlated many-electron dynamics occurring during and emerging from molecular ionisation by attosecond laser pulses. This is based on the time-dependent (TD) version of the B-spline restricted correlation space (RCS)-algebraic diagrammatic construction (ADC) method, with the full description of the photoelectron and inclusion of electron correlation effects, such as shakeup processes and inter-channel couplings. The nature of the ultrafast charge dynamics in the molecular ion is elucidated by quantitatively predicting the degree of electronic coherence and eigenstate content of the prepared molecular cationic state, beyond the commonly used sudden approximation. The results presented here for the acetylene and ethylene molecules show that even in the high photon energy regime the simulated hole dynamics is quantitatively different from the prediction of the sudden approximation. Moreover, for high-bandwidth ionising pulse, the residual interaction between the cation, in highly-excited shake-up states, and the emitted slow photoelectron gives rise to a loss of coherence in the ionic system which can persist for the first few femtoseconds after ionisation.

Restricted Correlation Space B-Spline ADC Approach to Molecular Ionization: Theory and Applications to Total Photoionization Cross-Sections

Herein is presented a new approach to the ab initio algebraic diagrammatic construction (ADC) schemes for the polarization propagator, which is explicitly designed to accurately and efficiently describe molecular ionization. The restricted correlation space (RCS) version of the ADC methods up to second order of perturbation theory is derived via the intermediate state representation (ISR) and implemented in the multicenter B-spline basis set for the electronic continuum. Remarkably a general close-coupling structure of the RCS-ADC many-electron wave function, connecting the N-particle to the ( N - 1)-particle ADC intermediate states, emerges naturally as a nontrivial result of the RCS ansatz. Moreover, the introduced RCS-ADC schemes prove to be particularly manageable from a computational point of view, overcoming the practical limitations of the conventional ADC approaches. The quality of the new RCS-ADC( n) approaches is verified by performing a series of total photoionization cross-section calculations on a test set of molecules. The excellent agreement of the results with existing accurate benchmarks demonstrates that the RCS versions of the ADC schemes are optimal and quantitatively accurate methods for studying multichannel molecular photoionization.

CHEERS: a tool for correlated hole-electron evolution from real-time simulations

We put forward a practical nonequilibrium Green's function (NEGF) scheme to perform real-time evolutions of many-body interacting systems driven out of equilibrium by external fields. CHEERS is a computational tool to solve the NEGF equation of motion in the so called generalized Kadanoff-Baym ansatz and it can be used for model systems as well as first-principles Hamiltonians. Dynamical correlation (or memory) effects are added to the Hartree-Fock dynamics through a many-body self-energy. Applications to time-dependent quantum transport, time-resolved photoabsorption and other ultrafast phenomena are discussed.