SEARCH · Search NASA

Results for “bootstrap”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Efficient ECCD non-inductive plasma current start-up, ramp-up, and sustainment for an ST fusion reactor

The elimination of the need for an Ohmic heating solenoid may be the most impactful design driver for the realization of economical compact fusion tokamak reactor systems. However, this would require fully non-inductive start-up and current ramp-up from zero plasma current and low electron temperature of sub-keV to the full plasma current of ~10–15 MA at 20–30 keV electron temperature. To address this challenge, an efficient solenoid-free start-up and ramp-up scenario utilizing a low-field-side-launched extraordinary mode at the fundamental electron cyclotron harmonic frequency (X–I) is proposed, which has more than two orders of magnitude higher electron cyclotron current drive (ECCD) efficiency than the conventional ECCD for the sub-keV start-up regime. A time dependent model was developed to simulate the start-up scenarios. For the Spherical Tokamak Advanced Reactor (STAR) (Menard et al 2023 Next-Step Low-Aspect-Ratio Tokamak Design Studies (IAEA)), it was found that to fully non-inductively ramp-up to 15 MA, it would take about 25 MW of EC power at 170 GHz. Because of the relatively large plasma volume of STAR, radiation losses must be considered. It is important to make sure that high Z impurities are kept sufficiently low during the early current start-up phase where the temperature is sub-keV range. Since the initial current ramp up takes place at a factor of ten lower density compared to the sustained regimes, it is important to transition into a higher bootstrap fraction discharge at lower density to minimize the ECCD power requirement during the densification. For the sustainment phase an array of eight gyrotron launchers with a total of about 60 MW of fundamental O-mode was found to be sufficient to provide the required axis-peaked external current drive. High efficiencies between 19–57 kA MW –1 were found with optimal aiming, and these were resilient to small changes in aiming angles and density and temperature profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Variable-spectrum mode control of high poloidal beta discharges

DIII-D experiments demonstrate that high pressure, broad current profile equilibria can be accessed in the high poloidal beta regime by optimizing the MHD mode control poloidal spectrum. A novel, variable spectrum (VS) magnetic feedback scheme implemented using the DIII-D internal non-axisymmetric coils (I-coils) facilitated access to reduced internal inductance $l$ i operation above the no-wall beta limit compared with both no feedback and fixed spectrum feedback. In addition, the VS feedback helped avoid beta collapses caused by marginally unstable resistive wall mode activity. The lower and upper I-coil rows were configured in two independent feedback loops, allowing the feedback field's poloidal spectrum to vary and track changes in the plasma mode structure as the edge safety factor q 95 varied from 11 to 6 during the discharges. The q 95 dependence of the measured phase difference between the lower and upper I-coil rows during VS feedback is qualitatively compatible with ideal MHD simulations of the least-stable plasma kink mode and with plasma response simulations that included kinetic modifications to ideal MHD. The VS feedback approach is a straightforward way to improve resilience to variations in mode structure that occur as plasma parameters change. The demonstrated expansion of the operating space to lower $l$ i is expected to improve the coupling of the plasma kink mode to external fields and beneficial wall eddy currents, and is compatible with high bootstrap fraction operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Impurity radiation seeding of neoclassical tearing mode growth

Abstract The physics of neoclassical tearing modes (NTMs) is of great concern to tokamak plasma stability and performance, especially in the burning plasma regime. Whereas in many situations the seed events can be clearly identified, such as sawteeth and edge localized modes, the potential seeding mechanism of NTMs due to the resistive tearing instability driven by the impurity radiation cooling still needs more study. Recent NIMROD simulations have demonstrated that local impurity radiation cooling can drive the seed island growth and trigger the subsequent onset of NTM instability. The seed island is mainly driven by the local helical perturbation of the diamagnetic current induced by the perturbed pressure gradient as a result of the impurity radiative cooling on the rational surface. A heuristic closure for the neoclassical viscosity is adopted, and the seed island is further driven by the perturbed bootstrap current induced from the neoclassical electron viscous stress in the extended Ohm’s law. The growth rate of the NTM in simulations is found proportional to the electron neoclassical viscosity, and a theoretical neoclassical driving term is adopted to account for the nonlinear neoclassical island growth in the simulations.

Zeng, Shiyong (ORCID:000000017098818X)↗

DOE OSTI2025

Identification of a window for quiescent H-mode operation in MHD stability diagram of DIII-D plasmas

A window for quiescent H-mode (QH-mode) operation in DIII-D plasmas in the diagram of pedestal magnetohydrodynamic (MHD) stability was identified through linear MHD stability analysis, accounting for the effects of plasma rotation and ion diamagnetic drift. The operation window lies between the stability boundaries of the kink/peeling mode (K/PM) identified with and without the inclusion of plasma rotation effects. In this region, the mode remains unstable unless rotation effects are considered alongside the ion diamagnetic drift, which is consistently included in the analysis. The stabilization of the MHD mode, facilitated by the coupled effects of plasma rotation and ion diamagnetic drift, plays a crucial role in widening the window, enabling the attainment of the QH-mode state. Furthermore, the suppression of edge localized modes (ELM) can be achieved by controlling the pedestal structure to ensure the plasma state remains within the operation window. The location of the operation window in the stability diagram depends on the K/PM stability properties. Therefore, optimizing conditions for QH-mode requires adjustments based on changes in stability characteristics. A pressure pedestal and its associated bootstrap current density near the last closed flux surface are advantageous for situating the plasma state within the window. However, excessive current density can trigger ELMs. This trend was confirmed through comparisons of MHD stability diagrams between QH-mode and ELMy H-mode plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Reconstruction of current profiles from external magnetic and MSE measurements on DIII-D

This work introduces a new current model into the EFIT code for plasma current profile reconstruction. Using this model, edge-localized current density profiles are reconstructed solely from external magnetic measurements. Quantitative comparisons show that edge current profiles obtained with the new current model align closely with those from kinetic equilibrium reconstruction. The new model enables current density reconstruction with fewer diagnostics than kinetic equilibrium reconstruction (Lao et al 1990 Nucl. Fusion 30 1035–49), demonstrating that edge-localized current density information can be obtained from external magnetic measurements only, without the help of kinetic profiles and internal diagnostics, and without using the bootstrap current model. Furthermore, with the constraints of external magnetic and MSE measurements, current profiles in full radius obtained from the new model are in good agreement with the results of kinetic equilibrium reconstruction, which not only improves the efficiency of current reconstruction but also helps avoid possible errors caused by different kinds of diagnoses. Possible applications are also discussed, and a brief analysis of the evolution of the full-radius current density is performed. Noticeable decreases in both the magnitude and gradient of the edge current around the time of ELM bursts are observed. Moreover, analysis of the full radius current density profiles ( ρ N = 0 − 1.0, ρN denotes to square root of the normalized toroidal flux) reveals that as the edge current changes, the inner current also changes following a sinusoidal-like pattern at almost all time slices during four type-I ELMs in discharge 161493. The current evolution in the inner region ( ρ N = 0.2 − 0.6) is synchronized with the boundary region ( ρ N = 0.89 − 1.0), while the intermediate region ( ρ N = 0.6 − 0.89) shows opposite behavior. This sinusoidal-like variations in full-radius current density profile indicates strong core-pedestal coupling, suggesting ELMs significantly affect core plasmas, though current changes manifest differently across radial regions.

EFIT↗

DOE OSTI2026

Spherical tokamak physics research in preparation for the operation of NSTX-U

The National Spherical Torus Experiment Upgrade (NSTX-U) is preparing to resume operation, representing a crucial step toward realizing compact, cost-effective fusion pilot plants. In advance of this, extensive modeling and data analysis have been conducted to advance the physics basis for low-aspect-ratio, high-performance plasma regimes, focusing on three core objectives: confinement and stability, power and particle handling, and steady-state operation. Significant progress has been made in understanding the electron temperature flattening in high-β plasmas, which is shown to be driven by a complex interplay of magnetohydrodynamic instabilities (e.g. non-resonant infernal modes), fast-ion-driven Alfvén eigenmodes, and electron and ion-scale micro-instabilities, particularly Kinetic Ballooning Modes (KBMs), whose destabilization is strongly dependent on parallel magnetic field fluctuations (δB ∥ ). Furthermore, a new gyrokinetic critical pedestal model was developed, accurately predicting pedestal structure by identifying KBMs as the primary stability limit, offering a critical constraint for future high-confinement scenarios. To address the challenge of high heat flux, novel liquid lithium plasma-facing components were modeled. The analysis confirmed that lithium vapor shielding is a self-regulating mechanism for heat mitigation, while also emphasizing that strong main ion parallel flow is essential to minimize core lithium contamination. Finally, progress toward steady-state operation was anchored by developing the required physics basis and control tools. This includes predictive modeling for reversed magnetic shear sustainment, demonstrating that magnetic island-induced bootstrap current reduction is negligible in STs, and advancing real-time control and disruption avoidance capabilities. The development of high-speed surrogate models (e.g. MMMNet) provides computationally efficient tools vital for non-inductive scenario optimization and integrated, low-disruptivity operations planned for NSTX-U.

NSTX-U↗

DOE OSTI2026

Macroscopic trends of neoclassical tearing stability in high-field H-mode tokamak pilot plants

The neoclassical tearing mode (NTM) stability metric—minimum marginally stable island width $w$$^{*}_{m}$—was compared across 14651 inductive high-field tokamak pilot plant equilibria. Larger devices with reduced elongation and/or increased minor radius demonstrated an order-of-magnitude increase in $w$$^{*}_{m}$, primarily due to a reduction in bootstrap drive. This work is part of an ongoing effort to ensure passive NTM-stability in the ARC tokamak, in which the technology to achieve active tearing-suppression with localised electron cyclotron current drive does not yet exist. The equilibrium scenarios in the database were Monte Carlo generated and normalised to the same >400MW fusion power, minimum pressure scenario at a range of plasma currents, before tearing analysis using the modified Rutherford equation was applied for all resonant poloidal and toroidal m, n modes up to n = 4. Single-helicity toroidal Δ' calculations in resistive DCON set the minimum marginally stable island width, and a simple modal scaling proportional to –m 2 n –1 was identified for high-m Δ' values. The dominant correlates of $w$$^{*}_{m}$ and Δ' across the database were analysed using interpretable machine learning techniques.

NTM seeding↗

DOE OSTI2026

Auxiliary heating and current drive physics for the ST-E1 fusion power plant

This work describes the physics basis for the proposed auxiliary heating and current drive system on the ST-E1 fusion power plant. The ST-E1 flattop plasma considered here is fully non-inductive with a bootstrap fraction of 0.9 and the remaining current driven by EC waves. Using the recently published physics-based optimization method for EC launchers (Lopez et al 2025 Plasma Phys. Control. Fusion 67 055012), we show that the target flattop ECCD can be achieved with a net efficiency of 52 kA MW −1 using fundamental O-mode (O1) with frequency range 160–200 GHz launched from the low-field side top half of the vacuum vessel (LFS top-launch). From considering two candidate rampup scenarios, we conclude that LFS top-launch O1 ECCD can be equally effective during the early stages of plasma operation, although poloidal steering might be needed. X-mode waves injected from the LFS midplane are also shown to be effective for rampup even when T e < 1 keV. We also present modeling results for the pre-conceptual design of an ICRH system proposed for ST-E1. Using TORIC, we find that an ICRH system aiming for 42–48 MHz and toroidal mode number n φ ~ 10 robustly achieves dominant ion damping via Helium-3 minority heating transitioning to second-harmonic Tritium heating. We then show that such waves can be efficiently generated by a 5-strap traveling-wave antenna (TWA) using the Petra-M code. The TWA has a 40–45 MHz passband within which ~60% of the power entering the TWA is coupled to the plasma with the remaining ~40% of the power being transmitted through the TWA and possibly recirculated; the power reflected back into the transmission lines is negligible. This passband structure persists even when the evanescent distance is increased by a factor of two, or when the magnetic-field angle is increased by 30°, demonstrating inherent load resilience that will be crucial for effective ICRH on ST-E1.

electron cyclotron↗

DOE OSTI2026

Achievement of a high-density, high-confinement, and high-beta tokamak plasma regime in DIII-D, and implications for a lower-current path for ITER and FPP

Experiments on DIII-D have demonstrated a density-confinement synergy that enables sustainment of high performance in a previously unattained parameter regime of simultaneous very high energy confinement quality (H 98y2 ≥ 1.5), very high line-average density Greenwald fraction (ƒ Gr = πa 2 < n >/I P ≥ 1.4), and high toroidal beta (β T ≥ 3%). Tokamak operation in this regime is essential for a compact steady-state FPP, as well as for Q=10 with 500 MW of fusion power in ITER at I P << 15 MA. These experiments leveraged the knowledge that, in the high-poloidal-beta (β P ) regime, impurity and density gradients can enhance turbulence stabilization caused by high α MHD (α MHD ~(dβ P )⁄dr). This was described by theoretical predictions and gyrokinetic transport simulations [M.T. Kotschenreuther et al, 2024 Nucl. Fusion, 64 076033], and later confirmed by experiments on DIII-D [S. Ding et al, 2024 Nature 629 555]. To increase both β P and β T , the new experiments increased the ideal-wall stability β N -limit by using a smaller plasma-outer wall distance and higher triangularity in the plasma cross section (top/bottom average δ~0.9), enabled by the recent “shape & volume rise” (SVR) modification to the DIII-D divertor. The higher triangularity also contributed to achieving higher ƒ Gr by enabling higher pedestal density. At high density, the pedestal is ballooning limited and exhibits small and frequent ELMs, while the divertor is near detachment even without any impurity seeding. High plasma performance was attained and sustained reproducibly, with the eventual terminations brought about by an MHD mode destabilized as the current profile slowly continued to evolve. A path to stationary fully noninductive operation might include ECH injection to reduce both core impurity accumulation and the electron collisionality, thus increasing the bootstrap current. These experiments provide the first experimental demonstration of the ƒ Gr , H 98y2 , and β T values required simultaneously for ITER Q = 10 at I P < 10 MA, pointing to practical ways to improve the energy confinement in a fusion reactor.

Garofalo, Andrea M. [General Atomics, San Diego, C↗

DOE OSTI2025

Unraveling the depth-dependent causal dynamics of methanogenesis and methanotrophy in a high-latitude fen peatland

The dynamics of methane (CH 4 ) cycling in high-latitude peatlands through different pathways of methanogenesis and methanotrophy are still poorly understood due to the spatiotemporal complexity of microbial activities and biogeochemical processes. Additionally, long-term in situ measurements within soil columns are limited and associated with large uncertainties in microbial substrates (e.g. dissolved organic carbon, acetate, hydrogen). To better understand CH 4 cycling dynamics, we first applied an advanced biogeochemical model, ecosys , to explicitly simulate methanogenesis, methanotrophy, and CH 4 transport in a high-latitude fen (within the Stordalen Mire, northern Sweden). Next, to explore the vertical heterogeneity in CH 4 cycling, we applied the PCMCI/PCMCI+ causal detection framework with a bootstrap aggregation method to the modeling results, characterizing causal relationships among regulating factors (e.g. temperature, microbial biomass, soil substrate concentrations) through acetoclastic methanogenesis, hydrogenotrophic methanogenesis, and methanotrophy, across three depth intervals (0–10 cm, 10–20 cm, 20–30 cm). Our results indicate that temperature, microbial biomass, and methanogenesis and methanotrophy substrates exhibit significant vertical variations within the soil column. Soil temperature demonstrates strong causal relationships with both biomass and substrate concentrations at the shallower depth (0–10 cm), while these causal relationships decrease significantly at the deeper depth within the two methanogenesis pathways. In contrast, soil substrate concentrations show significantly greater causal relationships with depth, suggesting the substantial influence of substrates on CH 4 cycling. CH 4 production is found to peak in August, while CH 4 oxidation peaks predominantly in October, showing a lag response between production and oxidation. Overall, this research provides important insights into the causal mechanisms modulating CH 4 cycling across different depths, which will improve carbon cycling predictions, and guide the future field measurement strategies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Downscaled CMIP5 projections of physical fire risk understate historical trends

Reliable projections of wildfire risk are important for multi-sector impacts analysis. Statistically downscaled and bias-corrected Earth system model ensemble products are routinely used to analyze regional physical wildfire risk, but evaluations of historical observed trends and variability are lacking. Here, we evaluate physical fire risk over the western United States using the Canadian Forest Fire Weather Index (FWI) by comparing model outputs from the Coupled Model Intercomparison Project Phase 5 (CMIP5), statistically downscaled via the Multivariate Adaptive Constructed Analogs (MACA) approach, against the observational target dataset gridMET, a gridded high-resolution surface meteorological product. We analyze multidecadal trends and interannual variability in seasonal average FWI for the historical period and future projections under two emissions scenarios, and we compare MACA-CMIP5 ensemble results with a simple time series model that generates historical and future projections of seasonal FWI based on bootstrapping observed historical trends and variability. Our findings indicate that MACA-CMIP5 accurately captures the magnitude and spatial patterns of seasonally averaged FWI but tends to underestimate historical decadal trends. We show that future increases in fire risk may be underestimated relative to the simple time series model that projects historical variability into the future. We also highlight that model biases in relative humidity contribute significantly to model-data differences. Our results underscore the importance of historical hindcasting exercises for informing broader multi-sector applications.

FWI↗

DOE OSTI2020

A neural network for determination of latent dimensionality in Nonnegative Matrix Factorization

Non-negative Matrix Factorization (NMF) has proven to be a powerful unsupervised learning method for uncovering hidden features in complex and noisy datasets with applications in data mining, text recognition, dimension reduction, face recognition, anomaly detection, blind source separation, and many other fields. An important input for NMF is the latent dimensionality of the data, that is, the number of hidden features, K, present in the explored dataset. Unfortunately, and this quantity is rarely known a priori. The existing methods for determining latent dimensionality, such as Automatic Relevance Determination (ARD), are mostly heuristic and utilize different characteristics to estimate the number of hidden features. However, all of them require human presence to make a final determination of K. Here we utilize a supervised machine learning approach in combination with a recent method for model determination, called NMFk, to determine the number of hidden features automatically. NMFk performs a set of NMF simulations on an ensemble of matrices, obtained by bootstrapping the initial dataset, and estimates which K produces stable groups of latent features that reconstruct the initial dataset well. We then train a Multi-Layter Perceptron (MLP) classifier network to determine the correct number of latent features utilizing the statistics and characteristics of the NMF solution, obtained from NMFk. In order to train the MLP classifier, a training set of 58,660 matrices with predetermined latent features were factorized with NMFk. The MLP classifier in conjunction with NMFk maintains a greater than 95% success rate when applied to a held out test set. Additionally, when applied to two well-known benchmark datasets, the swimmer and MIT face data, NMFk/MLP correctly recovers the established number of hidden features. Finally, we compare the accuracy of our method to the ARD, AIC and Stability-based methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

A measurement of the Hubble constant from Type II supernovae

ABSTRACT Progressive increases in the precision of the Hubble-constant measurement via Cepheid-calibrated Type Ia supernovae (SNe Ia) have shown a discrepancy of ∼4.4σ with the current value inferred from Planck satellite measurements of the cosmic microwave background radiation and the standard $\Lambda $cold dark matter (ΛCDM) cosmological model. This disagreement does not appear to be due to known systematic errors and may therefore be hinting at new fundamental physics. Although all of the current techniques have their own merits, further improvement in constraining the Hubble constant requires the development of as many independent methods as possible. In this work, we use SNe II as standardisable candles to obtain an independent measurement of the Hubble constant. Using seven SNe II with host-galaxy distances measured from Cepheid variables or the tip of the red giant branch, we derive H$_0= 75.8^{+5.2}_{-4.9}$ km s−1 Mpc−1 (statistical errors only). Our value favours that obtained from the conventional distance ladder (Cepheids + SNe Ia) and exhibits a difference of 8.4 km s−1 Mpc−1 from the Planck + ΛCDM value. Adding an estimate of the systematic errors (2.8 km s−1 Mpc−1) changes the ∼1.7σ discrepancy with Planck +ΛCDM to ∼1.4σ. Including the systematic errors and performing a bootstrap simulation, we confirm that the local H0 value exceeds the value from the early Universe with a confidence level of 95 per cent. As in this work, we only exchange SNe II for SNe Ia to measure extragalactic distances, we demonstrate that there is no evidence that SNe Ia are the source of the H0 tension.

de Jaeger, T.↗

DOE OSTI2023

Saccharomycopsis praedatoria sp. nov., a predacious yeast isolated from soil and rotten wood in an Amazonian rainforest biome

Three yeast isolates were obtained from soil and rotting wood samples collected in an Amazonian rainforest biome in Brazil. Comparison of the intergenic spacer 5.8S region and the D1/D2 domains of the large subunit rRNA gene showed that the isolates represent a novel species of the genus Saccharomycopsis. A tree inferred from the D1/D2 sequences placed the novel species near a subclade containing Saccharomycopsis lassenensis, Saccharomycopsis fermentans, Saccharomycopsis javanensis, Saccharomycopsis babjevae, Saccharomycopsis schoenii and Saccharomycopsis oosterbeekiorum, but with low bootstrap support. In terms of sequence divergence, the novel species had the highest identity in the D1/D2 domains with Saccharomycopsis capsularis, from which it differed by 36 substitutions. In contrast, a phylogenomic analysis based on 1061 single-copy orthologs for a smaller set of Saccharomycopsis species whose whole genome sequences are available indicated that the novel species represented by strain UFMG-CM-Y6991 is phylogenetically closer to Saccharomycopsis fodiens and Saccharomycopsis sp. TF2021a (=Saccharomycopsis phalluae). The novel yeast is homothallic and produces asci with one spheroidal ascospore with an equatorial or subequatorial ledge. The name Saccharomycopsis praedatoria sp. nov. is proposed to accommodate the novel species. The holotype of Saccharomycopsis praedatoria is CBS 16589 T . The MycoBank number is MB849369. S. praedatoria was able to kill cells of Saccharomyces cerevisiae by means of penetration with infection pegs, a trait common to most species of Saccharomycopsis.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2026

Third moments of nucleon unpolarized, polarized, and transversity parton distribution functions from physical-point lattice QCD

Using forward matrix elements of local leading-twist operators, we present a determination of the isovector third Mellin moments ⟨𝑥 2 ⟩ of nucleon unpolarized, polarized, and transversity parton distribution functions. Two lattice QCD ensembles at the physical pion mass are used, which were generated using a tree-level Symanzik-improved gauge action and 2+1 flavor tree-level improved Wilson Clover fermions coupling via 2-level HEX-smearing. Leveraging a wide set of operators, two extraction methods for the matrix elements, and the automatic inclusion of model uncertainties via bootstrapped model averages, we extract values of the third Mellin moments. Furthermore, this is the first direct calculation of these observables performed at the physical pion mass.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Lanczos Algorithm, the Transfer Matrix, and the Signal-to-Noise Problem

This Letter introduces a method for determining the energy spectrum of lattice quantum chromodynamics by applying the Lanczos algorithm to the transfer matrix and using a bootstrap generalization of the Cullum-Willoughby method to filter out spurious eigenvalues. Proof-of-principle analyses of the simple harmonic oscillator and the lattice quantum chromodynamics proton mass demonstrate that this method provides faster ground-state convergence than the “effective mass,” which is related to the power-iteration algorithm. Lanczos provides more accurate energy estimates than multistate fits to correlation functions with small imaginary times while achieving comparable statistical precision. Two-sided error bounds are computed for Lanczos results and guarantee that excited-state effects cannot shift Lanczos results far outside their statistical uncertainties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

X-ray-line coincidence photopumping in a potassium-chlorine mixed plasma

Exploiting the multiple long pulse capability and suite of x-ray diagnostics of the Orion laser, we have set out to explore line coincidence photopuming—the enhancement in population of an atomic level brought on by resonant absorption of x rays from a different emitting ion. Unlike previous work, the two ions are in the same plasma and so the experiment is an x-ray analog of the well-known Bowen resonance fluorescence mechanism that operates in astrophysical situations in the optical region. Herein, our measurements have shown enhanced fluorescence in a chlorine plasma, attributable to line coincident photopumping from co-mixed potassium ions. To detect this relatively low signal-to-noise phenomenon, the data from multiple shots are combined, and the statistical method of bootstrapping is used to assign a confidence value to the measured enhancement, resulting in an estimate of the enhancement of 39 ± 18 16 % compared to the null case, where no pumping occurs. The experimental results have been compared to coupled radiation-transport and radiation hydrodynamics simulations using the CRETIN code together with the NYM radiation hydrodynamics model and agreement has been found, with the simulations also predicting modest enhancement.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2021

Uncertainty quantification for deep learning in particle accelerator applications

With the advent of increased computational resources and improved algorithms, machine learning-based models are being increasingly applied to complex problems in particle accelerators. However, such data-driven models may provide overly confident predictions with unknown errors and uncertainties. For reliable deployment of machine learning models in high-regret and safety-critical systems such as particle accelerators, estimates of prediction uncertainty are needed along with accurate point predictions. In this investigation, we evaluate Bayesian neural networks (BNN) as an approach that can provide accurate predictions along with reliably quantified uncertainties for particle accelerator problems, and compare their performance with bootstrapped ensembles of neural networks. We select three accelerator setups for this evaluation: a storage ring, a photoinjector, and a linac. The problems span different data volumes and dimensionalities (e.g., scalar predictions as well as image outputs). It is found that BNN provide accurate predictions of the mean along with reliable estimates of predictive uncertainty across the test cases. In this vein, BNN may offer an attractive alternative to deterministic deep learning tools to generate accurate predictions with quantified uncertainties in particle accelerator applications.

43 PARTICLE ACCELERATORS↗