Next Article in Journal
Hepatocyte Proteome Alterations Induced by Individual and Combinations of Common Free Fatty Acids
Next Article in Special Issue
Electron Attachment to 5-Fluorouracil: The Role of Hydrogen Fluoride in Dissociation Chemistry
Previous Article in Journal
Mitochondrial Protein Homeostasis and Cardiomyopathy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulating the Feasibility of Using Liquid Micro-Jets for Determining Electron–Liquid Scattering Cross-Sections

1
College of Science & Engineering, James Cook University, Townsville, QLD 4811, Australia
2
Department of Medical Physics, Townsville University Hospital, Townsville, QLD 4814, Australia
3
Research School of Physics and Engineering, Australian National University, Canberra, ACT 0200, Australia
4
Synchronous Technologies PTE LTD, 6 Raffles Quay, #11-07, Singapore 048580, Singapore
5
College of Science & Engineering, Flinders University, Bedford Park, SA 5042, Australia
6
Institute of Actuarial Science and Data Analytics, Faculty of Business and Management, UCSI University, Kuala Lumpur 56000, Malaysia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(6), 3354; https://doi.org/10.3390/ijms23063354
Submission received: 1 March 2022 / Accepted: 17 March 2022 / Published: 20 March 2022

Abstract

:
The extraction of electron–liquid phase cross-sections (surface and bulk) is proposed through the measurement of (differential) energy loss spectra for electrons scattered from a liquid micro-jet. The signature physical elements of the scattering processes on the energy loss spectra are highlighted using a Monte Carlo simulation technique, originally developed for simulating electron transport in liquids. Machine learning techniques are applied to the simulated electron energy loss spectra, to invert the data and extract the cross-sections. The extraction of the elastic cross-section for neon was determined within 9% accuracy over the energy range 1–100 eV. The extension toward the simultaneous determination of elastic and ionisation cross-sections resulted in a decrease in accuracy, now to within 18% accuracy for elastic scattering and 1% for ionisation. Additional methods are explored to enhance the accuracy of the simultaneous extraction of liquid phase cross-sections.

1. Introduction

The interaction of a low-temperature plasma with liquids is fundamental for numerous new and emerging technologies, finding applications in important domains, including environmental remediation, and the synthesis of nanomaterials and medicine [1,2,3,4,5,6,7,8,9]. The goal of high-level optimisation, of the efficacy and selectivity, of these and future generation plasma–liquid applications depends on, among other things, a detailed understanding of the underlying fundamental nanoscale physics and associated predictive modelling, underpinned by accurate and complete transport theory and cross-sections. The key driver to these applications is the role of electrons (and other radical species) at the plasma–liquid interface, but despite their central role, electron-induced transport and processes at the interface are not well understood [10].
Developing our understanding of electron transport into and within liquids is critical for enhancing the predictive power of plasma–liquid models. Fundamental to transport theory, which governs the motion of electrons into and within such environments, are complete and accurate sets of electron impact cross-sections. While there is a wealth of knowledge of electron impact cross-sections in the gas phase (e.g., databases such as LXCat [11,12]), the same is not true for the liquid environments. The scattering and transport theory of pre-solvated electrons in non-polar liquids is reaching some level of maturity (see the ab initio treatment by the authors [13,14]), however this is not the case for polar liquids. An existing Monte Carlo simulation of low-energy electrons in liquid water [15] ignored many physical processes, that are relevant to developing appropriate scattering theory. Furthermore, current models use cross-sections calculated in the gas phase, or through electron reflection measurements from amorphous ice [16]. The absence of liquid phase data represents a large knowledge gap in the literature, and experiments proposed in this manuscript aim to address, at least in part, this knowledge gap.
Here, we propose a new experimental technique that extends gas phase beam experiment methodologies [17] to include electron scattering from liquid surfaces, through the use of liquid micro-jets (LµJs). A constant replenishment of liquid, a small surface area, and a laminar flow enable the scattering of particles from a smooth and stable liquid surface, whilst operating in a vacuum environment. Until now, LµJ investigations have focused primarily on photon and photoelectron scattering [18,19,20,21,22,23,24,25,26,27,28,29,30]. In this work, we introduce electron scattering from LµJs and investigate the feasibility of developing effective electron–liquid cross-sections from the measured electron energy loss spectra (EELS). Unlike gas phase experiments, which are designed to ensure single scattering processes only, the proposed experiment is necessarily multi-scattering, with electrons being scattered from the surface and bulk atoms/molecules. The connection of this with traditional (multiple collision) swarm experiments [31,32,33,34] is thus clear.
Swarm experiments have proven to be critical in the development of accurate cross-section sets [17]. In electron swarm experiments, electrons are driven through a gaseous (or liquid) medium by an applied electric field, and macroscopic descriptors such as current, drift velocity, and diffusion coefficients are determined. Simulation techniques then evaluate the accuracy and self-consistency of the cross-section sets through comparisons of the transport coefficients. The same techniques can be used, in principle, to iteratively improve or develop cross-section sets [35,36,37,38,39,40], although the question of degeneracy in the cross-section set remains open, and other means are usually required to minimise this. Recent studies have shown promising applications of neural networks (NN), toward improving cross-section sets through swarm transport data [41,42,43,44,45,46]. Utilising a similar NN methodology, along with energy loss spectra from a LµJ experiment, we propose the generation of cross-section sets for electrons within a liquid environment. To facilitate this, a benchmarked Monte Carlo simulation technique, for simulating electron transport in liquid environments [47], has been developed and implemented for this study, with extensions to calculate EELS arising from the scattering of electrons from an LµJ.
The paper is organised as follows. In Section 2, we provide a summary of LµJs and their implementation in the proposed experiment. In Section 3, we outline the Monte Carlo (MC) simulation method utilised here and discuss the characterisation of bulk liquid and interfacial effects. The characteristic signatures of the various multiple scattering elements on the EELS are also discussed. In Section 4, we use the developed Monte Carlo software to train a neural network to predict neon’s elastic and ionisation cross-sections as a proof-of-concept. Finally, in Section 5, we provide some concluding premarks.

2. The Proposed Electron–Liquid Micro-Jet Scattering Experiment

2.1. Liquid Micro-Jets

Liquid micro-jets were first developed to overcome the limitations of probing liquids under vacuum conditions, which are often required by scattering experiments [48]. Due to the high vapour pressures found in liquids, additional treatment is required to ensure a stable surface within vacuum conditions, and hence the developments of the liquid jets. In our proposed experiment, a replenishing and sufficiently thin liquid source, exhibiting laminar flow, facilitates a smooth and stable scattering surface within a vacuum environment [49].
Siegbahn and Siegbahn [48] first utilised a replenishing 0.2 mm diameter liquid jet, and scattered electrons from that source to measure binding energies in formamide ( HOCN H 2 ). In vacuo experiments were initially limited to low vapour pressure liquids, to avoid immediate freezing and/or evaporation. For high vapour pressure liquids such as water, a further reduction of exposed surface area is required to maintain a free vacuum surface [20]. To this end, significantly smaller (∼10 µ m ) liquid micro-jets were first developed to measure the velocity distributions of vapour molecules, while ensuring inter-molecular collisions were minimised [49]. Extending this, a focusing gas surrounding the nozzle was implemented to increase the length and decrease the diameter of the jet substantially [50,51,52]. Flat LµJs were also developed from the collision of two micro-jets, which results in a sheet-like scattering surface [53]. Additionally, cryogenic jet systems enable the use of super-cooled liquids with jet diameters reaching 1 µ m [54,55]. It is clear that a wide array of configurations are achievable, and can be tailored towards the desired experiment.
Until now, applications of the aforementioned LµJ designs include an array of photoelectron spectroscopy experiments, to measure properties such as binding energies and hydrogen bonding [18,19,20,21,22,23,24,25,26,27,28,29,30]. Other LµJ experiments include mass spectrometry [56], enhanced X-ray production [57,58,59,60], X-ray emission/absorption, to probe electronic structure and molecule orientation [61,62,63], to study the molecular dynamics of evaporation [64], electrokinetic power generation [65,66], and drug delivery alternatives [67]. For further general reading, one is directed to the summary papers cited here [23,68,69].
Scattering simulations, such as the MC code developed at James Cook University, rely on the characterisation of the LµJ features. Liquid dynamics studies of LµJs found jet length stability on the order of millimetres. In one study [70], source widths of 10, 20 and 50 µ m , along with velocities ranging between 25 and 150 ms 1 , were trialled for each diameter. A limited laminar flow was achieved with jet lengths between 1 and 15 mm . After emergence from a circular nozzle, the jet undergoes contraction to a final cross-sectional size due to surface tension forces [71]. Typical contraction values were found to range between 90 % and 60 % [70] of the original source width, after which capillary forces break the jet into droplets through a tendency to try and reduce surface energy [71].
Within the liquid, charge build-up, and the existence of a streaming potential, due to electrokinetic charging, can hinder the accurate simulation of charged particle transport. A streaming potential is produced when a pressurised liquid is forced through the nozzle, which disrupts the electric double layer created between the liquid and the inner wall of the nozzle [72,73]. For liquid water, the magnitude of this streaming potential exceeds 60 V depending on the jet velocity and diameter [19,73]. The addition of an electrolyte has since been shown to minimise this effect [73]. Additionally, the production of electrons within an insulating medium such as water will produce a current and subsequent surface potential [68]. The magnitude of this effect will depend on the application and relevant time scales.

2.2. Crossed-Beam Electron Scattering from Liquid Micro-Jets

In our proposed experimental work, we will extend existing crossed-beam methodologies from gas to liquid environments, utilising liquid micro-jets to measure differential electron energy loss spectra (EELS), as shown schematically in Figure 1. The liquid would be propelled into a vacuum environment through a 15 µ m circular nozzle to produce a continuous, stable, and smooth scattering surface. In this work, we assume that the jet diameter is the same as the nozzle diameter. As noted previously, this jet remains stable for several millimetres before decaying into droplets. Evaporation at the surface and the tapering of the jet’s diameter are neglected in the current work, as the characterisation of such parameters was not available, and their effects on the model are expected to be minimal. Electrons, sourced from the thermionic emission from a tungsten filament, and collimated and transported by a series of DC-potential elements, are then propelled through the vacuum normal to the jet’s direction, to undergo multiple surface and bulk scattering events before escaping (and being collected by one or more detectors) or solvating within the liquid. For the foreshadowed work, an electron source with an energy full width at half maximum (FWHM) of 0.5 eV will be used [74].
Once an electron escapes the jet, a scattered electron detector will resolve both its energy, ε , and scattering angle, χ , with a resolution of 0.8 eV (FWHM) and 1 , respectively [74]. Through the detection of both the energy and angular dependencies, EELS provides a high degree of insight into the underlying electron scattering cross-sections and their associated dynamics. This combination of the electron gun and scattered electron detector provides a combined energy resolution of 0.9 eV . We note that, with crossed-beam experiments, some restrictions will inevitably be placed on the available detection angles. This is simply due to the physical size of the electron source, the electron detector(s), and LµJ. Additionally, while an angle-resolved EELS provides sufficient information for the derivation of differential cross-sections within the gas phase, multiple scattering effects in liquids present significant challenges when deconstructing EELS into complete (differential) cross-section sets due to the high degrees of convolution arising. The extent of this convolution is investigated, in part, in Appendix B.

2.3. Extraction of Cross-Section Sets from EELS

