Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Milky Way Galaxy”

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

EDEN: Exploring Disks Embedded in N-body Simulations of Milky Way–mass Halos from Symphony

We investigate the impact of galactic disks on the tidal stripping of cold dark matter subhalos within Milky Way (MW)-mass halos (M vir ∼ 10 12 M ⊙ ) using a new simulation suite, EDEN. By re-simulating 45 MW-mass zoomed-in halos from the N-body Symphony compilation with embedded disk potentials, which evolve according to star formation histories predicted by the UniverseMachine model, we self-consistently tie disk growth to halo accretion rate and significantly expand the range of disk masses and formation histories studied. We use the particle-tracking-based subhalo finder Symfind to enhance the robustness of subhalo tracking. We find that disks near the median disk-to-halo mass ratio of our sample (M *,Disk /M vir,host = 2%) reduce subhalo peak mass functions within 100 kpc by about 10% for peak masses above 10 8 M ⊙ . Heavier, MW/M31-like disks (M *,Disk /M vir,host ≳ 5%) lead to a reduction of more than 40%. Subhalo abundance suppression is more pronounced near halo centers, with particularly enhanced stripping for subhalos accreted over 8 Gyr ago on orbits with pericenters <100 kpc. Suppression is further amplified when disk mass is increased within fixed halo and disk assembly histories. In all cases, the suppression we measure should be interpreted as stripping below the mass resolution limit rather than complete subhalo disruption. This study reshapes our understanding of the MW’s impact on its satellites, suggesting it strips subhalos more efficiently than typical MW-mass galaxies due to its larger disk-to-halo mass ratio and earlier disk formation.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1995

A far-infrared spectral line survey of 23 infrared-bright Galaxy nuclei

We present results from a KAO survey of fine-structure lines observed in 23 infrared-luminous galaxies. One or more of the following lines was observed and/or detected in each galaxy: (S III) 19, 33 microns, (Ne V) 24 microns, (O IV) 26 microns, (Fe II) 26 microns, (Si II) 35 microns, (O III) 52, 88 microns, (0 I) 63, 146 microns, (N III) 57 micro ns, (N II) 122, 205 microns, (C II) 158 microns. The galaxies span a wide range of morphologies (irregular to grand design), have varying metallicities, and include mergers, AGN's, and starburst systems. The observations were made beginning in 1988 using the facility Cryogenic Grating Spectrometer onboard the KAO at a typical resolution of approximately 60-140 km/s and with a 30-44 deg beam. We interpret the (C II) and (O I) fluxes, along with previous measurements of the IR continuum fluxes, in the context of photo dissociation region (PDR) models (Tielens & Hollenbach 1985; Wolfire et al. 1990). With these models, we obtain estimates of the typical interstellar UV fields incident on the line emitting regions (102-104 times the local interstellar radiation field) and the total masses (10(exp 7)-10(exp 8) Solar Mass), densities (10(exp 3)-10(exp 4)/cu cm), and temperatures (100-250 K) of the warm atomic gas. The (O III) (52/88) and (S III) (33/19) line flux ratios constrain the range of electron densities and pressures found within the ionized regions. The (O III) and (S III) lines also provide estimates of the effective temperature of the ionizing stars and elemental abundances within the ionized regions of these galactic nuclei. Our measurements imply typical gas pressures of nT approximately 5 x 10(exp 6)/cu cm K and typical upper mass cutoffs of 25-35 Solar Mass. The low-metallicity systems show high (C II)/CO and (O I)/CO flux ratios, 3-5 times the Milky Way value, indicating that they contain a larger fraction of photodissociated gas relative to the molecular material.

Lord, Steven D.↗
DOE OSTI2021

The kinematic properties of Milky Way stellar halo populations

