Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stellar astronomy”

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

An Educational Guide for 2D Stellar Structure Calculations of Rapidly Rotating Stars using the ESTER code

The Evolution STEllaire en Rotation (ESTER) code is the first 2D stellar structure code to be made open-source and freely available to the astronomy and astrophysics community. An important and novel advancement of this code is that it can reproduce the distorted shape and observable signatures (e.g., gravity darkening) of rapidly rotating stars. ESTER also calculates the steady-state large-scale flows within the star, namely their differential rotation and associated meridional circulation. In this report, we explore and document the physics implemented within version 1.1.0rc2 of the ESTER code, in a way that complements published descriptions. We illustrate this physics by plotting how stellar structure parameters vary through stellar interiors at a range of latitudes and at different angular velocities. We investigate how the thin convective envelopes of intermediate mass stars vary with latitude when rapidly rotating, becoming deeper and thicker near the equator. Simple comparisons of ESTER model predictions (e.g., central temperature and density, luminosity) with the output from the Modules for Experiments in Stellar Astrophysics (MESA) code [Paxton et al., 2010] shows generally good agreement. Additional comparisons provide important benchmarking and verification for ESTER as a comparatively young code. Finally, we provide a guide for installing and running the code on our local university cluster, aimed at helping students to begin work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Quenching time-scales in the IllustrisTNG simulation

ABSTRACT The time-scales for galaxy quenching offer clues to its underlying physical drivers. We investigate central galaxy quenching time-scales in the IllustrisTNG 100-1 simulation, their evolution over time, and the pre-quenching properties of galaxies that predict their quenching time-scales. Defining quenching duration τq as the time between crossing specific star formation rate (sSFR) thresholds, we find that ${\sim} 40{{\ \rm per\ cent}}$ of galaxies quench rapidly with τq < 1 Gyr, but a substantial tail of galaxies can take up to 10 Gyr to quench. Furthermore, 29 per cent of galaxies that left the star-forming main sequence (SFMS) more than 2 Gyr ago never fully quench by z = 0. While the median τq is fairly constant with epoch, the rate of galaxies leaving the SFMS increases steadily over cosmic time, with the rate of slow quenchers being dominant around z ∼ 2–0.7. Compared to fast quenchers (τq < 1 Gyr), slow-quenching galaxies (τq > 1 Gyr) were more massive, had more massive black holes, had larger stellar radii, and accreted gas with higher specific angular momentum (AM) prior to quenching. These properties evolve little by z = 0, except for the accreting gas AM for fast quenchers, which reaches the same high AM as the gas in slow quenchers. By z = 0, slow quenchers also have residual star formation in extended gas rings. Using the expected relationship between stellar age gradient and τq for inside-out quenching we find agreement with Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) integral field unit (IFU) observations. Our results suggest the accreting gas AM and potential well depth determine the quenching time-scale.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

An idealized model for collapse of a molecular core within a filament

Here, we present an analytic treatment of the self-similar collapse of a uniform density ellipsoid to linear order in the deviations from sphericity. First we obtain a self-consistent closed-form solution for the evolution of an isolated ellipsoid and then impose the effects of an external forcing. This model describes the evolution under gravity of a pre-stellar core of molecular gas embedded in a much larger and lower-density filament. We compare with numerical solutions for the collapse without the limitation of small deviations. These show how the external perturbing force producing the initial deviation from sphericity is eventually surpassed by the anisotropic forces generated by the collapsing ellipsoid itself. This model should be useful in interpreting the way in which environment shapes the evolution of pre-stellar cores.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

w22_stars: Turbulent Convection in a 19M ⊙ Star [Slide]

We find large scale turbulent convection within the O-shell region in our massive star model. Within this region we also observe a redistribution of angular momentum that differs from 1D stellar evolution models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The Pulsating White Dwarf G117-B15A: Still the Most Stable Optical Clock Known

The pulsating hydrogen atmosphere white dwarf star G 117-B15A has been observed since 1974. Its main pulsation period at 215.19738823(63) s, observed in optical light curves, varies by only (5.12 ± 0.82) × 10 -15 s s -1 and shows no glitches, as pulsars do. The observed rate of period change corresponds to a change of the pulsation period by 1 s in 6.2 million yr. Here, we demonstrate that this exceptional optical clock can continue to put stringent limits on fundamental physics, such as constraints on interaction from hypothetical dark matter particles, as well as to search for the presence of external substellar companions.

