Non-thermal evolution of dense plasmas driven by intense x-ray fields
Abstract
The advent of x-ray free-electron lasers (XFELs) has enabled a range of new experimental investigations into the properties of matter driven to extreme conditions via intense x-ray-matter interactions. The femtosecond timescales of these interactions lead to the creation of transient high-energy-density plasmas, where both the electrons and the ions may be far from local thermodynamic equilibrium (LTE). Predictive modelling of such systems remains challenging because of the substantially different timescales on which electrons and ions thermalize, and because of the vast number of atomic configurations that are required to describe the resulting highly-ionized plasmas. Here we explore the evolution of systems driven to high energy densities using CCFLY, a non-LTE, Fokker-Planck collisional-radiative code. We use CCFLY to investigate the evolution dynamics of a solid-density plasma driven by an XFEL, and explore the relaxation of the plasma to local thermodynamic equilibrium on femtosecond timescales in terms of the charge state distribution, electron density, and temperature.
pacs:
Valid PACS appear hereI Introduction
Experiments carried out at x-ray free-electron lasers (FELs) over the past decade have shown how intense x-ray fields can produce well-characterized plasmas at high temperatures and densities via isochoric x-ray photoabsortion Vinko et al. (2012); Fletcher et al. (2015); Levy et al. (2015); Ciricosta et al. (2016a). This approach has proven valuable in providing access to physical properties and atomic processes in warm dense matter, of importance to applications across plasma physics Ciricosta et al. (2012); Vinko et al. (2014), astrophysics and planetary studies Eggert et al. (2010); Smith et al. (2018), and fusion energy science Gaffney et al. (2018); Vinko et al. (2020). Experimentally, the systems produced are typically studied via x-ray emission or inelastic scattering Toleikis et al. (2010) spectroscopy, while the theoretical description is provided by time-dependent atomic kinetics calculations Chung et al. (2017); Ciricosta et al. (2016a).
Because of the short pulse duration of FEL x-ray sources, on the order of 100 fs, the time evolution of the ionization and ionic charge state populations must be modelled directly via collisional-radiative rate equations, and local thermodynamic equilibrium (LTE) cannot be assumed – such models are thus called non-LTE. In contrast, these same timescales are typically long for the purposes of electron equilibration, and the free electron distribution is often modelled using a Maxwell-Boltzmann distribution function at some evolving (but instantly equilibrating) electron temperature and density. Being able to replace the time evolution of the full electron distribution function by a time-dependent temperature and density is convenient computationally as it leads to a substantial simplification in modelling complexity, but it is also carries significant physical implications. For one, the shape of the electron distribution determines multiple properties that govern the evolution of the atomic kinetics, including photoemission and collisional ionization rates; these can be altered significantly by non-thermal electron distributions.
The importance of properly accounting for electron equilibration dynamics was demonstrated in early ultrafast experiments at the FLASH FEL with XUV pulses a few tens of fs in duration. There, the electron relaxation time was comparable to the typical lifetime of radiative decay, leading to fluorescence emission from a seemingly cold (non-equilibrated) distribution, alongside a bremsstrahlung emission spectrum arising from the non-thermal electrons consistent with a 20-times higher effective temperature Vinko et al. (2010); Medvedev et al. (2011). The contribution from non-thermal electrons to the total emission spectrum is amplified in such experiments due to the combination of low energy XUV photoexcitation ( eV) and short pulse lengths. In contrast, similar effects are expected to be smaller at harder x-ray photon energies (1 keV) and for typical SASE pulse durations of 60-80 fs Vinko et al. (2012); Preston et al. (2017). Nevertheless, experimental access to both high x-ray intensities and ultra-short pulse lengths (10 fs) will undoubtedly lead to novel x-ray matter interaction regimes where non-thermal electron dynamics dominate the plasma evolution, even on x-ray emission and scattering timescales. This poses a challenge to established modelling approaches, but also provides new opportunities to study ultrafast electron equilibration in a hot-dense plasmas experimentally in situ.
The effect of non-thermal electrons in extreme conditions can be modelled in various ways, and efforts have been made in both x-ray Thomson scattering Milder et al. (2019) and x-ray emission spectroscopy investigations via non-thermal atomic kinetics de la Varga et al. (2013). A kinetic modelling addition to the particle-in-cell (PIC) simulation suite PICLS Sentoku and Kemp (2008); Sentoku et al. (2014) allowed for the inclusion of self-consistent radiation transport, photoionization and Auger electron transport, enabling a full kinetic simulation of XFEL-driven plasmas Royle et al. (2017). While PIC-based models are very capable in capturing the evolution of the electron distribution in time and in space, they lack the full predictive capability for detailed x-ray emission spectroscopy. In contrast, non-thermal collisional-radiative codes deploying a Vlasov-Boltzmann-Fokker-Planck equation to treat non-Maxwellian electrons Le et al. (2019) are able to treat a vast number of ionic charge configurations, essential for predicting x-ray spectroscopy, and allow us to investigate the effect of non-thermal distributions on spectroscopic experiments at FEL facilities.
In this context we present CCFLY, a non-thermal, non-LTE collisional radiative code, designed to model the time-dependent evolution of of both electron distributions and ion states interacting with intense x-ray fields on ultra-short timescales. The code follows on from the SCFLY suite Chung et al. (2017, 2007), itself based on the FLYCHK code Chung et al. (2005), which has proven successful in modelling and interpreting x-ray FEL experiments over the past decade Ciricosta et al. (2016a).
CCFLY allows us to model the plasma evolution by tracking arbitrary electron distribution functions during the atomic-kinetics calculation, foregoing the need to assume a Maxwellian distribution for the free electrons. More precisely, the Fokker-Planck numerical framework is used to describe the evolution of the electron distribution during the elastic electron-electron collisions and account for the electrons released or absorbed during inelastic electron-ion interactions. A self-consistent method is achieved to compute the density and energy balance of the free-electron population together with the ion population dynamics in a partially ionized warm dense plasma. Thus, it becomes possible for us to investigate ionization dynamics in both equilibrated and non-equilibrated systems at high density. We note here that while CCFLY tracks the time-evolution of the electron distribution and of the ionic charge states, it does not explicitly take into account the motion of ions. As is normally the case in atomic-kinetics modelling, there is no molecular dynamics, and the ions have no associated kinetic energy (and thus no temperature). This is adequate to describe the evolution of x-ray-irradiated samples on femtosecond timescales, but starts to breakdown as the ions are heated. However, ion heating takes place on longer picosecond timescales, normally dictated by electron-phonon coupling dynamics Cho et al. (2011); Vinko et al. (2020); Lee et al. (2021); Zhang et al. (2022).
In this article use CCFLY to study the effect of non-thermal electrons in Al and Mg plasmas irradiated by x-ray pulses with characteristics achievable at current x-ray FEL facilities. We study both normal isochoric heating experimental setups, where the x-ray photon energy is chosen to be higher than the ionization threshold of several ionic charge states, and resonant pumping setups where the x-ray pulse is tuned to a specific bound-bound atomic resonance. In both cases we find that non-thermal electron processes become important for short pulse durations and at high intensities. We quantify the regimes for which instantaneous thermalization assumptions hold, and provide a more rigorous justification for the modelling assumptions employed in previous x-ray spectroscopy experiments at the Linac Coherent Light Source (LCLS) Andreasson et al. (2011); Vinko et al. (2012); Ciricosta et al. (2012); Cho et al. (2012); Brown et al. (2014); Rackstraw et al. (2014); Vinko et al. (2015); Ciricosta et al. (2016b, a); Preston et al. (2017); Kasim et al. (2018); van den Berg et al. (2018); Vinko et al. (2020).
II Theory
The atomic kinetics modelling approach aims to track the evolution of the charge states in the plasma as it interacts with an intense x-ray pulse. The model assumes we have ions in various ionization states, interacting with a classical free electron gas (the ionized free electrons), the density of which is constrained by the requirement for the overall plasma to be charge-neutral. We use a zero dimensional collisional-radiative model, where the transitions between various ion states are driven by electron collisional processes, or are mediated by the photon field. The evolution of the system is described in terms of the changes in the relative population density of each ionic charge state (labeled by the index ), and of the free electron distribution . A collisionally-induced transition between two ion states given by some cross section is described in the model by a transition rate per unit volume given by
(1) |
where is the kinetic energy of the electron and the electron mass. By computing the transition rates between all the ion states, we can construct a transition matrix . For each available ion state , its population density can then be determined by the following set of coupled rate equations
(2) |
Because our system is globally neural, the total electron density can be calculated from the ion density and the ionization as:
(3) |
where we have summed across all ion states and multiplied them by their ionization . From the perspective of energy balance, the ion transitions can either add energy or draw energy from the free electron population. Thus, due to the requirement of energy conservation, the energy change in the free electron distribution is given by:
(4) |
where is the energy difference in between ionic charge states and involved in the transition . This expression sums up all the energies that are released into, or drawn from, the free electron distribution due to photo-absorption, Auger decay, and collisional interactions. This includes bound-bound radiative transitions, bound-bound collisional transitions, bound-free radiative ionisation, bound-free radiative recombination, bound-free collisional ionisation, bound-free collisional recombination, autoionization, and electron capture.
We will call the model thermal if we assume that the electron distribution function is Maxwellian, uniquely described by an electron density and a temperature . For situations where energy is entering the system, for example via x-ray photoabsorption, the temperature and density are instantaneously equilibrated at each time step. Then,
(5) |
In this case, Eqs. (3) and (4) completely determine the free-electron distribution in terms of the changing temperature and density of the ionized electrons.
In contrast, if a thermal electron distribution cannot be assumed, then the change in the electron distribution needs to be modelled explicitly for all electron energies. This can be accomplished using the Fokker-Plank equation:
(6) |
The change in the electron distribution as a function of time is governed by four terms on the right hand side of this equation. The first two terms, containing and , describe the elastic part of the evolution. The third term denotes all contributions of inelastic collisional processes within the electron distribution, where the free electrons either gain or lose energy to the ionic charge state. The last term denotes the collection of all source interactions that do not depend on the electron distribution. By capturing the effects on the electron distribution in an inelastic collision operator that considers the photo-ionisation and Auger decay, a self-consistent approach for the ion populations and electrons can be achieved. Then the ionic charge state populations, the electron distribution and its time-derivatives can be fed into an ordinary differential equation (ODE) solver. We refer to this situation as the non-thermal case. Note that the ion population is never assumed to be thermal, and the time-dependent modelling is purely non-LTE for both the thermal and non-thermal simulation approaches.
II.1 Atomic data
CCFLY supports two levels of detail in the atomic data used in the atomic kinetics modelling: superconfigurations and relativistic configurations. In the super-configurational approach, we simply count how many electrons are bound with each quantum number. States are denoted in the notation which means that there are electrons in the K-shell (n = 1), electrons in the L-shell (n = 2), electrons in the M-shell (n = 3) and electrons in the N-shell (n = 4). The advantage of using superconfigurations is that the dataset can be kept relatively small compared with a more detailed configurational dataset. In contrast, relativistic configurations provide access to more detailed atomic physics modelling, but at a substantially greater computational cost. A relativistic configurations dataset can be generated using the Flexible Atomic Code Gu (2008) Gu (2004), via the CFACDB library that provides access to the cFAC-generated databases.
For the study presented in this work we will use the superconfigurational dataset for our atomic kinetics modelling. In contrast, a more detailed atomic model using relativistic configurations calculated by the DHS code Chen et al. (1980) was used for the spectral synthesis where energy levels, statistical weights, and radiative data are applied to produce frequency-dependent emissivities and opacities. Population distributions of relativistic configurations belonging to a super-configuration (or to a super-shell) were determined by Boltzmann statistics as a function of electron temperature, as described in more detail in ref. Chung et al. (2005, 2016); Ciricosta et al. (2016a).
II.2 Continuum lowering
The effect of continuum lowering, where the binding energies of atomic states are lowered due to high electron and ion density, is accounted for in the atomic kinetic modelling via an ionisation potential depression (IPD) model. The CCFLY code supports two main IPD models: the modified Ecker-Kroll IPD model (EK) Preston et al. (2013) and the Stewart-Pyatt IPD model Stewart and Pyatt (1966). All simulations presented here were carried out with the EK model, to make the results directly comparable to previously published related work.
II.3 Elastic Electron Collisions
The elastic part of the evolution of the electron distribution is given by the first two terms in the Fokker Planck equation in Eq. II. The elastic collision terms can alter the shape of the distribution, while preserving particle density and total kinetic energy. Additionally, the distribution function should equilibrate to the correct Maxwellian distribution. In CCFLY, the electron distribution is initialised and evolved in energy space but in the Fokker-Planck operator, the electrons are evolved in velocity space.

