We have recorded the diffraction patterns from individual xenon clusters irradiated with intense extreme ultraviolet pulses to investigate the influence of light-induced electronic changes on the scattering response. The clusters were irradiated with short wavelength pulses in the wavelength regime of different 4d inner-shell resonances of neutral and ionic xenon, resulting in distinctly different optical properties from areas in the clusters with lower or higher charge states. The data show the emergence of a transient structure with a spatial extension of tens of nanometers within the otherwise homogeneous sample. Simulations indicate that ionization and nanoplasma formation result in a light-induced outer shell in the cluster with a strongly altered refractive index. The presented resonant scattering approach enables imaging of ultrafast electron dynamics on their natural timescale.

Intense femtosecond short-wavelength pulses from free-electron lasers (FELs) open new avenues to investigate transient states and ultrafast processes with unprecedented spatial and temporal resolution.1–4 One prominent example is ultrafast x-ray diffraction methods such as femtosecond Coherent Diffraction Imaging (CDI), which have enabled the structure determination of individual nonperiodic nanoscale objects.5 The elastically scattered photons of a single-shot exposure form an interference pattern containing a snapshot of the object before it is quickly destroyed due to the large amount of deposited energy.6,7 The encoded structural information can be retrieved via phase retrieval methods8 or forward simulations,9 which allowed for the structural characterization of such fragile objects as single viruses,10,11 aerosols,12,13 atomic clusters,9,14,15 and even superfluid helium nanodroplets containing quantum vortices.16,17 Via pump-probe techniques, laser-induced processes in individual nanoparticles can also be studied in a time-resolved manner with unprecedented spatiotemporal resolution.18–22 Such studies are pivotal for understanding and mitigating the damage dynamics from ionization, plasma formation, and particle explosion, which limit the achievable resolution in CDI experiments.7,23

While it may be possible to outrun the structural damage induced in the nanoscale targets,6 the ultrafast changes in the electronic structure due to excitation, ionization, and plasma formation occur on a faster, few-femtosecond or subfemtosecond timescale, unseparable from the interaction with the intense pulse.24–27 The sensitivity of the diffraction process to the particle's electronic structure, on the other hand, holds tremendous opportunities to trace electronic structure changes with high spatial resolution in a time-resolved manner with ultrafast diffraction methods.28 In particular, near absorption resonances, the x-ray scattering cross sections depend sensitively on the energy of the incoming photons and the electronic structure of the target.27–30 

In this work, we demonstrate that ultrafast resonant scattering can be used to visualize the spatial distribution of transient charge states in an evolving nanoplasma. As samples, we use submicrometer-sized xenon clusters that are transformed to a highly excited nanoplasma during irradiation and imaged with the same intense femtosecond FEL pulse. On the timescale of the pulse, the position of the clusters is frozen in space, and the measured ion kinetic energies show that ionic motion in the generated nanoplasma can be neglected. Nevertheless, in the radial intensity profiles of single-shot single-cluster scattering patterns, we find an intensity dependent lobe structure corresponding to the appearance of an additional characteristic length in the otherwise homogeneous particle. At the FEL wavelength used in our study, neutral xenon and also its low charge state ions (Xe≤4+) are strongly absorbing, while higher charge states are almost transparent. The choice of a resonant wavelength allows us to discriminate between areas of different charge states. We carried out a one-dimensional Monte Carlo simulation of the photoionization process, suggesting the formation of a highly charged outer shell in the evolving nanoplasma with strongly altered optical properties. By grouping and averaging the experimental patterns obtained at similar FEL intensities, we suppress individual cluster effects, e.g. from a rough surface, while enhancing the relevant dynamic signature. The radial profiles of the grouped patterns are fitted with Mie-calculations for a concentric core–shell model to extract tendencies of the evolution as a function of illumination intensity. The fitting yields a sequence of core–shell structures with strongly altered refractive indices and increasing shell thickness. The experiments demonstrate the possibility to extract spatial information on transient plasma formation in isolated nanoparticles with resonant ultrafast x-ray scattering. The method provides the potential for imaging ultrafast excitation, ionization, and charge transfer dynamics in complex samples with femtosecond time and nanometer spatial resolution.

The experiments were performed at beamline BL2 at the soft x-ray free-electron laser FLASH.31,32 Intense extreme ultraviolet (XUV) pulses at 91 eV photon energy were produced with an electron bunch charge of 0.5 nC, yielding a pulse energy of 150 μJ as measured with the gas monitor detector.32 An estimate for the pulse duration of 100 fs is derived for these parameters from measurements33 carried out with the electron bunch length diagnostics LOLA.34,35 Considering a beamline transmission of 64% and a focal spot size of 20 μm (FWHM) at the beamline BL2,32 a power density of up to 3×1014 W/cm2 was reached. The pulses intersected a highly diluted jet of sub-micron xenon clusters.36 An adjustable piezo-skimmer slit ensured that only one single cluster was present in the focus volume per FEL shot.15 The scattering patterns were measured with a previously described14,28 large area scattering detector consisting of an MCP-phosphor stack with a center hole and an out-of-vacuum CCD camera. In addition to the diffraction images, coincident single-shot ion time-of-flight (tof) spectra were recorded.37,38 In the polarization plane of the FEL, the setup allowed us to measure diffraction patterns at scattering angles between 3° and 30°. In the perpendicular direction, the detection angle was limited to 10° due to a shadow of the spectrometer's electrodes. Prior to further analysis, the measured scattering intensities were corrected for a nonlinear detector response9 by taking each pixel's intensity to the power of 2.5 and for the flat detector geometry28 by multiplying with a factor of cos(θ)3.

Two representative examples of the single-shot single-particle diffraction patterns are shown in Figs. 1(a) and 1(b). The difference in brightness results from different irradiation intensities of the FEL due to the varying positions of the xenon clusters within the focal volume.38 As known from previous work, the basic structure of the diffraction images, with concentric but intermittent rings, indicates nearly spherical shapes with rough surfaces resulting from the coagulation-dominated growth process.36 The size of each single cluster (average radius) could be determined by comparing the spacing of the extrema in the diffraction patterns with Mie calculations.36 A total of 32 patterns within the size regime of R=(400±50) nm were selected (the raw data of all patterns were uploaded to the CXIDB39,40). Then, the patterns were radially averaged and plotted on a logarithmic scale vs scattering angle in Fig. 1(c).

FIG. 1.

Isolated xenon clusters were irradiated with intense XUV pulses (91 eV photon energy and 3×1014 W/cm2 peak intensity in the center of the focal spot). A total of 32 events with single clusters of (400±50) nm radius were selected for analysis by the characteristic spacing of the diffraction rings. (a) and (b) Representative diffraction images (second brightest and darkest image of 32 events). (c) Radial profiles of the 32 single-shot images (corrected for the flat detector and nonlinear response, averaged over the scattering angle ϕ; see text). The color coding indicates the binning of events with similar intensities (the least intense category only contains a single pattern). (d) Radial profiles of averaged patterns from bins A to D. For increasing scattered intensity, an upward shift of the profiles (linear response) and an additional modulation of the profiles (corresponding to the ionization and plasma formation) can be observed.

FIG. 1.

Isolated xenon clusters were irradiated with intense XUV pulses (91 eV photon energy and 3×1014 W/cm2 peak intensity in the center of the focal spot). A total of 32 events with single clusters of (400±50) nm radius were selected for analysis by the characteristic spacing of the diffraction rings. (a) and (b) Representative diffraction images (second brightest and darkest image of 32 events). (c) Radial profiles of the 32 single-shot images (corrected for the flat detector and nonlinear response, averaged over the scattering angle ϕ; see text). The color coding indicates the binning of events with similar intensities (the least intense category only contains a single pattern). (d) Radial profiles of averaged patterns from bins A to D. For increasing scattered intensity, an upward shift of the profiles (linear response) and an additional modulation of the profiles (corresponding to the ionization and plasma formation) can be observed.

Close modal

A number of observations can already be made when following the profiles' evolution in Fig. 1(c) from the lowest to the highest intensity: first, a high-frequency modulation can be observed in all profiles, originating from the ring structures in the patterns [see Figs. 1(a) and 1(b)] and reflecting the cluster size information. In the case of the least intense profile (yellow curve), the envelope agrees rather well with the expected curve for a homogeneous spherical xenon cluster, dropping linearly on a logarithmic scale. In the absence of light induced changes in the particle, all other profiles from the clusters irradiated with higher FEL intensity would follow a similar curve, just with a proportionally higher scattering signal. Instead, the envelopes of the more intense profiles show an additional large-scale structure.

In order to analyze only the intensity dependent signature in the patterns and to reduce effects from irregular shapes and slightly different sizes, the events were sorted for increasing detector brightness and binned into four different categories A to D [the bins are indicated by the color coding in Fig. 1(c)]. The diffraction patterns within each of the four bins A to D were averaged, and radial profiles were extracted from the averaged patterns. These averaged profiles are plotted in Fig. 1(d).

