Molecules in intense laser fields: Enhanced ionization in a one-dimensional model of H2

Dissociative multi-ionization of N2O molecules in strong femtosecond laser field

Multi-ionization and subsequent Coulomb explosion (CE) of the N2O molecule irradiated by linearly polarized 800 nm laser field is investigated by a reaction microscope, where a number of CE channels of N2Oq+ with q{less than or equal to}5 for two-body fragmentation and q{less than or equal to}8 for three-body fragmentation were observed. For two-body CE, by analyzing the internuclear separations extracted from kinetic energy releases (KERs), dissociation branching fractions, and laser intensity dependence, interestingly we found that fragmentation N2O5+→N3++NO2+ is produced directly from dissociating N2O3+ via non-sequential stairstep ionization whereas most of others result from the sequential stairstep ionization. For three-body CE, 25 fragmentation channels of N2Oq+ (q = 3-8) are distinguished in present charge-encoded multi-photoion coincidence plot and the concerted fragmentation mechanism is nominated in a typical Dalitz plot. With the help of the numerical computation with the measured KERs and momentum correlation angles, the geometric structures of molecular ions prior to fragmentation are reconstructed, which display the bending motion and simultaneous two-bond stretching before the CE. Increasing of bond length for high charged N2Oq+ indicates the dominating stairstep ionization in three-body fragmentation.

Electronic non-adiabatic dynamics in enhanced ionization of isotopologues of hydrogen molecular ions from the exact factorization perspective

An exact-factorization perspective of enhanced ionization in isotopologues of H₂⁺ demonstrates the concept of the exact potential driving the electrons in non-adiabatic motion of molecules in strong fields, and sets a new platform for introducing various approximations.

Time-dependent quantum simulation of coronene photoemission spectra

Photoelectron spectroscopy is usually described by a simple equation that relates the binding energy of the photoemitted electron, Ebinding, its kinetic energy, Ekinetic, the energy of the ionizing photon, Ephoton, and the work function of the spectrometer, ϕ, Ebinding = Ephoton - Ekinetic - ϕ. Behind this equation there is an extremely rich physics, which we describe here using as an example a relatively simple conjugated molecule, namely coronene. The theoretical analysis of valence band and C1s core level photoemission spectra showed that multiple excitations play an important role in determining the intensities of the final spectrum. An explicit, time-evolving model is applied, which is able to count all possible photo-excitations occurring during the photoemission process, showing that they evolve on a short time-scale, of about 10 fs. The method reveals itself to be a valid approach to reproduce photoemission spectra of polycyclic aromatic hydrocarbons (PAHs).

Ab Initio/RRKM Study of Dissociation Mechanism of Benzene Trication

Density functional B3LYP /6-31 G (d,p) calculations have been performed in order to investigate isomerization and dissociation of [Formula: see text] in the ground electronic state, which are relevant to the Coulomb explosion mechanism of benzene. The results demonstrate that the benzene-like isomer of [Formula: see text], 1, can decompose through various pathways leading to distinct fragmentation products. The most kinetically favorable channel involves ring opening in 1 accompanied with 1,2-shifts of two hydrogen atoms followed by rotation around the middle C–C bond and dissociation to H 2 CCCH 2+ + H 2 CCCH + with exothermicity and the highest barrier of 100.2 and 38.4 kcal/mol, respectively. Several other product channels, including [Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text], share the same highest barrier for the initial ring opening step with the pathway producing H 2 CCCH 2+ + H 2 CCCH + but exhibit higher second highest barriers. Elimination of [Formula: see text] to form [Formula: see text] is 137.9 kcal/mol exothermic but the highest barrier on the pathway leading to these products is 53.0 kcal/mol. A simple proton elimination to produce [Formula: see text] is computed to be 82–84 kcal/mol exothermic and to depict a 64–65 kcal/mol barrier. RRKM calculations of rate constants for individual reaction steps assuming that the initial internal energy of 1 is 110 kcal/mol and solving kinetic master equations to obtain relative branching ratios show that H 2 CCCH 2+ + H 2 CCCH + are the dominant products (81.5%) followed by [Formula: see text] (13.2%) and the other minor products include [Formula: see text] (2.6%), [Formula: see text] (1.1%), [Formula: see text] (0.55%), [Formula: see text] (0.49%), [Formula: see text] (0.20%), and [Formula: see text] (0.14%). The fragments are expected to be produced with high translational energy due to high Coulomb repulsion energy barriers.

