SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,081 records · Page 60

All sources NASA NTRS DOE OSTI

DOE OSTI2022

A Critical Review of the Circular Economy for Lithium-Ion Batteries and Photovoltaic Modules: Status, Challenges, and Opportunities

To meet net-zero emissions and cost targets for power production, recent analysis indicates that photovoltaic (PV) capacity in the United States could exceed 1 TW by 2050 alongside comparable levels of energy storage capacity, mostly from batteries. For comparison, the total U.S. utility-scale power capacity from all energy sources in 2020 was 1.2 TW (EIA 2022), of which solar satisfied approximately 3% (DOE 2021). With such massive scales of deployment, questions have arisen regarding issues of material supply for manufacturing, end-of-life management of technologies, environmental impacts across the life cycle, and economic costs to both individual consumers and society at large. A set of solutions to address these issues center on the development of a circular economy - shifting from a take-make-waste linear economic model to one that retains the value of materials and products as long as possible, recovering materials at end of life to recirculate back into the economy. With limited global experience, scholars and practitioners have begun to investigate circular economy pathways, focusing on applying novel technologies and analytical methods to fast-growing sectors like renewable energy. This critical review aims to synthesize the growing literature to identify key insights, gaps, and opportunities for research and implementation of a circular economy for two of the leading technologies that enable the transition to a renewable energy economy: solar PV and lithium-ion batteries (LIBs). We apply state-of-the-science systematic literature review procedures to critically analyze over 3,000 publications on the circular economy of solar PV and LIBs, categorizing those that pass a series of objective screens in ways that can illuminate the current state of the art, highlight existing impediments to a circular economy, and recommend future technological and analytical research. We conclude that while neither PV nor LIB industries have reached a circular economy, they are both on a path towards increased circularity. Based on our assessment of the state of current literature and scientific understanding, we recommend research move beyond its prior emphasis on recycling technology development to more comprehensively investigate other CE strategies, more holistically consider economic, environmental and policy aspects of CE strategies, increase leveraging of digital information systems that can support acceleration towards a CE, and to continue to study CE-related aspects of LIB and PV markets.

circular economy↗

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 OSTI2021

Preparing an Incompressible-Flow Fluid Dynamics Code for Exascale-Class Wind Energy Simulations

The U.S. Department of Energy has identified exascale-class wind farm simulation as critical to wind energy scientific discovery. A primary objective of the ExaWind project is to build high-performance, predictive computational fluid dynamics (CFD) tools that satisfy these modeling needs. GPU accelerators will serve as the computational thoroughbreds of next-generation, exascale-class supercomputers. Here, we report on our efforts in preparing the ExaWind unstructured mesh solver, Nalu-Wind, for exascale-class machines. For computing at this scale, a simple port of the incompressible-flow algorithms to GPUs is insufficient. To achieve high performance, one needs novel algorithms that are application aware, memory efficient, and optimized for the latest-generation GPU devices. The result of our efforts are unstructured-mesh simulations of wind turbines that can effectively leverage thousands of GPUs. In particular, we demonstrate a first-of-its-kind, incompressible-flow simulation using Algebraic Multigrid solvers that strong scales to more than 4000 GPUs on the Summit supercomputer.

algebraic multigrid↗

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 OSTI2021

Shotgun Immunoproteomic Approach for the Discovery of Linear B-Cell Epitopes in Biothreat Agents Francisella tularensis and Burkholderia pseudomallei

Peptide-based subunit vaccines are coming to the forefront of current vaccine approaches, with safety and cost-effective production among their top advantages. Peptide vaccine formulations consist of multiple synthetic linear epitopes that together trigger desired immune responses that can result in robust immune memory. The advantages of linear compared to conformational epitopes are their simple structure, ease of synthesis, and ability to stimulate immune responses by means that do not require complex 3D conformation. Prediction of linear epitopes through use of computational tools is fast and cost-effective, but typically of low accuracy, necessitating extensive experimentation to verify results. On the other hand, identification of linear epitopes through experimental screening has been an inefficient process that requires thorough characterization of previously identified full-length protein antigens, or laborious techniques involving genetic manipulation of organisms. In this study, we apply a newly developed generalizable screening method that enables efficient identification of B-cell epitopes in the proteomes of pathogenic bacteria. As a test case, we used this method to identify epitopes in the proteome of Francisella tularensis (Ft), a Select Agent with a well-characterized immunoproteome. Our screen identified many peptides that map to known antigens, including verified and predicted outer membrane proteins and extracellular proteins, validating the utility of this approach. We then used the method to identify seroreactive peptides in the less characterized immunoproteome of Select Agent Burkholderia pseudomallei (Bp). This screen revealed known Bp antigens as well as proteins that have not been previously identified as antigens. Although B-cell epitope prediction tools Bepipred 2.0 and iBCE-EL classified many of our seroreactive peptides as epitopes, they did not score them significantly higher than the non-reactive tryptic peptides in our study, nor did they assign higher scores to seroreactive peptides from known Ft or Bp antigens, highlighting the need for experimental data instead of relying on computational epitope predictions alone. The present workflow is easily adaptable to detecting peptide targets relevant to the immune systems of other mammalian species, including humans (depending upon the availability of convalescent sera from patients), and could aid in accelerating the discovery of B-cell epitopes and development of vaccines to counter emerging biological threats.