When following the averaged profiles of categories A to D, the evolution of the profile envelopes can be seen even more clearly. A lobe structure appears with a minimum roughly at a scattering angle of 15° that becomes more pronounced with increasing irradiation intensity. In general, a modulation in any diffracted intensity distribution corresponds to a characteristic length scale in the scattering object. We thus conclude that the evolving superstructure corresponds to the development of an additional scattering structure in the sample, an area with different optical properties. It is notable that the modulation structure appears to be a general feature because it survives the averaging over many single-cluster patters, which themselves incorporate the average scattering signal over the FEL pulse duration. This raises the question of the origin of the transient refractive structure with a spatial extension on the order of a few tens of nanometers, as estimated from an Airy pattern with a minimum at 15°.

The results from theoretical41–53 and experimental14,20,37,38,54–67 studies on clusters in intense XUV pulses provide basic knowledge about the interaction between an intense short wavelength pulse and a rare gas cluster. From this body of work, we can exclude ionic motion as the origin of the observed modulation feature. A general picture divides the dynamics in three phases.42,48,54,68 In the first phase, photoionization and Auger decay lead to the emission of electrons with residual kinetic energy from the cluster and therefore result in a charge-up of the cluster as a whole.48,69 As soon as the Coulomb attraction between the positively charged cluster and photoelectrons (or Auger electrons) exceeds their kinetic energy, the second phase starts, and subsequently, released electrons become trapped in the cluster potential.41,56 A nanoplasma is formed in which further electrons are released from the single atoms or ions, but they reside within the cluster. Within the following pico- to nanoseconds, in the third phase of the cluster dynamics, the electrons transfer their kinetic energy to the ions, expelling one surface layer after the other in a hydrodynamic expansion.20,37,38,48 Also, the net charge on the cluster leads to ion repulsion of the unshielded surface, referred to as Coulomb explosion.45,57,60,70 A theoretical study48 modeling argon clusters irradiated with 90 eV radiation predicts a motion of the outermost cluster shell of about 1.5 Å within the first 100 fs of the interaction for conditions comparable to the current experiment (the argon ions being a factor of 3 lighter and the intensity being a factor of 6 smaller).48 Furthermore, we can get a first-order estimate for the maximum motion of ions from the ion time-of-flight spectra measured in our experiment. An estimate for the acceleration at the surface can be derived from the final kinetic energies of the cluster ions of up to 600 eV per charge.38 Assuming a 400 nm sphere that accelerates a Xe5+ ion to its final kinetic energy of 3000 eV yields an effective charge of that sphere of 105e. Such a sphere drives Xe5+, starting at the sphere's surface from rest, less than 1 Å within the first 100 fs.

Having excluded ionic motion, we attribute the observed modulations to light-induced electronic structure changes (i) resulting from the ionization and plasma formation within the FEL-irradiated particle and (ii) visualized by the use of a wavelength of the FEL resonant with neutral xenon and low charge states while being nonresonant with higher charge states. In Sec. III B, based on a first-order model of the cluster ionization process, we develop a physical picture of the plasma formation and discuss the generation of an outer shell in the cluster with strongly altered optical properties. Subsequently, using a classical concentric core–shell Mie model, we fit the modulation features observed in Fig. 1 in order to extract tendencies for the evolution of the shell’s thickness and optical constants with increasing irradiation intensity (Sec. III C).

The appearance of the dynamic scattering features can be connected to the peculiar electronic properties of xenon atoms and ions in the vicinity of the photon energy of 91 eV. Absorption cross sections for gas-phase xenon atoms and atomic ions have been measured71–76 and are summarized in Fig. 2. A clear step from high to low absorption between Xe4+ and Xe5+ can be observed, with extremely high values for the charge state 4+, which exhibits a large ionic resonance at 91 eV. Correspondingly, the penetration depth, i.e., the depth into the material for the intensity to decay to 1/e, increases from about 30 to 300 nm.

FIG. 2.

Absorption of neutral xenon atoms and atomic ions at 91 eV. Total absorption cross sections σabs in Mbarn of neutral Xe,71 Xe+,72Xe2+,73Xe3+,74Xe4+,75 and Xe5,6,7+76 (colored points). Note that the value of 2 Mbarn for 5 to 7+ constitutes an upper bound. The corresponding penetration depth in nm (black crosses) is calculated using labs=1na·σabs, with na being the atomic density of solid xenon.77,78

FIG. 2.

Absorption of neutral xenon atoms and atomic ions at 91 eV. Total absorption cross sections σabs in Mbarn of neutral Xe,71 Xe+,72Xe2+,73Xe3+,74Xe4+,75 and Xe5,6,7+76 (colored points). Note that the value of 2 Mbarn for 5 to 7+ constitutes an upper bound. The corresponding penetration depth in nm (black crosses) is calculated using labs=1na·σabs, with na being the atomic density of solid xenon.77,78

Close modal

To gain a first-order model of the radial distributions of different charge states in an irradiated cluster, we investigate the photoionization of a one-dimensional chain of atoms using an atomistic Monte Carlo approach. The results are presented in Fig. 3. For the simulation, each photon is “propagated” along the chain, starting at x =0 from the first atom, by probing at every atomic position whether an absorption is taking place. If the random drawing dictates an absorption process, the photon is annihilated, and the charge state of the corresponding atom is increased by 1 (or 2 in the case of an Auger process), followed by the start of the propagation of the next photon, again at x =0. According to the number density of solid xenon77 of na=1.67×1028 particles per m3, we consider 824 atoms for a chain length of 400 nm. The charge state dependent absorption probabilities for each atom are derived from the absorption cross sections of different xenon charge states.71–76 Only linear photoabsorption and subsequent Auger processes are taken into account, while nonlinear effects, light scattering, and plasma processes such as collisional ionization are neglected. The derivation of the absorption probabilities and the number of photons impinging on the one-dimensional chain of atoms, as well as a benchmarking of the model by simulating the penetration depth for individual charge states, are given in the supplementary material.

FIG. 3.

(a) Simulation of the distributions of the relative charge state abundances ρq(x) for a one-dimensional chain of 824 atoms, i.e., 400 nm length. 870 photons (corresponding to 1014W/cm2) fall on the geometric cross section of one xenon atom and are propagated along the chain. Absorption cross sections of atomic xenon and its ions from Fig. 2 are used for calculating absorption probabilities. (b) The average charge state q¯(x) drops from around 6+ to neutral within 80 nm. The relative charge state abundances further allow us to determine an effective absorption index β(x), revealing a transition within 50 nm by an order of magnitude.

FIG. 3.

(a) Simulation of the distributions of the relative charge state abundances ρq(x) for a one-dimensional chain of 824 atoms, i.e., 400 nm length. 870 photons (corresponding to 1014W/cm2) fall on the geometric cross section of one xenon atom and are propagated along the chain. Absorption cross sections of atomic xenon and its ions from Fig. 2 are used for calculating absorption probabilities. (b) The average charge state q¯(x) drops from around 6+ to neutral within 80 nm. The relative charge state abundances further allow us to determine an effective absorption index β(x), revealing a transition within 50 nm by an order of magnitude.

Close modal

In Fig. 3(a), the simulated distributions of the relative charge state abundances ρq(x) with q =0 to 8+ are presented for an irradiation intensity of 1014 W/cm2, corresponding to 870 photons falling on the respective area of a single xenon atom (calculated via π·rXe2 with the atomic radius of xenon rXe being determined from the solid density value77 given above). From the relative charge state abundances, the distribution of the average charge state q¯(x) and an effective imaginary part of the refractive index β(x) along the chain are calculated,78 as given in Fig. 3(b). Note that the complex refractive index n=1δ+iβ is a dimensionless quantity with the absorption index β and the so-called refractive index decrement δ, which is related to the phase shift of light traveling through matter. The effective absorption index along the chain is derived as

using the atomic/ionic cross sections given in Fig. 2 and the wavelength λ. Both curves given in Fig. 3(b) indicate that after irradiation, an outer shell exists up to a propagation depth of about 120 nm. In this depth, the average charge state q¯(r) drastically drops in a transition region of about 80 nm thickness from about six to zero, while the effective absorption index β(r) jumps from 0.004 to 0.05, revealing an even more pronounced kink.

Our basic one-dimensional Monte Carlo simulation therefore indicates the formation of an outer shell in the cluster with only highly charged ions, which is rather transparent as compared to the opaque core. As mentioned above, the simulation omits nonlinear effects, the light scattering process itself, and all plasma related processes such as collisional ionization. While it yields simulated values of the imaginary part of the refractive index β, it does not allow for a prediction of possible changes in the refraction (δ) of the shell. Nevertheless, it provides a first-order explanation for the dynamic diffraction feature arising from the resonant interaction of 91 eV radiation with the xenon clusters.

