Search NASA⌕ Search

SEARCH · Search NASA

Results for “Astronomy data modeling”

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

Monte Carlo Simulations of Polarized Radiative Transfer in Neutron Star Atmospheres

Soft X-ray emission from neutron stars affords powerful diagnostic tools for uncovering their surface and interior properties, as well as their geometric configurations. In the atmospheres of neutron stars, the presence of magnetic fields alters the photon-electron scattering cross sections, resulting in nontrivial angular dependence of intensity and polarization of the emergent signals. This paper presents recent developments of our Monte Carlo simulation, MAGTHOMSCATT, which tracks the complex electric field vector for each photon during its transport. Our analysis encompasses the anisotropy and polarization characteristics of X-ray emission for field strengths ranging from nonmagnetic to extremely magnetized regimes that are germane to magnetars. In the very low field domain, we reproduced the numerical solution to the radiative transfer equation for nonmagnetic Thomson scattering, and provided analytical fits for the angular dependence of the intensity and the polarization degree. These fits can be useful for studies of millisecond pulsars and magnetic white dwarfs. By implementing a refined injection protocol, we show that, in the magnetar regime, the simulated intensity and polarization pulse profiles of emission from extended surface regions becomes invariant with respect to the ratio of photon (ω) and electron cyclotron (ω B ) frequencies once ω/ω B ≲ 0.01. This circumvents the need for simulations pertinent to really high magnetic field strengths, which are inherently slower. Our approach will be employed elsewhere to model observational data to constrain neutron star geometric parameters and properties of emitting hot spots on their surfaces.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Jipole: A Differentiable ipole-based Code for Radiative Transfer in Curved Spacetimes

Recent imaging of supermassive black holes by the Event Horizon Telescope has relied on exhaustive parameter-space searches, matching observations to large, precomputed libraries of theoretical models. As observational data become increasingly precise, the limitations of this computationally expensive approach grow more acute, creating a pressing need for more efficient methods. In this work, we present Jipole, an automatically differentiable (AD), ipole-based code for radiative transfer in curved spacetimes, designed to compute image gradients with respect to underlying model parameters. These gradients quantify how parameter changes—such as the black hole’s spin or the observer’s inclination—affect the image, enabling more efficient parameter estimation and reducing the number of required images. We validate Jipole against ipole in two analytical tests and then compare pixelwise intensity derivatives from AD with those from finite-difference methods. We then demonstrate the utility of these gradients by performing parameter recovery for an analytical model in three increasingly complex cases for the injected image: ideal, blurred, and blurred with added noise. In most cases, high-accuracy fits are obtained in only a few optimization steps, failing only in cases with extremely low signal-to-noise ratios. These results highlight the potential of AD-based methods to accelerate robust, high-fidelity model-data comparisons in current and future black hole imaging efforts.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2011

Contribution of X/Ka VLBI to Multi-Wavelength Celestial Frame Studies

This paper is an update of Sotuela et al. (2011) which improves their simulated Gaia frame tie precision by approximately 10% by adding three additional VLBI observing sessions. Astrometry at X/Ka-band (8.4/32 GHz) using NASAs Deep Space Network has detected 466 quasars with accuracies of 200-300 micro-arc seconds. A program is underway to reduce errors by a factor of 2-3. From our sample, 245 sources have optical magnitudes V less than 20 and should also be detectable by Gaia. A covariance study using existing X/Ka data and simulated Gaia uncertainties for the 345 objects yields a frame tie precision of 10-15 micro-arc seconds (1 - sigma). The characterization of wavelength dependent systematic from extended source morphology and core shift should benefit greatly from adding X/Ka-band measurements to S/X-band (2.3/8.4 GHz) measurements thus helping to constrain astrophysical models of the wavelength dependence of positions.

galaxies↗
DOE OSTI2025

Cosmic Ray Mediated Thermal Fronts in the Warm-Hot Circumgalactic Medium