60 APPLIED LIFE SCIENCES↗

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 OSTI2023

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

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 OSTI2021

Methane Mitigation Thermoelectric Generator (MMTEG) (Final Scientific/ Technical Report)

Gas Technology Institute (GTI) has completed a $1.815M (plus $500K cost share), 54-month Co-operative agreement with the Department of Energy’s National Energy Technology Laboratory to develop a Methane Mitigation Thermoelectric Generator (MMTEG) system for gas field applications. This novel system uses fugitive gas to produce electrical power and consists of Thermoelectric Generators (TEG) driven by a linear burner, an air compressor, an accumulator, valves, batteries, and power electronics. The electrical power generated by the system is used to compress air and, in turn, the air is used to operate the pneumatic valves at the well site instead of using natural gas (NG) which would then be vented. The team completed the project objectives which included (1) Design, fabricate and test an integrated 6We nominal MMTEG prototype system, (2) Design a “retrofit kit” MMTEG system capable of being field tested at a gas well production site, and (3) build and test the field MMTEG system in a laboratory environment. Multiple system options were developed prior to selecting a “Passive” system which meets cost, NG savings, and greenhouse gas (GHG) reduction targets although not as efficient as initially planned. The MMTEG system is built primarily from commercial off-the-shelf parts (in some cases re-purposed) including the heat exchanger, heat rejection, electronic components, and also the TEGs. The team developed and implemented a novel system including the control system developed by Morrison Applied Sciences (MAS). The team completed incremental demonstrations of the hardware prior to the MMTEG system demonstration. A one-year simulation of the air delivery and battery charging subsystems was completed prior to the integration into the MMTEG system. In addition, the team simulated two-years of thermal cycles for integrated Burner/TEG/heat rejection subsystem. Finally, the entire MMTEG system was assembled, and troubleshooting was completed over a two-week period. Next, the MMTEG system was tested to simulate over 15 weeks of entire system operation over approximately four weeks in an accelerated test fashion with minimal intervention (such as changing fuel tanks). The MMTEG system met the key goals of a unit cost of under $1500 while saving 97% of NG expended today (including leakage) on average by pneumatic systems venting to the atmosphere. The MMTEG system reduced GHG emissions by 1000X (using the methane intensification factor of 28 relative to CO2). GTI is currently pursuing a field test of the MMTEG system. Interfacing with producers has provided additional insight into system improvements. Planned improvements include additional weather protection and control system improvements including the implementation of a long- range radio capability to notify operators if there is a fault.

03 NATURAL GAS↗

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↗

DOE OSTI2025

Atacama Cosmology Telescope DR6 and DESI: Structure growth measurements from the cross-correlation of DESI legacy imaging galaxies and CMB lensing from ACT DR6 and 𝑃⁢𝑙⁢𝑎⁢𝑛⁢𝑐⁢𝑘 PR4

We measure the growth of cosmic density fluctuations on large scales and across the redshift range 0.3 < 𝑧 < 0.8 through galaxy clustering and the cross-correlation of the ACT data release 6 cosmic microwave background (CMB) lensing map and galaxies from the Dark Energy Spectroscopic Instrument Legacy Survey, using three galaxy samples spanning the redshifts of 0.3 ≲ 𝑧 ≲ 0.45, 0.45 ≲ 𝑧 ≲ 0.6, 0.6 ≲ 𝑧 ≲ 0.8. We adopt a scale cut where nonlinear effects are negligible, so that the cosmological constraints are derived from the linear regime. We determine the amplitude of matter fluctuations over all three redshift bins using Atacama Cosmology Telescope (ACT) data alone to be 𝑆 8 ≡ 𝜎 8 ⁢(Ω 𝑚 /0.3) 0.5 =0.772 ± 0.040 in a joint analysis combining the three redshift bins and ACT lensing alone. Using a combination of ACT and Planck data we obtain 𝑆 8 = 0.765 ± 0.032. The lowest redshift bin used is the least constraining and exhibits a ∼2⁢𝜎 tension with the other redshift bins; thus we also report constraints excluding the first redshift bin, giving 𝑆 8 = 0.785 ± 0.033 for the combination of ACT and Planck. This result is in excellent agreement at the 0.3⁢𝜎 level with measurements from galaxy lensing, but is 1.8⁢𝜎 lower than predictions based on Planck primary CMB data. Understanding whether this hint of discrepancy in the growth of structure at low redshifts arises from a fluctuation, from systematics in data, or from new physics is a high priority for forthcoming CMB lensing and galaxy cross-correlation analyses.

