Search NASA⌕ Search

SEARCH · Search NASA

Results for “ASTRONOMICAL MODEL”

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

Beyond Point Estimates: Benchmarking Uncertainty Quantification Methods on the AION-1 Astronomical Foundation Model

Foundation models for astronomical surveys offer powerful learned representations that can be transferred to downstream regression tasks such as galaxy property estimation. However, point predictions alone are insufficient for scientific inference; reliable uncertainty quantification (UQ) is essential. We compare seven UQ methods on galaxy property regression using frozen AION-1 foundation-model embeddings, predicting redshift, stellar mass, stellar-population age, gas-phase metallicity, and specific star-formation rate, from Legacy Survey photometry/imaging and DESI spectra, with PROVABGS-derived labels. Distribution-free conformal methods achieve marginal coverage within $\sim$1 pp of the nominal 90% across all properties, while non-conformal baselines (Deep Ensembles, MC~Dropout) fail to calibrate reliably. Among conformal approaches, Conformalized Quantile Regression (CQR) delivers the best coverage in the bin with the poorest model predictions. More importantly, only the Locally Valid and Discriminative (LVD) framework -- particularly when operating on AION-1 embeddings -- also provides finite-sample \emph{local validity}, producing intervals that adapt to each galaxy's local prediction difficulty rather than relying on marginal guarantees alone. These results establish conformal prediction, and LVD in particular, as the preferred UQ framework for uncertainty-aware inference on foundation-model embeddings in astrophysics.

