An Orbital Basis Set for Double Photoionization of Atoms and Molecules.

The ab initio theoretical treatment of one-photon double photoionization processes has been limited to atoms and diatomic molecules by the challenges posed by large grid-based representations of the double ionized continuum wave function. To provide a path for extensions to polyatomics, an energy-adapted orbital basis approach is demonstrated that reduces the dimensions of such representations and simultaneously allows larger time steps in time-dependent computational descriptions of double ionization. Additionally, an algorithm that exploits the diagonal nature of the two-electron integrals in the grid basis and dramatically accelerates the transformation between grid and orbital representations is presented. Excellent agreement between the present results and benchmark theoretical calculations is found for H- and Be atoms, as well as the hydrogen molecule, including for the triply differential cross sections that relate the angular distribution and energy sharing of all of the particles in the molecular frame.

Attosecond delays in X-ray molecular ionization.

The photoelectric effect is not truly instantaneous but exhibits attosecond delays that can reveal complex molecular dynamics1-7. Sub-femtosecond-duration light pulses provide the requisite tools to resolve the dynamics of photoionization8-12. Accordingly, the past decade has produced a large volume of work on photoionization delays following single-photon absorption of an extreme ultraviolet photon. However, the measurement of time-resolved core-level photoionization remained out of reach. The required X-ray photon energies needed for core-level photoionization were not available with attosecond tabletop sources. Here we report measurements of the X-ray photoemission delay of core-level electrons, with unexpectedly large delays, ranging up to 700 as in NO near the oxygen K-shell threshold. These measurements exploit attosecond soft X-ray pulses from a free-electron laser to scan across the entire region near the K-shell threshold. Furthermore, we find that the delay spectrum is richly modulated, suggesting several contributions, including transient trapping of the photoelectron owing to shape resonances, collisions with the Auger-Meitner electron that is emitted in the rapid non-radiative relaxation of the molecule and multi-electron scattering effects. The results demonstrate how X-ray attosecond experiments, supported by comprehensive theoretical modelling, can unravel the complex correlated dynamics of core-level photoionization.

Bond-length dependence of attosecond ionization delays in O2 arising from electron correlation to a shape resonance.

We experimentally and theoretically demonstrate that electron correlation can cause the bond-length sensitivity of a shape resonance to induce an unexpected vibrational state-dependent ionization delay in a nonresonant channel. This discovery was enabled by a high-resolution attosecond-interferometry experiment based on a 400-nm driving and dressing wavelength. The short-wavelength driver results in a 6.2-electron volt separation between harmonics, markedly reducing the spectral overlap in the measured interferogram. We demonstrate the promise of this method on O2, a system characterized by broad vibrational progressions and a dense photoelectron spectrum. We measure a 40-attosecond variation of the photoionization delays over the X2Πg vibrational progression. Multichannel calculations show that this variation originates from a strong bond-length dependence of the energetic position of a shape resonance in the [Formula: see text] channel, which translates to the observed effects through electron correlation. The unprecedented energy resolution and delay accuracies demonstrate the promise of visible-light-driven molecular attosecond interferometry.

Filming enhanced ionization in an ultrafast triatomic slingshot.

Filming atomic motion within molecules is an active pursuit of molecular physics and quantum chemistry. A promising method is laser-induced Coulomb Explosion Imaging (CEI) where a laser pulse rapidly ionizes many electrons from a molecule, causing the remaining ions to undergo Coulomb repulsion. The ion momenta are used to reconstruct the molecular geometry which is tracked over time (i.e., filmed) by ionizing at an adjustable delay with respect to the start of interatomic motion. Results are distorted, however, by ultrafast motion during the ionizing pulse. We studied this effect in water and filmed the rapid "slingshot" motion that enhances ionization and distorts CEI results. Our investigation uncovered both the geometry and mechanism of the enhancement which may inform CEI experiments in many other polyatomic molecules.

Quantum-Path-Resolved Attosecond High-Harmonic Spectroscopy.

Strong-field ionization of molecules releases electrons which can be accelerated and driven back to recombine with their parent ion, emitting high-order harmonics. This ionization also initiates attosecond electronic and vibrational dynamics in the ion, evolving during the electron travel in the continuum. Revealing this subcycle dynamics from the emitted radiation usually requires advanced theoretical modeling. We show that this can be avoided by resolving the emission from two families of electronic quantum paths in the generation process. The corresponding electrons have the same kinetic energy, and thus the same structural sensitivity, but differ by the travel time between ionization and recombination-the pump-probe delay in this attosecond self-probing scheme. We measure the harmonic amplitude and phase in aligned CO_{2} and N_{2} molecules and observe a strong influence of laser-induced dynamics on two characteristic spectroscopic features: a shape resonance and multichannel interference. This quantum-path-resolved spectroscopy thus opens wide prospects for the investigation of ultrafast ionic dynamics, such as charge migration.

