Strong-field photoelectron momentum imaging of OCS at finely resolved incident intensities

Photoelectron momentum distributions from strong-field ionization of carbonyl sulfide with 800 nm central-wavelength laser pulses at various peak intensities from $4.6$ to $13\times10^{13}$ W/cm$^2$ were recorded and analyzed regarding resonant Rydberg states and photoelectron orbital angular momentum. The evaluation of the differentials of the momentum distributions with respect to the peak intensity highly suppressed the impact of focal volume averaging and allowed for the unambiguous recognition of Freeman resonances. As a result, previously made assignments of photoelectron lines could be reassigned. An earlier reported empirical rule, which relates the initial state's orbital momentum and the minimum photon expense to ionize an ac Stark shifted atomic system to the observable dominant photoelectron orbital momentum, was confirmed for the molecular target.


I. INTRODUCTION
Studies in strong-field physics aim at the understanding and control of the electron wave packet emitted by an atomic or molecular target following the exposure to intense radiation. Vast theoretical and experimental efforts are focused on, e. g., the determination of intensitydependent ionization probabilities [1, 2] and cut-off energies for direct and rescattered electrons [2], the recognition of Freeman resonances [3][4][5], the observation of channelswitching and -closing effects [6,7], and the analysis of the photoelectron orbital angular momentum [5,7] in the multi-photon regime as well as the imaging of the initial state's electron density distribution [8][9][10]. In general for strong-field ionization, the electron wave packet's initial distribution in phase space -and thereby its subsequent dynamics in the field -is strongly shaped by the intensity of the electric field [11]. Within the framework of multi-photon ionization the intensity-dependent ac Stark shift alters the energies of initial and final target states as well as of any intermediate state that is resonantly passed through. In addition, the continuum the electron is born into is raised by the intensity-dependent ponderomotive energy. From the perspective of tunneling ionization the intensity dictates the shape of the potential barrier to be traversed in order to reach the continuum [12]. As a result, the outgoing electron wave packet is substantially dependent not only on the target system, but also on the intensity of the driving field.
Resolving the incident intensity in a strong-field experiment is as fundamental as knowing the incident photon energy in optical spectroscopy. Especially for the comparison of experimental data with predictions from strong-field ionization models intensity selectivity is highly beneficial, since intensity integration does not blur the results. How-ever, in experimental investigations employing typical laser-focus geometries the recorded target response results from intensities ranging from 0 to the peak intensity. As a consequence, the target information encoded in the experimental observables is often highly obscured by integration over all incident intensities. Various schemes were reported that allow the investigator to overcome this issue and enable access to the target response in an intensity-resolved fashion. An intensity-selective scanning method, which relies on the extraction of charged particles created within a restricted slice of the longitudinal intensity distribution -implemented through an aperture -in conjunction with an intensity deconvolution step in postprocessing was introduced [13] and later further refined [2]. For cw laser beams the acquisition of intensity-difference spectra was described [14]. The approach proposed is based on the evaluation of the derivative of the target response with respect to the peak intensity. Through analytical inversion by means of a power-series expansion in intensity, namely multi-photon expansion technique, the retrieval of distinct-intensity responses from a set of intensity-averaged measurements at different peak intensities was demonstrated [15]. While the intensity-selective scanning and the multi-photon expansion technique employ rather complex post-processing, they essentially eliminate the effect of focal intensity averaging depending on the peak-intensity step size of the data set. Although through the acquisition of intensity-difference spectra, using the most common focus geometry, the impact of intensity averaging is only highly suppressed yet not completely removed, it represents the by far simplest approach of the three.
In this manuscript, the advantages of intensitydifference spectra are studied for pulsed laser beams and utilized to investigate the intensity-dependent photoelectron momentum distributions from strong-field ionization of carbonyl sulfide (OCS) molecules in the intermediate regime between multi-photon and tunneling ionization.

