Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fixed target experiments”

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 145 records · Page 8

Multiphysics Simulations of Nanofibrous High Power Targets for the Inference of Failure Modes

Fermilab s High Power Targetry Research and Development (HPT R&D) group have been developing and studying an electrospun nanofiber target concept to support the need for robust targets in future fixed target facilities. These nanofiber mats have demonstrated resistance to radiation damage, and the free motion of the individual fibers is expected to mitigate the cyclic stresses induced by a pulsed, high-power beam. To evaluate the efficacy of this concept, nanofiber mat samples have been sent by the HPT R&D group to the HiRadMat facility at CERN for prototypic thermal shock testing; the outcomes of these experiments show that the survivability of a nanofiber target depends on its construction parameters, in particular the packing density of the fibers. Samples with higher packing densities have consistently been destroyed by exposure to the HiRadMat beam, with a hole at the beam center visible, and layers of nanofibers peeled away from the center hole, whereas samples with lower densities have survived with limited damage. The exact reason for the failure of the higher density targets was unclear at the time of the original experiments, but the results of our recent multiphysics simulations which recreate the experiments support the hypothesis that the expansion and pressurization of the air inside the target after being heated by the pulsed beam is the cause; in a high density nanofiber mat, the motion of air through the pores of the mat is much more restricted, and induces a larger pressure on the fibers, blowing the mat apart. In this talk, we ll share the results of these simulations and discuss how they support this hypothesis.

Asztalos, Will↗

LHCspin: a Polarized Gas Target for LHC

The goal of the LHCspin project is to develop innovative solutions for measuring the 3D structure of nucleons in high-energy polarized fixed-target collisions at LHC, exploring new processes and exploiting new probes in a unique, previously unexplored, kinematic regime. A precise multi-dimensional description of the hadron structure has, in fact, the potential to deepen our understanding of the strong interactions and to provide a much more precise framework for measuring both Standard Model and Beyond Standard Model observables. This ambitious task poses its basis on the recent experience with the successful installation and operation of the SMOG2 unpolarized gas target in front of the LHCb spectrometer. Besides allowing for interesting physics studies ranging from astrophysics to heavy-ion physics, SMOG2 provides an ideal benchmark for studying beam-target dynamics at the LHC and demonstrates the feasibility of simultaneous operation with beam-beam collisions. With the installation of the proposed polarized target system, LHCb will become the first experiment to simultaneously collect data from unpolarized beam-beam collisions at $\sqrt{s}$=14 TeV and polarized and unpolarized beam-target collisions at $\sqrt{s_{NN}}\sim$100 GeV. LHCspin has the potential to open new frontiers in physics by exploiting the capabilities of the world's most powerful collider and one of the most advanced spectrometers. This document also highlights the need to perform an R&D campaign and the commissioning of the apparatus at the LHC Interaction Region 4 during the Run 4, before its final installation in LHCb. This opportunity could also allow to undertake preliminary physics measurements with unprecedented conditions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measurement of the mass-changing, charge-changing, and production cross sections of C 11 , B 11 , and B 10 nuclei in C 12 + p interactions at 13.5 GeV / c per nucleon

A good knowledge of nuclear fragmentation cross sections is important to interpret the fluxes of secondary cosmic rays from the Galaxy. Here we report new measurements of nuclear fragmentation with the NA61/SHINE experiment at the CERN SPS. The specific focus is on cross sections important for the production of boron in the Galaxy from the interactions of C 12 nuclei with hydrogen in the interstellar medium, including the contribution from the decay of the short-lived C 11 fragments. The data were taken with the secondary C 12 beam at beam momentum of 13.5 GeV / c per nucleon and two fixed targets, polyethylene ( CH 2 ) and graphite (C), from which we derive the cross sections of carbon on hydrogen. We present the measurement of the fragmentation cross sections of C 11 , B 11 , and B 10 as well as the mass- and charge-changing cross sections.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

White Box Access to Quantum Testbeds for Co-Design

