Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complex Constraints”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Impact of host climate model on contrail cirrus effective radiative forcing estimates

Estimates of aviation effective radiative forcing (ERF) indicate that contrail cirrus is currently its largest contributor, although with a substantial associated uncertainty of ∼ 70 %. Here, we implement the contrail parameterisation developed for the Community Atmosphere Model (CAM) in the UK Met Office Unified Model (UM), allowing us to compare, for the first time, the impact of key features of the host climate model on contrail cirrus ERF. We find that differences in background humidity between the models result in the UM-simulated contrail fractions being 2 to 3 times larger than in CAM. Additionally, the models show contrasting responses in overall global cloud fraction, with contrails increasing the total cloud fraction in the UM and decreasing it in CAM. Differences in the complexity of the cloud microphysics schemes lead to significant differences in simulated changes to cloud ice water content due to aviation. After compensating for the unrealistically low contrail optical depth in the UM, we estimate the 2018 contrail cirrus ERF to be 40.8 mW m−2 in the UM, compared to 60.1 mW m−2 in CAM. These values highlight the substantial uncertainty in contrail cirrus ERF due to differences in microphysics and radiation schemes between the two models. We also find a factor-of-8 uncertainty in contrail cirrus ERF due to existing uncertainty in contrail cirrus optical depth. Future research should focus on better representing microphysical and radiative contrail characteristics in climate models and on improved observational constraints.