1799↗
DOE OSTI2023

Comparison of the stellar populations of bulges and discs using the MaNGA survey

We use the MaNGA integral field spectroscopic survey of low-redshift galaxies to compare the stellar populations of the bulge and disc components, identified from their Sérsic profiles, for various samples of galaxies. Bulge-dominated regions tend to be more metal-rich and have slightly older stellar ages than their associated disc-dominated regions. The metallicity difference is consistent with the deeper gravitational potential in bulges relative to discs, which allows bulges to retain more of the metals produced by stars. The age difference is due to star formation persisting longer in discs relative to bulges. Relative to galaxies with lower stellar masses, galaxies with higher stellar masses tend to have bulge-dominated regions that are more metal-rich and older (in light-weighted measurements) than their disc-dominated regions. This suggests high-mass galaxies quench from the inside out, while lower-mass galaxies quench across the whole galaxy simultaneously. Early-type galaxies tend to have bulge-dominated regions the same age as their disc-dominated regions, while late-type galaxies tend to have disc-dominated regions significantly younger than their bulge-dominated regions. Central galaxies tend to have a greater metallicity difference between their bulge-dominated regions and disc-dominated regions than satellite galaxies at similar stellar mass. This difference may be explained by central galaxies being subject to mergers or extended gas accretion bringing new, lower-metallicity gas to the disc, thereby reducing the average metallicity and age of the stars; quenching of satellite discs may also play a role.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

via machinae : Searching for stellar streams using unsupervised machine learning

ABSTRACT We develop a new machine learning algorithm, via machinae, to identify cold stellar streams in data from the Gaia telescope. via machinae is based on ANODE, a general method that uses conditional density estimation and sideband interpolation to detect local overdensities in the data in a model agnostic way. By applying ANODE to the positions, proper motions, and photometry of stars observed by Gaia, via machinae obtains a collection of those stars deemed most likely to belong to a stellar stream. We further apply an automated line-finding method based on the Hough transform to search for line-like features in patches of the sky. In this paper, we describe the via machinae algorithm in detail and demonstrate our approach on the prominent stream GD-1. Though some parts of the algorithm are tuned to increase sensitivity to cold streams, the via machinae technique itself does not rely on astrophysical assumptions, such as the potential of the Milky Way or stellar isochrones. This flexibility suggests that it may have further applications in identifying other anomalous structures within the Gaia data set, for example debris flow and globular clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Simulation of Stark-broadened Hydrogen Balmer-line Shapes for DA White Dwarf Synthetic Spectra

Abstract White dwarfs (WDs) are useful across a wide range of astrophysical contexts. The appropriate interpretation of their spectra relies on the accuracy of WD atmosphere models. One essential ingredient of atmosphere models is the theory used for the broadening of spectral lines. To date, the models have relied on Vidal et al., known as the unified theory of line broadening (VCS). There have since been advancements in the theory; however, the calculations used in model atmosphere codes have only received minor updates. Meanwhile, advances in instrumentation and data have uncovered indications of inaccuracies: spectroscopic temperatures are roughly 10% higher and spectroscopic masses are roughly 0.1 M ⊙ higher than their photometric counterparts. The evidence suggests that VCS-based treatments of line profiles may be at least partly responsible. Gomez et al. developed a simulation-based line-profile code Xenomorph using an improved theoretical treatment that can be used to inform questions around the discrepancy. However, the code required revisions to sufficiently decrease noise for use in model spectra and to make it computationally tractable and physically realistic. In particular, we investigate three additional physical effects that are not captured in the VCS calculations: ion dynamics, higher-order multipole expansion, and an expanded basis set. We also implement a simulation-based approach to occupation probability. The present study limits the scope to the first three hydrogen Balmer transitions (H α , H β , and H γ ). We find that screening effects and occupation probability have the largest effects on the line shapes and will likely have important consequences in stellar synthetic spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The MASSIVE Survey. XVI. The Stellar Initial Mass Function in the Center of MASSIVE Early-type Galaxies