We investigate the 1D plane-parallel front connecting the warm ($10^4$ K) and hot ($10^6$ K) phases of the circumgalactic medium (CGM), focusing on the influence of cosmic rays (CRs) in shaping these transition layers. We find that cosmic rays dictate the thermal balance while other fluxes (thermal conduction, radiative cooling, and gas flow) adjust to compensate. We compute column density ratios for selected transition temperature ions and compare them with observational data. While most of our models fail to reproduce the observations, a few are successful, although we make no claims for their uniqueness. Some of the discrepancies may indicate challenges in capturing the profiles in cooler, photoionized regions, as has been suggested for by previous efforts to model thermal transition layers.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2023

Primordial or Secondary? Testing Models of Debris Disk Gas with ALMA

The origin and evolution of gas in debris disks are still not well understood. Secondary gas production from cometary material or a primordial origin have been proposed. So far, observations have mostly concentrated on CO, with only a few C observations available. We overview the C and CO content of debris disk gas and test the state of the-art models. We use new and archival Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO and C I emission, complemented by C II data from Herschel, for a sample of 14 debris disks. This expands the number of disks with ALMA measurements of both CO and C I by 10 disks. We present new detections of C I emission toward three disks: HD 21997, HD 121191, and HD 121617. We use a simple disk model to derive gas masses and column densities. We find that current state-of-the-art models of secondary gas production overpredict the C 0 content of debris disk gas. This does not rule out a secondary origin, but might indicate that the models require an additional C removal process. Alternatively, the gas might be produced in transient events rather than a steady-state collisional cascade. We also test a primordial gas origin by comparing our results to a simplified thermochemical model. This yields promising results, but more detailed work is required before a conclusion can be reached. Our work demonstrates that the combination of C and CO data is a powerful tool to advance our understanding of debris disk gas.

Radiative transfer↗
DOE OSTI2024

A Novel Method to Train Classification Models for Structure Detection in In Situ Spacecraft Data

We present a method for creating spacecraft-like data which can be used to train Machine Learning (ML) models to detect and classify structures in in situ spacecraft data. First, we use the Grad-Shafranov equation to numerically solve for several magnetohydrostatic equilibria which are variations on a known analytic equilibrium. These equilibria are then used as the initial conditions for Particle-In-Cell simulations in which the structures of interest are observed and labeled. We then take one-dimensional slices through the simulations to replicate what a spacecraft collecting data from the simulation would observe. This sliced data then can be used as training data for the initial training of ML models intended for use on spacecraft data. We demonstrate the method applied to the problem of detecting small-scale plasmoids in the magnetotail, which is important for understanding complex magnetotail reconnection dynamics. The simple 1D classifier we train is able to detect more than 70% of the plasmoid points in the data set but also produces a large number of false positives. Our further work on this example problem is detailed, and further potential uses of the method are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1988

Contamination of grazing incidence EUV mirrors - An assessment

Contamination assessment for space optical systems requires an understanding of the sensitivity of component performance, e.g. mirror reflectance, to materials deposited on the mirror surface. In a previous study, the sensitivity of typical normal incidence mirror coatings to surface deposits of generic hydrocarbons was reported. Recent activity in the development of grazing incidence telescopes for extreme ultraviolet space astronomy has stimulated the need for a similar assessment in the spectral region extending from approximately 100 A to 1000 A. The model used for analysis treats the contamination layer as a continuous thin film deposited on the mirror surface. The mirror surfaces selected for this study are opaque vacuum deposited gold and the uncoated and polished Zerodur. Scatter caused by film irregularities or particulates are not included in this assessment. Parametric evaluations at 100, 500, and 1000 A determine the sensitivity of mirror reflectance to a range of optical constants selected for the generic contaminants. This sensitivity analysis combined with the limited amount of optical data in the EUV for hydrocarbons, is used to select representative optical constants for the three wavelength regions. Reflectance versus contamination layer thickness curves are then calculated and used to determine critical thickness limits based on allowable reflectance change. Initial observations indicate that thickness limits will be highly dependent on the real part of the complex index of refraction of the contaminant film being less than 1.0. Preliminary laboratory measurements of samples contaminated with some commonly encountered hydrocarbons confirm trends indicated in the analytical studies.

Osantowski, John F.↗
DOE OSTI2024

Light new physics in B → K ( * ) ν ν ¯ ?