Internuclear-distance dependence of electron correlation in nonsequential double ionization of H(2)

Using three-dimensional classical ensembles, we have investigated the internuclear distance dependence of nonsequential double ionization (NSDI) of H(2) molecules by an 800 nm, 1x10(14) W/cm(2) laser pulse. For the internuclear distance R ranging from 2 to 12 a.u., the NSDI of H(2) provides rich correlation patterns in the two-electron momentum distributions. These correlation patterns essentially reveal different microscopic dynamics in NSDI process. Moreover, our calculations show that R approximately 4 a.u. is the critical distance for double ionization yield of H(2). These results are qualitatively explained based on the classical barrier expression model and back analysis.

Space-time contours to treat intense field-dressed molecular states

In this article we consider a molecular system exposed to an intense short-pulsed external field. It is a continuation of a previous publication [A. K. Paul, S. Adhikari, D. Mukhopadhyay et al., J. Phys. Chem. A 113, 7331 (2009)] in which a theory is presented that treats quantum effects due to nonclassical photon states (known also as Fock states). Since these states became recently a subject of intense experimental efforts we thought that they can be treated properly within the existing quantum formulation of dynamical processes. This was achieved by incorporating them in the Born-Oppenheimer (BO) treatment with time-dependent coefficients. The extension of the BO treatment to include the Fock states results in a formidable enhancement in numerical efforts expressed, in particular, in a significant increase in CPU time. In the present article we discuss an approach that yields an efficient and reliable approximation with only negligible losses in accuracy. The approximation is tested in detail for the dissociation process of H(2) (+) as caused by a laser field.

Theoretical study of unimolecular decomposition of allene cations

Ab initio coupled clusters and multireference perturbation theory calculations with geometry optimization at the density functional or complete active space self-consistent-field levels have been carried out to compute ionization energies and to unravel the dissociation mechanism of allene and propyne cations, C(3)H(4)(n+) (n=1-3). The results indicate that the dominant decomposition channel of the monocation is c-C(3)H(3)(+) + H, endothermic by 37.9 kcal/mol and occurring via a barrier of 43.1 kcal/mol, with possible minor contributions from H(2)CCCH(+) + H and HCCCH(+) + H(2). For the dication, the competing reaction channels are predicted to be c-C(3)H(3)(+) + H(+), H(2)CCCH(+) + H(+), and CCCH(+) + H(3)(+), with dissociation energies of -20.5, 8.5, and 3.0 kcal/mol, respectively. The calculations reveal a H(2)-roaming mechanism for the H(3)(+) loss, where a neutral H(2) fragment is formed first, then roams around and abstracts a proton from the remaining molecular fragment before leaving the dication. According to Rice-Ramsperger-Kassel-Marcus calculations of energy-dependent rate constants for individual reaction steps, relative product yields vary with the available internal energy, with c-C(3)H(3)(+) + H(+) being the major product just above the dissociation threshold of 69.6 kcal/mol, in the energy range of 70-75 kcal/mol, and CCCH(+) + H(3)(+) taking over at higher energies. The C(3)H(4)(3+) trication is found to be not very stable, with dissociation thresholds of 18.5 and 3.7 kcal/mol for allene and propyne, respectively. Various products of Coulomb explosion of C(3)H(4)(3+), H(2)CCCH(2+) + H(+), CHCHCH(2+) + H(+), C(2)H(2)(2+) + CH(2)(+), and CCH(2)(2+) + CH(2)(+) are highly exothermic (by 98-185 kcal/mol). The tetracation of C(3)H(4) is concluded to be unstable and therefore no more than three electrons can be removed from this molecule before it falls apart. The theoretical results are compared to experimental observations of Coulomb explosions of allene and propyne.

Asymmetric molecular gating for supercontinuous high harmonic generation in the plateau

Molecular wave packets provide an alternative route to steer the high harmonic generation (HHG). Here we present a new method for isolated attosecond pulse generation from high harmonics in the plateau with asymmetric molecules. It is shown that the photoionization depends on the molecular structure. Through steering the ionization, HHG is controlled with the asymmetric molecule and supercontinuous high harmonics are produced in the plateau, from which an isolated 95-attosecond pulse is generated. In contrast to the cutoff which declines sharply, the plateau shows almost a constant intensity, indicating a higher yield. Moreover, the plateau covers a very broad bandwidth, thus is preferable to produce an isolated attosecond pulse with a rather short duration.

