Search NASASearch

SEARCH · Search NASA

Results for “m:N operations”

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 19 records

High photon-phonon pair generation rate in a two-dimensional optomechanical crystal

Integrated optomechanical systems are a leading platform for manipulating, sensing, and distributing quantum information, but are limited by residual optical heating. Here, we demonstrate a two-dimensional optomechanical crystal (OMC) geometry with increased thermal anchoring and a mechanical mode at 7.4 GHz, well aligned with the operation range of cryogenic microwave hardware and piezoelectric transducers. The eight times better thermalization than current one-dimensional OMCs, large optomechanical coupling rates, g 0 /2π ≈ 880 kHz, and high optical quality factors, Q opt = 2.4 × 10 5 , allow ground-state cooling (n m = 0.32) of the acoustic mode from 3 K and entering the optomechanical strong-coupling regime. In pulsed sideband asymmetry measurements, we show ground-state operation (n m < 0.45) at temperatures below 10 mK, with repetition rates up to 3 MHz, generating photon-phonon pairs at ≈ 147 kHz. Our results extend optomechanical system capabilities and establish a robust foundation for future microwave-to-optical transducers with entanglement rates exceeding state-of-the-art superconducting qubit decoherence rates.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Quantum Enhanced Dark-Matter Search with Entangled Fock States in High-Quality Cavities

We present a quantum-enhanced protocol for detecting wave-like dark matter using an array of $N$ entangled superconducting cavities initialized in an $m$-photon Fock state. By distributing and recollecting the quantum state with an entanglement-distribution operation, the scan rate scales as $N^2(m+1)$ while thermal excitation is the dominant background, significantly outperforming classical single-cavity methods under matched conditions. We evaluate the robustness of our scheme against additional noise sources, including decoherence and beamsplitter infidelity, through theoretical analysis and numerical simulations. In practice, the key requirements, namely high-Q superconducting radio-frequency cavities that support long integration times, high-fidelity microwave beamsplitters, and universal cavity control, are already available on current experimental platforms, making the protocol experimentally feasible.

Freiman, Benjamin [Fermilab]

Integrated modeling of advanced tokamaks for tearing mode avoidance on DIII-D

Understanding the physics of the onset and evolution of tearing modes (TMs) in tokamak plasmas is crucial for high-performance steady-state operations. The onset condition of (m,n) = (2,1) tearing stability, $Δ′$ > $Δ_c^′$> 0, has been studied in DIII-D steady-state hybrid discharges with accurate equilibrium reconstruction and well-measured plasma profiles. Here, m is the poloidal mode number and n is the toroidal mode number. The onset of n = 1 TMs is carefully identified using the Mirnov array to determine in-phase signals for the onset time and using electron cyclotron emission measurements to trace back from the saturated island width in the experiments for both initial island width and its growth rate. The tearing stability index $Δ′$ at TM onset is calculated in multiple ways, using the modified Rutherford equation and using the PEST-III or resistive DCON magnetohydrodynamic codes, and the results are in reasonable agreement with each other. The calculated $Δ′$ is compared with an analytical formula of the tearing stability threshold $Δ_c^′$ to determine the condition for the mode onset, $Δ′$ > $Δ_c^′$. Finally, the local gradient effect of the plasma current and the pressure profiles on the tearing stability is investigated in terms of the difference between $Δ′$ and $Δ_c^′$ through a series of modeled equilibria with the parametric variations. A database of experimental mode instability onsets is compatible with the simulated regions of instability obtained from the parametric variations.

Linear stability analysis

Insights into Probing the Effect of Molecular Weight of Poly(carboxybetaine methacrylate) on the Performance of Forward Osmosis Desalination

