Introduction

Atomic ionization processes in the presence of strong laser fields are well modeled by single active electron ionization approximations at moderate intensities1,2. In contrast, non-sequential double ionization (NSDI) is complicated by strong electron correlation and Coulombic attraction to the ion core3,4,5. Processes analogous to those occurring in atoms take place in molecules, albeit with added complexity introduced by the presence of other atoms, molecular structure, additional degrees of freedom, and higher density of electronic states. Furthermore, unlike atoms, molecules, especially nonlinear polyatomic molecules, can experience electron recollision processes at any of their atoms, making possible the observation of NSDI with elliptical and even circularly polarized fields6,7,8,9,10,11,12,13,14.

Understanding and controlling the behavior of molecules in strong fields requires the ability to discriminate between different double ionization mechanisms. The choice of a molecule with high ionization potential and long-wavelength short-pulse excitation has been shown to favor recollision-induced molecular NSDI3. Shorter wavelengths and longer pulse durations may cause NSDI in addition to sequential multiphoton ionization (MPI), wherein large fragment ions absorb additional photons from the field and undergo further fragmentation and ionization15,16.

Controlling the ionization process would be useful to simplify the analysis of molecular ionization in strong fields. Moreover, controlling electron recollision could impact a wide range of phenomena that depend on this process, for example: molecular fragmentation, high-harmonic generation (HHG), the generation of attosecond pulses via HHG, above threshold ionization (ATI), and even attosecond clocking5,17. Control of these processes via pulse shaping has been of interest for over a decade18. Most relevant to the findings here are calculations that have shown that large effects can be gained by ‘jumps’ in the optical cycles of a pulse. For example, a pulse resulting from joining two identical pulses with one having its carrier envelope phase shifted by π, was predicted theoretically to extend the HHG cutoff19. Similarly, a theoretical exploration of pulses with an instantaneous π phase jump in the time domain predicted a significant extension of the HHG cutoff energy and the energies achieve by ATI20,21. Unfortunately, these schemes have not been experimentally implemented. The latter case, for example, would require pulses spanning more than five octaves, a capability that is well beyond the present state of the art in ultrafast science.

The purpose of the present study is to enhance or suppress the yield of metastable dications via spectral-phase pulse shaping. Inspired by the elegant control experiments by Silberberg using phase steps to control two-photon excitation22, we use a phase step. However, in this work, laser-matter interactions are well outside the perturbation limit considered by Silberberg. In fact, perturbation theory only predicts a significant reduction in the yield of high-order (five photons or more) processes. The phase step in the frequency domain causes a jump in the time-dependent frequency of the pulse during the time when tunnel ionization and recollision occurs. The jump in frequency is relevant because recollision, NSDI, and HHG have been found to depend on the frequency of the laser elevated to a power of approximately five or eight23,24,25, stemming from quantum path interferences with contribution from multiple returning orbits. In addition to calculations of the kinetic energy acquired by a free electron in the shaped laser field, we develop a model exclusively based on the fact that the frequency of the shaped pulses varies through the temporal pulse. Thus, we assume the position of the phase step with respect to the center frequency of the pulse and the sign of the phase step in the frequency domain is proportional to the frequency of the pulse at the times when recollision takes place.

Here, using laser pulses with identical spectrum and peak intensity, but different phase characteristics, we observe control over the yield of ions produced via NSDI and MPI mechanisms. We find that the model predicts quite closely the observed enhancement or suppression of dications as well as other electron recollision processes such as high-harmonic generation. We test the control mechanism using circularly polarized pulses, which involve longer orbits for the recolliding electron wavepacket and lower recollision energies. We find that the contrast observed for doubly charged fragment ions is slightly larger than that measured for linearly polarized pulses. This indicates that the control mechanism is associated either with the recolliding electron wavepacket orbit period, which is here manipulated by the position and sign of the phase step, or the kinetic energy of the recolliding electron, which is above or below the double ionization threshold due to pulse shaping. In summary, we report the observation of significant control over the yield of metastable dications as we scan a ¾π step across the spectrum of an intense femtosecond laser pulse. The enhancement or suppression of electron wavepacket recollision appears to be universal based on similar phase-step measurements, presented here, of the yield of helium double ionization and the intensity of gas-phase third harmonic generation.

Results

Conceptual background