The stellar initial mass function (IMF) is a fundamental property in the measurement of stellar masses and galaxy star formation histories. In this work, we focus on the most massive galaxies in the nearby universe $\mathrm{log}({M}_{\star }/\,{M}_{\odot })\gt 11.2$. We obtain high-quality Magellan/LDSS-3 long-slit spectroscopy with a wide wavelength coverage of 0.4–1.01 μm for 41 early-type galaxies (ETGs) in the MASSIVE survey and derive high signal-to-noise spectra within an aperture of R e /8. Using detailed stellar synthesis models, we constrain the elemental abundances and stellar IMF of each galaxy through full spectral modeling. All the ETGs in our sample have an IMF that is steeper than a Milky Way (Kroupa) IMF. The best-fit IMF mismatch parameter, α IMF = (M/L)/(M/L) MW , ranges from 1.1 to 3.1, with an average of $\langle$α IMF $\rangle$ = 1.84, suggesting that on average, the IMF is more bottom heavy than Salpeter. Comparing the estimated stellar masses with the dynamical masses, we find that most galaxies have stellar masses that are smaller than their dynamical masses within the 1σ uncertainty. We complement our sample with lower-mass galaxies from the literature and confirm that $\mathrm{log}({\alpha }_{\mathrm{IMF}})$ is positively correlated with $\mathrm{log}(\sigma )$, $\mathrm{log}({M}_{\star })$, and $\mathrm{log}({M}_{\mathrm{dyn}})$. From the combined sample, we show that the IMF in the centers of more massive ETGs is more bottom heavy. In addition, we find that $\mathrm{log}({\alpha }_{\mathrm{IMF}})$ is positively correlated with both [Mg/Fe] and the estimated total metallicity [Z/H]. We find suggestive evidence that the effective stellar surface density Σ Kroupa might be responsible for the variation of α IMF . We conclude that σ, [Mg/Fe], and [Z/H] are the primary drivers of the global stellar IMF variation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

60 Microlensing Events from the Three Years of Zwicky Transient Facility Phase One

Abstract Microlensing events have historically been discovered throughout the Galactic bulge and plane by surveys designed solely for that purpose. We conduct the first multiyear search for microlensing events on the Zwicky Transient Facility (ZTF), an all-sky optical synoptic survey that observes the entire visible northern sky every few nights. We discover 60 high-quality microlensing events in the 3 yr of ZTF-I using the bulk lightcurves in the ZTF Public Data Release 5.19 of our events are found outside of the Galactic plane (∣ b ∣ ≥ 10°), nearly doubling the number of previously discovered events in the stellar halo from surveys pointed toward the Magellanic Clouds and the Andromeda galaxy. We also record 1558 ongoing candidate events as potential microlensing that can continue to be observed by ZTF-II for identification. The scalable and computationally efficient methods developed in this work can be applied to future synoptic surveys, such as the Vera C. Rubin Observatory’s Legacy Survey of Space and Time and the Nancy Grace Roman Space Telescope, as they attempt to find microlensing events in even larger and deeper data sets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Betelgeuse as a Merger of a Massive Star with a Companion

We investigate the merger between a 16M ⊙ star, on its way to becoming a red supergiant (RSG), and a 4M ⊙ main-sequence companion. Our study employs three-dimensional hydrodynamic simulations using the state-of-the-art adaptive mesh refinement code Octo-Tiger. The initially corotating binary undergoes interaction and mass transfer, resulting in the accumulation of mass around the companion and its subsequent loss through the second Lagrangian point (L2). The companion eventually plunges into the envelope of the primary, leading to its spin-up and subsequent merger with the helium core. We examine the internal structural properties of the post-merger star, as well as the merger environment and the outflow driven by the merger. Our findings reveal the ejection of approximately ∼0.6 M ⊙ of material in an asymmetric and somewhat bipolar outflow. We import the post-merger stellar structure into the MESA stellar evolution code to model its long-term nuclear evolution. In certain cases, the post-merger star exhibits persistent rapid equatorial surface rotation as it evolves in the H–R diagram toward the observed location of Betelgeuse. These cases demonstrate surface rotation velocities of a similar magnitude to those observed in Betelgeuse, along with a chemical composition resembling that of Betelgeuse. In other cases, efficient rotationally induced mixing leads to slower surface rotation. This pioneering study aims to model stellar mergers across critical timescales, encompassing dynamical, thermal, and nuclear evolutionary stages.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Quantifying Chemical and Kinematical Properties of Galactic Disks