II. INTENSITY-DIFFERENCE SIGNALS
The differential volume occupied by a distinct intensity in a cw laser beam can be readily described [14]. This focal intensity distribution is referred to as one implementation of a 2D configuration. In the following, the corresponding case of a pulsed laser beam will be investigated.
The intensity profile of an effective laser focus in a typical gas-phase experiment is assumed to have the following geometric properties: in space it is rotation-symmetric about the wave vector and its longitudinal distribution is given by a rectangular function with the width of the diameter of the skimmed gas-target beam D, relying on D to be smaller than twice the focus' Rayleigh length z R . The transverse spatial and the temporal profile follow normal distributions with standard deviations σ r and σ t . For a distinct time during the laser pulse's propagation through the focus volume the three-dimensional intensity distribution can be written as (1) Here, I 0 is the temporal peak intensity and σ z = cσ t is the spatial equivalent of the temporal standard deviation σ t . If D σ z , the full intensity spectrum, i. e., for all possible longitudinal positions of the pulse, will be quasiproportional to the intensity spectrum at a single position far away from the longitudinal borders given by the edges of the gas-target beam. Thus, the intensity distribution displayed above can be used to deduce the intensity spectrum of the complete space-and time-integrated focus interaction. In the following, a focus with such an intensity distribution is referred to as 3D configuration. After rearrangement of (1) to it is obvious that points of equal intensity are located on the surface of a spheroid. Accordingly, the volume function is with the differential volume −(∂V 3D /∂I) dI that is occupied by the iso-intensity shell dI around I following This provides a measure for the weight of any intensity between 0 and I 0 in the laser-target interaction. It complements the previous elaborations [14] by the additional consideration of the laser pulse's temporal intensity profile. Fig. 1 shows the resulting intensity spectra for two peak intensities I  difference spectrum. Evidently, in the difference spectrum the contributions from intensities < I (a) 0 are highly suppressed, albeit not completely eliminated. The pole point at I = 0 does not affect the determination of any intensity-dependent quantity, because the corresponding signal scales with ∆(−∂V 3D /∂I) · I n , n > 0. Here, n would be an integer for a pure single-channel multiphoton process but does rather represent an intensity-dependent quantity for a realistic strong-field interaction. The fullvolume and exemplary difference signals are shown in Fig. 2 for both, the 2D and the 3D configuration, as a function of intensity and intensity order n. Note that for a focal interaction with intensities normally distributed only along two dimensions (2D configuration), i. e., either a spatial laser sheet with a Gaussian temporal profile or a two-dimensional spatial profile with a flat-top temporal profile, the corresponding difference intensity spectrum carries non-zero contributions only for I (a) 0 Fig. 3 displays the mean intensities for various intensitydifference signals in dependence of the intensity order. In general, the use of a 2D focus configuration gives rise to mean intensities that are closer to the peak intensity compared to the 3D case. The black lines illustrate the change of the mean intensity with increasing intensity order for measurements at a single peak intensity. Especially for low intensity orders the corresponding mean intensities are far below the (upper) peak intensity. In contrast, in an intensity-difference setup the situation is much better, further improving with rising peak-intensity ratios Q = I 0 . The green graphs, Q = 0.99, provide insight into the achievable mean intensity for Q → 1: for the 2D configuration the mean intensity reflects, almost independently of the intensity order, the actual peak intensity. However, for the 3D case there is an appreciable deviation, especially for signals that arise at low intensity orders. Furthermore, as Q increases the fractional signal Normalized signals as a function of intensity and intensity order n for different volume differentials: a) full volume in 3D configuration, b) differential volume in 3D configuration, c) full volume in 2D configuration, d) differential volume in 2D configuration. These provide a measure for the relative contribution of a distinct intensity to an overall signal, that scales with intensity to the power of a fixed exponent. For the differential-volume signals a ratio of Q = I = 0.9 was used, which represents the lowest Q value employed in experiment. All signal maps depicted were row-wise normalized to their maxima. For the 2D configuration case the volume differential follows −∂V2D/∂I ∝ −1/I [14].
that is kept with respect to the full signal at peak intensity becomes smaller. This behavior is depicted in Fig. 4 and is identical for 2D and 3D configurations. In actual experiments one needs to find a compromise between intensity resolution and remaining signal.
In the following paragraph intensity-integrated observables and their derivatives with respect to peak intensity are elucidated. Any M -dimensional observable S(I 0 ), e. g., a two-dimensional projection of photoelectron momenta, that results from the full laser-target interaction at peak for the 3D (solid lines) and 2D configurations (dashed lines). The black lines correspond to signals at full volume integration, Q = 0, i. e., evaluated at a single peak intensity. Q values of 0.9 and 0.99 represent the lower and upper limits of peak-intensity ratios used in our experiments, vide infra.
intensity I 0 , obeys the relation with a normalized M -dimensional distribution s(I) that is created at the distinct intensity I with probability Ω(I) and the volume −(∂V /∂I) dI occupied by the iso-intensity shell dI. If −(∂V /∂I) dI is independent of I 0 (2D configuration), the evaluation of dS/ dI 0 enables the direct determination of Ω(I)s(I) by elimination of the integral sign in (5). This is not viable for the 3D configuration, since the volume differential is a function of I 0 . But also for this case dS/ dI 0 is proportional to an effective distributions, that arises from a restricted intensity spectrum as depicted in Fig. 1 and 2. Thus, the acquisition of dS/ dI 0 , or its discrete equivalent ∆S/∆I 0 , allows for the investigation of intensity-dependent changes in the signal distribution, while intensity-dependent occurrence probabilities, Ω(I), are not accessible.
Although a 2D focus configuration would enable the determination of intensity-dependent ionization rates as well as momentum distributions without contamination from intensities I < I (a) 0 , the experimental realization of such a setup suited for intensities appropriate for strongfield ionization studies proves to be challenging. However, the evaluation of ∆S/∆I 0 in a 3D configuration is easy to implement and provides the investigator with incidentintensity resolved signal distributionss. This enables analyses of the intensity-dependence of an M -dimensional signal, while the information about the integral signal remains obscured.