ABSTRACT In the Gaia era, stellar kinematics are extensively used to study Galactic halo stellar populations, to search for halo structures, and to characterize the interface between the halo and hot disc populations. We use distribution function-based models of modern data sets with 6D phase space data to qualitatively describe a variety of kinematic spaces commonly used in the study of the Galactic halo. Furthermore, we quantitatively assess how well each kinematic space can separate radially anisotropic from isotropic halo populations. We find that scaled action space (the ‘action diamond’) is superior to other commonly used kinematic spaces at this task. We present a new, easy to implement selection criterion for members of the radially anisotropic Gaia–Enceladus merger remnant. Assuming a 1:1 ratio of Gaia–Enceladus stars to more isotropic halo, we find our selection achieves a sample purity of 86 per cent in our models with respect to contamination from the more isotropic halo. We compare this criterion to literature criteria, finding that it produces the highest purity in the resulting samples, at the expense of a modest reduction in completeness. We also show that selection biases that underlie nearly all contemporary spectroscopic data sets can noticeably impact the E−Lz distribution of samples in a manner that may be confused for real substructure. We conclude by providing recommendations for how authors should use stellar kinematics in the future to study the Galactic stellar halo.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1998

Giant Molecular Clouds and High-Mass Star Formation in the Milky Way

We are conducting an extensive investigation of high-mass (OB) star formation within the dense cores of giant molecular clouds (GMCS) throughout the first Galactic quadrant of the Milky Way using enhanced resolution Infrared Astronomical Satellite (IRAS) images in combination with high-resolution ground-based observations in millimeter wave molecular transitions and radio continuum. As part of this investigation several resolution enhancement algorithms are applied to the IRAS data, including the HIgh RESolution (HIRES) algorithm developed at the IRAS Processing and Analysis Center (IPAC), as well as others ("pixon" image reconstruction). In addition, as part of a related study, we have completed a large survey of the CO emission in the first Galactic quadrant using the 15-element array detector (QUARRY) with the Five College Radio Astronomy Observatory (FCRAO) 14 m antenna, which provides sampling at an angular resolution of 50", comparable to that attained in the reprocessed IRAS data. Both of these data sets are compared with a sample of ultra-compact (UC) H II regions taken from a high-resolution multi-wavelength (6 and 20 cm) radio survey of the Galactic plane using the NRAO Very Large Array (VLA). Selected regions are observed in 1.3 mm continuum, which has proven to be particularly sensitive to the dust column density. Extensive observations of molecular clouds at high resolution in CO, CS and HCN are combined with the reprocessed IRAS high-resolution images to give a more complete picture of the physical conditions and kinematics of high-mass star forming GMCS. Our goals are to study in detail the morphology, structure, and rate of high-mass star formation within GMCs throughout the Galactic disk from the inner edge of the molecular ring to the outer Galaxy.

Source record↗
DOE OSTI2022

How Many Elements Matter?

Some studies of stars' multielement abundance distributions suggest at least 5–7 significant dimensions, but others show that many elemental abundances can be predicted to high accuracy from [Fe/H] and [Mg/Fe] (or [Fe/H] and age) alone. We show that both propositions can be, and are, simultaneously true. We adopt a machine-learning technique known as normalizing flow to reconstruct the probability distribution of Milky Way disk stars in the space of 15 elemental abundances measured by APOGEE. Conditioning on T eff and $\mathrm{log}\,g$ minimizes the differential systematics. After further conditioning on [Fe/H] and [Mg/Fe], the residual scatter for most abundances is σ [X/H] ≲ 0.02 dex, consistent with APOGEE's reported statistical uncertainties of ~0.01–0.015 dex and intrinsic scatter of 0.01–0.02 dex. Despite the small scatter, residual abundances display clear correlations between elements, which we show are too large to be explained by measurement uncertainties or by the finite sampling noise. We must condition on at least seven elements to reduce the correlations to a level consistent with the observational uncertainties. Our results demonstrate that cross-element correlations are a much more sensitive probe of a hidden structure than dispersion, and they can be measured precisely in a large sample even if the star-by-star measurement noise is comparable to the intrinsic scatter. We conclude that many elements have an independent story to tell, even for the mundane disk stars and elements produced by the core-collapse and Type Ia supernovae. The only way to learn these lessons is to measure the abundances directly, and not merely infer them.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Probing the Galactic halo with RR Lyrae stars – II. The substructures of the Milky Way

