Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MILKY WAY”

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 991 records · Page 55

popclass: A Python Package for Classifying Microlensing Events

popclass is a Python package that provides a flexible, probabilistic framework for classifying the lens of a gravitational microlensing event. Gravitational microlensing occurs when a massive foreground object (e.g., a star, white dwarf or black hole) passes in front of and lenses the light from a distant background source. This causes an apparent brightening, and shift in position, of the background source. In most cases, characteristics of the microlensing signal do not contain enough information to definitively identify the lens type. Different lens types lie in different but overlapping regions of the characteristics of the microlensing signal. For example, black holes tend to be more massive than stars and therefore cause microlensing signals that are longer. Current Galactic simulations enable us to predict where different lens types lie in the observational space and can therefore be used to classify events (e.g., Lam et al., 2020). popclass allows the user to match the characteristics of a microlensing signal with a simulation of the Galaxy to calculate lens type probabilities for the event (see Figure 1). Constraints on any microlensing signal properties and any Galactic model can be used. popclass comes with an interface to ArviZ (Kumar et al., 2019) and PyMultiNest (Buchner et al., 2014) for microlensing signal constraints, as well as pre-loaded Galactic models, plotting functionality, and methods to quantify the classification uncertainty. The probabilistic framework for popclass was developed in Perkins et al. (2024), used in Fardeen et al. (2024) and has been applied to classifying events in Kaczmarek et al. (2025).

97 MATHEMATICS AND COMPUTING↗

Efficient Detection of Emission-line Galaxies in the Cl0016+1609 and MACSJ1621.4+3810 Supercluster Filaments Using SITELLE

We observe a system of filaments and clusters around Cl0016+1609 and MACSJ1621.4+3810 using the SITELLE Fourier transform spectrograph at the Canada–France–Hawaii Telescope. For Cl0016+1609 (z = 0.546), the observations span an 11.8 Mpc × 4.3 Mpc region along an eastern filament that covers the main cluster core, as well as two 4.3 Mpc × 4.3 Mpc regions that each cover southern subclumps. For MACSJ1621.4+3810 (z = 0.465), 3.9 Mpc × 3.9 Mpc around the main cluster core is covered. We present the frequency and location of the emission-line galaxies and their emission-line images, and calculate the star formation rates, specific star formation rates, and merger statistics. In Cl0016+1609, we find 13 [O ii] 3727 Å emitting galaxies with star formation rates between 0.2 and 14.0 M {sub ⊙} yr{sup −1}. Of these, 91{sub −10}{sup +3}% are found in regions with moderate local galaxy density, avoiding the dense cluster cores. These galaxies follow the main filament of the superstructure and are mostly blue and disky, with several showing close companions and merging morphologies. In MACSJ1621.4+3810, we find 10 emission-line sources. All are blue (100{sub −15}{sup +0}%), with 40{sub −12}{sup +16}% classified as disky and 60{sub −16}{sup +12}% as merging systems. Eight avoid the cluster core (80{sub −17}{sup +7}%), but two (20{sub −7}{sup +17}%) are found near high-density regions, including the brightest cluster galaxy (BCG). These observations push the spectroscopic study of galaxies in filaments beyond z ∼ 0.3 to z ∼ 0.5. Their efficient confirmation is paramount to their usefulness as more galaxy surveys come online.

47 OTHER INSTRUMENTATION↗

The Open Cluster Chemical Abundances and Mapping Survey. VI. Galactic Chemical Gradient Analysis from APOGEE DR17

The goal of the Open Cluster Chemical Abundances and Mapping (OCCAM) survey is to constrain key Galactic dynamic and chemical evolution parameters by the construction and analysis of a large, comprehensive, uniform data set of infrared spectra for stars in hundreds of open clusters. This sixth contribution from the OCCAM survey presents analysis of SDSS/APOGEE Data Release 17 (DR17) results for a sample of stars in 150 open clusters, 94 of which we designate to be "high-quality" based on the appearance of their color–magnitude diagram. We find the APOGEE DR17-derived [Fe/H] values to be in good agreement with those from previous high-resolution spectroscopic open cluster abundance studies. Using a subset of the high-quality sample, the Galactic abundance gradients were measured for 16 chemical elements, including [Fe/H], for both Galactocentric radius ($R$ GC ) and guiding center radius ($R$ guide ). We find an overall Galactic [Fe/H] versus RGC gradient of -0.073 ± 0.002 dex kpc -1 over the range of 6 > $R$ GC < 11.5 kpc, and a similar gradient is found for [Fe/H] versus $R$ guide . Significant Galactic abundance gradients are also noted for O, Mg, S, Ca, Mn, Na, Al, K, and Ce. Our large sample additionally allows us to explore the evolution of the gradients in four age bins for the remaining 15 elements.