Doubly charged ions, especially if metastable, are more likely to be formed by NSDI than by stepwise MPI. The goal of our experiments is to find if temporal shaping of the laser pulses can influence NSDI. The link between temporal shaping and electron recollision, which causes NSDI, is given by the kinetic energy of the electron which is accelerated by the electromagnetic field. Using shaped pulses with phase jumps in the time domain we thus seek to influence the orbit of the recolliding electron wavepacket and hence the probability for NSDI. Figure 1 illustrates how two different shaped pulses result in different electron recollision orbits, and thus different electron kinetic energies. Our goal is then to find shaped pulses for which a small change, such as a sign in the phase, causes the kinetic energy of the recolliding electron to be above or below threshold for NSDI.

Fig. 1: Schematic representation of ethane molecules exposed to two different shaped laser pulses.
figure 1

The intense laser field (red and blue) drives an electron out of a sigma bond from the molecule through tunnel ionization accelerating away and then towards the molecule when the sign of the field changes. If the kinetic energy of the electron is sufficient to release the second electron it leads to the creation of metastable dications or, alternatively, the molecule undergoes multiphoton ionization during the subsequent optical cycles. We demonstrate that pulse shaping with a phase jump modifies the electric field and hence modifies the available electron kinetic energy. The position and sign of the phase jump thus control the nonsequential double ionization process and the resulting molecular fragments observed.

Preliminary experiments

We irradiate isolated gas phase ethane molecules and collect the resulting ions using a time-of-flight (TOF) mass spectrometer. Figure 2 compares the region of the mass spectrum containing the fragment ions with mass to charge ratio (m/Q) between 13 and 15.5 obtained for ethane with 42-fs transform-limited (TL) 800-nm pulses that are either (a) linearly polarized, (b) circularly polarized, or (c) chirped to seven times their TL duration. The corresponding complete mass spectra are provided in Supplementary Figs. 13, and Table 1 has the appearance energy for most of the fragment ions observed in our measurement8,26,27,28,29,30. We note that the metastable dications fragments C2H32+, C2H42+, C2H52+, and C2H62+ are observed for both linear and circular polarization. In contrast, for chirped pulses, the yield of doubly charged ions is negligible. Note that for some species, such as C2H42+, the metastable doubly charged fragment ions may dissociate with a lifetime of hundreds of nanoseconds on their way to the detector30. In place of the relatively narrow peaks observed in Fig. 2a, b, one observes pairs of peaks in Fig. 2c. These split time of flight (TOF) peaks are associated with fragment-ion pairs due the Coulomb repulsion between the ion pairs, such as CH2+ + CH2+ and CH3+ + CH3+. Specifically, the peak pairs appear in the TOF spectra because some fragments are ejected toward the detector and thus arrive earlier, while others are ejected in the opposite direction and thus arrive later in time. The peak separation is proportional to the kinetic energy release resulting from the repulsion between positively charged ions. The appearance of Coulomb exploding ion pairs has been found to increase with pulse duration, in agreement with Fig. 2c. This supports the observation that longer pulses lead to the breakage of multiple bonds31,32,33.

Fig. 2
figure 2

Portion of the mass spectra of ethane obtained when using transform-limited (TL) 42 fs (a) linearly, (b) circularly polarized laser pulses, and (c) chirped pulses. At m/Q = 14 and 15 there is an overlap between singly charged CH2+ and CH3+ with doubly charged C2H42+ and C2H62+, respectively.

Table 1 Principal fragments and ions resulting from the strong field ionization of ethane.

The peaks of the metastable dications C2H32+ and C2H52+ are easily identified by their unique half-mass unit position at m/Q of 13.5 and 14.5, respectively. The loss of three or one neutral hydrogen atoms, respectively, occurs with minimal kinetic energy such that these peaks are not broadened substantially. In contrast, both doubly and singly charged fragment ions can contribute to the TOF peaks at m/Q = 14 and 15. It is likely that the peak at m/Q = 14 has a significant contribution of the metastable C2H42+ dication, based on the fact that for linearly and circularly polarized light, the peak at m/Q = 14 is ten times larger than the peak at 14.5. Flanked by the peak at m/Q = 14 are the singly charged, CH2+ + CH2+ Coulomb explosion fragment-ion pairs. For chirped pulses, no doubly charged ions are detected and the single peak at m/Q = 14 becomes a double peak indicating Coulomb explosion among singly charged ion pairs. The chirped pulses are long enough to allow chemical bonds to lengthen and give rise to enhanced ionization31. A clear assignment cannot be made for m/Q = 15, therefore we shall focus the rest of our discussion on C2H42+ and C2H52+. The metastable fragment-ion C2H32+ has a lower yield, making it difficult to interpret quantitatively.