ABSTRACT We identify substructures of the Galactic halo using 3003 type ab RR Lyrae stars with six-dimensional position–velocity information from the Sloan Digital Sky Survey, the Large Sky Area Multi-Object Fiber Spectroscopic Telescope and the Gaia Early Data Release 3. Based on this information, we define the separation of any two of the stars in the integrals of motion space and we identify substructures by utilizing the friends-of-friends algorithm. We identify members belonging to several known substructures: the Sagittarius stream, the Gaia-Enceladus-Sausage (GES), the Sequoia and the Helmi streams. In addition to these known substructures, there are three other substructures possibly associated with globular clusters NGC 5272, 6656 and 5024, respectively. Finally, we also find three remaining unknown substructures, one of which has large angular momentum and a mean metallicity −2.13 dex, which may be a new substructure. As for the GES, we find that it accounts for a large part of substructures in the inner halo and the apocentre distance is in the range 10–34 kpc, which suggests that the GES is mainly distributed in the inner halo. The near one-third proportion of the GES and the peak value of 20 kpc of the apocentre distance suggest that the GES could account for the break in the density profile of the Galactic halo at the Galactocentric distance ~20–25 kpc. The similarity when comparing the kinematic properties of the GES with the Hercules–Aquila Cloud and Virgo Overdensity suggests that the three substructures may have similar origins.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Impact of bar resonances in the velocity–space distribution of the solar neighbourhood stars in a self-consistent N -body Galactic disc simulation

ABSTRACT The velocity–space distribution of the solar neighbourhood stars shows complex substructures. Most of the previous studies use static potentials to investigate their origins. Instead we use a self-consistent N-body model of the Milky Way, whose potential is asymmetric and evolves with time. In this paper, we quantitatively evaluate the similarities of the velocity–space distributions in the N-body model and that of the solar neighbourhood, using Kullback–Leibler divergence (KLD). The KLD analysis shows the time evolution and spatial variation of the velocity–space distribution. The KLD fluctuates with time, which indicates the velocity–space distribution at a fixed position is not always similar to that of the solar neighbourhood. Some positions show velocity–space distributions with small KLDs (high similarities) more frequently than others. One of them locates at $(R,\phi)=(8.2\,\,\rm{\mathrm{kpc}}, 30^\circ)$, where R and ϕ are the distance from the galactic centre and the angle with respect to the bar’s major axis, respectively. The detection frequency is higher in the inter-arm regions than in the arm regions. In the velocity maps with small KLDs, we identify the velocity–space substructures, which consist of particles trapped in bar resonances. The bar resonances have significant impact on the stellar velocity–space distribution even though the galactic potential is not static.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

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↗
DOE OSTI2020

Chromium Nucleosynthesis and Silicon–Carbon Shell Mergers in Massive Stars

In this work, we analyze the production of the element Cr in galactic chemical evolution (GCE) models using the NuGrid nucleosynthesis yields set. We show that the unusually large [Cr/Fe] abundance at [Fe/H] ≈ 0 reported by previous studies using those yields and predicted by our Milky Way model originates from the merging of convective Si-burning and C-burning shells in a 20 ${M}_{\odot }$ model at metallicity Z = 0.01, about an hour before the star explodes. This merger mixes the incomplete burning material in the Si shell, including 51 V and 52 Cr, out to the edge of the carbon/oxygen (CO) core. The adopted supernova model ejects the outer 2 ${M}_{\odot }$ of the CO core, which includes a significant fraction of the Cr-rich material. When including this 20 M⊙ model at Z = 0.01 in the yields interpolation scheme of our GCE model for stars between 15 and 25 ${M}_{\odot }$, we overestimate [Cr/Fe] by an order of magnitude at [Fe/H] ≈ 0 relative to observations in the Galactic disk. This raises a number of questions regarding the occurrence of Si–C shell mergers in nature, the accuracy of different simulation approaches, and the impact of such mergers on the presupernova structure and explosion dynamics. According to the conditions in this 1D stellar model, the substantial penetration of C-shell material into the Si shell could launch a convective–reactive global oscillation if a merger does take place. In any case, GCE provides stringent constraints on the outcome of this stellar evolution phase.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Spectral energy distributions of the brightest Palomar-Green quasars at intermediate redshifts

