Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stellar magnetic fields”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Star-spots and magnetism: testing the activity paradigm in the Pleiades and M67

We measure star-spot filling fractions for 240 stars in the Pleiades and M67 open star clusters using APOGEE high-resolution H-band spectra. For this work, we developed a modified spectroscopic pipeline which solves for star-spot filling fraction and star-spot temperature contrast. We exclude binary stars, finding that the large majority of binaries in these clusters (80 per cent) can be identified from Gaia DR3 and APOGEE criteria – important for field star applications. Our data agree well with independent activity proxies, indicating that this technique recovers real star-spot signals. In the Pleiades, filling fractions saturate at a mean level of 0.248 ± 0.005 for active stars with a decline at slower rotation; we present fitting functions as a function of Rossby number. In M67, we recover low mean filling fractions of 0.030 ± 0.008 and 0.003 ± 0.002 for main sequence GK stars and evolved red giants, respectively, confirming that the technique does not produce spurious spot signals in inactive stars. Star-spots also modify the derived spectroscopic effective temperatures and convective overturn time-scales. Effective temperatures for active stars are offset from inactive ones by −109 ± 11 K, in agreement with the Pecaut & Mamajek empirical scale. Star-spot filling fractions at the level measured in active stars changes their inferred overturn time-scale, which biases the derived threshold for saturation. Finally, we identify a population of stars statistically discrepant from mean activity–Rossby relations and present evidence that these are genuine departures from a Rossby scaling. Our technique is applicable to the full APOGEE catalogue, with broad applications to stellar, galactic, and exoplanetary astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Black hole to breakout: 3D GRMHD simulations of collapsar jets reveal a wide range of transients

We present a suite of the first 3D GRMHD collapsar simulations, which extend from the self-consistent jet launching by an accreting Kerr black hole (BH) to the breakout from the star. We identify three types of outflows, depending on the angular momentum, l, of the collapsing material and the magnetic field, B, on the BH horizon: (i) subrelativistic outflow (low l and high B), (ii) stationary accretion shock instability (SASI; high l and low B), (iii) relativistic jets (high l and high B). In the absence of jets, free-fall of the stellar envelope provides a good estimate for the BH accretion rate. Jets can substantially suppress the accretion rate, and their duration can be limited by the magnetization profile in the star. We find that progenitors with large (steep) inner density power-law indices (≳ 2), face extreme challenges as gamma-ray burst (GRB) progenitors due to excessive luminosity, global time evolution in the light curve throughout the burst and short breakout times, inconsistent with observations. Our results suggest that the wide variety of observed explosion appearances (supernova/supernova + GRB/low-luminosity GRBs) and the characteristics of the emitting relativistic outflows (luminosity and duration) can be naturally explained by the differences in the progenitor structure. Our simulations reveal several important jet features: (i) strong magnetic dissipation inside the star, resulting in weakly magnetized jets by breakout that may have significant photospheric emission and (ii) spontaneous emergence of tilted accretion disc-jet flows, even in the absence of any tilt in the progenitor.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Radiation GRMHD Models of Accretion onto Stellar-mass Black Holes. I. Survey of Eddington Ratios

We summarize results from a survey of radiation-dominated black hole accretion flows across a wide range of mass accretion rates, as well as two values of black hole spin and initial magnetic field geometry. These models apply an algorithm targeting direct solutions to the radiation transport equation in full general relativity and have been enabled by access to modern exascale computing systems. Super-Eddington accretion flows form geometrically thick radiation-pressure-supported disks that drive powerful equatorial outflows. A narrow funnel-shaped photosphere in the inner region results in very low radiative efficiencies in this regime. The structure of near- and sub-Eddington accretion depends on whether there is net vertical magnetic flux at the midplane of the disk. With net flux, the disk forms a thin, dense layer at the midplane surrounded by a magnetically dominated corona, whereas without net flux, the disk remains magnetically dominated everywhere. Although none of our models achieve the magnetically arrested disk regime, those with net vertical flux and a rapidly spinning black hole still produce powerful relativistic jets. Our calculations adopt simple opacity models (with scalings appropriate to stellar-mass black hole accretion). We discuss the application of our results to observations of X-ray binaries and ultraluminous X-ray sources such as Cyg X-3 and SS433. We also speculate on the application of our super-Eddington models to the interpretation of little red dots recently discovered by JWST.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Distinguishability of Radially Directed Energy and Particle Flux in the TJ-II Heliac: Controlling Transport by Feedback-Phasing Experiments