Zwitterionic polymers have proven to be a promising nonfouling material that can be applied in the design of selective layers of thin film composite (TFC) membranes. Extending the permeability and usage of TFC membranes have attracted increasing interest in membrane-based desalination processes since water-flux reduction associated with biofouling nowadays persists as a common challenge. By virtue of its strong hydration, this polymer category is very useful to counteract biofouling in marine and biomedical systems, but the benefits from their application in membrane technology are still emerging. The efficacy of the nonfouling property as a function of the polymer’s molecular weight remains unknown. In pursuit of that vision, this study fosters new scientific insights via probing different molecular weights of poly(carboxybetain methacrylate) (PCBMA) coated on the surface as a selective layer for the prepared TFC membranes. The coated zwitterionic membranes (zM) exhibited excellent performance in preventing water flux decay in a bench-scale forward osmosis system. The prepared zM membranes revealed enhanced hydrophilic properties and retained their operational water flux when compared to the control. Our results suggest that using an intermediate-size molecular weight (PCBMA M n 50,000) will result in the best operational performance. The intermediate size resulted in the lowest flux decline rate (R t ) of 0.01 ± 0.001 (zM-50) when compared to the unmodified control membrane 0.56 ± 0.071 (M0) after using a model BSA foulant solution. Furthermore, all coated membranes exhibited similar trends in the observed reverse salt flux profiles, as well. In conclusion, the constructed zM membranes will serve as a model to develop further selective layers in the construction of TFC membranes.

36 MATERIALS SCIENCE

Ferroelectric epitaxial Al(Sc/B)N/Mo/SiC heterostructures for high operating temperature devices

Advances in wurtzite nitride ferroelectrics of Al 1-x M x N (⁠ M = Sc or B) have led to novel capabilities, which must be integrated into existing fabrication processes. In the case of electronics operating > 200°C, a movement toward SiC-based platforms enables better performance over conventional Si, primarily due to the higher bandgap and lower intrinsic carrier concentration of SiC. Hence, the challenge is to develop a deposition process to integrate Al 1-x M x M x N ferroelectrics with elevated temperature-compatible materials for high temperature electronics such as the non-volatile memory component. We demonstrate epitaxial Al M N/Mo/SiC heterostructures, which provides both the crystalline and surface features that promote high-quality ferroelectric nitride film growth. Omega scans of the Mo (110) reflection exhibit a full width at half max of < 0.02° (40 arc sec) and the (0002) peak of the subsequently grown nitride film had a value of 1.1° for 160 nm thick Al 0.7 Sc 0.3 N and 1.3° for 400 nm thick Al 0.94 B 0.06 N. The crystallographic relationships found between the Al M N, Mo, and SiC layers indicate an advancement in sputter deposition of epitaxial films. Ferroelectric switching is also shown at 400 °C in both samples via polarization-electric field hysteresis and pulsed measurements, which exhibited P r values >100 μC cm -2 and E c between 3 and 4 M V cm -1 ⁠, despite a large presence of oxygen in both Al M N films ranging between 4 and 5 at.%, revealed by compositional analysis. This study demonstrates the process for synthesizing high-crystal quality ferroelectric nitride films, which can be used in extremely high temperature applications.

74 ATOMIC AND MOLECULAR PHYSICS

MHD-induced beta limits in the Large Helical Device

Using the extended-magnetohydrodynamics code, M3D-C1, we perform a systematic numerical study of the effect of externally applied heating on the achievable plasma beta in a ten field-period heliotron. Heat sources of varying intensity are applied to a vacuum magnetic field that is representative of the standard configuration of the Large Helical Device, with R 0 = 3.66 m, where R 0 is the radial position of the magnetic axis in vacuum. As the system is driven to a state that is unstable to low-n magnetohydrodynamic (MHD) modes, nonlinear mode interactions lead to the formation of chaotic magnetic fields. With sufficiently strong heating, a collapse of the electron temperature profile is observed. This demonstrates the necessity of simulating the self-consistent evolution of plasma profiles, without imposing assumptions on the structure of the magnetic field, to accurately determine transport properties in stellarator plasmas. It also highlights the value of these advanced simulation capabilities for accelerating the development of high-performance stellarator operating scenarios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Accelerating GNNs on GPU Sparse Tensor Cores through N:M Sparsity-Oriented Graph Reordering

