Dissociation, ionization, and Coulomb explosion of H2+ in an intense laser field by numerical integration of the time-dependent Schrödinger equation

Efficient Control of Electron Localization and Probability Modulation with Synthesized Two-Color Intense Laser Pulses

A coupled electron-nuclear dynamical study at attosecond time scale is performed on the HD+ and H2+ molecular ions under the influence of synthesized intense two-color electric fields. We have employed ω - 2ω and also, ω - 3ω two-color fields in the infrared/mid-infrared regime to study the different fragmentation processes originating from the interference of n - (n + i) (i = 1, 2) photon absorption pathways. The branching ratios corresponding to different photofragments are controlled by tuning the relative phase as well as intensity of the two-color pulses, while the effect of the initial nuclear wave function is also studied by taking an individual vibrational eigenstate or a coherent superposition of several eigenstates of HD+ and H2+. By comprehensive analysis, the efficacy of the two different types of synthesized two-color pulses (ω - 2ω and ω - 3ω) are analyzed with respect to one-color intense pulses in terms of controlling the probability modulation and electron localization asymmetry and compared with previous theoretical calculations and experimental findings. Through the detailed investigation, we have addressed which one is the major controlling knob to have better electron localization as well as probability modulation.

Strong Laser Field-Driven Coupled Electron-Nuclear Dynamics: Quantum vs Classical Description

We have performed a coupled electron-nuclear dynamics study of H2+ molecular ions under the influence of an intense few-cycle 4.5 fs laser pulse with an intensity of 4 × 1014 W/cm2 and a central wavelength of 750 nm. Both quantum and classical dynamical methods are employed in the exact similar initial conditions with the aim of head-to-head comparison of two methodologies. A competition between ionization and dissociation channel is explained under the framework of quantum and classical dynamics. The origin of the electron localization phenomena is elucidated by observing the molecular and electronic wave packet evolution pattern. By probing with different carrier envelope phase (CEP) values of the ultrashort pulse, the possibility of electron localization on either of the two nuclei is investigated. The effects of initial vibrational states on final dissociation and ionization probabilities for several CEP values are studied in detail. Finally, asymmetries in the dissociation probabilities are calculated and mutually compared for both quantum and classical dynamical methodologies, whereas Franck-Condon averaging over the initial vibrational states is carried out in order to mimic the existing experimental conditions.

Probing Resonant Photoionization Time Delay by Self-Referenced Molecular Attoclock

Attosecond time-resolved electron tunneling dynamics have been investigated by using attosecond angular streaking spectroscopy, where a clock reference to the laser field vector is required in atomic strong-field ionization and the situation becomes complicated in molecules. Here we reveal a resonant ionization process via a transient state by developing an electron-tunneling-site-resolved molecular attoclock in Ar-Kr^{+}. Two distinct deflection angles are observed in the photoelectron angular distribution in the molecular frame, corresponding to the direct and resonant ionization pathways. We find the electron is temporally trapped in the Coulomb potential wells of the Ar-Kr^{+} before finally releasing into the continuum when the electron tunnels through the internal barrier. By utilizing the direct tunneling ionization as a self-referenced arm of the attoclock, the time delay of the electron trapped in the resonant state is revealed to be 3.50±0.04  fs. Our results give an impetus to exploring the ultrafast electron dynamics in complex systems and also endow a semiclassical presentation of the electron trapping dynamics in a quantum resonant state.

Exact Factorization Adventures: A Promising Approach for Non-Bound States

Modeling the dynamics of non-bound states in molecules requires an accurate description of how electronic motion affects nuclear motion and vice-versa. The exact factorization (XF) approach offers a unique perspective, in that it provides potentials that act on the nuclear subsystem or electronic subsystem, which contain the effects of the coupling to the other subsystem in an exact way. We briefly review the various applications of the XF idea in different realms, and how features of these potentials aid in the interpretation of two different laser-driven dissociation mechanisms. We present a detailed study of the different ways the coupling terms in recently-developed XF-based mixed quantum-classical approximations are evaluated, where either truly coupled trajectories, or auxiliary trajectories that mimic the coupling are used, and discuss their effect in both a surface-hopping framework as well as the rigorously-derived coupled-trajectory mixed quantum-classical approach.

