Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model Inference”

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 667 records · Page 37

THE GRISM LENS-AMPLIFIED SURVEY FROM SPACE (GLASS). V. EXTENT AND SPATIAL DISTRIBUTION OF STAR FORMATION IN z ∼ 0.5 CLUSTER GALAXIES

We present the first study of the spatial distribution of star formation in z ∼ 0.5 cluster galaxies. The analysis is based on data taken with the Wide Field Camera 3 as part of the Grism Lens-Amplified Survey from Space (GLASS). We illustrate the methodology by focusing on two clusters (MACS 0717.5+3745 and MACS 1423.8+2404) with different morphologies (one relaxed and one merging) and use foreground and background galaxies as a field control sample. The cluster+field sample consists of 42 galaxies with stellar masses in the range 10{sup 8}–10{sup 11} M{sub ⊙} and star formation rates in the range 1–20 M{sub ⊙} yr{sup −1}. Both in clusters and in the field, Hα is more extended than the rest-frame UV continuum in 60% of the cases, consistent with diffuse star formation and inside-out growth. In ∼20% of the cases, the Hα emission appears more extended in cluster galaxies than in the field, pointing perhaps to ionized gas being stripped and/or star formation being enhanced at large radii. The peak of the Hα emission and that of the continuum are offset by less than 1 kpc. We investigate trends with the hot gas density as traced by the X-ray emission, and with the surface mass density as inferred from gravitational lens models, and find no conclusive results. The diversity of morphologies and sizes observed in Hα illustrates the complexity of the environmental processes that regulate star formation. Upcoming analysis of the full GLASS data set will increase our sample size by almost an order of magnitude, verifying and strengthening the inference from this initial data set.

79 ASTRONOMY AND ASTROPHYSICS↗

Deep-Learning-Based Koopman Modeling for Online Control Synthesis of Nonlinear Power System Transient Dynamics

Power system stability and control have become more challenging due to the increasing uncertainty associated with renewable generation. Here, the performance of conventional control is highly driven by the physics-based offline-developed dynamic models that can deviate from the actual system characteristics under different operating conditions and/or configurations. Data-driven approaches based on online measurements can be a better solution to addressing these issues by capturing real-time operation conditions. This article describes a novel fully data-driven probabilistic framework to derive a linear representation of postcontingency grid dynamics and online prescribe control based on the derived model to enhance transient stability. The complex nonlinear power system dynamics is approximated by a linear model by using multiple neural network modules that infer distributions of the observations and introducing a Koopman layer to sample possible Koopman linear models from the inferred distributions. The trained model features linearity that can be easily incorporated into the existing linear control design paradigm and ease the controller design process. The effectiveness of Koopman-based control designs is validated through comparative case studies, which demonstrate increased prediction accuracy and control performance when applied to a power system with heterogeneous generator dynamics.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Efficient kinetic thermal inverse modeling for organic material decomposition

The prevalent use of organic materials in manufacturing is a fire safety concern, and motivates the need for predictive thermal decomposition models. A critical component of predictive modeling is numerical inference of kinetic parameters from bench scale data. Currently, an active area of computational pyrolysis research focuses on identifying efficient, robust methods for optimization. This paper demonstrates that kinetic parameter calibration problems can successfully be solved using classical gradient-based optimization. We explore calibration examples that exhibit characteristics of concern: high nonlinearity, high dimensionality, complicated schemes, overlapping reactions, noisy data, and poor initial guesses. The examples demonstrate that a simple, non-invasive change to the problem formulation can simultaneously avoid local minima, avoid computation of derivative matrices, achieve a computational efficiency speedup of 10x, and make optimization robust to perturbations of parameter components. Techniques from the mathematical optimization and inverse problem communities are employed. By re-examining gradient-based algorithms, we highlight opportunities to develop kinetic parameter calibration methods that should outperform current methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling of Vapor Bubble Growth Under Nucleate Boiling Conditions in Reduced Gravity