Recent advancements in GPU hardware support have introduced the capability to leverage N:M sparse patterns for substantial performance gains. Graphs in Graph Neural Networks (GNNs) are typically sparse, but the sparsity is often irregular, not conforming to such sparse patterns. In this paper, we propose a novel graph reordering algorithm, the first of its kind, to reshape irregular graph data into the N:M structured sparse pattern at the tile level, allowing linear-algebra-based graph operations in GNNs to benefit from the N:M sparse hardware. The optimization is lossless, maintaining the accuracy of GNN. It can remove 98-100\% violations of the N:M sparse patterns at the vector level, and increase the proportion of conforming graphs in SuiteSparse collection from 5-9\% to 88.7-93.5\%. On A100 GPUs, the optimization accelerates Sparse Matrix Matrix (SpMM) by up to 43X (2.3X -- 7.5X on average) and speeds up the key graph operations in GNNs on real graphs by as much as 8.6X (3.5X on average).

artificial intelligence, graph neural networks

Trap-assisted Auger-Meitner recombination in GaN p-i-n diodes

Most properties of semiconductor devices are dominated by shallow impurities. However, deep defects often play an important role, for instance, in recombination processes or high field transport. While a variety of techniques are available to assess the density and energy levels of impurities, other properties, such as the recombination mechanisms of the defects, escape observation. We report on the direct measurement of hot electrons generated by trap-assisted Auger-Meitner recombination (TAAR) in GaN p-i-n diodes. By performing electron emission spectroscopy (EES) on diodes with surfaces activated to negative electron affinity by cesium, we observe the expected overflow electrons of p-i-n diodes under low current injection. However, when operating the devices at higher current densities, as low as ∼25 A/c⁢m 2 , we measure the emission of high-energy electrons. At variance with the observed hot electrons in light-emitting diodes (LEDs) using EES, the hot electrons generated in p-i-n diodes at our tested currents cannot be from eeh Auger-Meitner recombination due to the diodes' significantly lower carrier densities compared to those in LEDs. During our measurements, we observe the emission of accumulated electrons with energies ∼0.42 eV, ∼0.99 eV, ∼1.43 eV, and ∼2.32 eV above the conduction-band minimum (CBM) at various bias conditions, suggesting the existence of conduction-band features in GaN at these energies where electrons can be long-lived, such as satellite-valley minima and inflection points. We also measure incompletely relaxed hot electrons approaching energies 1.97 ± 0.13 eV and 2.94 ± 0.13 eV above the CBM, as the diodes are biased to high currents, suggesting at least some of the TAAR partaking defects have an energy level ≳1.97 eV and ≳2.94 eV away from either the conduction or valence band edges. Additionally, at our highest operating currents, we measure hot electrons with energies 3.28 ± 0.13 eV above the CBM, providing direct evidence of TAAR processes involving shallow impurities. This unexpected observation of TAAR in GaN p-i-n diodes spotlights the importance of further studies of defects in GaN and the necessity to incorporate the multi-phonon emission, radiative, and TAAR capture steps of defect-assisted recombination cycles into device modeling. Furthermore, this experiment demonstrates the applicability of the simplest semiconductor structures, p-i-n diodes, as a test bed to study the rich recombination physics of semiconductor materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Modelling and experiment to stabilize disruptive tearing modes in the ITER baseline scenario in DIII-D