This project supported a U.S. contribution to the TJ-II stellarator research program focused on distinguishing radially directed particle transport from energy transport in magnetized fusion plasmas. The motivating physics issue is that turbulent particle flux and turbulent heat or energy flux are not necessarily locked together: changes in density, electron temperature, plasma potential, and electric-field fluctuations can produce different phase relationships and therefore different radial fluxes. Resolving these relationships is important for understanding improved-confinement behavior, transport-barrier-like regimes, and the broader ability to predict and control confinement in stellarators and related toroidal devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Integrating novel stellarator single-stage optimization algorithms to design the Columbia stellarator experiment

Abstract The Columbia Stellarator eXperiment (CSX), currently being designed at Columbia University, aims to test theoretical predictions related to QA plasma behavior, and to pioneer the construction of an optimized stellarator using three-dimensional, non-insulated high-temperature superconducting (NI-HTS) coils. The magnetic configuration is generated by a combination of two circular planar poloidal field (PF) coils and two 3D-shaped interlinked (IL) coils, with the possibility to add windowpane coils to enhance shaping and experimental flexibility. The PF coils and vacuum vessel are repurposed from the former Columbia Non-Neutral Torus experiment, while the IL coils will be custom-wound in-house using NI-HTS tapes. To obtain a plasma shape that meets the physics objectives with a limited number of coils, novel single-stage optimization techniques are employed, optimizing both the plasma and coils concurrently, in particular targeting a tight aspect ratio QA plasma and minimized strain on the HTS tape. Despite the increased complexity due to the expanded degrees of freedom, these methods successfully identify optimized plasma geometries that can be realized by coils meeting engineering specifications. This paper discusses the derivation of the constraints and objectives specific to CSX, and describe how two recently developed single-stage optimization methodologies are applied to the design of CSX. A set of selected configurations for CSX is then described in detail.

Baillod, A. (ORCID:0000000303529180)↗
DOE OSTI2025

Tokamak Disruption Simulation (Final Technical Report)

The Tokamak Disruption Simulation collaboration was a SciDAC (Scientific Discovery through Advanced Computing) program led by Xianzhu Tang of the Los Alamos National Laboratory but with separate grants to each institution. I was the PI for Columbia University and the only person at Columbia directly funded by this SciDAC. Disruptions are a sudden loss of confinement in tokamaks, and a reliable method of preventing disruptions must be developed before tokamak fusion power plants become feasible. A related plasma confinement concept, the stellarator, is the only way known way of preventing disruptions with the required reliability. The largest devices, both existing and under construction, for confining plasmas for fusion using magnetic fields are tokamaks, so it is important to solve their disruption issues. The sudden loss of confinement in tokamaks can be from a bit of material falling into the plasma, which causes the plasma energy to be quickly lost by radiation, or by an instability causing a breakup of the confining surfaces formed by the magnetic field. This breakup can occur more than a million times faster than one would naively think possible. The sudden breakup of surfaces is called a fast magnetic reconnection and also results in a sudden cooling of the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Power-transfer and fixed-point analysis of sawtooth simulations of a current-carrying stellarator

Power-transfer and fixed-point analysis of previous NIMROD simulations (Roberds et al. , Phys. Plasmas , vol. 23, issue 9, 2016, 092513) improved the understanding of the effect of 3D (non-axisymmetric equilibrium) magnetic fields on sawtooth oscillations in the Compact Toroidal Hybrid (CTH) experiment. Computing the locations of order-1 fixed points, their Greene's residues, and local values for the rotational transform results in a description of CTH sawteeth consistent with Kadomtsev's model. A power-transfer analysis quantifies the distribution of energy among toroidal Fourier modes and their nonlinear interactions. The Lorentz power transfer drives sawtooth growth, and it is unambiguously interpreted as the flow of energy from toroidal mode $n'$ to mode $n$ , catalysed by $\boldsymbol {B}_{n-n'}$ . It has been reported previously that the CTH sawtooth frequency increases with the 3D field strength. This is attributed to an increased growth rate of the internal kink that drives sawtooth oscillations. Here, 3D fields remove energy from the kink, eliminating the possibility that these fields are an additional energy source that drives growth. Instead, 3D fields catalyse energy transfer from large-to-small scales, where magnetic reconnection is stronger. It is proposed that this energy transfer increases the reconnection rate at small scales, which is consistent with the increased growth rate observed at higher 3D field strengths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Stellarator coil optimization supporting multiple magnetic configurations