The study of the rare decays B → K ( * ) ν ν ¯ offers a window into the dynamics operating at the electroweak scale, allowing studies of the Standard Model and searches for heavy new physics. However, the analysis of these decays is also potentially sensitive to the on shell production of new light bosons X through the process B → K ( * ) X . In particular, Belle II has recently measured B + → K + ν ν ¯ , finding a 2.8 σ excess under the assumption of heavy new physics. Since this excess is rather localized in the kaon energy, a fit that includes the decay mode B + → K + X to the kinematic distributions prefers m X ≈ 2 GeV with branching fraction Br [ B → K X ] = ( 8.8 ± 2.5 ) × 10 − 6 and a significance of ≈ 3.6 σ . However, no excess was found in the measurements of B → K ( * ) ν ν ¯ , and a global analysis of the Belle II and data leads to Br [ B → K X ] = ( 5.1 ± 2.1 ) × 10 − 6 with a reduced significance of ≈ 2.4 σ . We then study various simplified dark-flavored models and present a possible UV completion based on a gauged B 3 − L 3 symmetry, highlighting the discovery potential of dedicated searches for B → K ( * ) X at Belle II. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
NASA NTRS1987

Starburst-driven superwinds from infrared galaxies

New data is presented that indicate that strong far infrared galaxies commonly have largescale emission line nebulae whose properties are suggestive of mass outflows (superwinds), presumably driven by the high supernova rate associated with the central starburst. These data include longslit spectra of M82 which show that the radial variation of the gas pressure in the emission line nebula is in excellent agreement with a previous wind model. The M82 nebula also has a LINER spectrum, consistent with shock heating. Morphologically and spectroscopically similar emission line nebulae were found in NGC253, and Arp 220 and NGC6240. A longslit spectroscopic investigation was conducted of 20 additional very powerful far-infrared galaxies and found that they generally have spatially extended emission line nebulae whose spectra closely resemble that of the M82 nebula. If the superwind interpretation is correct, it could have many important consequences in extragalactic astronomy.

Heckman, Timothy M.↗
NASA NTRS1999

Parameter Estimation in Astronomy with Poisson-Distributed Data: The (CHI)2(gamma) Statistic - 1

Applying the standard weighted mean formula, [Sigma (sub i)n(sub i)ssigma(sub i, sup -2)], to determine the weighted mean of data, n(sub i), drawn from a Poisson distribution, will, on average, underestimate the true mean by approx. 1 for all true mean values larger than approx.3 when the common assumption is made that the error of the i th observation is sigma(sub i) = max square root of n(sub i), 1).This small, but statistically significant offset, explains the long-known observation that chi-square minimization techniques which use the modified Neyman'chi(sub 2) statistic, chi(sup 2, sub N) equivalent Sigma(sub i)((n(sub i) - y(sub i)(exp 2)) / max(n(sub i), 1), to compare Poisson - distributed data with model values, y(sub i), will typically predict a total number of counts that underestimates the true total by about 1 count per bin. Based on my finding that weighted mean of data drawn from a Poisson distribution can be determined using the formula [Sigma(sub i)[n(sub i) + min(n(sub i), 1)](n(sub i) + 1)(exp -1)] / [Sigma(sub i)(n(sub i) + 1)(exp -1))], I propose that a new chi(sub 2) statistic, chi(sup 2, sub gamma) equivalent, should always be used to analyze Poisson- distributed data in preference to the modified Neyman's chi(exp 2) statistic. I demonstrated the power and usefulness of,chi(sub gamma, sup 2) minimization by using two statistical fitting techniques and five chi(exp 2) statistics to analyze simulated X-ray power - low 15 - channel spectra with large and small counts per bin. I show that chi(sub gamma, sup 2) minimization with the Levenberg - Marquardt or Powell's method can produce excellent results (mean slope errors approx. less than 3%) with spectra having as few as 25 total counts.

Mighell, Kenneth J.↗
DOE OSTI2024

Modeling the Enceladus dust plume based on in situ measurements performed with the Cassini Cosmic Dust Analyzer