We have combined low-dispersion International Ultraviolet Explorer (IUE) spectra with the optical/near-IR spectrophotometry of Neugebauer et al. (1987) in order to study the spectral energy distributions of seven of the brightest Palomar-Green (PG) quasars at intermediate redshifts (Z(sub em) greater than or equal to 0.9 and less than or equal to 1.5). Some of these PG quasars are barely detectable in long IUE exposures, so we have used the Gaussian Extraction (GEX) technique to maximize the signal-to-noise of the IUE data, and we have co-added all spectra available from the IUE archive for each QSO unless the ultraviolet spectra varied significantly from one exposure to the next. We have corrected the spectral energy distributions for Milky Way reddening using the observed neutral hydrogen column densities on each sight line and the gas-to-dust relation recently derived by Diplas & Savage. Six of the seven quasars are detected down to lambda much less than 700 A in the rest frame, and consequently continuum reddening due to dust in the immediate vicinity of the quasar can have a dramatic effect on the spectral energy distributions. In order to explore the possible importance of intrinsic continuum reddening, we have assembled a heuristic extinction curve which extends to lambda much less than 912 A. Using this heuristic extinction curve, we derive reasonable upper limits on the intrinsic E(B-V) for each quasar. We briefly discuss some of the implications of the derived intrinsic continuum reddening limits. We use geometrically thin accretion disk models to derive the black hole masses and accretion rates implied by the spectral energy distributions. Even if we neglect intrinsic reddening, we find that a large fraction of the quasars require super-Eddington accretion rates (which is not consistent with the thin disk assumption). Comparison of the data in this paper to a large body of data from the literature on the accretion disk M(sub BH) - M dot grid calculated by Wandel & Petrosian reveals that our quasars are among the brightest in the sky at 1450 A, and ostensibly suggests that the fraction of quasars which require super-Eddington accretion rates is much smaller than the fraction that we derive from our data alone. However, intrinsic continuum reddening has been ignored in this comparison, and a small amount of intrinsic reddening will push more of the quasars into the super-Eddington regime. We also plot the recent reverberation monitoring results on NGC 5548 and NGC 3783 on the Wandel & Petrosian grid, and we find that these Seyfert galaxies appear to vary along lines of constant M(sub BH). Continuum flux from two of the quasars in our main sample, PG 1338+416 and PG 1630+377, is detected at lambda(sub rest) less than 584 A. These quasars can in principle be used for the He I Gunn-Peterson test, but the S/N of IUE spectra of individual objects is usually too low to place interesting limits on the Gunn-Peterson optical depth. In order to improve the S/N, we have formed a composite spectrum from the spectra of five quasars detected with IUE at lambda(sub rest) less than 584 A, and we have used this composite spectrum to place a tighter limit on tau(sub GP, He I). We briefly discuss intermediate-redshift Lyman limit systems (Z(sub LL) greater than or equal to 0.5 and less than or equal to 1.6) detected in the IUE spectra of five quasars, including lower limits on N(H I) in each Lyman limit system.

Tripp, Todd M.↗
DOE OSTI2022

Reionization time of the Local Group and Local-Group-like halo pairs

ABSTRACT Patchy cosmic reionization resulted in the ionizing UV background asynchronous rise across the Universe. The latter might have left imprints visible in present-day observations. Several numerical simulation-based studies show correlations between the reionization time and overdensities and object masses today. To remove the mass from the study, as it may not be the sole important parameter, this paper focuses solely on the properties of paired haloes within the same mass range as the Milky Way. For this purpose, it uses CoDaII, a fully coupled radiation hydrodynamics reionization simulation of the local Universe. This simulation holds a halo pair representing the Local Group, in addition to other pairs, sharing similar mass, mass ratio, distance separation, and isolation criteria but in other environments, alongside isolated haloes within the same mass range. Investigations of the paired halo reionization histories reveal a wide diversity although always inside-out, given our reionization model. Within this model, haloes in a close pair tend to be reionized at the same time but being in a pair does not bring to an earlier time their mean reionization. The only significant trend is found between the total energy at z = 0 of the pairs and their mean reionization time: Pairs with the smallest total energy (bound) are reionized up to 50 Myr earlier than others (unbound). Above all, this study reveals the variety of reionization histories undergone by halo pairs similar to the Local Group, that of the Local Group being far from an average one. In our model, its reionization time is ∼625 Myr against 660 ± 4 Myr (z ∼ 8.25 against 7.87 ± 0.02) on average.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Diverse Chemo-Dynamical Properties of Nitrogen-Rich Stars Identified from Low-Resolution Spectra

