Linear dichroism in few-photon ionization of laser-dressed helium

Ionization of laser-dressed atomic helium is investigated with focus on photoelectron angular distributions stemming from two-color multi-photon excited states. The experiment combines extreme ultraviolet (XUV) with infrared (IR) radiation, while the relative polarization and the temporal delay between the pulses can be varied. By means of an XUV photon energy scan over several electronvolts, we get access to excited states in the dressed atom exhibiting various binding energies, angular momenta, and magnetic quantum numbers. Furthermore, varying the relative polarization is employed as a handle to switch on and off the population of certain states that are only accessible by two-photon excitation. In this way, photoemission can be suppressed for specific XUV photon energies. Additionally, we investigate the dependence of the photoelectron angular distributions on the IR laser intensity. At our higher IR intensities, we start leaving the simple multi-photon ionization regime. The interpretation of the experimental results is supported by numerically solving the time-dependent Schrödinger equation in a single-active-electron approximation.


Introduction
The fundamental interaction of photons and matter is omnipresent in nature. Scientists not only employ light to investigate matter to high precision, but also as a tool to modify it in a well-controlled manner. The continuous development and optimization of light sources and combining them in a variety of ways makes it possible to enter unprecedented realms in light-matter interaction. For example, intense infrared (IR) pulses focused into a gas target can trigger high-harmonic generation of extreme ultraviolet (XUV) photons [1,2]. The two intrinsically synchronized radiation fields can be used in a variety of experiments on atoms and molecules. The scheme of XUV excitation and subsequent IR ionization is employed in photoelectron measurements to determine angular-distribution parameters, relative population strengths, and delay-dependent ionization yields [3][4][5][6]. Complementary transient absorption measurements employed an energetically broad XUV spectrum in combination with the IR driving laser and revealed light-induced states (LIS) in helium [7]. These LIS a e-mail: severin.meister@mpi-hd.mpg.de (corresponding author) can be attributed to two-color multiphoton excitation, allowing transitions beyond the one-photon electricdipole selection rules. LIS were investigated in the context of transient absorption (see [8] for a review) with regard to their intensity dependence [9], quantum interference [10,11], and also for molecules [12]. Reduzzi et al. [13] brought up the aspect of polarization control, as the population of LIS depends on the relative orientation of the laser polarizations. We note that apart from the transient absorption context these states are dubbed dark states, as they are not accessible by singlephoton absorption. In the following we will refer to the basic helium eigenstates, but the notation obviously depends on the viewing angle.
XUV-IR photoelectron measurements employing synchrotron [14] or free-electron laser (FEL) [15] radiation in combination with an IR laser have also been reported. These experiments, however, focused on photoelectron angular distributions (PADs) stemming from single-photon excited states. Building upon our earlier work on photoelectron yields originating from multiphoton excited states [16], the present study comprises a detailed analysis by employing multidimensional data sets. Here, we focus on polarization-dependent or dichroic effects of PADs and their laserintensity dependence.
The measurement employs FEL-XUV radiation in combination with a synchronized IR laser. Both sources are linearly polarized with adjustable orientation relative to each other. For temporally separate (nonoverlapping) pulses, the preceding XUV pulse excites the helium atom to specific (1snp) 1 P states, while the delayed IR pulse ionizes the excited atom. For parallel orientation of polarizations and pulses in temporal overlap however, helium can be excited from its (1s 2 ) 1 S ground state to (1sns) 1 S and (1snd) 1 D states by simultaneous two-photon XUV-IR excitation. By absorbing additional IR photons, the excited atom can be ionized while details of the excited state can be inferred from the angular distribution and the kinetic energy of the emitted electrons. Measured PADs resemble hydrogen orbitals with varying angular momentum but conserved magnetic quantum number m = 0. For orthogonal orientation of the polarization axes, the PAD is altered fundamentally. The axial symmetry of the PAD along the polarization axis is broken, as the magnetic quantum number is m = ±1 with the IR defining the quantization axis. Furthermore, in this configuration excitation of 1 S states is completely suppressed, as predicted by selection rules for two-photon excitation [17]. Along with excitation, photoionization is suppressed at specific XUV energies, as the resonance enhancement is not available anymore.
We also report on a multiphoton feature just below the (1s2p) 1 P resonance involving excitation by one XUV and two IR photons. The analysis of the PADs allows to infer the excitation mechanism and an assignment to the (1s4f ) 1 F resonance.
Moreover, the IR intensity dependence of PADs is analyzed. Since the XUV pulse deposits most of the required excitation energy into the helium, IR ionization of excited states can be studied at lower laser intensities than in a single-color experiment. Our analysis suggests that the simple multi-photon description is valid except for high IR intensities, as our high-intensity data set shows. Finally, the interpretation of our experimental results is supported by numerical solutions of the time-dependent Schrödinger equation (TDSE) in the single-active electron (SAE) approximation.
This manuscript is organized as follows. Section 2 provides a summary of the experimental (Sect. 2.1) and theoretical (Sect. 2.2) methods. Our results are presented in Sect. 3, separated for ionization via the (1s3d) 1 D and (1s2s) 1 S excited states 3.1, the (1s4f ) 1 F excited state 3.2, and the IR intensity dependence of the PADs 3.3. A summary is given in Sect. 4.