First derived by Bobylev and Chuyanov Bobylev and Chuyanov (1976), and cast into a non-relativistic collision operator by Tzoufras et al. Tzoufras et al. (2011)Tzoufras et al. (2013), this formulation is given in integro-differential operator form as:
(7) |
where the function is given by:
(8) |
The function contains the time evolution of the electron distribution due to to electron-electron collisions that conserve particle number and kinetic energy. The is defined as:
(9) |
with the electron charge, and the Coulomb logarithm. The Coulomb logarithm for electron-electron interactions in hot plasmas is computed with a semi-empirical formula constructed by Huba Huba (2013):
(10) |
where the density and temperature have units of cm-3 and eV, respectively.
For large velocities, Eq. (8) can be evaluated through standard discretisation techniques such as the chained trapezoid method. For velocities that are small compared with the characteristic velocity of the distribution, the second and third term in Eq. (8) cancel up to . It is necessary to modify the numerical scheme to capture this, in order for the distribution to equilibrate to a Maxwellian. This can be done by writing as a Taylor series and . Setting , this results in:
(11) | ||||
Similarly, setting , we have
(12) | ||||
We show an example of how this elastic collisional operator relaxes a non-thermal electron distribution to a Maxwellian in Fig. 1. For simplicity, we assume here that only free electrons are involved in the collisional thermalization process, and all interactions with ions (e.g. collisional ionization and recombination) are ignored. We initialize the electrons to a thermal distribution at 40 eV, and add an additional non-thermal population of electrons with energies around 333 eV. The collisional operator equilibrated the distribution to a new overall temperature of 42 eV, the majority of which takes place on a timescale of a few femtoseconds.
III Non-thermal electrons in x-ray isochoric heating
The capabilities of the CCFLY code allow us to look into the detailed evolution of the electron distribution function in an isochoric heating experiment at an x-ray FEL. In particular, we are interested in exploring what role the non-thermal electrons play in the atomic kinetics of the system, and wish to assess the range of validity of the instantaneous electron thermalization assumption commonly used in the modelling and in the interpretation of experimental results.