79 ASTRONOMY AND ASTROPHYSICS↗

${S}^5$: The destruction of a bright dwarf galaxy as revealed by the chemistry of the Indus stellar stream

The recently discovered Indus stellar stream exhibits a diverse chemical signature compared to what is found for most other streams due to the abundances of two outlier stars, Indus$\_$0 and Indus$\_$13. Indus$\_$13, exhibits an extreme enhancement in rapid neutron-capture ($r$-)process elements with $\mathrm{[Eu/Fe]} = +1.81$. It thus provides direct evidence of the accreted nature of $r$-process enhanced stars. In this paper we present a detailed chemical analysis of the neutron-capture elements in Indus$\_$13, revealing the star to be slightly actinide poor. The other outlier, Indus$\_0$, displays a globular cluster-like signature with high N, Na, and Al abundances, while the rest of the Indus stars show abundances compatible with a dwarf galaxy origin. Hence, Indus$\_0$ provides the first chemical evidence of a fully disrupted dwarf containing a globular cluster. We use the chemical signature of the Indus stars to discuss the nature of the stream progenitor which was likely a chemically evolved system, with a mass somewhere in the range from Ursa Minor to Fornax.

79 ASTRONOMY AND ASTROPHYSICS↗

Spatial and Kinematic Clustering of Stars in the Galactic Disk

The Galactic disk is expected to be spatially and kinematically clustered on many scales due to both star formation and the Galactic potential. In this work we calculate the spatial and kinematic two-point correlation functions (TPCF) using a sample of 1.7 × 10 6 stars with radial velocities from Gaia DR2. Clustering is detected on spatial scales of 1–300 pc and a velocity scale of 15 km s –1 . After removing bound structures, the data have a power-law index of γ ≈ –1 for 1 pc < Δr < 100 pc and γ ≲ –1.5 for Δr > 100 pc. We interpret these results with the aid of a star-by-star simulation of the Galaxy, in which stars are born in clusters orbiting in a realistic potential that includes spiral arms, a bar, and giant molecular clouds. We find that the simulation largely agrees with the observations at most spatial and kinematic scales. In detail, the TPCF in the simulation is shallower than the data at ≲20 pc scales, and steeper than the data at ≳30 pc. We also find a persistent clustering signal in the kinematic TPCF for the data at large Δv (>5 km s –1 ) that is not present in the simulations. We speculate that this mismatch between observations and simulations may be due to two processes: hierarchical star formation and transient spiral arms. Furthermore, we also predict that the addition of ages and metallicities measured with a precision of 50% and 0.05 dex, respectively, will enhance the clustering signal beyond current measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Use of Field RR Lyrae As Galactic Probes: IV. New Insights Into and Around the Oosterhoff Dichotomy

We discuss the largest and most homogeneous spectroscopic data set of field RR Lyrae variables (RRLs) available to date. We estimated abundances using both high-resolution and low-resolution ($Δ$S method) spectra for fundamental (RRab) and first overtone (RRc) RRLs. The iron abundances for 7941 RRLs were supplemented with similar estimates that are available in the literature, ending up with 9015 RRLs (6150 RRab, 2865 RRc). The metallicity distribution shows a mean value of $\langle$[Fe/H]$\rangle$ = -1.51 ± 0.01, and $σ$(standard deviation) = 0.41 dex with a long metal-poor tail approaching [Fe/H] ≃ - 3 and a sharp metal-rich tail approaching solar iron abundance. The RRab variables are more metal-rich ( $\langle$[Fe/H]$\rangle$ ab = -1.48 ± 0.01, $σ$ = 0.41 dex) than RRc variables ($\langle$[Fe/H]$\rangle$ c = -1.58 ± 0.01, $σ$ = 0.40 dex). The relative fraction of RRab variables in the Bailey diagram (visual amplitude versus period) located along the short-period (more metal-rich) and the long-period (more metalpoor) sequences are 80% and 20%, while RRc variables display an opposite trend, namely 30% and 70%, respectively. We found that the pulsation period of both RRab and RRc variables steadily decreases when moving from the metal-poor to the metal-rich regime. Furthermore, the visual amplitude shows the same trend, but RRc amplitudes are almost two times more sensitive than RRab amplitudes to metallicity. We also investigated the dependence of the population ratio (N $c$ /N tot ) of field RRLs on the metallicity and we found that the distribution is more complex than in globular clusters. The population ratio steadily increases from ~0.25 to ~0.36 in the metal-poor regime, it decreases from ~0.36 to ~0.18 for -1.8 ≤ [Fe/H] ≤ -0.9 and it increases to a value of ~0.3 approaching solar iron abundance.