III. EXPERIMENTAL SETUP
The experimental setup was described elsewhere [16,17]. In brief, a cold molecular beam was created by supersonically expanding OCS seeded in 85 bar of helium through an Even-Lavie valve [18] operated at a repetition rate of 250 Hz. The skimmed molecular beam was intersected in the center of a velocity-map imaging spectrometer by a laser beam from a mode-locked Ti:sapphire system, which delivered Fourier-limited pulses at 800 nm central wavelength and 35 fs pulse duration (full width at half maximum) at a repetition rate of 1 kHz. Employing a 500 mm focal length biconvex lens resulted in a laser focus with a beam waist of 44 µm and a Rayleigh length of z R = 7.5 mm. The diameter of the molecular beam was found to be D = 2.2 mm and 2z R > D provides justification for treating the longitudinal intensity distribution of the effective focus as a rectangular function in the previous section. The pulse energy reaching the interaction region was defined using a half-wave plate and a thin-film polarizer. Ionization of OCS could be studied employing peak intensities of 4.6-13 × 10 13 W/cm 2 , respectively. Due to the chargedissipation limit of the in-vacuum detector components measurements at I 0 > 8 × 10 13 W/cm 2 were conducted at a reduced sample density -namely 1/3 of the density employed for measurements at I 0 ≤ 8×10 13 W/cm 2 . This was achieved by delaying the laser pulses with respect to the maximum of the molecular beam's temporal profile. Sample densities of 1 × 10 8 cm −3 and 1/3 × 10 8 cm −3 were utilized, respectively. Peak intensities as determined by a combination of beam-profiling, auto-correlation, and measurement of the pulse energy were found to be too low by 20 % compared to the intensity-dependent kinetic energy shift of distinct above-threshold ionization (ATI) peaks, see (8). In combination with the known energy spacing between adjacent ATI peaks, namely the photon energy hν, this approach allowed for an intrinsic calibration of both, the photoelectron momentum and the laser intensity. Two-dimensional projections of photoelectron momenta were recorded by mapping the electrons onto a position-sensitive detector, consisting of a microchannel plate (MCP), a phosphor screen and a high-frame-rate camera (Optronis CL600), by means of a velocity map imaging spectrometer (VMIS). This detection scheme allowed for counting and two-dimensional momentum-stamping of single electrons. The momentum resolution achieved -given by the voltage applied to the VMIS electrodes and the pixel size of the camera -was 5.6 × 10 −3 atomic units (electron velocity of 12 km/s). Momentum slices of the full three-dimensional momentum distributions S(I 0 ) were recovered through inverse Abel-transformations employing the "Onion-Bordas" onion-peeling algorithm [19,20] as implemented in the PyAbel software package [21].