The ionization model described above indicates the formation of an outer part in the cluster with altered optical properties, but it is important to note that the geometrical nanoplasma structure may considerably deviate from a concentric core–shell system. Instead, the considerations of the ionization process suggest a distribution that is asymmetric in the direction of the incident light, with a transparent part at the irradiated side of the cluster, while the back of the cluster remains neutral. However, assuming a concentric core–shell offers the advantage that the patterns can be further analyzed using classical core–shell Mie theory,79,80 i.e., the analytic solution of the Maxwell equations for the case of a concentric core–shell system. The input values and parameters of the Mie simulation—the size of particle and shell, the complex refractive index in both areas, and the wavelength and intensity of the incoming light—provide important handles to capture main tendencies of the nanoplasma formation. While the influence of the expected asymmetry needs to be tested with advanced theoretical models,70 we note that the contribution from the back of the cluster to the diffraction should be also small in a concentric Mie-model with a highly absorbing core, shadowing the back of the cluster from irradiation.

For the Mie analysis, not the averaged patterns A to D obtained at similar irradiation intensities [radial profiles are shown in Fig. 1(d)] are considered but their difference signal. This approach is conceptually similar to resonant x-ray imaging, e.g., of the ultrafast switching of magnetic domains,81 where the difference signal between patterns obtained at different helicities of circularly polarized light is analyzed, or of buried structures,30 where diffraction patterns just above and below an absorption edge are subtracted from each other to enhance the difference signal in the location of a certain element.

Using a difference signal approach in combination with Mie fitting is based on the following considerations: the measured patterns (and also averaged data) do not correspond to a single, intermediate plasma state that can be approximated by a single nanoplasma structure such as a single core–shell system. Always, the onset of the FEL pulse intercepts a cluster that is neutral and unchanged, and the last photons of the pulse interact with an evolved nanoplasma. By analyzing the difference spectra with Mie core–shell fits, we discretize this evolution and link the intensity-resolved information in different profiles with the idea of a common plasma dynamics (see also the mathematical derivation in the supplementary material).

This includes the simplifying assumption that the patterns belonging to groups A to D, obtained at different irradiation intensities, all result in principle from the same continuous evolution of nanoparticle ionization and plasma formation but up to different stages. In other words, profile A corresponds to an only weakly irradiated cluster, profile B contains the response of this initial phase and additionally the response of a more advanced nanoplasma, and so on. This perspective is equivalent to the statement that the number of impinged photons is the decisive factor for the plasma state reached. This assumption allows us to replace the variable of time by the variable of energy, but it can only be a rough approximation. Implicitly, this means that all nonlinear processes such as multiphoton absorption are neglected (for the instantaneous absorption of two photons, the intensity, i.e., the number of photons per time interval, is relevant, not only the number of photons in total).

Now, instead of a gradually evolving and changing system, we may approximate the nanoplasma formation process by a sequence of a few discrete steps. The profiles of the difference signals DC, CB, and so on are fitted with a concentric core–shell Mie model (note that profile A can be seen as the difference between A and 0). In Fig. 4(a), the radial profiles from the difference signal of the averaged patterns A − 0, BA, CB, and DC are given. The difference profiles reveal even more distinct features in the superstructure: with increasing irradiation intensity, a broad lobe appears, which becomes more and more pronounced, narrows, and shifts toward higher scattering angles.

FIG. 4.

(a) Difference profiles from the averaged profiles shown in Fig. 1(d). For better visibility, the upper curves were shifted by multiplication with a factor. (b) Fitted core–shell Mie profiles using the code from Shen.82 The refractive index of the core was kept constant to n=1.004+i·0.045 [values of neutral xenon at 91 eV (Ref. 84)]. See text for details. Analog to (a), the profiles II–IV were shifted by a multiplicative factor for better visibility. Dashed lines in (a) and (b) show the profiles below an angle of 10°, where the experimental data were excluded from the fitting process. (c) Parameters of the shell obtained from the fitting, i.e., shell thickness d (in nm), absorption index β, and refractive index decrement δ. (d) Visualization of the sequence of core–shell structures derived from the fitting with changing parameters of the shell [for the exact values of the refractive indices, compare with the 2D color map or with the graphs for β and δ given in (c)].

FIG. 4.

(a) Difference profiles from the averaged profiles shown in Fig. 1(d). For better visibility, the upper curves were shifted by multiplication with a factor. (b) Fitted core–shell Mie profiles using the code from Shen.82 The refractive index of the core was kept constant to n=1.004+i·0.045 [values of neutral xenon at 91 eV (Ref. 84)]. See text for details. Analog to (a), the profiles II–IV were shifted by a multiplicative factor for better visibility. Dashed lines in (a) and (b) show the profiles below an angle of 10°, where the experimental data were excluded from the fitting process. (c) Parameters of the shell obtained from the fitting, i.e., shell thickness d (in nm), absorption index β, and refractive index decrement δ. (d) Visualization of the sequence of core–shell structures derived from the fitting with changing parameters of the shell [for the exact values of the refractive indices, compare with the 2D color map or with the graphs for β and δ given in (c)].

Close modal

With a Mie code that was extended for spheres with a core–shell structure,79,80,82 profiles were simulated and fitted to the experimental curves [see Fig. 4(b)]. A global optimization was carried out using a differential genetic algorithm.83 This approach, even if it is among the slowest optimization methods, has the advantage of being flexible, and it does not require an initial guess for the solution. The scattering angle in a range between 10° and 30° was selected for the fitting because smaller angles are prone to detector saturation effects [see single profiles in Fig. 1(c)]. The optimization target was set to the L2-norm of the logarithm of the profiles (i.e., least squares minimization on the logarithmic scale). The shell thickness d, the refraction decrement δshell, and the absorption index βshell were varied within the parameter bounds given in Table I. For the core, the values of neutral xenon at 91 eV photon energy84 were used (see Table I). It is noted that the fitting results were rather insensitive to the choice of the core's refractive index, as long as the absorption of the core was considerably higher than that of the shell. Instead of using a single particle size, the average profiles from several spheres with different sizes, following a distribution given by the scaling parameter σ, were fitted to the experimental profiles. This allows us to make the fitting insensitive to influences from the high-frequency oscillations of the profiles corresponding to the cluster size by taking into account the individual particles' deviation from spherical shape and the difference in sizes. For further details and observations on the fitting procedure, see the supplementary material.

TABLE I.

Fixed values and parameter bounds of the input parameters for the fitting routine. The literature values of neutral solid xenon at 91 eV are taken from Henke tables.84 

Literature values δXe,91eV −0.004 
 βXe,91eV 0.045 
Parameter bounds dshell (0, 100) (nm) 
 δshell (−0.1, 0.1) 
 βshell (0, 0.1) 
 σshell (0, 50) (nm) 
Literature values δXe,91eV −0.004 
 βXe,91eV 0.045 
Parameter bounds dshell (0, 100) (nm) 
 δshell (−0.1, 0.1) 
 βshell (0, 0.1) 
 σshell (0, 50) (nm) 

The parameters of the shell resulting from the global optimization are given in Fig. 4(c). We find (i) a compared to the core low value of β around 0.02, which indicates a transparent shell, (ii) an increasing shell thickness of 20–55 nm, and (iii) a refractive index decrement δ increasing with intensity with positive values up to 0.08. The same fitting results are obtained for a broad range of starting values (cf. Table I). Systematic Mie simulations, in which the three parameters are varied separately, are shown in the supplementary material. They support the assumption of a single optimal parameter set, as they show that the three parameters d, β, and δ have a separable influence on different characteristics of the core–shell signature. The fitting results also provide a scaling factor for the incoming number of photons and thus enable an estimate of the incoming intensities for the average profiles A to D in Fig. 1(d). Assigning the maximum irradiation intensity, ID=3×1014 W/cm2 to group D, the values for the other average profiles can be determined according to IA=6×1013 W/cm2, IB=1×1014 W/cm2, and IC=2×1014 W/cm2.

Our one-dimensional first-order ionization model indicates that the transient characteristic length, revealed by the observed modulations, is connected to the formation and evolution of a highly charged outer shell with dramatically changed optical properties. The fact that we are able to fit the experimental profiles well with a sequence of concentric core–shell Mie models further supports the general physical picture, even though both models are clearly limited. For gaining a clearer picture, a full description of the light propagation via sophisticated many-particle simulations will be needed, which includes also other processes such as impact ionization, charge transfer, plasma shifts of the energy levels, and further nanoplasma dynamics.85–87 Nevertheless, our general considerations provide a first step toward understanding the observed results.

The results from the Mie-calculations are in good qualitative agreement with the atomistic ionization model. They show the same general trend of a rather transparent outer shell, whose thickness increases with the irradiation level. From the Mie-fitting, we further extract a tendency for the real part of the refractive index, indicating a strong change in the refraction between the core and the shell, which grows with the incoming intensity. This general agreement in combination with the robustness of the fitting results indicates that key trends are captured in the analysis.

On the other hand, the large differences in the absolute values of shell thickness and absorption index demonstrate the limitations of our modeling approaches. In addition, while the Mie approach assumes a discontinuous transition between the core and the shell, the ionization model indicates a transition range on the order of 50 nm thickness. This may partly explain the higher absorption values in the Mie fits since a higher absorption can in principle take over the effect of a gradual transition (see the variation of beta in the supplementary material). Nevertheless, a transition region of such broadness cannot actually explain the formation of an observable scattering feature. We have tested with a simple scattering simulation (see the supplementary material) that the modulation in the scattering profiles disappears already for a transition region of half that thickness. The required sharp transition is puzzling and cannot be explained by our simple models.