The data shown in Fig. 2a, b suggest that the double ionization mechanism observed for linear and for circular polarization is very similar, but different for chirped pulses. The Keldysh parameter γ,34 helps identify the ionization regime by comparing the ionization potential of the atom or molecule to the cycle-averaged quiver energy of a free electron. When γ < 1, tunnel ionization can be considered to take place within an optical cycle of the pulse, whereas the regime when γ > 1 is known as the MPI regime. Based on the laser parameters and the ionization potential of ethane, NSDI is expected for 42-fs TL pulses linearly polarized {(3.0 ± 0.5) × 1014 W cm−2, γ  0.57, 3.2Up ~ 57 eV}, and circularly polarized {(1.5 ± 0.3) × 1014 W cm−2, γ  0.8, 3.2Up ~ 28 eV}, but not for chirped pulses {(0.56 ± 0.10) × 1014 W cm−2, γ  1.3, 3.2Up  11 eV}, for the latter MPI followed by fragmentation is anticipated given that γ > 1.

The metastable dications C2H62+; C2H52+; C2H42+ and C2H32+ result from NSDI based on the experimentally measured ellipticity dependence (see Supplementary Fig. 4), the observation that C2H62+ exhibits the tale tell enhancement ‘knee’ of Y(A++)/Y(A+) under linear and circularly polarized pulses9. The observation of high-order ATI of ethane; and supported by the observation of NSDI from a number of saturated and unsaturated hydrocarbons9,11,12. The sigma bonding electron is responsible for the initial tunnel ionization, leaving a weakly bound M+, evidenced by its lower probability compared to C2H5+ and C2H4+, this results in a more pronounced yield of dications resulting from NSDI6,8.

Early work on electron recollision in noble atoms led to the conclusion that field ellipticity prevents electron recollision5,35. It was assumed that NSDI, and the characteristic knee structure in the Y(A++)/Y(A+) ratio associated with it, would similarly be suppressed in other systems exposed to elliptically-polarized fields. This assertion was contradicted soon after by the observation of NSDI under circularly polarized light in Mg atoms7, and then by similar observations in NO, O236, and ethane9. These observations are consistent with theoretical models and predictions37. The parameters required for an atom or molecule to exhibit NSDI have been summarized as a set of scaling laws and the presentation of a ‘phase diagram’ where it becomes clear that low ionization potential (below 13 eV for 800 nm laser pulses), and high peak intensity from relatively short pulses (20–50 fs) ensures the recolliding electron has sufficient energy to release the second electron38,39,40. Based on scaling laws proposed for photoelectrons, it has been possible to show the conditions under which one expects NSDI by circularly polarized light/pulses even for Ar and He atoms40. This conclusion is consistent with long recollision periodic orbits that are responsible for NSDI driven by circularly polarized light37.

Results from shaped laser pulses

Figure 3 shows the corresponding portion of the mass spectrum for fragment ions of ethane with m/Q between 13 and 15.5 obtained with spectrally shaped laser pulses. The results shown in Fig. 3a have a positive ¾π phase step in either the high energy (blue line) or low energy (red-dashed line) portion of the spectrum. Notice that the yield of doubly charged fragments is greatly enhanced (suppressed) when the phase step is in the higher (lower) energy portion of the laser spectrum. The situation is reversed when a negative ¾π phase step is used, as shown in Fig. 3b. Note that the pairs of shaped pulses have identical spectra, energy, polarization, and peak intensity. The only difference is the location of the spectral-phase step. The reason why pulse shaping, as being carried out here, leads to such different strong-field ionization processes can be explained by considering that the frequency becomes time dependent, as described below.

Fig. 3: Portion of the mass spectrum of ethane in the region showing the fragment ions with m/Q between 13 and 15.5 obtained using shaped pulses.
figure 3

