An official website of the United States government
Abstract Neutron stars provide a unique opportunity to study strongly interacting matter under extreme density conditions. The intricacies of matter inside neutron stars and their equation of state are not directly visible, but determine bulk properties, such as mass and radius, which affect the star's thermal X-ray emissions. However, the telescope spectra of these emissions are also affected by the stellar distance, hydrogen column, and effective surface temperature, which are not always well-constrained. Uncertainties on these nuisance parameters must be accounted for when making a robust estimation of the equation of state. In this study, we develop a novel methodology that, for the first time, can infer the full posterior distribution of both the equation of state and nuisance parameters directly from telescope observations. This method relies on the use of neural likelihood estimation, in which normalizing flows use samples of simulated telescope data to learn the likelihood of the neutron star spectra as a function of these parameters, coupled with Hamiltonian Monte Carlo methods to efficiently sample from the corresponding posterior distribution. Our approach surpasses the accuracy of previous methods, improves the interpretability of the results by providing access to the full posterior distribution, and naturally scales to a growing number of neutron star observations expected in the coming years.
more »
« less
Deducing neutron star equation of state from telescope spectra with machine-learning-derived likelihoods
The interiors of neutron stars reach densities and temperatures beyond the limits of terrestrial experiments, providing vital laboratories for probing nuclear physics. While the star's interior is not directly observable, its pressure and density determine the star's macroscopic structure which affects the spectra observed in telescopes. The relationship between the observations and the internal state is complex and partially intractable, presenting difficulties for inference. Previous work has focused on the regression from stellar spectra of parameters describing the internal state. We demonstrate a calculation of the full likelihood of the internal state parameters given observations, accomplished by replacing intractable elements with machine learning models trained on samples of simulated stars. Our machine-learning-derived likelihood allows us to performmaximum a posterioriestimation of the parameters of interest, as well as full scans. We demonstrate the technique by inferring stellar mass and radius from an individual stellar spectrum, as well as equation of state parameters from a set of spectra. Our results are more precise than pure regression models, reducing the width of the parameter residuals by 11.8% in the most realistic scenario. The neural networks will be released as a tool for fast simulation of neutron star properties and observed spectra.
more »
« less
- Award ID(s):
- 2206322
- PAR ID:
- 10484984
- Publisher / Repository:
- IOP Science
- Date Published:
- Journal Name:
- Journal of Cosmology and Astroparticle Physics
- Volume:
- 2023
- Issue:
- 12
- ISSN:
- 1475-7516
- Page Range / eLocation ID:
- 022
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Neutron stars provide a unique laboratory for studying matter at extreme pressures and densities. While there is no direct way to explore their interior structure, X-rays emitted from these stars can indirectly provide clues to the equation of state (EOS) of the superdense nuclear matter through the inference of the star's mass and radius. However, inference of EOS directly from a star's X-ray spectra is extremely challenging and is complicated by systematic uncertainties. The current state of the art is to use simulation-based likelihoods in a piece-wise method which relies on certain theoretical assumptions and simplifications about the uncertainties. It first infers the star's mass and radius to reduce the dimensionality of the problem, and from those quantities infer the EOS. We demonstrate a series of enhancements to the state of the art, in terms of realistic uncertainty quantification and a path towards circumventing the need for theoretical assumptions to infer physical properties with machine learning. We also demonstrate novel inference of the EOS directly from the high-dimensional spectra of observed stars, avoiding the intermediate mass-radius step. Our network is conditioned on the sources of uncertainty of each star, allowing for natural and complete propagation of uncertainties to the EOS.more » « less
-
Context. Cool stars, such as M giants, can only be analyzed in the near-infrared (NIR) regime due to the ubiquitous titanium oxide features in optical spectra of stars with T eff < 4000 K. In dust-obscured regions, the inner bulge and Galactic center region, the intrinsically bright M giants observed in the NIR are an optimal option for studying stellar abundances and the chemical evolution of stellar populations. Because of the uncertainties in photometric methods, a method for determining the stellar parameters for M giants from the NIR spectra themselves is needed. Aims. We develop a method for determining the stellar parameters for M giants from the NIR spectra. We validate the method by deriving the stellar parameters for nearby well-studied M giants with spectra from the spectral library of the Immersion GRating INfrared Spectrograph (IGRINS). We demonstrate the accuracy and precision of our method by determining the stellar parameters and α -element trends versus metallicity for solar neighborhood M giants. Methods. We carried out new observations of 44 M giant stars with IGRINS mounted on the Gemini South telescope. We also obtained the full H and K band IGRINS spectra of six nearby well-studied M giants at a spectral resolving power of R = 45 000 from the IGRINS spectral library. We used the tool called spectroscopy made easy in combination with one-dimensional (1D) model atmospheres in a radiative and convective scheme (MARCS) stellar atmosphere models to model the synthetic spectrum that fits the observed spectrum best. Results. The effective temperatures that we derive from our new method (tested for 3400 ≲ T eff ≲ 4000 K here) agree excellently with those of the six nearby well-studied M giants, which indicates that the accuracy is indeed high. For the 43 solar neighborhood M giants, our T eff , log g , [Fe/H], ξ micro , [C/Fe], [N/Fe], and [O/Fe] agree with APOGEE with mean differences and a scatter (our method – APOGEE) of −67±33 K, −0.31±0.15 dex, 0.02±0.05 dex, 0.22±0.13 km s −1 , −0.05±0.06 dex, 0.06±0.06 dex, and 0.02±0.09 dex, respectively. Furthermore, the tight offset with a small dispersion compared to the APOGEE T eff indicates a high precision in our derived temperatures and those derived from the APOGEE pipeline. The typical uncertainties in the stellar parameters are found to be ±100 K in T eff , ±0.2 dex in log g , ±0.1 dex in [Fe/H], and ±0.1 km s −1 in ξ micro . The α -element trends versus metallicity for Mg, Si, Ca, and Ti are consistent with the APOGEE DR17 trends for the same stars and with the GILD optical trends. We also find a clear enhancement in the abundances for thick-disk stars.more » « less
-
Abstract We present the third and final data release of the K2 Galactic Archaeology Program (K2 GAP) for Campaigns C1–C8 and C10–C18. We provide asteroseismic radius and mass coefficients,κRandκM, for ∼19,000 red giant stars, which translate directly to radius and mass given a temperature. As such, K2 GAP DR3 represents the largest asteroseismic sample in the literature to date. K2 GAP DR3 stellar parameters are calibrated to be on an absolute parallactic scale based on Gaia DR2, with red giant branch and red clump evolutionary state classifications provided via a machine-learning approach. Combining these stellar parameters with GALAH DR3 spectroscopy, we determine asteroseismic ages with precisions of ∼20%–30% and compare age-abundance relations to Galactic chemical evolution models among both low- and high-αpopulations forα, light, iron-peak, and neutron-capture elements. We confirm recent indications in the literature of both increased Ba production at late Galactic times as well as significant contributions tor-process enrichment from prompt sources associated with, e.g., core-collapse supernovae. With an eye toward other Galactic archeology applications, we characterize K2 GAP DR3 uncertainties and completeness using injection tests, suggesting that K2 GAP DR3 is largely unbiased in mass/age, with uncertainties of 2.9% (stat.) ± 0.1% (syst.) and 6.7% (stat.) ± 0.3% (syst.) inκRandκMfor red giant branch stars and 4.7% (stat.) ± 0.3% (syst.) and 11% (stat.) ± 0.9% (syst.) for red clump stars. We also identify percent-level asteroseismic systematics, which are likely related to the time baseline of the underlying data, and which therefore should be considered in TESS asteroseismic analysis.more » « less
-
Abstract We present the visual orbits of four spectroscopic binary stars, HD 61859, HD 89822, HD 109510, and HD 191692, using long baseline interferometry with the CHARA Array. We also obtained new radial velocities from echelle spectra using the APO 3.5 m, CTIO 1.5 m, and Fairborn Observatory 2.0 m telescopes. By combining the astrometric and spectroscopic observations, we solve for the full, three-dimensional orbits and determine the stellar masses to 1%–12% uncertainty and distances to 0.4%–6% uncertainty. We then estimate the effective temperature and radius of each component star through Doppler tomography and spectral energy distribution analyses. We found masses of 1.4–3.5 M ⊙ , radii of 1.5–4.7 R ⊙ , and temperatures of 6400–10,300 K. We then compare the observed stellar parameters to the predictions of the stellar evolution models, but found that only one of our systems fits well with the evolutionary models.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1088/1475-7516/2023/12/022