Abstract The achievement of high gain, stationary conditions in a tokamak scenario aimed at producing fusion energy in the ITER Project is crucial to the demonstration that this form of energy can be used in future reactors to provide cheap and clean energy globally. Disruptions are a challenge for the fusion energy field, in particular for the ‘ITER Baseline Scenario’ (IBS), as reproduced in the DIII-D tokamak. This work shows that a solution has been found for the m = 2/ n = 1 tearing modes that have consistently caused disruptions in the IBS: stable operation down to zero input torque was achieved by modifying the current density profile at the beginning of the pressure flattop and the ELM character later in the discharges, guided by previous results showing that the most likely cause of these instabilities is the current density profile. The coupling between sawteeth, n>2 modes and the 2/1 TMs is shown to not be statistically significant, nor the leading origin for the evolution towards instability. Ideal and resistive MHD modeling provide positive verification that a steeper ‘well’ in the region of the q = 2 rational surface leads to worse ideal stability, higher tearing index Δ’ and lower threshold Δ’ c for resistive instabilities, consistent with the experimental results. This provides confidence that the methods used in this work can be extrapolated to other devices and applied to avoid disruptions in ITER and pulsed fusion devices worldwide.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Operation above the Greenwald density limit in high performance DIII-D negative triangularity discharges

The density limit in strongly-shaped negative triangularity (NT) discharges is studied experimentally in the DIII-D tokamak. Record-high Greenwald fractions f G are obtained, using gas puff injection only, with values up to near 2, where f G is defined as the ratio of the line-averaged density over n G = I p /(π a 2 ), with I p [MA] the plasma current and a[m] the plasma minor radius. A clear higher operational limit with higher auxiliary power is also demonstrated, with the ohmic density limit about two times lower than with additional neutral beam injection heating. The evolution of the electron density, temperature and pressure profiles are analyzed as well. The core density can be up to twice the Greenwald density and keeps increasing, while the value at the separatrix remains essentially constant and slightly below n G . The edge temperature gradient collapses to near zero and NT plasmas are shown to be resilient to such profiles in terms of disruptivity. We also present the time evolution of the inverse electron pressure scale length with the value at the last closed flux surface (LCFS) decreasing below the value at the normalized radius 0.9 near the density limit, demonstrating the clear drop of confinement starting from the edge. This inverse scale length ‘collapse’ at the LCFS also defines well the characteristic behavior of the kinetic profiles approaching a density limit.

density limit

Determination of confinement regime boundaries via separatrix parameters on Alcator C-Mod based on a model for interchange-drift-Alfvén turbulence

The separatrix operational space (SepOS) model (Eich and Manz 2021 Nucl. Fusion 61 086017) is shown to predict the L–H transition, the L-mode density limit, and the ideal magnetohydrodynamic ballooning limit in terms of separatrix parameters for a wide range of Alcator C-Mod plasmas. The model is tested using Thomson scattering measurements across a wide range of operating conditions on C-Mod, spanning $\overline{n}_{\mathrm{e}}$ = 0.3–5.5 $\times 10^{20}\,$m −3 , $B_{\mathrm{t}} = 2.5$–8.0 T, and $B_{\mathrm{p}}$= 0.1–1.2 T. An empirical regression for the electron pressure gradient scale length, $\lambda_{{p}_{\mathrm{e}}}$, against a turbulence control parameter, $\alpha_{\mathrm{t}}$, and the poloidal fluid gyroradius, $\rho_{\mathrm{s,p}}$, is constructed for H-modes and found to require positive exponents for both regression parameters, indicating turbulence widening of near-scrape-off layer widths at high $\alpha_{\mathrm{t}}$ and an inverse scaling with $B_{\mathrm{p}}$, consistent with results on ASDEX Upgrade. The SepOS model is also tested in the unfavorable drift direction and found to apply well to all three boundaries, including the L–H transition as long as a correction to the Reynolds energy transfer term, $\alpha_\mathrm{RS} \lt 1$ is applied. I-modes typically exist in the unfavorable drift direction for values of $\alpha_{\mathrm{t}} \lesssim 0.3$. Finally, an experiment studying the transition between the Type-I ELMy and EDA H-mode is analyzed using the same framework. It is found that a recently identified boundary $\alpha_{\mathrm{t}} = 0.55$ at excludes most EDA H-modes but that the balance of wavenumbers responsible for the L-mode density limit, namely $k_\mathrm{EM} = k_\mathrm{RBM}$, may better describe the transition on C-Mod. The ensemble of boundaries validated and explored is then applied to project regime access and limit avoidance for the SPARC primary reference discharge parameters.

