Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Milky Way dynamics”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Substructure at High Speed. II. The Local Escape Velocity and Milky Way Mass with Gaia eDR3

Abstract Measuring the escape velocity of the Milky Way is critical in obtaining the mass of the Milky Way, understanding the dark matter velocity distribution, and building the dark matter density profile. In Necib & Lin, we introduced a strategy to robustly measure the escape velocity. Our approach takes into account the presence of kinematic substructures by modeling the tail of the stellar distribution with multiple components, including the stellar halo and the debris flow called the Gaia Sausage (Enceladus). In doing so, we can test the robustness of the escape velocity measurement for different definitions of the “tail” of the velocity distribution and the consistency of the data with different underlying models. In this paper, we apply this method to the Gaia eDR3 data release and find that a model with two components is preferred, although results from a single-component fit are also consistent. Based on a fit to retrograde data with two bound components to account for the relaxed halo and the Gaia Sausage, we find the escape velocity of the Milky Way at the solar position to be v esc = 445 − 8 + 25 km s −1 . A fit with a single component to the same data gives v esc = 472 − 12 + 17 km s −1 . Assuming a Navarro−Frenck−White dark matter profile, we find a Milky Way concentration of c 200 = 19 − 7 + 11 and a mass of M 200 = 4.6 − 0.8 + 1.5 × 10 11 M ⊙ , which is considerably lighter than previous measurements.

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

Effects of Dark Matter Self-interactions on Sagittarius and its Stream

This work explores how assumptions regarding the particle-physics nature of dark matter can alter the evolution of the Sagittarius (Sgr) dwarf spheroidal galaxy and its expansive stellar stream. We run a large suite of N-body simulations to model the infall of a spherically symmetric Sgr-like dwarf, exploring how the presence of dark matter self-interactions impacts its evolution. For a scattering cross section of σ/m χ = 30 cm 2 g −1 (at orbital velocity scales), these interactions result in significantly less stellar mass and little to no dark matter bound to the progenitor at the present day. To isolate the cause of this mass loss, we introduce a novel technique for controlling which pairs of dark matter simulation particles can interact. This enables us to identify ram-pressure evaporation—the scattering of satellite and host dark matter particles—as the primary source of the enhanced mass loss. The rapid disintegration of the Sgr progenitor when self-interactions are allowed alters some key properties of the resulting stellar stream, most dramatically suppressing the presence of a “spur” on the apocenter of the trailing stream arm that correlates with the mass of the satellite at last pericenter. We demonstrate how the effects on the Sgr system scale with the particular choice of self-interaction cross section, which affects the degree of ram-pressure evaporation. These findings generalize beyond the Sgr system, underscoring that dwarf stellar streams and dwarf galaxies with close passages may serve as sensitive probes for dark matter self-interactions.

Milky Way dark matter halo↗
DOE OSTI2023

Dynamically constraining the length of the Milky way bar

ABSTRACT We present a novel method for constraining the length of the Galactic bar using 6D phase-space information to directly integrate orbits. We define a pseudo-length for the Galactic bar, named RFreq, based on the maximal extent of trapped bar orbits. We find the RFreq measured from orbits is consistent with the RFreq of the assumed potential only when the length of the bar and pattern speed of said potential is similar to the model from which the initial phase-space coordinates of the orbits are derived. Therefore, one can measure the model’s or the Milky Way’s bar length from 6D phase-space coordinates by determining which assumed potential leads to a self-consistent measured RFreq. When we apply this method to ≈210 000 stars in APOGEE DR17 and Gaia eDR3 data, we find a consistent result only for potential models with a dynamical bar length of ≈3.5 kpc. We find the Milky Way’s trapped bar orbits extend out to only ≈3.5 kpc, but there is also an overdensity of stars at the end of the bar out to 4.8 kpc which could be related to an attached spiral arm. We also find that the measured orbital structure of the bar is strongly dependent on the properties of the assumed potential.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The Milky Way, coming into focus: Precision astrometry probes its evolution and its dark matter