Strongly coupled intermediate electronic states in one-color two-photon single valence ionization of O2.

We present an experimental and theoretical energy- and angle-resolved investigation on the non-dissociative photoionization dynamics of near-resonant, one-color, two-photon, single valence ionization of neutral O2 molecules. Using 9.3 eV femtosecond pulses produced via high harmonic generation and a 3-D momentum imaging spectrometer, we detect the photoelectrons and O2 + cations produced from one-color, two-photon ionization in coincidence. The measured and calculated photoelectron angular distributions show agreement, which indicates that a superposition of two intermediate electronic states is dominantly involved and that wavepacket motion on those near-resonantly populated intermediate states does not play a significant role in the measured two-photon ionization dynamics. Here, we find greater utility in the diabatic representation compared to the adiabatic representation, where invoking a single valence-character diabat is sufficient to describe the underlying two-photon ionization mechanism.

Two-Center Interference in the Photoionization Delays of Kr_{2}.

We present the experimental observation of two-center interference in the ionization time delays of Kr_{2}. Using attosecond electron-ion-coincidence spectroscopy, we simultaneously measure the photoionization delays of krypton monomer and dimer. The relative time delay is found to oscillate as a function of the electron kinetic energy, an effect that is traced back to constructive and destructive interference of the photoelectron wave packets that are emitted or scattered from the two atomic centers. Our interpretation of the experimental results is supported by solving the time-independent Schrödinger equation of a 1D double-well potential, as well as coupled-channel multiconfigurational quantum-scattering calculations of Kr_{2}. This work opens the door to the study of a broad class of quantum-interference effects in photoionization delays and demonstrates the potential of attosecond coincidence spectroscopy for studying weakly bound systems.

Step-by-step state-selective tracking of fragmentation dynamics of water dications by momentum imaging.

The double photoionization of a molecule by one photon ejects two electrons and typically creates an unstable dication. Observing the subsequent fragmentation products in coincidence can reveal a surprisingly detailed picture of the dynamics. Determining the time evolution and quantum mechanical states involved leads to deeper understanding of molecular dynamics. Here in a combined experimental and theoretical study, we unambiguously separate the sequential breakup via D+ + OD+ intermediates, from other processes leading to the same D+ + D+ + O final products of double ionization of water by a single photon. Moreover, we experimentally identify, separate, and follow step by step, two pathways involving the b 1Σ+ and a 1Δ electronic states of the intermediate OD+ ion. Our classical trajectory calculations on the relevant potential energy surfaces reproduce well the measured data and, combined with the experiment, enable the determination of the internal energy and angular momentum distribution of the OD+ intermediate.

Attosecond spectroscopy of size-resolved water clusters.

Electron dynamics in water are of fundamental importance for a broad range of phenomena1-3, but their real-time study faces numerous conceptual and methodological challenges4-6. Here we introduce attosecond size-resolved cluster spectroscopy and build up a molecular-level understanding of the attosecond electron dynamics in water. We measure the effect that the addition of single water molecules has on the photoionization time delays7-9 of water clusters. We find a continuous increase of the delay for clusters containing up to four to five molecules and little change towards larger clusters. We show that these delays are proportional to the spatial extension of the created electron hole, which first increases with cluster size and then partially localizes through the onset of structural disorder that is characteristic of large clusters and bulk liquid water. These results indicate a previously unknown sensitivity of photoionization delays to electron-hole delocalization and indicate a direct link between electronic structure and attosecond photoionization dynamics. Our results offer new perspectives for studying electron-hole delocalization and its attosecond dynamics.

Coupled nuclear-electronic decay dynamics of O2 inner valence excited states revealed by attosecond XUV wave-mixing spectroscopy.

Multiple Rydberg series converging to the O2+c4Σ-u state, accessed by 20-25 eV extreme ultraviolet (XUV) light, serve as important model systems for the competition between nuclear dissociation and electronic autoionization. The dynamics of the lowest member of these series, the 3sσg state around 21 eV, has been challenging to study owing to its ultra-short lifetime (<10 fs). Here, we apply transient wave-mixing spectroscopy with an attosecond XUV pulse to investigate the decay dynamics of this electronic state. Lifetimes of 5.8 ± 0.5 fs and 4.5 ± 0.7 fs at 95% confidence intervals are obtained for v = 0 and v = 1 vibrational levels of the 3s Rydberg state, respectively. A theoretical treatment of predissociation and electronic autoionization finds that these lifetimes are dominated by electronic autoionization. The strong dependence of the electronic autoionization rate on the internuclear distance because of two ionic decay channels that cross the 3s Rydberg state results in the different lifetimes of the two vibrational levels. The calculated lifetimes are highly sensitive to the location of the 3s potential with respect to the decay channels; by slight adjustment of the location, values of 6.2 and 5.0 fs are obtained computationally for the v = 0 and v = 1 levels, respectively, in good agreement with experiment. Overall, an intriguing picture of the coupled nuclear-electronic dynamics is revealed by attosecond XUV wave-mixing spectroscopy, indicating that the decay dynamics are not a simple competition between isolated autoionization and predissociation processes.