a Laser spectra for two different shaped pulses with a positive ¾π phase step in either the high energy (blue line) or low energy (red-dashed line) portion of the laser spectrum (see inset diagrams showing their location at ±σ/2 relative to ω0). b Spectra for negative ¾π phase step in either the high energy (blue line) or low energy (red-dashed line) portion of the laser spectrum. Notice the phase step serves to enhance or suppress the yield of doubly charged fragment ions associated with non-sequential double ionization, and that the effect is reversed upon changing the sign of the phase step. The four spectra, shown as measured (without normalization), were obtained within the span of 30 min under identical conditions except for the frequency where the phase step is located.

Simulation of the experimental data

A sharp π phase step in the frequency domain results in a phase jump in the time domain that breaks the pulse into two pulses. For a given electric field in the frequency domain, given by

$$E(\omega ) \equiv A(\omega )e^{i\varphi (\omega )} \propto \sqrt {I(\omega )} \;e^{i\varphi (\omega )},$$
(1)

where A(ω) is the amplitude and φ(ω) is the spectral phase, the shaped pulse in the time domain is obtained via the Fourier transformation

$$E(t) \equiv \frac{1}{{2\pi }} \int _{ - \infty }^{ + \infty } {\sqrt {I(\omega )} \;e^{i\varphi (\omega )}e^{ - i\omega\, t}\mathrm{d}\omega } .$$
(2)

Equations (1) and (2) can be used to obtain a time dependent description of the pulse I(t) = |E(t)|2. The temporal shape of the pulse, given by I(t), is identical for pairs of pulses with the same phase step sign and for which the position of the phase step in the frequency domain is detuned towards high/low energy symmetrically. We confirmed the shaped pulses were identical by comparing their second harmonic spectrum (Supplementary Fig. 5). Henceforth, we shall refer to these as “shaped-pulse pairs” that allow us to exclusively explore phase effects given that the pulses in the time domain are identical (see Supplementary Fig. 6). These types of pulses, with a positive or negative ¾π phase step, were employed in connection with control of the strong-field fragmentation of methanol41. The time-dependent intensity I(t) of the pair of shaped pulses used to obtain the TOF spectra shown in Fig. 3a is identical and can be described as a small shoulder preceding the main pulse. Clearly, the differences observed in TOF spectra shown in Fig. 3a indicate that I(t) alone is not responsible for the observed changes in fragment-ion yields. The position, magnitude, and sign of the phase step are the determining factors for the observed control over ionization. The same conclusions apply for the data shown in Fig. 3b, except that in this case I(t) is reversed in time when compared to the pulses in Fig. 3a. For Fig. 3b,I(t) can be described as a small shoulder following the main pulse.

Based on the aforementioned fifth- or sixth-order dependence of HHG on frequency42,43, and the observed enhancement of NSDI by higher frequencies25,44, we can assume that the difference in the yield of dications, observed here as an increase in the yield of doubly charged fragment ions, depends on the frequency of the field at the time of tunnel ionization. To simulate our data, we multiply the position of the phase step, expressed as ω − ω0, where ω0 is the central frequency of the pulse, by the spectrum of the pulse I(ω) = |E(ω)|2 elevated to the nth power given the nonlinearity of the process, where n is an odd number to maintain the sign dependence. Then, empirically we find that n= 5 fits the data better than n= 7. The fifth-order dependence is in agreement with the observed dependence of the ion yield on laser pulse intensity (linearly polarized) observed for Xe+ and C60+45. The product (ω − ω0)I5(ω − ω0) is plotted in Fig. 4a, b (solid line) as a function of the phase step position together with the relative difference (defined below) in the yield of ions. The ‘instantaneous frequency’46, describes how the frequency of the pulse changes as a function of time

$$\omega _{\mathrm{inst}}(t) \equiv \omega _0 + \frac{{{\mathrm{d}}\{{\mathrm{{Im}}} [\ln [E(t)]]\} }}{{{\mathrm{d}}t}},$$
(3)

where ω0 is the carrier frequency of the pulse, and the expression in the curly brackets corresponds to the time dependent phase of the pulse. This instantaneous frequency combines both frequency and time domains and can be unreliable when sharp changes in either frequency or time occur. Nevertheless, we calculated ωinst(t) when I(t) approaches its maximum amplitude. We then multiplied it by I5(ω − ω0) and obtained a function that is very similar to that plotted in Fig. 4a, allowing us to justify simply using the position of the phase step as a simple phenomenological model of the observed results.