The growing trove of precision astrometric observations from the Gaia space telescope and other surveys is revealing the structure and dynamics of the Milky Way in ever more exquisite detail. We summarize the current status of our understanding of the structure and the characteristics of the Milky Way, and we review the emerging picture: the Milky Way is evolving through interactions with the massive satellite galaxies that stud its volume, with evidence pointing to a cataclysmic past. It is also woven with stellar streams, and observations of streams, satellites, and field stars offer new constraints on its dark matter, both on its spatial distribution and its fundamental nature. The recent years have brought much focus to the study of dwarf galaxies found within our Galaxy’s halo and their internal matter distributions. In this review, we focus on the predictions of the cold dark matter paradigm at small mass scales through precision astrometric measurements, and we summarize the modern consensus on the extent to which small-scale probes are consistent with this paradigm. We note the discovery prospects of these studies, and also how they intertwine with probes of the dynamics and evolution of the Milky Way in various and distinct ways.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

ESCARGOT: Mapping Vertical Phase Spiral Characteristics Throughout the Real and Simulated Milky Way

Abstract The recent discovery of a spiral pattern in the vertical kinematic structure in the solar neighborhood provides a prime opportunity to study nonequilibrium dynamics in the Milky Way from local stellar kinematics. Furthermore, results from simulations indicate that even in a limited volume, differences in stellar orbital histories allow us to trace variations in the initial perturbation across large regions of the disk. We present ESCARGOT , a novel algorithm for studying these variations in both simulated and observed data sets. ESCARGOT automatically extracts key quantities from the structure of a given phase spiral, including the time since perturbation and the perturbation mode. We test ESCARGOT on simulated data and show that it is capable of accurately recovering information about the time since the perturbation occurred as well as subtle differences in phase spiral morphology due to stellar locations in the disk at the time of perturbation. We apply ESCARGOT to kinematic data from data release 3 of the Gaia mission in bins of guiding radius. We show that similar structural differences in morphology occur in the Gaia phase spirals as a function of stellar orbital history. These results indicate that the phase spirals are the product of a complex dynamical response in the disk with large-scale coupling between different regions of phase space.

Darragh-Ford, Elise (ORCID:0000000288005652)↗
DOE OSTI2023

A Blueprint for the Milky Way’s Stellar Populations. IV. A String of Pearls—the Galactic Starburst Sequence

We continue our series of papers on phase-space distributions of stars in the Milky Way based on photometrically derived metallicities and Gaia astrometry, with a focus on the halo-disk interface in the local volume. To exploit various photometric databases, we develop a method of empirically calibrating synthetic stellar spectra based on a comparison with observations of stellar sequences and individual stars in the Sloan Digital Sky Survey, the SkyMapper Sky Survey, and the Pan-STARRS1 surveys, overcoming band-specific corrections employed in our previous work. In addition, photometric zero-point corrections are derived to provide an internally consistent photometric system with a spatially uniform metallicity zero-point. Using our phase-space diagrams, we find a remarkably narrow sequence in the rotational velocity (v $\phi$ ) versus metallicity ([Fe/H]) space for a sample of high proper-motion stars (>25 mas yr -1 ), which runs along Gaia Sausage/Enceladus (GSE) and the Splash substructures and is linked to the disk, spanning nearly 2 dex in [Fe/H]. Notably, a rapid increase of v $\phi$ from a nearly zero net rotation to ~180 km s -1 in a narrow metallicity interval (-0.6 ≲ [Fe/H] ≲ -0.4) suggests that some of these stars emerged quickly on a short gas-depletion timescale. Through measurements of a scale height and length, we argue that these stars are distinct from those heated dynamically by mergers. This chain of high proper-motion stars provides additional support for recent discoveries suggesting that a starburst took place when the young Milky Way encountered the gas-rich GSE progenitor, which eventually led to the settling of metal-enriched gas onto the disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

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

ABSTRACT 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$\lt 6\times 10^{-4} \, \mathrm{M}_{\odot }$ yr−1, and no X-rays are detected down to a sensitivity of 1038 erg s−1. This rules out a hot corona with a $M \gt 10^8 \, \mathrm{M}_{\odot }$ 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.