The derivation of cross-sections from a reflected EELS measurement is akin to the ‘Inverse Milne Problem’, which consists of deriving the reflected EELS from a ‘half-plane’ medium, of finite or infinite depth from z = 0 to z = L (where L is the half-plane depth), and extends to infinity in all other directions [75]. Outside the half-plane, a vacuum exists from which particles are propelled into the half-plane. The presence of high multiple-scattering at large L, relative to the scattering length, results in a complex combinatorial problem with a restricted set of observables, thanks to a reliance on external measurements [75].
The vast majority of treatments involve energy-independent radiative transfer [75,76] and neutron scattering [77,78]. The inclusion of energy dependencies was conceptually treated [79], with strict delta function requirements placed on both the source and detector, in order to derive total and differential cross-sections. In other studies, sufficiently thin (10–100 nm ) solid films were used to derive surface and bulk inelastic mean free paths [80,81,82], assuming a knowledge of the elastic cross-section. Likewise, cross-sections for amorphous ice films were derived under a two-stream assumption, along with a variation of film thickness between 0 and 5.456 nm [16,83,84]. While the extension of these techniques towards LµJs may be feasible, other indirect approaches, such as Monte Carlo simulations, together with machine learning, could simplify the methodology significantly.
In multiple studies [41,42,43,45,46], machine learning was applied to the inversion of transport coefficients found through swarm experiments. While swarm transport coefficient measurements differ in nature to EELS, they each represent an ill-posed inversion problem associated with high degrees of multiple scattering. An extension of the current methodology towards liquid phase EELS requires a sufficiently large training dataset, which is facilitated through a highly configurable non-equilibrium Monte Carlo simulation of electron transport developed for this project. In this study, we thus seek an implicit cross-section extraction technique utilising a neural network trained through Monte Carlo simulations of an LµJ. In what follows, we outline the Monte Carlo simulation employed, and its application to the proposed LµJ measurements, before applying machine learning to this ill-posed inverse problem of determining cross-sections from LµJ-measured EELS.

3. Simulation of Electron Transport through Liquid Micro-Jets

Machine learning requires a substantial volume of training data to ensure a robust fitting process. To generate these data efficiently, a Monte Carlo simulation method was developed. Monte Carlo (MC) simulations provide flexible environments in which the manipulation of spatial and temporal parameters is both efficient and straightforward. LµJs inherently require specific and precise spatial variation in terms of the neutral density, and hence, an MC simulation is well suited for the task. The simulation technique is discussed in detail in [47], hence, here, we detail its application to LµJs, while a section of benchmark results is presented in Appendix A.

3.1. Liquid Dynamics

Relatively low-density gaseous environments involve an inherent assumption that each scattering event involves instantaneous localised interactions. In liquid environments, however, the de Broglie wavelength of the electron can be smaller than, or of the same order as, the inter-particle spacing and the mean free path, which necessitates modified coherent collision dynamics. To account for coherent scattering, a new method was recently developed which incorporated a structure factor modified cross-section [47]. Through this, the macroscopic effects of coherent scattering are realised, while retaining the efficiency of effective single particle scattering events. Therefore, the relative simplicity of modelling single particle-particle collisions is maintained, while replicating the macroscopic effects of multiple scattering. Through an integral structure factor [47,85], elastic collisions are split into three collisional processes, which are weighted by an angle-integrated structure factor, Γ ( ε ) , produced through either theoretical or independent experimental techniques.
As a result of these microscopic processes, electron scattering through liquids is accurately simulated, given an appropriate structure factor. For real systems, experimentally measured structure factors are available through X-ray or neutron scattering [86,87,88]. While its functional form will vary, the methodology remains the same. Previously [47], it was shown that these processes match the required momentum and energy transfers derived using Boltzmann’s equation [47]. While isotropic scattering was assumed, the proof of principle was extended by the present authors to show its validity for the use of differential cross-sections. In the current work, we investigate structure effects upon measured spectra, while its implementation in our machine learning model is left for future study.

3.2. Interfacial Dynamics

In addition to liquid dynamics, interfacial effects must be considered in the simulation of LµJs. A recent study investigated density-dependent inter-facial effects on electron transport [89]. Due to a change in de-localised electron energy across the interface, a density-dependent ( n 0 ), and hence spatially-dependent potential V 0 ( n 0 ( r ) ) , exists, which produces an effective electric field E 0 ( r ˜ ) = V 0 ( r ) . A spatial density dependency was thus implemented within the MC simulation, along with the resulting potential, and hence electric field. In that study [89], it was found that a discrete change in density, and hence potential, did not significantly change the path and associated transport dynamics of such particles when compared to a realistic functional form, which in addition to computational considerations, motivated the use of a step function form for the change in density in the current work. Additionally, any variation in the potential resulting from the density change is left for future extensions of this work.

3.3. Experimental Parameters

Under realistic experimental conditions, detector resolutions and uncertainties will impact the prominence of correlations between EELS and their corresponding cross-sections. A critical factor for extracting important information is the experimental energy resolution of both the electron source and detector. In the models discussed above, an approximation to the energy spread of the electron source was incorporated by sampling a Gaussian distribution of FWHM of 0.5 eV . In Figure 2, an experimentally feasible combined EELS energy resolution of 0.9 eV [74], which is incorporated into the simulations in what follows in Section 4, alongside a hypothetical resolution of 0.1 eV , are shown for comparison. Distinct information loss occurs in the 0.9 eV energy-resolved spectra, for inelastic peaks with similar threshold energies, when compared to the 0.1 eV spectra. This clearly indicates that the foreshadowed experimental measurements should be conducted with as narrow an energy resolution as possible.

3.4. Liquid Micro-Jet Parameters

To develop training data, simulated electrons were scattered through a cylindrical liquid micro-jet using the benchmarked MC simulation. Each electron, with initial energy ε 0 , was fired along the z-axis towards a 15 µ m jet, which was directed along the x-axis. It was assumed that the motion of the jet (∼100 ms 1 ) has a negligible effect on the measured EELS, as deflection angles were measured within the z-y plane and integrated over the x-axis.
Scattered electrons are detected once they reach a radial distance sufficiently far from the jet, such that the angle measured was relative to the jet’s centre in the z-y plane. The EELS was then integrated over the x-axis, such that there existed a χ , ε dependence, where ε is the electron’s final energy, before integrating over χ . Each electron was simulated until it was detected, or until it lost sufficient energy, such that the probability of escape was negligible.
As a proof-of-concept for this work, we assume that the jet is gaseous neon at a liquid density of 2.13 × 10 28 m 3 , and at a temperature of 0 K . Based on the proposed experiment [74], the energy variance of the initial electrons was sampled from a Gaussian with an FWHM of 0.5 eV , centred around the initial energy ε 0 , while its spatial extent was assumed to be a delta function directed along the surface normal. This serves as a first approximation to the experimental work, where there will be an electron beam, with a finite diameter of 1 mm , incident on the LµJ.

4. Determining Cross-Sections from Electron Energy Loss Spectra Using Machine Learning

Through a sensitivity analysis, which is presented in Appendix B, we show that there exists a highly complex and degenerate correlation between the energy loss spectra and their underlying cross-sections. It was found that, while absolute magnitudes had little effect on the EELS, the effects of both relative magnitudes and energy dependencies of the cross-section were relatively significant. Deriving cross-section sets directly from EELS, without prior knowledge, is therefore a formidable task.
In this section, we apply machine learning to this ill-posed inverse problem of determining cross-sections from LµJ-measured EELS. As a proof-of-concept, we initially consider the task of determining electron–Ne cross-sections from EELS that are calculated using our Monte Carlo simulation. In what follows, we assume gas phase isotropic scattering, due to the present limited availability of liquid phase electron scattering and differential cross-sections, the inclusion of which remains a critical step towards the determination of liquid phase electron cross-sections. We note that, given appropriate liquid phase test data, the methodology that we propose here could be used to produce effective liquid phase electron cross-sections.

4.1. Machine Learning Methodology

To obtain a solution to the “inverse EELS problem”, we apply the same general machine learning approach as proposed by Stokes et al. [41,42,43,44], for the analysis of electron swarm transport data, and utilise an artificial neural network of the form:
y x = A 4 mish A 3 mish A 2 mish A 1 x ,
where A n x W n x + b n are affine mappings defined by dense weight matrices W n and bias vectors b n , and mish x = x tanh ln 1 + e x is a nonlinear activation function [91] that is applied element-wise. The output vector, y , contains each cross-section of interest:
y = σ 1 ε σ 2 ε ,
all of which are a function of energy, ε , which becomes an input to the neural network alongside the available EELS data:
x = ε i n e I 1 I 2 ,
where I 1 , I 2 , are the electron intensities accumulated in each of the considered EELS energy bins. Note that we apply suitable logarithmic transformations to ensure that all inputs and outputs of the network are dimensionless and lie within 1 , 1 . In what immediately follows, we specify that each bias vector contains 64 parameters, with the exception of b 4 , of which the size must match the number of cross-sections in y . The weight matrices are sized accordingly.
In order to train the neural network, Equation (1), we require an appropriate set of example solutions to the inverse EELS problem. To ensure cross-sections provided by the network are physically plausible, we train on cross-sections from the LXCat project [11,12]. Specifically, we train with cross-sections of the form [42]:
σ ε = σ 1 1 r ε + ε 1 ε 1 1 r ε 2 r σ 2 r ε + ε 2 ε 1 1 r ε 2 r ,
where σ 1 ε and σ 2 ε are a random pair of LXCat electron scattering cross-sections, sampled from the available targets, of a given type (e.g., excitation, ionisation, etc.), r is a pseudo-random number uniformly distributed between 0 and 1, and ε 1 and ε 2 are their respective threshold energies. Once suitable cross-sections are found for the training, corresponding LµJ EELS can be determined using our Monte Carlo simulation. In total, we consider 10,000 such training exemplars.
We implement and train the neural network, Equation (1), using the Flux.jl machine learning framework [92]. The network is initialised such that its biases are zero and its weights are uniform random numbers, as described by Glorot and Bengio [93]. Training is performed using the AdaBelief optimiser [94] with Nesterov momentum [95,96], step size α = 10 3 , exponential decay rates β 1 = 0.9 and β 2 = 0.999 , and the small parameter ϵ = 10 8 . At each iteration, the optimiser is provided with a different batch of 4096 training examples, where each batch consists of 32 training cross-sections each evaluated at 128 random energies of the form ε = 10 s , where s 0 , 2 is sampled from a continuous uniform distribution. For each batch, the optimiser adjusts the neural network weights and biases with the aim of further minimising the mean absolute error in solving the inverse EELS problem for that batch. Training is continued for 250,000 iterations, providing an equal number of potential solutions to the inverse problem. We then select every 10 for the last 100,000 cross-sections, and the quality of each of these solutions is subsequently assessed by simulating their corresponding EELS and comparing those to Ne’s EELS to find the ‘best’ regression.

4.2. Cross-Section Regression Given the EELS