Experimental Setup
The experiment was carried out with the reaction microscope (REMI) endstation [18,19] at the free-electron laser in Hamburg (FLASH2) [20,21] and a short-pulse IR laser based on optical parametric chirped-pulse amplification (OPCPA) [22]. The wavelength tunability of FLASH2 [23] was employed by scanning the XUV photon energy between 20.4 eV and 23.2 eV in steps of 0.1 eV with a bandwidth of 0.2 eV full-width at half maximum (FWHM). The FEL pulse duration was estimated to be about 40 fs FWHM in intensity and the FEL pulse energy was reduced to < 10 nJ, a level where two-XUV-photon absorption is suppressed. The IR laser was stretched to pulse durations of 90 fs (FWHM) to maintain, temporal overlap with the XUV pulses despite the timing jitter of several ten femtoseconds. The IR focal spot diameter of about 30 µm is significantly larger than the FEL focus (≈ 10 µm) and reached an intensity on the order of 10 12 W/cm 2 . The effective repetition rate in the measurement was 500 Hz.
The IR and FEL pulses enter the REMI collinearly while their foci overlap in the REMI center. Both the XUV and the IR radiation are linearly polarized, and the IR polarization can be rotated with respect to XUV polarization. A dilute supersonic gas jet of atomic helium crosses the laser foci, triggering single photoionization events. Charged fragments are separated by a homogeneous electric field and guided onto two timeand position-sensitive detectors. In order to achieve high momentum resolution, we selected a low electric field strength of 1.2 V/cm. A magnetic field, collinear with the electric field, forces charged particles on spiral trajectories ensuring a detection in the entire solid angle of 4π [24]. The chosen settings allow a 4π detection of electrons with kinetic energy up to 5 eV. Time-of-flight and position information is used to retrieve the particles' momentum vectors at the time of emission. Imposing momentum conservation on ions and electrons produced during the same laser pulse makes it possible to sort out fragments stemming from the same parent particle.