The following considerations may allow us to hypothesize on the origin of a narrow reflective layer tens of nanometers deep in the nanoplasma. In general, a high reflectivity is connected to a strong change in the real part of the refractive index. The above discussed first-order ionization simulation only describes the absorption of the nanoplasma, i.e., the imaginary part β, and a model of its radial dependence. For the optical response of the cluster, both the real and imaginary parts of the refractive index, 1δ and β, are relevant, which are interrelated through the Kramers–Kronig dispersion relations. In this context, the optical properties of the charge state Xe3+ are worth a closer look. Plasma calculations of the atomic scattering factors of Xe3+ indicate that between 90 eV and 98 eV, the real part of the atomic scattering factor f1 (proportional to the refractive index decrement δ) rapidly changes from strongly positive to negative values and back several times.88 Now, we have to take into account that in the environment of the nanoplasma, the Xe3+ ions are not isolated (as in the gas-phase measurements carried out to determine the atomic absorption cross sections74) but instead surrounded by other ionic species and quasi-free electrons in the nanoplasma. By comparing the Xe3+-distribution [orange curve in Fig. 3(a)], peaking between 150 and 200 nm in depth with the average charge state at the same x-positions [magenta curve in Fig. 3(b)], we find that the neighborhood of the Xe3+ ions drastically changes. In a plasma environment, atomic or ionic resonances can be shifted in energy up to several eV.89 It is to be expected that the change in the environment as a function of the propagation depth translates into an energy shift, possibly from just below a sharp resonance to just above the resonance. This would result in a drastic change of the real part of the refractive index within a short distance, acting like a transient plasma mirror. A similar argument could be made for other xenon charge states Xe4+ to Xe6+, which also exhibit narrow and very strong absorption resonances in the vicinity of 91 eV.74,75 For testing this hypothesis, more sophisticated theoretical approaches will be required.

In summary, we have presented scattering patterns of single large xenon clusters resonantly excited with intense XUV pulses. The patterns reveal strong intensity dependent modulations in the radial distribution of the diffracted light, indicating the formation of an additional characteristic scattering length scale in an otherwise homogeneous nanoscale particle. Using a first-order modeling of the ionization in combination with Mie-based simulations, we assigned the transient diffraction signal to light induced electronic core–shell structures with an increasingly thick outer shell of low absorption and high refraction. An abrupt change in refraction, needed to explain the prominent diffraction feature observed, may be correlated with the radially changing plasma environment of higher charge states, translating into a radially changing shift of the electronic resonances. Our work shows that ultrafast resonant light scattering can map the transient spatial charge distributions in laser-excited nanoscale matter. The method can be employed to develop a deeper understanding of nanoplasma formation and charge transfer dynamics, which play a key role in many areas ranging from single-shot x-ray imaging to fusion and warm dense matter research and condensed matter physics. In the future, the approach provides an avenue to resolve ultrafast electron dynamics in extended systems on their natural timescale with intense attosecond pulses currently under development at FELs and lab-based sources.90–95 

See the supplementary material for a mathematical expression of the discretization approach, further details on the fitting procedure, a systematic investigation of the influence of each input parameter on the modulation features, simulations for sharp and smooth interfaces, and the derivation of the radial absorption index.

The authors thank Thomas Fennel and Christian Peltz for enlightening and helpful discussions. Excellent support from IOAP and DESY machine shops is acknowledged. The experiments received funding from BMBF (Grant No. 05K13KT2/05K16KT2) and DFG (Grant Nos. MO 719/13-1, 14-1, and 14-2). Further support from the Leibniz Grant (No. SAW/2017/MBI4; D.R. and M.S.), the NCCR MUST of the SNF (D.R. and A.C.), the Australian Research Council Centre of Excellence in Advanced Molecular Imaging (No. CE140100011; L.F.), the Volkswagen Foundation via a Peter Ewald fellowship (T.G.)., the U.S. Department of Energy, Office of Science, Basic Energy Sciences, Chemical Sciences, Geosciences, and Biosciences Division, the project Structural Dynamics of Biomolecular Systems (No. CZ.02.1.01/0.0/0.0/15_003/0000447) from European Regional Development Fund, and the project Advanced Research Using High Intensity Laser Produced Photons and Particles (No. CZ.02.1.01/0.0/0.0/16_019/0000789) from European Regional Development Fund (M.K.) is acknowledged.

The data that support the findings of this study are openly available at the CXI data bank at http://cxidb.org/id-146.html (reference number, 146).

1.
J.
Feldhaus
,
M.
Krikunova
,
M.
Meyer
,
T.
Möller
,
R.
Moshammer
,
A.
Rudenko
,
T.
Tschentscher
, and
J.
Ullrich
, “
AMO science at the FLASH and European XFEL free-electron laser facilities
,”
J. Phys. B
46
,
164002
(
2013
).
2.
C.
Bostedt
,
S.
Boutet
,
D. M.
Fritz
,
Z.
Huang
,
H. J.
Lee
,
H. T.
Lemke
,
A.
Robert
,
W. F.
Schlotter
,
J. J.
Turner
, and
G. J.
Williams
,
Rev. Mod. Phys.
88
,
015007
(
2016
).
3.
F.
Bencivenga
,
F.
Capotondi
,
E.
Principi
,
M.
Kiskinova
, and
C.
Masciovecchio
, “
Coherent and transient states studied with extreme ultraviolet and x-ray free electron lasers: Present and future prospects
,”
Adv. Phys.
63
,
327
404
(
2014
).
4.
M.
Yabashi
,
H.
Tanaka
,
T.
Tanaka
,
H.
Tomizawa
,
T.
Togashi
,
M.
Nagasono
,
T.
Ishikawa
,
J. R.
Harries
,
Y.
Hikosaka
,
A.
Hishikawa
,
K.
Nagaya
,
N.
Saito
,
E.
Shigemasa
,
K.
Yamanouchi
, and
K.
Ueda
, “
Compact XFEL and AMO sciences: SACLA and SCSS
,”
J. Phys. B
46
,
164001
(
2013
).
5.
H. N.
Chapman
,
A.
Barty
,
M. J.
Bogan
,
S.
Boutet
,
M.
Frank
,
S.
Hau-Riege
,
S.
Marchesini
,
B.
Woods
,
S.
Bajt
,
W. H.
Benner
,
R. A.
London
,
E.
Plonjes
,
M.
Kuhlmann
,
R.
Treusch
,
S.
Dusterer
,
T.
Tschentscher
,
J. R.
Schneider
,
E.
Spiller
,
T.
Möller
,
C.
Bostedt
,
M.
Hoener
,
D. A.
Shapiro
,
K. O.
Hodgson
,
D.
van der Spoel
,
F.
Burmeister
,
M.
Bergh
,
C.
Caleman
,
G.
Huldt
,
M.
Seibert
,
F. R. N. C.
Maia
,
R. W.
Lee
,
A.
Szoke
,
N.
Timneanu
, and
J.
Hajdu
, “
Femtosecond diffractive imaging with a soft-x-ray free-electron laser
,”
Nat. Phys.
2
,
839
843
(
2006
).
6.
R.
Neutze
 et al, “
Potential for biomolecular imaging with femtosecond X-ray pulses
,”
Nature
406
,
752
757
(
2000
).
7.
H. N.
Chapman
,
C.
Caleman
, and
N.
Timneanu
, “
Diffraction before destruction
,”
Philos. Trans. R. Soc. B
369
,
20130313
(
2014
).
8.
J. R.
Fienup
, “
Phase retrieval algorithms: A comparison
,”
Appl. Opt.
21
,
2758
2769
(
1982
).
9.
I.
Barke
,
H.
Hartmann
,
D.
Rupp
,
L.
Flückiger
,
M.
Sauppe
,
M.
Adolph
,
S.
Schorb
,
C.
Bostedt
,
R.
Treusch
,
C.
Peltz
,
S.
Bartling
,
T.
Fennel
,
K.
Meiwes-Broer
, and
T.
Möller
,
Nat. Commun.
6
,
6187
(
2015
).
10.
M.
Seibert
 et al, “
Single mimivirus particles intercepted and imaged with an X-ray laser
,”
Nature
470
,
78
81
(
2011
).
11.
T.
Ekeberg
,
M.
Svenda
,
C.
Abergel
,
F. R. N. C.
Maia
,
V.
Seltzer
,
J.-M.
Claverie
,
M.
Hantke
,
O.
Jönsson
,
C.
Nettelblad
,
G.
van der Schot
,
M.
Liang
,
D. P.
DePonte
,
A.
Barty
,
M. M.
Seibert
,
B.
Iwan
,
I.
Andersson
,
N. D.
Loh
,
A. V.
Martin
,
H.
Chapman
,
C.
Bostedt
,
J. D.
Bozek
,
K. R.
Ferguson
,
J.
Krzywinski
,
S. W.
Epp
,
D.
Rolles
,
A.
Rudenko
,
R.
Hartmann
,
N.
Kimmel
, and
J.
Hajdu
, “
Three-dimensional reconstruction of the giant mimivirus particle with an x-ray free-electron laser
,”
Phys. Rev. Lett.
114
,
098102
(
2015
).
12.
M.
Bogan
 et al, “