We aim to quantify the chemical and kinematical properties of Galactic disks with a sample of 119,558 giant stars having abundances and 3D velocities taken or derived from the APOGEE DR17 and Gaia EDR3 catalogs. A Gaussian mixture model is employed to distinguish the high-α and low-α sequences along the metallicity by simultaneously using chemical and kinematical data. Four disk components are identified and quantified; they are named the hαmp, hαmr, lαmp, and lαmr disks and correspond to the high-α or low-α, and metal-poor or metal-rich properties. Combined with the spatial and stellar-age information, we confirm that they are well interpreted by the two-infall formation model. The first infall of turbulent gas quickly forms the hot and thick hαmp disk with consequent thinner hαmr and lαmr disks. Then the second gas accretion forms a thinner and outermost lαmp disk. We find that the inside-out and upside-down scenario does not only satisfy the overall Galactic disk formation of these two major episodes but is also presented in the formation sequence of the three inner disks. Importantly, we reveal the inverse age–[M/H] trend of the lαmr disk, which means its younger stars are more metal-poor, indicating that the rejuvenated gas from the second accretion gradually dominates later star formation. Meanwhile, the recently formed stars converge to [M/H] ~ –0.1 dex, demonstrating a sufficient mixture of gas from two infalls.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Planets Across Space and Time (PAST). II. Catalog and Analyses of the LAMOST–Gaia–Kepler Stellar Kinematic Properties

The Kepler telescope has discovered over 4000 planets (candidates) by searching ~200,000 stars over a wide range of distance (order of kpc) in our Galaxy. Characterizing the kinematic properties (e.g., Galactic component membership and kinematic age) of these Kepler targets (including the planet candidate hosts) is the first step toward studying Kepler planets in the Galactic context, which will reveal fresh insights into planet formation and evolution. In this paper, the second part of the Planets Across the Space and Time (PAST) series, by combining the data from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and Gaia and then applying the revised kinematic methods from PAST I, we present a catalog of kinematic properties (i.e., Galactic positions, velocities, and the relative membership probabilities among the thin disk, thick disk, Hercules stream, and the halo) as well as other basic stellar parameters for 35,835 Kepler stars. Further analyses of the LAMOST–Gaia–Kepler catalog demonstrate that our derived kinematic age reveals the expected stellar activity-age trend. Furthermore, we find that the fraction of thin (thick) disk stars increases (decreases) with the transiting planet multiplicity (N {sub p} = 0, 1, 2 and 3+) and the kinematic age decreases with N {sub p}, which could be a consequence of the dynamical evolution of planetary architecture with age. The LAMOST–Gaia–Kepler catalog will be useful for future studies on the correlations between the exoplanet distributions and the stellar Galactic environments as well as ages.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Emission line predictions for mock galaxy catalogues: a new differentiable and empirical mapping from DESI

ABSTRACT We present a simple, differentiable method for predicting emission line strengths from rest-frame optical continua using an empirically determined mapping. Extensive work has been done to develop mock galaxy catalogues that include robust predictions for galaxy photometry, but reliably predicting the strengths of emission lines has remained challenging. Our new mapping is a simple neural network implemented using the JAX Python automatic differentiation library. It is trained on Dark Energy Spectroscopic Instrument Early Release data to predict the equivalent widths (EWs) of the eight brightest optical emission lines (including H α, H β, [O ii], and [O iii]) from a galaxy’s rest-frame optical continuum. The predicted EW distributions are consistent with the observed ones when noise is accounted for, and we find Spearman’s rank correlation coefficient ρs > 0.87 between predictions and observations for most lines. Using a non-linear dimensionality reduction technique, we show that this is true for galaxies across the full range of observed spectral energy distributions. In addition, we find that adding measurement uncertainties to the predicted line strengths is essential for reproducing the distribution of observed line-ratios in the BPT diagram. Our trained network can easily be incorporated into a differentiable stellar population synthesis pipeline without hindering differentiability or scalability with GPUs. A synthetic catalogue generated with such a pipeline can be used to characterize and account for biases in the spectroscopic training sets used for training and calibration of photo-z’s, improving the modelling of systematic incompleteness for the Rubin Observatory LSST and other surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Substructure at High Speed. I. Inferring the Escape Velocity in the Presence of Kinematic Substructure