79 ASTRONOMY AND ASTROPHYSICS↗

The Atacama Cosmology Telescope: Microwave Intensity and Polarization Maps of the Galactic Center

Herein we present arcminute-resolution intensity and polarization maps of the Galactic center made with the Atacama Cosmology Telescope. The maps cover a 32 deg 2 field at 98, 150, and 224 GHz with |l| ≤ 4°, |b| ≤ 2°. We combine these data with Planck observations at similar frequencies to create coadded maps with increased sensitivity at large angular scales. With the coadded maps, we are able to resolve many known features of the Central Molecular Zone (CMZ) in both total intensity and polarization. We map the orientation of the plane-of-sky component of the Galactic magnetic field inferred from the polarization angle in the CMZ, finding significant changes in morphology in the three frequency bands as the underlying dominant emission mechanism changes from synchrotron to dust emission. Selected Galactic center sources, including Sgr A*, the Brick molecular cloud (G0.253+0.016), the Mouse pulsar wind nebula (G359.23-0.82), and the Tornado supernova remnant candidate (G357.7-0.1), are examined in detail. These data illustrate the potential for leveraging ground-based cosmic microwave background polarization experiments for Galactic science.

79 ASTRONOMY AND ASTROPHYSICS↗

From the Fire: A Deeper Look at the Phoenix Stream

We use 6 yr of data from the Dark Energy Survey to perform a detailed photometric characterization of the Phoenix stellar stream, a 15° long, thin, dynamically cold, low-metallicity stellar system in the Southern Hemisphere. We use natural splines, a nonparametric modeling technique, to simultaneously fit the stream track, width, and linear density. This updated stream model allows us to improve measurements of the heliocentric distance (17.4 ± 0.1 (stat.) ±0.8 (sys.) kpc) and distance gradient (−0.009 ± 0.006 kpc deg$^{−1}$) of Phoenix, which corresponds to a small change of 0.13 ± 0.09 kpc in heliocentric distance along the length of the stream. We measure linear intensity variations on degree scales, as well as deviations in the stream track on ∼2° scales, suggesting that the stream may have been disturbed during its formation and/or evolution. We recover three peaks and one gap in linear intensity along with fluctuations in the stream track. Compared to other thin streams, the Phoenix stream shows more fluctuations and, consequently, the study of Phoenix offers a unique perspective on gravitational perturbations of stellar streams. We discuss possible sources of perturbations to Phoenix, including baryonic structures in the Galaxy and dark matter subhalos.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetic Misalignment of Interstellar Dust Filaments

We present evidence for scale-independent misalignment of interstellar dust filaments and magnetic fields. We estimate the misalignment by comparing millimeter-wave dust-polarization measurements from Planck with filamentary structures identified in neutral-hydrogen (H I ) measurements from H I 4PI. We find that the misalignment angle displays a scale independence (harmonic coherence) for features larger than the H I 4PI beamwidth (16$\acute{.}$2). We additionally find a spatial coherence on angular scales of ${ \mathcal O }(1^\circ )$. We present several misalignment estimators formed from the auto- and cross-spectra of dust-polarization and H I -based maps, and we also introduce a map-space estimator. Applied to large regions of the high-Galactic-latitude sky, we find a global misalignment angle of ~2°, which is robust to a variety of masking choices. By dividing the sky into small regions, we show that the misalignment angle correlates with the parity-violating TB cross-spectrum measured in the Planck dust maps. The misalignment paradigm also predicts a dust EB signal, which is of relevance in the search for cosmic birefringence but as yet undetected; the measurements of EB are noisier than those of TB, and our correlations of EB with misalignment angle are found to be weaker and less robust to masking choices. We also introduce an H I -based dust-polarization template constructed from the Hessian matrix of the H I intensity, which is found to correlate more strongly than previous templates with Planck dust B modes.

79 ASTRONOMY AND ASTROPHYSICS↗

A Swing of the Pendulum: The Chemodynamics of the Local Stellar Halo Indicate Contributions from Several Radial Merger Events