cosmic microwave background↗

DOE OSTI2024

Reducing beta functions in the NSLS-II long straight section for NEXT-III beamlines

The NSLS-II Experimental Tools III (NEXT-III) Project with the goal of building 8-12 new beamlines is ongoing. We were requested to explore the possibility of modifying the long straight section in Cell 26 to accommodate two new undulators for Advanced Nanoscale Imaging (ANI) Beamline and Tender X-ray Nanoprobe (TXN) Beamline. The scope of work includes: linear lattice matching; confirmation of the minimum ID gaps consistent with beam stay-clear limitations; calculation of the X-ray brightness; optimization of the dynamic aperture and momentum acceptance; multi-particle tracking to confirm the injection efficiency does not deteriorate significantly. The beam simulations have to be done using the most recent NSLS-II lattice model with all present insertion devices and realistic magnet errors, misalignments, and apertures. Originally two lattice options were considered for the assessment:1) with single -minimum beta functions and 2) NSLS-II with double-minimum beta functions. However, since the beamlines are not interested in the reduction of the horizontal beta function, just the vertical one, we decided to stop considering the single-minimum option because the current vertical beta is already low (3.3m) and further reducing it at the straight section center can increase it too much at the edge of a long ID, and limit the beam stay-clear. So, this report is focused on the double mini-beta scheme only.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

Entanglement, trace anomaly, and confinement in QCD

We formulate confinement in quantum chromodynamics (QCD) as an entropic surface phenomenon. Quark and gluon quantum information is localized on a transverse, entangling two-sphere of radius 𝑅 𝐸⁢𝐸 ; at this radius the QCD vacuum—partitioned by a hadron into interior and exterior regions—reaches its maximal entanglement entropy. Lattice-QCD determinations of the scalar (trace) gravitational form factors fix both 𝑅 𝐸⁢𝐸 and the transverse trace-anomaly density 𝜌 ℎ ⁡(𝑅 𝐸⁢𝐸 ), yielding a parameter-free slope 𝑐 ℎ = 8⁢𝜋 2 ⁢𝑅𝑆$^2_{𝐸⁢𝐸}$𝜌 ℎ ⁡(𝑅 𝐸⁢𝐸 ) and a mechanical entropy 𝑆 𝐸⁢𝐸 ⁡(𝑦) = 𝑐 ℎ ⁢𝑦 that grows linearly with rapidity 𝑦. The entropy gradient ∂ 𝑅 𝑆 𝐸⁢𝐸 changes sign at 𝑅 𝐸⁢𝐸 : it pushes colored degrees of freedom outward for 𝑟 <𝑅 𝐸⁢𝐸 and pulls them inward for 𝑟 >𝑅 𝐸⁢𝐸 , thereby localizing them on the codimension-2 entangling two-sphere Σ ⊥ = 𝑆$^2_{𝑅_{𝐸⁢𝐸}}$ (which, in the infinite-momentum frame (IMF), projects onto the transverse plane)—the “information wall.” This provides a high-energy (large-𝑦) entropic confinement diagnostic that complements—rather than replaces—Wilson’s area-law criterion, which probes long-distance dynamics near the rest frame (𝑦 → 0). Imposing unitarity on an entropic ansatz for the amplitude yields 𝜎⁡(𝑠) ∝ 𝑦 𝛿 . World data favor 𝛿 = 2 for elastic 𝑝⁡𝑝⁡($𝑝\bar{⁡𝑝}$) scattering and heavy-quark photoproduction, whereas 𝜙 photoproduction favors a softer 𝛿 = 0.387. All extracted cross sections remain well below the Froissart-Martin bound. These results provide a confinement criterion quantified directly from nonperturbative QCD inputs, unifying the trace anomaly, entanglement entropy, and high-energy scattering within a single quantitative framework.

Color confinement↗