Controlling the multi-electron dynamics in the high harmonic spectrum from N2O molecule using TDDFT

In this study, high harmonic generation from a multi-atomic nitrous oxide molecule was investigated. A comprehensive three-dimensional calculation of the molecular dynamics and electron trajectories through an accurate time-dependent density functional theory was conducted to efficiently explore a broad harmonic plateau. The effects of multi-electron and inner orbitals on the harmonic spectrum and generated coherent attosecond pulses were analyzed. The role of the valence electrons in controlling the process and extending the harmonic plateau was investigated. The main issue of producing a super-continuum harmonic spectrum via a frequency shift was considered. The time-frequency representation by means of a wavelet transform of the induced dipole acceleration provided a good insight into the distorted effects from the nonlinear processes in high harmonic emission. The effect of the chirped laser pulse on the production of broadband amplitude was justified in this model. By adjusting the optimal laser parameters to an input intensity of 2.5 × 10¹⁴ W cm^-2, an isolated 68 as pulse was generated.

A redshift mechanism of high-order harmonics: Change of ionization energy

We theoretically study the high-order harmonic generation of H₂⁺ and its isotopes beyond the Born-Oppenheimer dynamics. It is surprising that the spectral redshift can still be observed in high harmonic spectra of H₂⁺ driven by a sinusoidal laser pulse in which the trailing (leading) edge of the laser pulse is nonexistent. The results confirm that this spectral redshift originates from the reduction in ionization energy between recombination time and ionization time, which is obviously different from the nonadiabatic spectral redshift induced by the falling edge of the laser pulse. Additionally, the improved instantaneous frequency of harmonics by considering the changeable ionization energy can deeply verify our results. Therefore, this new mechanism must be taken into account when one uses the nonadiabatic spectral redshift to retrieve the nuclear motion.

Alignment dependent ultrafast electron-nuclear dynamics in molecular high-order harmonic generation

We investigated the high-order harmonic generation (HHG) process of diatomic molecular ion H₂⁺ in non-Born-Oppenheimer approximations (NBOA). The corresponding three-dimensional time-dependent Schrödinger equation is solved with arbitrary alignment angles. It is found that the nuclear motion can lead to spectral modulation of HHG in both the tunneling and multiphoton ionization regimes. The universal redshifts of the whole spectrum are unique in molecular HHG. The spectral width of HHG increases in NBOA. We calculated possible influences on redshifts of HHG in real experimental conditions and found that redshifts decrease with the increase of alignment angles of the molecules and are sensitive to the initial vibrational states. It can be used to extract the ultrafast electron-nuclear dynamics and image molecular structure. It will be instructive to related experiments.

Dissociation and Ionization of Quasi-Periodically Vibrating H2+ in Intense Few-Cycle Mid-Infrared Laser Fields

Using quantum mechanics calculations, we theoretically study the dissociation and ionization dynamics of the hydrogen-molecule ion in strong laser fields. Having prepared the nuclear wave packet of H₂⁺ in a specific vibrational state, a pump laser is used to produce a vibrational excitation, leading to quasi-periodical vibration without ionization. Then, a time-delayed few-cycle laser is applied to trigger the dissociation or ionization of H₂⁺. Both the time delay and the intensity of the probe laser alter the competition between dissociation and ionization. We also explore the dependence of kinetic-energy release spectra of fragments on the time delay, showing that the channels of above-threshold dissociation and below-threshold dissociation are opened and closed periodically. Also, dissociation from different channels is influenced by nuclear motion. The dissociation mechanism has been described in detail using the Floquet picture. This work provides a useful method for steering the electronic and nuclear dynamics of diatomic molecules in intense laser fields.