A. Integral photoelectron yield
Extensive strong-field ionization studies [1] employing various laser peak intensities demonstrated that the asymptotic slope of an integral ionization signal Σ(I 0 ) with respect to ln I 0 , arising from the full intensity spectrum of the laser-target interaction at peak intensity I 0 , follows α is the instrument sensitivity, σ r the standard deviation of the transverse intensity distribution, D the length of the focus volume and the sample density. The only requirements for the validity of this equation are the predominance of single ionization and the cylindrical symmetry of the spatial focus geometry. Due to the onset of saturation, starting at intensity I sat , the slope of Σ(ln I 0 ) converges to a constant value for increasing peak intensity. Within the notation of this manuscript Σ(I 0 ) equals  6. Photoelectron momentum slices a) for the full-intensity distribution S at a peak intensity of 5.7 × 10 13 W/cm 2 and b) for the differential-intensity distribution S , which is dominated by intensities of 5.4-5.7 × 10 13 W/cm 2 . The black-line wedge sketches the polar angles used to evaluate the intensity-dependent radial distributions shown in Fig. 7 and 8; the black circles separate contributions from distinct photon numbers N and correspond to kinetic energies T (N ± 0.5,Ĩ0) as obtained from (8).
the momentum-space integral of the distribution S(I 0 ), which for the present case is just the total number of electrons detected. Fig. 5 shows the resulting values of Σ plotted versus I 0 on a semi-logarithmic scale. With increasing values of lnI 0 the slope of the total electron yield converges to a constant value. In the following analysis, this behavior is used as implicit evidence for a negligible contribution from multiple ionization. The asymptotic slope and the saturation onset were deduced by modeling a line through the highest peak-intensity points providing a saturation onset of I sat = 7.2 × 10 13 W/cm 2 . The saturation intensity is a function of both, the pulse duration and the orientation-dependent ionization probability [1]. Using the MCP's open area ratio of 0.64 as instrument sensitivity α, a transverse focal standard deviation of σ r = 22 µm (beam waist of 44 µm) and a molecular beam diameter of D = 2.2 mm resulted in a sample density of = 1 × 10 8 W/cm 3 . Since the intensity could be intrinsically calibrated from the ponderomotive shift of the ATI peaks (see next subsection), the deduced sample density is as accurate as ασ 2 r D.