We now implement and train neural networks of the form of Equation (1), to determine a selection of Ne’s cross-sections given the corresponding EELS, while assuming full knowledge of Ne’s remaining cross-sections. In total, we use eleven EELS, corresponding to initial electron energies of 100, 79.43, 63.1, 50.12, 39.81, 31.62, 25.12, 19.95, 15.85, 12.59 , and 10.0 eV . In each spectra, a combined energy resolution of 0.9 eV was assumed based on the detector performance that we currently find at Flinders University. Additionally, for computational considerations, each electron was simulated until 90 % of its energy was lost. For the initial energies considered, no electron was simulated below 1 eV , and thus we restrict the prediction domain to 1 eV , 100 eV .
We first determine only Ne’s elastic momentum transfer cross-section (MTCS) from the considered EELS. Figure 3a shows a reasonable level of agreement between the elastic MTCS for Ne [90] and those found by the neural network, with the best regression seen to be accurate to within 9 % . Figure 3b shows that, despite the large range in cross-sections for the 100 best regressions depicted in Figure 3a, the corresponding range of the EELS is much smaller, and in good agreement with those for Ne. This highlights the nonuniqueness of the inverse EELS problem, and suggests that there is not much room for improvement to the best cross-section fit plotted in Figure 3a unless; additional EELS were included, the energy resolution of the spectra was improved, or additional information about the unknown cross-sections was provided as input to the network.
We follow the above by now simultaneously determining both the elastic MTCS and ionisation cross-sections for Ne, using the same set of EELS. In this model, the bias vector was increased to 128, while 300,000 iterations were performed. As expected, the uncertainty in the elastic MTCS has increased here, with a larger elastic MTCS envelope found in Figure 4a, with the best regression seen to be accurate to within 18 % . In contrast, the corresponding ionisation cross-section envelope is particularly small, with the best regression being accurate here to within 1 % . We attribute this high accuracy to the prominent ionisation “shoulder” present in six out of eleven of the EELS. Figure 4b shows that the corresponding EELS are now close in line with those for Ne, despite having a range that is slightly larger than their counterparts in Figure 3b.
Using the MC simulation, every combination of the best 100 elastic and ionisation cross-section pairs was compared through their resulting EELS to find a further improved cross-section set. While the maximum error in the elastic MTCS fit increased to 20 % , the average error decreased resulting in a closer fit. For ionisation, the maximum error decreased to 0.2 % , which further emphasises the remarkable predictive capability of this model for ionisation cross-sections.
The extension of this methodology to include the prediction of excitation cross-sections, or more specifically, the simultaneous prediction of three or more cross-sections, results in a substantial decrease in accuracy within the current methodology. Thus, the simultaneous prediction of elastic, ionisation, and excitation cross-sections is left for future iterations of this method, which should include advances in either the experimental resolution, or the fitting process.

5. Conclusions

We have developed a joint Monte Carlo and machine learning solution to the inverse Milne problem, that extracts electron cross-sections based on electron energy loss spectra from a micro-jet of dense gas. Machine learning was conducted using similar techniques outlined in the literature [41,42,43,44], while a new Monte Carlo simulation of a liquid micro-jet was developed. As a proof-of-concept, we found that the neural network determined neon’s gas phase elastic momentum transfer cross-sections to within 9 % .
The extension towards the simultaneous determination of neon’s ionisation cross-sections, in addition to elastic scattering, decreased the accuracy to within 18 % for elastic scattering, but was accurate within 1 % for the ionisation cross-section. A combinatorial search was conducted using the 100 best elastic and ionisation cross-section pairs, which resulted in an accuracy to within 20 % for the elastic MTCS and 0.2 % for ionisation. While the maximum error for neon’s elastic MTCS increased, overall, the fit was improved, with only the higher energy regime suffering in accuracy. The determination of three or more simultaneous cross-sections resulted in a substantial decrease in accuracy, and is thus left for future iterations of this method.
As expected, absolute magnitudes were difficult to determine, although the prediction of the energy dependencies showed promise. Ionisation cross-sections were remarkably well determined, while the elastic cross-section prediction saw only quite slight discrepancies. Considering a relatively well-formed energy dependence, theoretical values at the higher energies could be used to ‘fix’ each cross-section to improve the current prediction. Alternatively (or in addition), the relative flow technique used for gas phase scattering experiments to achieve absolute cross-sections could be adapted to these liquid micro-jet experiments.
The main limitation with our current technique revolves around information density. Spectra are comprised of several ‘dead zones’, in which little to no scattering information exists. Additionally, a coarse energy resolution was shown to hide important information within the spectra, especially around inelastic peaks with similar threshold energies. To increase information density, a feature extraction algorithm might be utilised for each spectrum along with an improved experimental detector resolution. In other spectral studies, Gaussian de-convolution algorithms are employed to extract peaks from spectra. Utilising a similar approach, along with an increased detector resolution, one could provide further clarity to the network, assuming an appropriate treatment of asymmetry resulting from both multiple scattering and the ionisation energy sharing profiles.
Currently, a fitting algorithm is in development which, when applied to the method outlined in this study, aims to improve the predictive capability of the neural network. Additionally, the inclusion of appropriate liquid structure factors, surface potentials and space charge effects are necessary for real-world derivations of effective electron–liquid cross-sections. Overall, through a proof-of-concept model, we have shown that utilising machine learning, along with a significant wealth of Monte Carlo training data, one can reasonably predict individual and two simultaneous electron scattering cross-sections from EELS. The extension towards predicting full, self-consistent cross-section sets first begins with improvements on the methodology to ensure an improved accuracy in determining the cross-section magnitude.

Author Contributions

Conceptualization, M.J.B. and R.D.W.; Data curation, D.L.M. and P.W.S.; Formal analysis, D.L.M. and P.W.S.; Funding acquisition, M.J.B. and R.D.W.; Investigation, D.L.M. and P.W.S.; Methodology, D.L.M., P.W.S., J.R.G., D.B.J. and R.D.W.; Project administration, M.J.B. and R.D.W.; Software, D.L.M., P.W.S. and D.G.C.; Supervision, R.D.W.; Visualization, D.L.M. and P.W.S.; Writing—original draft, D.L.M. and P.W.S.; Writing—review and editing, D.L.M., P.W.S., D.G.C., J.R.G., D.B.J., M.J.B. and R.D.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors acknowledge the financial assistance of the Australian Research Council (ARC) through its Discovery program (DP180101655, DP190100696, DP220101480). DBJ is the recipient of an Australian Research Council Future Fellowship (FT210100264) funded by the Australian Government.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Benchmarking of the Liquid Micro-Jet Simulation

Throughout the development of the present MC simulation, systematic benchmarking was conducted through comparisons of transport coefficients for swarm experiments under highly non-equilibrium conditions (Appendix A.1). Furthermore, we have developed a relatively simple analytic model for the EELS, and have systematically benchmarked the simulation against this (Appendix A.2).

Appendix A.1. Swarm Benchmarks

In this section, we first present Lucas and Saelee’s electron impact ionisation model to evaluate the accuracy of elastic, inelastic, and ionisation collisions [97]. We then present a modified Percus–Yevic benchmark [98], to benchmark the liquid scattering dynamics discussed in Section 3.1. Lucas and Saelee’s electron impact ionisation model is defined as the following:
σ e l a s t i c = 4 U 1 2 Å 2 , σ i n e l a s t i c = 0.1 ( 1 F ) ( U 15.6 ) Å 2 , U 15.6 eV 0 , U < 15.6 eV σ i o n i s a t i o n = 0.1 F ( U 15.6 ) Å 2 , U 15.6 eV 0 , U < 15.6 e V E / n 0 = 10 Td , m / m 0 = 10 3 , T 0 = 0 K ,
where U is energy in eV, E is the electric field, m is the electron mass, n 0 and m 0 are the background density and mass, respectively, while F is a parameter introduced to control the ratio of ionisation to inelastic collisions and reduced electric field, E / n 0 , is given in Townsend units of Td = 10 21 V m 2 . This model’s strength comes from the fact that the total cross-section is independent of F. Therefore, variation of the ratio between ionisation and inelastic collisions allows for the isolation of the effects produced by each collision type. The results of this application are shown in Table A1 for F values of 0, 0.5 and 1, with excellent agreement being found between the Boltzmann equation and present MC results for each transport coefficient.
Table A1. Swarm benchmark results for Lucas and Saelee’s electron impact ionisation using the current MC code, presented alongside an independent solution to the Boltzmann equation [98]. Provided, for comparison, are the mean energy ( ε ), drift velocity (W) and the transverse ( D T ) and longitudinal ( D L ) diffusion coefficients for each value of F. Note that error estimates are given by the standard deviation of each transport coefficient with respect to time. For this benchmark 10 6 particles were simulated for n 0 t = 2 × 10 16 sm 3 , where averages were conducted over the latter half of the simulation time.
Table A1. Swarm benchmark results for Lucas and Saelee’s electron impact ionisation using the current MC code, presented alongside an independent solution to the Boltzmann equation [98]. Provided, for comparison, are the mean energy ( ε ), drift velocity (W) and the transverse ( D T ) and longitudinal ( D L ) diffusion coefficients for each value of F. Note that error estimates are given by the standard deviation of each transport coefficient with respect to time. For this benchmark 10 6 particles were simulated for n 0 t = 2 × 10 16 sm 3 , where averages were conducted over the latter half of the simulation time.
F ϵ W n 0 D T n 0 D L
[ eV ] [ 10 4 ms 1 ] [ 10 24 m 1 s 1 ] [ 10 24 m 1 s 1 ]
0[98]5.5657.31927.2626.54
Current5.5637.32727.2826.64
Uncertainty0.00010.0030.030.04
0.5[98]5.2248.59327.2628.65
Current5.2238.59427.2328.62
Uncertainty0.0010.0050.040.07
1[98]4.9699.47427.2329.33
Current4.9689.48727.2529.42
Uncertainty0.0020.0080.050.01
We also present results from an adaptation of the Percus–Yevick model liquid benchmark, where the elastic cross-section is held constant and an additional, low energy inelastic cross-section is included:
σ e l a s t i c = 6 Å 2 , σ i n e l a s t i c = 0 , U < 2 eV 0.1 Å 2 , U 2 eV ϕ = 0 , 0.2 , 0.3 , 0.4 m 0 = 4 amu , E / n 0 = 3 Td , T = 0 K ,
where ϕ is the so-called packing factor. For this model, a representative structure factor, detailed in [98], was implemented to verify coherent scattering. The time-averaged values are presented in Table A2, for various values of the packing factor ϕ , which is proportional to the density of the liquid used. Again, excellent agreement between the Boltzmann equation and present MC results, for each transport coefficient, is observed.
Table A2. Swarm benchmark results for our modified Percus–Yevik benchmark (Equation (A2)), from the current MC code, presented alongside an independent solution to the Boltzmann equation [98]. Provided, for comparison, are the mean energies ( ε ), drift velocities (W), and the transverse ( D T ) and longitudinal ( D L ) diffusion coefficients for each value of ϕ . Note that error estimates are given by the standard deviation of each transport coefficient with respect to time. For this benchmark 10 5 particles were simulated for varying times, where averages were conducted over the latter half of the simulation time.
Table A2. Swarm benchmark results for our modified Percus–Yevik benchmark (Equation (A2)), from the current MC code, presented alongside an independent solution to the Boltzmann equation [98]. Provided, for comparison, are the mean energies ( ε ), drift velocities (W), and the transverse ( D T ) and longitudinal ( D L ) diffusion coefficients for each value of ϕ . Note that error estimates are given by the standard deviation of each transport coefficient with respect to time. For this benchmark 10 5 particles were simulated for varying times, where averages were conducted over the latter half of the simulation time.
ϕ ϵ W n 0 D T n 0 D L
[ eV ] [ 10 4 ms 1 ] [ 10 24 m 1 s 1 ] [ 10 24 m 1 s 1 ]
0[98]0.8331.385 2.38
Current0.8341.3842.8252.38
Uncertainty0.00010.0010.0170.02
0.2[98]0.9763.397 6.32
Current0.9773.3889.096.35
Uncertainty0.00010.0030.040.05
0.3[98]1.0805.929 11.2
Current1.0805.91517.9211.1
Uncertainty0.00010.0010.0090.01
0.4[98]1.23310.52 19.51
Current1.23410.5034.9719.51
Uncertainty0.00010.0040.0130.013