We start by modelling the heating process of a thin Al foil irradiated by x-rays at a photon energy of 1590 eV. This energy was deliberately chosen to be above the Al K-edge (1560 eV) so that the dominant absorption mechanism is K-shell (1s) photoionization. The x-ray pulse is chosen to be a Gaussian in time with a FWHM duration of 40 fs, and has a spectral bandwidth of 0.4%, typical of a SASE FEL pulse. Four different intensities are compared, between Wcm-2 and Wcm-2, depositing an ever increasing amount of energy into the target. Experimentally, these intensities would be realized by focusing the x-rays to various spot sizes: focusing the beam to a spot of a few micron in diameter typically yields intensities in region of Wcm-2 and can be readily achieved on most facilities.
We show the evolution of the electron density in Fig. 2, where we compare the results of the thermal and non-thermal models side by side. The x-ray pulse peaks in intensity in the middle of the simulation window, at 40 fs. The thermal simulations simply show a Maxwellian distribution that gets increasingly hotter in time over the duration of the excitation pulse. More deposited energy leads to increasing ionization, and thus the free electron density also increases over time. The temperatures and densities are higher for higher irradiation intensities, since more energy is deposited into smaller volumes, as expected. The electron distribution function for the non-thermal simulation is visibly different, and can be divided into two components. The first is the bulk part of the free electron distribution, lying at lower energies. This assumes a near-Maxwellian shape because the electron thermalization time is short even on the 100 fs timescales used here (see Fig. 1 and related discussion). This bulk component is largely equilibrated, and grows in prominence as the x-rays are absorbed and the sample heats. At higher energies, however, we can see a large energy spread over which non-thermal electrons can be found. These are mostly created by the equilibrating hot electrons generated in the x-ray-matter interaction process, and subsequent atomic kinetics. Two main sources of quasi-mono-energetic hot electrons can be identified: a feature at due to the photoionization of the Al L-shell electrons (at 1520 eV), and a feature due to the Auger KLL process (at 1420 eV) that efficiently fills a K-shell core hole. Both features can be readily seen in the figure, but are clearest at the lowest intensities. They are harder to discern at high intensity because their relative populations, compared with the bulk free-electron distribution, decreases with intensity.
Figure 2 shows the density distribution on a logarithmic scale, and thus amplifies the non-thermal signature. But for the underlying physics, what matters most is the relative contribution of the quasi-thermalized component of the electron distribution in driving the atomic kinetics, compared with the non-thermal hot electron population. To examine this, we plot the populations for ionic charge states to in Fig. 3, in both the thermal and non-thermal case. We immediately note that the results are comparable: the charge state populations do not change dramatically under the instantaneous thermalization approximation, validating this common assumption. Nevertheless, we do see a lag in the time at which certain charge states peak in the population distribution between the thermal and non-thermal calculations. In particular, we note that lower charge states in the more detailed non-thermal calculations tend to appear later in time compared with the thermal prediction, whereas higher charge states appear earlier in time. This is primarily due to the mechanism through which these charge states are created. The lower charge states are more efficiently generated via collisional ionization processes driven by the bulk free electrons; these electrons heat up faster in the thermal calculation as all energy is placed into the thermal distribution instantaneously, whereas in the non-thermal case the bulk temperature rises more slowly as the energy needs to be transferred from the highly energetic photoionization and Auger electrons. In contrast, the higher charge states require substantially higher temperatures to be collisionally ionized, and are thus more efficiently created via collisional interactions with the non-thermal electrons. Even if there are fewer of them available, they have sufficient energy to play an important, and possibly dominant role, in the dynamics of collisional ionization.
In an experiment one cannot normally measure individual plasma charge state populations. However, these populations can be inferred via x-ray emission spectroscopy of Kα satellite lines, when the sample is excited by x-rays above the relevant ionization thresholds Vinko et al. (2012); Vinko (2015). For charge states photo-pumped below their ionization threshold, the appearance of emission lines can be instead indicative of collisional ionization processes, and thus used to measure collisional rates Vinko et al. (2015). For this measurement to be reliable, however, the bulk free electrons must be the dominant contributor to collisional ionization, and the contributions from non-thermal electrons should be small.
We investigate the contribution of non-thermal electrons to the total Kα emission spectrum in Fig. 4. The x-ray photon energy is set just above the K-edges of the first two charge states, but below the ionization threshold of higher charge states. Thus only the first two charge states can be revealed directly via Kα emission spectroscopy. While higher charge states are still present in the plasma, as illustrated in Fig. 3, their states cannot be photo-ionized and they therefore cannot relax emitting a Kα photon. Despite this, emission from higher charge states can still be observed, because of electron collisional ionization of the L-shell in the presence of a K-shell hole van den Berg et al. (2018). Non-thermal electrons affect the emission spectrum in two key ways. Firstly, because the equilibration of photoionization and Auger electrons is the main heating channel of the bulk free-electron distribution, non-thermal electrons help determine the equilibration timescales and drive the evolution of the charge state population towards LTE. These evolving populations ultimately determine the emission intensity of charge states pumped above their K-edge thresholds. The second mechanism applies to charge states that are pumped below their K-shell ionization threshold. These are primarily driven by collisional ionization of the L-shell electrons, which can be either due to the bulk electrons at some average temperature eV, or to the non-thermal electrons at a much higher energy (typically ). Higher electron energies are more efficient at L-shell ionization, and so a system with a large non-thermal electron population and relatively cold bulk distribution will lead to more collisional ionization, and consequently to more intense emission from higher charge states. We see this play out in the panels of Fig. 4. Because the x-ray pulse length is relatively long compared with the electron relaxation timescale, the thermal and non-thermal modelling predict very similar intensities for the first two emission lines. The charge state population is very close to LTE for intensities up to around 1018 Wcm-2, and the intensity of the first two lines is simply a reflection of the population of those charge states modulated by the temporal intensity of the x-ray pulse: the second charge state is most prominent at irradiation intensities of around 1016 Wcm-2, and decreases for lower (insufficient heating) and higher (too much heating) intensities. We note a slight difference in the intensity of the second line at the highest intensities: here the charge state is swept through very quickly, and the instantaneous thermalization assumption starts to breakdown as it underestimates the time needed to establish LTE in the plasma distribution. The emission intensity of higher lines (VI and above) shows a considerably stronger dependence on the x-ray intensity. For intensities up to several times 1016 Wcm-2 the thermal and non-thermal predictions match, indicating that the L-shell ionization process described above is indeed dominated by quasi-thermal bulk free electrons. Above this, the number of non-thermal electrons remains sufficiently high to more than double the effective ionization rate, indicating a stark departure from LTE collisional dynamics.