Aerosol imaging with a soft X-ray free electron laser
,”
Aerosol Sci. Technol.
44
,
1
6
(
2010
).
13.
N.
Loh
,
C. Y.
Hampton
,
A. V.
Martin
,
D.
Starodub
,
R. G.
Sierra
,
A.
Barty
,
A.
Aquila
,
J.
Schulz
,
L.
Lomb
,
J.
Steinbrener
 et al, “
Fractal morphology, imaging and mass spectrometry of single aerosol particles in flight
,”
Nature
486
,
513
517
(
2012
).
14.
C.
Bostedt
,
M.
Adolph
,
E.
Eremina
,
M.
Hoener
,
D.
Rupp
,
S.
Schorb
,
H.
Thomas
,
A.
de Castro
, and
T.
Möller
, “
Clusters in intense flash pulses: Ultrafast ionization dynamics and electron emission studied with spectroscopic and scattering techniques
,”
J. Phys. B
43
,
194011
(
2010
).
15.
D.
Rupp
,
M.
Adolph
,
T.
Gorkhover
,
S.
Schorb
,
D.
Wolter
,
R.
Hartmann
,
N.
Kimmel
,
C.
Reich
,
T.
Feigl
,
A.
de Castro
,
R.
Treusch
,
L.
Strüder
,
T.
Möller
, and
C.
Bostedt
, “
Identification of twinned gas phase clusters by single-shot scattering with intense soft x-ray pulses
,”
New J. Phys.
14
,
055016
(
2012
).
16.
L.
Gomez
,
K.
Ferguson
,
J.
Cryan
,
C.
Bacellar
,
R.
Tanyag
,
C.
Jones
,
S.
Schorb
,
D.
Anielski
,
A.
Belkacem
,
C.
Bernando
,
R.
Boll
,
J.
Bozek
,
S.
Carron
,
G.
Chen
,
T.
Delmas
,
L.
Englert
,
S. W.
Epp
,
B.
Erk
,
L.
Foucar
,
R.
Hartmann
,
A.
Hexemer
,
M.
Huth
,
J.
Kwok
,
S. R.
Leone
,
J. H. S.
Ma
,
F. R. N. C.
Maia
,
E.
Malmerberg
,
S.
Marchesini
,
D. M.
Neumark
,
B.
Poon
,
J.
Prell
,
D.
Rolles
,
B.
Rudek
,
A.
Rudenko
,
M.
Seifrid
,
K. R.
Siefermann
,
F. P.
Sturm
,
M.
Swiggers
,
J.
Ullrich
,
F.
Weise
,
P.
Zwart
,
C.
Bostedt
,
O.
Gessner
, and
A. F.
Vilesov
, “
Shapes and vorticities of superfluid helium nanodroplets
,”
Science
345
,
906
909
(
2014
).
17.
B.
Langbehn
,
K.
Sander
,
Y.
Ovcharenko
,
C.
Peltz
,
A.
Clark
,
M.
Coreno
,
R.
Cucini
,
M.
Drabbels
,
P.
Finetti
,
M.
Di Fraia
,
L.
Giannessi
,
C.
Grazioli
,
D.
Iablonskyi
,
A. C.
LaForge
,
T.
Nishiyama
,
V.
Oliver Álvarez de Lara
,
P.
Piseri
,
O.
Plekan
,
K.
Ueda
,
J.
Zimmermann
,
K. C.
Prince
,
F.
Stienkemeier
,
C.
Callegari
,
T.
Fennel
,
D.
Rupp
, and
T.
Möller
, “
Three-dimensional shapes of spinning helium nanodroplets
,”
Phys. Rev. Lett.
121
,
255301
(
2018
).
18.
J. N.
Clark
,
L.
Beitra
,
G.
Xiong
,
A.
Higginbotham
,
D. M.
Fritz
,
H. T.
Lemke
,
D.
Zhu
,
M.
Chollet
,
G. J.
Williams
,
M.
Messerschmidt
,
B.
Abbey
,
R. J.
Harder
,
A. M.
Korsunsky
,
J. S.
Wark
, and
I. K.
Robinson
, “
Ultrafast three-dimensional imaging of lattice dynamics in individual gold nanocrystals
,”
Science
341
,
56
59
(
2013
).
19.
T.
Gorkhover
,
S.
Schorb
,
R.
Coffee
,
M.
Adolph
,
L.
Foucar
,
D.
Rupp
,
A.
Aquila
,
J. D.
Bozek
,
S. W.
Epp
,
B.
Erk
,
L.
Gumprecht
,
L.
Holmegaard
,
A.
Hartmann
,
R.
Hartmann
,
G.
Hauser
,
P.
Holl
,
A.
Hömke
,
P.
Johnsson
,
N.
Kimmel
,
K.-U.
Kühnel
,
M.
Messerschmidt
,
C.
Reich
,
A.
Rouzée
,
B.
Rudek
,
C.
Schmidt
,
J.
Schulz
,
H.
Soltau
,
S.
Stern
,
G.
Weidenspointner
,
B.
White
,
J.
Küpper
,
L.
Strüder
,
I.
Schlichting
,
J.
Ullrich
,
D.
Rolles
,
A.
Rudenko
,
T.
Möller
, and
C.
Bostedt
, “
Femtosecond and nanometre visualization of structural dynamics in superheated nanoparticles
,”
Nat. Photonics
10
,
93
(
2016
).
20.
L.
Flückiger
,
D.
Rupp
,
M.
Adolph
,
T.
Gorkhover
,
M.
Krikunova
,
M.
Müller
,
T.
Oelze
,
Y.
Ovcharenko
,
M.
Sauppe
,
S.
Schorb
,
C.
Bostedt
,
S.
Düsterer
,
M.
Harmand
,
H.
Redlin
,
R.
Treusch
, and
T.
Möller
, “
Time-resolved x-ray imaging of a laser-induced nanoplasma and its neutral residuals
,”
New J. Phys.
18
,
043017
(
2016
).
21.
K. R.
Ferguson
,
M.
Bucher
,
T.
Gorkhover
,
S.
Boutet
,
H.
Fukuzawa
,
J. E.
Koglin
,
Y.
Kumagai
,
A.
Lutman
,
A.
Marinelli
,
M.
Messerschmidt
,
K.
Nagaya
,
J.
Turner
,
K.
Ueda
,
G. J.
Williams
,
P. H.
Bucksbaum
, and
C.
Bostedt
, “
Transient lattice contraction in the solid-to-plasma transition
,”
Sci. Adv.
2
,
e1500837
(
2016
).
22.
T.
Nishiyama
,
Y.
Kumagai
,
A.
Niozu
,
H.
Fukuzawa
,
K.
Motomura
,
M.
Bucher
,
Y.
Ito
,
T.
Takanashi
,
K.
Asa
,
Y.
Sato
,
D.
You
,
Y.
Li
,
T.
Ono
,
E.
Kukk
,
C.
Miron
,
L.
Neagu
,
C.
Callegari
,
M.
Di Fraia
,
G.
Rossi
,
D. E.
Galli
,
T.
Pincelli
,
A.
Colombo
,
T.
Kameshima
,
Y.
Joti
,
T.
Hatsui
,
S.
Owada
,
T.
Katayama
,
T.
Togashi
,
K.
Tono
,
M.
Yabashi
,
K.
Matsuda
,
C.
Bostedt
,
K.
Nagaya
, and
K.
Ueda
, “
Ultrafast structural dynamics of nanoparticles in intense laser fields
,”
Phys. Rev. Lett.
123
,
123201
(
2019
).
23.
B.
Ziaja
,
H.
Chapman
,
R.
Fäustlin
,
S.
Hau-Riege
,
Z.
Jurek
,
A.
Martin
,
S.
Toleikis
,
F.
Wang
,
E.
Weckert
, and
R.
Santra
, “
Limitations of coherent diffractive imaging of single objects due to their damage by intense x-ray radiation
,”
New J. Phys.
14
,
115015
(
2012
).
24.
S.-K.
Son
,
L.
Young
, and
R.
Santra
, “
Impact of hollow-atom formation on coherent x-ray scattering at high intensity
,”
Phys. Rev. A
83
,
033402
(
2011
).
25.
H. M.
Quiney
and
K. A.
Nugent
, “
Biomolecular imaging and electronic damage using x-ray free-electron lasers
,”
Nat. Phys.
7
,
142
146
(
2011
).
26.
K.
Nass
,
L.
Foucar
,
T. R. M.
Barends
,
E.
Hartmann
,
S.
Botha
,
R. L.
Shoeman
,
R. B.
Doak
,
R.
Alonso-Mori
,
A.
Aquila
,
S.
Bajt
,
A.
Barty
,
R.
Bean
,
K. R.
Beyerlein
,
M.
Bublitz
,
N.
Drachmann
,
J.
Gregersen
,
H. O.
Jönsson
,
W.
Kabsch
,
S.
Kassemeyer
,
J. E.
Koglin
,
M.
Krumrey
,
D.
Mattle
,
M.
Messerschmidt
,
P.
Nissen
,
L.
Reinhard
,
O.
Sitsel
,
D.
Sokaras
,
G. J.
Williams
,
S.
Hau-Riege
,
N.
Timneanu
,
C.
Caleman
,
H. N.
Chapman
,
S.
Boutet
, and
I.
Schlichting
, “
Indications of radiation damage in ferredoxin microcrystals using high-intensity X-FEL beams
,”
J. Synchrotron Radiat.
22
,
225
(
2015
).
27.
P. J.
Ho
,
B. J.
Daurer
,
M. F.
Hantke
,
J.
Bielecki
,
A.
Al Haddad
,
M.
Bucher
,
G.
Doumy
,
K. R.
Ferguson
,
L.
Flückiger
,
T.
Gorkhover
,
B.
Iwan
,
C.
Knight
,
S.
Moeller
,
T.
Osipov
,
D.
Ray
,
S. H.
Southworth
,
M.
Svenda
,
N.
Timneanu
,
A.
Ulmer
,
P.
Walter
,
J.
Hajdu
,
L.
Young
,
F. R. N. C.
Maia
, and
C.
Bostedt
, “
The role of transient resonances for ultra-fast imaging of single sucrose nanoclusters
,”
Nat. Commun.
11
,
167
(
2020
).
28.
C.
Bostedt
,
E.
Eremina
,
D.
Rupp
,
M.
Adolph
,
H.
Thomas
,
M.
Hoener
,
A. R. B.
de Castro
,
J.
Tiggesbäumker
,
K.-H.
Meiwes-Broer
,
T.
Laarmann
,
H.
Wabnitz
,
E.
Plönjes
,
R.
Treusch
,
J. R.
Schneider
, and
T.
Möller
, “
Ultrafast x-ray scattering of xenon nanoparticles: Imaging transient states of matter
,”
Phys. Rev. Lett.
108
,
093401
(
2012
).
29.
S.-K.
Son
,
H. N.
Chapman
, and
R.
Santra
, “
Multiwavelength anomalous diffraction at high x-ray intensity
,”
Phys. Rev. Lett.
107
,
218102
(
2011
).
30.
C.
Song
,
R.
Bergstrom
,
D.
Ramunno-Johnson
,
H.
Jiang
,
D.
Paterson
,
M. D.
de Jonge
,
I.
McNulty
,
J.
Lee
,
K. L.
Wang
, and
J.
Miao
, “
Nanoscale imaging of buried structures with elemental specificity using resonant x-ray diffraction microscopy
,”
Phys. Rev. Lett.
100
,
025504
(
2008
).
31.
W.
Ackermann
,
G.
Asova
,
V.
Ayvazyan
,
A.
Azima
,
N.
Baboi
,
J.
Bahr
,
V.
Balandin
,
B.
Beutner
,
A.
Brandt
,
A.
Bolzmann
,
R.
Brinkmann
,
O.
Brovko
,
M.
Castellano
,
P.
Castro
,
L.
Catani
,
E.
Chiadroni
,
S.
Choroba
,
A.
Cianchi
,
J.
Costello
,
D.
Cubaynes
,
J.
Dardis
,
W.
Decking
,
H.
Delsim-Hashemi
,
A.
Delserieys
,
G.
DiPirro
,
M.
Dohlus
,
S.
Dusterer
,
A.
Eckhardt
,
H.
Edwards
,
B.
Faatz
,
J.
Feldhaus
,
K.
Flottmann
,
J.
Frisch
,
L.
Frohlich
,
T.
Garvey
,
U.
Gensch
,
C.
Gerth
,
M.
Gorler
,
N.
Golubeva
,
H.-J.
Grabosch
,
M.
Grecki
,
O.
Grimm
,
K.
Hacker
,
U.
Hahn
,
J.
Han
,
K.
Honkavaara
,
T.
Hott
,
M.
Huning
,
Y.
Ivanisenko
,
E.
Jaeschke
,
W.
Jalmuzna
,
T.
Jezynski
,
R.
Kammering
,
V.
Katalev
,
K.
Kavanagh
,
E.
Kennedy
,
S.
Khodyachykh
,
K.
Klose
,
V.
Kocharyan
,
M.
Korfer
,
M.
Kollewe
,
W.
Koprek
,
S.
Korepanov
,
D.
Kostin
,
M.
Krassilnikov
,
G.
Kube
,
M.
Kuhlmann
,
C.
Lewis
,
L.
Lilje
,
T.
Limberg
,
D.
Lipka
,
F.
Lohl
,
H.
Luna
,
M.
Luong
,
M.
Martins
,
M.
Meyer
,
P.
Michelato
,
V.
Miltchev
,
W.
Moller
,
L.
Monaco
,
W.
Muller
,
O.
Napieralski
,
O.
Napoly
,
P.
Nicolosi
,
D.
Nolle
,
T.
Nunez
,
A.
Oppelt
,
C.
Pagani
,
R.
Paparella
,
N.
Pchalek
,
J.
Pedregosa-Gutierrez
,
B.
Petersen
,
B.
Petrosyan
,
G.
Petrosyan
,
L.
Petrosyan
,
J.
Pfluger
,
E.
Plonjes
,
L.
Poletto
,
K.
Pozniak
,
E.
Prat
,
D.
Proch
,
P.
Pucyk
,
P.
Radcliffe
,
H.
Redlin
,
K.
Rehlich
,
M.
Richter
,
M.
Roehrs
,
J.
Roensch
,
R.
Romaniuk
,
M.
Ross
,
J.
Rossbach
,
V.
Rybnikov
,
M.
Sachwitz
,
E.
Saldin
,
W.
Sandner
,
H.
Schlarb
,
B.
Schmidt
,
M.
Schmitz
,
P.
Schmuser
,
J.
Schneider
,
E.
Schneidmiller
,
S.
Schnepp
,
S.
Schreiber
,
M.
Seidel
,
D.
Sertore
,
A.
Shabunov
,
C.
Simon
,
S.
Simrock
,
E.
Sombrowski
,
A.
Sorokin
,
P.
Spanknebel
,
R.
Spesyvtsev
,
L.
Staykov
,
B.
Steffen
,
F.
Stephan
,
F.
Stulle
,
H.
Thom
,
K.
Tiedtke
,
M.
Tischer
,
S.
Toleikis
,
R.
Treusch
,
D.
Trines
,
I.
Tsakov
,
E.
Vogel
,
T.
Weiland
,
H.
Weise
,
M.
Wellhofer
,
M.
Wendt
,
I.
Will
,
A.
Winter
,
K.
Wittenburg
,
W.
Wurth
,
P.
Yeates
,
M.
Yurkov
,
I.
Zagorodnov
, and
K.
Zapfe
, “
Operation of a free-electron laser from the extreme ultraviolet to the water window
,”
Nat. Photonics
1
,
336
342
(
2007
).
32.
K.
Tiedtke
 et al,