Time-dependent density-functional theory with a self-interaction correction

We discuss an implementation of the self-interaction correction for the local-density approximation to time-dependent density-functional theory. A variational formulation is given, taking care of the necessary constraints. A manageable and transparent propagation scheme using two sets of wave functions is proposed and applied to laser excitation with subsequent ionization of a dimer molecule.

Strong-field induced vibrational coherence in the ground electronic state of hot I2

We observe large amplitude coherent vibrations in vibrationally hot ground state neutral I2 molecules created through "Lochfrass," or "R-selective ionization." We directly measure the phase and amplitude of the vibrations. Our results support the notion of enhanced ionization at large internuclear separation over recent theoretical predictions for heavy molecules. Furthermore, simulations of the vibrational motion show that for Lochfrass the vibrational coherence, contrary to most coherent control schemes, is stronger for hot molecules than for cold molecules.

Numerical simulation of nonadiabatic electron excitation in the strong-field regime. 3. Polyacene neutrals and cations

The electron optical response for a series of linear polyacenes and their molecular ions (mono and dications) in strong laser fields was studied using time-dependent Hartree-Fock theory. The interactions of benzene, naphthalene, anthracene, and tetracene with pulsed fields at a frequency of 1.55 eV and intensities of 8.77 x 10(13), 3.07 x 10(13), 1.23 x 10(13), and 2.75 x 10(12) W/cm2, respectively, were calculated using the 6-31G(d,p) basis set. Nonadiabatic processes, including nonadiabatic time evolution of the dipole moment, Löwden charges, and occupation numbers, were studied. The nonadiabatic response increased with the length of the molecule and was greatest for the molecular monocations. The only exception was tetracene, in which the very strong response of the dication was due to a near resonance with the applied field. The intensity and frequency dependence of the dipole moment response for the monocations of naphthalene and anthracene was also calculated. As the intensity increased, the population of higher-energy excited-states increased, and as the frequency increased, the excitation volume increased in good agreement with the Dykhne approximation.

Space-time contours to treat intense field-dressed molecular states. I. Theory

A molecular system exposed to an intense external field is considered. The strength of the field is measured by the number L of electronic states that become populated during this process. In the present article the authors discuss a rigorous way, based on the recently introduced space-time contours [R. Baer, et al., J. Chem. Phys. 119, 6998 (2003)], to form N coupled Schrodinger equations where N<L, which maintains the effects due to the remaining (L-N) populated states. It is shown that whereas the size of L is unlimited, the main requirement concerning N is that the original group of N field-free states forms a Hilbert subspace in the spatial region of interest. From previous studies it is known that a group of states forms a Hilbert subspace if and only if the corresponding topological D matrix is diagonal [M. Baer, et al., Farad, Discuss 127, 337 (2004)].

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.

Electronic and nuclear responses of fixed-in-space H2S to ultrashort intense laser fields

The Coulomb explosion dynamics of H2S, H2S3+-->H+ +S+ + H+, in ultrashort intense laser fields (12 fs, approximately 2 x 10(14) W/cm2) is studied by the coincidence momentum imaging of the three fragment ions. Different electronic and nuclear responses are identified depending on the direction of laser polarization epsilon in the molecular frame. The dependence can be interpreted in terms of the electronic and bonding characters of charge transfer states of H2S coupled to the electronic ground state.

Numerical Simulation of Nonadiabatic Electron Excitation in the Strong Field Regime. 2. Linear Polyene Cations

Time-dependent Hartree-Fock theory has been used to study the electronic optical response of a series of linear polyene cations (+1 and +2) in strong laser fields. The interaction of ethylene, butadiene, and hexatriene, with pulsed and CW fields corresponding to 8.75 x 10(13) W/cm(2) and 760 nm, have been calculated using the 6-31G(d,p) basis set. Nonadiabatic processes including nonlinear response of the dipole moment to the field and non-resonant energy deposition into excited states were more pronounced for the monocations in comparison with dications. For a given charge state and geometry, the nonadiabatic effects in the charge distribution and instantaneous dipole increased with the length of the polyene. For pulsed fields, the instantaneous dipole continued to oscillate after the field returned to zero and corresponded to a non-resonant electronic excitation involving primarily the lowest electronic transition. For a given molecule and fixed charge state, the degree of nonadiabatic coupling and excitation was greater for geometries with lower excitation energies.