Strong-field ionization and the resulting electronic dynamics are important for a range of processes such as high harmonic generation, photodamage, charge resonance enhanced ionization, and ionization-triggered charge migration. Modeling ionization dynamics in molecular systems from first-principles can be challenging due to the large spatial extent of the wavefunction which stresses the accuracy of basis sets, and the intense fields which require non-perturbative time-dependent electronic structure methods. In this paper, we develop a time-dependent density functional theory approach which uses a Gaussian-type orbital (GTO) basis set to capture strong-field ionization rates and dynamics in atoms and small molecules. This involves propagating the electronic density matrix in time with a time-dependent laser potential and a spatial non-Hermitian complex absorbing potential which is projected onto an atom-centered basis set to remove ionized charge from the simulation. For the density functional theory (DFT) functional we use a tuned range-separated functional LC-PBE*, which has the correct asymptotic 1/r form of the potential and a reduced delocalization error compared to traditional DFT functionals. Ionization rates are computed for hydrogen, molecular nitrogen, and iodoacetylene under various field frequencies, intensities, and polarizations (angle-dependent ionization), and the results are shown to quantitatively agree with time-dependent Schrödinger equation and strong-field approximation calculations. This tuned DFT with GTO method opens the door to predictive all-electron time-dependent density functional theory simulations of ionization and ionization-triggered dynamics in molecular systems using tuned range-separated hybrid functionals.

We have studied, experimentally and theoretically, the ionization probability of carbonyl sulfide (OCS) molecules in intense linearly-polarized 800 nm laser pulses as a function of the angle between the molecular axis and the laser polarization. Experimentally, the molecules are exposed to two laser pulses with a relative time delay. The first, weaker pulse induces a nuclear rotational wave packet within each molecule such that the ensemble exhibits preferential alignment in the laboratory frame at specific times. The second, stronger pulse induces ionization, and the variation in single and double ionization yields is measured as a function of the delay between the two pulses. The angular-dependence of the ionization yield is extracted by fitting the delay-dependent yields to a sum of delay-dependent moments of the rotational wave packet's angular distribution. We compute these same angular-dependent strong-field ionization yields for OCS using time-dependent density functional theory (TDDFT). For the single ionization case, our measurements agree well with TDDFT calculations and with previous experiments. Furthermore, analysis of the simulated one-body density reveals that, when averaged over a laser cycle, the resulting hole is delocalized across the molecule for light polarized perpendicular to the molecular axis, and mostly localized on the sulfur for parallel polarization. This suggests that preferential molecular alignment is a key parameter for controlling charge migration dynamics initiated by strong-field ionization. For double ionization, the agreement between experiment and theory is less compelling, reflecting the substantial challenges of computing double ionization yields using TDDFT methods.

We study, experimentally and theoretically, the ionization probability of singly halogenated methane molecules, CH3Cl and CH3Br, in intense linearly polarized 800 nm laser pulses as a function of the angle between the molecular axis and the laser polarization. Experimentally, the molecules are exposed to two laser pulses with a relative time delay. The first, weaker pulse induces a nuclear rotational wave packet within the molecules, which are then ionized by the second, stronger pulse. The angle-dependent ionization yields are extracted from fits of the measured delay-dependent ionization signal to a superposition of moments of the rotational wave packet’s angular distribution. Angle-dependent strong-field ionization (SFI) yields are also calculated using time-dependent density functional theory. Good agreement between measurements and theory is obtained. Interestingly, we find a marked difference between the angle-dependence of the ionization yields for these two halomethane species despite the similar structure of their highest occupied molecular orbitals. Calculations reveal that these differences are a result of multichannel (CH3Cl) vs single-channel (CH3Br) ionization and of increased hole localization on Br vs Cl. By adding calculations for CH3F, we can discern clear trends in the ionization dynamics with increasing halogen mass. These results are illustrative, as chemical functionalization and molecular alignment are likely to be important parameters for initiating and controlling charge migration dynamics via SFI.

Real-time methods are convenient for simulating core-level absorption spectra but suffer from nonphysical intruder peaks when using atom-centered basis sets. In transient absorption spectra, these peaks exhibit highly nonphysical time-dependent modulations in their energies and oscillator strengths. In this paper, we address the origins of these intruder peaks and propose a straightforward and effective solution based on a filtered dipole operator. In combination with real-time time-dependent density functional theory (RT-TDDFT), we demonstrate how to compute intruder-free attosecond transient X-ray absorption spectra for the aminophenol (C6H7NO) oxygen and nitrogen K-edges and the α-quartz (SiO2) silicon L-edge. Without filtering, the computed spectra are qualitatively wrong. This procedure is suitable for both static and transient inner-shell spectroscopy studies and can easily be implemented in a range of real-time methodologies.

We present a theoretical investigation of charge migration following strong-field ionization in a multi-electron system. We study a model homo-nuclear molecule with two electrons, each restricted to one dimension (1+1D), interacting with a strong, static electric field. We show that in this system charge migration results from the interplay between multiple ionization channels that overlap in space, creating a coherent electron-hole wave packet in the cation. We also find that, in our case, charge migration following the first ionization manifests as a modulation of the subsequent double ionization signal. We derive a parametrized semiclassical model from the full multi-electron system and we discuss the importance of the choice of cation electronic-structure basis for the efficacy of the semiclassical representation. We use the ab-initio solution of the full 1+1D system as a reference for the qualitative and quantitative results of the parametrized semiclassical model. We discuss the extension of our model to long wavelength time-dependent fields with full dimension, many-electron targets.