Fig. 4: Dependence of ion yields on the position of a ¾ π phase step.
figure 4

a Relative difference between measured ion yields of C2H42+, C2H52+, and H3+, which are associated with non-sequential double ionization, as well as CH2+ Coulomb explosion fragment. b Relative difference between measured yield for doubly charged He and third harmonic generation (THG) photon yield in air for pulses similarly shaped as for the data in panel (a). The thin black line shown in both panels (a) and (b) corresponds to our phenomenological model, i.e. (ω − ω0)I5(ω − ω0), and the measured relative difference is shown “as measured” without normalization. Error bars, defined as the standard deviation of integrated ion count, are within the symbol size.

In order to test this phenomenological model, we acquired TOF spectra as a function of the position of a ¾π phase step with respect to the center of the spectrum ω − ω0, for both positive and negative phase steps. The results were plotted as the relative difference, a quantity that expresses the difference between two ion yields normalized by the average of the two, specifically (Y+ − Y)/[(Y+ + Y)/2], where the superscript sign indicates that the yield is measured for either a positive or negative phase step. For example, when the yield of some ion for a positive phase is three times greater than the yield for a negative phase at the same phase position, then the relative difference at that position is 1.0.

We find that the relative difference for doubly charged, C2H42+ and C2H52+, fragment ions is in very good agreement with the model (thin black line). Interestingly, H3+, a minor fragment ion resulting from double ionization of ethane30,47, follows the same trend, see Fig. 4a. In alcohols, NSDI most likely involves the non-bonding electrons at the oxygen atom, causes a reduction in electron density from the nearest carbon atom and leads to formation of neutral H2, which roams and extracts an additional proton to produce H3+48,49. Apparently, a similar process occurs in ethane. This observation is consistent with experiments on CH3CD3 for which the yields of H2+ and D2+ amounted to 93% of the diatomic ions, and the yields of H2D+ and HD2+ amounts to 65% of the triatomic ions and the rest to H3+ and D3+30. The yield for Coulomb explosion CH2+ ions, also shown in Fig. 3a, identified as shoulders flanking the m/Q = 14 peak, are found to have a weaker dependence on the phase step. Note that in Fig. 3a, b the amplitude of the CH2+ fragment ion is less affected by the position or sign of the phase step than the C2H42+ ions.

The agreement between the model and experiment in Fig. 4 suggests that the control mechanism is related to the position of the phase step, which causes a change in the instantaneous frequency and hence the period of the optical cycle of the pulse at the time of tunnel ionization. We observe maximum differences for step positions ±0.015 eV relative to the photon energy 1.55 eV. The pulses used are well described by a Gaussian function and the step position corresponds to ±σ/2 relative to ω0, where the intensity of the pulse is reduced by 1/e at ±σ. Because in each pulse pair we compare symmetric blue and red shifted frequency positions, the change in frequency is ~0.03 eV. While this difference seems to be small, as pointed above, recollision, NSDI, and HHG have been found to depend on the frequency of the pulse to the fifth or sixth power23,24,25. In terms of the frequency difference, we would anticipate a ~23% relative difference in the yield I6(ωblue) − I6(ωred) divided by the average of those quantities. This estimate is much smaller than the experimentally achieved control factor of three. Possible reasons for the much greater control achieved include the fact that multiple returning orbits of the electron can cause NSDI, and the experiments are carried out near the threshold for highly nonlinear processes; requiring 11.5 eV for first ionization (equivalent to 8 photons), and 32 ± 2 eV (equivalent to 21 photons) for double ionization. Note that the appearance energy of doubly charged fragment ions and H3+ fragment ions following ethane double ionization is 33 ± 1 eV (equivalent to 21 photons)28,29. Preliminary one-dimensional kinetic energy calculations for an unbound electron in the presence of the shaped laser field based on classical mechanics were carried out. An electron with mass m at rest interacting with a strong laser field E(t) experiences a force that is proportional to the electric field. Therefore, the kinetic energy of the electron at time τ is given by \(E_{{{\mathrm{kin}}}}(\tau ) = (1/2)m( { - {\smallint }_{t_i}^\tau E(t){{{\mathrm{d}}t}}} )^2\). The results from those calculations are shown in Supplementary Fig. 7. We hope that our results will motivate theorists to simulate our observations with more advanced methods.

