Development of Machine Learned Interatomic Potentials for Modeling Plasma-Material Interactions
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Abstract not provided.
Final performance report for our DOE grant
Self-focusing instability is a well-known phenomenon of nonlinear optics, which is of great importance in the field of laser–plasma interactions. Self-focusing instability leads to beam focusing and, consequently, breakup into multiple laser filaments. The majority of applications tend to avoid a laser filamentation regime due to its detrimental role on laser spot profile and peak intensity. In our work, using nonlinear Schrödinger equation solver and particle-in-cell simulations, we address the problem of interaction of multiple parallel beams in plasmas. We consider both non-relativistic and moderately relativistic regimes and demonstrate how the physics of parallel beam interaction transitions from the familiar self- and mutual-focusing instabilities in the non-relativistic regime to a moderately relativistic regime, where an analytical description of filament interaction is not available.
Reliable modeling of laser-plasma accelerators, where a short and intense laser pulse propagates in an underdense plasma over long distances, is a computationally challenging task. This is due to the great disparity among the scales involved in the modeling, ranging from the micrometer scale of the laser wavelength to, for instance, the meter scale of the laser-plasma interaction length for a multi-GeV-class laser-plasma accelerator. To reduce such imbalance, the time-averaged ponderomotive approximation may be used, where the plasma particle dynamics is analytically averaged over the laser frequency, and only spatiotemporal scales associated with the laser envelope are retained in the calculations, resulting in significant computational savings. Here, we characterize the accuracy and robustness of the time-averaged ponderomotive approximation for a range of laser parameters of interest for present and future laser-plasma accelerators, and we show that the error introduced by the averaging process is small in all relevant cases.
Laser–plasma interaction instabilities can be detrimental for direct-drive inertial confinement fusion by generating high-energy electrons that preheat the target. Moreover, an experimental platform has been developed and fielded on the National Ignition Facility to investigate hot-electron production from laser–plasma instabilities at direct-drive ignition-relevant conditions. The radiation-hydrodynamic code DRACO has been used to design planar-target experiments that generate plasma and interaction conditions comparable to direct-drive ignition designs: I L ~10 15 W/cm 2 , T e > 3 keV, and density-gradient scale lengths of L n ~ 600 μ m in the quarter-critical density region. The hot-electron properties were inferred by comparing the experimentally observed hard x-ray spectra to Monte Carlo simulations of hard x-ray emission from hot electrons depositing energy in the target. Hot-electron temperatures of ~40 keV to 60 keV and the fraction of laser energy converted to hot electrons of ~0.5% to 5% were inferred in plastic targets for laser intensities at the quarter-critical density surface of (~4 to 14) × 10 14 W/cm 2 . The use of silicon ablators was found to mitigate the hot-electron preheat by increasing the threshold laser intensity for hot-electron generation from ~3.5 × 10 14 W/cm 2 in plastic to ~6 × 10 14 W/cm 2 in silicon. The overall hot-electron production is further reduced in silicon ablators when the intensity threshold is exceeded.
Ultra-high intensity laser–plasma interactions can produce ultra-relativistic electrons via direct laser acceleration, assisted by quasi-static plasma magnetic and electric fields. These fields transversely confine electron motion and induce betatron oscillations. The net energy gain is strongly influenced by the interplay between two frequencies: the betatron frequency and the frequency of laser field oscillations experienced by the electron. Prior work has shown that energy gain is enabled by a resonance between the betatron oscillations and the oscillations of the laser field. In particular, higher-order resonances occur when the laser field completes multiple cycles during one betatron oscillation, allowing additional regimes of energy transfer beyond the fundamental (betatron) resonance. In this work, we demonstrate that such resonances become particularly effective when the laser's phase velocity is superluminal. Although the two frequencies generally evolve differently with increasing electron energy—leading to detuning—a superluminal phase velocity introduces a non-monotonic frequency ratio with a global minimum. This minimum allows sustained frequency matching over a broad energy range, thereby enabling enhanced energy gain. As the phase velocity increases, the betatron resonance becomes ineffective due to premature frequency detuning. At the same time, higher-order resonances become increasingly effective, emerging as the dominant mechanisms for enhanced energy gain in this regime of direct laser acceleration.
Here, an extended thermal particle boundary condition is devised to more efficiently and accurately model laser–plasma interactions in overdense plasmas. Particle-in-cell simulations of such interactions require many particles per cell, and a large region of background plasma is often necessary to correctly mimic a semi-infinite plasma and avoid electron refluxing from a truncated plasma. For long-pulse lasers of many picoseconds, such constraints can become prohibitively expensive. Here, an extended particle boundary condition (absorber) is designed that instantaneously stops and re-emits energetic particles streaming toward the simulation boundary over a defined region, allowing sufficient time and space for a suitably cool return current to develop in the background plasma. Tunable parameters of the absorber are explained, and simulations using the absorber with a 3-ps laser are shown to accurately reproduce those of a causally separated boundary while requiring only 20% the number of particles.
Plasma treatment of per- and polyfluoroalkyl substances (PFAS) contaminated water is a potentially energy efficient remediation method. In this treatment, an atmospheric pressure plasma interacts with surface-resident PFAS molecules. Developing a reaction mechanism and modeling of plasma–PFAS interactions requires fundamental data for electron–molecule reactions. In this paper, we present results of electron scattering calculations, potential energy landscapes and their implications for plasma modelling of a dielectric barrier discharge in PFAS contaminated gases, a first step towards modelling of plasma–water–PFAS intereactions. It is found that the plasma degradation of PFAS is dominated by dissociative electron attachment with the importance of other contributing processes varying depending on the molecule. All molecules posses a large number of shape resonances – transient negative ion states – from near-threshold up to ionization threshold. These states lie in the region of the most probable electron energies in the plasma (4–5 eV) and consequently are expected to further enhance the fragmentation dynamics in both dissociative attachment and dissociative excitation.
Since the invention of chirped pulse amplification, which was recognized by a Nobel Prize in physics in 2018, there has been a continuing increase in available laser intensity. Combined with advances in our understanding of the kinetics of relativistic plasma, studies of laser-plasma interactions are entering a new regime where the physics of relativistic plasmas is strongly affected by strong-field quantum electrodynamics (QED) processes, including hard photon emission and electron-positron (e⁻-e⁺) pair production. Additionally, this coupling of quantum emission processes and relativistic collective particle dynamics can result in dramatically new plasma physics phenomena, such as the generation of dense e⁻-e⁺ pair plasma from near vacuum, complete laser energy absorption by QED processes, or the stopping of an ultra-relativistic electron beam, which could penetrate a cm of lead, by a hair's breadth of laser light. In addition to being of fundamental interest, it is crucial to study this new regime to understand the next generation of ultra-high intensity laser-matter experiments and their resulting applications, such as high energy ion, electron, positron, and photon sources for fundamental physics studies, medical radiotherapy, and next generation radiography for homeland security and industry.
In this LDRD, we made significant advances on multiple fronts on the investigation and design of plasma-based optical structures to manipulate light at extreme intensities. We have unveiled new fundamental properties of laser-plasma interaction that were leveraged to introduce transformative applications. We first investigated new ways to generate high-amplitude plasma-based optical structures, based on the auto-resonance mechanism in order to produce refractive structures in plasmas. New insights into the interaction of such plasma structures with light waves allowed us to propose new methods for polarization manipulation and control for high-intensity lasers in the nonlinear regime, i.e. allowing the manipulation of an intense laser via its interaction with a less intense “control” laser beam in plasmas. Another application that we demonstrated experimentally is the capability to manipulate the velocity of light waves in plasmas by tailoring its refractive index; we reported the first demonstration of “slow light” and supra-luminal light propagation in plasmas. We also developed a new technological concept that would allow the probing of the type of plasma structures investigated in such applications, using high-bandwidth laser probing [Ludwig Phys. Plasmas 2019]. Finally, first steps were taken towards the demonstration of a plasma-based amplifier and compressor for the National Ignition Facility, in order to achieve powers equal of greater than those delivered by the Advanced Radiographic Capability (ARC). This technique uses a plasma-based grating acting like a compressor and amplifier, which could potentially compress and amplify an 88 ps, 56 J beam from the National Ignition Facility to <10 ps and hundreds of J.
The three-dimensional (3D) quasi-static particle-in-cell (PIC) algorithm is a very efficient method for modeling short-pulse laser or relativistic charged particle beam–plasma interactions. In this algorithm, the plasma response, i.e., plasma wave wake, to a non-evolving laser or particle beam is calculated using a set of Maxwell’s equations based on the quasi-static approximate equations that exclude radiation. The plasma fields are then used to advance the laser or beam forward using a large time step. The algorithm is many orders of magnitude faster than a 3D fully explicit relativistic electromagnetic PIC algorithm. It has been shown to be capable to accurately model the evolution of lasers and particle beams in a variety of scenarios. Additionally, at the same time, an algorithm in which the fields, currents and Maxwell equations are decomposed into azimuthal harmonics has been shown to reduce the algorithmic complexity of a 3D explicit PIC algorithm to that of a 2D algorithm when the expansion is truncated while maintaining accuracy for problems with near azimuthal symmetry. This hybrid algorithm uses a PIC description in r–z and a gridless description in . We describe a novel method that combines the quasi-static and hybrid PIC methods. This algorithm expands the fields, charge and current density into azimuthal harmonics. A set of the quasi-static field equations is derived for each harmonic. The complex amplitudes of the fields are then solved using the finite difference method. The beam and plasma particles are advanced in Cartesian coordinates using the total fields. Details on how this algorithm was implemented using a similar workflow to an existing quasi-static code, QuickPIC, are presented. The new code is called QPAD for QuickPIC with Azimuthal Decomposition. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), a full 3D quasi-static code (QuickPIC), and the new quasi-static PIC code with azimuthal decomposition (QPAD) are also presented.
Next step fusion devices that will operate in steady-state will require complex plasma-facing components (PFCs) that can survive the harsh environment over long timescales not common in current devices. This will require robust plasma facing surfaces that are integrated with active cooling systems. In a collaboration between CEA and ORNL, a plasma-interacting diagnostic is being designed for the W Environment in Steady-state Tokamak (WEST) in Cadarache, France which requires plasma-facing protection like those needed for steady-state PFCs. This integrated diagnostic studies edge transport and impurity migration within WEST, and is planned to include imbedded temperature Langmuir probe sensors, as well as removeable sample slots for ex-situ surface analysis of plasma-material interactions. The entire assembly is expected to move into the plasma edge for periods up to 1000 s having an energy removeable capability of ∼6 kW and seeing a peak heat flux of more than 7 MW/m2. The assembly compliments WEST high fluence campaigns that plan for multiple 1000 s pulses. Within these specifications, the assembly will require a refractory metal plasma facing surface and integral cooling in order to function within the limited space allotted for such diagnostics. Because of the limited space and linear actuator needs, additive manufacturing of the high heat flux working end of the assembly is being considered which could allow for precision cooling-channels and lighter weight designs. Conceptual design along with simulation and analysis results will be presented for this complex diagnostic with novel PFCs.
Understanding how atoms interact with hot dense matter is essential for astrophysical and laboratory plasmas. Interactions in high-density plasmas broaden spectral lines, providing a rare window into interactions that govern, for example, radiation transport in stars. However, up to now, spectral line-shape theories employed at least one of three common approximations: second-order Taylor treatment of broadening operator, dipole-only interactions between atom and plasma, and classical treatment of perturbing electrons. In this Letter, we remove all three approximations simultaneously for the first time and test the importance for two applications: neutral hydrogen and highly ionized magnesium and oxygen. Finally, we found 15%–50% change in the spectral line widths, which are sufficient to impact applications including white-dwarf mass determination, stellar-opacity research, and laboratory plasma diagnostics.
A third generation flowing liquid lithium (FLiLi) limiter with a substrate made of TZM, an alloy with >99% Mo, was fabricated by conventional manufacturing techniques. TZM has a high corrosion resistance, a high sputtering threshold, and a good wettability to Li, as compared to stainless steel (SS), which had been used as an FLiLi substrate surface in 2014 and 2016. The third generation FLiLi was inserted into the edge in EAST H-mode plasmas in an upper single-null configuration with an ion grad-B drift toward the upper divertor with a limiter temperature of 330–380 °C and an auxiliary heating power of about 2–8 MW. Analysis has shown that by using TZM FLiLi, fuel particle recycling continuously decreased and near-complete edge localized mode elimination was achieved in H-mode plasmas with RF-only heating. The main impurities during the initial FLiLi discharges were Mo, Fe, and W resulting from strong plasma interaction at the Mo plate side, the SS collector, and the upper W divertor. Plasma stored energy increased by about 10 kJ in subsequent FLiLi discharges due to decreased impurity radiation. Engineering analysis shows a uniform lithium flow with an ~80% Li coverage ratio on the limiter surface, similar to the second FLiLi, even though only one of two J × B pumps was functioning. Despite technical difficulties, the FLiLi gen.3 improved the overall plasma performance, providing support for flowing liquid Li plasma facing component applications in present and future devices.
Magnetized plasma interactions are ubiquitous in astrophysical and laboratory plasmas. Various physical effects have been shown to be important within colliding plasma flows influenced by opposing magnetic fields, however, experimental verification of the mechanisms within the interaction region has remained elusive. Here we discuss a laser-plasma experiment whereby experimental results verify that Biermann battery generated magnetic fields are advected by Nernst flows and anisotropic pressure effects dominate these flows in a reconnection region. These fields are mapped using time-resolved proton probing in multiple directions. Various experimental, modelling and analytical techniques demonstrate the importance of anisotropic pressure in semi-collisional, high- β plasmas, causing a reduction in the magnitude of the reconnecting fields when compared to resistive processes. Anisotropic pressure dynamics are crucial in collisionless plasmas, but are often neglected in collisional plasmas. We show pressure anisotropy to be essential in maintaining the interaction layer, redistributing magnetic fields even for semi-collisional, high energy density physics (HEDP) regimes.
This project investigated how plasmas interact with energetic particles and solid materials under extreme conditions relevant to fusion energy, space plasmas, and planetary environments. Using experiments on the DIII-D National Fusion Facility, the research first examined how high-energy electrons move, become trapped, and are released in plasmas containing magnetic islands—structures commonly found in fusion reactors and Earth’s magnetosphere—providing new insight into particle transport and acceleration processes. The project also explored plasma-driven chemical reactions that can occur during meteoroid entry into planetary atmospheres, demonstrating that simple molecules such as ammonia can be produced and survive in high-temperature plasma conditions. Together, these results improve understanding of plasma behavior across laboratory, space, and planetary systems while informing fusion plasma control and plasma–material interaction studies. The project additionally contributed to workforce development by training graduate students, undergraduates, and early-career researchers and by disseminating results through peer-reviewed publications and international scientific conferences.
We explore the applications of machine learning techniques in relativistic laser-plasma experiments beyond optimization purposes. We predict the beam charge of electrons produced in a laser wakefield accelerator given the laser wavefront change caused by a deformable mirror. Machine learning enables feature analysis beyond merely searching for an optimal beam charge, showing that specific aberrations in the laser wavefront are favored in generating higher beam charges. Supervised learning models allow characterizing the measured data quality as well as recognizing irreproducible data and potential outliers. Furthermore, we also include virtual measurement errors in the experimental data to examine the model robustness under these conditions. This work demonstrates how machine learning methods can benefit data analysis and physics interpretation in a highly nonlinear problem of relativistic laser-plasma interaction.