We analyzed data recorded by the Cosmic Dust Analyzer on board the Cassini spacecraft during Enceladus dust plume traversals. Our focus was on profiles of relative abundances of grains of different compositional types derived from mass spectra recorded with the Dust Analyzer subsystem during the Cassini flybys E5 and E17. The E5 profile, corresponding to a steep and fast traversal of the plume, has already been analyzed. In this paper, we included a second profile from the E17 flyby involving a nearly horizontal traversal of the south polar terrain at a significantly lower velocity. Additionally, we incorporated dust detection rates from the High Rate Detector subsystem during flybys E7 and E21. We derived grain size ranges in the different observational data sets and used these data to constrain parameters for a new dust plume model. This model was constructed using a mathematical description of dust ejection implemented in the software package DUDI. Further constraints included published velocities of gas ejection, positions of gas and dust jets, and the mass production rate of the plume. Our model employs two different types of sources: diffuse sources of dust ejected with a lower velocity and jets with a faster and more colimated emission. From our model, we derived dust mass production rates for different compositional grain types, amounting to at least 28 kg s –1 . Previously, salt-rich dust was believed to dominate the plume mass based on E5 data alone. The E17 profile shows a dominance of organic-enriched grains over the south polar terrain, a region not well constrained by E5 data. By including both E5 and E17 profiles, we find the salt-rich dust contribution to be at most 1% by mass. This revision also results from an improved understanding of grain masses of various compositional types that implies smaller sizes for salt-rich grains. Our new model can predict grain numbers and masses for future mission detectors during plume traversals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Chimera F-Series Gravitational Wave Emission Sourced from Matter and Neutrino Anisotropy

This dataset contains gravitational wave data sourced from the the time-dependent fluid quadrupole motion as well as neutrino anisotropy, in the Chimera F-Series two-dimensional core collapse supernova simulations. Data from two models initiated from two different progenitors are presented: F15.78 and F15.79. Please see the README for more information about the data structure and progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

First detailed calculation of atmospheric neutrino foregrounds to the diffuse supernova neutrino background in Super-Kamiokande

The diffuse supernova neutrino background (DSNB)—a probe of the core-collapse mechanism and the cosmic star-formation history—has not been detected, but its discovery may be imminent. A significant obstacle for DSNB detection in Super-Kamiokande (Super-K) is detector backgrounds, especially due to atmospheric neutrinos (more precisely, these are foregrounds), which are not sufficiently understood. We perform the first detailed theoretical calculations of these foregrounds in the range 16–90 MeV in detected electron energy, taking into account several physical and detector effects, quantifying uncertainties, and comparing our predictions to the 15.9 live time years of pre-gadolinium data from Super-K stages I–IV. We show that our modeling reasonably reproduces this low-energy data as well as the usual high-energy atmospheric-neutrino data. To accelerate progress on detecting the DSNB, we outline key actions to be taken in future theoretical and experimental work. In a forthcoming paper, we use our modeling to detail how low-energy atmospheric-neutrino events register in Super-K and suggest new cuts to reduce their impact. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Strong Lensing by Galaxies

Strong gravitational lensing at the galaxy scale is a valuable tool for various applications in astrophysics and cosmology. Some of the primary uses of galaxy-scale lensing are to study elliptical galaxies’ mass structure and evolution, constrain the stellar initial mass function, and measure cosmological parameters. Since the discovery of the first galaxy-scale lens in the 1980s, this field has made significant advancements in data quality and modeling techniques. In this review, we describe the most common methods for modeling lensing observables, especially imaging data, as they are the most accessible and informative source of lensing observables. We then summarize the primary findings from the literature on the astrophysical and cosmological applications of galaxy-scale lenses. We also discuss the current limitations of the data and methodologies and provide an outlook on the expected improvements in both areas in the near future.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Advanced radiometric and interferometric milimeter-wave scene simulations