Appendix A.2. Beam Benchmark

In this section, we derive a simplified analytic expression for results from reflected EELS measurements from an infinite half-plane given the underlying scattering cross-sections. We first note that, due to threshold excitation collisions, distinct peaks are expected to be observed in the EELS. The relationship between each scattering process (i.e., elastic, excitation and ionisation collisions) and their associated cross-sections, and the area of the corresponding EELS peak can be represented through an infinite series. Given the probability, p n , of an electron escaping the jet after n collisions, the integral of an observed elastic scattering peak over the EELS, I ε , is given by the infinite series:
p 1 σ e l ε c σ t ε c + p 2 σ e l ε c σ t ε c 2 + p 3 σ e l ε c σ t ε c 3 + . . . = 1 N e l I ε d ϵ ,
where the integration bound e l is the peak width associated with pure elastic collision paths, ε is electron energy, σ e l and σ t are the elastic and total cross-sections, respectively, and N is the initial number of electrons. For notation purposes, the left-hand side of Equation (A3) has been simplified as each subsequent collision will result in small energy losses that are encompassed by ε c . Assuming a single collision, Equation (A3) reduces to:
p 1 σ e l ε 0 σ t ε 0 = 1 I 0 e l I ε d ϵ ,
where ε 0 is the initial energy. If we denote the unknown background medium with the subscript u n , and introduce a reference medium r e f , we can then use the ratio of their peaks to arrive at the following relationship:
σ e l a s t i c u n ε 0 σ t o t a l u n ε 0 = p 1 r e f p 1 u n I 0 r e f e l I u n ε d ϵ I 0 u n e l I r e f ε d ϵ σ t o t a l r e f ε 0 σ e l a s t i c r e f ε 0 .
Given the above relationship and assuming that p 1 u n = p 1 r e f = 1 , due to a restriction to single scattering events, such as with beam experiments in the gas phase, we can then derive unknown cross-section sets given the reference cross-section sets and the experimental limitations being considered. Attempts to extend this to include multiple collisions results in an ill posed set of equations [99]. However, under certain assumptions, limited analytical expressions can be derived, which are subsequently used for benchmarking the developed MC simulation package.
Equation (A3) can be used to derive an expression for each peak within the EELS, assuming we know the cross-section set and each probability p n . Using a known cross-section, we thus seek an analytical or numerical solution for p n . First, we consider the case for p 1 . To escape, each electron must travel a distance s 1 without colliding, and then undergo a collision before travelling a distance s towards the surface without undergoing another collision. If we assume a constant cross-section, the probability that a collision will occur at some distance s is given by the following:
p s d s = 1 l exp s l d s ,
where l is the mean free path. Assuming now isotropic scattering, the probability of scattering with azimuthal angle χ and polar angle ψ , with respect to the initial velocity, is given by:
p ( χ ) sin χ d χ = 1 2 sin χ d χ ,
where the ψ dependence has been integrated over and normalised. Each detected electron then escapes the jet by travelling a distance s , defined through:
s = s cos χ .
The probability to not collide in some distance s , is given by the cumulative distribution function (CDF):
P ( s > s ) = 1 0 s 1 l exp s l d s ,
P ( s > s ) = exp s l cos χ .
To couple the spatial and angular probability densities, we multiply Equations (A6), (A7) and (A10) together to obtain the probability that an electron will escape after a collision at distance s and angle χ :
P χ , s sin χ d χ d s = 1 2 l exp s l exp s l cos χ sin χ d χ d s .
After performing a variable change to energy loss u f , and utilising the following:
d u f d χ = 2 m m 0 m + m 0 2 sin χ ,
where m and m 0 are the electron mass and target mass, respectively, then integrating Equation (A11) with respect to u f , s we arrive at the expression for p 1 :
p 1 = 1 2 m + m 0 2 2 m m 0 1 u f d u f .
We then extend the method to calculate p 2 , which requires the use of numerical integration techniques. Each electron now initially travels a distance s 1 and s 2 , before each collision. The detected electron must then travel a distance s , towards the surface to escape. Note that the scattering angles χ and ψ , are now relative to the direction of motion prior to each respective collision. The probability that an electron will undergo each collision at s 1 and s 2 , with scattering angles of χ 1 , ψ 1 and χ 2 , ψ 2 , before escaping after some distance s , is given by:
P s 1 , s 2 , s , χ 1 , χ 2 , ψ 1 , ψ 2 = 1 16 π 2 l 2 exp s 1 l exp s 2 l exp s l ,
where s is represented through a similar CDF as Equation (A10). Assuming an initial position of r = [ 0 , 0 , 0 ] , velocity v = 0 , 0 , v z , and a jet surface along r z = 0 , each detected electron must return to the surface defined by r z = 0 . To achieve an expression for this, we rotate the reference frame of the second collision by χ 1 , about an appropriate unit vector, to the initial frame of reference according to Rodrigues’ rotation formula. The z co-ordinate after two collisions is then given by:
r z = s 1 + s 2 cos χ 1 + s cos χ 1 cos χ 2 sin χ 1 sin χ 2 cos ψ 1 ψ 2 .
Setting r z = 0 , we arrive at an expression for s :
s = s 1 + s 2 cos χ 1 cos χ 1 cos χ 2 sin χ 1 sin χ 2 cos ψ 1 ψ 2 .
Substituting this into Equation (A14) we find:
P s 1 , s 2 , χ 1 , χ 2 , ψ 1 , ψ 2 = 1 16 π 2 l 2 exp s 1 l exp s 2 l × exp 1 l s 1 + s 2 cos χ 1 cos χ 1 cos χ 2 sin χ 1 sin χ 2 cos ψ 1 ψ 2 .
The integration of Equation (A15) over the appropriate limits, along with a transformation of variables to energy loss, gives p 2 . Any non-physical limits were removed using step functions, and they are detailed in what follows. First, bounds for which s < 0 were removed as if s < 0 the electron has already escaped after a single collision and is accounted for in p 1 . Additionally, any divisions by 0 were also removed. Extending this methodology to higher terms is theoretically possible. However, in this investigation we include only p 1 and p 2 due to limitations in the available computation power.
For illustrative purposes and a comparison, a MC simulation of a 15 µ m LµJ was then considered with the following model:
σ e l a s t i c = 1 Å 2 , σ i n e l a s t i c , 1 = 0 , ϵ < 0.3 eV 1 Å 2 , ϵ > 0.3 eV , σ i n e l a s t i c , 2 = 0 , ϵ < 0.5 eV 1 Å 2 , ϵ > 0.5 eV , n 0 = 2.13 × 10 28 m 3 .
Additionally, a numerical integration (using a Monte Carlo technique) was conducted for up to two collisions, and compared with the simulated results in Figure A1. MC simulations that measure both all and up to 2 collisions are included as separate series. Also included as shaded regions are error estimations based upon standard deviations between successive simulations. Excellent agreement was found for up to two collisions. While being somewhat limited in scope, this benchmarking further lends confidence to the accuracy of the simulation.
Figure A1. Comparison of the elastic peak in the EELS using a MC simulation including all collisions (red line), a MC simulation with 2 collisions (black line) and a Monte Carlo integration of 2 collisions (black dots). Shaded areas indicate the standard deviation of successive simulations. Resulting from single collisions, a large clear initial hump is present and consistent between all three models. A sharp drop-off followed by a decaying tail, then represents electrons that underwent multiple collisions. Excellent agreement was found between both the numerical and Monte Carlo method.
Figure A1. Comparison of the elastic peak in the EELS using a MC simulation including all collisions (red line), a MC simulation with 2 collisions (black line) and a Monte Carlo integration of 2 collisions (black dots). Shaded areas indicate the standard deviation of successive simulations. Resulting from single collisions, a large clear initial hump is present and consistent between all three models. A sharp drop-off followed by a decaying tail, then represents electrons that underwent multiple collisions. Excellent agreement was found between both the numerical and Monte Carlo method.
Ijms 23 03354 g0a1

Appendix B. Electron Energy Loss Spectra: Sensitivity to the Scattering Dynamics and Experimental Parameters

In this study we seek an implicit relationship between the EELS and its underlying cross-section set through the use of demonstrative cross-section sets. In what follows we now systematically investigate the ramifications of individual changes to the cross-section sets upon the EELS. For the conditions under investigation, no electron transmission through the beam of 15 µ m was observed.

Appendix B.1. Magnitude and Energy Dependence of Elastic and Excitation Cross-Sections on the EELS

Appendix B.1.1. Effect of the Elastic Cross-Section Dependence on the EELS

In order to de-convolute features present within the EELS and their relationship with the underlying cross-section set, each parameter within the model was isolated and tested for its effect upon the EELS. The first model considered was a constant elastic cross-section:
σ e l a s t i c = σ 0 Å 2 ,
where σ 0 was varied to determine the impact of the cross-section magnitude on the simulated EELS. As expected, an invariance with respect to magnitude was observed as no transmission occurred and thus cross-section magnitude changes will scale the problem without affecting the reflected intensities. The second elastic model considers the effect of the energy dependence of the cross-section. Specifically, in Figure A2 we compare the effect of various power laws, through the elastic only model:
σ e l a s t i c = ε l Å 2 ,
where ε is energy in eV and l is varied in accordance with the legend of Figure A2. Each spectrum showed an invariance with respect to the power law at low energy losses, and exhibited intensity decay rates at higher energy loss that were dependent upon the power law used. The magnitude of this effect depends on the energy dependence of the cross-section’s derivative, and gradual changes were found to have a minimal impact upon the spectra.
Figure A2. Electron Energy Loss Spectra (a), and their corresponding elastic cross-sections (b), for various power laws (see legend of (a)). An initial energy of 50 eV was used. Each spectrum initially remains the same between each power law, before a divergence is observed at higher energy loss. Each spectrum for a single cross-section remains independent of the cross-section magnitude; however evident here is a dependence upon each cross-section derivative.
Figure A2. Electron Energy Loss Spectra (a), and their corresponding elastic cross-sections (b), for various power laws (see legend of (a)). An initial energy of 50 eV was used. Each spectrum initially remains the same between each power law, before a divergence is observed at higher energy loss. Each spectrum for a single cross-section remains independent of the cross-section magnitude; however evident here is a dependence upon each cross-section derivative.
Ijms 23 03354 g0a2
To explain this physically, one can imagine the swarm of electrons as an expanding sphere in Cartesian space. Figure A3 shows the electron swarm within the jet at time 2.83 × 10 11 s after impact with a LµJ for l 1 , 0 , 1 . A constant mean free path with energy loss ( l = 0 ) results in considerable back scattering after jet surface impact, and the remaining electron cloud travels further inward with a constant expansion rate. Nothing inherently changes during their flight and the back-scattering flux is proportional to the number of remaining electrons within the jet. A decreasing mean free path with energy loss ( l = 1 ) results in a diminishing rate of expansion, that results in a relatively smaller back-scattered intensity. In contrast to this, an increasing mean free path with energy loss ( l = 1 ) accelerates the rate of expansion, resulting in an increased likelihood of back scattering out of the jet.
Figure A3. Heatmap of the electron cloud at 2.83 × 10 11 s after impact, with the LµJ surface for l = 1 (a), l = 0 (b) and l = 1 (c) according to Equation (A17). For this model, 10 7 1 eV electrons were simulated through a 15 µ m LµJ with the jet directed into the page. The electron counts for each corresponding depth and horizontal displacement were then recorded. A log 10 transform of non-zero values was then performed for visualisation purposes. Initially, in each scenario, the elastic cross-section is equal due to the initial energy of 1 eV . In (a), the electron swarm loses energy through collisions with the background gas and experiences an increasing cross-section (decreasing mean free path) with energy loss, which results in a decreased expansion of the electron swarm. In (b), the electron swarm experiences a constant cross-section (constant mean free path) with energy loss, which results in a steady expansion of the electron swarm. Finally, in (c), the electron swarm experiences a decreasing cross-section (increasing mean free path) with energy loss, which results in an increased expansion of the electron swarm.
Figure A3. Heatmap of the electron cloud at 2.83 × 10 11 s after impact, with the LµJ surface for l = 1 (a), l = 0 (b) and l = 1 (c) according to Equation (A17). For this model, 10 7 1 eV electrons were simulated through a 15 µ m LµJ with the jet directed into the page. The electron counts for each corresponding depth and horizontal displacement were then recorded. A log 10 transform of non-zero values was then performed for visualisation purposes. Initially, in each scenario, the elastic cross-section is equal due to the initial energy of 1 eV . In (a), the electron swarm loses energy through collisions with the background gas and experiences an increasing cross-section (decreasing mean free path) with energy loss, which results in a decreased expansion of the electron swarm. In (b), the electron swarm experiences a constant cross-section (constant mean free path) with energy loss, which results in a steady expansion of the electron swarm. Finally, in (c), the electron swarm experiences a decreasing cross-section (increasing mean free path) with energy loss, which results in an increased expansion of the electron swarm.
Ijms 23 03354 g0a3