We find that the chemical abundances and dynamics of APOGEE and GALAH stars in the local stellar halo are inconsistent with a scenario in which the inner halo is primarily composed of debris from a single massive, ancient merger event, as has been proposed to explain the Gaia-Enceladus/Gaia Sausage (GSE) structure. The data contain trends of chemical composition with energy that are opposite to expectations for a single massive, ancient merger event, and multiple chemical evolution paths with distinct dynamics are present. We use a Bayesian Gaussian mixture model regression algorithm to characterize the local stellar halo, and find that the data are fit best by a model with four components. We interpret these components as the Virgo Radial Merger (VRM), Cronus, Nereus, and Thamnos; however, Nereus and Thamnos likely represent more than one accretion event because the chemical abundance distributions of their member stars contain many peaks. Although the Cronus and Thamnos components have different dynamics, their chemical abundances suggest they may be related. We show that the distinct low- and high-α halo populations from Nissen & Schuster are explained by VRM and Cronus stars, as well as some in situ stars. Because the local stellar halo contains multiple substructures, different popular methods of selecting GSE stars will actually select different mixtures of these substructures, which may change the apparent chemodynamic properties of the selected stars. We also find that the Splash stars in the Solar region are shifted to higher v $\phi$ and slightly lower [Fe/H] than previously reported.

79 ASTRONOMY AND ASTROPHYSICS↗

A Deep, High-angular-resolution 3D Dust Map of the Southern Galactic Plane

We present a deep, high-angular-resolution 3D dust map of the southern Galactic plane over 239° < l < 6° and ∣ b∣ < 10° built on photometry from the DECaPS2 survey, in combination with photometry from VISTA Variables in the Via Lactea, the Two Micron All Sky Survey, and “Unofficial” Wide-field Infrared Survey Explorer and parallaxes from Gaia Data Release 3 where available. To construct the map, we first infer the distance, extinction, and stellar types of over 700 million stars using the brutus stellar inference framework with a set of theoretical MESA Isochrone and Stellar Tracks (MIST) stellar models. Our resultant 3D dust map has an angular resolution of 1′ , roughly an order of magnitude finer than existing 3D dust maps and comparable to the angular resolution of the Herschel 2D dust emission maps. We detect complexes at the range of distances associated with the Sagittarius-Carina and Scutum-Centaurus arms in the fourth quadrant, as well as more distant structures out to a maximum reliable distance of d ≈ 10 kpc from the Sun. The map is sensitive up to a maximum extinction of roughly A V ≈ 12 mag. We publicly release both the stellar catalog and the 3D dust map, the latter of which can easily be queried via the Python package dustmaps. When combined with the existing Bayestar19 3D dust map of the northern sky, the DECaPS 3D dust map fills in the missing piece of the Galactic plane, enabling extinction corrections over the entire disk ∣b∣ < 10°. Our map serves as a pathfinder for the future of 3D dust mapping in the era of LSST and Roman, targeting regimes accessible with deep optical and near-infrared photometry but often inaccessible with Gaia.

Milky Way galaxy↗

Reconstructing the Stripping History of the Sagittarius Stream with Neural Networks

The Sagittarius (Sgr) Stream is produced by the ongoing disruption of the Sgr dwarf spheroidal (dSph) galaxy and is thought to contain multiple wraps that were stripped during different pericentric passages. In this study, we introduce a neural-network–based method trained on N-body simulations to infer the stripping time of Sgr Stream stars directly from their phase-space coordinates. We combine spectroscopic data from SEGUE, APOGEE DR17, and LAMOST DR7 low-resolution spectroscopic (LRS) survey with Gaia EDR3 astrometry and distance estimates from the latest StarHorse catalog to identify high-quality Sgr Stream members. Applying our method to these stars, we measure a clear metallicity gradient with stripping time, well described by a linear relation with slope ∼0.3 dex Gyr −1 . We further predict the stripping times of globular clusters previously suggested to originate from the Sgr dSph. M 54, Terzan 7, Terzan 8, and Arp 2 exhibit stripping times consistent with being currently bound to the Sgr remnant. Pal 12, Whiting 1, and NGC 2419 are inferred to have been stripped 0.9 ± 0.1, 1.1 ± 0.2, and 2.1 ± 0.2 Gyr ago, respectively. For NGC 4147 and NGC 5634, whose membership in the Sgr system remains uncertain, our analysis suggests stripping times of 1.1 ± 0.4 and 1.1 ± 0.1 Gyr, respectively, if they are ultimately confirmed as genuine Sgr members. These results demonstrate that data-driven models of dynamical stripping histories offer a promising approach for reconstructing the formation and chemical evolution of the Sgr Stream.