Abstract We present a technique that can be used to design stellarators with a high degree of experimental flexibility. For our purposes, flexibility is defined by the range of values the rotational transform can take on the magnetic axis of the vacuum field while maintaining satisfactory quasisymmetry. We show that accounting for configuration flexibility during the modular coil design improves flexibility beyond that attained by previous methods. Careful placement of planar control coils and the incorporation of an integrability objective enhance the quasisymmetry and nested flux surface volume of each configuration. We show that it is possible to achieve flexibility, quasisymmetry, and nested flux surface volume to reasonable degrees with a relatively simple coil set through an NCSX-like example. This example coil design is optimized to achieve three rotational transform targets and nested flux surface volumes in each magnetic configuration larger than the NCSX design plasma volume. Our work suggests that there is a tradeoff between flexibility, quasisymmetry, and volume of nested flux surfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Improvement in the spatial resolution of heavy ion beam probe measurements through application of ion optics

A technique for more accurately modeling and improving the spatial resolution of heavy ion beam probe measurements is described. We use a set of particle trajectories to numerically determine the focusing properties of a complicated three-dimensional magnetic field and characterize those properties with a transfer matrix. We then modify the transfer matrix approach of traditional ion optics to include a parameter that describes the ionization location of the detected ions. The ion optics model calculated using this technique enables a more accurate description of the particle trajectories than previously feasible. The model also allows one to easily determine an initial beam focus that could be used during experimental operation to optimize the spatial resolution of measurements. The technique has been applied to the design of a heavy ion beam probe diagnostic for the Wendelstein 7-X stellarator and improvements in the modeled spatial resolution by a factor of about 2 over previous estimates are possible. As a result, the improved spatial resolution will enable measurements of plasma fluctuations with smaller wavelengths than would otherwise be possible.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Construction of quasisymmetric stellarators using a direct coordinate approach

In this work, optimized stellarator configurations and their analytical properties are obtained using a near-axis expansion approach. Such configurations are associated with good confinement as the guiding center particle trajectories and neoclassical transport are isomorphic to those in a tokamak. This makes them appealing as fusion reactor candidates. Using a direct coordinate approach, where the magnetic field and flux surface functions are found explicitly in terms of the position vector at successive orders in the distance to the axis, the set of ordinary differential equations for first and second order quasisymmetry is derived. Examples of quasi-axisymmetric shapes are constructed using a pseudospectral numerical method. Finally, the direct coordinate approach is used to independently verify two hypotheses commonly associated with quasisymmetric magnetic fields, namely that the number of equations exceeds the number of parameters at third order in the expansion and that the near-axis expansion does not prohibit exact quasisymmetry from being achieved on a single flux surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A tutorial on radio frequency sheath physics for magnetically confined fusion devices

Radio frequency (RF) sheaths occur under a wide variety of conditions when RF waves, material surfaces and plasma coexist. RF sheaths are of special importance in describing the interaction of ion cyclotron range of frequency (ICRF) waves with the boundary plasma in tokamaks, stellarators and other magnetic confinement devices. In this article the basic physics of RF sheaths is discussed in the context of magnetic fusion research. Techniques for modelling RF sheaths, their interaction with RF wave fields and the resulting consequences are highlighted. The article is intended as a guide for the early-career ICRF researcher, but it may equally well serve to provide an overview of basic RF sheath concepts and modelling directions for any interested fusion scientist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effects of drifts on scrape-off layer transport in W7-X