ELM suppression

Use of differential plasma rotation to prevent disruptive tearing mode onset from 3-wave coupling

Abstract Plasma differential rotation is found to be capable of preventing disruptive neoclassical tearing modes (NTMs) seeded by nonlinear three-wave coupling. As tearing modes degrade confinement and can lead to disruptions, stabilization strategies are crucial to the successful operation of future devices. In ITER-relevant scenarios on DIII-D, rotationally coupled m / n = 1/1 and 3/2 modes have been observed to drive 2/1 islands through three-wave coupling. The frequency of the driven 2/1 mode is set by matching conditions and the frequencies of the driving modes. When the driven mode frequency matches the local plasma rotation frequency, e.g. at low differential rotation, the driven 2/1 island can grow into a disruptive NTM. Using neutral beam torque as an actuator to scan the differential rotation, these experiments demonstrate that a sufficiently large frequency mismatch prevents destabilization of disruptive 2/1 NTMs by three-wave coupling. This work indicates that differential rotation can be used as an actuator to prevent NTMs seeded by three-wave coupling.

Richner, N. J. (ORCID:0000000155443915)

The root cause of disruptive NTMs and paths to stable operation in DIII-D ITER baseline scenario plasmas

Analyses of the DIII-D ITER Baseline Scenario database support that the disruptive m,n=2,1 magnetic islands are pressure gradient driven, non-linear instabilities seeded in a sequence of stochastic transient magnetic perturbations, and that the current profile relaxation does not affect the m,n=2,1 island onset rate. At low torque, these Neoclassical Tearing Modes are most commonly seeded by non-linear 3-wave coupling when the differential rotation between the q=1 & q=2 rational surfaces approaches zero. Lack of statistically significant difference between the current profiles of stable and unstable states, as well as lack of correlation between the tearing mode onset rate and the current profile relaxation both reject causality between the current profile evolution and the 2,1 magnetic island onsets in these plasmas. These support that preserving the differential rotation between the q=1 and q=2 rational surfaces is key to long pulse stable operation in the plasma scenario planned for ITER, while optimization of the current profile within the explored parameter space may lead to much weaker improvements than sustaining the differential rotation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Nuclear Data Sensitivity/Uncertainty Studies of Tantalum-Reflected Systems

Tantalum (Ta) is a refractory metal with high melting point and high corrosion resistance. This makes it useful as a mold material for plutonium casting operations. Unfortunately, there are few criticality benchmarks with high nuclear data sensitivities. Given that, it is desirable to design new experiments to provide additional validation data for nuclear criticality safety, which is the objective of the Thales project. Multiple application models were used to represent process conditions (either normal conditions or credible upset conditions) both at Los Alamos National Laboratory (LANL) and Savannah River Site (SRS). A total of 40 application models were obtained. Given this, it was useful to down-select the number of needed application models. This work documents the down-selection based on several nuclear data sensitivity/uncertainty metrics. In addition to showing the result of the specific down-selection used for the Thales project, this work discusses potential ways to represent a large number of application models with a smaller number of proposed experiment designs.

42 ENGINEERING

BSEC flux towers: CSAT3B and TRH