The results shown in Fig. 4 are for a single x-ray intensity, and are thus not directly comparable with experimental results that are measured over a broad range of x-ray intensities. However, the average experimental intensity (rather than the often quoted peak intensity) is indicative of the overall irradiation regime, and can be linked to the discussion above. Average intensities typically achieved at the LCLS FEL, such as those fielded for the work on collisional ionization presented in Refs. Vinko et al. (2015); van den Berg et al. (2018), is on the order of Wcm-2.
III.1 The effect of x-ray pulse length
We have seen that for a given x-ray pulse duration, the effect of focusing the x-ray beam to achieve increasing intensities is to increase the relative contribution of non-thermal electrons to the evolution of the plasma, and thus to the observed emission spectrum. We now turn our attention to the effect of pulse length. Changing the pulse length is another way to change the intensity on target, but it adds an additional layer of complexity because the pulse duration also gates the observation timescale of the system. This is primarily because we only observe x-ray emission in the Kα window when electrons can be ionized, and the overall temperatures reached in a typical microfocusing experiment ( eV) are insufficient to drive collisional K-shell ionization to any notable extent. Thus, emission is only observed during the x-ray pulse, and for a short time after depending on the femtosecond core-hole lifetime. By changing the duration of the x-ray pulse, while maintaining a constant pulse energy, we thus not only modify the x-ray intensity seen by the plasma, but also change the period over which we view the evolution of the charge state distribution. We show the simulated results for a Al plasma in Figs. 5 and 6. For these simulations we have assumed a pulse energy of 1 mJ, and a spotsize on-target of 12 m2, so that a pulse length of 80 fs corresponds to an x-ray intensity of about 1017 Wcm-2. These are typical values achievable at FEL facilities today Vinko et al. (2012). The pulse energy is then kept constant, and the pulse duration lowered down to 8 fs.