Appendix B.1.2. Effect of the Excitation Cross-Section Dependence on the EELS

Excitation cross-sections are introduced through the model:
σ e l a s t i c = 10 Å 2 , σ e x c i t a i t o n = 0 , ε < 5 eV σ 0 Å 2 , ε > 5 eV ,
where σ 0 is varied to investigate the relative elastic/excitation cross-section magnitudes and their effects on the EELS. As shown in Figure A4, it is clear that the intensity peak at 0 eV energy loss is proportional to the elastic cross-section. Inelastic peaks are then observed at intervals equal to the 5 eV threshold energy. Once the energy of the electrons falls below that threshold, energy loss is minimal due to the restriction to elastic collisions. This results in a high energy loss peak that is proportional to the elastic mean free path. Additionally, as the inelastic cross-section magnitude increases, the below threshold peak at 45 eV energy loss also increases due to an increased energy loss rate. Note that while the initial energy in all cases is sampled from a Gaussian, energy loss is defined as the difference between the average initial energy and final energy, which can result in an apparent negative energy loss in some cases.
As expected, the relative elastic/excitation cross-section magnitude is correlated with each peak found in the EELS. As shown previously, in a single collision regime and within an ideal environment, direct relationships exist between the two quantities. However, the introduction of multiple scattering convolutes their dependence. While a clear trend also exists here, excitation processes within a multiple scattering environment introduce peak dependencies upon relative cross-section magnitudes at multiple energies and results in a large combinatorial problem.
Figure A4. EELS for elastic and a single inelastic cross-section, with varying relative magnitude. An initial energy of 50 eV was used, with an inelastic threshold of 5 eV . The inelastic cross-section magnitude is shown in the legend. At 0 eV energy loss an elastic peak exists, followed by threshold separated sequential inelastic peaks. Additionally, after 45 eV of energy loss, electrons exclusively undergo elastic collisions which results in an increased intensity due to their minimal energy loss. Evident in this figure is a relative elastic to inelastic ratio dependence of the EELS upon the relative magnitude of the respective cross-sections.
Figure A4. EELS for elastic and a single inelastic cross-section, with varying relative magnitude. An initial energy of 50 eV was used, with an inelastic threshold of 5 eV . The inelastic cross-section magnitude is shown in the legend. At 0 eV energy loss an elastic peak exists, followed by threshold separated sequential inelastic peaks. Additionally, after 45 eV of energy loss, electrons exclusively undergo elastic collisions which results in an increased intensity due to their minimal energy loss. Evident in this figure is a relative elastic to inelastic ratio dependence of the EELS upon the relative magnitude of the respective cross-sections.
Ijms 23 03354 g0a4
Shown in Figure A5, we expand the previous model to include energy dependencies with multiple cross-sections through the following model:
σ e l a s t i c = ε x Å 2 , σ i n e l a s t i c = 0 , ε < 5 eV ε x Å 2 , ε > 5 eV ,
where x is varied according to the legend of Figure A5. No change in the elastic peak at 0 eV is observed, however each subsequent peak decays at a rate dependent upon the power law used.
Figure A5. Simulated EELS (a), and area under each peak (b), as a function of the electron energy loss for a single elastic and inelastic cross-section, but both now with an equal power law functional form defined by the legend in (b). A distinct trend is observed, where the power law is proportional to the each cross-section’s slope.
Figure A5. Simulated EELS (a), and area under each peak (b), as a function of the electron energy loss for a single elastic and inelastic cross-section, but both now with an equal power law functional form defined by the legend in (b). A distinct trend is observed, where the power law is proportional to the each cross-section’s slope.
Ijms 23 03354 g0a5
For comparison, integration over the fitted exponential modified Gaussians for each peak are shown in the right hand graph of Figure A5 (note that for x = 1 , the 40 eV peak was not found by the fitting algorithm). In contrast to the single elastic cross-section result outlined in Appendix B.1.1, a distinct correlation was found. The introduction of an inelastic process dramatically increases the energy loss rate, which in turn will amplify the effects of an energy-dependent cross-section. A direct relationship remains convoluted, however, as the magnitude and energy dependence of the individual cross-sections, within a given cross-section set, will vary, resulting in significant degeneracy issues.
Of particular interest is the below threshold peak at 45 eV in the EELS. Here it was found that while their areas remain relatively constant, there exists a relationship between their asymmetric nature and the power law used. The exact nature of the relationship depends upon each cross-section’s, and hence mean free path’s, derivative with respect to energy. As discussed in Appendix B.1.1, when compared to scattering with a constant mean free path, an increasing mean free path with each subsequent collision increases the expansion rate for the electron swarm due to the reduced collisional impediment, resulting in a larger scattered intensity. As the swarm falls below the threshold energies, where energy loss is minimal, fewer electrons remain in the liquid when compared with the result for the constant mean free path model and thus a decreased scattered intensity was observed. A decreasing mean free path with energy loss slows the rate at which the swarm expands, which results in a reduced intensity when compared to a constant mean free path. As the mean free path decreases even further at lower energies, the probability of escape reduces even further, which effectively traps the electron swarm within the liquid.

Appendix B.1.3. Effect of the Ionisation Cross-Section Dependence on the EELS

To investigate ionisation collisions, and their effect upon the EELS, the following model was used:
σ e l a s t i c = 1 Å 2 , σ i o n i s a t i o n = 0 , ε < 25 eV 1 Å 2 , ε > 25 eV ,
with varying ionisation energy-sharing profiles between the primary and secondary electrons as shown in Figure A6. At this point, it is important to consider the effects of fractional energy loss on the multiple scattering features of the EELS. Energy loss is largely driven by threshold excitations, energy partitioning losses (i.e., distribution of the excess energy between the scattered and ejected electrons) and kinetic energy exchange between the target and the electron (see Equation (A12)). Energy exchange losses are proportional to the incident energy of the electron (and the mass ratio with the target atom/molecules), and result in a narrowing of the features (i.e., a reduction in average energy loss) at low energies.
While excitation collisions introduce relatively simple peaks within the spectra, ionisation collisions introduce unique features depending upon the energy-sharing profile between the primary and secondary electrons. Four model scenarios were considered (see Figure A6; equal, biased, uniform and skewed sharing). Equal energy sharing (i.e., a 50–50% energy partition) between the primary and secondary electrons resulted in a single peak at an energy loss of the ionisation threshold plus half the ionisation threshold energy. Once below the energy threshold, electrons undergo elastic collisions at minimal energy loss resulting in a relatively large ionisation peak. Biased energy sharing was represented through a 90–10% sharing partition, resulting in separate peaks that exhibit broadening due to energy exchange losses. While the observed amplitudes and full width half maximums vary according to the aforementioned narrowing effect, the total peak area was equal in accordance with the equal cross-sectional magnitudes used (see Figure A6).
Figure A6. EELS for an elastic scattering and ionisation model. Initial energy was set to 50 eV , with an ionisation threshold of 25 eV . The energy sharing ratios, between the primary and secondary electrons, were set to 50–50% (a), 90–10% (b), uniformly randomised (c), and skewed randomisation in accordance to a normalised exponential cumulative distribution function (d).
Figure A6. EELS for an elastic scattering and ionisation model. Initial energy was set to 50 eV , with an ionisation threshold of 25 eV . The energy sharing ratios, between the primary and secondary electrons, were set to 50–50% (a), 90–10% (b), uniformly randomised (c), and skewed randomisation in accordance to a normalised exponential cumulative distribution function (d).
Ijms 23 03354 g0a6
As expected, a uniform probability of all sharing ratios resulted in a shoulder beginning at the ionisation threshold energy loss. A gradual increase was observed in the shoulder intensity once again, due to the energy exchange losses discussed above. Finally, a normalised exponential cumulative distribution function was utilised to emulate a skewed energy sharing distribution for the ionisation process. A semi-symmetric saddle shape was observed, with an asymmetry again driven through energy exchange losses.

Appendix B.2. Anisotropy in the Scattering Dynamics

Appendix B.2.1. Coherent Scattering Effects

As discussed previously, modelling coherent scattering in liquid environments was achieved through the addition of liquid structure factors [47,85]. Figure A7 compares how those coherent liquid scattering structure factors affect the present EELS utilising the model:
σ e l a s t i c = 6 Å 2 , σ i n e l a s t i c = 0 , ε < 2 eV 0.1 Å 2 , ε > 2 eV ,
with the Percus-Yevick structure factor [47]. Coherent scattering effects are isolated to the low energy regime (<5 eV), and hence its effects were limited outside that range. At low energies, electrons tend to coherently scatter through multiple targets, with little angle change, and hence will travel further in a particular direction when compared to isotropic scattering binary collisions. It was found that initial energies below 5 eV resulted in a decreased surface scattering intensity, due to enhanced forward scattering. Multiple scattering then randomises the preferred direction, such that the swarm will expand in size faster than the isotropic equivalent, which resulted in an increased intensity observed at the higher energy loss.
Figure A7. An elastic and inelastic EELS, with the Percus-Yevick liquid structure factor result (Black) compared with the equivalent gas phase result (Red). In each model, an inelastic threshold of 2 eV was used with incident energies of 5 eV (a), and 2.5 eV (b). As coherent scattering is predominantly restricted to low energies (<5 eV), its effects are minimal when the incident energies exist outside that domain.
Figure A7. An elastic and inelastic EELS, with the Percus-Yevick liquid structure factor result (Black) compared with the equivalent gas phase result (Red). In each model, an inelastic threshold of 2 eV was used with incident energies of 5 eV (a), and 2.5 eV (b). As coherent scattering is predominantly restricted to low energies (<5 eV), its effects are minimal when the incident energies exist outside that domain.
Ijms 23 03354 g0a7