79 ASTRONOMY AND ASTROPHYSICS↗

Nonthermal galactic emission below 10 megahertz

The Radio Astronomy Explorer 2 (RAE 2) lunar-orbiting satellite has provided measurements of the nonthermal galactic radio emission at frequencies below 10 MHz. Measurements of the emission spectra are presented for the center, anticenter, north polar, and south polar directions at 22 frequencies between 0.25 and 9.18 MHz. Survey maps of the spatial distribution of the observed low-frequency Galactic emission at 1.31, 2.20, 3.93, 4.70, 6.55, and 9.18 MHz are presented. The observations were obtained with the 229-m traveling-wave V antenna on this lunar-orbiting spacecraft. The improved frequency coverage offers additional insights into the structure of the local Galactic neighborhood.

Milky Way↗

Testing the Universality of the Stellar IMF with Chandra and HST

The stellar initial mass function (IMF), which is often assumed to be universal across unresolved stellar populations, has recently been suggested to be bottom-heavy for massive ellipticals. In these galaxies, the prevalence of gravity-sensitive absorption lines (e.g., Na I and Ca II) in their near-IR spectra implies an excess of low-mass (m < or approx. = 0.5 Stellar Mass) stars over that expected from a canonical IMF observed in low-mass ellipticals. A direct extrapolation of such a bottom-heavy IMF to high stellar masses (m > or approx. = 8 Stellar Mass) would lead to a corresponding deficit of neutron stars and black holes, and therefore of low-mass X-ray binaries (LMXBs), per unit near-IR luminosity in these galaxies. Peacock et al. searched for evidence of this trend and found that the observed number of LMXBs per unit K-band luminosity (N/LK) was nearly constant. We extend this work using new and archival Chandra X-ray Observatory and Hubble Space Telescope observations of seven low-mass ellipticals where N/LK is expected to be the largest and compare these data with a variety of IMF models to test which are consistent with the observed N/LK. We reproduce the result of Peacock et al., strengthening the constraint that the slope of the IMF at m > or approx. = 8 Stellar Mass must be consistent with a Kroupa-like IMF. We construct an IMF model that is a linear combination of a Milky Way-like IMF and a broken power-law IMF, with a steep slope (alpha1 = 3.84) for stars < 0.5 Stellar Mass (as suggested by near-IR indices), and that flattens out (alpha2 = 2.14) for stars > 0.5 Stellar Mass, and discuss its wider ramifications and limitations.

galaxies: elliptical and lenticular↗

Black Holes and Neutron Stars in Nearby Galaxies: Insights from NuSTAR

Nearby galaxy surveys have long classified X-ray binaries (XRBs) by the mass category of their donor stars (high-mass and low-mass). The NuSTAR (Nuclear Spectroscopic Telescope Array) observatory, which provides imaging data at E greater than 10 keV (kiloelectronvolts), has enabled the classification of extragalactic XRBs by their compact object type: neutron star (NS) or black hole (BH). We analyzed NuSTAR/Chandra/XMM-Newton observations from a NuSTAR-selected sample of 12 galaxies within 5 Mpc (megaparsecs) having stellar masses (M (sub asterisk)) 10 (sup 7-11) solar masses and star formation rates (SFRs) approximately 0.01-15 solar masses per year. We detected 128 NuSTAR sources to a sensitivity of 10 (sup 38) ergs per second. Using NuSTAR color-intensity and color-color diagrams we classified 43 of these sources as candidate NSs and 47 as candidate BHs. We further subdivide BHs by accretion states (soft, intermediate, and hard) and NSs by weak (Z/Atoll) and strong (accreting pulsar) magnetic field. Using eight normal (Milky Way-type) galaxies in the sample, we confirm the relation between the SFR and galaxy X-ray point source luminosity in the 4-25 and 12-25 keV energy bands. We also constrained galaxy X-ray point source luminosity using the relation L (sub X) equals M (sub asterisk) plus Beta SFR, finding agreement with previous work. The X-ray luminosity function (XLF) of all sources in the 4-25 and 12-25 keV energy bands matches the alpha = 1.6 slope for high-mass XRBs. We find that NS XLFs suggest a decline beginning at the Eddington limit for a 1.4 solar mass NS, whereas the BH fraction shows an approximate monotonic increase in the 4-25 and 12-25 keV energy bands. We calculate the overall ratio of BH to NS to be approximately 1 for 4-25 keV and approximately 2 for 12-25 keV.