Lee, Chris H.↗
NASA NTRS2010

The Nuclear Spectroscopic Telescope Array (NuSTAR)

The Nuclear Spectroscopic Telescope Array (NuSTAR) is a NASA Small Explorer mission that will carry the first focusing hard X-ray (5 - 80 keV) telescope to orbit. NuSTAR will offer a factor 50 - 100 sensitivity improvement compared to previous collimated or coded mask imagers that have operated in this energy band. In addition, NuSTAR provides sub-arcminute imaging with good spectral resolution over a 12-arcminute field of view. After launch, NuSTAR will carry out a two-year primary science mission that focuses on four key programs: studying the evolution of massive black holes through surveys carried out in fields with excellent multiwavelength coverage, understanding the population of compact objects and the nature of the massive black hole in the center of the Milky Way, constraining explosion dynamics and nucleosynthesis in supernovae, and probing the nature of particle acceleration in relativistic jets in active galactic nuclei. A number of additional observations will be included in the primary mission, and a. guest observer program will be proposed for an extended mission to expand the range of scientific targets. The payload consists of two co-aligned depth-graded multilayer coated grazing incidence optics focused onto solid state CdZnTe pixel detectors. To be launched in early 2012 on a Pegasus rocket into a low-inclination Earth orbit. NuSTAR largely avoids SAA passages, and will therefore have low and stable detector backgrounds. The telescope achieves a 10.15-meter focal length through on-orbit deployment of all mast. An aspect and alignment metrology system enable reconstruction of the absolute aspect and variations in the telescope alignment resulting from mast flexure during ground data processing. Data will be publicly available at GSFC's High Energy Astrophysics Science Archive Research Center (HEASARC) following validation at the science operations center located at Caltech.

Harrison, Fiona A.↗
NASA NTRS2020

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

DESI Early Data Release Milky Way Survey value-added catalogue

We present the stellar value-added catalogue based on the Dark Energy Spectroscopic Instrument (DESI) Early Data Release. The catalogue contains radial velocity and stellar parameter measurements for $\simeq$ 400 000 unique stars observed during commissioning and survey validation by DESI. These observations were made under conditions similar to the Milky Way Survey (MWS) currently carried out by DESI but also include multiple specially targeted fields, such as those containing well-studied dwarf galaxies and stellar streams. The majority of observed stars have $16\lt r\lt 20$ with a median signal-to-noise ratio in the spectra of $\sim$ 20. In the paper, we describe the structure of the catalogue, give an overview of different target classes observed, as well as provide recipes for selecting clean stellar samples. We validate the catalogue using external high-resolution measurements and show that radial velocities, surface gravities, and iron abundances determined by DESI are accurate to 1 km s −1 , 0.3 dex, and $\sim$ 0.15 dex respectively. We also demonstrate possible uses of the catalogue for chemo-dynamical studies of the Milky Way stellar halo and Draco dwarf spheroidal. The value-added catalogue described in this paper is the very first DESI MWS catalogue. The next DESI data release, expected in less than a year, will add the data from the first year of DESI survey operations and will contain approximately 4 million stars, along with significant processing improvements.

Galaxy: general↗
DOE OSTI2021

Uncovering fossils of the distant Milky Way with UNIONS: NGC 5466 and its stellar stream

ABSTRACT We examine the spatial clustering of blue horizontal branch (BHB) stars from the u-band of the Canada–France Imaging Survey (CFIS, a component of the Ultraviolet Near-Infrared Optical Northern Survey, or UNIONS). All major groupings of stars are associated with previously known satellites, and among these is NGC 5466, a distant (16 kpc) globular cluster. NGC 5466 reportedly possesses a long stellar stream, although no individual members of the stream have previously been identified. Using both BHBs and more numerous red giant branch stars cross-matched to Gaia Data Release 2, we identify extended tidal tails from NGC 5466 that are both spatially and kinematically coherent. Interestingly, we find that this stream does not follow the same path as the previous detection at large distances from the cluster. We trace the stream across 31° of sky and show that it exhibits a very strong distance gradient in the range 10 < Rhelio < 30 kpc. We compare our observations to simple dynamical models of the stream and find that they are able to broadly reproduce the overall path and kinematics. The fact that NGC 5466 is so distant, traces a wide range of Galactic distances, has an identified progenitor, and appears to have recently had an interaction with the Galaxy’s disc makes it a unique test-case for dynamical modelling of the Milky Way.