The second generation of stars in the globular clusters (GCs) of the Milky Way (MW) exhibit unusually high N, Na, or Al, compared to typical Galactic halo stars at similar metallicities. The halo field stars enhanced with such elements are believed to have originated in disrupted GCs or escaped from existing GCs. We identify such stars in the metallicity range –3.0 < [Fe/H] < 0.0 from a sample of ~ 36,800 giant stars observed in the Sloan Digital Sky Survey and Large Sky Area Multi-Object Fiber Spectroscopic Telescope survey, and present their dynamical properties. The N-rich population (NRP) and N-normal population (NNP) among our giant sample do not exhibit similarities in either in their metallicity distribution function (MDF) or dynamical properties. We find that, even though the MDF of the NRP looks similar to that of the MW’s GCs in the range of [Fe/H] < –1.0, our analysis of the dynamical properties does not indicate similarities between them in the same metallicity range, implying that the escaped members from existing GCs may account for a small fraction of our N-rich stars, or the orbits of the present GCs have been altered by the dynamical friction of the MW. We also find a significant increase in the fraction of N-rich stars in the halo field in the very metal-poor (VMP; [Fe/H] < –2.0) regime, comprising up to ~ 20% of the fraction of the N-rich stars below [Fe/H] = –2.5, hinting that partially or fully destroyed VMP GCs may have in some degree contributed to the Galactic halo. A more detailed dynamical analysis of the NRP reveals that our sample of N-rich stars do not share a single common origin. Although a substantial fraction of the N-rich stars seem to originate from the GCs formed in situ, more than 60% of them are not associated with those of typical Galactic populations, but probably have extragalactic origins associated with Gaia Sausage/Enceladus, Sequoia, and Sagittarius dwarf galaxies, as well as with presently unrecognized progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Robust clustering of the local Milky Way stellar kinematic substructures with Gaia eDR3

Understanding local stellar kinematic substructures in the solar neighbourhood helps build a complete picture of the formation of the Milky Way, as well as an empirical phase space distribution of dark matter that would inform detection experiments. We apply the clustering algorithm HDBSCAN on the Gaia early third data release to identify a list of stable clusters in velocity space and action-angle space by taking into account the measurement uncertainties and studying the stability of the clustering results. We find 1405 (497) stars in 23 (6) robust clusters in velocity space (action-angle space) that are consistently not associated with noise. We discuss the kinematic properties of these structures and study whether many of the small clusters belong to a similar larger cluster based on their chemical abundances. They are attributed to the known structures: the Gaia Sausage-Enceladus, the Helmi Stream, and globular cluster NGC 3201 are found in both spaces, while NGC 104 and the thick disc (Sequoia) are identified in velocity space (action-angle space). Although we do not identify any new structures, we find that the HDBSCAN member selection of already known structures is unstable to input kinematics of the stars when resampled within their uncertainties. We therefore present the stable subset of local kinematic structures, which are consistently identified by the clustering algorithm, and emphasize the need to take into account error propagation during both the manual and automated identification of stellar structures, both for existing ones as well as future discoveries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

S 5 : Probing the Milky Way and Magellanic Clouds potentials with the 6D map of the Orphan–Chenab stream