A dynamic model is developed to describe the evolution of a vapor bubble growing at a nucleation site on a superheated surface under arbitrary gravity. The bubble is separated from the surface by a thin microlayer and grows due to the evaporation from the microlayer interface. The average thickness of the microlayer increases as the bubble expands along the surface if the evaporation rate is lower than some critical value. The corresponding threshold value of the surface temperature has to be associated with the burn-out crisis. Two main reasons make for bubble separation, which are the buoyancy force and a force caused by the vapor momentum that comes to the bubble with vapor molecules. The latter force is somewhat diminished if condensation takes place at the upper bubble surface in subcooled liquids. The action of the said forces is opposed by inertia of the additional mass of liquid as the bubble center rises above the surface and by inertia of liquid being expelled by the growing bubble in radial directions. An extra pressure force arises due to the liquid inflow into the microlayer with a finite velocity. The last force helps in holding the bubble close to the surface during an initial stage of bubble evolution. Two limiting regimes with distinctly different properties can be singled out, depending on which of the forces that favor bubble detachment dominates. Under conditions of moderately reduced gravity, the situation is much the same as in normal gravity, although the bubble detachment volume increases as gravity diminishes. In microgravity, the buoyancy force is negligible. Then the bubble is capable of staying near the surface for a long time, with intensive evaporation from the microlayer. It suggests a drastic change in the physical mechanism of heat removal as gravity falls below a certain sufficiently low level. Inferences of the model and conclusions pertaining to effects caused on heat transfer processes by changes in bubble hydrodynamics induced by gravity are discussed in connection with experimental evidence, both available in current and in as yet unpublished literature.

Buyevich, Yu A.↗

Application of X-Ray Pulsar Navigation: A Characterization of the Earth Orbit Trade Space

The potential for pulsars as a navigation source has been studied since their discovery in 1967. X-ray pulsar navigation (XNAV) is a celestial navigation system that uses the consistent timing nature of x-ray photons from milli-second pulsars (MSP) to perform space navigation. By comparing the detected arrival of x-ray photons to a reference database of expected pulsar lightcurve timing models, one can infer a range and range rate measurement based on light time delay. Much of the challenge of XNAV comes from the faint signal, availability, and distant nature of pulsars. This is a study of potential pulsar XNAV measurements to measure extended Kalman filter (EKF) tracking performance with a wide trade space of bounded Earth orbits, using a simulation of existing x-ray detector space hardware. An example of an x-ray detector for XNAV is the NASA Station Explorer for X-ray Timing and Navigation (SEXTANT) mission, a technology demonstration of XNAV set to perform on the International Space Station (ISS) in late 2016early 2017. XNAV hardware implementation is driven by trajectory and environmental influences which add noise to the x-ray pulse signal. In a closed Earth orbit, the radiation environment can exponentially increase the signal noise from x-ray pulsar sources, decreasing the quality and frequency of measurements. The SEXTANT mission in particular improves on the signal to noise ratio by focusing an array of 56 x-ray silicon drift detectors at one pulsar target at a time. This reduces timing glitches and other timing noise contributions from ambient x-ray sources to within a 100 nanosecond resolution. This study also considers the SEXTANT scheduling challenges inherent in a single target observation. Finally, as the navigation sources are now relatively inertial targets, XNAV measurements are also subject to periods of occultation from various celestial bodies. This study focuses on the characterization of these drivers in closed Earth orbits and is not a tuning analysis of the EKF. The study shows that the closed Earth orbit for XNAV performance is reliant on the orbit semi-major axis and eccentricity as well as orbit inclination. These parameters are the primary drivers of pulsar measurement availability and significantly influence the natural spacecraft orbit dynamics. Sensitivity to initial orbit determination error growth due to the scarcity of XNAV measurements within an orbital period require appropriate timing of initial XNAV measurements. The orbit angles of argument of perigee and right ascension of the ascending node, alongside the other orbit parameters, complete the initial cadence of XNAV measurements. The performance of initial XNAV measurements then propagates throughout the experimental period. The study provides a basis to missions who wish to consider XNAV as a potential navigation source in a closed Earth orbit design. It provides a family of orbit trajectories as well as other modeling considerations needed to effectively evaluate if XNAV is an effective navigation source for a potential mission. As an EKF is sensitive to a linearized estimated state, this study has a direct benefit of providing effective XNAV measurements to maintain spacecraft tracking, independent of other navigation sources. In the particular use case of the SEXTANT mission, it also provides a novel scheduling algorithm which addresses the need to prioritize and manage pulsar observations for effective navigation.

PULSAR↗

Application of X-Ray Pulsar Navigation: A Characterization of the Earth Orbit Trade Space