Distinguishing resonance symmetries with energy-resolved photoion angular distributions from ion-pair formation in O2 following two-photon absorption of a 9.3 eV femtosecond pulse.

We present a combined experimental and theoretical study on the photodissociation dynamics of ion-pair formation in O2 following resonant two-photon absorption of a 9.3 eV femtosecond pulse, where the resulting O+ ions are detected using 3D momentum imaging. Ion-pair formation states of Σg-3 and 3Πg symmetry are accessed through predissociation of optically dark continuum Rydberg states converging to the B Σg-2 ionic state, which are resonantly populated via a mixture of both parallel-parallel and parallel-perpendicular two-photon transitions. This mixture is evident in the angular distribution of the dissociation relative to the light polarization and varies with the kinetic energy release (KER) of the fragmenting ion pair. The KER-dependent photoion angular distribution reveals the underlying two-photon absorption dynamics involved in the ion-pair production mechanism and indicates the existence of two nearly degenerate continuum resonances possessing different symmetries, which can decay by coupling to ion-pair states of the same total symmetry through internal conversion.

Probing molecular bond-length using molecular-frame photoelectron angular distributions.

The molecular-frame photoelectron angular distributions (MFPADs) in O 1s photoemission from CO2 molecule were measured. Patterns due to photoelectron diffractions were observed in the MFPADs. The polarization-averaged MFPADs were compared with theoretical calculation and were found to be useful in determining the molecular bond-length, which is a component to determine molecular structures.

Canonical Approach To Generate Multidimensional Potential Energy Surfaces.

A new force-based canonical approach for the accurate generation of multidimensional potential energy surfaces is demonstrated. Canonical transformations previously developed for diatomic molecules are used to construct accurate approximations to the 3-dimensional potential energy surface of the water molecule from judiciously chosen (adopting the right coordinate system) 1-dimensional planar slices that are shown to have the same canonical shape as the classical Lennard-Jones potential curve. Spline interpolation is then used to piece together the 1-dimensional canonical potential curves, to obtain the full 3-dimensional potential energy surface of a water molecule with a relative error less than 0.01. This work provides an approach to greatly reduce the computational cost of constructing potential energy surfaces in molecules from ab initio calculations. The canonical transformation techniques developed in this work illuminate a pathway to deepening our understanding of chemical bonding.

Is there any fundamental difference between ionic, covalent, and others types of bond? A canonical perspective on the question.

The concept of chemical bonding is normally presented and simplified through two models: the covalent bond and the ionic bond. Expansion of the ideal covalent and ionic models leads chemists to the concepts of electronegativity and polarizability, and thus to the classification of polar and non-polar bonds. In addition, the intermolecular interactions are normally viewed as physical phenomena without direct correlation to the chemical bond in any simplistic model. Contrary to these traditional concepts of chemical bonding, recently developed canonical approaches demonstrate a unified perspective on the nature of binding in pairwise interatomic interactions. This new canonical model, which is a force-based approach with a basis in fundamental molecular quantum mechanics, confirms much earlier assertions that in fact there are no fundamental distinctions among covalent bonds, ionic bonds, and intermolecular interactions including the hydrogen bond, the halogen bond, and van der Waals interactions.

Morse, Lennard-Jones, and Kratzer Potentials: A Canonical Perspective with Applications.

Canonical approaches are applied to classic Morse, Lennard-Jones, and Kratzer potentials. Using the canonical transformation generated for the Morse potential as a reference, inverse transformations allow the accurate generation of the Born-Oppenheimer potential for the H2+ ion, neutral covalently bound H2, van der Waals bound Ar2, and the hydrogen bonded one-dimensional dissociative coordinate in a water dimer. Similar transformations are also generated using the Lennard-Jones and Kratzer potentials as references. Following application of inverse transformations, vibrational eigenvalues generated from the Born-Oppenheimer potentials give significantly improved quantitative comparison with values determined from the original accurately known potentials. In addition, an algorithmic strategy based upon a canonical transformation to the dimensionless form applied to the force distribution associated with a potential is presented. The resulting canonical force distribution is employed to construct an algorithm for deriving accurate estimates for the dissociation energy, the maximum attractive force, and the internuclear separations corresponding to the maximum attractive force and the potential well.