TDSE Calculation
The theoretical part of this study is based on numerically solving the TDSE in the SAE approximation. Although the model was described already in [16], we repeat it here for completeness and describe the extensions to implement the polarization. The ground state is effectively described as a (1s1s ) 1 S state, where the 1s is close to the He + orbital and the 1s is treated like a valence orbital, similar to the n orbitals in all the other 1sn Rydberg states. As always in theory, the binding energies of the n valence electrons are not exact (see Fig. 1 for a quantitative overview), especially since we used a single-electron model potential rather than fully accounting for the electron-electron interaction and correlation in an ab initio manner. While 1sn Rydberg states with orbital angular momenta ≥ 2 have very accurate binding energies, this is not quite the case for p electrons and particularly for s electrons, due to the small or missing centrifugal barrier. Since excita- Fig. 1 Comparison of helium excitation energies between NIST [25] and the model of this study. With the ground state set at 0.00 eV, this indicates that the excitation energy of the 2s state is underestimated, while all the other excitation energies are slightly overestimated tion energies are measured from the ground state, much of the remaining discrepancies are due to the binding energy of the 1s orbital.
Specifically, we used the one-electron potential, where r is the distance from the nucleus, to calculate the valence orbitals. The deviations in the excitation energies, compared to the recommended values from the NIST database [25], are less than 0.2 eV even in the worst-case scenario (see Fig. 1) and do not alter the essential conclusions presented below. In the following we will usually omit the inner 1s electron and the total spin state (it would be a singlet in the two-electron description) to simplify the notation. Doubly excited states are not accessible with the photon energies used in the present experiment, and spin-forbidden transitions (i.e., excitation of triplet states) are negligible. The laser intensity was chosen according to the available knowledge regarding the actual experimental conditions. The XUV pulse duration was taken as 40 fs (FWHM value of a peak intensity of 10 10 W/cm 2 with a Gaussian envelope) and the IR pulse duration as 80 fs (FWHM). While the XUV photon energy was varied over a range in steps of 0.05 eV, the central IR photon energy was held fixed at 1.55 eV (800 nm).
Due to the fact that the lasers are linearly polarized along different directions, the cylindrical symmetry of the problem that we could take advantage of in the previous work [16] is no longer valid. The complications, however, are manageable in a straightforward way by employing an updated version of the code described by Douguet et al. [26]. Specifically, we used the velocity gauge and a variable radial grid step starting at 0.1 a.u. near the nucleus and gradually increasing to 0.2 a.u. far away, where we only need to describe rel-

Fig. 2
Excitation and ionization scheme of the XUV energy scan in helium. For preceding XUV pulses there is only single XUV-photon absorption inducing p-state excitation (hatched orange region). For excitation into states with = 1 a combination of an XUV photon and IR photons is required. Two-color multi-photon excitation is only possible for temporally overlapping pulses atively slow electrons. The timestep was set constant at 0.05 a.u., partial waves up to max = 12 were coupled, and the numerical convergence of the predictions was checked by varying the parameters. Any apparent shortcomings of the theoretical results presented below are expected to be due to either the deficiency of the SAE model itself and/or the remaining uncertainties in the knowledge and, consequently, the description of the laser parameters.