The potential for pulsars as a navigation source has been studied since their discovery in 1967. X-ray pulsar navigation (XNAV) is a celestial navigation system that uses the consistent timing nature of x-ray photons from millisecond pulsars (MSP) to perform space navigation. By comparing the detected arrival of x-ray photons to a reference database of expected pulsar light-curve timing models, one can infer a range and range rate measurement based on light time delay. Much of the challenge of XNAV comes from the faint signal, availability, and distant nature of pulsars. This is a study of potential pulsar XNAV measurements to measure extended Kalman filter (EKF) tracking performance with a wide trade space of bounded Earth orbits, using a simulation of existing x-ray detector space hardware. An example of an x-ray detector for XNAV is the NASA Station Explorer for X-ray Timing and Navigation (SEXTANT) mission, a technology demonstration of XNAV set to perform on the International Space Station (ISS) in late 2016early 2017. XNAV hardware implementation is driven by trajectory and environmental influences which add noise to the x-ray pulse signal. In a closed Earth orbit, the radiation environment can exponentially increase the signal noise from x-ray pulsar sources, decreasing the quality and frequency of measurements. The SEXTANT mission in particular improves on the signal to noise ratio by focusing an array of 56 x-ray silicon drift detectors at one pulsar target at a time. This reduces timing glitches and other timing noise contributions from ambient x-ray sources to within a 100 nanosecond resolution. This study also considers the SEXTANT scheduling challenges inherent in a single target observation. Finally, as the navigation sources are now relatively inertial targets, XNAV measurements are also subject to periods of occultation from various celestial bodies. This study focuses on the characterization of these drivers in closed Earth orbits and is not a tuning analysis of the EKF. The study shows that the closed Earth orbit for XNAV performance is reliant on the orbit semi-major axis and eccentricity as well as orbit inclination. These parameters are the primary drivers of pulsar measurement availability and significantly influence the natural spacecraft orbit dynamics. Sensitivity to initial orbit determination error growth due to the scarcity of XNAV measurements within an orbital period require appropriate timing of initial XNAV measurements. The orbit angles of argument of perigee and right ascension of the ascending node, alongside the other orbit parameters, complete the initial cadence of XNAV measurements. The performance of initial XNAV measurements then propagates throughout the experimental period. The study provides a basis to missions who wish to consider XNAV as a potential navigation source in a closed Earth orbit design. It provides a family of orbit trajectories as well as other modeling considerations needed to effectively evaluate if XNAV is an effective navigation source for a potential mission. As an EKF is sensitive to a linearized estimated state, this study has a direct benefit of providing effective XNAV measurements to maintain spacecraft tracking, independent of other navigation sources. In the particular use case of the SEXTANT mission, it also provides a novel scheduling algorithm which addresses the need to prioritize and manage pulsar observations for effective navigation.

x-ray↗

Applying the principles of isotope analysis in plant and animal ecology to forensic science in the Americas

Forensic science is an applied science that helps investigators answer questions at the heart of a crime: Who, What, When, Where, and How. As part of that process, forensic practitioners use basic laws of chemistry and physics to examine evidence and theories of origin to interpret said evidence. Ecological research has had an extraordinary impact on forensic science through the development and application of a specialized measurement technique useful for examining evidence that appears chemically and physically identical: stable isotope analysis. In this review, we examine the proven utility of stable isotope analysis in forensic settings from an ecological perspective, concentrating on work from the Americas completed within the last three decades. Our primary focus is on the interpretation of stable isotope ratios using plant and animal physiological models for “source inference,” although we also recognize that sample-to-sample comparison of stable isotope ratios can be useful in many forensic settings. Examples of the forensic application of stable isotope analysis cover the gamut, from cotton used in counterfeit bills to anthrax shipped through the U.S. Postal Service, from beer adulterated with cheap adjuncts to human remains discovered in shallow graves and never identified. Recent methodological developments—including isotope analysis of specific compounds, multiple elements, and radioisotopes—and the generation of isotope landscapes, or isoscapes, approaches for data interpretation promises that stable isotope analysis will continue to be used in forensic settings for another three decades… and beyond.

Chesson, Lesley A.↗

MADLens, a python package for fast and differentiable non-Gaussian lensing simulations