Vulic, N.↗

Ultraviolet and X-ray properties of Coma’s ultra-diffuse galaxies

Many ultra-diffuse galaxies (UDGs) have been discovered in the Coma cluster, and there is evidence that some, notably Dragonfly 44, have Milky Way-like dynamical masses despite dwarf-like stellar masses. We used X-ray, ultraviolet (UV), and optical data to investigate the star formation and nuclear activity in the Coma UDGs, and we obtained deep UV and X-ray data (Swift and XMM–Newton) for Dragonfly 44 to search for low-level star formation, hot circumgalactic gas, and the integrated emission from X-ray binaries. Among the Coma UDGs, we find UV luminosities consistent with quiescence but NUV−r colours indicating star formation in the past Gyr. This indicates that the UDGs were recently quenched. The r-band luminosity declines with projected distance from the Coma core. The Dragonfly 44 UV luminosity is also consistent with quiescence, with SFR<6×10^−4M(ʘ)/yr,and no X-rays are detected down to a sensitivity of 10^38erg/s. This rules out a hot corona with a M>10^8M(ʘ) within the virial radius, which would normally be expected for a dynamically massive galaxy. The absence of bright, low-mass X-ray binaries is consistent with the expectation from the Galaxy total stellar mass, but it is unlikely if most low-mass X-ray binaries form in globular clusters, as Dragonfly 44 has a very large population. Based on the UV and X-ray analysis, the Coma UDGs are consistent with quenched dwarf galaxies, although we cannot rule out a dynamically massive population.

Chris H. Lee↗

HaloSat Observations of Heliospheric Solar Wind Charge eXchange

X-ray emission from solar wind charge exchange (SWCX) produced in interplanetary space contaminates every astrophysical observation, regardless of the line-of-sight. Unfortunately, the primary SWCX emission lines also happen to be important diagnostics of astrophysical plasmas. Models of SWCX emission are limited by two main uncertainties: the local solar wind fluxes along the line-of-sight, and the charge exchange cross sections. The He cone, a localized density enhancement of helium neutrals, is the only heliospheric SWCX emission feature that is small enough and bright enough to be observationally isolated from the X-ray background and the broader SWCX emission. HaloSat, an X-ray CubeSat mission, has recently completed two series of specialized observations, near and far from the ecliptic plane, during two Earth transits of the He cone. These observations were used to test the predictions of a SWCX emission model against the emission observed at low ecliptic latitudes, where the solar wind data are monitored, and at high ecliptic latitudes, where the solar wind data are extrapolated. The measured SWCX emission for the set of observations near the ecliptic plane were consistent with the line intensities predicted by the model, but were underpredicted for the set of observations at high ecliptic latitude near the south ecliptic pole. Additionally, high temperature Galactic Halo emission components are reported for both spectral sets.

Solar Wind↗

Joint Modeling of GD-1 and C-19 as Old Streams

DESI observational data for the GD-1 and C-19 streams are compared to stream simulations in a common evolving multi-halo potential of a Milky Way-like galaxy based on a cosmological simulation. The goal is to find the best match of the stream velocity spread and the density power spectrum stream density to simulations having either CDM or WDM subhalos. The cocoon velocity width integrated over the length of the stream is independent of orbital blurring along the stream and the power spectrum integrates over the width of the stream, sidestepping the geometric details of the streams. Streams develop from star clusters inserted at $\simeq$1 Gyr after the Big Bang and evolved for 13 Gyr to their current orbital positions. Streams in a CDM subhalo population provide the best match to the velocity width, with streams younger than 10 Gyr ruled out as insufficiently hot. The progenitor star cluster masses, which determine the fraction of stars released at late times which comprise the stream core, are found to be $\simeq 8\times 10^4 M_\odot$ for GD-1 and $\simeq 4\times 10^4 M_\odot$ for C-19, although the mass depends on the star cluster half mass radius. Stream heating leads to stream lumpiness which is measurable for the relatively large and clean GD-1 dataset. The stream density power spectrum measured along the length of the DESI GD-1 sample is in good agreement with CDM simulations, with 1.7 to 1.9 times more power than WDM 7 keV and 5.5 keV simulations.

Carlberg, Raymond G. [Toronto U.] (ORCID:000000027↗