Canonical Force Distributions in Pairwise Interatomic Interactions from the Perspective of the Hellmann-Feynman Theorem.

Force-based canonical approaches have recently given a unified but different viewpoint on the nature of bonding in pairwise interatomic interactions. Differing molecular categories (covalent, ionic, van der Waals, hydrogen, and halogen bonding) of representative interatomic interactions with binding energies ranging from 1.01 to 1072.03 kJ/mol have been modeled canonically giving a rigorous semiempirical verification to high accuracy. However, the fundamental physical basis expected to provide the inherent characteristics of these canonical transformations has not yet been elucidated. Subsequently, it was shown through direct numerical differentiation of these potentials that their associated force curves have canonical shapes. However, this approach to analyzing force results in inherent loss of accuracy coming from numerical differentiation of the potentials. We now show that this serious obstruction can be avoided by directly demonstrating the canonical nature of force distributions from the perspective of the Hellmann-Feynman theorem. This requires only differentiation of explicitly known Coulombic potentials, and we discuss how this approach to canonical forces can be used to further explain the nature of chemical bonding in pairwise interatomic interactions. All parameter values used in the canonical transformation are determined through explicit physical based algorithms, and it does not require direct consideration of electron correlation effects.

Canonical Approaches to Applications of the Virial Theorem.

Canonical approaches are applied for investigation of the extraordinarily accurate electronic ground state potentials of H2(+), H2, HeH(+), and LiH using the virial theorem. These approaches will be dependent on previous investigations involving the canonical nature of E(R), the Born-Oppenheimer potential, and F(R), the associated force of E(R), that have been demonstrated to be individually canonical to high accuracy in the case of the systems investigated. Now, the canonical nature of the remaining functions in the virial theorem [the electronic kinetic energy T(R), the electrostatic potential energy V(R), and the function W(R) = RF(R)] are investigated and applied to H2, HeH(+), and LiH with H2(+) chosen as reference. The results will be discussed in the context of a different perspective of molecular bonding that goes beyond previous direct applications of the virial theorem.

Canonical Potentials and Spectra within the Born-Oppenheimer Approximation.

A generalized formulation of canonical transformations and spectra are used to investigate the concept of a canonical potential strictly within the Born-Oppenheimer approximation. Data for the most accurate available ground electronic state pairwise intermolecular potentials in H2, HD, D2, HeH(+), and LiH are used to rigorously evaluate such transformations. The corresponding potentials are generated explicitly using parameters calculated with algebraic functions from that of the single canonical potential of the simplest molecule, H2(+). The efficacy of this approach is further tested by direct comparison of the predicted eigenvalues of all vibrational states in the selected molecular systems considered with the corresponding most accurately known Born-Oppenheimer eigenvalues currently available. Deviations are demonstrated to be less than 2 cm(-1) for all vibrational states in H2, HD, D2, HeH(+), and LiH, with an average standard deviation of 0.27 cm(-1) for the 87 states considered. The implications of these results for molecular quantum chemistry are discussed.

A general transformation to canonical form for potentials in pairwise interatomic interactions.

A generalized formulation of explicit force-based transformations is introduced to investigate the concept of a canonical potential in both fundamental chemical and intermolecular bonding. Different classes of representative ground electronic state pairwise interatomic interactions are referenced to a chosen canonical potential illustrating application of such transformations. Specifically, accurately determined potentials of the diatomic molecules H2, H2(+), HF, LiH, argon dimer, and one-dimensional dissociative coordinates in Ar-HBr, OC-HF, and OC-Cl2 are investigated throughout their bound potentials. Advantages of the current formulation for accurately evaluating equilibrium dissociation energies and a fundamentally different unified perspective on nature of intermolecular interactions will be emphasized. In particular, this canonical approach has significance to previous assertions that there is no very fundamental distinction between van der Waals bonding and covalent bonding or for that matter hydrogen and halogen bonds.

Vibrationally specific photoionization cross sections of acrolein leading to the X̃²A' ionic state.

The vibrational branching ratios in the photoionization of acrolein for ionization leading to the X̃²A' ion state were studied. Computed logarithmic derivatives of the cross section and the corresponding experimental data derived from measured vibrational branching ratios for several normal modes (ν9, ν10, ν11, and ν12) were found to be in relatively good agreement, particularly for the lower half of the 11-100 eV photon energy range considered. Two shape resonances have been found near photon energies of 15.5 and 23 eV in the photoionization cross section and have been demonstrated to originate from the partial cross section of the A' scattering symmetry. The wave functions computed at the resonance complex energies are delocalized over the whole molecule. By looking at the dependence of the cross section on the different normal mode displacements together with the wave function at the resonant energy, a qualitative explanation is given for the change of the cross sections with respect to changing geometry.