Abstract Drifts affect particle, momentum, and energy transport in the scrape-off layer (SOL) of tokamaks and stellarators, altering plasma flows and creating asymmetries between divertors. To understand how drifts affect SOL transport in the W7-X island divertor, an experiment was performed to compare plasmas with matched core parameters but opposite magnetic field directions, and therefore opposite drift transport directions. Parallel flow measurements made with coherence imaging spectroscopy are interpreted with the aid of a diagnostic forward model and a 1D simple SOL model that includes the E × B drift. In low-density plasmas ( n ‾ e < 2 × 10 19 m −3 ), the poloidal E × B drift induces a large poloidal density asymmetry within the island SOL, as measured by divertor Langmuir probes. This in turn causes the parallel flow stagnation point to shift from the position halfway between targets to the X-point in the drift direction, leading to near-unidirectional flow throughout the SOL. As density increases, the effects of the poloidal E × B drift decrease substantially, resulting in a smaller density asymmetry and the development of a counter-streaming flow pattern. For the entire density range probed in this experiment ( n ‾ e = 1.5 − 6 × 10 19 m −3 ), the experimental observations are more consistent with the effects of the poloidal E × B drift than the radial E × B drift.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Physics Basis, Optimization, and Control for Integrated 3D Edge Long-pulse Tokamak Scenarios

The aim of this project with international collaborations (PI: J.-K. Park of Princeton Plasma Physics Laboratory) is to develop a unified physics basis and predictive capability for the control of edge-localized modes (ELMs) with optimized non-axisymmetric (3D) fields, by leveraging the unique research capabilities of international KSTAR tokamak in Korea. At the University of California, Irvine, we use a fusion simulation code GTC to study the effects of the 3D magnetic fields (including both magnetic islands and stochastic magnetic field lines) on microturbulence. We focus on the effects of 3D equilibrium on turbulent and neoclassical transport that leads to ELM suppression. During the project period, the UCI team of this collaboration has made significant progress in understanding turbulent transport in 3D geometry including both tokamaks with RMP and stellarators. Four UCI postdoctoral researchers, Drs. Javier Nicolau, Gyungjin Choi, Pengfei Liu, and Xishuo Wei were partially supported by this project and collaborating with the entire project team. We have published 4 journal papers and gave several invited talks. This project advance physics understand of 3D effects on microturbulence, specially effects of magnetic islands. It will contribute to the DOE SciDAC ISEP project in the area of GTC development for general 3D geometry capability, rigorous verification and validation, and readiness for the emerging exascale computing. Through collaborations with computational scientists, the project helps to advance high performance computing in the area of GPU acceleration, linear solver, and data management.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Mapping the space of quasisymmetric stellarators using optimized near-axis expansion

A method is demonstrated to rapidly calculate the shapes and properties of quasi-axisymmetric and quasi-helically symmetric stellarators. In this approach, optimization is applied to the equations of magnetohydrodynamic equilibrium and quasisymmetry, expanded in the small distance from the magnetic axis, as formulated by Garren & Boozer [ Phys. Fluids B, vol. 3, 1991, p. 2805]. Due to the reduction of the equations by the expansion, the computational cost is significantly reduced, to times of the order of 1 cpu second, enabling wide and high-resolution scans over parameter space. In contrast to traditional stellarator optimization, here, the cost function serves to maximize the volume in which the expansion is accurate. A key term in the cost function is $\| \boldsymbol {\nabla }\boldsymbol B \|$ , the norm of the magnetic field gradient, to maximize scale lengths in the field. Using this method, a database of $5\times 10^5$ optimized configurations is calculated and presented. Quasisymmetric configurations are observed to exist in continuous bands, varying in the ratio of the magnetic axis length to average major radius. Several qualitatively new types of configuration are found, including quasi-helically symmetric fields in which the number of field periods is two or more than six.

Landreman, Matt (ORCID:000000027233577X)↗
DOE OSTI2021

White paper on NIF Discovery Science and frontier regimes of HEDP