We present a 6D map of the Orphan–Chenab (OC) stream by combining the data from Southern Stellar Stream Spectroscopic Survey (S 5 ) and Gaia. We reconstruct the proper motion, radial velocity, distance, on-sky track, and stellar density along the stream with spline models. The stream has a total luminosity of M V = -8.2 and metallicity of [Fe/H] = -1.9, similar to classical Milky Way (MW) satellites like Draco. The stream shows drastic changes in its physical width varying from 200 pc to 1kpc, but a constant line-of-sight velocity dispersion of 5 km s -1 . Despite the large apparent variation in the stellar number density along the stream, the flow rate of stars along the stream is remarkably constant. We model the 6D stream track by a Lagrange-point stripping method with a flexible MW potential in the presence of a moving extended Large Magellanic Cloud (LMC). This allows us to constrain the mass profile of the MW within the distance range 15.6 < r < 55∙5 kpc, with the best measured enclosed mass of (2.85 ± 0.1) X 10 11 M ⊙ within 32.4 kpc. Our stream measurements are highly sensitive to the LMC mass profile with the most precise measurement of its enclosed mass made at 32.8 kpc, 7.02 ± 0.9) X 10 11 M ⊙ . We also detect that the LMC dark matter halo extends to at least 53 kpc. The fitting of the OC stream allows us to constrain the past LMC trajectory and the degree of dynamical friction it experienced. We demonstrate that the stars in the OC stream show large energy and angular momentum spreads caused by LMC perturbation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Photometric calibration in u -band using blue halo stars

ABSTRACT We develop a method to calibrate u-band photometry based on the observed colour of blue Galactic halo stars. The Galactic halo stars belong to an old stellar population of the Milky Way and have relatively low metallicity. The ‘blue tip’ of the halo population – the main sequence turn-off (MSTO) stars – is known to have a relatively uniform intrinsic edge u-g colour with only slow spatial variation. In SDSS data, the observed variation is correlated with Galactic Latitude, which we attribute to contamination by higher metallicity disc stars and fit with an empirical curve. This curve can then be used to calibrate u-band imaging if g-band imaging of matching depth is available. Our approach can be applied to single-field observations at |b| > 30°, and removes the need for standard star observations or overlap with calibrated u-band imaging. We include in our method the calibration of g-band data with ATLAS-Refcat2. We test our approach on stars in KiDS DR 4, ATLAS DR 4, and DECam imaging from the NOIRLab Source Catalog (NSC DR2), and compare our calibration with SDSS. For this process, we use synthetic magnitudes to derive the colour equations between these data sets, in order to improve zero-point accuracy. We find an improvement for all data sets, reaching a zero-point precision of 0.016 mag for KiDS (compared to the original 0.033 mag), 0.020 mag for ATLAS (originally 0.027 mag), and 0.016 mag for DECam (originally 0.041 mag). Thus, this method alone reaches the goal of 0.02 mag photometric precision in u-band for the Rubin Observatory’s Legacy Survey of Space and Time (LSST).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Decoding the age–chemical structure of the Milky Way disc: an application of copulas and elicitable maps

In the Milky Way, the distribution of stars in the [α/Fe] versus [Fe/H] and [Fe/H] versus age planes holds essential information about the history of star formation, accretion, and dynamical evolution of the Galactic disc. We investigate these planes by applying novel statistical methods called copulas and elicitable maps to the ages and abundances of red giants in the Apache Point Observatory Galactic Evolution Experiment survey. We find that the high- and low-α disc stars have a clean separation in copula space and use this to provide an automated separation of the α sequences using a purely statistical approach. This separation reveals that the high-α disc ends at the same [α/Fe] and age at high [Fe/H] as the low-[Fe/H] start of the low-α disc, thus supporting a sequential formation scenario for the high- and low-α discs. We then combine copulas with elicitable maps to precisely obtain the correlation between stellar age τ and metallicity [Fe/H] conditional on Galactocentric radius R and height z in the range 0 < R < 20 kpc and |z| < 2 kpc. The resulting trends in the age–metallicity correlation with radius, height, and [α/Fe] demonstrate a ≈0 correlation wherever kinematically cold orbits dominate, while the naively expected negative correlation is present where kinematically hot orbits dominate. This is consistent with the effects of spiral-driven radial migration, which must be strong enough to completely flatten the age–metallicity structure of the low-α disc.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

The ISM in the M82 starburst