Toward the Full Quantum Dynamical Description of Photon Induced Processes in D2. J Phys Chem A 2016;120:9411-9421. [PMID: 27934332 DOI: 10.1021/acs.jpca.6b09623]

The dissociative ionization (multiphoton regime) of the D₂⁺ ion by ultrashort laser pulses has been studied theoretically using ab initio calculations. The combined ionization and dissociation spectrum was explored for fixed molecular axis orientations. In accordance with previous investigations, the dominant features in the obtained joint energy spectrum were multiphoton peaks. In addition to this, in the present work, photoelectron angular distributions were analyzed as well. By performing a partial wave analysis for each multiphoton peak, we have identified the number of absorbed photons. Moreover, we also found that the angular distribution can significantly change inside a multiphoton peak as a function of electron and nuclear kinetic energy.

Coupled Electron-Nuclear Dynamics on H2+ within Time-Dependent Born-Oppenheimer Approximation

Quantum dynamical behavior of H₂⁺ in the presence of a linearly polarized, ultrashort, intense, infrared laser pulse has been studied by numerically solving the time-dependent Schrödinger equation with nuclear motion restricted in one-dimension along the direction of laser polarization and electronic motion in three-dimensions. On the basis of the time-dependent Born-Oppenheimer approximation, we have constructed time-dependent potentials for the ground electronic state (1sσ_g) of H₂⁺. Subsequent nuclear dynamics is then carried out on these field-dressed potential energy surfaces, and the dissociation dynamics is investigated. Our analyses reveal that although the electronic longitudinal degree of freedom plays the major role in governing the dissociation dynamics, contributions from the electronic transverse degree of freedom should also have to be taken into account to obtain accurate results. Also, modeling electron-nuclei Coulomb interactions in a one-dimensional calculation with an artificially chosen constant softening parameter leads to a discrepancy with the exact results. Comparing our results with other quantum and classical dynamical studies showed a good agreement with exact results.

Role of electron heating in efficient interaction of a nanosecond laser with the cluster media: a case study on tetrahydrofuran cluster system

Interaction of tetrahydrofuran clusters with nanosecond laser pulses has been investigated at 532 and 1064 nm, using a time-of-flight mass spectrometer and home-built electron analyzer setup.

Direct Signature of Light-Induced Conical Intersections in Diatomics

Nonadiabatic effects are ubiquitous in physics, chemistry, and biology. They are strongly amplified by conical intersections (CIs), which are degeneracies between electronic states of triatomic or larger molecules. A few years ago it was revealed that CIs in molecular systems can be formed by laser light, even in diatomics. Because of the prevailing strong nonadiabatic couplings, the existence of such laser-induced conical intersections (LICIs) may considerably change the dynamical behavior of molecular systems. By analyzing the photodissociation process of the D2+ molecule carefully, we found a robust effect in the angular distribution of the photofragments that serves as a direct signature of the LICI, providing undoubted evidence of its existence.

Probing nuclear motion by frequency modulation of molecular high-order harmonic generation

Molecular high-order harmonic generation (MHOHG) in a non-Born-Oppenheimer treatment of H(2)(+), D(2)(+), is investigated by numerical simulations of the corresponding time-dependent Schrödinger equations in full dimensions. As opposed to previous studies on amplitude modulation of intracycle dynamics in MHOHG, we demonstrate redshifts as frequency modulation (FM) of intercycle dynamics in MHOHG. The FM is induced by nuclear motion using intense laser pulses. Compared to fixed-nuclei approximations, the intensity of MHOHG is much higher due to the dependence of enhanced ionization on the internuclear distance. The width and symmetry of the spectrum of each harmonic in MHOHG encode rich information on the dissociation process of molecules at the rising and falling parts of the laser pulses, which can be used to retrieve the nuclear dynamics. Isotope effects are studied to confirm the FM mechanism.

Quantum wave-packet dynamics in spin-coupled vibronic states