Here, we present MADLens a python package for producing non-Gaussian lensing convergence maps at arbitrary source redshifts with unprecedented precision. MADLens is designed to achieve high accuracy while keeping computational costs as low as possible. A MADLens simulation with only particles produces convergence maps whose power agrees with theoretical lensing power spectra up to within the accuracy limits of HaloFit. This is made possible by a combination of a highly parallelizable particle-mesh algorithm, a sub-evolution scheme in the lensing projection, and a machine-learning inspired sharpening step. Further, MADLens is fully differentiable with respect to the initial conditions of the underlying particle-mesh simulations and a number of cosmological parameters. These properties allow MADLens to be used as a forward model in Bayesian inference algorithms that require optimization or derivative-aided sampling. Another use case for MADLens is the production of large, high resolution simulation sets as they are required for training novel deep-learning-based lensing analysis tools. We make the MADLens package publicly available under a Creative Commons License

79 ASTRONOMY AND ASTROPHYSICS↗

Multichannel Analysis of Surface Waves Accelerated (MASWAccelerated): Software for efficient surface wave inversion using MPI and GPUs

Multichannel Analysis of Surface Waves (MASW) is a technique frequently used in geotechnical engineering and engineering geophysics to infer 1D layered models of seismic shear wave velocities in the top tens to hundreds of meters of the subsurface. We aim to accelerate MASW calculations by capitalizing on modern computer hardware available in the workstations of most engineers: multiple cores and graphics processing units (GPUs). We propose new parallel and GPU accelerated algorithms for computing 1D MASW inversion, and provide software implementations in C using Message Passing Interface (MPI) and CUDA. These algorithms take advantage of sparsity that arises in the problem, and the work balance between processes considers typical data trends. We compare our methods to an existing open source Matlab MASW tool. Our serial C implementation achieves a 2x speedup over the Matlab software, and we continue to see improvements by parallelizing the problem with MPI. Here we see nearly perfect strong and weak scaling for uniform data, and improve strong scaling for realistic data by repartitioning the problem to process mapping. By utilizing GPUs available on most modern workstations, we observe an additional 1.3x speedup over the serial C implementation on the first use of the method. We typically repeatedly evaluate theoretical dispersion curves as part of an optimization procedure, and on the GPU the kernel can be cached for faster reuse on later runs. We observe a 3.2x speedup on the cached GPU runs compared to the serial C runs. This work is the first open-source parallel or GPU-accelerated software tool for MASW imaging, and should enable geotechnical engineers to fully utilize all computer hardware at their disposal.

58 GEOSCIENCES↗

Explainable machine learning of the underlying physics of high-energy particle collisions

We present an implementation of an explainable and physics-aware machine learning model capable of inferring the underlying physics of high-energy particle collisions using the information encoded in the energy-momentum four-vectors of the final state particles. We demonstrate the proof-of-concept of our White Box AI approach using a Generative Adversarial Network (GAN) which learns from a DGLAP-based parton shower Monte Carlo event generator. The constrained generator network architecture mimics the structure of a parton shower exhibiting similarities with Recurrent Neural Networks (RNNs). We show, for the first time, that our approach leads to a network that is able to learn not only the final distribution of particles, but also the underlying parton branching mechanism, i.e. the Altarelli-Parisi splitting function, the ordering variable of the shower, and the scaling behavior. While the current work is focused on perturbative physics of the parton shower, we foresee a broad range of applications of our framework to areas that are currently difficult to address from first principles in QCD. Examples include nonperturbative and collective effects, factorization breaking and the modification of the parton shower in heavy-ion, and electron-nucleus collisions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Fragmentation of the edge of a terminated Cu nanolayer within a Nb matrix upon annealing

Here, we investigate the evolution of the tip of a terminated Cu nano-layer embedded within a Nb matrix upon annealing at 800°C. The tip of the terminated Cu layer fragments progressively into a series of single-crystal Cu particles upon continued annealing. This finding shows that layer fragmentation is the second step, after layer pinchoff, in the thermal coarsening of nanolayered metal composites. Comparison to prior phase-field modeling permits us to infer that the mobility at 800°C of Cu-Nb, as modeled by the Cahn-Hillard equation, is M(c)=(M 0 = 7.3 ± 3.6 nm 5 / s • eV) |1 – c 2 |, where the order parameter, c, is 1 for Cu and -1 for Nb.

36 MATERIALS SCIENCE↗

Towards fast and accurate predictions of radio frequency power deposition and current profile via data-driven modelling: applications to lower hybrid current drive

