Zero-energy photoelectric effect

Sajad Azizi, Ulf Saalmann and Jan M. Rost Max-Planck-Institut für Physik komplexer Systeme, Nöthnitzer Str. 38, 01187 Dresden, Germany
(July 23, 2024)
Abstract

We predict a near-threshold (“zero energy”) peak in multi-photon ionization for a dynamical regime where the photon frequency is large compared to the binding energy of the electron. The peak position does not depend on the laser frequency, but on the binding energy and the pulse duration. The effect originates from the fact that bound-continuum dipole transitions are stronger than continuum-continuum ones. To clearly observe this zero-energy photoelectric effect, the spectral width of the laser pulse should be comparable to the binding energy of the ionized orbital, and the second ionization potential should be larger than the photon energy. This suggests negative ions as ideal candidates for corresponding experiments.

A threshold for an observable A𝐴Aitalic_A indicates the transition of the system from one regime to another upon change of the relevant parameter ε𝜀\varepsilonitalic_ε and provides therefore important information about the system. Often, the observable changes near threshold at ε0subscript𝜀0\varepsilon_{0}italic_ε start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT with a certain power of the parameter, i.e., A(ε+ε0)(ε/ε01)αproportional-to𝐴𝜀subscript𝜀0superscript𝜀subscript𝜀01𝛼A(\varepsilon{\to}{+}\varepsilon_{0})\propto\left(\varepsilon/\varepsilon_{0}{% -}1\right)^{\alpha}italic_A ( italic_ε → + italic_ε start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ∝ ( italic_ε / italic_ε start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - 1 ) start_POSTSUPERSCRIPT italic_α end_POSTSUPERSCRIPT. Probably best known are thermodynamical variables near phase transitions [1], but also quantum critical points in condensed matter [2] or fragmentation/ionization thresholds in atomic and molecular physics [3] are examples. In the latter case the so-called Wigner threshold law [4] for fragmentation cross sections of particles under short-range forces is a universal property, similar as its counterpart for long-range (Coulomb) forces, the Wannier law [5, 6].

The Wigner law has been verified in fragmentation scenarios including multi-photon detachment of a negative ions [7], which were a popular target for theoretical considerations [8] regarding above-threshold ionization (ATI). However, so far it has gone unnoticed that the combination of the Wigner threshold behavior with short intense pulses can give rise to a peculiar zero-energy photoelectric effect (ZEPE) with a characteristic maximum, whose position at very low photo-electron energy does not depend on the photon energy, but on the duration of the (short) laser pulse. That short intense pulses can lead to unusual electron dynamics has been pointed out in the context of non-adiabatic photo-ionization [9, 10] where short pulses can even be used for coherent control [11].

The ZEPE effect requires pulses short enough such that their spectral bandwidth ΔEΔ𝐸\Delta Eroman_Δ italic_E is larger than the binding energy EEAsubscript𝐸EAE_{\rm EA}italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT of the detached electron, typically of the order of a few femtoseconds. This was clearly not the case in the ATI detachment experiments with 800 nm light. ZEPE is enabled by a two-photon process, where a bound electron absorbs a photon and emits a photon, ending at the same (bound) energy as it started from. Therefore, this mechanism does typically not contribute to the photo-electron spectrum. Only, if the laser pulse is short enough that the spectral peak at the binding energy “leaks” into the continuum, see Fig. 1, this process becomes visible in form of photo-electrons and can lead, together with the Wigner power law, to a pronounced ZEPE maximum.

Negative ions are an ideal target to clearly identify the ZEPE, since they combine a low electron affinity (EA), i.e., the ionization potential of the negative ion, with a typically large gap to the next ionization potential (IP), i.e., the ionization potential of the (neutral) atom. Therefore, a relatively large photon frequency, still fulfilling EEA<ω<EIPsubscript𝐸EA𝜔subscript𝐸IPE_{\rm EA}<\omega<E_{\rm IP}italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT < italic_ω < italic_E start_POSTSUBSCRIPT roman_IP end_POSTSUBSCRIPT, can be chosen, which on the one hand energetically prevents ionization of more deeply bound electrons and on the other hand lets the first ATI peak appear at relatively high energy, keeping an energy interval just above threshold pristine for a clean signature of ZEPE.

Refer to caption
Figure 1: Multi-photon detachment of a an electron, weakly bound by electron affinity EEAsubscript𝐸EAE_{\rm EA}italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT to an ion, will give rise to peaks in the photo-electron spectrum at energies E=EEA+nω𝐸subscript𝐸EA𝑛𝜔E\,{=}\,{-}E_{\rm EA}{+}n\omegaitalic_E = - italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT + italic_n italic_ω, shown here for n=1,2𝑛12n{=}1,2italic_n = 1 , 2. However, following absorption of the 1st photon it is more likely that a 2nd photon is not absorbed but emitted (ωPlanck-constant-over-2-pi𝜔{-}\hbar\omega- roman_ℏ italic_ω) since the dipole for this transition is much stronger. This process with net-zero energy absorption becomes only visible in the spectrum for very short pulses (blue) leading to a broad peak at EEAsubscript𝐸EA-E_{\rm EA}- italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT whose tail reaches into the continuum (lightblue-shaded), but not for longer pulses (red).

We will demonstrate the effect in the following with negative ions of hydrogen and oxygen in a simple and transparent fashion. To this end, we use an effective potential for the bound electron of the negative ion which is designed to reproduce the EA well [12, 13],

V(r)=Zrb1+c[exp(r/r0)1],𝑉𝑟𝑍𝑟𝑏1𝑐delimited-[]𝑟subscript𝑟01V(r)=-\frac{Z}{r}\frac{b}{1+c\,[\exp(r/r_{0}){-}1]}\,,italic_V ( italic_r ) = - divide start_ARG italic_Z end_ARG start_ARG italic_r end_ARG divide start_ARG italic_b end_ARG start_ARG 1 + italic_c [ roman_exp ( start_ARG italic_r / italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) - 1 ] end_ARG , (1)

where Z𝑍Zitalic_Z is the nuclear charge, Z= 1𝑍1Z\,{=}\,1italic_Z = 1 and Z= 8𝑍8Z\,{=}\,8italic_Z = 8 for hydrogen and oxygen, respectively. We use atomic units unless stated otherwise. With parameter values b= 1.1𝑏1.1b\,{=}\,1.1italic_b = 1.1, c= 1𝑐1c\,{=}\,1italic_c = 1, and r0= 0.5292subscript𝑟00.5292r_{0}\,{=}\,0.5292italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.5292 Å the computed E1s=0.75subscript𝐸1s0.75E_{\rm 1s}\,{=}\,{-}0.75italic_E start_POSTSUBSCRIPT 1 roman_s end_POSTSUBSCRIPT = - 0.75 eV for H- closely matches the experimental value [14]. Similarly, with b= 1𝑏1b\,{=}\,1italic_b = 1, c= 1.9607𝑐1.9607c\,{=}\,1.9607italic_c = 1.9607, and r0= 0.4689subscript𝑟00.4689r_{0}\,{=}\,0.4689italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.4689 Å, the computed E2p=1.464subscript𝐸2p1.464E_{\rm 2p}\,{=}\,{-}1.464italic_E start_POSTSUBSCRIPT 2 roman_p end_POSTSUBSCRIPT = - 1.464 eV for O- is in good agreement with the experimental value of 1.4611.461-1.461- 1.461 eV [15].

With the potential (1) the Hamilton operator reads

H^(t)=𝒑^2/2+V(r)+𝒑^𝒆zA0g(t)cos(ωt),^𝐻𝑡superscript^𝒑22𝑉𝑟^𝒑subscript𝒆𝑧subscript𝐴0𝑔𝑡𝜔𝑡\widehat{H}(t)=\widehat{\bm{p}}{\,}^{2}\!/2+V(r)+\widehat{\bm{p}}\cdot\bm{e}_{% z}\,A_{0}g(t)\cos(\omega t),over^ start_ARG italic_H end_ARG ( italic_t ) = over^ start_ARG bold_italic_p end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 2 + italic_V ( italic_r ) + over^ start_ARG bold_italic_p end_ARG ⋅ bold_italic_e start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_A start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_g ( italic_t ) roman_cos ( start_ARG italic_ω italic_t end_ARG ) , (2)

where g(t)=exp(t2/T2)𝑔𝑡superscript𝑡2superscript𝑇2g(t)\,{=}\,{\exp}({-}t^{2}/T^{2})italic_g ( italic_t ) = roman_exp ( - italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_T start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) is the laser-pulse envelope with the pulse duration T=Tfwhm/2ln2𝑇subscript𝑇fwhm22T\,{=}\,T_{\rm fwhm}/\sqrt{2\ln\!2}italic_T = italic_T start_POSTSUBSCRIPT roman_fwhm end_POSTSUBSCRIPT / square-root start_ARG 2 roman_ln 2 end_ARG, linearly polarized along the z𝑧zitalic_z-axis and dipole coupled in velocity gauge. In single-active-electron approximation (2), the propagation of the 3D time-dependent Schrödinger equation (TDSE), for processes discussed here, is straightforward, since the number of photons involved is moderate. What is challenging, however, is the numerically accurate description of the very small detachment probabilities close to threshold. To this end we calculate eigenstates and dipole matrix elements for the potentials given in (1) for = 0max 0subscriptmax\ell\,{=}\,0\ldots\ell_{\rm max}roman_ℓ = 0 … roman_ℓ start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT in a box r= 0R𝑟 0𝑅r\,{=}\,0\ldots Ritalic_r = 0 … italic_R up to a cutoff energy Ecutsubscript𝐸cutE_{\rm cut}italic_E start_POSTSUBSCRIPT roman_cut end_POSTSUBSCRIPT. The parameters used are max= 4subscriptmax4\ell_{\rm max}\,{=}\,4roman_ℓ start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT = 4, R= 3×103𝑅3superscript103R\,{=}\,3{\times}10^{3}italic_R = 3 × 10 start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT Å and Ecut= 3subscript𝐸cut3E_{\rm cut}\,{=}\,3italic_E start_POSTSUBSCRIPT roman_cut end_POSTSUBSCRIPT = 3 keV. The TDSE is propagated numerically with this field-free basis. From the final amplitudes obtained we calculate the photo-electron spectrum. Attaching to every eigenstate a normalized Gaussian, whose height is the absolute square of the amplitude and whose width is given by the level spacing at the corresponding eigenenergy, produces a continuous spectrum. Thereby, the density of states and the detachment process are correctly included.

Refer to caption
Figure 2: Spectra from H- for Tfwhm= 1subscript𝑇fwhm1T_{\rm fwhm}\,{=}\,1italic_T start_POSTSUBSCRIPT roman_fwhm end_POSTSUBSCRIPT = 1 fs, I= 1014𝐼superscript1014I\,{=}\,10^{14}italic_I = 10 start_POSTSUPERSCRIPT 14 end_POSTSUPERSCRIPTW/cm2 and 4 different photon frequencies (ω= 3,6,9,12𝜔36912\omega\,{=}\,3,6,9,12italic_ω = 3 , 6 , 9 , 12 eV). Note, while higher-order peaks (at energies E 6greater-than-or-equivalent-to𝐸6E\,{\gtrsim}\,6italic_E ≳ 6 eV) shift with increasing ω𝜔\omegaitalic_ω to higher energies, the zero-peak position (Es 0.3subscript𝐸s0.3E_{\rm s}\,{\approx}\,0.3italic_E start_POSTSUBSCRIPT roman_s end_POSTSUBSCRIPT ≈ 0.3 eV and Ed 1.4subscript𝐸d1.4E_{\rm d}\,{\approx}\,1.4italic_E start_POSTSUBSCRIPT roman_d end_POSTSUBSCRIPT ≈ 1.4 eV) is independent of ω𝜔\omegaitalic_ω. Spectra are scaled such that they agree for E= 0.1𝐸0.1E\,{=}\,0.1italic_E = 0.1 eV.

Figure 2 shows the resulting photo-electron spectra for H- exposed to a laser pulse of 1 fs duration and a peak intensity of I= 1014W/cm2𝐼superscript1014Wsuperscriptcm2I\,{=}\,10^{14}\,\mathrm{W/cm^{2}}italic_I = 10 start_POSTSUPERSCRIPT 14 end_POSTSUPERSCRIPT roman_W / roman_cm start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT for different photon frequencies. The unique feature of the zero-energy photo-effect is evident: There is a pronounced maximum close to threshold (Es0.3subscript𝐸s0.3E_{\rm s}{\approx}0.3italic_E start_POSTSUBSCRIPT roman_s end_POSTSUBSCRIPT ≈ 0.3 eV and Ed1.4subscript𝐸d1.4E_{\rm d}{\approx}1.4italic_E start_POSTSUBSCRIPT roman_d end_POSTSUBSCRIPT ≈ 1.4 eV) which is independent of the photon frequency, clearly different in shape from the two-photon ATI peaks which can be identified at the respective energies E= 2ωEEA𝐸2𝜔subscript𝐸EAE\,{=}\,2\omega\,{-}\,E_{\rm EA}italic_E = 2 italic_ω - italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT. If the ZEPE maximum is much smaller than and too close to the ATI peaks, typically for higher final partial waves >00\ell>0roman_ℓ > 0, then it may get buried under the rise to the two-photon ATI peak, which is the case for H{}^{-}{\to}start_FLOATSUPERSCRIPT - end_FLOATSUPERSCRIPT → d detachment with 3 eV photons, as visible in Fig. 2b.

Refer to caption
Figure 3: Photo-electron spectrum (colored solid lines) of H- for a Gaussian pulse with different pulse duration Tfwhm=1,2,3subscript𝑇fwhm123T_{\rm fwhm}=1,2,3italic_T start_POSTSUBSCRIPT roman_fwhm end_POSTSUBSCRIPT = 1 , 2 , 3 fs, carrier frequency ω=9𝜔9\omega=9\,italic_ω = 9eV, and intensity I=1014W/cm2𝐼superscript1014Wsuperscriptcm2I=10^{14}\,\mathrm{W/cm^{2}}italic_I = 10 start_POSTSUPERSCRIPT 14 end_POSTSUPERSCRIPT roman_W / roman_cm start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT for the s (a) and the d-channel (b), respectively, and fits (black dashed lines) according to Eq. (4).

To understand how ZEPE comes about qualitatively and quantitatively, we take a closer look at the near-threshold energy range with the spectra of Fig. 3 on a linear scale for different pulse durations T𝑇Titalic_T. One immediately sees that the spectra depend on T𝑇Titalic_T, as already anticipated. The Wigner threshold law [4] states that for break-up of two fragments under short-range forces (which is the case for electron detachment), the ionization probability near threshold is given by the available continuum states in momentum space, i.e.,

P(E0)dpp2(+1)δ(p2/2E)E+1/2,proportional-to𝑃𝐸0differential-d𝑝superscript𝑝21𝛿superscript𝑝22𝐸proportional-tosuperscript𝐸12P(E{\to}0)\propto\mbox{\large$\int$}\!\mathrm{d}p\>p^{2(\ell+1)}\delta(p^{2}\!% /2{-}E)\propto E^{\ell+1/2},italic_P ( italic_E → 0 ) ∝ ∫ roman_d italic_p italic_p start_POSTSUPERSCRIPT 2 ( roman_ℓ + 1 ) end_POSTSUPERSCRIPT italic_δ ( italic_p start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 2 - italic_E ) ∝ italic_E start_POSTSUPERSCRIPT roman_ℓ + 1 / 2 end_POSTSUPERSCRIPT , (3)

where \ellroman_ℓ is the angular momentum of the fragment pair. Following the intuition that ZEPE originates in the combination of the Wigner law and a Gaussian energy distribution induced by the short laser pulse, located at the binding energy EEAsubscript𝐸EA-E_{\rm EA}- italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT due to the two-photon zero-energy process, the detachment probability should be given by

PZEPE,β(E)superscriptsubscript𝑃ZEPE𝛽𝐸\displaystyle P_{\rm ZEPE}^{\ell,\beta}(E)italic_P start_POSTSUBSCRIPT roman_ZEPE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_ℓ , italic_β end_POSTSUPERSCRIPT ( italic_E ) =P(E)s,β(E/EEA),absentsubscript𝑃𝐸subscript𝑠𝛽𝐸subscript𝐸EA\displaystyle=P_{*}(E)\,s_{\ell,\beta}(E/E_{\rm EA}),= italic_P start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_E ) italic_s start_POSTSUBSCRIPT roman_ℓ , italic_β end_POSTSUBSCRIPT ( italic_E / italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT ) , (4)
P(E)subscript𝑃𝐸\displaystyle P_{*}(E)italic_P start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_E ) =[1+E/E]1absentsuperscriptdelimited-[]1𝐸subscript𝐸1\displaystyle=[1+E/E_{*}]^{-1}= [ 1 + italic_E / italic_E start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT (4a)
s,β(x)subscript𝑠𝛽𝑥\displaystyle s_{\ell,\beta}(x)italic_s start_POSTSUBSCRIPT roman_ℓ , italic_β end_POSTSUBSCRIPT ( italic_x ) =β4x+12exp(β2[x+1]2).absentsuperscript𝛽4superscript𝑥12superscript𝛽2superscriptdelimited-[]𝑥12\displaystyle=\beta^{4}x^{\ell+\frac{1}{2}}\exp(-\beta^{2}[x{+}1]^{2}).= italic_β start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_x start_POSTSUPERSCRIPT roman_ℓ + divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT roman_exp ( start_ARG - italic_β start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [ italic_x + 1 ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) . (4b)

The universal shape s,β(x)subscript𝑠𝛽𝑥s_{\ell,\beta}(x)italic_s start_POSTSUBSCRIPT roman_ℓ , italic_β end_POSTSUBSCRIPT ( italic_x ) of the zero-energy photo-electron spectrum is fully determined by the angular momentum \ellroman_ℓ of the photo-electron and β=EEA/ΔE𝛽subscript𝐸EAΔ𝐸\beta\,{=}\,E_{\rm EA}/\Delta Eitalic_β = italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT / roman_Δ italic_E the ratio of electron affinity EEAsubscript𝐸EAE_{\rm EA}italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT to spectral pulse width ΔEΔ𝐸\Delta Eroman_Δ italic_E. From the spectral representation of the pulse and the fact that ZEPE is a two-photon process follows ΔE= 2/TΔ𝐸2𝑇\Delta E\,{=}\,2/Troman_Δ italic_E = 2 / italic_T, which is confirmed by 2nd-order perturbation theory discussed below. Furthermore, P(E)subscript𝑃𝐸P_{*}(E)italic_P start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ( italic_E ) takes care of the slowly-varying background, being different for each ion.

Refer to caption
Figure 4: Photo-electron spectra for Hsuperscript𝐻H^{-}italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT (a,b) as in Fig. 3 but also for the two final channels of O- (c,d) in double-logarithmic scale. Note the different vertical scale for each panel. The gray line in panel (c) is for a lower intensity of I=1012𝐼superscript1012I{=}10^{12}italic_I = 10 start_POSTSUPERSCRIPT 12 end_POSTSUPERSCRIPTW/cm2 and therefore scaled by a factor of 104superscript10410^{4}10 start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT. The dotted lines in panel (d) are computed with a basis of states up to Ecut=30subscript𝐸cut30E_{\rm cut}{=}30italic_E start_POSTSUBSCRIPT roman_cut end_POSTSUBSCRIPT = 30 eV which would allow to describe 3-photon processes, the solid curves are for Ecut=3subscript𝐸cut3E_{\rm cut}{=}3italic_E start_POSTSUBSCRIPT roman_cut end_POSTSUBSCRIPT = 3 keV. Arrows mark the maxima according to Eq. (5).

The dashed lines in Fig. 3 show PZEPE,β(E)subscriptsuperscript𝑃𝛽ZEPE𝐸P^{\ell,\beta}_{\rm ZEPE}(E)italic_P start_POSTSUPERSCRIPT roman_ℓ , italic_β end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ZEPE end_POSTSUBSCRIPT ( italic_E ) according to Eq. (4). The only fit parameter Esubscript𝐸E_{\ast}italic_E start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT assumes the values E=0.896subscript𝐸0.896E_{\ast}{=}0.896italic_E start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT = 0.896 eV and E=0.684subscript𝐸0.684E_{\ast}{=}0.684italic_E start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT = 0.684 eV for the “s” and “p” photo-electrons of hydrogen and oxygen, respectively, while E=13.4subscript𝐸13.4E_{\ast}{=}13.4italic_E start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT = 13.4 eV and E>1000subscript𝐸1000E_{\ast}{>}1000italic_E start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT > 1000 eV for the electrons detached with angular momentum “d” and “f” reveal that the background is nearly constant in these cases over the energy interval considered. Obviously, Eq. (4) represents the numerical results very well, giving confidence in the interpretation and description of ZEPE provided. Deviations for O- at “longer” pulses can be attributed to higher-order effects occurring for intense pulses, since the agreement with the analytic description for a weaker pulse with I=1012𝐼superscript1012I{=}10^{12}italic_I = 10 start_POSTSUPERSCRIPT 12 end_POSTSUPERSCRIPTW/cm2, as shown by the gray line in Fig. 4c, is excellent.

Peculiar at second glance is the variation of the maximal detachment with the pulse length T𝑇Titalic_T or ΔEΔ𝐸\Delta Eroman_Δ italic_E, respectively. Note, that the highest maximum is achieved in Fig.3 with Tfwhm= 2subscript𝑇fwhm2T_{\rm fwhm}{=}\,2italic_T start_POSTSUBSCRIPT roman_fwhm end_POSTSUBSCRIPT = 2 fs for the =00\ell{=}0roman_ℓ = 0 case, while for all other \ellroman_ℓ the maxima seem to increase monotonically. To elucidate their behavior systematically, we determine the maximum of PZEPE,βsuperscriptsubscript𝑃ZEPE𝛽P_{\rm ZEPE}^{\ell,\beta}italic_P start_POSTSUBSCRIPT roman_ZEPE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_ℓ , italic_β end_POSTSUPERSCRIPT with respect to the two-dimensional parameter space {β,x}𝛽𝑥\{\beta,x\}{ italic_β , italic_x }, i. e., the (scaled) pulse width β=EEA/ΔE𝛽subscript𝐸EAΔ𝐸\beta\,{=}\,E_{\rm EA}/\Delta Eitalic_β = italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT / roman_Δ italic_E and the (scaled) excess energy x=E/EEA𝑥𝐸subscript𝐸EAx\,{=}\,E/E_{\rm EA}italic_x = italic_E / italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT. The corresponding optimal parameters are given by analytical but lengthy expressions. Inserting β(x)=2/(1+x)𝛽𝑥21𝑥\beta(x)\,{=}\,\sqrt{2}/(1{+}x)italic_β ( italic_x ) = square-root start_ARG 2 end_ARG / ( 1 + italic_x ), the solution of PZEPE,β/β=0superscriptsubscript𝑃ZEPE𝛽𝛽0\partial\!P_{\rm ZEPE}^{\ell,\beta}/\partial\!\beta{=}0∂ italic_P start_POSTSUBSCRIPT roman_ZEPE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_ℓ , italic_β end_POSTSUPERSCRIPT / ∂ italic_β = 0 into PZEPE,β/x=0superscriptsubscript𝑃ZEPE𝛽𝑥0\partial\!P_{\rm ZEPE}^{\ell,\beta}/\partial\!x{=}0∂ italic_P start_POSTSUBSCRIPT roman_ZEPE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_ℓ , italic_β end_POSTSUPERSCRIPT / ∂ italic_x = 0 gives an analytical but lengthy expression for βMax(x,)subscript𝛽Maxsuperscript𝑥\beta_{\rm Max}(x^{*},\ell)italic_β start_POSTSUBSCRIPT roman_Max end_POSTSUBSCRIPT ( italic_x start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT , roman_ℓ ), i. e. the pulse width that gives the highest peak of the detachment probability. Ignoring the background modification (4a) by taking the limit xsuperscript𝑥x^{*}{\to}\inftyitalic_x start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT → ∞ one obtains with the simpler expressions βMax=(72)/42subscript𝛽Max7242\beta_{\rm Max}\,{=}\,(7{-}2\ell)/{4\sqrt{2}}italic_β start_POSTSUBSCRIPT roman_Max end_POSTSUBSCRIPT = ( 7 - 2 roman_ℓ ) / 4 square-root start_ARG 2 end_ARG and xMax=(1+2)/(72)subscript𝑥Max1272x_{\rm Max}\,{=}\,(1{+}2\ell)/(7{-}2\ell)italic_x start_POSTSUBSCRIPT roman_Max end_POSTSUBSCRIPT = ( 1 + 2 roman_ℓ ) / ( 7 - 2 roman_ℓ ) directly the optimal pulse duration Tfwhm=8log2βMax/EEAsubscript𝑇fwhm82subscript𝛽Maxsubscript𝐸EAT_{\rm fwhm}{=}\sqrt{8\log 2}\,\beta_{\rm Max}/E_{\rm EA}italic_T start_POSTSUBSCRIPT roman_fwhm end_POSTSUBSCRIPT = square-root start_ARG 8 roman_log 2 end_ARG italic_β start_POSTSUBSCRIPT roman_Max end_POSTSUBSCRIPT / italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT. They values are 2.55 fs and 1.09 fs for H- (s- and d-channel) and 0.92 fs and 0.18 fs for O- (p- and f-channel), consistent with the data shown in Fig. 4. Within the same approximation one can determine the location of the maxima of PZEPE,βsuperscriptsubscript𝑃ZEPE𝛽P_{\rm ZEPE}^{\ell,\beta}italic_P start_POSTSUBSCRIPT roman_ZEPE end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_ℓ , italic_β end_POSTSUPERSCRIPT for a given pulse durations β𝛽\betaitalic_β. They read

Emax=EEA2(1+(2+1)/β21)subscript𝐸maxsubscript𝐸EA2121superscript𝛽21E_{\rm max}=\frac{E_{\rm EA}}{2}\big{(}\sqrt{1{+}(2\ell{+}1)/\beta^{2}}-1\big{)}italic_E start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT = divide start_ARG italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG ( square-root start_ARG 1 + ( 2 roman_ℓ + 1 ) / italic_β start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - 1 ) (5)

and are shown in Fig. 4 for the respective pulses with vertical arrows. Although for =00\ell{=}0roman_ℓ = 0 (H-) and =11\ell{=}1roman_ℓ = 1 (O-) the background is relevant, since Esubscript𝐸E_{*}italic_E start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT is in the energy range of interest, the values (5) explain the maxima there quite well.

Hence, the overview of ZEPE for H- and O- in Fig. 4, highlighting the influence of different partial waves =0,,303\ell=0,...,3roman_ℓ = 0 , … , 3 of the fragments reveal that essentially all features discussed so far can be identified and are confirmed. For completeness we estimate the number of ZEPE electrons per shot one could detect in an experiment,

Nexp=Nion[Iexp/1014W/cm2]2PT,subscript𝑁expsubscript𝑁ionsuperscriptdelimited-[]subscript𝐼expsuperscript1014Wsuperscriptcm22subscript𝑃𝑇N_{\rm exp}=N_{\rm ion}\big{[}I_{\rm exp}\big{/}10^{14}{\rm W/cm}^{2}\big{]}^{% 2}P_{T},italic_N start_POSTSUBSCRIPT roman_exp end_POSTSUBSCRIPT = italic_N start_POSTSUBSCRIPT roman_ion end_POSTSUBSCRIPT [ italic_I start_POSTSUBSCRIPT roman_exp end_POSTSUBSCRIPT / 10 start_POSTSUPERSCRIPT 14 end_POSTSUPERSCRIPT roman_W / roman_cm start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT , (6)

where Nionsubscript𝑁ionN_{\rm ion}italic_N start_POSTSUBSCRIPT roman_ion end_POSTSUBSCRIPT is the number of ions in the target volume, Iexpsubscript𝐼expI_{\rm exp}italic_I start_POSTSUBSCRIPT roman_exp end_POSTSUBSCRIPT the experimental laser intensity and PTsubscript𝑃𝑇P_{T}italic_P start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT is given in Tab. 1.

Table 1: Total ionization probability PTsubscript𝑃𝑇P_{T}italic_P start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT for H- and O- for accessible channels and three pulse durations Tfwhmsubscript𝑇fwhmT_{\rm fwhm}italic_T start_POSTSUBSCRIPT roman_fwhm end_POSTSUBSCRIPT.
H{}^{-}{\to}\,start_FLOATSUPERSCRIPT - end_FLOATSUPERSCRIPT →s H{}^{-}{\to}\,start_FLOATSUPERSCRIPT - end_FLOATSUPERSCRIPT →d O{}^{-}{\to}\,start_FLOATSUPERSCRIPT - end_FLOATSUPERSCRIPT →s O{}^{-}{\to}\,start_FLOATSUPERSCRIPT - end_FLOATSUPERSCRIPT →d
Tfwhmsubscript𝑇fwhmT_{\rm fwhm}italic_T start_POSTSUBSCRIPT roman_fwhm end_POSTSUBSCRIPT=1fs 8.9×1048.9superscript1048.9{\times}10^{-4}8.9 × 10 start_POSTSUPERSCRIPT - 4 end_POSTSUPERSCRIPT 1.5×1061.5superscript1061.5{\times}10^{-6}1.5 × 10 start_POSTSUPERSCRIPT - 6 end_POSTSUPERSCRIPT 1.8×1041.8superscript1041.8{\times}10^{-4}1.8 × 10 start_POSTSUPERSCRIPT - 4 end_POSTSUPERSCRIPT 7.1×1087.1superscript1087.1{\times}10^{-8}7.1 × 10 start_POSTSUPERSCRIPT - 8 end_POSTSUPERSCRIPT
2fs 4.4×1044.4superscript1044.4{\times}10^{-4}4.4 × 10 start_POSTSUPERSCRIPT - 4 end_POSTSUPERSCRIPT 9.1×1089.1superscript1089.1{\times}10^{-8}9.1 × 10 start_POSTSUPERSCRIPT - 8 end_POSTSUPERSCRIPT 3.1×1063.1superscript1063.1{\times}10^{-6}3.1 × 10 start_POSTSUPERSCRIPT - 6 end_POSTSUPERSCRIPT 9.7×10119.7superscript10119.7{\times}10^{-11}9.7 × 10 start_POSTSUPERSCRIPT - 11 end_POSTSUPERSCRIPT
3fs 1.1×1041.1superscript1041.1{\times}10^{-4}1.1 × 10 start_POSTSUPERSCRIPT - 4 end_POSTSUPERSCRIPT 6.4×1096.4superscript1096.4{\times}10^{-9}6.4 × 10 start_POSTSUPERSCRIPT - 9 end_POSTSUPERSCRIPT 3.9×1083.9superscript1083.9{\times}10^{-8}3.9 × 10 start_POSTSUPERSCRIPT - 8 end_POSTSUPERSCRIPT 6.8×10136.8superscript10136.8{\times}10^{-13}6.8 × 10 start_POSTSUPERSCRIPT - 13 end_POSTSUPERSCRIPT

The double-logarithmic scale in Fig. 4 was chosen to emphasize the threshold behavior and one can see that the quasi-analytical formula (4) performs in all cases extremely well compared to the numerical TDSE calculations, particularly with regard to the large dynamic ranges considered. In fact, the TDSE results in panel 4d contain a surprise: They appear not to be converged for the = 33\ell\,{=}\,3roman_ℓ = 3 spectrum of oxygen. Indeed, despite very small photo-electron energies of the order of 1 meV, numerically one has to include continuum electrons up to 3 keV to achieve convergence (see also Fig. 5). An explanation is provided by 2nd-order perturbation theory for the two-photon zero-energy process comprising the two different events of absorption followed by emission (η=``+"𝜂``"\eta{=}``+"italic_η = ` ` + "), typically much stronger than emission followed by absorption (η=``"𝜂``"\eta{=}``-"italic_η = ` ` - "). The ionization amplitude to energy E𝐸Eitalic_E for these processes reads [16]

a(E)𝑎𝐸\displaystyle a(E)italic_a ( italic_E ) =kdEkdkEA+dtA(t)ei[EEk]ttdtA(t)ei[Ek+EEA]tabsentsubscript𝑘subscript𝑑𝐸𝑘subscript𝑑𝑘EAsuperscriptsubscriptdifferential-d𝑡𝐴𝑡superscripteidelimited-[]𝐸subscript𝐸𝑘𝑡superscriptsubscript𝑡differential-dsuperscript𝑡𝐴superscript𝑡superscripteidelimited-[]subscript𝐸𝑘subscript𝐸EAsuperscript𝑡\displaystyle=\sum_{k}d_{Ek}d_{k\,\rm EA}\int_{-\infty}^{+\infty}\!\!\!\mathrm% {d}t\,A(t)\mathrm{e}^{{\mathrm{i}[E-E_{k}]t}}\int_{-\infty}^{t}\!\!\!\mathrm{d% }t^{\prime}\,A(t^{\prime})\mathrm{e}^{{\mathrm{i}[E_{k}+E_{\rm EA}]t^{\prime}}}= ∑ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_E italic_k end_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_k roman_EA end_POSTSUBSCRIPT ∫ start_POSTSUBSCRIPT - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + ∞ end_POSTSUPERSCRIPT roman_d italic_t italic_A ( italic_t ) roman_e start_POSTSUPERSCRIPT roman_i [ italic_E - italic_E start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ] italic_t end_POSTSUPERSCRIPT ∫ start_POSTSUBSCRIPT - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT roman_d italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_A ( italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) roman_e start_POSTSUPERSCRIPT roman_i [ italic_E start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT + italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT ] italic_t start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT
=π8T2A02kdEkdkEAη=±[e[[Δkη+EEA]2+[EΔkη]2]T2/42iπe[E+EEA]2T2/8F([2Δkη+EEAE]T/8)]absent𝜋8superscript𝑇2superscriptsubscript𝐴02subscript𝑘subscript𝑑𝐸𝑘subscript𝑑𝑘EAsubscript𝜂plus-or-minusdelimited-[]superscriptedelimited-[]superscriptdelimited-[]subscriptsuperscriptΔ𝜂𝑘subscript𝐸EA2superscriptdelimited-[]𝐸subscriptsuperscriptΔ𝜂𝑘2superscript𝑇242i𝜋superscriptesuperscriptdelimited-[]𝐸subscript𝐸EA2superscript𝑇28𝐹delimited-[]2subscriptsuperscriptΔ𝜂𝑘subscript𝐸EA𝐸𝑇8\displaystyle=\frac{\pi}{8}T^{2}A_{0}{\!}^{2}\sum_{k}d_{Ek}d_{k\,\rm EA}\,\sum% _{\eta=\pm}\Big{[}\mathrm{e}^{{-[[\Delta^{\eta}_{k}+E_{\rm EA}]^{2}+[E-\Delta^% {\eta}_{k}]^{2}]T^{2}\!/4}}-\tfrac{2\mathrm{i}}{\sqrt{\pi}}\mathrm{e}^{{-[E+E_% {\rm EA}]^{2}T^{2}\!/8}}\,F\big{(}[2\Delta^{\eta}_{k}{+}E_{\rm EA}{-}E]T/\sqrt% {8}\big{)}\Big{]}= divide start_ARG italic_π end_ARG start_ARG 8 end_ARG italic_T start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_A start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_E italic_k end_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_k roman_EA end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_η = ± end_POSTSUBSCRIPT [ roman_e start_POSTSUPERSCRIPT - [ [ roman_Δ start_POSTSUPERSCRIPT italic_η end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT + italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + [ italic_E - roman_Δ start_POSTSUPERSCRIPT italic_η end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] italic_T start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 4 end_POSTSUPERSCRIPT - divide start_ARG 2 roman_i end_ARG start_ARG square-root start_ARG italic_π end_ARG end_ARG roman_e start_POSTSUPERSCRIPT - [ italic_E + italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_T start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 8 end_POSTSUPERSCRIPT italic_F ( [ 2 roman_Δ start_POSTSUPERSCRIPT italic_η end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT + italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT - italic_E ] italic_T / square-root start_ARG 8 end_ARG ) ] (7)

with the dipole matrix elements djk𝒆zϕj|𝒅|ϕksubscript𝑑𝑗𝑘subscript𝒆𝑧quantum-operator-productsubscriptitalic-ϕ𝑗𝒅subscriptitalic-ϕ𝑘d_{jk}\,{\equiv}\,\bm{e}_{z}\,{\cdot}\,\langle\phi_{j}|\bm{d}|\phi_{k}\rangleitalic_d start_POSTSUBSCRIPT italic_j italic_k end_POSTSUBSCRIPT ≡ bold_italic_e start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ⋅ ⟨ italic_ϕ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT | bold_italic_d | italic_ϕ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ⟩ and the detunings Δk±EkωsubscriptsuperscriptΔplus-or-minus𝑘minus-or-plussubscript𝐸𝑘𝜔\Delta^{\pm}_{k}\,{\equiv}\,E_{k}{\mp}\omegaroman_Δ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ≡ italic_E start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ∓ italic_ω. We have not explicitly specified that |ϕ0ketsubscriptitalic-ϕ0|\phi_{0}\rangle| italic_ϕ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ⟩, |ϕkketsubscriptitalic-ϕ𝑘|\phi_{k}\rangle| italic_ϕ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ⟩, |ϕEketsubscriptitalic-ϕ𝐸|\phi_{E}\rangle| italic_ϕ start_POSTSUBSCRIPT italic_E end_POSTSUBSCRIPT ⟩ have different angular momenta \ellroman_ℓ, dictated by selection rules in the dipole matrix elements, and for simplicity we have assumed that all virtual states |ϕkketsubscriptitalic-ϕ𝑘|\phi_{k}\rangle| italic_ϕ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ⟩ are discrete in accordance with our numerical treatment of the continuum. Finally, due to the necessary cutoff at Ecutsubscript𝐸cutE_{\rm cut}italic_E start_POSTSUBSCRIPT roman_cut end_POSTSUBSCRIPT, the sum over virtual states k𝑘kitalic_k is finite. From the real part of (7) one sees that continuum states, which are resonant with initial one-photon absorption Δk= 0subscriptΔ𝑘 0\Delta_{k}\,{=}\,0roman_Δ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = 0, are dominant as well as the final energy equal to the initial energy, E=EEA𝐸subscript𝐸EAE\,{=}\,{-}E_{\rm EA}italic_E = - italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT, as expected for the zero-energy photoelectric effect.

Refer to caption
Figure 5: Photo-electron spectra in linear (left panel) and double-logarithmic (right) scale for the f-channel of O- from TDSE calculations (ω𝜔\omegaitalic_ω and I𝐼Iitalic_I as in Fig. 3, Tfwhm=1subscript𝑇fwhm1T_{\rm fwhm}{=}1italic_T start_POSTSUBSCRIPT roman_fwhm end_POSTSUBSCRIPT = 1 fs) with different cut-off energies Ecutsubscript𝐸cutE_{\rm cut}italic_E start_POSTSUBSCRIPT roman_cut end_POSTSUBSCRIPT, see text.

Coming back to the requirement of including very high energies Eksubscript𝐸𝑘E_{k}italic_E start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT in the calculation to reach convergence for small finite energy E𝐸Eitalic_E, we note that the Dawson function F𝐹Fitalic_F [17] in the imaginary part of Eq. (7) falls off very slowly F(x) 1/(2x)similar-to𝐹𝑥12𝑥F(x)\,{\sim}\,1/(2x)italic_F ( italic_x ) ∼ 1 / ( 2 italic_x ) for large x𝑥xitalic_x. That implies that non-resonant states EkEmuch-greater-thansubscript𝐸𝑘𝐸E_{k}\gg Eitalic_E start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ≫ italic_E can contribute significantly. Since the detachment probability (E/EEA)+1/2proportional-toabsentsuperscript𝐸subscript𝐸EA12{\propto}\,(E/E_{\rm EA})^{\ell+1/2}∝ ( italic_E / italic_E start_POSTSUBSCRIPT roman_EA end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT roman_ℓ + 1 / 2 end_POSTSUPERSCRIPT is very small near threshold, it must be determined with high absolute accuracy which explains why unexpectedly high virtual energies need to be taken into account. Note, that for the very same reason the two-photon process of first emitting and then absorbing a photon, in standard situations never considered, has a non-negligible contribution.

In summary, we have discussed a universal two-photon process optimally observable in negative ions, where the weakly-bound electron interacts with a short laser pulse of spectral width comparable to the electronic binding energy. Together with the typical Wigner power law of the electron-detachment cross section near threshold, a characteristic maximum forms at very-low electron energies, which does not dependent on the photon frequency, as accurate numerical calculations confirm. They require the inclusion of surprisingly energetic electronic continuum states to converge, highlighting subtle features of this unconventional process which also surface in 2nd-order time-dependent perturbation theory. Following physical intuition, we have derived an analytical form of the spectrum. It describes the relative peak height and location as well as the dependence on pulse length and electron affinity of this peculiar zero-energy photoelectric effect very well which will facilitate its experimental realization.

S. A. acknowledges discussion with Jonathan Dubois in the early stage of this work.

References

  • [1] H.E. Stanley, Introduction to phase transitions and critical phenomena. International series of monographs on physics. Oxford University Press 1971.
  • [2] S. Sachdev, Quantum phase transitions. Cambridge University Press 2000.
  • [3] J. M. Rost, Critical phenomena in atomic physics Physica E 9, 467 (2001).
  • [4] E. P. Wigner, On the behavior of cross sections near thresholds Phys. Rev. 73, 1002 (1948).
  • [5] G. H. Wannier, The threshold law for single ionization of atoms or ions by electrons Phys. Rev. 90, 817 (1953).
  • [6] J. M. Rost, Two-electron escape near threshold: A classical process? Phys. Rev. Lett. 72, 1998 (1994).
  • [7] R. Reichle, H. Helm, and I Yu. Kiyan, Photodetachment of H- in a strong infrared laser field Phys. Rev. Lett. 87, 243001 (2001).
  • [8] K. Krajewska, I. I. Fabrikant, and A. F. Starace, Threshold effects in strong-field detachment of H- and F-: Plateau enhancements and angular distribution variations Phys. Rev. A 74, 053407 (2006).
  • [9] K. Toyota, O. I. Tolstikhin, T. Morishita, and S. Watanabe, Slow electrons generated by intense high-frequency laser pulses Phys. Rev. Lett. 103, 153003 (2009).
  • [10] Q.-C. Ning, U. Saalmann, and J. M. Rost, Electron dynamics driven by light-pulse derivatives Phys. Rev. Lett. 120, 033203 (2018).
  • [11] S. Azizi, U. Saalmann, and Jan M. Rost, Non-adiabatic ionization with tailored laser pulses J. Phys. B 54, 134001 (2021).
  • [12] P. S. Ganas, J. D. Talman, and A. E. S. Green, Independent-particle models for light negative atomic ions Phys. Rev. A 22, 336 (1980).
  • [13] R. Shakeshaft and X. Tang, Integral-equation approach to multiphoton ionization by intense fields. II. Application to H and H- Phys. Rev. A 36, 3193 (1987).
  • [14] K. R. Lykke, K. K. Murray, and W. C. Lineberger, Threshold photodetachment of H- Phys. Rev. A 43, 6104 (1991).
  • [15] M. K. Kristiansson, K. Chartkunchand, G. Eklund, O. M. Hole, E. K. Anderson, N. de Ruette, M. Kamińska, N. Punnakayathil, J. E. Navarro-Navarrete, St. Sigurdsson, J. Grumer, A. Simonsson, M. Björkhage, St. Rosén, P. Reinhed, M. Blom, A. Källberg, J. D. Alexander, H. Cederquist, H. Zettergren, H. T. Schmidt, and D. Hanstorp, High-precision electron affinity of oxygen Nat. Commun. 13, 5906 (2022).
  • [16] S. Azizi, Three aspects of photo-ionization in ultrashort pulses PhD thesis, Technische Universität Dresden 2023.
  • [17] The Dawson function F(z)𝐹𝑧F(z)italic_F ( italic_z ) [18] is related to the (complex) error function, F(z)=π2iez2erf(iz)𝐹𝑧𝜋2isuperscriptesuperscript𝑧2erfi𝑧F(z)=\frac{\sqrt{\pi}}{2\mathrm{i}}\mathrm{e}^{{-z^{2}}}\mbox{erf}(\mathrm{i}z)italic_F ( italic_z ) = divide start_ARG square-root start_ARG italic_π end_ARG end_ARG start_ARG 2 roman_i end_ARG roman_e start_POSTSUPERSCRIPT - italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT erf ( roman_i italic_z ).
  • [18] H. G. Dawson, On the numerical value of 0dxex2superscriptsubscript0differential-d𝑥superscriptesuperscript𝑥2\int_{0}^{\infty}\mathrm{d}x\,\mathrm{e}^{x^{2}}∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT roman_d italic_x roman_e start_POSTSUPERSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT Proc. London Math. Soc. XXIX, 519 (1897).