Dichroism in Ionization via the 3d and 2s Excited States
A straightforward interpretation of the presented excitation and ionization mechanism is provided for XUV pulses preceding IR pulses. At matching XUV photon energies in the scanning range, the helium atom can be excited from the 1s ground state to np states. By absorbing additional photons of the moderately intense (≈ 4 × 10 12 W/cm 2 ) delayed IR pulse, the excited atom can be ionized. The energy of the measured photoelectron depends on the binding potential of the excited state and the number of absorbed IR photons. ization mechanism. The hatched regions in orange indicate the pathways that are accessible in the case of preceding XUV pulses. For example, the hatched region on the left shows ionization of helium via the 2p (21.2 eV) state. One XUV photon excites the atom to the 2p state and the additional absorption of at least three IR photons (1.55 eV each) from the delayed IR pulse results in ionization. The emitted photoelectron is detected and its measured kinetic energy can be reproduced by calculating the excess energy as follows: The atom can also absorb more photons than the minimum number required for ionization. The effect of above-threshold ionization (ATI) [27,28] manifests itself as regions of high yield separated by 1.55 eV in the electron energy.
The gray "continuum" region in the top of Fig. 2 illustrates the experimentally accessible part of the scheme. Respectively, Fig. 3 depicts the measured photoelectron yield as a function of the photoelectron kinetic energy and the scanned XUV photon energy (Eγ-plot) for temporally delayed IR pulses. Here, the same XUV energy range as in Fig. 2 is depicted, while Fig. 4 and Sect. 3.2 will focus on regions B and A (cf. Fig. 2), respectively. According to the mechanism described, Fig. 3 shows a high yield of photoelectrons at XUV photon energies of the 2p (21.2 eV) and the 3p (23.1 eV) states. ATI peaks are found along the electron-energy axis, separated by the energy of one IR photon.
The diagonal line, starting at around 21.8 eV in Fig. 3, can be attributed to the third harmonic of the FEL. At this fundamental frequency, the third harmonic has an energy of 65.4 eV, just sufficient to ionize the helium atom (I p = 24.6 eV) and leaving the He + ion in its first excited state (n = 2, 40.8 eV). The diagonal has a slope of three, as the photoelectron takes the excess energy.
For XUV and IR pulses in temporal overlap, the situation changes substantially, as various excitation pathways open up. In this case the helium atom is dressed in the IR radiation field during the XUV pulse and can absorb or emit IR photons, allowing access to transitions beyond the one-photon electric dipole selection rules. The simplest case of two-color two-photon excitation allows to couple the 1s ground state to ns and nd states. By absorbing additional IR photons, the excited atom can again get ionized, and the emitted photoelectron is measured.
In the case where both radiation fields are simultaneously present, all ionization pathways depicted in Fig. 2 are accessible. For excitation into the 3d state, the helium atom conjointly absorbs one XUV photon and one IR photon from the laser field. We also observe excitation into the 2s state, where conjointly to the XUV-photon absorption one IR photon is emitted to the laser field. As reported in reference [16], these two states dominate the light-induced features for the chosen IR and XUV parameters between the 2p and the 3p states. A detailed measurement in the XUV energy region B (cf. Fig. 2) is presented in Fig. 4. Subplot (a) is recorded for parallel polarization and shows two regions of high yield around 21.8 eV and 22.3 eV, which can be attributed to the 3d and the 2s states, respectively. The distribution of electron energies is broad due to the IR laser bandwidth and shortened lifetimes of the rapidly ionizing states. The broadness of the peaks along the XUV energy axis can be attributed to the XUV bandwidth and the AC Stark shift. The latter effect becomes Fig. 4 Eγ-plot of region B (cf. Fig. 2). Experiment: a parallel, b orthogonal polarization. Calculation: c parallel, d orthogonal polarization even more pronounced as the timing jitter between pulses leads to varying IR intensities when the XUV pulses are present. This effect was accounted for in the TDSE calculation of Fig. 4c by summing over three IR intensities (1, 2 and 3 × 10 12 W/cm 2 ). Here, the two dominant multiphoton excitation features of the experiment are well reproduced, even though they are located at slightly different XUV photon energies. As mentioned in Sect. 2.2, the binding and, consequently, the excitation energies of our model deviate from the ones recommended by NIST. In fact our model underestimates the excitation energy of the 2s state and overestimates that of the 3d state, which explains the remaining differences between calculation and experiment.
Two-color multi-photon excitation includes the additional degree of freedom to change the relative orientation of the polarizations. By exploiting this handle, specific excitation pathways, and with it photoionization, can be deliberately suppressed. The measurement depicted in Fig. 4b is recorded with the same parameters as in (a), except for the orthogonal orientation of the polarization axes. While the 3d feature remains, the signal of the 2s state vanishes. This finding is predicted by the "two-photon electric-dipole selection rules" [17] for photons of differing frequency. According to this rule, the transition ΔJ : 0 → 0 of the total electronic angular momentum is forbidden, while the magnetic quantum number m has to change by Δm = ±1. Since spin-flips during the excitation process are negligible, the ground-state spin (S = 0) is conserved and L = J.
The elongated region of high yield from the 3d state along the XUV photon-energy axis can be explained by the AC Stark shift. A pronounced shift of this state is described in Chen et al. [7] and Meister et al. [16]. In combination with the temporal jitter of the two radiation fields, it leads to the elongated feature. The calculated data for the orthogonal case are depicted in Fig. 4d. Similar to the experimental observation, the 2s state disappears while the 3d feature remains.
In addition to the kinetic energy of the photoelectrons, their angular distribution is measured. This allows us to deduce the quantum numbers of the excited state from which the electron was emitted. Figure 5 shows the measured three-dimensional photoelectron momentum distribution, which contains the PAD. Bins in momentum space are color-coded by the number of counts; each dot represents a photoelectron count. The two-dimensional plots on the walls depict the momentum distribution integrated along one axis.
Electrons in both subplots stem from the 3d resonance and were recorded at 21.8 eV XUV photon energy with their kinetic energy in the range from 0.07 eV to 0.16 eV. The weak contribution of the third harmonic is negligible and shows an isotropic PAD.
The PAD in Fig. 5a was measured for parallel polarization axes. In this configuration, each contributing photon changes the angular momentum by Δ = ±1 while the magnetic quantum number m = 0 is conserved from its ground-state value. The change of and m upon photon absorption is depicted in Fig. 6. For the measured "f -like" ( = 3) continuum distribution emission is found along the polarization axes and in two rings orthogonal to them. This suggests that the ionization proceeds via a d state. The PAD changes substantially in the case of orthogonal polarizations. The axial symmetry along the IR polarization is broken, and the distribution resembles that of a hydrogen orbital with = 3 and m = ±1. This can be explained by taking the IR polarization as the quantization axis, so the linear polarization of the XUV is expressed as a superposition of right-hand and left-hand circularly polarized radiation fields. In this way, it is straightforward to see that the magnetic quantum number changes according to Δm = ±1 upon absorption of an XUV photon. The absorption of IR photons, on the other hand, does not change m. This is illustrated in Fig. 6.