For the longer pulse lengths of 40 and 80 fs we see little difference in the computed emission spectrum for the thermal and non-thermal cases, and a lag in the creation of various charge states analogous with that described in the previous section. In these cases, the observation window is sufficiently long for the evolution of the change state populations to agree well with the assumption of instantaneous electron thermalization. The overall population of the L-shell states is collisionally driven by the equilibrated bulk free electron distribution, and the charge states evolve as if in near-LTE for each point in time. For pulse durations of 20 fs and less, these assumptions start to breakdown, and the calculated spectrum assuming instantaneous electron thermalization underestimates the emission intensity.


III.2 Non-thermal electrons in resonant spectroscopy

Recently, van den Berg et al. proposed driving a resonant transition in Mg to measure collisional ionization rates in dense plasmas van den Berg et al. (2018). The principle is straightforward: an x-ray pulse form an FEL is tuned to resonantly drive a specific inner-shell Kα transition, and a spectroscopic measurement is performed to look for emission from neighbouring Kα emission lines. Because of the narrow bandwidth of the FEL beam, the x-ray pulse can only excite a single Kα line, and so the observation of multiple lines is indicative of L-shell collisional processes taking place during the lifetime of a K-shell core hole. This method thus allows for a measurement of an ultra-fast collisional ionization rate. However, the real objective is not so much to infer a rate, but rather to measure a collisional cross section. This can be accomplished using the expression in Eq. (1), under the caveat that the collisions are predominantly due to an interaction of the ion with the bulk free-electron distribution at some known temperature and density. But if the electron distribution is non-thermal, then it is easy to see that the cross section cannot uniquely be extracted from the integral expression of Eq. (1). Estimating the effect of non-thermal electrons is then important to guide the reliable interpretation of such experiments, especially since we have shown that the contribution of non-thermal electrons can be significant with increasing x-ray intensity or decreasing x-ray pulse length.
For this investigation we set out simulation parameters to match the conditions fo Ref. van den Berg et al. (2018): an optically thin Mg plasma is pumped at a photon energy of 1304 eV, corresponding to the resonance energy of the transition energy in a 7+ Mg ion. At the nominal pulse duration of 60 fs the intensity on target is Wcm-2. We show the evolution of the ionic populations as a function of time, using the thermal and a non-thermal models, in Fig. 7. The related emission spectra are shown in Fig. 8. Here we see that the assumption of instantaneous thermalization holds well for nominal pulse durations, above 30 fs, providing a more rigorous justification of the assumptions made in the work on collisional ionization rates. For shorter pulse durations, however, the contribution of non-thermal electrons starts to become significant, and even dominant for pulse lengths below 10 fs. At these short pulse lengths, modelling the full evolution of the electron distribution function becomes essential to produce atomic kinetics simulations with predictive capability.
IV Conclusions
We have presented the atomic kinetics simulation package CCFLY, which can model the non-thermal (electron) and non-LTE (ion) evolution of matter interacting with intense x-ray radiation on ultra-short timescales. We have used these capabilities to investigate the evolution of solid-density systems, heated isochorically to high temperatures via femtosecond FEL pulses, and assess the role of non-thermal electron distributions. We find that in mid-Z systems such as Al and Mg, non-thermal electrons become important in determining the short-time plasma kinetics and electron dynamics for intensities above Wcm-2, and for pulse lengths below 10 fs.
To date, the use of Fokker-Planck models within collisional-radiative simulations to interpret experiments at free-electron lasers has been limited, largely due to the substantial computational cost the method requires. Whilst we have shown that this is often justifiable for experiments using a typical SASE FEL pulse on the timescales of order 100 fs, ongoing developments enabling the routine generation of energetic few-fs and sub-fs pulses, and of ultra-high intensities via narrow bandwidth self-seeding and x-ray nanofocusing, are increasingly accessing x-ray regimes where the direct modelling of non-thermal electron relaxation will be essential to retain predictive power in the modelling of extreme x-ray matter interactions.
V Acknowledgments
S.R., M.F.K., J.S.W. and S.M.V. acknowledge support from the UK EPSRC grant EP/P015794/1 and the Royal Society. S.M.V. is a Royal Society University Research Fellow.
References
- Vinko et al. (2012) S. M. Vinko et al., Nature 482, 59 (2012).
- Fletcher et al. (2015) L. B. Fletcher et al., Nature Photonics 9, 274 (2015).
- Levy et al. (2015) A. Levy et al., Physics of Plasmas 22, 030703 (2015).
- Ciricosta et al. (2016a) O. Ciricosta et al., Physics of Plasmas 23, 022707 (2016a).
- Ciricosta et al. (2012) O. Ciricosta et al., Phys. Rev. Lett. 109, 065002 (2012).
- Vinko et al. (2014) S. M. Vinko, O. Ciricosta, and J. S. Wark, Nat. Commun. 5, 3533 (2014).
- Eggert et al. (2010) J. H. Eggert et al., Nature Physics 6, 40 (2010).
- Smith et al. (2018) R. F. Smith et al., Nature Astronomy 2, 452 (2018).
- Gaffney et al. (2018) J. Gaffney et al., High Energy Density Physics 28, 7 (2018).
- Vinko et al. (2020) S. M. Vinko et al., Phys. Rev. Lett. 124, 225002 (2020).
- Toleikis et al. (2010) S. Toleikis et al., Journal of Physics B: Atomic, Molecular and Optical Physics 43, 194017 (2010).
- Chung et al. (2017) H.-K. Chung et al., AIP Conference Proceedings 1811, 020001 (2017).
- Vinko et al. (2010) S. M. Vinko et al., Phys. Rev. Lett. 104, 225001 (2010).
- Medvedev et al. (2011) N. Medvedev, U. Zastrau, E. Förster, D. O. Gericke, and B. Rethfeld, Phys. Rev. Lett. 107, 165003 (2011).
- Preston et al. (2017) T. R. Preston et al., Phys. Rev. Lett. 119, 085001 (2017).
- Milder et al. (2019) A. L. Milder, S. T. Ivancic, J. P. Palastro, and D. H. Froula, Physics of Plasmas 26, 022711 (2019).
- de la Varga et al. (2013) A. G. de la Varga, P. Velarde, F. de Gaufridy, D. Portillo, M. Cotelo, A. Barbas, A. González, and P. Zeitoun, High Energy Density Physics 9, 542 (2013).
- Sentoku and Kemp (2008) Y. Sentoku and A. Kemp, Journal of Computational Physics 227, 6846 (2008).
- Sentoku et al. (2014) Y. Sentoku, I. Paraschiv, R. Royle, R. C. Mancini, and T. Johzaki, Phys. Rev. E 90, 051102 (2014).
- Royle et al. (2017) R. Royle, Y. Sentoku, R. C. Mancini, I. Paraschiv, and T. Johzaki, Phys. Rev. E 95, 063203 (2017).
- Le et al. (2019) H. P. Le, M. Sherlock, and H. A. Scott, Phys. Rev. E 100, 013202 (2019).
- Chung et al. (2007) H.-K. Chung, M. Chen, and R. Lee, High Energy Density Physics 3, 57 (2007).
- Chung et al. (2005) H.-K. Chung, M. Chen, W. Morgan, Y. Ralchenko, and R. Lee, High Energy Density Physics 1, 3 (2005).
- Cho et al. (2011) B. I. Cho et al., Phys. Rev. Lett. 106, 167601 (2011).
- Lee et al. (2021) J.-W. Lee et al., Phys. Rev. Lett. 127, 175003 (2021).
- Zhang et al. (2022) J. Zhang, R. Qin, W. Zhu, and J. Vorberger, Materials 15 (2022), 10.3390/ma15051902.
- Andreasson et al. (2011) J. Andreasson et al., Phys. Rev. E 83, 016403 (2011).
- Cho et al. (2012) B. I. Cho et al., Phys. Rev. Lett. 109, 245003 (2012).
- Brown et al. (2014) C. R. D. Brown et al., Scientific Reports 4, 5214 (2014).
- Rackstraw et al. (2014) D. Rackstraw et al., High Energy Density Physics 11, 59 (2014).
- Vinko et al. (2015) S. M. Vinko et al., Nat Commun 6, 6397 (2015).
- Ciricosta et al. (2016b) O. Ciricosta et al., Nat. Commun. 7, 11713 (2016b).
- Kasim et al. (2018) M. F. Kasim, J. S. Wark, and S. M. Vinko, Scientific Reports 8, 6276 (2018).
- van den Berg et al. (2018) Q. Y. van den Berg et al., Phys. Rev. Lett. 120, 055002 (2018).
- Gu (2008) M. F. Gu, Canadian Journal of Physics 86, 675 (2008).
- Gu (2004) M. F. Gu, AIP Conference Proceedings 730, 127 (2004).
- Chen et al. (1980) M. H. Chen, B. Crasemann, and H. Mark, Phys. Rev. A 21, 436 (1980).
- Chung et al. (2016) H.-K. Chung, S. B. Hansen, and H. A. Scott, “Generalized collisional radiative model using screened hydrogenic levels,” in Modern Methods in Collisional-Radiative Modeling of Plasmas, edited by Y. Ralchenko (Springer International Publishing, Cham, 2016) pp. 51–79.
- Preston et al. (2013) T. R. Preston, S. M. Vinko, O. Ciricosta, H.-K. Chung, R. W. Lee, and J. S. Wark, High Energy Density Physics 9, 258 (2013).
- Stewart and Pyatt (1966) J. C. Stewart and J. Pyatt, Kedar D., Astrophys. J. 144, 1203 (1966).
- Bobylev and Chuyanov (1976) A. Bobylev and V. Chuyanov, USSR Computational Mathematics and Mathematical Physics 16, 121 (1976).
- Tzoufras et al. (2011) M. Tzoufras, A. R. Bell, P. A. Norreys, and F. S. Tsung, J. Comput. Phys. 230, 6475 (2011).
- Tzoufras et al. (2013) M. Tzoufras, A. Tableman, F. S. Tsung, W. B. Mori, and A. R. Bell, Physics of Plasmas 20, 056303 (2013).
- Huba (2013) J. D. Huba, Plasma Physics (Naval Research Laboratory, Washington, DC, 2013) pp. 1–71.
- Vinko (2015) S. M. Vinko, Journal of Plasma Physics 81, 365810501 (2015).