New J. Phys.
11
,
023029
(
2009
).
33.
C.
Behrens
,
N.
Gerasimova
,
C.
Gerth
,
B.
Schmidt
,
E.
Schneidmiller
,
S.
Serkez
,
S.
Wesch
, and
M.
Yurkov
, “
Constraints on photon pulse duration from longitudinal electron beam diagnostics at a soft x-ray free-electron laser
,”
Phys. Rev. Spec. Top.-Accel. Beams
15
,
030707
(
2012
).
34.
S.
Düsterer
,
M.
Rehders
,
A.
Al-Shemmary
,
C.
Behrens
,
G.
Brenner
,
O.
Brovko
,
M.
DellAngela
,
M.
Drescher
,
B.
Faatz
,
J.
Feldhaus
,
U.
Frühling
,
N.
Gerasimova
,
N.
Gerken
,
C.
Gerth
,
T.
Golz
,
A.
Grebentsov
,
E.
Hass
,
K.
Honkavaara
,
V.
Kocharian
,
M.
Kurka
,
T.
Limberg
,
R.
Mitzner
,
R.
Moshammer
,
E.
Plönjes
,
M.
Richter
,
J.
Rönsch-Schulenburg
,
A.
Rudenko
,
H.
Schlarb
,
B.
Schmidt
,
A.
Senftleben
,
E. A.
Schneidmiller
,
B.
Siemer
,
F.
Sorgenfrei
,
A. A.
Sorokin
,
N.
Stojanovic
,
K.
Tiedtke
,
R.
Treusch
,
M.
Vogt
,
M.
Wieland
,
W.
Wurth
,
S.
Wesch
,
M.
Yan
,
M. V.
Yurkov
,
H.
Zacharias
, and
S.
Schreiber
, “
Development of experimental techniques for the characterization of ultrashort photon pulses of extreme ultraviolet free-electron lasers
,”
Phys. Rev. Spec. Top.-Accel. Beams
17
,
120702
(
2014
).
35.
The name LOLA is for historic reasons, since G. A. Loew and R. R. Larsen and O. H. Altenmueller used such a tranverse deflecting structure in the 1960's to separate particles as described in .
36.
D.
Rupp
,
M.
Adolph
,
L.
Flückiger
,
T.
Gorkhover
,
J. P.
Müller
,
M.
Müller
,
M.
Sauppe
,
D.
Wolter
,
S.
Schorb
,
R.
Treusch
,
C.
Bostedt
, and
T.
Möller
, “
Generation and structure of extremely large clusters in pulsed jets
,”
J. Chem. Phys.
141
,
044306
(
2014
).
37.
T.
Gorkhover
,
M.
Adolph
,
D.
Rupp
,
S.
Schorb
,
S. W.
Epp
,
B.
Erk
,
L.
Foucar
,
R.
Hartmann
,
N.
Kimmel
,
K.-U.
Kühnel
,
D.
Rolles
,
B.
Rudek
,
A.
Rudenko
,
R.
Andritschke
,
A.
Aquila
,
J. D.
Bozek
,
N.
Coppola
,
T.
Erke
,
F.
Filsinger
,
H.
Gorke
,
H.
Graafsma
,
L.
Gumprecht
,
G.
Hauser
,
S.
Herrmann
,
H.
Hirsemann
,
A.
Hömke
,
P.
Holl
,
C.
Kaiser
,
F.
Krasniqi
,
J.-H.
Meyer
,
M.
Matysek
,
M.
Messerschmidt
,
D.
Miessner
,
B.
Nilsson
,
D.
Pietschner
,
G.
Potdevin
,
C.
Reich
,
G.
Schaller
,
C.
Schmidt
,
F.
Schopper
,
C. D.
Schröter
,
J.
Schulz
,
H.
Soltau
,
G.
Weidenspointner
,
I.
Schlichting
,
L.
Strüder
,
J.
Ullrich
,
T.
Möller
, and
C.
Bostedt
, “
Nanoplasma dynamics of single large xenon clusters irradiated with superintense x-ray pulses from the linac coherent light source free-electron laser
,”
Phys. Rev. Lett.
108
,
245005
(
2012
).
38.
D.
Rupp
,
L.
Flückiger
,
M.
Adolph
,
T.
Gorkhover
,
M.
Krikunova
,
J. P.
Müller
,
M.
Müller
,
T.
Oelze
,
Y.
Ovcharenko
,
B.
Röben
,
M.
Sauppe
,
S.
Schorb
,
D.
Wolter
,
R.
Mitzner
,
M.
Wöstmann
,
S.
Roling
,
M.
Harmand
,
R.
Treusch
,
M.
Arbeiter
,
T.
Fennel
,
C.
Bostedt
, and
T.
Möller
, “
Recombination-enhanced surface expansion of clusters in intense soft x-ray laser pulses
,”
Phys. Rev. Lett.
117
,
153401
(
2016
).
39.
F. R.
Maia
, “
The coherent x-ray imaging data bank
,”
Nat. Methods
9
,
854
(
2012
).
40.
See http://cxidb.org/id-146.html for “Web Link for the Raw Data on CDI Data Bank,
2020
.”
41.
C.
Siedschlag
and
J.-M.
Rost
, “
Small rare-gas clusters in soft x-ray pulses
,”
Phys. Rev. Lett.
93
,
043402
(
2004
).
42.
U.
Saalmann
,
C.
Siedschlag
, and
J.
Rost
, “
Mechanisms of cluster ionization in strong laser pulses
,”
J. Phys. B
39
,
R39
R77
(
2006
).
43.
M. R.
Islam
,
U.
Saalmann
, and
J. M.
Rost
, “
Kinetic energy of ions after coulomb explosion of clusters induced by an intense laser pulse
,”
Phys. Rev. A
73
,
041201
(
2006
).
44.
B.
Ziaja
,
H.
Wabnitz
,
E.
Weckert
, and
T.
Möller
, “
Atomic clusters of various sizes irradiated with short intense pulses of VUV radiation
,”
Europhys. Lett.
82
,
24002
(
2008
).
45.
B.
Ziaja
,
T.
Laarmann
,
H.
Wabnitz
,
F.
Wang
,
E.
Weckert
,
C.
Bostedt
, and
T.
Möller
, “
Emission of electrons from rare gas clusters after irradiation with intense VUV pulses of wavelength 100 nm and 32 nm
,”
New J. Phys.
11
,
103012
(
2009
).
46.
T.
Fennel
 et al, “