Results from shaped circularly polarized laser pulses

Based on recent simulations aimed at understanding the observed NSDI from Mg atoms7,10,37,50,51,52,53,54,55 and some molecules45,56,57,58,59,60,61, circularly polarized pulses are associated with orbits having longer recollision times. These findings allow us to test our hypothesis that the observed control over the yield of doubly charged ions is caused by the instantaneous frequency at the time of ionization. Although longer orbits will lead to greater electron wavepacket dispersion, for a given orbit, higher frequencies imply shorter periods and less time for electron wavepacket dispersion. Therefore, longer recollision orbits should be even more sensitive to changes in instantaneous frequency of the shaped pulses.

Figure 5 shows the mass spectrum of ethane fragment ions with m/Q between 13 and 15.5 obtained with spectrally shaped laser pulses that were circularly polarized. A reduction in overall yield of all ions was found upon changing the polarization of the pulses given that the laser pulse energy was kept constant, and thus the peak intensity was reduced by a factor of 2 (recall that the peak field for circularly polarized pulses is decreased by √2 compared to linearly polarized pulses). The yield of dications for circularly polarized light was only 10% of that found for linearly polarized light. The ion yield dependence on polarization ellipticity is included in (Supplementary Fig. 4). We focus our analysis on the difference (blue/red lines) observed for the doubly charged C2H42+ and C2H52+, ions which is ~3 and >10, respectively, as determined by integrating the respective peaks in the TOF spectra. These ratios are contrasted with those observed for linearly polarized pulses, which are ~2 and ~6, respectively. The ~50% greater contrast observed for circularly polarized light seems to be consistent with the hypothesis that longer orbits are more sensitive to changes in the driving frequency. Note, these ratios are obtained from Figs. 3 and 5, and are not the relative differences shown in Fig. 4, which compare positive and negative phases at a given phase-step position. The lower contrast observed for C2H42+ indicates that this peak has a contribution of CH2+ fragment ions paired with neutral counterparts30. We note that, as in Fig. 3, the yield for the CH2+ Coulomb explosion channel is less affected by pulse shaping. We find that in the case of ethane, Coulomb explosion fragment ions are more prevalent than doubly charged fragments when linearly chirped pulses are used, see Fig. 2c.

Fig. 5: Portion of the mass spectrum of ethane in region showing the fragment ions with m/Q between 13 and 15.5 obtained using shaped circularly polarized pulses.
figure 5

a Results for a positive ¾π phase step in either the high energy (blue line) or low energy (red-dashed line) portion of the laser spectrum (see inset diagrams showing their location at ±σ/2 relative to ω0). b Spectra for negative ¾π phase step in either the high energy (blue line) or low energy (red-dashed line) portion of the laser spectrum. Notice the phase step serves to enhance or suppress doubly charged fragment ions associated with non-sequential double ionization, and that the effect is reversed upon changing the sign of the phase step. The four spectra were obtained under identical conditions, except for the phase step location. The data is shown as measured without normalization.

Considerations for future experiments

The degree of control reported here, a factor of three to six (corresponding to relative differences of 1–1.4) in the yield of doubly charged fragments resulting from pulse shaping, warrants extensive scrutiny. The explanation given is based on the chirp rate at the time of double ionization, albeit taking into consideration multiple electron orbits being accelerated by a few optical cycles. We tested the generality of the recollision control obtained by ¾π phase step shaped pulses and found that the yield of He2+ and of third harmonic generation in air have a similar dependence on the observed yield of ethane dications found here, see Fig. 3(b). Furthermore, our findings are consistent with enhanced H3+ formation from methanol observed under similar experimental conditions41. Further confirmation using shorter pulses, where the greater bandwidth will allow a greater frequency change, will be valuable. Electron correlation measurements62 will help provide a definite assignment of the NSDI process; in particular, if it is associated with recollision excitation followed by double ionization (RESI) or recollision-induced ionization (RII)13,62. Future experiments capable of detecting the momenta of the two electrons and the associated ions in coincidence, may help further identify NSDI dissociative pathways. Such a task is very demanding given that in some of the channels of interest there is a neural fragment or overlapping competing processes. Moreover, unlike photodissociation experiments carried out with single photon, the amount of energy (i.e., number of photons) absorbed by the molecule from the strong field used in this experiment is unknown. Simulations based on simpler systems, such as helium could also shed light on the role of the shaped pulses. Finally, other pulse shaping approaches may yield even greater changes in frequency while maintaining short-pulse duration, and therefore may lead to even higher level of control. For example, we have found that a very small (almost unmeasurable) chirp, causing <1.005 times pulse broadening, affects the control observed (see Supplementary Fig. 8). Comparison with recent methods aimed at tailoring the excitation field via the addition of other frequencies will also be instructive63.