Zhang, Weiyu [University of Leeds (United Kingdom)↗
DOE OSTI2024

Resonant multiscalar production in the generic complex singlet model in the multi-TeV region

We develop benchmarks for resonant discalar production in the generic complex singlet scalar extension of the Standard Model (SM) with no additional symmetries, which contains two new scalars. These benchmarks maximize discalar resonant production modes at future 𝑝⁢𝑝 colliders: 𝑝⁢𝑝 → ℎ 2 → ℎ 1⁡ ℎ 1 , 𝑝⁢𝑝 → ℎ 2 → ℎ 1 ⁡ℎ 3 , and 𝑝⁢𝑝 → ℎ 2 → ℎ 3 ⁡ℎ 3 , where ℎ 1 is the observed SM-like Higgs boson and ℎ 2,3 are new scalars. The decays ℎ 2 → ℎ 1 ⁡ℎ 3 and ℎ 2 → ℎ 3 ⁡ℎ 3 may be the only way to discover ℎ 3 , leading to a discovery of two new scalars at once. Current LHC and projected future collider (HL-LHC, FCC-ee+HL-LHC, IL⁢C⁢500+HL-LHC) constraints on this model are used to produce benchmarks at the HL-LHC for ℎ 2 masses between 250 GeV and 1 TeV and a future 𝑝⁢𝑝 collider (FCC-hh) for ℎ 2 masses between 250 GeV and 12 TeV. We update the current LHC bounds on the singlet-Higgs boson mixing angle for these benchmarks. As the mass of ℎ 2 approaches the multi-TeV region, certain limiting behaviors of the maximum rates are uncovered due to theoretical constraints on the parameters. These limits, which can be derived analytically, are BR⁡(ℎ 2 → ℎ 1⁡ ℎ 1 ) → 0.25, BR⁡(ℎ 2 → ℎ 3 ⁡ℎ 3 ) → 0.5, and BR⁡(ℎ 2 → ℎ 1⁡ ℎ 3 ) → 0. It can also be shown that the maximum rates of 𝑝⁢𝑝 → ℎ 2 → ℎ 1 ⁡ℎ 1 and 𝑝⁢𝑝 → ℎ 2 → ℎ 3⁡ ℎ 3 approach the same value. Hence, all three ℎ 2 → ℎ 𝑖 ⁡ℎ 𝑗 decays are promising discovery modes for ℎ 2 masses at and below 𝒪⁡(1 TeV), while above 𝒪⁡(1 TeV) the decays ℎ 2 → ℎ 1 ⁡ℎ 1 and ℎ 2 → ℎ 3 ⁡ℎ 3 are more encouraging. We choose benchmark masses for ℎ 3 to produce a large range of decay signatures including multi-𝑏, multivector boson, and multi-SM-like Higgs production. As we will show, the behavior of the maximum rates leads to the surprising conclusion that in the multi-TeV region this model may be discovered in the Higgs quartet production mode via ℎ 2 → ℎ 3 ⁡ℎ 3 → 4⁢ℎ 1 decays before Higgs triple production is observed. The maximum di- and four Higgs production rates are similar in the multi-TeV range.

Extensions of Higgs sector↗
DOE OSTI2026

Mitigating Transition in the Fermilab Booster Using a Triple Phase Jump

The PIP-II project will significantly enhance neutrino production for DUNE, Fermilab’s flagship long-baseline neutrino oscillation experiment, by doubling the beam power delivered by the accelerator complex. To achieve this, the total charge injected from the new PIP-II linac into the Booster rapid cycling synchrotron (RCS) will increase from $4.5\times 10^{12}$ to $6.5\times 10^{12}$ protons per pulse. Simultaneously, the Booster’s ramp rate will rise from 15 to 20 Hz, while the injection energy will go from 400 MeV to 800 MeV. The Booster accelerates the beam to 8 GeV and crosses transition at $\gamma_t=5.45$. In current operations, no dedicated $\gamma_t$ jump system is employed; instead, longitudinal emittance growth is controlled via an active quadrupole-mode feedback system. Downstream, the Recycler Ring uses slip-stacking to accumulate beam and increase bunch intensity. However, this process imposes a constraint on the longitudinal emittance at Booster extraction, limiting it to 0.1 eV-s (95\%) to avoid excessive particle loss. Because collective effects scale with intensity, additional measures to mitigate transition crossing may be required to stay below this limit. One potential approach is the so-called triple phase-jump technique. While the method has known limitations, it offers some advantages: it can be implemented using the existing digital low-level RF (LLRF) system, requires no additional magnets or pulsed power supplies, and remains compatible with quadrupole-mode feedback.

Ostiguy, J.-F. [Fermilab] (ORCID:0000000290883681)↗
DOE OSTI2025

Mitigating Transition in the Fermilab Booster Using a Triple Phase Jump

The PIP-II project will significantly enhance neutrino production for DUNE, Fermilab s flagship long-baseline neutrino oscillation experiment, by doubling the beam power delivered by the accelerator complex. To achieve this, the total charge injected from the new PIP-II linac into the Booster rapid cycling synchrotron (RCS) will increase from $4.5\times 10^{12}$ to $6.5\times 10^{12}$ protons per pulse. Simultaneously, the Booster s ramp rate will rise from 15 to 20 Hz, while the injection energy will go from 400 MeV to 800 MeV. The Booster accelerates the beam to 8 GeV and crosses transition at $\gamma_t=5.45$. In current operations, no dedicated $\gamma_t$ jump system is employed; instead, longitudinal emittance growth is controlled via an active quadrupole-mode feedback system. Downstream, the Recycler Ring uses slip-stacking to accumulate beam and increase bunch intensity. However, this process imposes a constraint on the longitudinal emittance at Booster extraction, limiting it to 0.1 eV-s (95\%) to avoid excessive particle loss. Because collective effects scale with intensity, additional measures to mitigate transition crossing may be required to stay below this limit. One potential approach is the so-called triple phase-jump technique. While the method has known limitations, it offers some advantages: it can be implemented using the existing digital low-level RF (LLRF) system, requires no additional magnets or pulsed power supplies, and remains compatible with quadrupole-mode feedback.

Ostiguy, Jean-Francois [Fermilab] (ORCID:000000029↗
DOE OSTI2026

A Decomposition-Based Learn-To-Optimize Approach with Feasibility Layer Assistance for Sub-Hourly Unit Commitment

Sub-hourly unit commitment (UC) with 15-min intervals is gaining significant attention as a way to respond rapidly to the fluctuations in electricity supply and demand introduced by renewable resources. However, the increased temporal resolution and complex inter-temporal dependencies pose substantial computational challenges for traditional optimization methods. To this end, this paper explores a decomposition-based learn-to-optimize approach. Building on recent advances in machine learning, our method revisits the long- overlooked Lagrangian relaxation framework, which is a classical decomposition technique that enables tractable subproblem solving. These smaller subproblems are inherently well-suited for machine learning, as their reduced dimensionality and structural regularity allow predictive models to efficiently learn and generalize solution patterns. We thus propose a generic predictive model, which embeds Gated Recurrent Units (GRUs) and Attention in the encoder-decoder structure, and integrate a rule-based feasibility layer to capture temporal dependencies, reduce training effort, and improve feasibility w.r.t. unit-level constraints. Our method has been validated on the IEEE 118-bus system, demonstrating promising performance in solving sub-hourly UC problems efficiently and feasibly.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Mauka Energy FEVER Tool Dataset

Mauka Energy’s dataset, developed under the Forestry Electric Vehicle Energy Routing (FEVER) project and funded by the U.S. Department of Energy’s Small Business Innovation Research program, is a high-resolution geospatial resource designed to support energy modeling for electric log trucks in complex forestry environments. The dataset integrates detailed spatial and road network data to enable accurate simulation of vehicle performance across varied terrain. At its core, the dataset incorporates lidar-derived elevation models, road alignments, and surface classifications from Oregon State University’s McDonald-Dunn Research Forest. These data capture fine-scale variations in slope, curvature, and surface conditions across forest road systems, allowing for vehicle-level analysis of energy consumption and recovery. The dataset also includes data collected on the surrounding public and private road networks in Benton County, Oregon, used in real-world haul routes. These connecting segments provide critical context for modeling transitions between forest operations and regional transportation infrastructure, incorporating attributes such as grade profiles, elevation change, and speed constraints. This combined dataset underpins the development of Mauka Energy’s rolldown tool, which quantifies energy use and regenerative braking potential on downhill and variable-grade segments. By leveraging high-resolution terrain and road data, the FEVER project enables more accurate assessment of electric vehicle feasibility and performance in forestry applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

PhaseT3M: 3D imaging at 1.6 Å resolution via electron cryo-tomography with nonlinear phase retrieval

Electron cryo-tomography (cryo-ET) enables 3D imaging of complex, radiation-sensitive structures with molecular detail. However, image contrast from the interference of scattered electrons is nonlinear with atomic density and multiple scattering further complicates interpretation. These effects degrade resolution, particularly in conventional reconstruction algorithms, which assume linearity. Particle averaging can reduce such issues but is unsuitable for heterogeneous or dynamic samples ubiquitous in biology, chemistry, and materials sciences. Here, we develop a phase retrieval-based cryo-ET method, PhaseT3M. We experimentally demonstrate its application to an approximately 7 nm Co3O4 nanoparticle on an approximately 30 nm carbon substrate, achieving a maximum resolution of 1.6 Å, surpassing conventional limits using standard cryo-TEM equipment. PhaseT3M uses a multislice model for multiple scattering and Bayesian optimization for alignment and computational aberration correction, with a positivity constraint to recover ‘missing wedge’ information. Applied directly to biological particles, it enhances reconstruction quality and reduces artifacts, establishing a standard for routine 3D imaging with phase contrast.

Biophysics↗
DOE OSTI2024

Embedded symmetric positive semi-definite machine-learned elements for reduced-order modeling in finite-element simulations with application to threaded fasteners

Here, we present a machine-learning strategy for finite element analysis of solid mechanics wherein we replace complex portions of a computational domain with a data-driven surrogate. In the proposed strategy, we decompose a computational domain into an “outer” coarse-scale domain that we resolve using a finite element method (FEM) and an “inner” fine-scale domain. We then develop a machine-learned (ML) model for the impact of the inner domain on the outer domain. In essence, for solid mechanics, our machine-learned surrogate performs static condensation of the inner domain degrees of freedom. This is achieved by learning the map from displacements on the inner-outer domain interface boundary to forces contributed by the inner domain to the outer domain on the same interface boundary. We consider two such mappings, one that directly maps from displacements to forces without constraints, and one that maps from displacements to forces by virtue of learning a symmetric positive semi-definite (SPSD) stiffness matrix. We demonstrate, in a simplified setting, that learning an SPSD stiffness matrix results in a coarse-scale problem that is well-posed with a unique solution. We present numerical experiments on several exemplars, ranging from finite deformations of a cube to finite deformations with contact of a fastener-bushing geometry. We demonstrate that enforcing an SPSD stiffness matrix drastically improves the robustness and accuracy of FEM–ML coupled simulations, and that the resulting methods can accurately characterize out-of-sample loading configurations with significant speedups over the standard FEM simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The Galaxy–Galaxy Strong Lensing Cross Section and the Internal Distribution of Matter in ΛCDM Substructure

Abstract Strong gravitational lensing offers a powerful probe of the detailed distribution of matter in lenses, while magnifying and bringing faint background sources into view. Observed strong lensing by massive galaxy clusters, which are often in complex dynamical states, has also been used to map their dark matter (DM) substructures on smaller scales. Deep high-resolution imaging has revealed the presence of strong lensing events associated with these substructures, namely galaxy-scale sub-halos. However, an inventory of these observed galaxy–galaxy strong lensing (GGSL) events is noted to be discrepant with state-of-the-art ΛCDM simulations. Cluster sub-halos appear to be over-concentrated compared to their simulated counterparts yielding an order-of-magnitude higher value of GGSL. In this paper, we explore the possibility of resolving this observed discrepancy by redistributing the mass within observed cluster sub-halos in ways that are consistent within the ΛCDM paradigm of structure formation. Lensing mass reconstructions from data provide constraints on the mass enclosed within apertures and are agnostic to the detailed mass profile within them. Therefore, as the detailed density profile within cluster sub-halos currently remains unconstrained by data, we are afforded the freedom to redistribute the enclosed mass. We investigate if rearranging the mass to a more centrally concentrated density profile helps alleviate the GGSL discrepancy. We report that refitting cluster sub-halos to the ubiquitous ΛCDM-motivated Navarro–Frenk–White profile, and further modifying them to include significant baryonic components, does not resolve this tension. A resolution to this persisting GGSL discrepancy may require more careful exploration of alternative DM models.

Astronomy & Astrophysics↗
DOE OSTI2025

A neural master equation framework for multiscale modeling of molecular processes: application to atomic-scale plasma processes

Plasma-surface interactions (PSI) play a crucial role in microelectronics fabrication; however, their multiscale nature and array of complex, often unknown interactions make computational modeling of PSIs extremely difficult. To this end, we propose a general neural master equation (NME) framework that uses master equations to describe the dynamics of a molecular process, wherein neural networks learned from atomistic simulations represent unknown transitions between different system states. By leveraging the physics-based structure of master equations and data-driven state transitions, the NME framework promotes generalizability and physics interpretability, and can bridge disparate length and time scales. The framework is demonstrated for multiscale modeling of Si atomic layer etching and reactive ion etching, where the learned NME-based surface kinetic models exhibit good predictive and extrapolative capabilities for predicting experimentally relevant observables as a function of process parameters. The NME-based surface kinetic models obey physical constraints, which are violated in models based on neural ordinary differential equations. The proposed NME framework for multiscale modeling of molecular processes can pave the way for the discovery of new chemistries and materials in atomic-scale plasma processes.

Chemical engineering↗
DOE OSTI2025

Even Higher-Level Synthesis: An Exploration of AI Hardware Accelerators using HLS4ML

With the rise of artificial intelligence, the popularization of deep learning, and a constantly evolving industry, the demand for flexible and efficient tools has never been greater. As algorithms grow more complex, their runtime and energy consumption increase exponentially. Customized hardware accelerators, long used for specific mathematical operations, remain essential for managing modern applications' computational and power demands. Hardware accelerators can speed up complex computations by orders of magnitude, but their manual design and verification processes are often challenging and time-consuming. High-Level Synthesis (HLS) provides a solution by transforming high-level algorithm descriptions, typically written in C++ or SystemC, into synthesizable RTL suitable for hardware implementation. This approach reduces development time for RTL engineers while offering flexibility beyond what traditional handwritten RTL can provide. We extended this capability to the machine-learning domain with the open-source framework hls4ml, which allows neural networks trained in Python frameworks like Tensorflow or PyTorch to be synthesized into efficient hardware representations for the traditional FPGA and ASIC flows. This breakthrough addresses the growing need for reduced design turnaround and easy verification of ML hardware accelerators with low latency and power efficiency constraints. During this tutorial, we will demonstrate how Python complements HLS by simplifying the ML design process, bridging the gap between software and hardware development. Attendees will explore how we translate neural networks modeled in Python into fixed-point C++ models suitable for HLS workflows. We will dive into strategies like Value-Range Analysis and Quantization-Aware Training, which optimize these designs for deployment and evaluate their accuracy, power consumption, and energy efficiency. To exemplify these concepts, experts from Fermilab will share their experiences applying this technology to high-energy physics experiments, where real-time, low-latency processing is critical. Over the years, Fermilab engineers have demonstrated how deep neural networks, optimized for hardware using hls4ml, can meet the stringent requirements of trigger systems at the CERN Large Hadron Collider. These systems rely on rapid decision-making to process immense data volumes while retaining only the most relevant events for further analysis. The application of hls4ml has also been extended to innovative technologies like smart pixel arrays. These smart pixels integrate ML inference capabilities directly into sensor devices, enabling localized data processing at the pixel level. This approach drastically reduces the need to transmit raw data to external processing units, significantly decreasing power consumption and latency. By embedding neural networks within the pixel architecture, the smart pixels can identify and prioritize relevant data in real time, providing a highly efficient solution for edge computing in scenarios such as particle detectors and imaging systems. Fermilab's work highlights the potential of hardware-accelerated ML in scenarios where both speed and power efficiency are mission-critical. Through this tutorial, attendees will gain valuable insights into the challenges and solutions of deploying ML in hardware. Understanding how HLS and hls4ml streamline the development of neural network-based hardware accelerators is fundamental for the industry's future. Participants will learn how these technologies are shaping the future of AI and scientific computing.

Di Guglielmo, Giuseppe [Fermilab]↗
DOE OSTI2025

Parallel derivative-free optimization for simulation-based design of behind-the-meter energy systems

In this work, the integrated design and dispatch of behind-the-meter or distributed resources (e.g. stationary battery storage and solar PV generation) is considered. A simulation-based framework is employed, generating high-fidelity results with closed-loop predictive control at a fine resolution, at the expense of high computational cost (several minutes to a few hours per design point). To address this challenge, parallel derivative-free design methods are considered. Four methods are compared, including state-of-the-art surrogate-based methods (Radial-Basis Functions and Gaussian processes) and sampling strategies, an evolutionary-based method, and a simple sequential grid refinement method. As a case study, two types of design problem with increasing complexity are considered, namely, the design of behind-the-meter resources (three design variables) and the inclusion of grid capacity (four design variables). The second yields a constrained design problem for which violations can only be determined after solving the computationally expensive simulation. For the three-dimensional case, all methods present a good performance, achieving a solution within 1% of the optimum after the first iteration, with the sequential grid refinement exhibiting the fastest convergence and achieving the best final objective value. This indicates that the parallel evaluation of multiple sampling points may be more important than the choice of method for small decision spaces. For the four-dimensional constrained case, the Genetic Algorithm presents the best tradeoff between performance and computational effort, while the rough objective function terrain generated by constraint violation penalties reduces the performance of surrogate-based methods. Contour plots with flat regions indicate flexibility in the optimal design and highlight the importance of characterizing the solution space.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Discovery of Powerful Multivelocity Ultrafast Outflows in the Starburst Merger Galaxy IRAS 05189–2524 with XRISM

We observed the X-ray-bright ultraluminous infrared galaxy IRAS 05189−2524 with XRISM during its performance verification phase. The unprecedented energy resolution of the onboard X-ray microcalorimeter revealed complex spectral features at ∼7–9 keV, which can be interpreted as blueshifted Fe XXV/XXVI absorption lines with various velocity dispersions, originating from ultrafast outflow (UFO) components with multiple bulk velocities of ∼0.076c, ∼0.101c, and ∼0.143c. In addition, a broad Fe–K emission line was detected around ∼7 keV, forming a P Cygni profile together with the absorption lines. The onboard X-ray CCD camera revealed a 0.4–12 keV broadband spectrum characterized by a neutrally absorbed power-law continuum with a photon index of ∼2.3 and intrinsic flare-like variability on timescales of ∼10 ks, both of which are likely associated with near-Eddington accretion. We also found potential variability of the UFO parameters on a timescale of ∼140 ks. Using these properties, we propose new constraints on the outflow structure and suggest the presence of multiple outflowing regions on scales of about tens to 100 Schwarzschild radii, located within roughly 2000 Schwarzschild radii. Since both the estimated momentum and energy outflow rates of the UFOs exceed those of galactic molecular outflows, our results indicate that powerful, multivelocity UFOs are already well developed during a short-lived evolutionary phase following a major galaxy merger, characterized by intense starburst activity and likely preceding the quasar phase. This system is expected to evolve into a quasar, sustaining strong UFO activity and suppressing star formation in the host galaxy.

Astronomy and AstroPhysics↗
DOE OSTI2026

Measurements of three-flavor neutrino oscillations from a PISCES two-detector fit to the NOvA Experiment data

NOvA is a long-baseline neutrino oscillation experiment with two functionally identical detectors: a Near Detector (ND) at Fermilab, placed 1 km from the neutrino source, and a Far Detector (FD) located 810 km away from the ND in Minnesota. NOvA s primary physics goals are to measure the neutrino oscillation parameters $\theta_{23}$ and $\Delta m^2_{32}$ with high precision, determine the neutrino mass hierarchy, and constrain the value of $\delta_{CP}$, primarily via the study of muon neutrino to electron neutrino oscillation. Extracting values for oscillation parameters from fits to data usually relies on treating systematic uncertainties as nuisance parameters, a strategy that suffers from poor scalability as the number of uncertainties becomes larger. This work introduces PISCES (Parameter Inference with Systematic Covariance and Exact Statistics), a novel method that circumvents this scalability problem by encoding systematic uncertainties into a covariance matrix. PISCES utilizes a nested minimization in which optimal systematic pulls are first computed using the covariance matrix in an inner minimization step, then the oscillation parameters are profiled over in the outer minimization. PISCES also uses a Poisson Likelihood term, making it ideal for the inclusion of low-statistic samples in the fits. PISCES is a flexible framework that also supports complex fits, such as a joint Near and Far detector fit. In the standard NOvA analysis, oscillation parameters are extracted using an extrapolation technique in which the ND data indirectly constrain the FD prediction via a ratio method. PISCES, on the other hand, enables a simultaneous ND+FD fit, allowing the high-statistics ND data to directly constrain systematic uncertainties across all samples. This thesis presents the full PISCES joint ND+FD fit for the NOvA three-flavor analysis, details its implementation, and evaluates its performance through extensive robustness tests and fake data studies. It also provides a comparison between the PISCES joint ND+FD results and the standard NOvA extrapolation method using the full NOvA 10-year data set. The results demonstrate that PISCES can successfully fit NOvA data while incorporating the constraints from the ND detectors consistently, using physically motivated systematic uncertainties to account for data/MC discrepancies.

Rajaoalisoa, Miriama [Cincinnati U.]↗
DOE OSTI2024

A Gaussian process based surrogate approach for the optimization of cylindrical targets

Simulating direct-drive inertial confinement experiments presents significant computational challenges, both due to the complexity of the codes required for such simulations and the substantial computational expense associated with target design studies. Machine learning models, and in particular, surrogate models, offer a solution by replacing simulation results with a simplified approximation. In this study, we apply surrogate modeling and optimization techniques that are well established in the existing literature to one-dimensional simulation data of a new cylindrical target design containing deuterium–tritium fuel. These models predict yields without the need for expensive simulations. We find that Bayesian optimization with Gaussian process surrogates enhances sampling efficiency in low-dimensional design spaces but becomes less efficient as dimensionality increases. Nonetheless, optimization routines within two-dimensional and five-dimensional design spaces can identify designs that maximize yield, while also aligning with established physical intuition. Optimization routines, which ignore constraints on hydrodynamic instability growth, are shown to lead to unstable designs in 2D, resulting in yield loss. However, routines that utilize 1D simulations and impose constraints on the in-flight aspect ratio converge on novel cylindrical target designs that are stable against hydrodynamic instability growth in 2D and achieve high yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Machine Learning for Mapping Multipactor Susceptibility in RF Systems: Capabilities and Generalization Constraints

Multipactor is a surface-driven electron avalanche phenomenon that degrades the performance and reliability of radio-frequency (RF) systems in particle accelerator and vacuum electronics applications. Multipactor behavior in a given device structure is conventionally assessed through susceptibility charts, which provide a parameter-space characterization of the instability. In this work, we assess the capabilities of machine-learning (ML) models to learn and predict such susceptibility charts and analyze the constraints governing their generalization across materials. Using a simulation-derived dataset spanning six distinct secondary-electron-yield material profiles in a canonical two-surface planar geometry, we train supervised regression models and artificial neural networks to predict the time-averaged electron growth rate, δavg, across the relevant parameter space. Model performance is evaluated using metrics that explicitly probe the structure of susceptibility charts, including Intersection over Union, Structural Similarity Index, and correlation analysis. Tree-based ensemble models outperform neural-network models in reconstructing susceptibility regions and in generalizing across material domains. Principal-component analysis reveals disjoint material feature distributions, indicating that the piecewise mode structure of multipactor susceptibility is difficult to represent with a single global model and that generalization is constrained by data coverage rather than by model complexity. An exhaustive reduced-coverage study further shows that sparse material-space coverage can yield mean performance in the same general range but producing large variability in the susceptibility-region overlap. These results clarify the capabilities of ML-based surrogate models for parameter-space characterization of multipactor discharge. They also provide guidance for their appropriate use in RF system design.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Scan‐Path‐ and Initial‐State‐Dependent Superdomain Switching in (111)‐Oriented PZT

Polarization switching in ferroelectric materials arises from the collective evolution of complex domain hierarchies, yet deterministic control over these processes remains challenging. Here, we investigate scan-path- and initial-state-dependent switching in epitaxial (111)-oriented PbZr 0.2 Ti 0.8 O 3 thin films using automated AFM-based writing combined with quantitative 3D piezoresponse force microscopy. We show that the scan trajectory acts as an experimentally accessible control parameter for superdomain formation. Box-in-box raster scans reproducibly stabilize ordered stripe superdomains with a reduced subset of symmetry-allowed variants, whereas spiral trajectories generate frustrated mixed-variant states with a broader distribution of final microstructures. Automated pulsing experiments further show that the local superdomain configuration at the nucleation site strongly influences the final written morphology. Phase-field modeling qualitatively reproduces the contrast between representative initial-state geometries and supports the role of compatibility constraints among competing ferroelastic pathways. These findings establish scan-path and initial-state engineering as practical handles to program ferroic order in hierarchical ferroelectric domain structures.

Vasudevan, Rama K. [Oak Ridge National Laboratory ↗
DOE OSTI2025

Requirements on the gain calibration for LiteBIRD polarisation data with blind component separation

The detection of primordial B modes of the cosmic microwave background (CMB) could provide information about the early stages of the Universe's evolution. The faintness of this signal requires exquisite calibration accuracy and control of instrumental systematic effects which otherwise could bias the measurements. In this work, we study the impact of an imperfect relative polarisation gain calibration on the recovered value of the tensor-to-scalar ratio r for the LiteBIRD experiment, through the application of the blind Needlet Internal Linear Combination (NILC) foreground-cleaning method. We derive requirements on the relative calibration accuracy of the overall polarisation gain (Δg ν ) for each LiteBIRD frequency channel. Our results show that minimum variance techniques, as NILC, are less sensitive to systematic gain calibration uncertainties compared to a parametric approach, if the latter is not equipped with a proper modelling of these instrumental effects. In this study, the most stringent requirements are found in the channels where the CMB signal is relatively brighter, with the tightest constraints at 166 GHz (Δg ν ≈ 0.16%). This differs from the outcome of an analogous analysis performed with a parametric method, where the tightest requirements are obtained for the foreground-dominated channels. Gain calibration uncertainties, corresponding to the derived requirements, are then simultaneously propagated into all frequency channels. By doing so, we find that the overall impact on estimated r is lower than the total gain systematic budget for LiteBIRD approximately by a factor 5, due to the correlations of the impacts of gain calibration uncertainties in different frequency channels. In order to decouple the systematic effect from the specific choice of the model, we derive the requirements assuming constant spectral parameters for the foreground emission. To assess the robustness of the obtained results against more realistic scenarios, we repeat the analysis assuming sky models of intermediate and high complexity. In these further cases, we adopt an optimised NILC pipeline, called the Multi-Clustering NILC (MC-NILC). We find that the impact of gain calibration uncertainties on r is lower than the LiteBIRD gain systematics budget for the intermediate-complexity sky model. For the high-complexity case, instead, it would be necessary to tighten the requirements by a factor 1.8.

79 ASTRONOMY AND ASTROPHYSICS↗