Smart munitions and weapons utilize various imaging sensors (including passive IR, active and passive millimeter-wave, and visible wavebands) to detect/identify targets at short standoff ranges and in varied terrain backgrounds. In order to design and evaluate these sensors under a variety of conditions, a high-fidelity scene simulation capability is necessary. Such a capability for passive millimeter-wave scene simulation exists at TRW. TRW's Advanced Radiometric Millimeter-Wave Scene Simulation (ARMSS) code is a rigorous, benchmarked, end-to-end passive millimeter-wave scene simulation code for interpreting millimeter-wave data, establishing scene signatures and evaluating sensor performance. In passive millimeter-wave imaging, resolution is limited due to wavelength and aperture size. Where high resolution is required, the utility of passive millimeter-wave imaging is confined to short ranges. Recent developments in interferometry have made possible high resolution applications on military platforms. Interferometry or synthetic aperture radiometry allows the creation of a high resolution image with a sparsely filled aperture. Borrowing from research work in radio astronomy, we have developed and tested at TRW scene reconstruction algorithms that allow the recovery of the scene from a relatively small number of spatial frequency components. In this paper, the TRW modeling capability is described and numerical results are presented.

Hauss, B. I.↗
DOE OSTI2025

A Probabilistic Model for Global EMIC Wave Activity Using Van Allen Probes Observations

Electromagnetic ion cyclotron (EMIC) waves play a key role in radiation belt dynamics through resonant interactions. However, their low occurrence probability, high variability, and spatial intermittency pose challenges for accurate modeling. In this study, we present a machine learning (ML)-based global EMIC wave model built on the entire data set from the Van Allen Probes mission. To capture the distinct statistical characteristics of wave occurrence and amplitude, the model is separated into two modules: an occurrence model trained using ML techniques, and a wave amplitude model sampled from observed probability distributions. The input parameters are limited to real-time or predictable variables to ensure practical applicability. Our model shows strong performance across the entire test set and demonstrates improved predictive capability over a baseline random occurrence model, particularly during quiet geomagnetic conditions. Evaluation during both quiet and active periods confirms the model's ability to represent the clustered and intermittent nature of EMIC wave activity. Furthermore, the model provides global estimates of wave power, enabling integration with radiation belt electron data and showing signatures consistent with wave-induced scattering. We found a good correlation between the global wave activity from the model and relativistic electron observation by Van Allen Probes, regardless of the availability of in situ wave observations. The modular structure of the model also allows for straightforward expansion for additional wave properties, such as wave frequency, which can be modeled independently. This flexible, event-sensitive approach offers a promising framework for data-driven radiation belt simulations and space weather applications.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Constraining nuclear mass models using 𝑟-process observables with multiobjective optimization

Modeling nuclear masses, particularly for nuclei far from stability, remains a key objective in nuclear physics. One contemporary approach is machine learning (ML), which trains on experimental data, but can suffer large errors when extrapolating toward neutron-rich species. In nature, such masses shape observables for the rapid neutron capture process (𝑟 process), which in principle could inform ML models. Here, we introduce a multiobjective optimization approach using the Pareto front algorithm. We show that this technique, capable of identifying models that generate 𝑟-process abundances aligning with both solar and stellar data, is a promising method to select ML models with reliable extrapolation power.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2012

A Study of Planetary Nebulae using the Faint Object Infrared Camera for the SOFIA Telescope

A planetary nebula is formed following an intermediate-mass (1-8 solar M) star's evolution off of the main sequence; it undergoes a phase of mass loss whereby the stellar envelope is ejected and the core is converted into a white dwarf. Planetary nebulae often display complex morphologies such as waists or torii, rings, collimated jet-like outflows, and bipolar symmetry, but exactly how these features form is unclear. To study how the distribution of dust in the interstellar medium affects their morphology, we utilize the Faint Object InfraRed CAmera for the SOFIA Telescope (FORCAST) to obtain well-resolved images of four planetary nebulae--NGC 7027, NGC 6543, M2-9, and the Frosty Leo Nebula--at wavelengths where they radiate most of their energy. We retrieve mid infrared images at wavelengths ranging from 6.3 to 37.1 micron for each of our targets. IDL (Interactive Data Language) is used to perform basic analysis. We select M2-9 to investigate further; analyzing cross sections of the southern lobe reveals a slight limb brightening effect. Modeling the dust distribution within the lobes reveals that the thickness of the lobe walls is higher than anticipated, or rather than surrounding a vacuum surrounds a low density region of tenuous dust. Further analysis of this and other planetary nebulae is needed before drawing more specific conclusions.

infrared sensors↗