Extending the Shin-Metiu two-electron Hamiltonian, we construct a new Hamiltonian with effective singlet-triplet couplings. The Born-Oppenheimer electronic potentials and couplings are obtained for different parameters, and the laser-free dynamics is calculated with the full Hamiltonian and in the adiabatic limit. We compare the dynamics of the system using nuclear wave packets for different numbers of Born-Oppenheimer potentials and vibronic wave packets on a full 3-dimensional (two electron coordinates plus one nuclear coordinate) grid. Using strong fields, we show that it is possible to dynamically lock the spin state of the system by decoupling the singlet-triplet transition via a nonresonant dynamic Stark effect in the adiabatic limit. Although a similar spin-locking mechanism is observed in the dynamics of vibronic wave packets, multiphoton ionization cannot be neglected leading to the breakdown of the control scheme.

Multiphoton above threshold effects in strong-field fragmentation

We present a study of multiphoton dissociative ionization from molecules. By solving the time-dependent Schrödinger equation for H(2)(+) and projecting the solution onto double continuum scattering states, we observe the correlated electron-nuclear ionization dynamics in detail. We show-for the first time-how multiphoton structure prevails as long as one accounts for the energies of all the fragments. Our current work provides a new avenue to analyze strong-field fragmentation that leads to a deeper understanding of the correlated molecular dynamics.

THE APPLICATION OF BESSEL DISCRETE VARIABLE REPRESENTATION TO ATOMIC HYDROGEN IN INTENSE LASER FIELDS

The Bessel discrete variable representation (DVR) method is tested to describe the interaction of atomic hydrogen with intense laser fields by numerically solving the time-dependent Schrödinger equation. Using the Bessel functions of the first kind, the singular terms, r-2 or r-1 at the origin, in the kinetic energy operators are analytically solved. As an illustration example, the high-order harmonic generation (HOHG) spectra in atomic hydrogen is calculated in length and acceleration forms. From the numerical results, it is concluded that this simple Bessel DVR may be a useful method for describing the interaction of atomic hydrogen with intense laser fields.

Classical dynamics of laser-driven D₃⁺

A classical model of the triatomic D₃⁺ molecule subjected to an intense, few-cycle laser pulse is introduced. The model is capable of describing the laser-induced correlated motion of both electrons and nuclei in three dimensions, and allows us to follow the motion of the two electrons and three deuterons from the initial field-free state, during the pulse, and until the bond breaking into the final fragments. By averaging over many trajectories, we calculate the relative yields of the ionization and dissociation channels, as well as the kinetic energy release (KER) from the fragment ions. A comparison with recent experimental KER spectra shows good qualitative agreement. In addition, we find a pathway in which an emitted electron recombines into a high-lying Rydberg state, resulting in D + D⁺ + D⁺ fragments with the same KER as in the D⁺ + D⁺ + D⁺ channel.

Theoretical study of the photodissociation of Li(2)+ in one-color intense laser fields

A theoretical treatment of the photodissociation of the molecular ion Li(2) (+) in one-color intense laser fields, using the time-dependent wave packet approach in a Floquet Born-Oppenheimer representation, is presented. Six electronic states 1,2 (2)Σ(g)(+), 1,2 (2)Σ(u)(+), 1 (2)Π(g), and 1 (2)Π(u) are of relevance in this simulation and have been included. The dependences of the fragmental dissociation probabilities and kinetic energy release (KER) spectra on pulse width, peak intensity, polarization angle, wavelength, and initial vibrational level are analyzed to interpret the influence of control parameters of the external field. Three main dissociation channels, 1 (2)Σ(g)(+) (m = -1), 2 (2)Σ(g)(+) (m = -2), and 2 (2)Σ(u)(+) (m = -3), are seen to dominate the dissociation processes under a wide variety of laser conditions and give rise to well separated groups of KER features. Different dissociation mechanisms for the involved Floquet channels are discussed.

Exact factorization of the time-dependent electron-nuclear wave function