Laser-driven nonlinear cluster dynamics
,”
Rev. Mod. Phys.
82
,
1793
1842
(
2010
).
47.
M.
Arbeiter
and
T.
Fennel
, “
Ionization heating in rare-gas clusters under intense XUV laser pulses
,”
Phys. Rev. A
82
,
013201
(
2010
).
48.
M.
Arbeiter
and
T.
Fennel
, “
Rare-gas clusters in intense VUV, XUV and soft x-ray pulses: Signatures of the transition from nanoplasma-driven cluster expansion to Coulomb explosion in ion and electron spectra
,”
New J. Phys.
13
,
053022
(
2011
).
49.
E.
Ackad
,
N.
Bigaouette
,
K.
Briggs
, and
L.
Ramunno
, “
Clusters in intense XUV pulses: Effects of cluster size on expansion dynamics and ionization
,”
Phys. Rev. A
83
,
063201
(
2011
).
50.
E.
Ackad
,
N.
Bigaouette
,
S.
Mack
,
K.
Popov
, and
L.
Ramunno
, “
Recombination effects in soft-x-ray cluster interactions at the xenon giant resonance
,”
New J. Phys.
15
,
053047
(
2013
).
51.
M.
Arbeiter
,
C.
Peltz
, and
T.
Fennel
, “
Electron-relocalization dynamics in xenon clusters in intense soft-x-ray fields
,”
Phys. Rev. A
89
,
043428
(
2014
).
52.
I.
Last
and
J.
Jortner
, “
Electron–ion impact energy transfer in nanoplasmas of coulomb exploding clusters
,”
Advances in Quantum Chemistry
(
Elsevier
,
2017
), Vol.
75
, pp.
27
52
.
53.
I.
Last
and
J.
Jortner
, “
Microscopic and hydrodynamic impact energy transfer from nanoplasma electrons to ions in exploding clusters
,”
Mol. Phys.
116
,
2461
2473
(
2018
).
54.
H.
Wabnitz
 et al,
Nature
420
,
482
485
(
2002
).
55.
M.
Hoener
,
C.
Bostedt
,
H.
Thomas
,
L.
Landt
,
E.
Eremina
,
H.
Wabnitz
,
T.
Laarmann
,
R.
Treusch
,
A.
De Castro
, and
T.
Möller
, “
Charge recombination in soft x-ray laser produced nanoplasmas
,”
J. Phys. B
41
,
181001
(
2008
).
56.
C.
Bostedt
 et al, “
Multistep ionization of argon clusters in intense femtosecond extreme ultraviolet pulses
,”
Phys. Rev. Lett.
100
,
133401
(
2008
).
57.
H.
Thomas
 et al,
J. Phys. B
42
,
134018
(
2009
).
58.
C.
Bostedt
,
H.
Thomas
,
M.
Hoener
,
T.
Möller
,
U.
Saalmann
,
I.
Georgescu
,
C.
Gnodtke
, and
J.-M.
Rost
, “
Fast electrons from multi-electron dynamics in xenon clusters induced by inner-shell ionization
,”
New J. Phys.
12
,
083004
(
2010
).
59.
K.
Nagaya
 et al, “
Investigation of the interaction of xenon cluster with intense EUV–FEL pulses using pulsed cluster beam source and momentum imaging spectrometer
,”
J. Electron Spectrosc. Relat. Phenom.
181
,
125
128
(
2010
).
60.
H.
Iwayama
,
A.
Sugishima
,
K.
Nagaya
,
M.
Yao
,
H.
Fukuzawa
,
K.
Motomura
,
X.
Liu
,
A.
Yamada
,
C.
Wang
,
K.
Ueda
 et al, “
Inhomogeneous charge redistribution in Xe clusters exposed to an intense extreme ultraviolet free electron laser
,”
J. Phys. B
43
,
161001
(
2010
).
61.
S.
Schorb
 et al,
Phys. Rev. Lett.
108
,
233401
(
2012
).
62.
M.
Krikunova
,
M.
Adolph
,
T.
Gorkhover
,
D.
Rupp
,
S.
Schorb
,
C.
Bostedt
,
S.
Roling
,
B.
Siemer
,
R.
Mitzner
,
H.
Zacharias
 et al, “
Ionization dynamics in expanding clusters studied by XUV pump–probe spectroscopy
,”
J. Phys. B
45
,
105101
(
2012
).
63.
H.
Iwayama
,
K.
Nagaya
,
M.
Yao
,
H.
Fukuzawa
,
X.
Liu
,
G.
Prümper
,
K.
Motomura
,
K.
Ueda
,
N.
Saito
,
A.
Rudenko
 et al, “