We have observed (O 1) (63 microns) and (Si 2) (35 microns) in the central 700 pc of the starburst galaxy M82. The luminosities in these transitions are 7.1 x 10(exp 7) solar luminosity and 6.2 x 10(exp 7) solar luminosity, respectively, which are each approx. 0.15% of the bolometric luminosity from this region. The ratios of (O 1) line luminosity to (O 3), (Si 2) (35 microns), and to bolometric luminosities in M82 are similar to those in M42, M17, and Sgr A. These similarities, and the association of the bulk of the (O 1) and (Si 2) emission with the ionized emission, suggest that the dominant emission mechanism for (O 1) and (Si 2) in M82 is the same as in these Galactic regions, namely warm gas photodissociated by UV flux from the OB stars responsible for the nearby H 2 regions. We argue that shock or x ray heated gas or H 2 plasma is a minor contributor to the intensities of these fine structure lines. Both the (O 1) (53 microns) and the (Si 2) (35 microns) spectrum show an asymmetric line profile indistinguishable in shape from those of the (O 3) (52 and 88 microns) and (N 3) (57 microns) lines and similar to that of the more extended (C 2) 158 micron line measured previously in M82. We detect two distinct velocity components, which we attribute to emission from two regions at either end of the central bar, where the bar connects to an orbiting torus of neutral gas seen in H 1 and CO J = 1-0. We model separately the two velocity components and derive the physical conditions in these two regions. The clouds in these regions are small, R approx. 1-2 pc, have warm neutral gas surfaces, T approx. 200 K, and are concentrated with volume filling factors of approx. 0.02 and area filling factors of 1-5. The entire central region (R approx. 700 pc) is characterized by a large number, approx. 5 x 10(exp 4), of 2 x 10(exp 3) solar mass clouds with surface densities of approx. 3 x 10(exp 4) cm(exp -3), illuminated by FUV fluxes 10(exp 4) times the average local interstellar value for the Milky Way. These clouds reside in the harsh conditions of a starburst nucleus, with photoevaporation times of 10(exp 6) yr, and collision timescales only about an order of magnitude longer.

Lord, Steven↗
DOE OSTI2022

Substructure in the stellar halo near the Sun: II. Characterisation of independent structures

In an accompanying paper, we present a data-driven method for clustering in ‘integrals of motion’ space and apply it to a large sample of nearby halo stars with 6D phase-space information. The algorithm identified a large number of clusters, many of which could tentatively be merged into larger groups. The goal here is to establish the reality of the clusters and groups through a combined study of their stellar populations (average age, metallicity, and chemical and dynamical properties) to gain more insights into the accretion history of the Milky Way. To this end, we developed a procedure that quantifies the similarity of clusters based on the Kolmogorov–Smirnov test using their metallicity distribution functions, and an isochrone fitting method to determine their average age, which is also used to compare the distribution of stars in the colour–absolute magnitude diagram. Also taking into consideration how the clusters are distributed in integrals of motion space allows us to group clusters into substructures and to compare substructures with one another. We find that the 67 clusters identified by our algorithm can be merged into 12 extended substructures and 8 small clusters that remain as such. The large substructures include the previously known Gaia-Enceladus, Helmi streams, Sequoia, and Thamnos 1 and 2. We identify a few over-densities that can be associated with the hot thick disc and host a small metal-poor population. Especially notable is the largest (by number of member stars) substructure in our sample which, although peaking at the metallicity characteristic of the thick disc, has a very well populated metal-poor component, and dynamics intermediate between the hot thick disc and the halo. We also identify additional debris in the region occupied by Sequoia with clearly distinct kinematics, likely remnants of three different accretion events with progenitors of similar masses. Although only a small subset of the stars in our sample have chemical abundance information, we are able to identify different trends of [Mg/Fe] versus [Fe/H] for the various substructures, confirming our dissection of the nearby halo. We find that at least 20% of the halo near the Sun is associated to substructures. When comparing their global properties, we note that those substructures on retrograde orbits are not only more metal-poor on average but are also older. We provide a table summarising the properties of the substructures, as well as a membership list that can be used for follow-up chemical abundance studies for example.

79 ASTRONOMY AND ASTROPHYSICS↗