Discussion

We found that scanning the position of a ¾π phase step across the spectrum of an intense femtosecond laser pulse leads to significant (factors of 3–6) changes in the yield of doubly charged fragment ions, as well as fragment ions associated with NSDI. A phenomenological model based on the position of the phase step relative to the central frequency of the pulse was found to be in excellent agreement with the experimental results. We found the maximum degree of control, enhancement/suppression of doubly charged species, is observed when the phase step is near ±σ/2 of the Gaussian laser spectrum. We experimentally confirmed the general ability of ¾π phase step shaped pulses to control electron recollision by quantifying the yield of metastable dications from ethane, as well as for double ionization of helium and third harmonic generation in air. Finally, we found that the degree of control increased for ethane under circularly polarized laser pulses. All processes that depend on electron recollision, such as HHG, generation of attosecond pulses via HHG, ATI, attochirp measurements, and attosecond clocking, should thus be controllable by tailoring the phase of the incident femtosecond laser pulses as indicated here without the need for adding pulses with different frequencies.

We predict that the magnitude of the control reported here using shaped pulses will increase as the duration of the laser pulses decreases. Methods to introduce a phase step causing minimum or no loss64 of pulse energy could represent significant enhancements in harmonic generation intensities, modest gains in cutoff energies, and significant increases in ATI yields. In terms of molecular processes, this type of pulse shaping could enable the double ionization of non-bonding electron pairs within a molecule, given the fact that linearly polarized light will favor electron recollision back to the orbital from which the electron tunnel ionized17,65. This degree of localized double ionization, as opposed to ionization out of delocalized valence orbitals, involving for example oxygen and nitrogen atoms, could allow unprecedented control of chemical reactivity, such as being able to cleave strong chemical bonds and not neighboring weaker bonds66.

Methods

Experiments are carried out employing a Wiley-McLaren TOF mass spectrometer operated at 10−7 Torr. The laser pulses are generated by a 1 kHz titanium sapphire chirped pulse amplified regenerative amplifier and attenuated to achieve the desired intensities (between 0.2 and 0.5 mJ per pulse). The pulses are first compressed to within 99.9% of the theoretical transform limit (40 fs), eliminating higher order dispersion, using the multiphoton intrapulse interference phase scan (MIIPS) method67,68, and then shaped using a commercial pulse shaper (model MIIPS-HD, Biophotonic Solutions, now part of IPG Photonics). Achieving pairs of shaped pulses that agree with theoretical predictions based on the second harmonic spectrum of the shaped pulses (see Supplementary Fig. 4)69 required using TL pulses with <40-fs2 of chirp and no measurable high-order dispersion. The compressed pulses are focused into the chamber by an achromatic 300 mm focal length lens. Peak intensities are calibrated in situ by comparing the ratio between singly and doubly charged oxygen, nitrogen, and argon with experimental values reported in the literature36,70. The reported peak intensities are accurate within ~20%. Ions formed by the laser pulse are detected using a microchannel plate, the resulting TOF spectrum was then digitized. Background electronic noise signature from the laser was acquired while the laser was blocked and subtracted from the raw time-of flight transients. Ion yields are obtained by integrating the appropriate peaks in the TOF spectrum. The mass spectra in Figs. 2, 3 and 5 are averages of 100,000 pulses, and the mass to charge range shown is obtained from 205 TOF points which are converted to m/Q. The data points in Fig. 4 represent averaged integrated intensity for each ion shown from 5000 laser pulses for each phase step position. Each such scan consisted of 200 phase step positions. To avoid overlapping of multiple traces, only half of the points are shown after 3-point smoothing of the original scans (for H3+ we only show one quarter of the points obtained). The results for He were obtained using 40 fs TL pulses with a peak intensity of 6 × 1014 W cm−2.