Frustration of photoionization of Ar nanoplasma produced by extreme ultraviolet FEL pulses
,”
J. Phys. B
46
,
164019
(
2013
).
64.
L.
Schroedter
 et al,
Phys. Rev. Lett.
112
,
183401
(
2014
).
65.
B.
Schütte
,
M.
Arbeiter
,
T.
Fennel
,
M. J.
Vrakking
, and
A.
Rouzée
, “
Rare-gas clusters in intense extreme-ultraviolet pulses from a high-order harmonic source
,”
Phys. Rev. Lett.
112
,
073003
(
2014
).
66.
B.
Schütte
,
F.
Campi
,
M.
Arbeiter
,
T.
Fennel
,
M.
Vrakking
, and
A.
Rouzée
, “
Tracing electron-ion recombination in nanoplasmas produced by extreme-ultraviolet irradiation of rare-gas clusters
,”
Phys. Rev. Lett.
112
,
253401
(
2014
).
67.
B.
Schütte
,
M.
Arbeiter
,
A.
Mermillod-Blondin
,
M. J.
Vrakking
,
A.
Rouzée
, and
T.
Fennel
, “
Ionization avalanching in clusters ignited by extreme-ultraviolet driven seed electrons
,”
Phys. Rev. Lett.
116
,
033001
(
2016
).
68.
I.
Last
and
J.
Jortner
, “
Dynamics of the coulomb explosion of large clusters in a strong laser field
,”
Phys. Rev. A
62
,
013201
(
2000
).
69.
I.
Last
and
J.
Jortner
, “
Quasiresonance ionization of large multicharged clusters in a strong laser field
,”
Phys. Rev. A
60
,
2215
(
1999
).
70.
C.
Peltz
,
C.
Varin
,
T.
Brabec
, and
T.
Fennel
, “
Time-resolved x-ray imaging of anisotropic nanoplasma expansion
,”
Phys. Rev. Lett.
113
,
133401
(
2014
).
71.
D.
Ederer
, “
Photoionization of the 4d electrons in xenon
,”
Phys. Rev. Lett.
13
,
760
(
1964
).
72.
Y.
Itoh
,
A.
Ito
, and
M.
Kitajima
, “
Absolute photoionization cross section measurements of Xe+ ions in the 4d threshold energy region
,”
J. Phys. B
4075
,
3493
3499
(
2001
).
73.
P.
Andersen
,
T.
Andersen
,
F.
Folkmann
,
V.
Ivanov
,
H.
Kjeldsen
, and
J.
West
, “
Absolute cross sections for the photoionization of 4d electrons in Xe+ and Xe2+ ions
,”
J. Phys. B
4075
,
2009
2019
(
2001
).
74.
E. D.
Emmons
,
A.
Aguilar
,
M. F.
Gharaibeh
,
S. W. J.
Scully
,
R. A.
Phaneuf
,
A. L. D.
Kilcoyne
,
A. S.
Schlachter
,
I.
Álvarez
,
C.
Cisneros
, and
G.
Hinojosa
, “
Photoionization and electron-impact ionization of Xe3+
,”
Phys. Rev. A
71
,
042704
(
2005
).
75.
A.
Aguilar
,
J. D.
Gillaspy
,
G. F.
Gribakin
,
R. A.
Phaneuf
,
M. F.
Gharaibeh
,
M. G.
Kozlov
,
J. D.
Bozek
, and
A. L. D.
Kilcoyne
, “
Absolute photoionization cross sections for Xe4+, Xe5+, and Xe6+ near 13.5 nm: Experiment and theory
,”
Phys. Rev. A
73
,
1
10
(
2006
).
76.
J.
Bizau
,
J.-M.
Esteva
,
D.
Cubaynes
,
F. J.
Wuilleumier
,
C.
Blancard
,
A.
Compant La Fontaine
,
C.
Couillaud
,
J.
Lachkar
,
R.
Marmoret
,
C.
Remond
,
J.
Bruneau
,
D.
Hitz
,
P.
Ludwig
, and
M.
Delaunay
, “
Photoionization of highly charged ions using an ECR ion source and undulator radiation
,”
Phys. Rev. Lett.
84
,
435
438
(
2000
).
77.
D.
Sears
and
H.
Klug
, “
Density and expansivity of solid xenon
,”
J. Chem. Phys.
37
,
3002
(
1962
).
78.
D.
Attwood
,
Soft X-Ray and Extreme Ultraviolet Radiation
(
Cambridge University Press
,
2007
).
79.
G.
Mie
, “
Beiträge zur Optik trüber Medien, speziell kolloidaler Metallösungen
,”
Ann. Phys.
330
,
377
(
1908
).
80.
C.
Bohren
and
D.
Huffman
,
Absorption and Scattering of Light by Small Particles
(
Wiley
,
1983
).
81.
C.
von Korff Schmising
,
B.
Pfau
,
M.
Schneider
,
C. M.
Günther
,
M.
Giovannella
,
J.
Perron
,
B.
Vodungbo
,
L.
Müller
,
F.
Capotondi
,
E.
Pedersoli
,
N.
Mahne
,
J.
Lüning
, and
S.
Eisebitt
, “
Imaging ultrafast demagnetization dynamics after a spatially localized optical excitation
,”
Phys. Rev. Lett.
112
,
217203
(
2014
).
82.
See https://scattport.org/files/jianqi_shen/Mie%20CoatedParticleIntensity.rar for “Web Link for the Code Package by Jianqi Shen for Coated Spheres,
2014
.”
83.
R.
Storn
and
K.
Price
, “
Differential evolution: A simple and efficient heuristic for global optimization over continuous spaces
,”
J. Global Optim.
11
,
341
359
(
1997
).
84.
B. L.
Henke
,
E. M.
Gullikson
, and
J. C.
Davis
, “
X-ray interactions: Photoabsorption, scattering, transmission and reflection E = 50–30,000 eV, z = 1–92
,”
Atomic Data and Nuclear Data Tables
54
,
181
342
(
1993
).
85.
U.
Saalmann
, “
Cluster nanoplasmas in strong flash pulses: Formation, excitation and relaxation
,”
J. Phys. B
43
,
194012
(
2010
).
86.
C.
Peltz
,
C.
Varin
,
T.
Brabec
, and
T.
Fennel
, “
Fully microscopic analysis of laser-driven finite plasmas using the example of clusters
,”
New J. Phys.
14
,
065011
(
2012
).
87.
P. J.
Ho
and
C.
Knight
, “
Large-scale atomistic calculations of clusters in intense x-ray pulses
,”
J. Phys. B
50
,
104003
(
2017
).
88.
J.
Nilsen
,
W.
Johnson
, and
K.
Cheng
, “
Searching for plasmas with anomalous dispersion in the soft X-ray regime
,”
Proc. SPIE
6702
,
67020N
(
2007
).
89.
A.
Gets
and
V.
Krainov
, “
The ionization potentials of atomic ions in laser-irradiated Ar, Kr and Xe clusters
,”
J. Phys. B
39
,
1787
1795
(
2006
).
90.
A. A.
Zholents
and
G.
Penn
, “
Obtaining attosecond x-ray pulses using a self-amplified spontaneous emission free electron laser
,”
Phys. Rev. Spec. Top.-Accel. Beams
8
,
050704
(
2005
).
91.
E. L.
Saldin
,
E. A.
Schneidmiller
, and
M. V.
Yurkov
, “
Self-amplified spontaneous emission FEL with energy-chirped electron beam and its application for generation of attosecond x-ray pulses
,”
Phys. Rev. Spec. Top.-Accel. Beams
9
,
050702
(
2006
).
92.
R. N.
Coffee
,
J. P.
Cryan
,
J.
Duris
,
W.
Helml
,
S.
Li
, and
A.
Marinelli
, “
Development of ultrafast capabilities for x-ray free-electron lasers at the linac coherent light source
,”
Philos. Trans. R. Soc. A
377
,
20180386
(
2019
).
93.
J.
Duris
,
S.
Li
,
T.
Driver
,
E. G.
Champenois
,
J. P.
MacArthur
,
A. A.
Lutman
,
Z.
Zhang
,
P.
Rosenberger
,
J. W.
Aldrich
,
R.
Coffee
 et al, “
Tunable isolated attosecond x-ray pulses with gigawatt peak power from a free-electron laser
,”
Nat. Photonics
14
,
30
36
(
2020
).
94.
P.
Tzallas
,
E.
Skantzakis
,
L.
Nikolopoulos
,
G. D.
Tsakiris
, and
D.
Charalambidis
, “
Extreme-ultraviolet pump–probe studies of one-femtosecond-scale electron dynamics
,”
Nat. Phys.
7
,
781
784
(
2011
).
95.
E. J.
Takahashi
,
P.
Lan
,
O. D.
Mücke
,
Y.
Nabekawa
, and
K.
Midorikawa
, “
Attosecond nonlinear optics using gigawatt-scale isolated attosecond pulses
,”
Nat. Commun.
4
,
2691
(
2013
).

Supplementary Material