Jensen, Jaclyn↗
DOE OSTI2022

Combined $\mathrm{APOGEE}$-$\mathrm{GALAH}$ stellar catalogues using the Cannon

APOGEE and GALAH are two high resolution multi-object spectroscopic surveys that provide fundamental stellar parameters and multiple elemental abundance estimates for about half a million stars in the Milky Way. Both surveys observe in different wavelength regimes and use different data reduction pipelines leading to significant offsets and trends in stellar parameters and abundances for the common stars observed in both surveys. Such systematic differences/offsets in stellar parameters and abundances make it difficult to effectively utilize them to investigate Galactic abundance trends in spite of the unique advantage provided by their complementary sky coverage and different Milky Way components they observe. Hence, we use the Cannon data-driven method selecting a training set of 4418 common stars observed by both surveys. This enables the construction of two catalogues, one with the APOGEE-scaled and the other with the GALAH-scaled stellar parameters. Using repeat observations in APOGEE and GALAH, we find high precision in metallicity (~0.02–0.4 dex) and alpha abundances (~0.02–0.03 dex) for spectra with good signal-to-noise ratio (SNR >; 80 for APOGEE and SNR > 40 for GALAH). We use open and globular clusters to validate our parameter estimates and find small scatter in metallicity (0.06 dex) and alpha abundances (0.03 dex) in APOGEE-scaled case. The final catalogues have been cross-matched with the Gaia EDR3 catalogue to enable their use to carry out detailed chemo-dynamic studies of the Milky Way from perspectives of APOGEE and GALAH.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A search for stellar structures around nine outer halo globular clusters in the Milky Way

We use deep imaging from the Dark Energy Camera to explore the peripheral regions of nine globular clusters in the outer halo of the Milky Way. Apart from Whiting 1 and NGC 7492, which are projected against the Sagittarius stream, we see no evidence for adjacent stellar populations to indicate any of these clusters is associated with coherent tidal debris from a destroyed host dwarf. We also find no evidence for tidal tails around any of the clusters in our sample; however, both NGC 1904 and 6981 appear to possess outer envelopes. Here, motivated by a slew of recent Gaia -based discoveries, we compile a sample of clusters with robust detections of extra-tidal structure, and search for correlations with orbital properties. While we observe that clusters with tidal tails are typically on moderately or very eccentric orbits that are highly inclined to the Galactic plane and often retrograde, these are neither necessary nor sufficient conditions for the formation of extra-tidal structure. That many objects with tidal tails appear to be accreted leads us to speculate that this lack of consistency may stem from the inhomogeneous dynamical history of the Milky Way globular cluster system. Finally, we note that clusters with prominent stellar envelopes detected in ground-based imaging (such as NGC 1851 and 7089) are now all known from Gaia to possess long tidal tails – experimental confirmation that the presence of an extended envelope is indicative of tidal erosion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Functional Data Analysis for Extracting the Intrinsic Dimensionality of Spectra: Application to Chemical Homogeneity in the Open Cluster M67

High-resolution spectroscopic surveys of the Milky Way have entered the Big Data regime and have opened avenues for solving outstanding questions in Galactic archeology. However, exploiting their full potential is limited by complex systematics, whose characterization has not received much attention in modern spectroscopic analyses. In this work, we present a novel method to disentangle the component of spectral data space intrinsic to the stars from that due to systematics. Using functional principal component analysis on a sample of 18,933 giant spectra from APOGEE, we find that the intrinsic structure above the level of observational uncertainties requires ≈10 functional principal components (FPCs). Our FPCs can reduce the dimensionality of spectra, remove systematics, and impute masked wavelengths, thereby enabling accurate studies of stellar populations. To demonstrate the applicability of our FPCs, we use them to infer stellar parameters and abundances of 28 giants in the open cluster M67. We employ Sequential Neural Likelihood, a simulation-based Bayesian inference method that learns likelihood functions using neural density estimators, to incorporate non-Gaussian effects in spectral likelihoods. By hierarchically combining the inferred abundances, we limit the spread of the following elements in M67: Fe ≲ 0.02 dex; C ≲ 0.03 dex; O, Mg, Si, Ni ≲ 0.04 dex; Ca ≲ 0.05 dex; N, Al ≲ 0.07 dex (at 68% confidence). Our constraints suggest a lack of self-pollution by core-collapse supernovae in M67, which has promising implications for the future of chemical tagging to understand the star formation history and dynamical evolution of the Milky Way.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Dynamics of stellar disc tilting from satellite mergers

