Abstract The physical processes responsible for shaping how galaxies form and quench over time leave imprints on both the spatial (galaxy morphology) and temporal (star formation history; SFH) ...tracers that we use to study galaxies. While the morphology–SFR connection is well studied, the correlation with past star formation activity is not as well understood. To quantify this, we present Katachi (形), an interpretable convolutional neural network framework that learns the connection between the factors regulating star formation in galaxies on different spatial and temporal scales. Katachi is trained on 9904 galaxies at 0.02 < z < 0.1 in the SDSS-IV MaNGA DR17 sample to predict stellar mass ( M * ; root mean square error (RSME) 0.22 dex), current star formation rate (SFR; RMSE 0.31 dex), and half-mass time ( t 50 ; RMSE 0.23 dex). This information allows us to reconstruct nonparametric SFHs for each galaxy from gri imaging alone. To quantify the morphological features informing the SFH predictions, we use SHAP (SHapley Additive exPlanations). We recover the expected trends of M * governed by the growth of galaxy bulges, as well as those of SFR correlating with spiral arms and other star-forming regions. We also find that the SHAP maps of D4000 are more complex than those of M * and SFR, and that morphology is correlated with t 50 even at fixed mass and SFR. Katachi serves as a scalable public framework to predict galaxy properties from large imaging surveys including Rubin, Roman, and Euclid, with large data sets of high signal-to-noise ratio imaging across limited photometric bands.
Context. High-resolution 3D maps of interstellar dust are critical for probing the underlying physics shaping the structure of the interstellar medium, and for foreground correction of astrophysical ...observations affected by dust. Aims. We aim to construct a new 3D map of the spatial distribution of interstellar dust extinction out to a distance of 1.25 kpc from the Sun. Methods. We leveraged distance and extinction estimates to 54 million nearby stars derived from the Gaia BP/RP spectra. Using the stellar distance and extinction information, we inferred the spatial distribution of dust extinction. We modeled the logarithmic dust extinction with a Gaussian process in a spherical coordinate system via iterative charted refinement and a correlation kernel inferred in previous work. In total, our posterior has over 661 million degrees of freedom. We probed the posterior distribution using the variational inference method MGVI. Results. Our 3D dust map has an angular resolution of up to 14′ ( N side = 256), and we achieve parsec-scale distance resolution, sampling the dust in 516 logarithmically spaced distance bins spanning 69 pc to 1250 pc. We generated 12 samples from the variational posterior of the 3D dust distribution and release the samples alongside the mean 3D dust map and its corresponding uncertainty. Conclusions. Our map resolves the internal structure of hundreds of molecular clouds in the solar neighborhood and will be broadly useful for studies of star formation, Galactic structure, and young stellar populations. It is available for download in a variety of coordinate systems online and can also be queried via the publicly available dustmaps Python package.
ABSTRACT
The North Polar Spur (NPS) is one of the largest structures observed in the Milky Way in both the radio and soft X-rays. While several predictions have been made regarding the origin of the ...NPS, modelling the structure is difficult without precise distance constraints. In this paper, we determine accurate distances to the southern terminus of the NPS and towards latitudes ranging up to 55°. First, we fit for the distance and extinction to stars towards the NPS using optical and near-infrared photometry and Gaia Data Release 2 astrometry. We model these per-star distance–extinction estimates as being caused by dust screens at unknown distances, which we fit for using a nested sampling algorithm. We then compare the extinction to the Spur derived from our 3D dust modelling with integrated independent measures from XMM–Newton X-ray absorption and H i column density measures. We find that we can account for nearly 100 per cent of the total column density of the NPS as lying within 140 pc for latitudes >26° and within 700 pc for latitudes <11°. Based on the results, we conclude that the NPS is not associated with the Galactic Centre or the Fermi bubbles. Instead, it is likely associated, especially at higher latitudes, with the Scorpius–Centaurus association.
ABSTRACT
We present a hierarchical Bayesian inference approach to estimating the structural properties and the phase-space centre of a globular cluster (GC) given the spatial and kinematic ...information of its stars based on lowered isothermal cluster models. As a first step towards more realistic modelling of GCs, we built a differentiable, accurate emulator of the lowered isothermal distribution function using interpolation. The reliable gradient information provided by the emulator allows the use of Hamiltonian Monte Carlo methods to sample large Bayesian models with hundreds of parameters, thereby enabling inference on hierarchical models. We explore the use of hierarchical Bayesian modelling to address several issues encountered in observations of GC including an unknown GC centre, incomplete data, and measurement errors. Our approach not only avoids the common technique of radial binning but also incorporates the aforementioned uncertainties in a robust and statistically consistent way. Through demonstrating the reliability of our hierarchical Bayesian model on simulations, our work lays out the foundation for more realistic and complex modelling of real GC data.
Abstract
With an eye towards the computational requirements of future large-scale surveys such as Euclid and Large Synoptic Survey Telescope (LSST) that will require photometric redshifts (photo-z’s) ...for ≳ 109 objects, we investigate a variety of ways that ‘fuzzy archetypes’ can be used to improve photometric redshifts and explore their respective statistical interpretations. We characterize their relative performance using an idealized LSST ugrizY and Euclid YJH mock catalogue of 10 000 objects spanning z = 0–6 at Y = 24 mag. We find most schemes are able to robustly identify redshift probability distribution functions that are multimodal and/or poorly constrained. Once these objects are flagged and removed, the results are generally in good agreement with the strict accuracy requirements necessary to meet Euclid weak lensing goals for most redshifts between 0.8 ≲ z ≲ 2. These results demonstrate the statistical robustness and flexibility that can be gained by combining template-fitting and machine-learning methods and provide useful insights into how astronomers can further exploit the colour–redshift relation.
ABSTRACT We present the first detailed chemical-abundance analysis of stars from the dwarf-galaxy stellar stream Wukong/LMS-1 covering a wide metallicity range ($-3.5 \lt \rm Fe/H \lesssim -1.3$). We ...find abundance patterns that are effectively indistinguishable from the bulk of Indus and Jhelum, a pair of smaller stellar streams proposed to be dynamically associated with Wukong/LMS-1. We confirmed a carbon-enhanced metal-poor star ($\rm C/Fe \gt +0.7$ and $\rm Fe/H \sim -2.9$) in Wukong/LMS-1 with strong enhancements in Sr, Y, and Zr, which is peculiar given its solar-level Ba/Fe. Wukong/LMS-1 stars have high abundances of α elements up to $\rm Fe/H \gtrsim -2$, which is expected for relatively massive dwarfs. Towards the high-metallicity end, Wukong/LMS-1 becomes α-poor, revealing that it probably experienced fairly standard chemical evolution. We identified a pair of N- and Na-rich stars in Wukong/LMS-1, reminiscent of multiple stellar populations in globular clusters. This indicates that this dwarf galaxy contained at least one globular cluster that was completely disrupted in addition to two intact ones previously known to be associated with Wukong/LMS-1, which is possibly connected to similar evidence found in Indus. From these ≥3 globular clusters, we estimate the total mass of Wukong/LMS-1 to be ${\approx }10^{10} \, \mathrm{M}_\odot$, representing ∼1 per cent of the present-day Milky Way. Finally, the Eu/Mg ratio in Wukong/LMS-1 continuously increases with metallicity, making this the first example of a dwarf galaxy where the production of r-process elements is clearly dominated by delayed sources, presumably neutron-star mergers.
Photometric Biases in Modern Surveys Portillo, Stephen K. N.; Speagle, Joshua S.; Finkbeiner, Douglas P.
The Astronomical journal,
04/2020, Letnik:
159, Številka:
4
Journal Article
Recenzirano
Odprti dostop
Many surveys use maximum-likelihood (ML) methods to fit models when extracting photometry from images. We show that these ML estimators systematically overestimate the flux as a function of the ...signal-to-noise ratio and the number of model parameters involved in the fit. This bias is substantially worse for resolved sources: while a 1% bias is expected for a 10 point source, a 10 resolved galaxy with a simplified Gaussian profile suffers a 2.5% bias. This bias also behaves differently depending how multiple bands are used in the fit: simultaneously fitting all bands leads the flux bias to become roughly evenly distributed between them, while fixing the position in "non-detection" bands (i.e., forced photometry) gives flux estimates in those bands that are biased low, compounding a bias in derived colors. We show that these effects are present in idealized simulations, outputs from the Hyper Suprime-Cam fake-object pipeline (SynPipe), and observations from Sloan Digital Sky Survey Stripe 82. Prescriptions to correct for the ML bias in flux, and its uncertainty, are provided.
The tidal disruption of the Sagittarius dwarf galaxy has generated a spectacular stream of stars wrapping around the entire Galaxy. We use data from Gaia and the H3 Stellar Spectroscopic Survey to ...identify high-quality Sagittarius members based on their angular momenta. The H3 Survey is largely unbiased in metallicity, and so our sample of Sagittarius members is similarly unbiased. Stream stars span a wide range in Fe/H from −0.2 to −3.0, with a mean overall metallicity of . We identify a strong metallicity dependence to the kinematics of the stream members. At Fe/H > −0.8 nearly all members belong to the well-known cold ( ) leading and trailing arms. At intermediate metallicities (−1.9 < Fe/H < −0.8) a significant population (24%) emerges of stars that are kinematically offset from the cold arms. These stars also appear to have hotter kinematics. At the lowest metallicities (Fe/H −2), the majority of stars (69%) belong to this kinematically offset diffuse population. Comparison to simulations suggests that the diffuse component was stripped from the Sagittarius progenitor at earlier epochs, and therefore resided at larger radius on average than the colder metal-rich component. We speculate that this kinematically diffuse, low-metallicity population is the stellar halo of the Sagittarius progenitor system.
ABSTRACT
This work presents the Globular cluster Extra-tidal Mock Star (GEMS) catalogue of extra-tidal stars and binaries created via three-body dynamical encounters in globular cluster cores. Using ...the particle-spray code Corespray, we sample $N=50\, 000$ extra-tidal stars and escaped recoil binaries for 159 Galactic globular clusters. Sky positions, kinematics, stellar properties, and escape information are provided for all simulated stars. Stellar orbits are integrated in seven different static and time-varying Milky Way gravitational potential models where the structure of the disc, perturbations from the Large Magellanic Cloud and the mass and sphericity of the Milky Way’s dark matter halo are all investigated. We find that the action coordinates of the mock extra-tidal stars are largely Galactic model independent, where minor offsets and broadening of the distributions between models are likely due to interactions with substructure. Importantly, we also report the first evidence for stellar stream contamination by globular cluster core stars and binaries for clusters with pericentre radii larger than five kiloparsecs. Finally, we provide a quantitative tool that uses action coordinates to match field stars to host clusters with probabilities. Ultimately, combining data from the GEMS catalogue with information of observed stars will allow for association of extra-tidal field stars with any Galactic globular cluster; a requisite tool for understanding population-level dynamics and evolution of clusters in the Milky Way.
Here, we present high-fidelity cosmology results from a blinded joint analysis of galaxy-galaxy weak lensing (ΔΣ) and projected galaxy clustering (wp) measured from the Hyper Suprime-Cam Year-1 ...(HSC-Y1) data and spectroscopic Sloan Digital Sky Survey (SDSS) galaxy catalogs in the redshift range 0.15 < z < 0.7. We define luminosity-limited samples of SDSS galaxies to serve as the tracers of wp in three spectroscopic redshift bins, and as the lens samples for ΔΣ. For the ΔΣ measurements, we select a single sample of 4×106 source galaxies over 140 deg2 from HSC-Y1 with photometric redshifts (photo z) greater than 0.75, enabling a better handle of photo- z errors by comparing the ΔΣ amplitudes for the three lens redshift bins. The deep, high-quality HSC-Y1 data enable significant detections of the ΔΣ signals, with integrated signal-to-noise ratio S/N ~ 15 in the range 3 ≤ R/h–1 Mpc ≤ 30 for the three lens samples, despite the small area coverage. For cosmological parameter inference, we use an input galaxy-halo connection model built on the dark emulator package (which uses an ensemble set of high-resolution N-body simulations and enables fast, accurate computation of the clustering observables) with a halo occupation distribution that includes nuisance parameters to marginalize over modeling uncertainties. We model the ΔΣ and wp measurements on scales from R≃3 and 2h–1 Mpc , respectively, up to 30 h–1 Mpc (therefore excluding the baryon acoustic oscillations information) assuming a flat Λ CDM cosmology, marginalizing over about 20 nuisance parameters and demonstrating the robustness of our results to them. With various tests using mock catalogs described in Miyatake et al. preceding paper, Phys. Rev. D 106, 083519 (2022), we show that any bias in the clustering amplitude S8 ≡ σ8(Ωm/0.3)0.5 due to uncertainties in the galaxy-halo connection is less than ~ 50 % of the statistical uncertainty of S8, unless the assembly biaseffect is unexpectedly large. Our best-fit models have S8 = 0.795$_{-0.042}^{+0.049}$ (mode and 68% credible interval) for the flat Λ CDM model; we find tighter constraints on the quantity S8(α = 0.17)≡σ8(Ωm/0.3)0.17 = 0.745$_{-0.031}^{+0.039}$.