Appendix B.2.2. Differential Cross-Section Effects

Anisotropic scattering was included through a differential cross-section (DCS) in the following model scenario:
σ e l a s t i c = 6 Å 2 , σ i n e l a s t i c = 0 , ε < 2 eV 0.1 Å 2 , ε > 2 eV ,
where varying degrees of forward scattering were considered and represented through a normalised powers of the cos χ function, where χ is relative to the incident direction of motion. Each modified DCS was applied to both the elastic and inelastic cross-sections.
Similar to coherent scattering effects, an initial decrease in intensity at 0 eV energy loss was observed in Figure A8 along with an increased signal at higher energy loss that is proportional to the degree of forward scattering. The explanation of that latter observation is the same as in Appendix B.2.1, except now it is applicable to all incident energies.
Figure A8. Elastic and inelastic EELS, with a varying differential cross-section. An initial energy of 10 eV is used with an inelastic threshold of 2 eV . Each DCS is represented using a normalised power of cos χ , as indicated in the legend. For higher degrees of forward scattering, a decrease/increase was observed at lower/higher energy loss.
Figure A8. Elastic and inelastic EELS, with a varying differential cross-section. An initial energy of 10 eV is used with an inelastic threshold of 2 eV . Each DCS is represented using a normalised power of cos χ , as indicated in the legend. For higher degrees of forward scattering, a decrease/increase was observed at lower/higher energy loss.
Ijms 23 03354 g0a8