At Lawrence Livermore National Laboratory (LLNL), we operate and maintain the Quantum Device and Integration Testbed (QuDIT) facility, a small quantum testbed that supports about 10 active research teams (including our own) and over 50 internal and external collaborators. This testbed is designed to give remote white box access to users for research, training, and outreach. A guiding principle behind the development of our testbed infrastructure, software and user interfaces is to empower users to perform experiments at the cutting edge of quantum information science at any level of abstraction, from materials studies, device physics and control and characterization techniques to algorithm development and quantum operating system design. Our testbed targets a multilevel quantum system (qudit) to expand the accessible Hilbert space of a simple-to-manufacture quantum device and focuses on quantum simulation, typically implemented through custom gates designed with quantum optimal control methods, rather than on a universal computing framework with a fixed gate set. We leverage the Lab’s high-performance computing (HPC) program and related expertise to simulate quantum systems, develop hybrid algorithms, and generate gates optimized for given simulations. Additionally, we have adopted a co-design philosophy from the HPC community in designing new hardware, so that the systems we develop are optimized for the specific physics simulations we plan to use them for.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

National Ignition Facility. Facility and Infrastructure Systems Maintenance Plan

Ensuring the reliability of the NIF, including its support systems, laser systems, target diagnostic systems, and utilities, is essential to the availability of the NIF in its support of NNSA missions. NIF is a key capability in the DOE Stockpile Stewardship Program and supports high energy physics experiments for nuclear weapons, energy, and astrophysics applications. High system reliability provides opportunities for shots and scientific discoveries with opportunities to enhance and upgrade capabilities. This Maintenance Plan (MP) identifies the policies and procedures used to perform and support asset management of the NIF Facility and Infrastructure Systems (FInS), NIF Lasers & Alignment (LASE), NIF Target Experimental Operations (TOPS), NIF Target Area Science and Engineering (TASE), and NIF&PS Control Systems (NCS). The FInS systems include the facility, HVAC, contamination control, beampath, and Line Replaceable Units (LRUs) as well as utilities which create the beampath environments, such as vacuum, argon, or clean dry air. The LASE systems are Programmatic systems which include laser diagnostics, alignment, power conditioning, pulsed power, and input laser systems. The TOPS and TASE systems are also Programmatic systems which include target and diagnostic delivery systems and positioners, many different insertable and fixed target diagnostics, and cryogenic and target gas fill systems. Finally, the NCS systems include both software and hardware for industrial and shot operation control systems. Policies governing administrative and operational practices related to maintenance of FInS, LASE, TOPS, TASE, and NCS systems are described in this plan. In addition, the plan provides processes and procedures for managing, tracking, and documenting the work. This document, the NIF Operations Management Plan, NIF-5020544 (Ref. 1), and NIF Shot Operations Plan, NIF-5018506 (Ref. 2), together satisfy the requirements of the Conduct of Operations. Duties, responsibilities, and reporting requirements of the various positions associated with FInS, LASE, and TOPS maintenance are detailed in this plan. The FInS systems include both Real Property systems with asset management requirements specified in DOE Order 430.1C (Ref. 3) and Programmatic systems. In addition, for FInS, there is a list of the System Level Maintenance Plans (SLMPs) in NIF-1007419198 (Ref. 4) which provide the system descriptions and maintenance plan and schedule. In addition, the list includes the Reliability Centered Maintenance (RCM) and Experience Centered Maintenance (ECM) evaluations that have been performed for applicable FInS systems as well as reliability criticality per Section 3.5. The -AM version of the NIF Maintenance Plan focuses on the reliability program for FInS, LASE, TOPS, TASE, and NCS within the context of the overall NIF Reliability, Availability, and Maintainability (RAM) program and incorporates changes since the -AL version from August 2011 and has been updated to be fully consistent with the updates to Ref. 1. It also includes asset management considerations, updates to the Work Order (WO) process within the NIF Computerized Maintenance Management System (CMMS) which is EAM infor® System Maintenance and Reliability Tracking (SMaRT) (Ref. 5), and updates to metrics and key performance indicators (KPIs).

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Timelike Compton Scattering from a longitudinally polarised target with CLAS12 at Jefferson Lab