The data were collected as part of the BSEC project, during the period from June 2025 to May 2026. Directory "broadway" contains data collected on a multi-level flux tower (US-BWf) in the Broadway East neighborhood (1808 North Patterson Park Ave., Baltimore City, MD 21213; LAT: 39o18'40.31'' N; LONG: 76o35'12.43'' W). At each of the four measurement heights (8.5 m, 11.1 m, 13.4 m, 15.9 m), a Campbell Scientific CSAT3B sonic anemometer was operated at 50 Hz to measure virtual temperature (tc) and three velocity components (u: 270 degrees; v: 180 degrees; w: vertical), and a RM Young temperature sensor (model 41382VC) was operated at 1 Hz inside a compact aspirated radiation shield (model 43502) to measure absolute temperature (T) and relative humidity (RH). Inside directory "broadway", directory "netcdf" contains data collected each day in 5-minute chunks that have been converted to NetCDF format (before quality checking), while "4hr" contains data arranged into 4-hour chunks (also in NetCDF format) that have been through basic quality checking steps (treating data points with nonzero diagnostic codes as missing data; fixing six or fewer consecutive missing data points using linear interpolation). Users are recommended to start with data in directory "4hr", while data in directory "netcdf" can be used for reference purposes.

Baltimore

Strong interaction physics at the luminosity frontier with 22 GeV electrons at Jefferson Lab

Here, the purpose of this document is to outline the developing scientific case for pursuing an energy upgrade to 22 GeV of the Continuous Electron Beam Accelerator Facility (CEBAF) at the Thomas Jefferson National Accelerator Facility (TJNAF, or JLab). This document was developed with input from a series of workshops held in the period between March 2022 and April 2023 that were organized by the JLab user community and staff with guidance from JLab management (see Sec. 10). The scientific case for the 22 GeV energy upgrade leverages existing or already planned Hall equipment and world-wide uniqueness of CEBAF high-luminosity operations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Initial measurement of reactor antineutrino oscillation at SNO+

The SNO+ collaboration reports its first spectral analysis of long-baseline reactor antineutrino oscillation using 114 tonne-years of data. Fitting the neutrino oscillation probability to the observed energy spectrum yields constraints on the neutrino mass-squared difference Δ$m^2_{21}$. In the ranges allowed by previous measurements, the best-fit Δ$m^2_{21}$ is ($8.85^{+1.10}_{-1.33}$) × 10 -5 eV 2 . This measurement is continuing in the next phases of SNO+ and is expected to surpass the present global precision on Δ$m^2_{21}$ with about three years of data. In the following, we first describe the configuration of the SNO+ detector when it was partially filled with scintillator, and a characterization of the detector response using intrinsic radioactivity. Next, we detail the event selection and expectations for reactor IBDs and (α, n) reactions. Then, we present the results of an energy spectrum analysis using 114 tonne-years of data. We conclude with prospects of future results from the SNO+ detector, which has been operating fully-filled with 780 tonnes of scintillator

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Efficient First-Order Algorithms for Large-Scale, Non-Smooth Maximum Entropy Models with Application to Wildfire Science

Maximum entropy (MaxEnt) models are a class of statistical models that use the maximum entropy principle to estimate probability distributions from data. Due to the size of modern data sets, MaxEnt models need efficient optimization algorithms to scale well for big data applications. State-of-the-art algorithms for MaxEnt models, however, were not originally designed to handle big data sets; these algorithms either rely on technical devices that may yield unreliable numerical results, scale poorly, or require smoothness assumptions that many practical MaxEnt models lack. In this paper, we present novel optimization algorithms that overcome the shortcomings of state-of-the-art algorithms for training large-scale, non-smooth MaxEnt models. Our proposed first-order algorithms leverage the Kullback–Leibler divergence to train large-scale and non-smooth MaxEnt models efficiently. For MaxEnt models with discrete probability distribution of n elements built from samples, each containing m features, the stepsize parameter estimation and iterations in our algorithms scale on the order of O(mn) operations and can be trivially parallelized. Moreover, the strong ℓ1 convexity of the Kullback–Leibler divergence allows for larger stepsize parameters, thereby speeding up the convergence rate of our algorithms. To illustrate the efficiency of our novel algorithms, we consider the problem of estimating probabilities of fire occurrences as a function of ecological features in the Western US MTBS-Interagency wildfire data set. Our numerical results show that our algorithms outperform the state of the art by one order of magnitude and yield results that agree with physical models of wildfire occurrence and previous statistical analyses of wildfire drivers.

Physics