References

  1. Shimizu, T.; Steffes, B.; Pompl, R.; Jamitzky, F.; Bunk, W.; Ramrath, K.; Georgi, M.; Stolz, W.; Schmidt, H.; Urayama, T.; et al. Characterization of microwave plasma torch for decontamination. Plasma Process. Polym. 2008, 5, 577–582. [Google Scholar] [CrossRef]
  2. Tanaka, H.; Mizuno, M.; Ishikawa, K.; Nakamura, K.; Kajiyama, H.; Kano, H.; Kikkawa, F.; Hori, M. Plasma-activated medium selectively kills glioblastoma brain tumor cells by down-regulating a survival signaling molecule, AKT kinase. Plasma Med. 2011, 1, 265–277. [Google Scholar] [CrossRef] [Green Version]
  3. Taylor, D.W.; Petrera, M.; Hendry, M.; Theodoropoulos, J.S. A systematic review of the use of platelet-rich plasma in sports medicine as a new treatment for tendon and ligament injuries. Clin. J. Sport Med. 2011, 21, 344–352. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Schlegel, J.; Köritzer, J.; Boxhammer, V. Plasma in cancer treatment. Clin. Plasma Med. 2013, 1, 2–7. [Google Scholar] [CrossRef]
  5. Keidar, M. Plasma for cancer treatment. Plasma Sources Sci. Technol. 2015, 24, 033001. [Google Scholar] [CrossRef]
  6. Bekeschus, S.; Schmidt, A.; Weltmann, K.; von Woedtke, T. The plasma jet kINPen—A powerful tool for wound healing. Clin. Plasma Med. 2016, 4, 19–28. [Google Scholar] [CrossRef]
  7. Bernhardt, T.; Semmler, M.L.; Schäfer, M.; Bekeschus, S.; Emmert, S.; Boeckmann, L. Plasma medicine: Applications of cold atmospheric pressure plasma in dermatology. Oxidative Med. Cell. Longev. 2019, 2019, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. von Woedtke, T.; Schmidt, A.; Bekeschus, S.; Wende, K.; Weltmann, K. Plasma medicine: A field of applied redox biology. In Vivo 2019, 33, 1011–1026. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Liu, D.; Szili, E.J.; Ostrikov, K. Plasma medicine: Opportunities for nanotechnology in a digital age. Plasma Process. Polym. 2020, 17, 2000097. [Google Scholar] [CrossRef]
  10. Adamovich, I.; Baalrud, S.D.; Bogaerts, A.; Bruggeman, P.J.; Cappelli, M.; Colombo, V.; Czarnetzki, U.; Ebert, U.; Eden, J.G.; Favia, P.; et al. The 2017 Plasma Roadmap: Low temperature plasma science and technology. J. Phys. Appl. Phys. 2017, 50, 323001. [Google Scholar] [CrossRef]
  11. Pancheshnyi, S.; Biagi, S.; Bordage, M.C.; Hagelaar, G.J.; Morgan, W.L.; Phelps, A.V.; Pitchford, L.C. The LXCat project: Electron scattering cross sections and swarm parameters for low temperature plasma modeling. Chem. Phys. 2012, 398, 148–153. [Google Scholar] [CrossRef]
  12. Pitchford, L.C.; Alves, L.L.; Bartschat, K.; Biagi, S.F.; Bordage, M.C.; Bray, I.; Brion, C.E.; Brunger, M.J.; Campbell, L.; Chachereau, A.; et al. LXCat: An Open-Access, Web-Based Platform for Data Needed for Modeling Low Temperature Plasmas. Plasma Process. Polym. 2017, 14, 1600098. [Google Scholar] [CrossRef]
  13. Boyle, G.J.; McEachran, R.P.; Cocks, D.G.; White, R.D. Electron scattering and transport in liquid argon. J. Chem. Phys. 2015, 142, 154507. [Google Scholar] [CrossRef] [PubMed]
  14. Boyle, G.J.; McEachran, R.P.; Cocks, D.G.; Brunger, M.J.; Buckman, S.J.; Dujko, S.; White, R.D. Ab initio electron scattering cross-sections and transport in liquid xenon. J. Phys. Appl. Phys. 2016, 49, 355201. [Google Scholar] [CrossRef] [Green Version]
  15. Meesungnoen, J.; Jay-Gerin, J.; Filali-Mouhim, A.; Mankhetkorn, S. Low-energy electron penetration range in liquid water. Radiat. Res. 2002, 158, 657–660. [Google Scholar] [CrossRef]
  16. Michaud, M.; Wen, A.; Sanche, L. Cross sections for low-energy (1–100 eV) electron elastic and inelastic scattering in amorphous ice. Radiat. Res. 2003, 159, 3–22. [Google Scholar] [CrossRef]
  17. Brunger, M.J.; Buckman, S.J. Electron–molecule scattering cross-sections. I. Experimental techniques and data for diatomic molecules. Phys. Rep. 2002, 357, 215–458. [Google Scholar] [CrossRef]
  18. Faubel, M.; Steiner, B. Photoelectron spectroscopy at liquid water surfaces. In Linking the Gaseous and Condensed Phases of Matter; Springer: Berlin/Heidelberg, Germany, 1994; pp. 517–523. [Google Scholar]
  19. Faubel, M.; Steiner, B.; Toennies, J.P. Photoelectron spectroscopy of liquid water, some alcohols, and pure nonane in free micro jets. J. Chem. Phys. 1997, 106, 9013–9031. [Google Scholar] [CrossRef]
  20. Faubel, M.; Steiner, B.; Toennies, J.P. Measurement of He I photoelectron spectra of liquid water, formamide and ethylene glycol in fast-flowing microjets. J. Electron Spectrosc. Relat. Phenom. 1998, 95, 159–169. [Google Scholar] [CrossRef]
  21. Winter, B.; Weber, R.; Hertel, I.V.; Faubel, M.; Jungwirth, P.; Brown, E.C.; Bradforth, S.E. Electron binding energies of aqueous alkali and halide ions: EUV photoelectron spectroscopy of liquid solutions and combined ab initio and molecular dynamics calculations. J. Am. Chem. Soc. 2005, 127, 7203–7214. [Google Scholar] [CrossRef]
  22. Winter, B.; Aziz, E.F.; Hergenhahn, U.; Faubel, M.; Hertel, I.V. Hydrogen bonds in liquid water studied by photoelectron spectroscopy. J. Chem. Phys. 2007, 126, 124504. [Google Scholar] [CrossRef] [PubMed]
  23. Winter, B. Liquid microjet for photoelectron spectroscopy. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 2009, 601, 139–150. [Google Scholar] [CrossRef]
  24. Ottosson, N.; Faubel, M.; Bradforth, S.E.; Jungwirth, P.; Winter, B. Photoelectron spectroscopy of liquid water and aqueous solution: Electron effective attenuation lengths and emission-angle anisotropy. J. Electron Spectrosc. Relat. Phenom. 2010, 177, 60–70. [Google Scholar] [CrossRef]
  25. Thürmer, S.; Seidel, R.; Faubel, M.; Eberhardt, W.; Hemminger, J.C.; Bradforth, S.E.; Winter, B. Photoelectron angular distributions from liquid water: Effects of electron scattering. Phys. Rev. Lett. 2013, 111, 173005. [Google Scholar] [CrossRef]
  26. Brown, M.A.; Redondo, A.B.; Jordan, I.; Duyckaerts, N.; Lee, M.; Ammann, M.; Nolting, F.; Kleibert, A.; Huthwelker, T.; Mächler, J.; et al. A new endstation at the Swiss Light Source for ultraviolet photoelectron spectroscopy, X-ray photoelectron spectroscopy, and X-ray absorption spectroscopy measurements of liquid solutions. Rev. Sci. Instrum. 2013, 84, 073904. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Riley, J.W.; Wang, B.; Parkes, M.A.; Fielding, H.H. Design and characterization of a recirculating liquid-microjet photoelectron spectrometer for multiphoton ultraviolet photoelectron spectroscopy. Rev. Sci. Instrum. 2019, 90, 083104. [Google Scholar] [CrossRef] [Green Version]
  28. Buttersack, T.; Mason, P.E.; McMullen, R.S.; Martinek, T.; Brezina, K.; Hein, D.; Ali, H.; Kolbeck, C.; Schewe, C.; Malerz, S.; et al. Valence and Core-Level X-ray Photoelectron Spectroscopy of a Liquid Ammonia Microjet. J. Am. Chem. Soc. 2019, 141, 1838–1841. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Mudryk, K.D.; Seidel, R.; Winter, B.; Wilkinson, I. The electronic structure of the aqueous permanganate ion: Aqueous-phase energetics and molecular bonding studied using liquid jet photoelectron spectroscopy. Phys. Chem. Chem. Phys. 2020, 22, 20311–20330. [Google Scholar] [CrossRef] [PubMed]
  30. Nishitani, J.; Karashima, S.; West, C.W.; Suzuki, T. Surface potential of liquid microjet investigated using extreme ultraviolet photoelectron spectroscopy. J. Chem. Phys. 2020, 152, 144503. [Google Scholar] [CrossRef] [PubMed]
  31. Huxley, L.G.H.; Crompton, R.W. Diffusion and Drift of Electrons in Gases; Wiley Series in Plasma Physics: Hoboken, NJ, USA, 1974. [Google Scholar]
  32. Crompton, R.W. Benchmark measurements of cross sections for electron collisions: Electron swarm methods. Adv. At. Mol. Opt. Phys. 1994, 33, 97–148. [Google Scholar]
  33. Petrović, Z.L.; Šuvakov, M.; Nikitović, Ž.; Dujko, S.; Šašić, O.; Jovanović, J.; Malović, G.; Stojanović, V. Kinetic phenomena in charged particle transport in gases, swarm parameters and cross section data. Plasma Sources Sci. Technol. 2007, 16, S1. [Google Scholar] [CrossRef] [Green Version]
  34. Petrović, Z.L.; Banković, A.; Dujko, S.; Marjanović, S.; Malović, G.; Sullivan, J.P.; Buckman, S.J. Data for modeling of positron collisions and transport in gases. AIP Conf. Proc. 2013, 1545, 115–131. [Google Scholar]
  35. Ness, K.F.; Robson, R.E. Transport properties of electrons in water vapor. Phys. Rev. A 1988, 38, 1446. [Google Scholar] [CrossRef] [PubMed]
  36. Yousfi, M.; Benabdessadok, M.D. Boltzmann equation analysis of electron-molecule collision cross sections in water vapor and ammonia. J. Appl. Phys. 1996, 80, 6619–6630. [Google Scholar] [CrossRef]
  37. Robson, R.E.; White, R.D.; Ness, K.F. Transport coefficients for electrons in water vapor: Definition, measurement, and calculation. J. Chem. Phys. 2011, 134, 064319. [Google Scholar] [CrossRef]
  38. Ness, K.F.; Robson, R.E.; Brunger, M.J.; White, R.D. Transport coefficients and cross sections for electrons in water vapour: Comparison of cross section sets using an improved Boltzmann equation solution. J. Chem. Phys. 2012, 136, 024318. [Google Scholar] [CrossRef]
  39. De Urquijo, J.; Basurto, E.; Juárez, A.M.; Ness, K.F.; Robson, R.E.; Brunger, M.J.; White, R.D. Electron drift velocities in He and water mixtures: Measurements and an assessment of the water vapour cross-section sets. J. Chem. Phys. 2014, 141, 014308. [Google Scholar] [CrossRef]
  40. Casey, M.J.E.; De Urquijo, J.; Serkovic Loli, L.N.; Cocks, D.G.; Boyle, G.J.; Jones, D.B.; Brunger, M.J.; White, R.D. Self-consistency of electron-THF cross sections using electron swarm techniques. J. Chem. Phys. 2017, 147, 195103. [Google Scholar] [CrossRef]
  41. Stokes, P.W.; Casey, M.J.E.; Cocks, D.G.; de Urquijo, J.; García, G.; Brunger, M.J.; White, R.D. Self-consistent electron–THF cross sections derived using data-driven swarm analysis with a neural network model. Plasma Sources Sci. Technol. 2020, 29, 105008. [Google Scholar] [CrossRef]
  42. Stokes, P.W.; Cocks, D.G.; Brunger, M.J.; White, R.D. Determining cross sections from transport coefficients using deep neural networks. Plasma Sources Sci. Technol. 2020, 29, 055009. [Google Scholar] [CrossRef] [Green Version]
  43. Stokes, P.W.; Foster, S.P.; Casey, M.J.E.; Cocks, D.G.; González-Magaña, O.; de Urquijo, J.; García, G.; Brunger, M.J.; White, R.D. An improved set of electron-THFA cross sections refined through a neural network-based analysis of swarm data. J. Chem. Phys. 2021, 154, 084306. [Google Scholar] [CrossRef] [PubMed]
  44. Stokes, P.W.; White, R.D.; Campbell, L.; Brunger, M.J. Toward a complete and comprehensive cross section database for electron scattering from NO using machine learning. J. Chem. Phys. 2021, 155, 084305. [Google Scholar] [CrossRef] [PubMed]
  45. Nam, J.; Yong, H.; Hwang, J.; Choi, J. Training an artificial neural network for recognizing electron collision patterns. Phys. Lett. A 2021, 387, 127005. [Google Scholar] [CrossRef]
  46. Jetly, V.; Chaudhury, B. Extracting Electron Scattering Cross Sections from Swarm Data using Deep Neural Networks. Mach. Learn. Sci. Technol. 2021, 2, 035025. [Google Scholar] [CrossRef]
  47. Tattersall, W.J.; Cocks, D.G.; Boyle, G.J.; Buckman, S.J.; White, R.D. Monte Carlo study of coherent scattering effects of low-energy charged particle transport in Percus-Yevick liquids. Phys. Rev. E 2015, 91, 043304. [Google Scholar] [CrossRef] [Green Version]
  48. Siegbahn, H.; Siegbahn, K. ESCA applied to liquids. J. Electron Spectrosc. Relat. Phenom. 1973, 2, 319–325. [Google Scholar] [CrossRef]
  49. Faubel, M.; Schlemmer, S.; Toennies, J.P. A molecular beam study of the evaporation of water from a liquid jet. Z. Phys. D Atoms Mol. Clust. 1988, 10, 269–277. [Google Scholar] [CrossRef]
  50. Zahoor, R.; Bajt, S.; Šarler, B. Influence of gas dynamic virtual nozzle geometry on micro-jet characteristics. Int. J. Multiph. Flow 2018, 104, 152–165. [Google Scholar] [CrossRef]
  51. Zahoor, R.; Bajt, S.; Šarler, B. Numerical investigation on influence of focusing gas type on liquid micro-jet characteristics. Int. J. Hydromechatronics 2018, 1, 222–237. [Google Scholar] [CrossRef]
  52. Zahoor, R.; Belšak, G.; Bajt, S.; Šarler, B. Simulation of liquid micro-jet in free expanding high-speed co-flowing gas streams. Microfluid. Nanofluidics 2018, 22, 1–20. [Google Scholar] [CrossRef]
  53. Ekimova, M.; Quevedo, W.; Faubel, M.; Wernet, P.; Nibbering, E.T. A liquid flatjet system for solution phase soft-x-ray spectroscopy. Struct. Dyn. 2015, 2, 054301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Grisenti, R.E.; Fraga, R.A.C.; Petridis, N.; Dörner, R.; Deppe, J. Cryogenic microjet for exploration of superfluidity in highly supercooled molecular hydrogen. EPL (Europhys. Lett.) 2006, 73, 540. [Google Scholar] [CrossRef]
  55. Kim, J.B.; Göde, S.; Glenzer, S.H. Development of a cryogenic hydrogen microjet for high-intensity, high-repetition rate experiments. Rev. Sci. Instrum. 2016, 87, 11E328. [Google Scholar] [CrossRef] [PubMed]
  56. Holstein, W.L.; Hammer, M.R.; Metha, G.F.; Buntine, M.A. Time-of-flight mass spectrometric detection of mono-and di-substituted benzenes at parts per million concentrations by way of liquid microjet injection and laser ionisation. Int. J. Mass Spectrom. 2001, 207, 1–12. [Google Scholar] [CrossRef]
  57. Hemberg, O.; Otendal, M.; Hertz, H.M. Liquid-metal-jet anode electron-impact x-ray source. Appl. Phys. Lett. 2003, 83, 1483–1485. [Google Scholar] [CrossRef] [Green Version]
  58. Higashiguchi, T.; Hamada, M.; Kubodera, S. Development of a liquid tin microjet target for an efficient laser-produced plasma extreme ultraviolet source. Rev. Sci. Instrum. 2007, 78, 036106. [Google Scholar] [CrossRef] [PubMed]
  59. Ueno, Y.; Ariga, T.; Soumagne, G.; Higashiguchi, T.; Kubodera, S.; Pogorelsky, I.; Pavlishin, I.; Stolyarov, D.; Babzien, M.; Kusche, K.; et al. Efficient extreme ultraviolet plasma source generated by a CO2 laser and a liquid xenon microjet target. Appl. Phys. Lett. 2007, 90, 191503. [Google Scholar] [CrossRef]
  60. Lar’Kin, A.; Uryupina, D.; Ivanov, K.; Savel’Ev, A.; Bonnet, T.; Gobet, F.; Hannachi, F.; Tarisien, M.; Versteegen, M.; Spohr, K.; et al. Microjet formation and hard x-ray production from a liquid metal target irradiated by intense femtosecond laser pulses. Phys. Plasmas 2014, 21, 093103. [Google Scholar] [CrossRef] [Green Version]
  61. Lange, K.M.; Könnecke, R.; Ghadimi, S.; Golnak, R.; Soldatov, M.A.; Hodeck, K.F.; Soldatov, A.; Aziz, E.F. High resolution X-ray emission spectroscopy of water and aqueous ions using the micro-jet technique. Chem. Phys. 2010, 377, 1–5. [Google Scholar] [CrossRef]
  62. Lange, K.M.; Soldatov, M.; Golnak, R.; Gotz, M.; Engel, N.; Könnecke, R.; Rubensson, J.; Aziz, E.F. X-ray emission from pure and dilute H2O and D2O in a liquid microjet: Hydrogen bonds and nuclear dynamics. Phys. Rev. B 2012, 85, 155104. [Google Scholar] [CrossRef] [Green Version]
  63. Dierker, B.; Suljoti, E.; Atak, K.; Lange, K.M.; Engel, N.; Golnak, R.; Dantz, M.; Hodeck, K.; Khan, M.; Kosugi, N.; et al. Probing orbital symmetry in solution: Polarization-dependent resonant inelastic soft X-ray scattering on liquid micro-jet. New J. Phys. 2013, 15, 093025. [Google Scholar] [CrossRef] [Green Version]
  64. Maselli, O.J.; Gascooke, J.R.; Lawrance, W.D.; Buntine, M.A. The dynamics of evaporation from a liquid surface. Chem. Phys. Lett. 2011, 513, 1–11. [Google Scholar] [CrossRef]
  65. Duffin, A.M.; Saykally, R.J. Electrokinetic power generation from liquid water microjets. J. Phys. Chem. C 2008, 112, 17018–17022. [Google Scholar] [CrossRef] [Green Version]
  66. Li, C.; Meng, P.; Jiang, H.; Hu, X. Power generation from microjet array of liquid water. J. Phys. D Appl. Phys. 2018, 51, 285501. [Google Scholar] [CrossRef]
  67. Jang, H.; Yu, H.; Lee, S.; Hur, E.; Kim, Y.; Lee, S.; Kang, N.; Yoh, J.J. Towards clinical use of a laser-induced microjet system aimed at reliable and safe drug delivery. J. Biomed. Opt. 2014, 19, 058001. [Google Scholar] [CrossRef] [PubMed]
  68. Winter, B.; Faubel, M. Photoemission from liquid aqueous solutions. Chem. Rev. 2006, 106, 1176–1211. [Google Scholar] [CrossRef]
  69. Seidel, R.; Thurmer, S.; Winter, B. Photoelectron spectroscopy meets aqueous solution: Studies from a vacuum liquid microjet. J. Phys. Chem. Lett. 2011, 2, 633–641. [Google Scholar] [CrossRef]
  70. Faubel, M. Photoelectron spectroscopy at liquid surfaces. In Photoionization and Photodetachment: In 2 Parts; World Scientific: Singapore, 2000; pp. 634–690. [Google Scholar]
  71. Rayleigh, J.W.S.; Rayleigh, J.W.S.B. The Theory of Sound; Courier Corporation: London, UK, 1945; Volume 2. [Google Scholar]
  72. Bockris, J.; Devanathan, M.; Müller, K. On the structure of charged interfaces. In Electrochemistry; Elsevier: Amsterdam, The Netherlands, 1965; pp. 832–863. [Google Scholar]
  73. Kurahashi, N.; Karashima, S.; Tang, Y.; Horio, T.; Abulimiti, B.; Suzuki, Y.; Ogi, Y.; Oura, M.; Suzuki, T. Photoelectron spectroscopy of aqueous solutions: Streaming potentials of NaX (X = Cl, Br, and I) solutions and electron binding energies of liquid water and X. J. Chem. Phys. 2014, 140, 174506. [Google Scholar] [CrossRef] [PubMed]
  74. Cavanagh, S.; Lohmann, B. Coplanar asymmetric (e, 2e) measurements of ionization of N2O. J. Phys. At. Mol. Opt. Phys. 1999, 32, L261. [Google Scholar] [CrossRef]
  75. McCormick, N.J. Inverse radiative transfer problems: A review. Nucl. Sci. Eng. 1992, 112, 185–198. [Google Scholar] [CrossRef]
  76. Case, K.M. Inverse problem in transport theory. Phys. Fluids 1973, 16, 1607–1611. [Google Scholar] [CrossRef]
  77. Davison, B.; Sykes, J.B. Neutron Transport Theory; Clarendon Press: Oxford, UK, 1957. [Google Scholar]
  78. Larsen, E.W.; Keller, J.B. Asymptotic solution of neutron transport problems for small mean free paths. J. Math. Phys. 1974, 15, 75–81. [Google Scholar] [CrossRef]
  79. Larsen, E.W. Solution of three dimensional inverse transport problems. Transp. Theory Stat. Phys. 1988, 17, 147–167. [Google Scholar] [CrossRef]
  80. Vos, M. Extracting detailed information from reflection electron energy loss spectra. J. Electron Spectrosc. Relat. Phenom. 2013, 191, 65–70. [Google Scholar] [CrossRef]
  81. Vos, M. A model dielectric function for low and very high momentum transfer. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. Atoms 2016, 366, 6–12. [Google Scholar] [CrossRef]
  82. Afanas’, V.P.; Gryazev, A.S.; Efremenko, D.S.; Kaplya, P.S. Differential inverse inelastic mean free path and differential surface excitation probability retrieval from electron energy loss spectra. Vacuum 2017, 136, 146–155. [Google Scholar] [CrossRef] [Green Version]
  83. Michaud, M.; Sanche, L. Absolute vibrational excitation cross sections for slow-electron (1–18 eV) scattering in solid H2O. Phys. Rev. A 1987, 36, 4684. [Google Scholar] [CrossRef] [PubMed]
  84. Michaud, M.; Sanche, L. Total cross sections for slow-electron (1–20 eV) scattering in solid H2O. Phys. Rev. A 1987, 36, 4672. [Google Scholar] [CrossRef]
  85. Wojcik, M.; Tachiya, M. Electron transport and electron–ion recombination in liquid argon: Simulation based on the Cohen–Lekner theory. Chem. Phys. Lett. 2002, 363, 381–388. [Google Scholar] [CrossRef]
  86. Winter, R.; Hensel, F.; Bodensteiner, T.; Gläser, W. The static structure factor of cesium over the whole liquid range up to the critical point. Berichte Bunsenges. Phys. Chem. 1987, 91, 1327–1330. [Google Scholar] [CrossRef]
  87. Greenfield, A.J.; Wellendorf, J.; Wiser, N. X-ray determination of the static structure factor of liquid Na and K. Phys. Rev. A 1971, 4, 1607. [Google Scholar] [CrossRef]
  88. Yarnell, J.L.; Katz, M.J.; Wenzel, R.G.; Koenig, S.H. Structure factor and radial distribution function for liquid argon at 85 K. Phys. Rev. A 1973, 7, 2130. [Google Scholar] [CrossRef]
  89. Garland, N.A.; Simonović, I.; Boyle, G.J.; Cocks, D.G.; Dujko, S.; White, R.D. Electron swarm and streamer transport across the gas–liquid interface: A comparative fluid model study. Plasma Sources Sci. Technol. 2018, 27, 105004. [Google Scholar] [CrossRef]
  90. Biagi, S.F. Biagi-v7.1 Database. 2004. Available online: https://us.lxcat.net/contributors/ (accessed on 28 February 2022).
  91. Misra, D. Mish: A self regularized non-monotonic activation function. arXiv 2019, arXiv:1908.08681. [Google Scholar]
  92. Innes, M. Flux: Elegant machine learning with Julia. J. Open Source Softw. 2018, 3, 602. [Google Scholar] [CrossRef] [Green Version]
  93. Glorot, X.; Bengio, Y. Understanding the difficulty of training deep feedforward neural networks. In Proceedings of the Thirteenth International Conference on Artificial Intelligence and Statistics, Sardinia, Italy, 13–15 May 2010; Volume 9, pp. 249–256. [Google Scholar]
  94. Zhuang, J.; Tang, T.; Ding, Y.; Tatikonda, S.C.; Dvornek, N.; Papademetris, X.; Duncan, J. AdaBelief Optimizer: Adapting Stepsizes by the Belief in Observed Gradients. In Advances in Neural Information Processing Systems; Larochelle, H., Ranzato, M., Hadsell, R., Balcan, M.F., Lin, H., Eds.; Curran Associates, Inc.: New York, NY, USA, 2020; Volume 33, pp. 18795–18806. [Google Scholar]
  95. Nesterov, Y. A method of solving a convex programming problem with convergence rate O1/k2. Sov. Math. Dokl. 1983, 27, 372–376. [Google Scholar]
  96. Dozat, T. Incorporating Nesterov Momentum into Adam. In Proceedings of the 4th International Conference on Learning Representations, San Juan, PR, USA, 2–4 May 2016. [Google Scholar]
  97. Ness, K.F.; Robson, R.E. Velocity distribution function and transport coefficients of electron swarms in gases. II. Moment equations and applications. Phys. Rev. A 1986, 34, 2185. [Google Scholar] [CrossRef] [PubMed]
  98. Boyle, G.J. The Modelling of Non-Equilibrium Light Lepton Transport in Gases and Liquids. Ph.D. Thesis, James Cook University, North Queensland, Australia, 2015. [Google Scholar]
  99. McCormick, N.J.; Kuščer, I. On the inverse problem in radiative transfer. J. Math. Phys. 1974, 15, 926–927. [Google Scholar] [CrossRef]