Three machine learning techniques (multilayer perceptron, random forest and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. A single GENRAY/CQL3D simulation without radial diffusion of fast electrons requires several minutes of wall-clock time to complete, which is acceptable for many purposes, but too slow for integrated modelling and real-time control applications. The machine learning models use a database of more than 16 000 GENRAY/CQL3D simulations for training, validation and testing. Latin hypercube sampling methods ensure that the database covers the range of nine input parameters ( $n_{e0}$ , $T_{e0}$ , $I_p$ , $B_t$ , $R_0$ , $n_{\|}$ , $Z_{{\rm eff}}$ , $V_{{\rm loop}}$ and $P_{{\rm LHCD}}$ ) with sufficient density in all regions of parameter space. The surrogate models reduce the inference time from minutes to $\sim$ ms with high accuracy across the input parameter space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bayesian Active Learning for Scanning Probe Microscopy: From Gaussian Processes to Hypothesis Learning

Recent progress in machine learning methods and the emerging availability of programmable interfaces for scanning probe microscopes (SPMs) have propelled automated and autonomous microscopies to the forefront of attention of the scientific community. However, enabling automated microscopy requires the development of task-specific machine learning methods, understanding the interplay between physics discovery and machine learning, and fully defined discovery workflows. This, in turn, requires balancing the physical intuition and prior knowledge of the domain scientist with rewards that define experimental goals and machine learning algorithms that can translate these to specific experimental protocols. Here, we discuss the basic principles of Bayesian active learning and illustrate its applications for SPM. We progress from the Gaussian process as a simple data-driven method and Bayesian inference for physical models as an extension of physics-based functional fits to more complex deep kernel learning methods, structured Gaussian processes, and hypothesis learning. These frameworks allow for the use of prior data, the discovery of specific functionalities as encoded in spectral data, and exploration of physical laws manifesting during the experiment. Here, the discussed framework can be universally applied to all techniques combining imaging and spectroscopy, SPM methods, nanoindentation, electron microscopy and spectroscopy, and chemical imaging methods and can be particularly impactful for destructive or irreversible measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Effects of Localized Thermal Pressure on Equilibrium Magnetic Fields and Particle Drifts in The Inner Magnetosphere

Localized injections of hot anisotropic plasma sheet particles into the inner magnetosphere can significantly deform the quiet time dipole–like magnetic field and thus disturb electron and ion's drift paths and scattering rates. Although many details of magnetic field deformation can be inferred from empirical models, roles of different characteristics of injected plasma on the structure of such deformation require further investigation. In this study, we use the 2–D axisymmetric equilibrium model to calculate self–consistent magnetic field in force balance with a Gaussian thermal pressure distribution characterized by four input parameters: the ratio between plasma pressure and magnetic pressure ( β ) at the pressure peak β 0 , the radial location of the pressure peak L 0 , the width of the half peak pressure σ 0 , and the equatorial pressure anisotropy A e . Using the modeled magnetic field, we find that the magnetic field perturbation increases with increasing β 0 and decreasing σ 0 while the magnetic curvature perturbation increases with increasing A e , β 0 , and σ 0 and decreasing L 0 . For energetic particles the change of magnetic gradient drift motion is much greater than that of curvature drift motion. The magnetic dip structure formation requires a critical β value that increases with increasing σ 0 and decreasing L 0 . Despite the unavailability of observations in the existing literatures to check the condition of magnetic dip formation, such condition will be checked against observations as a future study. Lastly, we also use 3–D ring current–atmosphere interactions model with self–consistent magnetic field model to illustrate the effect of azimuthal pressure distribution, which is relevant to asymmetric ring current.

79 ASTRONOMY AND ASTROPHYSICS↗

Observational Evidence of State‐Dependent Positive and Negative Land Surface Feedback on Afternoon Deep Convection Over the Southern Great Plains

Abstract The detection of soil moisture‐precipitation feedback is complicated by the dependence of feedback sign on atmospheric state and cloud types. We evaluated responses of convective clouds to changes in surface evaporative fraction (EF) from satellite and ground‐based measurements over the Southern Great Plains in summer months. Strong responses to EF were inferred only after accounting for tropospheric states supporting systematically positive or negative responses. State dependence (inferred from controlled model experiments) explains the weak response found in previous studies. The probability of transitions from shallow to deep cumulus, and subsequent precipitation onset, increases by 58% over drier surfaces in weakly stratified environments having deeper boundary layers. Boundary layer dynamics play a more important role than low‐level relative humidity in these environments. Wetter surfaces yield increased deep convection in moist and thermally stratified environments. The distinct pathways of convective initiation over dry and wet land surfaces can play important roles in drought onset and termination.

Qiu, Shaoyue↗

A Robust Estimate of Continental-Scale Terrestrial Carbon Sinks Using GOSAT XCO 2 Retrievals

Satellite XCO 2 retrievals could improve the estimates of surface carbon fluxes, but it remains unknown on what scales these estimates are robust. Here, we use the time-dependent Bayesian synthesis top-down method and prior net ecosystem exchanges (NEEs) from 12 terrestrial biosphere models (TBMs) to infer the monthly carbon fluxes of 51 land regions with constraints by GOSAT XCO 2 retrievals. We find that the uncertainty (standard deviation of 12 TBMs) reduction rates (uncertainty reduction rate (URR)) decrease significantly at decreasing spatial scales. On the continental-scale, the mean URR is about 57%, and the annual and seasonal cycle estimates of NEE are rather robust. The evaluation shows that the posterior CO 2 concentrations are significantly improved at the continental scale. Our study suggests that the GOSAT XCO 2 can only promise a robust continental-scale NEE estimate and improving the XCO 2 accuracy is an effective way to achieve robust estimates on smaller scales under current spatial coverage.

58 GEOSCIENCES↗

Bioinformatic and experimental evidence for suicidal and catalytic plant THI4s

Like fungi and some prokaryotes, plants use a thiazole synthase (THI4) to make the thiazole precursor of thiamin. Fungal THI4s are suicide enzymes that destroy an essential active-site Cys residue to obtain the sulfur atom needed for thiazole formation. In contrast, certain prokaryotic THI4s have no active-site Cys, use sulfide as sulfur donor, and are truly catalytic. The presence of a conserved active-site Cys in plant THI4s and other indirect evidence implies that they are suicidal. To confirm this, we complemented the Arabidopsistz-1 mutant, which lacks THI4 activity, with a His-tagged Arabidopsis THI4 construct. LC–MS analysis of tryptic peptides of the THI4 extracted from leaves showed that the active-site Cys was predominantly in desulfurated form, consistent with THI4 having a suicide mechanism in planta. Unexpectedly, transcriptome data mining and deep proteome profiling showed that barley, wheat, and oat have both a widely expressed canonical THI4 with an active-site Cys, and a THI4-like paralog (non-Cys THI4) that has no active-site Cys and is the major type of THI4 in developing grains. Transcriptomic evidence also indicated that barley, wheat, and oat grains synthesize thiamin de novo, implying that their non-Cys THI4s synthesize thiazole. Structure modeling supported this inference, as did demonstration that non-Cys THI4s have significant capacity to complement thiazole auxotrophy in Escherichia coli. There is thus a prima facie case that non-Cys cereal THI4s, like their prokaryotic counterparts, are catalytic thiazole synthases. Bioenergetic calculations show that, relative to suicide THI4s, such enzymes could save substantial energy during the grain-filling period.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Can local white dwarfs constrain dark matter interactions?

We investigate whether nearby white dwarfs (WDs) can constrain dark matter (DM) interactions with ordinary matter. As experimental sensitivity improves, driven by the Gaia mission, the sample volume of nearby WDs has been increasing over recent years. We carefully select a sample of ten cold, isolated, non-magnetic WDs within 13 pc of the Sun. We model their carbon-oxygen interior using a finite temperature relativistic equation of state and model atmospheres to infer their core temperatures. This enables us to perform a thorough estimation of the DM capture rate and evaporation mass using actual astrophysical observations. Given the low local DM density, we focus on DM that annihilates into long-lived mediators, which escape the WD and later decay into photons. While Fermi-LAT data shows no significant gamma-ray excess, future telescopes, CTA North & South, LHAASO, SWGO, could probe DM-nucleon cross sections down to ∼ 10 -41 cm 2 for DM masses above the TeV scale. Our results are competitive with current direct detection bounds (e.g., LZ) in the multi-TeV regime. This work underscores the importance of systematic WD studies in the broader landscape of DM detection and demonstrates the synergy between astrophysical and terrestrial searches in exploring DM interactions.

79 ASTRONOMY AND ASTROPHYSICS↗