B. Radial and angular distributions
Due to their rotational symmetry the three-dimensional momentum distributions S(I 0 ) can be completely described by two-dimensional slices of momenta parallel (p ) and perpendicular (p ⊥ ) to the polarization axis. In the following, the slope of S with respect to I 0 , is analyzed regarding its radial and angular distributions.
Here, the mean intensityĨ 0 corresponds to the arithmetic mean of every two peak intensities I that were used to measure S. Since S (Ĩ 0 ) arose from a more narrow distribution of incident intensities than S(I 0 ), it allows for the detailed investigation of intensity dependent effects.
Both S(I 0 ) and S (Ĩ 0 ) exhibit rich angular and radial structure. Fig. 6 displays a typical momentum slice for the two cases, here for intensities around 5.5 × 10 13 W/cm 2 . Both momentum distributions show broad radial patterns that repeat with integer multiples of the photon energy -referred to as ATI structure -, sharp radial features on top of the ATI structures -referred to as Freeman resonances [22] -and a carpet-like pattern perpendicular to the polarization axis [23].

Radial distributions
In Fig. 7 the peak-intensity dependent radial distributions for polar angles between −20 • and +20 • with respect to the polarization axis are illustrated for S(I 0 ). There are two kinds of features to be noted: (i ) broad peaks with a spacing of hν that experience a shift in kinetic energy proportional to the negative value of the peak intensity and (ii ) sharp peaks at constant kinetic energies. The radial distributions of S(I 0 ) are dominated by feature (ii ). Fig. 8 displays the equivalent radial distributions for S (Ĩ 0 ), which exhibit predominant contributions from feature (i ). Note that there is a change in sample density and intensity bin size at 8.3 × 10 13 W/cm 2 , c. f. Fig. 5.
In the following, feature (i ), which corresponds to nonresonant above-threshold ionization, is analyzed employ- ing the radial distributions of S (Ĩ 0 ). Dashed black lines in Fig. 8 correspond to final kinetic energies [24] T (N, with the number of absorbed photons N , the photon energy hν = 1.55 eV, the field-free ionization potential of OCS, I (0) p = 11.19 eV [25], and the ponderomotive potential U p =Ĩ 0 /4ω 2 (for atomic units). Clearly, (8) describes the kinetic energies of the radial local maxima in Fig. 8 quite well except for minor deviations towards lower kinetic energies at T < 3 eV, which are ascribed to the Coulomb interaction between photoelectron and cation. Thus, the total ionization energy appears to be dictated by the ponderomotive continuum shift. For the intensity range of the current study the relative Stark shift between the neutral and the cationic species is dominated by the permanent-dipole interaction [8] and, hence, equals zero in average over the isotropic distribution of molecular orientations.
Within the peak-intensity range investigated all radial distributions of S(I 0 ) carry series of sharp peaks, feature (ii ), that are also reproduced periodically with the photon energy hν. Close inspection of Fig. 8 reveals that the kinetic energies of those peak series are independent of the incident intensity, i. e., they appear as vertical lines imprinted onto the non-resonant ATI pattern, feature (i ). Based on this intensity-independence and the energy spacing between individual series, these peaks can unambiguously be assigned to Freeman resonances [22], i. e., photoelectron emission from resonance-enhanced multiphoton ionization through Rydberg states. Since those Rydberg states experience the same ac Stark shift as their corresponding continua, namely the ponderomotive energy, the resulting final kinetic energies T F are independent of the incident intensity [24,26]: with the number of photons k that ionize a resonant Rydberg state with principal quantum number n and the Rydberg constant R ∞ . The subtrahend reflects the energy difference between the Rydberg state and its asymptotic unbound state. Resonant transitions for n = 4, 5, 7 could be observed together with a combination band for n = 8-10, with k = 1-3. These assignments agree with previous single-photon absorption studies on OCS [27][28][29]. In the present study individual angular momentum states could not be resolved: at a kinetic energy of 1.2 eV, where the k = 1 Freeman manifold occurred, an energy resolution of 0.05 eV was achieved, which is roughly ten times larger than what would be required to distinguish between different l states.
Previous results on OCS for comparable peak intensities of 2.8-9.4 × 10 13 W/cm 2 , analyzed without the intensitydifference spectra, lacked information about the incident ponderomotive shift and, hence, assigned intensityindependent photoelectron lines to valence-valence transitions [4]. Three sharp peaks in the photoelectron kinetic energy spectra, which follow T 1,3 = 0.73 eV +khν, with k = (0, 1), and T 2 = 1.17 eV, were attributed to the valence transitions 1 ∆ ←X 1 Σ + and 1 Π ←X 1 Σ + . The corresponding excited valence states have field-free vertical ionization energies of roughly 5.8 eV ( 1 ∆) and 3.9 eV ( 1 Π) [28]. Both, the 1 ∆ and the 1 Π state, were considered to experience ac Stark shifts equal to the ponderomotive potential, resulting in intensity-independent kinetic energy releases [4]. Since U p is the kinetic energy imparted to an unbound electron in an oscillating electric field, it represents the asymptotic ac Stark shift of an electronic state at the onset of the continuum. However, treating the 1 ∆ and 1 Π states as quasi-unbound states is in conflict with their large binding energies. Thus, in our present investigation the very same three peaks are, instead, assigned to ionization of resonant Rydberg states leading to final kinetic energies of T 1,3 = T F (k = (1, 2), n = 4) and T 2 = T F (k = 1, n = 7).
While S (Ĩ 0 ) mainly carries contributions from nonresonant above-threshold ionization, S(I 0 ) is dominated by resonance-enhanced multi-photon ionization. Through averaging over all intensities present in the focal volume the ponderomotively shifted non-resonant contributions get washed out in the final momentum distribution S(I 0 ), while resonant contributions with constant final momentum, e. g., Freeman resonances, gain intensity [30]. By inspection of S(I 0 ) only, the contributions from resonant ionization pathways are highly overestimated and final kinetic energies cannot be related to distinct numbers of absorbed photons, because the ponderomotive continuum shift U p cannot be resolved. In the momentum distributions S(I 0 ) the closing of ionization channels with increasing peak intensities is highly obscured, because contributions from all intensities between 0 and I 0 are contained. Thus, even if a channel is already closed at the peak intensity, S(I 0 ) will embody all contributions from I < I 0 with even larger weight, i. e., larger differential volume, at which the channel is still open. In principle, S(I 0 ) contains just more and more open ionization channels with increasing peak intensity. Nevertheless, experimental evidence for channel-closing was reported for xenon at fully focal-volume averaged momentum distributions [5,6]. In the present study, three channel-closings could be clearly observed, visible in Fig. 8  (P5(cos θ)) 2 (P6(cos θ)) 2 (P5(cos θ)) 2 FIG. 9. Angular distributions of the lowest order nonvanishing ATI channels for the minimum numbers of photons of 10, 11, 12 required to ionize the OCS molecule. The three distributions were retrieved within the intensity intervals [4.6, 7.2[, [7.2, 9.8[, [9.8, 12.4[×10 13 W/cm 2 , which represent the channel-closing intensities obtained by setting T (N,Ĩ0) = 0 in (8). Experimental data are shown by circles, calculated traces of the squares of Legendre polynomials P l (cos θ) of order l with the polar angle θ to the laser polarization axis are shown by dashed lines. The distributions are plotted with an arbitrary vertical offset for clarity.

Angular distributions
In the following, the non-resonant ATI channels for ionization with a minimum number of photons are analyzed with respect to their photoelectron angular distributions. These channels give rise to the lowest radial layers of the carpet-like pattern [23] visible in Fig. 6a. The minimum number of photons needed to ionize the target is a staircase function proceeding upward with increasing intensity. Accordingly, S (Ĩ 0 ) was sectioned by the channel-closing intensities, which were obtained by evaluating the roots of (8). Within the intensity range investigated the closing of the 9, 10, 11, 12 photon channels were encountered at intensities of 4.6, 7.2, 9.8, 12.4 × 10 13 W/cm 2 , respectively. In Fig. 9, circles depict the resulting angular distributions for ionization by 10, 11, 12 photons, each representing the lowest order non-vanishing ATI channel in its corresponding intensity regime. The angular distributions associated with the absorption of 10, 11, 12 photons exhibit 5, 6, 5 local minima, respectively.
As shown by a theoretical study of the strong-field ionization of atomic hydrogen [31] and by experimental investigations of xenon [5], the dominant orbital angular momentum quantum number l of the emitted electrons can be deduced from fitting the angular distributions by squared Legendre polynomials P l (cos θ) 2 of order l. According to the number of local minima the absorption of N = 10, 11, 12 photons gives rise to dominant orbital angular momentum quantum numbers of l = 5, 6, 5. Best matches of P l (cos θ) 2 to the experimental angular distributions are shown as dashed lines in Fig. 9. An empirical rule was given that predicts the dominant photoelectron orbital angular momentum quantum number l based only on the minimum number of photons to reach the continuum and the orbital quantum number l 0 of the initial bound electronic state [32]. For intensities implying minimum numbers of 10, 11, 12 photons needed to ionize argon, ground state configuration [Ne]3s 2 3p 6 (l 0 = 1), dominant orbital-angular-momentum quantum numbers of l = 5, 6, 5 were found. These values for l agree with those obtained for OCS, ground state configuration . . . 9σ 2 3π 4 (l 0 = 1), providing clear evidence of the empirical rule and its extension to molecules.
The time-dependent Schrödinger equation for the strong-field ionization of H and Ar pictured the trajectory of the outgoing electron as a quiver motion along the laser-polarization axis superimposed to a drift motion following a Kepler hyperbola [33]. Identifying the hyperbola's pericenter as the electron's quiver amplitude in the field provided an analytical expression for an effective final orbital momentum that -for single ionization -only depends on the electron's final linear momentum and its quiver amplitude in the field. For the ionization of OCS through 10, 11, 12 photons effective orbital quantum numbersl = 5.5, 6.5, 7.1 were obtained accordingly [33]. Although this model considers neither the initial state's orbital momentum nor a parity selection rule, the values forl reflect our experimental observations remarkably well. This further supports the statement that the momentum distribution of an electron wave packet from non-resonant strong-field ionization is dictated by the electric field and its Coulomb interaction with the cation. The initial bound electronic state imprints onto the freeelectron wave packet through its vertical ionization energy, its orbital angular momentum, and the symmetry of its electron density distribution. While the impact of the former two properties was demonstrated within the scope of this manuscript, the investigation of the latter would require referencing to the molecular frame.

V. CONCLUSIONS
Three-dimensional photoelectron momentum distributions from strong-field ionization of OCS at differential peak intensities were recorded and analyzed in matters of resonant Rydberg states and photoelectron angular momentum. The evaluation of the derivatives of the momentum distributions with respect to the peak intensity allowed for the unambiguous recognition of Freeman resonances. As a result, sharp photoelectron lines, previously ascribed to valence-valence transitions [4], were reassigned to progressions of Freeman resonances corresponding to Rydberg states with principal quantum numbers n = 4, 5,7,8,9,10. Furthermore, the empirical rule that relates the initial state's angular momentum and the minimum photon expense to ionize an ac Stark shifted atomic system to the observable dominant photoelectron momentum [32] could be confirmed for a molecular target. In addition, the clear closing of ionization channels with increasing peak intensity is shown for the OCS molecule, which could eventually enable the probing of a distinct strong-field ionization channel in a coherent-control experimental setup. It is demonstrated that in order to gain insight into strong-field ionization processes on a quantum-mechanical level it is essential to relate distinct incident intensities to the resulting response of the system under investigation. Especially for the exploration of the transition regime between multi-photon and tunneling ionization a well-defined and narrow intensity distribution is crucial. Studies employing a clear relation between incident intensity and target response could finally shed light onto the importance of resonant states and, thus, show the shortcomings of present molecular strong-field ionization theories to describe the outgoing electron wave packet.