Abstract The local escape velocity provides valuable inputs to the mass profile of the galaxy, and requires understanding the tail of the stellar speed distribution. Following Leonard & Tremaine, various works have since modeled the tail of the stellar speed distribution as ∝ ( v esc − v ) k , where v esc is the escape velocity, and k is the slope of the distribution. In such studies, however, these two parameters were found to be largely degenerate and often a narrow prior is imposed on k in order to constrain v esc . Furthermore, the validity of the power-law form can breakdown in the presence of multiple kinematic substructures or other mis-modeled features in the data. In this paper, we introduce a strategy that for the first time takes into account the presence of kinematic substructure. We model the tail of the velocity distribution as a sum of multiple power laws as a way of introducing a more flexible fitting framework. Using mock data and data from FIRE simulations of Milky Way-like galaxies, we show the robustness of this method in the presence of kinematic structure that is similar to the recently discovered Gaia Sausage. In a companion paper, we present the new measurement of the escape velocity and subsequently the mass of the Milky Way using Gaia eDR3 data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The stellar mass of the Gaia -Sausage/Enceladus accretion remnant

The Gaia-Sausage/Enceladus (GS/E) structure is an accretion remnant that comprises a large fraction of the Milky Way’s stellar halo. We study GS/E using high-purity samples of kinematically selected stars from APOGEE DR16 and Gaia. Employing a novel framework to account for kinematic selection biases using distribution functions, we fit density profiles to these GS/E samples and measure their masses. We find that GS/E has a shallow density profile in the inner Galaxy, with a break between 15 and 25 kpc beyond which the profile steepens. We also find that GS/E is triaxial, with axis ratios 1:0.55:0.45 (nearly prolate), and the major axis is oriented about 80° from the Sun–Galactic centre line and 16° above the plane. We measure a stellar mass for GS/E of $1.45\, ^{+0.92}_{-0.51}\, \mathrm{(stat.)}\, ^{+0.13}_{-0.37} \mathrm{(sys.)}\ \times 10^{8}$ M ⊙ . Our mass estimate is lower than others in the literature, a finding we attribute to the excellent purity of the samples we work with. We also fit a density profile to the entire Milky Way stellar halo, finding a mass in the range of 6.7–8.4 × 10 8 M ⊙ , and implying that GS/E could make up as little as 15–25 per cent of the mass of the Milky Way stellar halo. Our lower stellar mass combined with standard stellar mass-to-halo mass relations implies that GS/E constituted a minor 1:8 mass-ratio merger at the time of its accretion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Estimate of the Mass and Radial Profile of the Orphan–Chenab Stream's Dwarf-galaxy Progenitor Using MilkyWay@home

We fit the mass and radial profile of the Orphan–Chenab Stream's (OCS) dwarf-galaxy progenitor by using turnoff stars in the Sloan Digital Sky Survey and the Dark Energy Camera to constrain N-body simulations of the OCS progenitor falling into the Milky Way on the 1.5 PetaFLOPS MilkyWay@home distributed supercomputer. We infer the internal structure of the OCS's progenitor under the assumption that it was a spherically symmetric dwarf galaxy composed of a stellar system embedded in an extended dark matter halo. We optimize the evolution time, the baryonic and dark matter scale radii, and the baryonic and dark matter masses of the progenitor using a differential evolution algorithm. The likelihood score for each set of parameters is determined by comparing the simulated tidal stream to the angular distribution of OCS stars observed in the sky. We fit the total mass of the OCS's progenitor to (2.0 ± 0.3) × 10 7 M ⊙ with a mass-to-light ratio of γ = 73.5 ± 10.6 and (1.1 ± 0.2) × 10 6 M ⊙ within 300 pc of its center. Within the progenitor's half-light radius, we estimate a total mass of (4.0 ± 1.0) × 10 5 M ⊙ . We also fit the current sky position of the progenitor's remnant to be (α, δ) = ((166.0 ± 0.9)°, (–11.1 ± 2.5)°) and show that it is gravitationally unbound at the present time. The measured progenitor mass is on the low end of previous measurements and, if confirmed, lowers the mass range of ultrafaint dwarf galaxies. Our optimization assumes a fixed Milky Way potential, OCS orbit, and radial profile for the progenitor, ignoring the impact of the Large Magellanic Cloud.

79 ASTRONOMY AND ASTROPHYSICS↗