The Milky Way’s stellar disc can tilt in response to torques that result from infalling satellite galaxies and their associated tidal debris. In this work, we explore the dynamics of disc tilting by running N-body simulations of mergers in an isolated, isotropic Milky Way-like host galaxy, varying over satellite virial mass, initial position, and orbit. We develop and validate a first-principles understanding of the dynamics that govern how the host galaxy’s stellar disc responds to the satellite’s dark matter (DM) debris. We find that the degree of disc tilting can be large for cosmologically motivated merger histories. In particular, our results suggest that the Galactic disc may still be tilting in response to Gaia-Sausage-Enceladus, one of the most significant recent mergers in the Milky Way’s history. In conclusion, these findings have implications for terrestrial direct detection experiments as disc tilting changes the relative location of the Sun with respect to DM substructure left behind by a merging galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Slant, Fan, and Narrow: The Response of Stellar Streams to a Tilting Galactic Disk

Stellar streams are sensitive tracers of the gravitational potential, which is typically assumed to be static in the inner Galaxy. However, massive mergers like Gaia-Sausage-Enceladus can impart torques on the stellar disk of the Milky Way that result in the disk tilting at rates of up to 10°–20° Gyr –1 . Here, we demonstrate the effects of disk tilting on the morphology and kinematics of stellar streams. Through a series of numerical experiments, we find that streams with nearby apocenters (r apo ≲ 20 kpc) are sensitive to disk tilting, with the primary effect being changes to the stream's on-sky track and width. Interestingly, disk tilting can produce both more diffuse streams and more narrow streams, depending on the orbital inclination of the progenitor and the direction in which the disk is tilting. Our model of Pal 5's tidal tails for a tilting rate of 15° Gyr –1 is in excellent agreement with the observed stream's track and width, and reproduces the extreme narrowing of the trailing tail. We also find that failure to account for a tilting disk can bias constraints on shape parameters of the Milky Way's local dark matter distribution at the level of 5%–10%, with the direction of the bias changing for different streams. Disk tilting could therefore explain discrepancies in the Milky Way's dark matter halo shape inferred using different streams.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Cosmic-Ray Feedback on Bistable Interstellar Medium Turbulence

Abstract While cosmic rays ( E ≳ 1 GeV) are well coupled to a galaxy’s interstellar medium (ISM) at scales of L > 100 pc, adjusting stratification and driving outflows, their impact on small scales is less clear. Based on calculations of the cosmic-ray diffusion coefficient from observations of the grammage in the Milky Way, cosmic rays have little time to dynamically impact the ISM on those small scales. Using numerical simulations, we explore how more complex cosmic-ray transport could allow cosmic rays to couple to the ISM on small scales. We create a two-zone model of cosmic-ray transport, with the cosmic-ray diffusion coefficient set at the estimated Milky Way value in cold gas but smaller in warm gas. We compare this model to simulations with a constant diffusion coefficient. Quicker diffusion through cold gas allows more cold gas to form compared to a simulation with a constant, small diffusion coefficient. However, slower diffusion in warm gas allows cosmic rays to take energy from the turbulent cascade anisotropically. This cosmic-ray energization comes at the expense of turbulent energy which would otherwise be lost during radiative cooling. Finally, we show our two-zone model is capable of matching observational estimates of the grammage for some transport paths through the simulation.

79 ASTRONOMY AND ASTROPHYSICS↗