We present an exact decomposition of the complete wave function for a system of nuclei and electrons evolving in a time-dependent external potential. We derive formally exact equations for the nuclear and electronic wave functions that lead to rigorous definitions of a time-dependent potential energy surface (TDPES) and a time-dependent geometric phase. For the H(2)(+) molecular ion exposed to a laser field, the TDPES proves to be a useful interpretive tool to identify different mechanisms of dissociation.

Photodissociation of H2(+) upon exposure to an intense pulsed photonic Fock state

Producing and controlling nonclassical light states are now the subject of intense experimental efforts. In this paper we consider the interaction of such a light state with a small molecule. Specifically, we develop the theory and apply it numerically to calculate in detail how a short pulse of nonclassical light, such as the high intensity Fock state, induces photodissociation in H(2)(+). We compare the kinetic energy distributions and photodissociation yields with the analogous results of quasi-classical light, namely a coherent state. We find that Fock-state light decreases the overall probability of dissociation for low vibrational states of H(2)(+) as well as the location of peaks and line shapes in the kinetic energy distribution of the nuclei.

Attosecond-resolution quantum dynamics calculations for atoms and molecules in strong laser fields

A parallel quantum electron and nuclei wave packet computer code, LZH-DICP, has been developed to study laser-atom-molecule interaction in the nonperturbative regime with attosecond resolution. The nonlinear phenomena occurring in that regime can be studied with the code in a rigorous way by numerically solving the time-dependent Schrödinger equation of electrons and nuclei. Time propagation of the wave functions is performed using a split-operator approach, and based on a sine discrete variable representation. Photoelectron spectra for hydrogen and kinetic-energy spectra for molecular hydrogen ion in linearly polarized laser fields are calculated using a flux operator scheme, which testifies to the validity and the high efficiency of LZH-DICP.

Nonadiabatic chemical dynamics in an intense laser field: electronic wave packet coupled with classical nuclear motions

Dynamics of molecules in an intense laser field is studied in terms of the quantum electronic wave packet coupled with classical nuclear motions. The equations of motion are derived taking a proper account of molecular interactions with the vector potential of a classical electromagnetic field, along with the nonadiabatic interaction due to the breakdown of the Born-Oppenheimer approximation. With the aid of electronic structure calculations, the present method enables us to track, in an ab initio manner, the dynamics of polyatomic molecules in an intense field. Preliminary calculations are carried out for the vibrational state of LiF and a collision of Li+F under an intense laser pulse, which are limited to the domain of no ionization.

Electronic optical response of molecules in intense fields: Comparison of TD-HF, TD-CIS, and TD-CIS(D) approaches

Time-dependent Hartree-Fock (TD-HF) and time-dependent configuration interaction (TD-CI) methods with Gaussian basis sets have been compared in modeling the response of hydrogen molecule, butadiene, and hexatriene exposed to very short, intense laser pulses (760 nm, 3 cycles). After the electric field of the pulse returns to zero, the molecular dipole continues to oscillate due to the coherent superposition of excited states resulting from the nonadiabatic excitation caused by the pulse. The Fourier transform of this residual dipole gives a measure of the nonadiabatic excitation. For low fields, only the lowest excited states are populated, and TD-CI simulations using singly excited states with and without perturbative corrections for double excitations [TD-CIS(D) and TD-CIS, respectively] are generally in good agreement with the TD-HF simulations. At higher field strengths, higher states are populated and the methods begin to differ significantly if the coefficients of the excited states become larger than approximately 0.1. The response of individual excited states does not grow linearly with intensity because of excited state to excited state transitions. Beyond a threshold in the field strength, there is a rapid increase in the population of many higher excited states, possibly signaling an approach to ionization. However, without continuum functions, the present TD-HF and TD-CI calculations cannot model ionization directly. The TD-HF and TD-CIS simulations are in good accord because the excitation energies obtained by linear response TD-HF [also known as random phase approximation (RPA)] agree very well with those obtained from singly excited configuration interaction (CIS) calculations. Because CIS excitation energies with the perturbative doubles corrections [CIS(D)] are on average lower than the CIS excitation energies, the TD-CIS(D) response is generally stronger than TD-CIS.