Ionization via the 4f Excited State
For high IR intensities (≈ 8 × 10 12 W/cm 2 ), we find a further peak at 20.8 eV, below the 2p resonance peak . A LIS for similar IR parameters in the same XUV energy range was observed by Chini et al. [10] in a transient-absorption (TA) measurement. It was tentatively assigned to excitation of np states by the absorption of one XUV photon and two IR photons. The present measurement however, suggests the feature to be a resonant enhancement of ionization via the 4f state as explained in the following. Figure 7 shows an Eγ-plot of region A (cf. Fig. 2) for parallel polarization axes. In subplot (a), the XUV pulse precedes the IR pulse, and thus ionization is restricted to resonant enhancement via p-state excitation. Within the depicted XUV energy range (scanned  Fig. 2. For pulses in temporal overlap, the dressing IR-field opens up an additional excitation pathway. A new peak appears in Fig. 7b at slightly lower XUV photon energies and electron kinetic energies (E XUV ≈ 20.8 eV, E kin ≈ 0.6 eV) compared to the first ATI peak from the 2p state.
The following energetic considerations rule out certain excitation pathways and lead to the final assignment of the involved state. As the peak is only present for overlapping pulses, two-color excitation has to be considered. For excitation by one XUV photon plus one IR photon, even the lowest s and d states are out of reach (20.8 ± 1.55 eV) to be connected to the peak. The 2s (20.6 eV) state is too close to be populated, while the next 3s (22.9 eV) and 3d (23.1 eV) states are already too high to be populated. By considering an excitation involving absorption of one XUV photon plus two IR photons, there are nf and np states available. The angular momenta of the three photons either add up, changing the ground-state angular momentum by Δ = +3, or two angular momenta cancel each other resulting in Δ = +1. The case of negative Δ is omitted, because the angular momentum of the helium ground state is already the lowest possible. The total available excitation energy is E EXC = E XUV (20.8 eV) + 2 · E IR (1.55 eV) = 23.9 eV.
Possible excited states are 4p and 4f , both with a binding energy of about 23.75 eV, and the 5p and 5f states, both with a binding energy of about 24.05 eV. In order to decide which states are more likely to be populated (either n = 4 or n = 5), we consider the AC Stark shift. To estimate the shift of these higher excited states, we take the shift between the 2p peak at 1.05 eV in Fig. 7b and the field-free value of Eq. 2 at 1.2 eV. By including the shift of about 0.15 eV, the binding energy of the 4p and 4f states is found at the matching excitation energy of 23.9 eV. This estimate is reasonable, as the highly excited states (n = 4, 5) experience a similar energetic shift with IR intensity as the continuum level. The possible resonances involved in the laser-dressing peak in Fig. 7b are narrowed down to the 4p and the 4f states.
The final identification is performed by analyzing the angular distribution of the photoelectrons stemming from the excited state. Figure 8 shows the measured and calculated photoelectron momentum distributions for parallel orientation of the polarization axes. The three-dimensional data are integrated along the mutual propagation direction y of the radiation fields and normalized to the total yield to provide a convenient comparison. The two plots in the upper row show the measured PAD of ionization through the 2p state and the nearby light-induced 4f state. In both cases emission is found along the polarization axis x and in three rings around the axes, appearing as three vertical lines. This distribution is "g-like" and corresponds to a continuum state with angular momentum = 4. Consequently, this PAD pinpoints the excited state to be the 4f state rather than the 4p state, supported by monitoring the excited-state population during and after the pulse, which is possible in the calculation.
The two plots in the bottom row of Fig. 8 were generated by the numerical calculation described in Sect. 2.2. The electron kinetic energy was cut off at 1.8 eV to avoid an overlay of higher ATIs in the projection. The IR intensity was set to 3 × 10 12 W/cm 2 and the XUV photon energy range was adapted to select the peak of interest. We introduced a Gaussian uncertainty in the electron energy of 0.2 eV to mimic the experimental conditions, including the IR bandwidth and the Stark broadening due to the temporal jitter between the pulses.
The 4f peak is also observed for orthogonal polarization axes. The measured Eγ-plot in Fig. 9 samples region A (cf. Fig. 2) in steps of 0.1 eV. Similar to Fig. 7, subplot (a) is recorded for preceding XUV pulses, while (b) is recorded in temporal overlap. The peak appears at the same electron energy and XUV photon energy as in the case of parallel polarization.
Again, by analyzing the PAD of the peak, the ionization-enhancing resonance can be identified. Subplots in the top row of Fig. 10 depict the measured photoelectron momentum distribution for orthogonal polarization. In this case the magnetic quantum number increases to one, following the same argumentation as given in Sect laser propagation axis y. The node count corresponds to an = 4 continuum state and, therefore, to an = 3 excited f -state the electron was emitted from. There is good agreement in Fig. 10 between the experimental data in the top row and the theoretical predictions in the bottom row. Characteristics like the number of nodes and their relative intensity are reproduced as well as the diameter corresponding to the photoelectron kinetic energy.

IR Intensity Dependence of PADs
In the present XUV-IR scheme, the main energy contribution to ionization is provided by one XUV photon. This allows ionization of helium from excited states where electrons are less strongly bound. Electrons can be lifted over the remaining potential barrier by absorbing just a few additional IR photons, in contrast to ionization by an IR laser only, where strong fields are required. In the extreme case of tunnel ionization, electrons are predominantly emitted along the polarization axis, giving rise to the classical picture of electrons being accelerated along the electric field lines. On the contrary, at lower laser intensities within the multiphoton regime, complex PADs emerge, corresponding to discrete number of absorbed photons and angular momenta.
The present measurement exhibits PADs that indicate the photoionization process leaving the simple multi-photon regime. The Keldysh parameter γ [29] can be used to roughly quantify the localization within the continuous transition from the tunneling regime (γ < 1) to the multi-photon regime (γ 1). The parameter is defined as Here I p is the ionization potential, ω the frequency, and E 0 the maximum field strength (all in atomic units).
In the present experiment, helium is ionized from the 2p state, yielding γ high = 1.9 in the high IR intensity case (≈ 8 × 10 12 W/cm 2 ) and γ low = 5.3 in the low IR intensity case (≈ 1 × 10 12 W/cm 2 ). While γ low is well within the multi-photon regime, γ high enters the transition region. Figure 11 shows calculated PADs for XUV and IR pulses in temporal overlap with parallel polarization axes. The angle Θ is defined as the angle between the polarization axis and the emission direction of the photoelectron. Photoelectrons in Fig. 11 are selected to cover the first ATI peak of ionization via the 2p state (XUV energy from 21.1 eV to 21.5 eV, electron energy from 0.6 eV to 1.4 eV). The data are normalized to the total integral yield. The number of maxima in Fig. 11 corresponds to the four absorbed photons (compare Fig. 2) and thus to the accumulated angular momentum. At low IR intensity, the PAD exhibits very pronounced minima and maxima, while at higher IR intensity the distribution changes distinctly. The minima at  Corresponding to the calculation, Fig. 12 shows the measured PAD of electrons emitted from the first ATI at E XUV ≈ 21.2 eV and E kin ≈ 1.2 eV. Similar to theory, we find the PAD to change with IR intensity. This is not obvious in first place as the number of absorbed photons stays the same. The two minima around cos(Θ) = ±0.8 smear out for increasing IR intensity, as predicted by the calculation. However, in contrast to theory, the center maximum increases in magnitude.
The solid lines in Fig. 12 depict fitted results to the data with their corresponding parameters listed in Table 1. The PAD can be described by the finite sum of even-rank Legendre polynomials P n scaled with β n Table 1 Beta parameters of the fit in Fig. 12 parameters [30,31]: Here σ 0 is the total emission cross section, and the angle Θ is defined as before. As explained in Ref. [5], the β 8 parameter corresponds to g-waves generated by fourphoton absorption from the = 0 ground state, which is dominating in the presented case. For high intensities β 8 decreases as well as lower order β-parameters (cf. Table 1). To explain the qualitative and quantitative observations at high IR intensities, mechanisms of higher order have to be considered. In the simplest case, the emitted electron can absorb one IR photon and emit one IR photon while conserving the photoelectron energy. The PAD, however, which reflects the angular momenta, is altered by this process. This is taken into account by extending the sum of equation 4 up to β 12 , which allows to describe a six-photon process. However, the products β 2n P 2n for different n are not disentangled, which makes a strict quantitative statement about the contribution of higher orders tentative.
Mayer et al. [5] investigated the laser-intensity dependence of PADs for similar parameters in an HHG-IR laser scheme. However, their findings are only comparable with ours to a limited extent, as the periodic HHG-based XUV spectrum couples to multiple excitation pathways.
An explanation for the slight discrepancy between theory and experiment (Figs. 11,12 ), but also inconsistencies in Table 1, can be found in additional resonances being involved in the ionization process. The energy difference between the 2p state and higher-lying np states (n = 6, 7, 8) is about twice the IR photon energy. Theory might populate or miss some of these states due to its narrow IR bandwidth in contrast to experiment. Slight differences in resonance energies (imperfect model energies) and intensities (experimental temporal jitter) between the calculation and experiment lead to different contributions from the Rydberg np states, thus resulting in a discrepancy between the measured and calculated PADs.

Summary
We investigated photoelectron emission via excited states in laser-dressed atomic helium. The experimental results are in good agreement with TDSE calculations based on an SAE model. The simultaneous interplay of XUV and IR radiation enabled populating 1 S, 1 D and 1 F states, which are not accessible by single-photon absorption from the ground state. Two-color multiphoton excitation has an additional degree of freedom in the relative orientation of the polarizations. We analyzed the dichroic effect of different relative orientations on the PAD of electrons stemming from the (1s3d) 1 D state. For parallel orientation we found a PAD resembling a hydrogen orbital with = 3 and m = 0, while for orthogonal orientation there is = 3 and m = ±1. Furthermore, by changing the relative polarization direction, the photoionization signal can be suppressed at specific XUV energies. This can be explained by the fact that for orthogonal polarization axes the transition from the 1 S ground state to an 1 S excited state is forbidden by the two-photon selection rules, in contrast to the case of parallel polarization.
Finally, we investigated the IR laser intensity dependence of PADs obtained from photoelectrons of fixed energy. The electrons were emitted from the 2p state by the absorption of three IR-photons. At high IR intensities, we found the PAD to blur and higher-order terms to come into play. This observation can be regarded as the onset of the process to transit from the multiphoton to the tunneling regime. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.