Tame-Narvaez, Karla [Fermilab] (ORCID:000000022249↗

Towards an astronomical foundation model for stars with a transformer-based model

ABSTRACT Rapid strides are currently being made in the field of artificial intelligence using transformer-based models like Large Language Models (LLMs). The potential of these methods for creating a single, large, versatile model in astronomy has not yet been explored. In this work, we propose a framework for data-driven astronomy that uses the same core techniques and architecture as used by LLMs. Using a variety of observations and labels of stars as an example, we build a transformer-based model and train it in a self-supervised manner with cross-survey data sets to perform a variety of inference tasks. In particular, we demonstrate that a single model can perform both discriminative and generative tasks even if the model was not trained or fine-tuned to do any specific task. For example, on the discriminative task of deriving stellar parameters from Gaia XP spectra, we achieve an accuracy of 47 K in Teff, 0.11 dex in log g, and 0.07 dex in [M/H], outperforming an expert XGBoost model in the same setting. But the same model can also generate XP spectra from stellar parameters, inpaint unobserved spectral regions, extract empirical stellar loci, and even determine the interstellar extinction curve. Our framework demonstrates that building and training a single foundation model without fine-tuning using data and parameters from multiple surveys to predict unmeasured observations and parameters is well within reach. Such ‘Large Astronomy Models’ trained on large quantities of observational data will play a large role in the analysis of current and future large surveys.

Leung, Henry W. (ORCID:0000000200362752)↗

Parameters of cosmological models and recent astronomical observations

For different gravitational models we consider limitations on their parameters coming from recent observational data for type Ia supernovae, baryon acoustic oscillations, and from 34 data points for the Hubble parameter H(z) depending on redshift. We calculate parameters of 3 models describing accelerated expansion of the universe: the ΛCDM model, the model with generalized Chaplygin gas (GCG) and the multidimensional model of I. Pahwa, D. Choudhury and T.R. Seshadri. In particular, for the ΛCDM model 1σ estimates of parameters are: H{sub 0}=70.262±0.319 km {sup -1}Mp {sup -1}, Ω{sub m}=0.276{sub -0.008}{sup +0.009}, Ω{sub Λ}=0.769±0.029, Ω{sub k}=-0.045±0.032. The GCG model under restriction 0α≥ is reduced to the ΛCDM model. Predictions of the multidimensional model essentially depend on 3 data points for H(z) with z≥2.3.

79 ASTRONOMY AND ASTROPHYSICS↗

The FluxCompensator: Making Radiative Transfer Models of Hydrodynamical Simulations Directly Comparable to Real Observations

When modeling astronomical objects throughout the universe, it is important to correctly treat the limitations of the data, for instance finite resolution and sensitivity. In order to simulate these effects, and to make radiative transfer models directly comparable to real observations, we have developed an open-source Python package called the FluxCompensator that enables the post-processing of the output of 3D Monte Carlo radiative transfer codes, such as Hyperion. With the FluxCompensator, realistic synthetic observations can be generated by modeling the effects of convolution with arbitrary point-spread functions, transmission curves, finite pixel resolution, noise, and reddening. Pipelines can be applied to compute synthetic observations that simulate observatories, such as the Spitzer Space Telescope or the Herschel Space Observatory . Additionally, this tool can read in existing observations (e.g., FITS format) and use the same settings for the synthetic observations. In this paper, we describe the package as well as present examples of such synthetic observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling the Altitude Distribution of Meteor Head Echoes Observed with HPLA Radars: Implications for the Radar Detectability of Meteoroid Populations

The altitude distribution of meteors detected by a radar is sensitive to the instrument’s response function and can thus provide insight into the physical processes involved in radar measurements. This, in turn, can be used to determine the rate of ablation and ionization of the meteoroids and ultimately the input flux on Earth. In this work, we model the radar meteor head echo altitude distribution for three High Power and Large Aperture radar systems, by considering meteoroid populations from the main cometary family sources. In this simulation, we first use the results of a dynamical model of small meteoroids impacting Earth’s upper atmosphere to model the incoming mass, velocity, and entry angular distributions. We then combine these with the Chemical Ablation Model and establish the meteoroid ionization rates as a function of mass, velocity, and entry angle in order to determine the altitude at which these radars should detect the produced meteors and the portion of produced meteors from each population that are detected by these radars. We explore different sizes of head plasma as well as the possible effects on radar scattering of the head echo aspect sensitivity. We find that the modeled altitude distributions are generally in good agreement with measurements, particularly for ultra-high-frequency radars. In addition, our results indicate that the number of particles from Jupiter Family Comets (JFCs) required to fit the observations is lower than predicted by astronomical models. It is not clear yet if this discrepancy is due to the overprediction of JFC meteoroids by dynamical models or due to unaccounted physical processes in the treatment of ablation, ionization, and detections of meteoroids as they pass through Earth’s atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Poisson_CCD: A dedicated simulator for modeling CCDs

A dedicated simulator, Poisson_CCD, has been constructed which models astronomical CCDs by solving Poisson’s equation numerically and simulating charge transport within the CCD. The potentials and free carrier densities within the CCD are self-consistently solved for, giving realistic results for the charge distribution within the CCD storage wells. The simulator has been used to model the CCDs that are being used to construct the LSST digital camera. The simulator output has been validated by comparing its predictions with several different types of CCD measurements, including astrometric shifts, brighter-fatter induced pixel–pixel covariances, saturation effects, and diffusion spreading. The code is open source and freely available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Toward accurate thermochemistry of the {sup 24}MgH, {sup 25}MgH, and {sup 26}MgH molecules at elevated temperatures: Corrections due to unbound states

The total partition functions Q(T) and their first two moments Q{sup ′}(T) and Q{sup ″}(T), together with the isobaric heat capacities C{sub p}(T), are computed a priori for three major MgH isotopologues on the temperature range of T = 100–3000 K using the recent highly accurate potential energy curve, spin-rotation, and non-adiabatic correction functions of Henderson et al. [J. Phys. Chem. A 117, 13373 (2013)]. Nuclear motion computations are carried out on the ground electronic state to determine the (ro)vibrational energy levels and the scattering phase shifts. The effect of resonance states is found to be significant above about 1000 K and it increases with temperature. Even very short-lived states, due to their relatively large number, have significant contributions to Q(T) at elevated temperatures. The contribution of scattering states is around one fourth of that of resonance states but opposite in sign. Uncertainty estimates are given for the possible error sources, suggesting that all computed thermochemical properties have an accuracy better than 0.005% up to 1200 K. Between 1200 and 2500 K, the uncertainties can rise to around 0.1%, while between 2500 K and 3000 K, a further increase to 0.5% might be observed for Q{sup ″}(T) and C{sub p}(T), principally due to the neglect of excited electronic states. The accurate thermochemical data determined are presented in the supplementary material for the three isotopologues of {sup 24}MgH, {sup 25}MgH, and {sup 26}MgH at 1 K increments. These data, which differ significantly from older standard data, should prove useful for astronomical models incorporating thermodynamic properties of these species.

74 ATOMIC AND MOLECULAR PHYSICS↗

Fully anharmonic infrared cascade spectra of polycyclic aromatic hydrocarbons

The infrared (IR) emission of polycyclic aromatic hydrocarbons (PAHs) permeates our universe; astronomers have detected the IR signatures of PAHs around many interstellar objects. The IR emission of interstellar PAHs differs from their emission as seen under conditions on Earth as they emit through a collisionless cascade down through their excited vibrational states from high internal energies. The difficulty in reproducing interstellar conditions in the laboratory results in a reliance on theoretical techniques. However, the size and complexity of PAHs require careful consideration when producing the theoretical spectra. In this work, we outline the theoretical methods necessary to lead to fully theoretical IR cascade spectra of PAHs including: an anharmonic second order vibrational perturbation theory treatment, the inclusion of Fermi resonances through polyads, and the calculation of anharmonic temperature band shifts and broadenings (including resonances) through a Wang-Landau approach. We also suggest a simplified scheme to calculate vibrational emission spectra that retain the essential characteristics of the full IR cascade treatment and can directly transform low temperature absorption spectra in IR cascade spectra. Finally, we show that past astronomical models were in error in assuming a 15 cm -1 correction was needed to account for anharmonic emission effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An Investigation of LED Street Lighting's Impact on Sky Glow

A significant amount of public attention has recently focused on perceived impacts of converting street lighting from incumbent lamp-based products to LED technology. Much of this attention pertains to the higher content of short wavelength light (commonly referred to as "blue light") of LEDs and its attendant influences on sky glow (a brightening of the night sky that can interfere with astronomical observation and may be associated with a host of other issues). The complexity of this topic leads to common misunderstandings and misperceptions among the public, and for this reason the U.S. Department of Energy Solid-State Lighting Program embarked on a study of sky glow using a well-established astronomical model to investigate some of the primary factors influencing sky glow. This report details the results of the investigation and attempts to present those results in terms accessible to the general lighting community. The report also strives to put the results into a larger context, and help educate interested readers on various topics relevant to the issues being discussed.

99 GENERAL AND MISCELLANEOUS↗

Blind Denoising by Self-Supervised Neural Networks in Astronomical Datasets (Noise2Self4Astro)

the same part of the sky, with different shot, read-out, cosmic ray, and detector noise. Traditionally, in order to gain denoised model for astronomical objects, researchers take the average of multiple images of the same object or apply pre-processing filters such as Gaussian filter or Median filter, feeding them into “Galfit”, which is a software for astronomical object model fitting. Instead of that, we may train a network to predict one noisy image from the other based on the Noise2Self algorithm. The resulting output will be a semi-denoised image, which would later be put into galfit to complete the denoising process.

79 ASTRONOMY AND ASTROPHYSICS↗

The Origin of Chemical Elements: Connecting Laboratory Nuclear Astrophysics with Astronomical Observations through Nucleosynthesis Modeling

Nuclear astrophysics is a field at the frontier of nuclear science which continues to enjoy significant progress and exciting prospects for the future. It addresses the fundamental question of where and how the chemical elements are made. Except for hydrogen and helium, all other chemical elements are made either in stars or in cosmic explosions such as supernovae or neutron star mergers. To understand the formation of elements we need to understand the relevant nuclear physics processes that govern the interaction between nuclei. We need to understand the conditions in the star or the explosion under which these nuclear processes take place. And we need to understand how individual events of element formation contribute to the overall chemical enrichment of our Galaxy over time. Experimental facilities allow for terrestrial experiments on the nuclear physics governing the element synthesis. Earth-based and space-based telescopes probe the Galaxy around us and provide information about the conditions in stars and supernovae, and about the chemical element enrichment of these different environments. Numerical simulations of stars, supernovae, and other astrophysical phenomena allow us to probe the conditions of extreme densities and temperatures required for element synthesis. Nucleosynthesis research is ideally posed to integrate advances from all three of these areas and advance our understanding of how and where elements are made. Nucleosynthesis research provides theory support for experiments and observations. And it makes predictions that can be tested in nuclear physics experiments and astronomical observations. The goals of this award were in three main areas: 1. Develop and implement an effective core-collapse supernova model that is computationally efficient and that captures the relevant physics for the scientific questions of interest 2. Advance our understanding of the synthesis of the lightest heavy elements 3. Increase the realism of supernova nucleosynthesis and determine its dependence on the properties of the exploding star As the projects progressed, an additional focus was added to this list: 1. Study the explosion properties and the signals from the most massive stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Joint Image Framework

The software implements probabilistic modeling of astronomical images of stars and galaxies as seen by multiple telescopes and/or at multiple observation times. The image models are rendered with the external software package GalSim (https://github.com GalSim- developers GalSim ). JIF primarily handles the parsing of image model parameters, Monte Carlo sampling, and diagnostic plots.

Dawson, William↗

A COMPREHENSIVE DUST MODEL APPLIED TO THE RESOLVED BETA PICTORIS DEBRIS DISK FROM OPTICAL TO RADIO WAVELENGTHS

We investigate whether varying the dust composition (described by the optical constants) can solve a persistent problem in debris disk modeling—the inability to fit the thermal emission without overpredicting the scattered light. We model five images of the β Pictoris disk: two in scattered light from the Hubble Space Telescope ( HST )/Space Telescope Imaging Spectrograph at 0.58 μ m and HST /Wide Field Camera 3 (WFC 3) at 1.16 μ m, and three in thermal emission from Spitzer /Multiband Imaging Photometer for Spitzer (MIPS) at 24 μ m, Herschel /PACS at 70 μ m, and Atacama Large Millimeter/submillimeter Array at 870 μ m. The WFC3 and MIPS data are published here for the first time. We focus our modeling on the outer part of this disk, consisting of a parent body ring and a halo of small grains. First, we confirm that a model using astronomical silicates cannot simultaneously fit the thermal and scattered light data. Next, we use a simple generic function for the optical constants to show that varying the dust composition can improve the fit substantially. Finally, we model the dust as a mixture of the most plausible debris constituents: astronomical silicates, water ice, organic refractory material, and vacuum. We achieve a good fit to all data sets with grains composed predominantly of silicates and organics, while ice and vacuum are, at most, present in small amounts. This composition is similar to one derived from previous work on the HR 4796A disk. Our model also fits the thermal spectral energy distribution, scattered light colors, and high-resolution mid-IR data from T-ReCS for this disk. Additionally, we show that sub-blowout grains are a necessary component of the halo.

79 ASTRONOMY AND ASTROPHYSICS↗

SIMULATION OF ASTRONOMICAL IMAGES FROM OPTICAL SURVEY TELESCOPES USING A COMPREHENSIVE PHOTON MONTE CARLO APPROACH

We present a comprehensive methodology for the simulation of astronomical images from optical survey telescopes. We use a photon Monte Carlo approach to construct images by sampling photons from models of astronomical source populations, and then simulating those photons through the system as they interact with the atmosphere, telescope, and camera. We demonstrate that all physical effects for optical light that determine the shapes, locations, and brightnesses of individual stars and galaxies can be accurately represented in this formalism. By using large scale grid computing, modern processors, and an efficient implementation that can produce 400,000 photons s{sup −1}, we demonstrate that even very large optical surveys can be now be simulated. We demonstrate that we are able to (1) construct kilometer scale phase screens necessary for wide-field telescopes, (2) reproduce atmospheric point-spread function moments using a fast novel hybrid geometric/Fourier technique for non-diffraction limited telescopes, (3) accurately reproduce the expected spot diagrams for complex aspheric optical designs, and (4) recover system effective area predicted from analytic photometry integrals. This new code, the Photon Simulator (PhoSim), is publicly available. We have implemented the Large Synoptic Survey Telescope design, and it can be extended to other telescopes. We expect that because of the comprehensive physics implemented in PhoSim, it will be used by the community to plan future observations, interpret detailed existing observations, and quantify systematics related to various astronomical measurements. Future development and validation by comparisons with real data will continue to improve the fidelity and usability of the code.

79 ASTRONOMY AND ASTROPHYSICS↗

Red Noise–based False Alarm Thresholds for Astrophysical Periodograms via Whittle’s Approximation to the Likelihood

Astronomers who search for periodic signals using Lomb–Scargle periodograms rely on false alarm level (FAL) estimates to identify statistically significant peaks. Although FALs are often calculated from white noise models, many astronomical time series suffer from red noise. Prewhitening is a statistical technique in which a continuum model is subtracted from the log power spectrum estimate, after which the observer can proceed with a white-noise treatment. Here we present a prewhitening-based method of calculating frequency-dependent FALs. We fit power laws and autoregressive models of order 1 to each Lomb–Scargle periodogram by minimizing the Whittle approximation to the negative log-likelihood (NLL), then calculate FALs based on the best-fit model power spectrum. Our technique is a novel extension of the Whittle NLL to datasets with uneven time sampling. We demonstrate FAL calculations using observations of α Cen B, GJ 581, HD 192310, synthetic data from the radial velocity (RV) fitting challenge, and Kepler observations of a differential rotator. The Kepler data analysis shows that only true rotation signals are detected by red noise FALs, while white noise FALs suggest all spurious peaks in the low-frequency range are significant. A high-frequency sinusoid injected into α Cen B logR$'$ HK observations exceeds the 1% red noise FAL despite having only 8.9% of the power of the dominant rotation signal. In a periodogram of HD 192310 RVs, peaks associated with differential rotation and planets are detected against the 5% red noise FAL without iterative model fitting or subtraction. The software for calculating red noise–based FALs is available on GitHub.

Astrostatistics (1882)↗

Mapping Stellar Surfaces. II. An Interpretable Gaussian Process Model for Light Curves

The use of Gaussian processes (GPs) as models for astronomical time series data sets has recently become almost ubiquitous, given their ease of use and flexibility. In particular, GPs excel at marginalization over the stellar signal when the variability due to starspots is treated as a nuisance, as in exoplanet transit modeling. However, these effective models are less useful in cases where the starspot signal is of primary interest, since it is not obvious how the parameters of the GP relate to physical parameters like the spot size, contrast, and latitudinal distribution. Instead, it is common practice to explicitly model the effect of individual starspots on the light curve and attempt to infer their properties via optimization or posterior inference. Unfortunately, this process is ill-posed and often computationally intractable when applied to stars with more than a few spots and/or to ensembles of many stars. Here we derive a closed-form expression for a GP that describes the light curve of a rotating, evolving stellar surface conditioned on a given distribution of starspot sizes, contrasts, and latitudes. We demonstrate that this model is correctly calibrated, allowing one to robustly infer physical parameters of interest from one or more light curves, including the typical spot radii and latitudes. Our GP has far-ranging implications for understanding the variability and magnetic activity of stars from light curves and radial velocity measurements, as well as for modeling correlated noise in exoplanet searches. Our implementation is efficient, user-friendly, and open-source, available in the package starry-process.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