Explorations into the internal dynamics of hadrons are constantly evolving, and the requirement for experimental results to verify theoretical models of hadron structure is paramount. A key area in this field is the study of Generalised Parton Distributions (GPDs), which are functions used to model the momenta of quarks and gluons within hadrons, and the methods to access GPDs experimentally. One such scattering process that allows access to these is Timelike Compton Scattering. TCS complements existing Deeply Virtual Compton Scattering experiments and allows investigation into the universality of GPDs through access to the real and imaginary parts of the parton helicity independent GPD Hq via beam spin asymmetries (BSA), and it provides novel access to the real and imaginary parts of the parton helicity dependent GPD ˜Hq through target polarisation asymmetries (TSA). This thesis work presents a comparative study with the first published BSA for TCS at the Thomas Jefferson National Accelerator Facility (JLab), alongside a first time extraction of a Target Spin Asymmetry with the Summer 2022 data taking run. JLab hosts the Continuous Electron Beam Accelerator Facility (CEBAF) which provides a 12 GeV electron beam to four experimental halls. Hall-B contains the CEBAF Large Acceptance Spectrometer, which took data across three run periods on a longitudinally polarised NH3 and ND3 fixed target from 2022-2023, to extract measurements of electron-proton scattering, from which a TCS signal could be extracted. The thesis discusses work done to understand and eliminate contributions from the non-/low-polarised nuclear background, testing pre-established cuts to eliminate pion background from a dilepton (e+e-) final state and modifying them as needed for the new experimental run, and attempts to hone in on a clean TCS signal from which to extract the two asymmetry observables. A comparison with existing BSA results was performed; however, the statistical errors are too large to draw a significant conclusion as to whether there is agreement across each bin. More data is needed for a multidimensionally binned extraction. A proof of principle was achieved in the TSA measurements, with two out of four kinematic bins showing preliminary agreement in shape with theoretical values. Again, the errors are significant due to the contributions from the nuclear background. To support these conclusions, a further study was done, which takes into account an estimate of the asymmetries with the full available dataset (this thesis is based only on data taken in the summer set; at the time of writing processing was still being conducted for the final two datasets), as well as an estimate including additional future experiment days that were awarded in July 2024. Additional work was done on a secondary project exploring the feasibility of measuring TCS at the upcoming Electron Ion Collider, supporting the design proposal for the detector for the first interaction region and giving a positive outlook for the future of these types of measurements beyond JLab.