Figure 1. Diagram of crossed-beam electron scattering from a liquid micro-jet. The liquid is propelled into the page, while detectors are positioned normal to the liquid surface.
Figure 1. Diagram of crossed-beam electron scattering from a liquid micro-jet. The liquid is propelled into the page, while detectors are positioned normal to the liquid surface.
Ijms 23 03354 g001
Figure 2. An example EELS comparison between a 0.1 eV energy resolution (a) and a 0.9 eV energy resolution (b), for an initial electron energy of 50 eV scattering from neon [90]. Significant information is lost at 0.9 eV , especially when the peaks are in close proximity.
Figure 2. An example EELS comparison between a 0.1 eV energy resolution (a) and a 0.9 eV energy resolution (b), for an initial electron energy of 50 eV scattering from neon [90]. Significant information is lost at 0.9 eV , especially when the peaks are in close proximity.
Ijms 23 03354 g002
Figure 3. Neural network regression of Ne’s elastic MTCS [90], (a), and the corresponding agreement, (b), between its simulated EELS and the EELS used to perform the fit. The vertical lines in (b) denote the initial energies for each of the eleven EELS considered (from 10 eV to 100 eV incident energy, inclusive). Additionally, a shaded area in (a), along with vertical error bars in (b), are provided to indicate the range of the 100 best regressions. See also the figure legend for further details, noting that black is used to represent each colour of that marker shape and where colours correspond to each initial electron energies of 100, 79.43, 63.1, 50.12, 39.81, 31.62, 25.12, 19.95, 15.85, 12.59 and 10.0 eV , indicated by vertical lines. In this model, a bias vector size of 64 was used, along with 250,000 iterations.
Figure 3. Neural network regression of Ne’s elastic MTCS [90], (a), and the corresponding agreement, (b), between its simulated EELS and the EELS used to perform the fit. The vertical lines in (b) denote the initial energies for each of the eleven EELS considered (from 10 eV to 100 eV incident energy, inclusive). Additionally, a shaded area in (a), along with vertical error bars in (b), are provided to indicate the range of the 100 best regressions. See also the figure legend for further details, noting that black is used to represent each colour of that marker shape and where colours correspond to each initial electron energies of 100, 79.43, 63.1, 50.12, 39.81, 31.62, 25.12, 19.95, 15.85, 12.59 and 10.0 eV , indicated by vertical lines. In this model, a bias vector size of 64 was used, along with 250,000 iterations.
Ijms 23 03354 g003
Figure 4. Simultaneous neural network regression of Ne’s elastic MTCS and ionisation cross-sections [90], (a), and the corresponding agreement, (b), between its simulated EELS and the EELS used to perform the fit. The vertical lines in (b) denote the initial energies for each of the eleven EELS considered (from 10 eV to 100 eV incident energy, inclusive). Additionally, a shaded area in (a), along with vertical error bars in (b), are provided to indicate the range of the 100 best regressions. See also the figure legends for further details, noting that black is used to represent each colour of that marker shape, where colours correspond to each initial electron energy of 100, 79.43, 63.1, 50.12, 39.81, 31.62, 25.12, 19.95, 15.85, 12.59 , and 10.0 eV , indicated by vertical lines. In this model, a bias vector size of 128 was used along with 300,000 iterations.
Figure 4. Simultaneous neural network regression of Ne’s elastic MTCS and ionisation cross-sections [90], (a), and the corresponding agreement, (b), between its simulated EELS and the EELS used to perform the fit. The vertical lines in (b) denote the initial energies for each of the eleven EELS considered (from 10 eV to 100 eV incident energy, inclusive). Additionally, a shaded area in (a), along with vertical error bars in (b), are provided to indicate the range of the 100 best regressions. See also the figure legends for further details, noting that black is used to represent each colour of that marker shape, where colours correspond to each initial electron energy of 100, 79.43, 63.1, 50.12, 39.81, 31.62, 25.12, 19.95, 15.85, 12.59 , and 10.0 eV , indicated by vertical lines. In this model, a bias vector size of 128 was used along with 300,000 iterations.
Ijms 23 03354 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Muccignat, D.L.; Stokes, P.W.; Cocks, D.G.; Gascooke, J.R.; Jones, D.B.; Brunger, M.J.; White, R.D. Simulating the Feasibility of Using Liquid Micro-Jets for Determining Electron–Liquid Scattering Cross-Sections. Int. J. Mol. Sci. 2022, 23, 3354. https://doi.org/10.3390/ijms23063354

AMA Style

Muccignat DL, Stokes PW, Cocks DG, Gascooke JR, Jones DB, Brunger MJ, White RD. Simulating the Feasibility of Using Liquid Micro-Jets for Determining Electron–Liquid Scattering Cross-Sections. International Journal of Molecular Sciences. 2022; 23(6):3354. https://doi.org/10.3390/ijms23063354

Chicago/Turabian Style

Muccignat, Dale L., Peter W. Stokes, Daniel G. Cocks, Jason R. Gascooke, Darryl B. Jones, Michael J. Brunger, and Ronald D. White. 2022. "Simulating the Feasibility of Using Liquid Micro-Jets for Determining Electron–Liquid Scattering Cross-Sections" International Journal of Molecular Sciences 23, no. 6: 3354. https://doi.org/10.3390/ijms23063354

APA Style

Muccignat, D. L., Stokes, P. W., Cocks, D. G., Gascooke, J. R., Jones, D. B., Brunger, M. J., & White, R. D. (2022). Simulating the Feasibility of Using Liquid Micro-Jets for Determining Electron–Liquid Scattering Cross-Sections. International Journal of Molecular Sciences, 23(6), 3354. https://doi.org/10.3390/ijms23063354

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop