Search NASA⌕ Search

SEARCH · Search NASA

Results for “Linear accelerators”

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.

At least 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Enhancing lifetime, forecasting, and economic benefits of photovoltaic technologies undergoing UV-induced degradation with optical filtering

Ultraviolet-induced degradation (UV-ID) of various PV cell types was analyzed under optical UV filters with different cutoff wavelengths. Cell types studied included interdigitated back contact (IBC), passivated emitter and rear totally diffused (PERT), and heterojunction technology (HJT) based on crystalline Si (c-Si), and metal halide perovskite (MHP) cells. Analyzing degradation rates in two distinct regimes proved beneficial for all cell types. We used empirical linearizing functions ln(t) for c-Si technologies and 2 √t for MHP samples where t is time. These were applied to extrapolate UV-induced degradation over the lifetime of PV modules under various levels of optical UV filtering and used to predict the relative economic benefits for PV power plants. Degradation rates for all technologies were generally faster under the long pass optical filters having shorter cutoff wavelengths transmitting more UV irradiation and at elevated temperatures when testing MHP samples in the range between 60 °C and 90 °C.

14 SOLAR ENERGY↗
DOE OSTI2026

Unpolarized GPDs at small x and non-zero skewness

We study the small-x asymptotics of unpolarized generalized parton distributions (GPDs) and generalized transverse momentum distributions (GTMDs). Unlike the previous works in the literature, we consider the case of non-zero (but small) skewness while allowing for non-linear contributions to the evolution equations. We first show that unpolarized GPDs and GTMDs at small x are related to the eikonal dipole amplitude N, whose small-x evolution is given by the BK/JIMWLK evolution equations, and to the odderon amplitude $\mathscr{O}$, whose evolution is also known in the literature. We then show that the effect of non-zero skewness ξ ≠ 0 is to modify the value of the evolution parameter (rapidity) in the arguments for the dipole amplitudes N and $\mathscr{O}$ from Y = ln(1/x) to Y = ln min{1/|x|, 1/|ξ|}.

Generalized parton distributions (GPDs)↗
DOE OSTI2024

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗
DOE OSTI2026

Magic Diamond: Covalent Bond Formation of Melamine and Other Amines on Nanodiamond Surfaces

High-temperature, high-pressure (HPHT) nanodiamond (ND) hosts nitrogen-vacancy (NV) centers, solid-state qubits that enable room-temperature quantum sensing by all-optical magnetometry, electrometry, and thermometry. However, the covalent surface functionalization of nanoscale diamond remains largely limited to carboxylate-based chemistries. Amine termination is particularly attractive because theoretical studies predict suppression of midgap states and extended electron-spin coherence times. Recently, chemical activation of alcohol-terminated NDs to alkyl bromides (ND-Br) using SOBr2 has enabled nucleophilic substitution through a carbocation intermediate, allowing formation of simple amine terminations. Here, we evaluate whether sterically demanding amines can form covalent diamond−nitrogen bonds on ND-Br surfaces. ND-Br was reacted with branched, linear, and cyclic amines, including polyethylenimine, diethylenetriamine, and melamine. X-ray spectroscopies were used to confirm successful and to probe the resulting electronic structure at the diamond−amine interface. These results expand the chemical toolbox for tuning diamond surface dipoles and electron affinity, providing new pathways for engineering nanodiamond surfaces for quantum sensing and photocatalysis applications.

Amines↗
DOE OSTI2025

Probing the Higgs boson CP properties in vector-boson fusion production in the H → τ + τ − channel with the ATLAS detector

The CP properties of the Higgs boson are studied in the vector-boson fusion production mode. The analysis exploits the decay mode of the Higgs boson into two τ-leptons using 140 fb −1 of proton-proton collision data at $\sqrt{s}$ = 13 TeV collected by the ATLAS experiment at the Large Hadron Collider. Results are obtained using the Optimal Observable method. CP-violating interactions between the Higgs boson and electroweak gauge bosons are considered in the effective field theory framework, with the interaction strength described in the HISZ basis by $\tilde{d}$, and in the Warsaw basis by c H$\tilde{W}$ , c H$\tilde{B}$ , and c H$\tilde{W}$B . No deviations relative to the Standard Model are observed, and limits are obtained on the strength parameters. The $\tilde{d}$ parameter is constrained to the interval [−0.012, 0.044] at the 95% confidence level while c H$\tilde{W}$ is constrained to [−0.24, 0.83], when considering both linear and quadratic effects of physics beyond the Standard Model.

Hadron-Hadron Scattering↗
DOE OSTI2025

Using Separation-Enhanced Isotope Ratio Mass Spectrometry to Enable Increased Renewable Carbon Content in Transportation Fuels (CRADA 525)

Stable isotope ratio measurements of carbon atoms using isotope ratio mass spectrometry (IRMS) can be an effective tool for quantifying biogenic carbon in co-processed fuels, with results approaching the precision and accuracy of accelerator mass spectrometry (AMS). The lower cost of an IRMS may enable deployment to refineries, improving access and analysis turnaround times (≤2 hours), and, by extension, provide data that can allow process optimization to maximize renewable carbon in desired refinery products. This project explored the integration of chemical separation with IRMS analyses to enable highly detailed tracking of biogenic carbon into fuel product streams separated by boiling point range, chemical class, or specific compound. Forty-nine fuels and fuel components of fossil and biogenic origin, spanning gasoline and diesel boiling point ranges, were received from three refiners and were analyzed for their δ 13 C values via IRMS. Results spanned a 13 C range from ca. 10‰ to 44‰ and reflect materials derived from sustainable sources (e.g., C4 or C3 plants, animal-based pathways, syngas) or from fossil-derived fuels. Common ranges are approximately 18‰ to 9‰ and approximately 30‰ to 20‰ for C4 and C3 plants, respectively, and approximately 34‰ to 24‰ and approximately 70‰ to 33‰ for petroleum-derived fuels and methane, respectively. Fuel-like standards were developed and tested using direct-injection elemental analyzer (EA) IRMS for liquid fuels. This method was compared with the published methods, yielding statistically similar results. Four blend curve sets were produced ranging from 0% to 100% of a fuel containing biogenic carbon, focusing on 0% to 10% biogenic carbon. Linear fits were the most applicable for two of the four blend curve sets; however, two sets were found to exhibit slightly quadratic behavior, which was more pronounced in low biogenic blend samples, necessitating second-order fits. The origin of the slight quadratic behavior remains unclear; however, the discussion points to possible interpretations. CanmetENERGY thoroughly characterized a majority of the samples using one- and two-dimensional gas chromatography (GC and GC×GC, respectively) and other analyses. Selected samples were subjected to solid phase extraction (SPE) for saturate, olefin, aromatic, and polar (SOAP) analysis, and the resulting solvent-diluted fractions containing saturates and aromatics were returned to Pacific Northwest National Laboratory (PNNL), where the solvent was removed via evaporation or physical separation using GC techniques. Characterization and separations provided an understanding of saturate and aromatic content, as well as boiling point ranges for each sample and sample fraction. Samples resulting from SPE were examined using EA-IRMS and gas chromatography combustion IRMS (GC-C-IRMS) analyses. Both approaches suggest that the range in values between end-members can be increased by selecting the paraffinic or aromatic fraction of the end-member or by selecting among individual compounds resulting from GC separation of the paraffinic fractions. Considerable work remains to put these approaches into practice and statistically validate the benefit for using a fraction or individual compound over bulk analysis of a sample. However, initial results suggest that separations provide advantages for samples having blend ratios of less than 10% biogenic blendstocks. 13 C results showed statistically similar biofuel blend results to those obtained at PNNL, although additional work is needed to obtain better reproducibility. Select samples were sent to Los Alamos National Laboratory (LANL) for IRMS measurements and Beta Analytics for AMS measurements. This work suggests that IRMS and AMS yield closely comparable results and in some circumstances, IRMS could serve as a surrogate for AMS. While additional work is needed to better resolve statistical advantages for separations and better show the comparable nature of IRMS and AMS in both the biogenic carbon analysis of bulk chemical classes, initial results from this study suggest that these should be pursued in order to proliferate this approach for quantifying biogenic carbon in transportation fuels to the refinery level, thereby potentially enabling process optimization in co-processing scenarios.

09 BIOMASS FUELS↗
DOE OSTI2025

AI-Optimized Polarization at Jefferson Lab

The AI-Optimized Polarization project seeks to develop experimental control applications for polarized targets and beams at Jefferson Lab using AI/ML. This paper will focus on two ongoing efforts involving a cryogenic polarized target and a linearly-polarized photon beam. Firstly, cryogenic targets, such as those used in Halls B and C (and approved for Hall D), are complex systems that are sensitive to a number of factors, including the temperature, beam currents, and the microwave and NMR apparatus. Secondly, the Hall D photon beam polarization depends on the optimal orientation of a diamond radiator, which produces coherent bremsstrahlung radiation from the electron beam incident upon it. Manual operation of both systems is tedious and error prone; implementing well-designed, interpretable control systems that incorporate AI is expected to lead to improved real-time polarization. AI optimization of nuclear physics experiments will lead, not just to cost-savings, but also to more efficient and higher-quality data, and this project will help to lay the foundation for future autonomous experiments.

Moran, Patrick [College of William and Mary, Willi↗
DOE OSTI2025

Geometric GNNs for charged particle tracking at GlueX

Nuclear physics experiments are aimed at uncovering the fundamental building blocks of matter. The experiments involve high-energy collisions that produce complex events with many particle trajectories. Tracking charged particles resulting from collisions in the presence of a strong magnetic field is critical to enable the reconstruction of particle trajectories and precise determination of interactions. It is traditionally achieved through combinatorial approaches that scale worse than linearly as the number of hits grows. Since particle hit data naturally form a point cloud and can be structured as graphs, graph neural networks (GNNs) emerge as an intuitive and effective choice for this task. In this study, we evaluate the GNN model for track finding on the data from the GlueX experiment at Jefferson Lab. We use simulation data to train the model and test on both simulation and real GlueX measurements. We demonstrate that GNN-based track finding outperforms the currently used traditional method at GlueX in terms of segment-based efficiency at a fixed purity while providing faster inferences. We show that the GNN model can achieve significant speedup by processing multiple events in batches, which exploits the parallel computation capability of graphical processing units (GPUs). Finally, we compare the GNN implementation on GPU and field-programmable gate array and describe the trade-off.

batched GNN pipeline↗
DOE OSTI2025

Phonon lifetimes and mode softening in cubic Cs 2 AgBiBr 6

Metal halide perovskites (MHPs) have emerged as noteworthy candidates for photovoltaic applications in recent years. Their high power-conversion efficiency is largely attributed to the structure-function relationship, which is not well understood. In this study, we examine the low-frequency phonons of the double perovskite Cs 2 AgBiBr 6 using neutron inelastic scattering. We find that the acoustic phonon lifetimes decrease from 16 to 3 ps along Γ to X, which is indicative of significant anharmonicity that contributes to the ultralow thermal conductivity. Additionally, we observe a linear temperature dependence of the square of the zone-center soft optical phonon energy, which is consistent with a weakly first-order displacive cubic-tetragonal structural phase transition. These results provide a deeper understanding of the lattice dynamics and phase transitions in MHPs as well as the effects of anharmonicity by comparison with prototypical hybrid organic-inorganic MHPs, (CH 3 NH 3 ) PbX 3 (X = Br, I).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Proton Quenching in Rare-Earth Inorganic Scintillators: GAGG:Ce and YSO:Ce

Scintillator detectors are an integral component of radiation detection systems for a variety of applications such as medical imaging, accelerator diagnostics, and space science. Typically, a scintillator detector’s response is characterized using gamma sources to understand the detection response to different types of radiation, including charged particle detection. However, there exists a nonlinearity of the amount of light produced from an incident gamma ray of specific energy and the light produced from an incident charged particle of the same energy. This important effect, known as quenching, must be accounted for to interpret energies from charged particles incident on detectors. In this article, we present results of quenching parameterization for two types of cerium-doped inorganic scintillators, Y2SiO5:Ce (YSO:Ce) and Gd3Al2Ga3O12:Ce (GAGG:Ce). We measured the light output from incident proton energies from 1 to 25 MeV using a 3-MV tandem accelerator and two reactions: Au(p,p)Au and 3He(d,p)⁴He. Using gamma-ray sources to calibrate the detectors, we compared the measured electron-equivalent energy versus the incident energy expected. Using an adaptation of the Birks semi-empirical formula, we extracted the Birks parameter (kB) to understand quenching. For one of the GAGG:Ce samples, the kB parameter of 0.0072 [g cm-2 MeV-1] is comparable to a similar study where the value of kB was 0.0065 [g cm-2 MeV-1]. For YSO:Ce, no other kB values were found in the literature. Three different types of GAGG:Ce were used to collect measurements of kB as a function of dopant concentration.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Twisted Tin‐Chloride Perovskite Single‐Crystal Heterostructures

Self-assembly affords simpler synthetic routes to heterostructures compared with manual layer-by-layer stacking, yet controlling interlayer twist angles in a bulk solid remains an outstanding challenge. We report two new single-crystal heterostructures: (Sn 2 Cl 2 )(CYS) 2 SnCl 4 (CYS = + NH 3 (CH 2 ) 2 S – ; Sn_CYS) and (Sn 2 Cl 2 )(SeCYS) 2 SnCl 4 (SeCYS = + NH 3 (CH 2 ) 2 Se – ; Sn_SeCYS) synthesized in solution, with alternating perovskite and intergrowth layers. Notably, compared to the recently reported lead analog, (Pb 2 Cl 2 )(CYS) 2 PbCl 4 (Pb_CYS), the tin heterostructures feature a twist between the perovskite and intergrowth layers. We trace this twist to local distortions at the Sn centers, which change the interfacial lattice-matching requirements compared to those of the Pb analog. Electronic band structure calculations show that the striking differences in the relative energies of perovskite- and intergrowth-derived bands in Sn_CYS and Pb_CYS arise from structural and not compositional differences. The structural anisotropy of Sn_CYS is also reflected in a large in-plane photoluminescence linear anisotropy ratio. Interfacial strain further affords differential incorporation of Pb into the perovskite and intergrowth layers of the Sn heterostructures, resulting in redshifted optical absorption onsets. Thus, we posit that local structural distortions may be exploited to manipulate the twist angle and interfacial strain in bulk heterostructures, providing a new handle for tuning the band alignments of bulk quantum-well electronic structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Synthesis of Correct Digital Controller Models from Specifications by Model Transformation (21-0320)

The design of high consequence controllers (in weapons systems, autonomy, etc.) that do what they are supposed to do is a significant challenge. Testing simply does not come close to meeting the requirements for assurance. Today circuit designers at Sandia (and elsewhere) typically capture the core behavior of their components using state models in tools such as STATEFLOW. They then check that their models meet certain requirements (e.g. “The system bus must not deadlock” or “both traffic lights at an intersection must not be green at the same time”) using tools called model checkers. If the model checker returns “yes” then the property is guaranteed to be satisfied by the model. However, there are several drawbacks to this industry practice: (1) there is a lot of detail to get right, this is particularly challenging when there are multiple components requiring complex coordination (2) any errors returned by the model checker have to be traced back through the design and fixed, necessitating rework, (3) there are severe scalability problems with this approach, particularly when dealing with concurrency. All this places high demands on the designers who now face not only an accelerated schedule but also controllers of increasing complexity. This report describes a new and fundamentally different approach to the construction of safety-critical digital controllers. Instead of directly constructing a complete model and then trying to verify it, the designer can start with an initial abstract (think “sketch”) model plus the requirements, from which a correct concrete model is automatically synthesized. There is no need for post-hoc verification of required functional properties. Having tool to carry this out will significantly impact the nation’s ability to ensure the safety of high-consequence digital systems. The approach has been implemented in a prototype tool, along with a suite of examples, including ones that reflect actual problems faced by designers. Our approach operates on a variant of Statecharts developed at Sandia called Qspecs. Statecharts are a widely used formalism for developing concurrent reactive systems, supporting scalability through allowing state models containing composite states, which are the serial or parallel composition of substates which can themselves contain statecharts. Statecharts enable an incremental style of development, in which states are progressively refined to incorporate greater detail in an incremental model of software development. Our approach formulates a set of constraints from the structure of the models and the requirements and propagates these constraints to a fixpoint. The solution to the constraints is an inductive invariant along with guards on the transitions. We also show how our approach extends to implementation refinement, decomposition, composition, and elaboration. We currently handle safety requirements written in LTL (Linear Temporal Logic)

42 ENGINEERING↗
DOE OSTI2025

CGSim: A Simulation Framework for Large Scale Distributed Computing Environment

Large-scale distributed computing infrastructures such as the Worldwide LHC Computing Grid (WLCG) require comprehensive simulation tools for evaluating performance, testing new algorithms, and optimizing resource allocation strategies. However, existing simulators suffer from limited scalability, hardwired algorithms, lack of real-time monitoring, and inability to generate datasets suitable for modern machine learning approaches. We present CGSim, a simulation framework for large-scale distributed computing environments that addresses these limitations. Built upon the validated SimGrid simulation framework, CGSim provides high-level abstractions for modeling heterogeneous grid environments while maintaining accuracy and scalability. Key features include a modular plugin mechanism for testing custom workflow scheduling and data movement policies, interactive real-time visualization dashboards, and automatic generation of event-level datasets suitable for AI-assisted performance modeling. We demonstrate CGSim’s capabilities through a comprehensive evaluation using production ATLAS PanDA workloads, showing significant calibration accuracy improvements across WLCG computing sites. Scalability experiments show near-linear scaling for multi-site simulations, with distributed workloads achieving 6 × better performance compared to single-site execution. The framework enables researchers to simulate WLCG-scale infrastructures with hundreds of sites and thousands of concurrent jobs within practical time budget constraints on commodity hardware.

Vatsavai, Sairam Sri [Brookhaven National Laborato↗
DOE OSTI2024

Calcium-organic matter fouling in nanofiltration: Synchrotron-based X-ray fluorescence and absorption near-edge structure spectroscopy for speciation

Calcium (Ca)-enhanced organic matter (OM) fouling of nanofiltration (NF) membranes leads to reduced flux during desalination and requires frequent cleaning. Fouling mechanisms are not fully understood, which limits the development of targeted fouling control methods. This study employed synchrotron-based X-ray fluorescence (XRF) and X-ray absorption near-edge structure (XANES) spectroscopy to quantify the spatial distribution and mass of Ca deposition as well as changes in the Ca coordination environment characteristic of specific fouling mechanisms, respectively. Bench-scale filtration experiments were performed using feed solutions containing Ca and ten different types of organic matter (OM), as well as the common scalants, calcium carbonate (CaCO 3 ) and calcium sulfate (CaSO 4 ). Osmotic backwash (OB) was performed at regular intervals for fouling control. Ca-OM aggregation resulted in greater flux decline and lower flux recovery during OB than Ca conditioning of membranes followed by filtration of feed solution with OM. Linear combination fitting (LCF) of XANES absorption spectra from fouled membranes indicated that Ca-OM aggregation preferentially occurred for OM types that exhibited both high carboxylic group and negative charge density. Consequently, these OM types exhibited greater deposition of Ca and TOC on the membrane surface when compared to other OM types. For the coexistence of scalants and OM, Ca speciation within the fouling layer was characteristic of both Ca bound to the membrane (i.e. potential bridging, charge screening) as well as Ca-OM aggregation and deposition mechanisms, while a range of crystal polymorphs were observed to occur simultaneously. XRF and XANES represent powerful tools for the elucidation of NF fouling mechanisms by quantification of Ca deposition as well as Ca speciation. Fouling control methods should target OM types with high carboxyl group density and negative charge to neutralize or eliminate interactions with Ca.

42 ENGINEERING↗
DOE OSTI2025

Maximizing quantum enhancement in axion dark matter experiments

We provide a comprehensive comparison of linear amplifiers and microwave photon counters in axion dark matter experiments. The study is done assuming a range of realistic operating conditions and detector parameters, over the frequency range between 1 and 30 GHz. As expected, photon counters are found to be advantageous under low background, at high frequencies (𝜈 >5 GHz), if they can be implemented with robust wide-frequency tuning or a very low dark count rate. Additional noteworthy observations emerging from this study include: (1) an expanded applicability of off-resonance photon background reduction, including the single-quadrature state squeezing, for scan rate enhancements; (2) a much broader appeal for operating the haloscope resonators in the overcoupling regime, up to 𝛽 ∼10; (3) the need for a detailed investigation into the cryogenic and electromagnetic conditions inside haloscope cavities to lower the photon temperature for future experiments; (4) the necessity to develop a distributed network of coupling ports in high-volume axion haloscopes to utilize these potential gains in the scan rate.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Compensating slice emittance growth in high brightness photoinjectors using sacrificial charge

Achieving maximum electron beam brightness in photoinjectors requires detailed control of the 3D bunch shape and precise tuning of the beam focusing. Even in state-of-the-art designs, slice emittance growth due to nonlinear space charge forces and partial nonlaminarity often remains non-negligible. In this work, we introduce a new means to linearize the transverse slice phase space: a sacrificial portion of the bunch’s own charge distribution, formed into a wavebroken shock front by highly nonlinear space charge forces within the gun, whose downstream purpose is to dynamically linearize the desired bunch core. We show that linearization of an appropriately prepared bunch can be achieved via strongly nonlaminar focusing of the sacrificial shock front, while the inner core focuses laminarly. This leads to a natural spatial separation of the two distributions: a dense core surrounded by a diffuse halo of sacrificial charge that can be collimated. Multiobjective genetic algorithm optimizations of the ultracompact x-ray free electron laser injector employ this concept, and we interpret it with an analytic model that agrees well with the simulations. In simulation, we demonstrate a final bunch charge of 100 pC, peak current ∼ 30 A, and a sacrificial charge of 150 pC (250 pC total emitted from cathode) with normalized emittance growth of < 20 nm rad due to space charge. This implies a maximum achievable brightness approximately an order of magnitude greater than existing free electron laser injector designs. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

MAGIC: M arching Cubes Isosurface Uncertainty Visualization for G auss i an Uncertain Data With Spatial C orrelation

Here, in this paper, we study the propagation of data uncertainty through the marching cubes algorithm for isosurface visualization for correlated uncertain data. Consideration of correlation has been shown paramount for avoiding errors in uncertainty quantification and visualization in multiple prior studies. Although the problem of isosurface uncertainty with spatial data correlation has been previously addressed, there are two major limitations to prior treatments. First, there are no analytical formulations for uncertainty quantification of isosurfaces when the data uncertainty is characterized by a Gaussian distribution with spatial correlation. Second, as a consequence of the lack of analytical formulations,existing techniques resort to a Monte Carlo sampling approach, which is expensive and difficult to integrate into visualization tools. To address these limitations, we present a closed-form framework to efficiently derive uncertainty in marching cubes level-sets for Gaussian uncertain data with spatial correlation (MAGIC). To derive closed-form solutions, we leverage the Hinkley's derivation on the ratio of Gaussian distributions. With our analytical framework, we achieve a significant speed-up and enhanced accuracy of uncertainty quantification over classical Monte Carlo methods. We further accelerate our analytical solutions using many-core processors to achieve speed-ups up to 585× and integrability with production visualization tools for broader impact. We demonstrate the effectiveness of our correlation-aware uncertainty framework through experiments on meteorology, urban flow, and astrophysics simulation datasets.

Gaussian↗
DOE OSTI2024

Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity

A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f⁡(vy) or f⁡(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