Over the past two decades, experimental capabilities at high energy density (HED) facilities such as the National Ignition Facility (NIF), Omega, and Omega EP lasers, and the Z pulsed power facility are able to make precision measurements of fundamental quantities such as equations of state (EOS), material phase, opacities, and ionization levels Z. Flow based quantities are also being studied, such as hydrodynamic instabilities, evolution into turbulence, turbulent dynamo magnetic field amplification, magnetic reconnection; plasma instabilities, plasma wakefield particle acceleration; and relativistic plasma phenomena, such as (e + ,e - ) pair creation. Experimental regimes can be created that reproduce conditions relevant to planetary and exoplanet interiors, stellar interiors, the interiors of brown dwarfs, and the envelope conditions of white dwarf stars. Diagnostic techniques have been developed hand in hand with these new experimental regimes that can be accessed. In this white paper, a brief summary of key experimental achievements from the NIF Discovery Science (basic science) program will be presented; followed by a discussion of current challenges, needs, and potential paths forward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

MPEX AI Digital Twins

All magnetically confined plasma fusion power plant concepts (Tokamak, Spherical Tokamak, Stellarator, Mirror, ...) must exhaust the heat and plasma from the core confinement region to the material walls. The primary channel for this exhaust is through a plasma divertor which directs plasma along open magnetic field lines to a material target. The Material Plasma Exposure eXperiment (MPEX) illustrated in Figure 1, is a high-power, steady-state linear plasma device designed to produce the plasma material interaction (PMI) conditions of the divertor of future magnetic confinement fusion power plants: energy flux 20MW/m 2 , ion fluence 1031/m 2 , pulse duration 106 sec. These goals of plasma exposure in MPEX are well beyond those achieved in magnetic fusion experimental devices. Successfully achieving these high power steady state conditions for long pulses requires operational control of the heating and particle sources and the plasma flux to the walls and target. The MPEX AI Hot Spot Controller, proposed in this project, will help achieve the operational milestones of MPEX. The MPEX device will begin commissioning at the end of FY26. A smaller proto-MPEX was operated for 14,666 plasma discharges and will resume operation in September of 2025 as proto-MPEX-lite, with reduced capability, to test a new window for the Helicon plasma source. The proto-MPEX data has undergone surrogate modeling with machine learning methods (R. Archibald, 2022 IEEE International Conference on Big Data). This proto-MPEX data will be used to begin development of the AI digital twins described in this white paper. The scientific mission of MPEX is to qualify materials of different composition for use in the high energy and plasma flux conditions of a fusion power plant. The materials exposed in MPEX will in some cases be exposed to high neutron fluxes at other ORNL facilities to measure the changes to their PMI properties. The targets exposed in MPEX will be transported under vacuum to a Surface Analysis Station (SAS). The SAS will be equipped with the following diagnostics: Focused Ion Beam (FIB) for trench milling, 100-400 angstrom resolution scanning electron microscope (SEM), surface mapping x-ray spectrometer, high resolution camera, and a future upgrade to a laser induced breakdown spectroscopy quadruple mass spectrometer (LIBS-QMS). The MPEX experiments will generate diverse pre- and post-exposure measurement data of detailed material properties down to the crystal grain level in 3D for post-exposure assessment of PMI damage (e.g. cracking, melting, erosion and redeposition of the material). Physics models for the PMI, and how the material composition and manufacturing impact its performance under high energy plasma exposure, need to be validated with MPEX data to guide the selection of new candidate materials. Our vision for the MPEX AI Digital Twins project is to supply experimental and physics model simulation data to train Artificial Intelligence (AI) models for data processing, analysis, operational control, PMI and materials simulation to maximize the scientific output of the MPEX device. Ultimately, an AI digital twin of MPEX material assessment metrics for tested and synthetic material types with simulated PMI will be trained by the AI Modeling Teams on the experimental and physics simulation data submitted to the American Science Cloud by this project. A purely empirical search for the best material is inefficient given the finite number of samples that can be tested on MPEX. In order to expand the material properties database for training the MPEX Material Assessment AI Digital Twin, and to gain physics understanding of the PMI processes, physics models of the material properties and PMI processes are required. The physics simulations provide detailed simulation data, like impact angles for plasma ions, sputtering yields, transport of the ionized sputtered target material in the plasma, and redeposition locations. This simulation data expands the measurement data for deeper physics understanding. The experimental data is essential to validate the PMI and material structure simulation models. The validated models can then be used to generate new simulation data of MPEX material assessments for synthetic material compositions that have not been exposed in MPEX. These predictive simulations, plus the whole experimental dataset, will be used to train the MPEX Material Assessment AI Digital Twin allowing a rapid generative AI search for new materials with reduced PMI damage by interpolating the domain of the training set. These new optimum materials can be simulated with the physics codes and/or tested in MPEX. The ability of AI neural networks to interpolate multi-dimensional parameter spaces and generate virtual data is exploited for a more efficient search for optimum materials. The advent of the Transformational AI Models Consortium (TAIMC) is an opportunity to engage with state of the art private and public AI developers to achieve the goals of the AI digital twins and AI accelerated physics models proposed in this project. Our partners at ORNL from the Advance Scientific Computing Research (ASCR) organization will collaborate in accelerating the integrated plasma material interaction simulation framework. This simulation framework will provide a platform for generating simulation data across a range of physical fidelities, including hybrid methods that produce multi-fidelity results. This data will be leveraged for AI model development, both for generation of surrogates and the automation of simulation campaigns. A part of the research below will include collaborative efforts with the TAIMC to (i) adapt data storage approaches to ensure AI-readiness, (ii) provide a protypical exemplar to inform and exercise constructed workflows, and (iii) generate and share data, using the TAIMC unified AI data standard, for foundational models that will be trained from multiple sources across the DOE complex. We will also collaborate with the TAIMC, as well as the planned AI modeling teams, to develop approaches for reducing the cost of data generation. These include tailored multi-fidelity approaches as well as fine-tuning strategies to augment general, large-scale foundational models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Core plasma fueling by fast inward particle transport after hydrogen pellet injection in Wendelstein 7-X