Gates, Kayleigh [Univ. of Glasgow, Scotland (Unite↗

The ePix10k 2-megapixel hard X-ray detector at LCLS

The ePix10ka2M (ePix10k) is a new large area detector specifically developed for X-ray free-electron laser (XFEL) applications. The hybrid pixel detector was developed at SLAC to provide a hard X-ray area detector with a high dynamic range, running at the 120 Hz repetition rate of the Linac Coherent Light Source (LCLS). The ePix10k consists of 16 modules, each with 352 × 384 pixels of 100 µm × 100 µm distributed on four ASICs, resulting in a 2.16 megapixel detector, with a 16.5 cm × 16.5 cm active area and ∼80% coverage. The high dynamic range is achieved with three distinct gain settings (low, medium, high) as well as two auto-ranging modes (high-to-low and medium-to-low). Here the three fixed gain modes are evaluated. The resulting dynamic range (from single photon counting to 10000 photons pixel −1 pulse −1 at 8 keV) makes it suitable for a large number of different XFEL experiments. The ePix10k replaces the large CSPAD in operation since 2011. The dimensions of the two detectors are similar, making the upgrade from CSPAD to ePix10k straightforward for most setups, with the ePix10k improving on experimental performance. The SLAC-developed ePix cameras all utilize a similar platform, are tailored to target different experimental conditions and are designed to provide an upgrade path for future high-repetition-rate XFELs. Here the first measurements on this new ePix10k detector are presented and the performance under typical XFEL conditions evaluated during an LCLS X-ray diffuse scattering experiment measuring the 9.5 keV X-ray photons scattered from a thin liquid jet.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Measurements of Electron Beam-Induced Spin-Relaxation in Frozen-Spin Hydrogen-Deuteride (HD)

Targets of solid hydrogen-deuteride (HD) can be prepared in a frozen-spin state with spin relaxation times (T1) in excess of a year. The present work has studied their potential use in experiments with electron beams. Polarization relaxation rates have been measured during exposure to sub-nanoAmp currents of 9.7 MeV electron beams from the recently commissioned Upgraded Injector Test Facility (UITF) at Jefferson Lab (JLab). These UITF measurements can be used to anticipate the expected performance at the GeV energies typical of JLab experiments, since the energy deposition in a target is almost independent of beam energy. An in-beam dilution refrigerator equipped with superconducting solenoids has been used to maintain solid HD samples at ?0.1 K and ?1 tesla, and internal NMR coils have been used to monitor hydrogen polarization. Spin-relaxation in two frozen-spin targets, with initial H-polarizations of 40% and 34%, have been tracked while exposed to beams under varying conditions of current, dose, beam duty factor and temperature. At a fixed accumulated dose, the spin-relaxation rates drop with current, suggesting depolarization by the charge cloud of the beam. Imposing a duty factor on the otherwise continuous UITF beam, with forced millisecond scale gaps, has shown no obvious correlation to polarization loss rates. After an accumulated dose of ?2 ?C/cm2, beam-off spin-relaxation rates drop from their immeasurably long pre-irradiation values to the order of weeks (with beam-on T1 values typically an order of magnitude shorter), reflecting a buildup of paramagnetic charge centers within the HD lattice. The accumulated polarization loss was approximately proportional to dose in both targets, dropping to 1/e of their initial values after ?6 ?C/cm2. Thermal equilibrium polarizations of targets not in the frozen-spin state (with intentionally short T1) have been used to deduce the in situ temperature of solid HD while under electron bombardment. A model for depolarization by beam-associated paramagnetic impurities largely accounts for the data, and suggests that improvements in heat removal could lead to significant increases in the in-beam T1.

O'Connell, Thomas↗

Power handling in a highly-radiative negative triangularity pilot plant

Abstract This work explores detailed power handling solutions for a class of high-field, highly-radiative negative triangularity (NT) reactors based around the MANTA concept (Rutherford et al 2024 Plasma Phys. Control. Fusion ). The divertor design is kept as simple as possible, opting for a standard divertor with standard leg length. FreeGS is used to create an equilibrium for the boundary region, prioritizing a short outer leg length of only ∼50 cm (∼40% of the minor radius). The UEDGE code package is used for the boundary plasma solution to track plasma temperatures and fluxes to the divertor targets. It is found that for P SOL = 25 MW and n sep = 0.96 × 10 20 m −3 , conditions consistent with initial core transport modeling, little additional power mitigation is necessary. For a fixed impurity fraction of just 0.13% Ne in the plasma, the peak heat flux density at the more heavily loaded outer targets falls to 7.8 MW m −2 , while the electron temperature T e remains just under 5 eV. Scans around the parameter space reveal that even at densities lower than in the primary operating scenario, P SOL can be increased up to 50 MW, so long as a slightly higher fraction of extrinsic radiator is used. With less than 1% neon (Ne) impurity content, the divertor still experiences less than 10 MW m −2 at the outer target. Design of the plasma-facing components includes a close-fitting vacuum vessel with a tungsten inner surface as well as FLiBe-carrying cooling channels fashioned into the VV wall directly behind the divertor targets. For the seeded heat flux profile, Ansys Fluent heat transfer simulations estimate that the outer target temperature remains at just below 1550 ∘ C. Initial scoping of advanced divertor designs shows that for an X-divertor, detachment of the outer target becomes much simpler, and plasma fluxes to the targets drop considerably with only 0.01% Ne content.

Miller, M. A. (ORCID:0000000265406533)↗

TAMU TRACER: Targeted Mobile Measurements to Isolate the Impacts of Aerosols and Meteorology on Deep Convection

Difficulty in using observations to isolate the impacts of aerosols from meteorology on deep convection often stems from inability to resolve the spatiotemporal variations in the environment serving as the storm’s inflow region. During the DOE TRacking Aerosol Convection interactions ExpeRiment (TRACER) in June-September 2022, a Texas A&M University (TAMU) team conducted a mobile field campaign to characterize the meteorological and aerosol variability in airmasses that serve as inflow to convection across the ubiquitous mesoscale boundaries associated with the sea- and bay-breezes in the Houston, Texas, region. These boundaries propagate inland over the fixed DOE Atmospheric Radiation Measurement (ARM) sites. However, convection occurs on either or both the continental or maritime sides or along the boundary. The maritime and continental airmasses serving as convection inflow may be quite distinct, with different meteorological and aerosol characteristics that fixed-site measurements cannot simultaneously sample. Thus, a primary objective of TAMU TRACER was to provide mobile measurements similar to those at the fixed sites, but in the opposite airmass across these moving mesoscale boundaries. TAMU TRACER collected radiosonde, lidar, aerosol, cloud condensation nuclei (CCN), and ice nucleating particle (INP) measurements on 29 enhanced operations days covering a variety of maritime, continental, outflow, and pre-frontal airmasses. This paper summarizes the TAMU TRACER deployment and measurement strategy, instruments, available datasets, and provides sample cases highlighting differences between these mobile measurements and those made at the ARM sites. We also highlight the exceptional TAMU TRACER undergraduate student participation in high impact learning activities through forecasting and field deployment opportunities.

54 ENVIRONMENTAL SCIENCES↗

Integrated performance of large HDC-capsule implosions on the National Ignition Facility

In this paper, we report on eight, indirect-drive, deuterium–tritium-layered, inertial-confinement-fusion experiments at the National Ignition Facility to determine the largest capsule that can be driven symmetrically without relying on cross-beam energy transfer or advanced Hohlraum designs. Targets with inner radii of up to 1050 μm exhibited controllable P 2 symmetry, while larger capsules suffered from diminished equatorial drive. Reducing the Hohlraum gas-fill-density from 0.45 mg/cm 3 to 0.3 mg/cm 3 did not result in a favorable shift of P 2 amplitude as observed in preceding tuning experiments. Reducing the laser-entrance-hole diameter from 4 mm to 3.64 mm decreased polar radiation losses as expected, resulting in an oblate symmetry. The experiments exhibited the expected performance benefit from increased experimental scale, with yields at a fixed implosion velocity roughly following the predicted 1D dependence. With an inner radius of 1050 μm and a case-to-capsule-ratio of 3.0, experiment N181104 is the lowest implosion-velocity experiment to exceed a total neutron yield of 10 16 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiscale Richtmyer-Meshkov instability experiments to isolate the strain rate dependence of strength

Theoretical analysis of Richtmyer-Meshkov instability (RMI) experiments for solid strength shows that the strain rate for a given shock should be inversely proportional to the length scale of the sine wave perturbations when η 0 k , the nondimensional amplitude to wavelength ratio, is held fixed. To isolate the effect of strain rate on strength, free-surface RMI specimens of annealed copper were prepared with three perturbation regions with the same η 0 k but different length scales, characterized by the wavelength λ varying by a factor of 4.9 from 65 to 130 to 320 µ m . Three such targets with different fixed η 0 k ′ s were impacted to a shock pressure of 25 GPa, and the instability evolution was measured with photon Doppler velocimetry. Strengths estimated by comparing hydrocode simulation to the data increased from 700 to 1200 MPa as λ decreased. The different η 0 k targets exercised increasing amounts of plastic strain yet showed no evidence of strain hardening. Physical regime sensitivity analysis determined that for 320 − 65 µ m wavelength perturbations, the effective strain rates increased from 8.7 × 10 6 to 3.3 × 10 7 s − 1 , a factor of 3.8. Thus, the predicted strain rate scaling was mostly achieved but slightly suppressed by increased strength at higher rates. The RMI strength estimates were plotted against constitutive testing data on copper from the literature to show striking evidence of the strength upturn at higher strain rates. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗

Stability Constants of Lanthanide-nitrate Complexes in Aqueous Solutions: A Theoretical Study

Calculating stability constants for lanthanide--nitrate complexes in aqueous solution is challenging due to the complex free-energy landscapes of the participating species. In this work, we compare cluster-continuum solvation and condensed-phase approaches using universal machine learning interatomic potentials (MLIPs) for determining the stability constants of lanthanide--nitrate complexes in aqueous solutions. Within the cluster--continuum solvation framework at the B3LYP level of theory, reactions involving lanthanide coordination numbers of both 8 and 9 are found to be relevant. After an empirical linear free-energy correction, the cluster--continuum results fall on the same order-of-magnitude scale as the experimental stability constants. By contrast, MACE-MP0 MLIP underestimates lanthanide hydration numbers and gives PMF-derived stability constants with large deviations from experiment, whereas MACE-MATPES-R2SCAN improves both hydration structure and the stability-constant scale but still does not quantitatively reproduce the detailed lanthanide trend. Across both approaches, nitrate binding is best viewed as a labile coordination motif rather than a fixed mono- or bidentate structure, with hydration-shell structure influencing which configurations are favored. Overall, the cluster--continuum calculations provide a practical semi-quantitative baseline for the experimental stability-constant scale, while the explicit-solvent MLIP benchmarks show clear progress from MACE-MP0 to MACE-MATPES-r2SCAN but also highlight the need for lanthanide-targeted training, fine-tuning, or improved long-range and polarization treatments to obtain predictive thermodynamics for complex aqueous lanthanide chemistry.

Dinpajooh, Mohammadhasan↗

Advanced manufacturing provides tailor-made solutions for crystallography with x-ray free-electron lasers

Serial crystallography at large facilities, such as x-ray free-electron lasers and synchrotrons, evolved as a powerful method for the high-resolution structural investigation of proteins that are critical for human health, thus advancing drug discovery and novel therapies. However, a critical barrier to successful serial crystallography experiments lies in the efficient handling of the protein microcrystals and solutions at microscales. Microfluidics are the obvious approach for any high-throughput, nano-to-microliter sample handling, that also requires design flexibility and rapid prototyping to deal with the variable shapes, sizes, and density of crystals. Here, we discuss recent advances in polymer 3D printing for microfluidics-based serial crystallography research and present a demonstration of emerging, large-scale, nano-3D printing approaches leading into the future of 3D sample environment and delivery device fabrication from liquid jet gas-dynamic virtual nozzles devices to fixed-target sample environment technology.

36 MATERIALS SCIENCE↗

Autonomous Nanoparticle Synthesis Guided by In Situ Multiscale Structural Characterization

Autonomous synthesis platforms promise rapid exploration of vast parameter spaces; yet, integrating in situ structural characterization in closed-loop synthesis optimization remains challenging. We demonstrate a realization of such a closed-loop platform coupled with a droplet-flow microreactor, in situ X-ray scattering methods (SAXS/WAXS), and Gaussian process optimization to synthesize citrate-reduced Au nanoparticles with targeted characteristics. The system efficiently explored ∼19,000 synthesis recipes through 365 experiments, achieving precise control over size (4–60 nm) and polydispersity (σ < 0.11) across large citrate/gold ratios, exceeding traditional synthesis boundaries (1–10). Beyond confirming classical Turkevich–Frens trends, partial-dependence analysis revealed strong nonlinear coupling among precursor, citrate, and pH effects. Combining quantitative SAXS/WAXS analysis with electron microscopy characterization, we uncovered that crystallite size (d c ) and particle size (d) follow d c = 0.18d + β, where synthesis chemistry controls the intercept β while maintaining a universal slope. This parallel-band structure enables independent tuning of crystallite domain size at fixed particle diameter through a combination of chloride, gold precursor, citrate, and pH contributions (cross-validated Spearman ρ = 0.7 ± 0.1). High-resolution electron microscopy shows multiple lattice-fringe orientations within single particles, directly confirming polycrystalline domains and the ability to tune d c at the fixed d. The platform’s validation includes indistinguishable static versus flowing measurements, stable droplet transport at 100 °C, and <5% run-to-run variation, establishing a robust framework for mapping and controlling multiscale nanoparticle structure across expansive chemical spaces. In conclusion, the developed closed-loop platform can be applied to a borad range of nanosyntheis processes.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Low- n stability and plasma response to RMP in various STEP scenarios

The low-n (n is the toroidal mode number) magnetohydrodynamic (MHD) stability and plasma response are numerically investigated for various scenarios designed for STEP, that are relevant for the H-mode pedestal analysis. Control of the edge-localized modes (ELMs) with externally applied resonant magnetic perturbations (RMPs) is considered. Optimization of the ELM control coil current configuration, based on the computed plasma MHD response and well-established figures of merit validated on present-day experiments, finds reasonable robustness of a fixed coil phasing (for a given n-number) to control ELMs in all five STEP plasmas considered. Based on certain semi-empirical criteria, the required coil current to achieve ELM suppression is estimated to be about 10–20 kAt with the n = 1 or 2 RMP configuration and about 100–200 kAt for the n = 4 RMP. Systematic linear stability calculations are used to map out stability windows for the low-n kink-peeling modes, in terms of the ideal-wall location and variation of the edge safety factor q 95 with respect to the target design. The kink-peeling stability boundary is found to be generally sensitive to the q 95 variation, which has implications for achieving the quiescent H-mode regime in STEP. Full toroidal quasilinear initial-value simulations for these STEP plasmas find that generation of the edge-harmonic oscillations (EHOs) depends sensitively on the plasma scenario, the initial linear stability of the kink-peeling modes, the initial plasma toroidal flow and q 95 . In general, it is easier (more robust) to access the EHO-regime for two of the cases considered with smaller plasma volume and higher on-axis safety factor. Finally, quasilinear simulations find robust density pumpout due to applied RMPs in these STEP plasmas, but the effect on the plasma toroidal flow varies among different cases.

EHO↗

Modeling Equilibrium Solid–Liquid Interfaces under Effective Constant Chemical Potential Using Machine Learning Interatomic Potentials

The chemical potential (μ) of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot maintain constant chemical potential with reference to a targeted concentration or chemical potential under nonequilibrium or dynamic conditions, as solute species can migrate to the interface and deplete (or enrich) the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative quasi-constant chemical potential molecular dynamics (iqCμMD) simulation, which helps simulate targeted molar concentrations of species in solution. iqCμMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target bulk concentration (chemical potential), which allows simulation of the interface under the bulk conditions comparable to experiment. We demonstrate our approach using machine learning interatomic potential (MLIP)-based MD simulations of the Na 2 SO 4,aq –graphene interface, and to show the transferability of our approach, we also perform classical force field-based MD simulations of NaCl aq –air and NaCl aq –graphite interfaces, which produce comparable results to previous CμMD simulations. Our results also show that the iqCμMD approach efficiently achieves the desired bulk ion concentration within two iterations, and by utilizing MLIPs, we can achieve converged results using relatively small-scale simulations compared to previous CμMD simulations. By combining iqCμMD with MLIP-driven simulations, solid–liquid interfaces can be modeled under an effective constant chemical potential with DFT-level accuracy. Here, we show that iqCμMD offers a robust and simple computational framework for constant chemical potential simulations, as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.

Chemical structure↗

Probabilistic partition of unity networks for high–dimensional regression problems

We explore the probabilistic partition of unity network (PPOU-Net) model in the context of high-dimensional regression problems and propose a general framework focusing on adaptive dimensionality reduction. With the proposed framework, the target function is approximated by a mixture of experts model on a low-dimensional manifold, where each cluster is associated with a fixed-degree polynomial. We present a training strategy that leverages the expectation maximization (EM) algorithm. During the training, we alternate between (i) applying gradient descent to update the DNN coefficients; and (ii) using closed-form formulae derived from the EM algorithm to update the mixture of experts model parameters. Under the probabilistic formulation, step (ii) admits the form of embarrassingly paralleliazable weighted least-squares solves. The PPOU-Nets consistently outperform the baseline fully-connected neural networks of comparable sizes in numerical experiments of various data dimensions. Here, we also explore the proposed model in applications of quantum computing, where the PPOU-Nets act as surrogate models for cost landscapes associated with variational quantum circuits.

97 MATHEMATICS AND COMPUTING↗