A large database of more than 1000 individual cryogenic hydrogen pellets injected into Wendelstein 7-X for plasma fueling was analyzed to improve the understanding of the three phases of the process: the ablation, deposition and transport of the pellet material. Kilohertz-sampled electron density and temperature measurements revealed a more complex drift behavior than predicted by numerical code simulation. It could be explained by the poloidal plasma E r x B- drift rotation, which plays a significant role in stellarators, but was not previously considered in pellet injection codes like HPI2. The drift results in a fast poloidal rotation of the pellet material around the plasma core, leading to an almost homogeneous deposition over the involved flux surfaces regardless of magnetic high and low field side injection geometry. Additionally, a novel fast inward directed transport mechanism (‘FIT-effect’) was observed. The effect occurs on timescales of tens of milliseconds and cannot be explained by neoclassical transport or diffusion. It might be linked to the turbulence pinch recently found in Wendelstein 7-X. When the FIT-effect occurs, the pellet particles are rapidly transferred from the deposition flux surfaces to the plasma core, causing the plasma density profile to peak, which is beneficial for confinement in Wendelstein 7-X. The large pellet injection database was statistical analyzed with regard to pellet and plasma parameters, which delivered some starting points towards developing an understanding of the physics behind the FIT-effect. The results indicate, that plasma core fueling via pellet injection is largely independent of the injection geometry in stellarators under certain conditions, reducing the technical complexity of the injection system.

Wendelstein 7-X↗
DOE OSTI2026

Probing the axion-electron coupling with NuSTAR observations of galaxies

We search for the existence of ultralight axions coupling to electrons and photons using data from the NuSTAR telescope directed toward the galaxies M82, M87, and M31. We focus on electron bremsstrahlung and Compton scattering for axion production in stars, summing over the stellar populations found in the target galaxies when computing the axion luminosity. We then compute the hard x-ray signal that arises from the conversion of these axions to photons in each galaxy’s magnetic fields, inferred from analog galaxies in cosmological magnetohydrodynamic simulations. Analyzing NuSTAR data toward these galaxies between roughly 20 to 70 keV, we find no evidence for axions and set leading constraints on the combined axion-electron and axion-photon coupling at the level of |𝑔 𝑎⁢𝑒⁢𝑒 ×𝑔 𝑎⁢𝛾⁢𝛾 | ≲ 8.3 ×10 −27 GeV −1 for 𝑚 𝑎 ≲ 10 −10 eV at 95% confidence, with M82 providing the most stringent constraints.

Axion-like particles↗