Search

Metals in the diffuse, ionized gas at the boundary between the Milky Way's interstellar medium (ISM) and circumgalactic medium (CGM), known as the disk-halo interface (DHI), are valuable tracers of the feedback processes that drive the Galactic fountain. However, metallicity measurements in this region are challenging due to obscuration by the Milky Way ISM and uncertain ionization corrections that affect the total hydrogen column density. In this work, we constrain the ionization corrections to neutral hydrogen column densities using precisely measured electron column densities from the dispersion measure of pulsars that lie in the same globular clusters as UV-bright targets with high-resolution absorption spectroscopy. We address the blending of absorption lines with the ISM by jointly fitting Voigt profiles to all absorption components. We present our metallicity estimates for the DHI of the Milky Way based on detailed photoionization modeling to the absorption from ionized metal lines and ionization-corrected total hydrogen columns. Generally, the gas clouds show a large scatter in metallicity, ranging between $0.04-3.2\ Z_{\odot}$, implying that the DHI consists of a mixture of gaseous structures having multiple origins. We estimate the inflow and outflow timescales of the DHI ionized clouds to be $6 - 35$ Myr. We report the detection of an infalling cloud with super-solar metallicity that suggests a Galactic fountain mechanism, whereas at least one low-metallicity outflowing cloud ($Z < 0.1\ Z_{\odot}$) poses a challenge for Galactic fountain and feedback models., Comment: 28 pages, 15 figures, accepted for publication in ApJ

We present an analysis of \textit{HST}/COS/G160M observations of CIV in the inner circumgalactic medium (CGM) of a novel sample of eight z$\sim$0, L$\approx$L$^{\star}$ galaxies, paired with UV-bright QSOs at impact parameters ($R_\mathrm{proj}$) between 25-130 kpc. The galaxies in this stellar-mass-controlled sample (log$_{10}$M$_{\star}$/M$_{\odot}$ $\sim$ 10.2-10.9 M$_{\odot}$) host super-massive black holes (SMBHs) with dynamically-measured masses spanning log$_{10}$M$_\mathrm{BH}$/M$_{\odot}$ $\sim$ 6.8-8.4; this allows us to compare our results with models of galaxy formation where the integrated feedback history from the SMBH alters the CGM over long timescales. We find that the \ion{C}{IV} column density measurements (N$_{\rm C IV}$) (average log$_{10}$N$_{\rm C IV, CH}$ = 13.94$\pm$0.09 cm$^{-2}$) are largely consistent with existing measurements from other surveys of N$_{\rm C IV}$ in the CGM (average log$_{10}$N$_{\rm C IV, Lit}$ = 13.90$\pm$0.08 cm$^{-2}$), but do not show obvious variation as a function of the SMBH mass. In contrast, specific star-formation rate (sSFR) is highly correlated with the ionized content of the CGM. We find a large spread in sSFR for galaxies with log$_{10}$M$_\mathrm{BH}$/M$_{\odot}$ $>$ 7.0, where the CGM \ion{C}{IV} content shows clear dependence on galaxy sSFR but not M$_\mathrm{BH}$. Our results do not indicate an obvious causal link between CGM CIV and the mass of the galaxy's SMBH; however through comparisons to the EAGLE, Romulus25, $\&$ IllustrisTNG simulations, we find that our sample is likely too small to constrain such causality., Comment: Accepted to ApJ, 28 pages, 7 figures

We present a sample of 305 QSO candidates having $|b| < 30^{\circ}$, the majority with GALEX magnitudes NUV < 18.75. To generate this sample, we apply UV-IR color selection criteria to photometric data from the Ultraviolet GAlactic Plane Survey (UVGAPS) as part of GALEX-CAUSE, the Million Quasars Catalog, Gaia DR2, and Pan-STARRS DR1. 165 of these 305 candidate UV-bright AGN (54%) have published spectroscopic redshifts from 45 different surveys, confirming them as AGN. We further obtained low-dispersion, optical, longslit spectra with the APO 3.5-m, MDM 2.4-m, and MDM 1.3-m telescopes for 84 of the candidates, and confirm 86% (N = 72) as AGN, generally with z < 0.6. These sources fill a gap in the Galactic latitude coverage of the available samples of known UV-bright QSO background probes. Along with a description of the confirmed QSO properties, we provide the fully-reduced, flux and wavelength-calibrated spectra of 84 low-latitude QSOs through the Mikulski Archive for Space Telescopes. Future HST/COS spectroscopy of these low-Galactic-latitude QSOs has the potential to transform our view of the Milky Way and Local Group circumgalactic medium., Comment: 18 pages, 4 figures, 2 tables. Submitted to ApJS, comments welcome. If MAST DOI linked in text is not yet published, please email first author for full table of Plane QSOs

The cycling of metals between interstellar gas and dust is a critical aspect of the baryon cycle of galaxies, yet our understanding of this process is limited. This study focuses on understanding dust depletion effects in the low metallicity regime (< 20% Zo) typical of cosmic noon. Using medium-resolution UV spectroscopy from the COS onboard the Hubble Space Telescope, gas-phase abundances and depletions of iron and sulfur were derived toward 18 sightlines in local dwarf galaxies IC 1613 and Sextans A. The results show that the depletion of Fe and S is consistent with that found in the Milky Way, LMC and SMC. The depletion level of Fe increases with gas column density, indicating dust growth in the interstellar medium (ISM). The level of Fe depletion decreases with decreasing metallicity, resulting in the fraction of iron in gas ranging from 3% in the MW to 9% in IC 1613 and ~19% in Sextans A. The dust-to-gas and dust-to-metal ratios (D/G, D/M) for these dwarf galaxies were estimated based on the MW relations between the depletion of Fe and other elements. The study finds that D/G decreases only slightly sub-linearly with metallicity, with D/M decreasing from 0.41 +/- 0.05 in the MW to 0.11 +/- 0.11 at 0.10 Zo (at log N(H) = 21 cm-2). The trend of D/G vs. metallicity using depletion in local systems is similar to that inferred in Damped Ly-alpha systems from abundance ratios but lies higher than the trend inferred from FIR measurements in nearby galaxies.

The Large Magellanic Cloud (LMC) is home to many HII regions, which may lead to significant outflows. We examine the LMC's multiphase gas ($T\sim10^{4-5}$ K) in HI, SII, SiIV, and CIV using 110 stellar sight lines from the HST's Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) program. We develop a continuum fitting algorithm based on the concept of Gaussian Process regression and identify reliable LMC interstellar absorption over $v_{\rm helio}=175-375$ km s$^{-1}$. Our analyses show disk-wide ionized outflows in SiIV and CIV across the LMC with bulk velocities of $|v_{\rm out, bulk}|\sim20-60$ km s$^{-1}$, which indicates that most of the outflowing mass is gravitationally bound. The outflows' column densities correlate with the LMC's star formation rate surface densities ($\Sigma_{\rm SFR}$), and the outflows with higher $\Sigma_{\rm SFR}$ tend to be more ionized. Considering outflows from both sides of the LMC as traced by CIV, we conservatively estimate a total outflow rate of $\dot{M}_{\rm out}\gtrsim 0.03~M_\odot {\rm yr}^{-1}$ and a mass loading factor of $\eta\gtrsim 0.15$. We compare the LMC's outflows with those detected in starburst galaxies and simulation predictions, and find a universal scaling relation of $|v_{\rm out, bulk}|\propto \Sigma_{\rm SFR}^{0.23}$ over a wide range of star-forming conditions ($\Sigma_{\rm SFR}\sim10^{-4.5}-10^{2}~M_\odot {\rm yr}^{-1} {\rm kpc}^{-2}$). Lastly, we find that the outflows are co-rotating with the LMC's young stellar disk and the velocity field does not seem to be significantly impacted by external forces; we thus speculate on the existence of a bow shock leading the LMC, which may have shielded the outflows from ram pressure as the LMC orbits the Milky Way., Comment: Accepted for publication in ApJ. Key findings can be found in Figures 9-12. Main updates from the previous version include new estimates on star formation rate surface densities. Normalized SII, SiIV, and CIV line spectra derived for this work are published as a High Level Science Product called LMC-FLOWS (doi: 10.17909/hz0m-np43), available on website: https://archive.stsci.edu/hlsp/lmc-flows

We combine datasets from the CGM$^{2}$ and CASBaH surveys to model a transition point, $R_{\rm cross}$, between circumgalactic and intergalactic media (CGM and IGM, respectively). In total, our data consist of 7244 galaxies at z < 0.5 with precisely measured spectroscopic redshifts, all having impact parameters of 0.01 - 20 comoving Mpc from 28 QSO sightlines with high-resolution UV spectra that cover H I Ly$\alpha$. Our best-fitting model is an exclusionary two-component model that combines a 3D absorber-galaxy cross correlation function with a simple Gaussian profile at inner radii to represent the CGM. By design, this model gives rise to a determination of $R_{\rm cross}$ as a function of galaxy stellar mass, which can be interpreted as the boundary between the CGM and IGM. For galaxies with $10^8 \leq M_{\star}/M_{\odot} \leq 10^{10.5}$, we find that $R_{\rm cross}(M_{\star}) \approx 2 \pm 0.6 R_{\rm vir}$. Additionally, we find excellent agreement between $R_{\rm cross}(M_{\star})$ and the theoretically-determined splashback radius for galaxies in this mass range. Overall, our results favor models of galaxy evolution at z < 0.5 that distribute $T \approx 10^{4}$K gas to distances beyond the virial radius., Comment: Submitted to ApJ, 17 pages, 8 figures

This study addresses how the incidence rate of strong O VI absorbers in a galaxy's circumgalactic medium (CGM) depends on galaxy mass and, independently, on the amount of star formation in the galaxy. We use HST/COS absorption spectroscopy of quasars to measure O VI absorption within 400 projected kpc and 300 km s$^{-1}$ of 52 $M_{*}\sim 10^{10}$ $M_\odot$ galaxies. The galaxies have redshifts $0.12

The evolution of the metal content of the universe can be tracked through rest-frame UV spectroscopy of damped Ly-$\alpha$ systems (DLAs). Gas-phase abundances in DLAs must be corrected for dust depletion effects, which can be accomplished by calibrating the relation between abundance ratios such as [Zn/Fe] and depletions (the fraction of metals in gas, as opposed to dust). Using samples of gas-phase abundances and depletions in the Milky Way (MW), LMC, and SMC, we demonstrate that the relation between [Zn/Fe] and other abundance ratios does not change significantly between these local galaxies and DLAs, indicating that [Zn/Fe] should trace depletions of heavy elements in those systems. The availability of photospheric abundances in young massive stars, a proxy for the total (gas+dust) metallicity of neutral gas, in the MW LMC, and SMC allows us to calibrate the relation between [Zn/Fe] and depletions in these nearby galaxies. We apply the local calibrations of depletions to DLA systems. We find that the fraction of metals in dust, the dust-to-gas-ratio, and total abundances are 2-5 times lower than inferred from previous depletion calibrations based on MW measurements and a different formalism. However, the trend of dust abundance vs. metallicity remains only slightly sub-linear for all existing depletion calibrations, contrary to what is inferred from FIR, 21 cm, and CO emission in nearby galaxies and predicted by chemical evolution models. Observational constraints on the FIR dust opacity and depletions at metallicities lower than 20\% solar will be needed to resolve this tension., Comment: 37 pages; 22 figures

This work presents a new detection of H$_2$ absorption arising in a high velocity cloud (HVC) associated with either the Milky Way or the Large Magellanic Cloud (LMC). The absorber was found in an archival Far Ultraviolet Spectroscopic Explorer spectrum of the LMC star Sk-70$^\circ$32. This is the fifth well-characterized H$_2$ absorber to be found in the Milky Way's halo and the second such absorber outside the Magellanic Stream and Bridge. The absorber has a local standard of rest central velocity of $+$140 km s$^{-1}$ and a H$_2$ column density of $10^{17.5}$ cm$^{-2}$. It is most likely part of a cool and relatively dense inclusion ($T\approx 75$ K, $n_{\rm H}\sim 100$ cm$^{-3}$) in a warmer and more diffuse halo cloud. This halo cloud may be part of a still-rising Milky Way Galactic fountain flow or an outflow from the Large Magellanic Cloud., Comment: Accepted to ApJ. 13 pages, 6 figures

In this work we demonstrate that the Perseus Arm is not a continuous structure of molecular gas in the second quadrant. We first show that the observed, distanced-resolved velocity structure of the Galaxy in the outer disk is capable of creating illusory spiral arms, as was first proposed by Burton (1971). Second, we measure the distances to a collection of CO clouds at velocities consistent with the Perseus arm with $135^\circ < l < 160^\circ$. We find these distances using 3D dust maps from Green et al. (2019). We determine that these molecular cloud do not preferentially lie at the distance of a purported Perseus arm, but rather extend over 3 kpc in distance, with some evidence for a closer, high pitch angle structure between 1 and 1.5 kpc away. Finally, we demonstrate that velocity perturbations of the amplitude found near the Perseus arm can wreak havoc on our interpretation of the longitude-velocity diagram for more than half of the Milky Way disk., Comment: 14 pages, 7 figures. Accepted to the Astrophysical Journal

For over 60 years, the scientific community has studied actively growing central super-massive black holes (active galactic nuclei -- AGN) but fundamental questions on their genesis remain unanswered. Numerical simulations and theoretical arguments show that black hole growth occurs during short-lived periods ($\sim$ 10$^{7}$ -10$^{8}$ yr) of powerful accretion. Major mergers are commonly invoked as the most likely dissipative process to trigger the rapid fueling of AGN. If the AGN-merger paradigm is true, we expect galaxy mergers to coincide with black hole accretion during a heavily obscured AGN phase (N$_H$ $ > 10^{23}$ cm$^{-2}$). Starting from one of the largest samples of obscured AGN at 0.5 $<$ $z$ $<$ 3.1, we select 40 non-starbursting lower-luminosity obscured AGN. We then construct a one-to-one matched redshift- and near-IR magnitude-matched non-starbursting inactive galaxy control sample. Combining deep color \textit{Hubble Space Telescope} imaging and a novel method of human classification, we test the merger-AGN paradigm prediction that heavily obscured AGN are strongly associated with galaxies undergoing a major merger. On the total sample of 80 galaxies, we estimate each individual classifier's accuracy at identifying merging galaxies/post-merging systems and isolated galaxies. We calculate the probability of each galaxy being in either a major merger or isolated system, given the accuracy of the human classifiers and the individual classifications of each galaxy. We do not find statistically significant evidence that obscured AGN at cosmic noon are predominately found in systems with evidence of significant merging/post-merging features., Comment: 19 pages, 11 figures, accepted in ApJ

Significant galaxy mergers throughout cosmic time play a fundamental role in theories of galaxy evolution. The widespread usage of human classifiers to visually assess whether galaxies are in merging systems remains a fundamental component of many morphology studies. Studies that employ human classifiers usually construct a control sample, and rely on the assumption that the bias introduced by using humans will be evenly applied to all samples. In this work, we test this assumption and develop methods to correct for it. Using the standard binomial statistical methods employed in many morphology studies, we find that the merger fraction, error, and the significance of the difference between two samples are dependent on the intrinsic merger fraction of any given sample. We propose a method of quantifying merger biases of individual human classifiers and incorporate these biases into a full probabilistic model to determine the merger fraction and the probability of an individual galaxy being in a merger. Using 14 simulated human responses and accuracies, we are able to correctly label a galaxy as ''merger'' or ''isolated'' to within 1\% of the truth. Using 14 real human responses on a set of realistic mock galaxy simulation snapshots our model is able to recover the pre-coalesced merger fraction to within 10\%. Our method can not only increase the accuracy of studies probing the merger state of galaxies at cosmic noon, but also can be used to construct more accurate training sets in machine learning studies that use human classified data-sets., Comment: 23 pages, 16 figures, accepted for publication in ApJ, https://github.com/elambrid/merger_or_not

We use hydrodynamical simulations of two Milky Way-mass galaxies to demonstrate the impact of cosmic-ray pressure on the kinematics of cool and warm circumgalactic gas. Consistent with previous studies, we find that cosmic-ray pressure can dominate over thermal pressure in the inner 50 kpc of the circumgalactic medium (CGM), creating an overall cooler CGM than that of similar galaxy simulations run without cosmic rays. We generate synthetic sightlines of the simulated galaxies' CGM and use Voigt profile fitting methods to extract ion column densities, Doppler-b parameters, and velocity centroids of individual absorbers. We directly compare these synthetic spectral line fits with HST/COS CGM absorption-line data analyses, which tend to show that metallic species with a wide range of ionization potential energies are often kinematically aligned. Compared to the Milky-Way simulation run without cosmic rays, the presence of cosmic-ray pressure in the inner CGM creates narrower OVI absorption features and broader SiIII absorption features, a quality which is more consistent with observational data. Additionally, because the cool gas is buoyant due to nonthermal cosmic-ray pressure support, the velocity centroids of both cool and warm gas tend to align in the simulated Milky Way with feedback from cosmic rays. Our study demonstrates that detailed, direct comparisons between simulations and observations, focused on gas kinematics, have the potential to reveal the dominant physical mechanisms that shape the CGM., Comment: 26 pages including 19 figures and appendix. Accepted by ApJ. Updated version includes additional figures, tests, and expanded discussion. Supplemental figures at irynabutsky.me/data

A key component of the baryon cycle in galaxies is the depletion of metals from the gas to the dust phase in the neutral ISM. The METAL (Metal Evolution, Transport and Abundance in the Large Magellanic Cloud) program on the Hubble Space Telescope acquired UV spectra toward 32 sightlines in the half-solar metallicity LMC, from which we derive interstellar depletions (gas-phase fractions) of Mg, Si, Fe, Ni, S, Zn, Cr, and Cu. The depletions of different elements are tightly correlated, indicating a common origin. Hydrogen column density is the main driver for depletion variations. Correlations are weaker with volume density, probed by CI fine structure lines, and distance to the LMC center. The latter correlation results from an East-West variation of the gas-phase metallicity. Gas in the East, compressed side of the LMC encompassing 30 Doradus and the Southeast HI over-density is enriched by up to +0.3dex, while gas in the West side is metal-deficient by up to -0.5dex. Within the parameter space probed by METAL, no correlation with molecular fraction or radiation field intensity are found. We confirm the factor 3-4 increase in dust-to-metal and dust-to-gas ratios between the diffuse (logN(H)~20 cm-2) and molecular (logN(H)~22 cm-2) ISM observed from far-infrared, 21 cm, and CO observations. The variations of dust-to-metal and dust-to-gas ratios with column density have important implications for the sub-grid physics of chemical evolution, gas and dust mass estimates throughout cosmic times, and for the chemical enrichment of the Universe measured via spectroscopy of damped Lyman-alpha systems.

We present initial results from the \textit{COS and Gemini Mapping the Circumgalactic Medium} (\mbox{CGMCGM} $\equiv$ CGM$^{2}$) survey. The CGM$^{2}$ survey consists of 1689 galaxies, all with high-quality Gemini GMOS spectra, within 1 Mpc of twenty-two $z \lesssim 1$ quasars, all with S/N$\sim$10 {\emph{HST/COS}} G130M$+$G160M spectra. For 572 of these galaxies having stellar masses $10^{7} M_{\odot} < M_{\star} < 10^{11} M_{\odot}$ and $z \lesssim 0.5$, we show that the \ion{H}{1} covering fraction above a threshold of \NHI$>10^{14} $cm$^{-2}$ is $\gtrsim 0.5$ within 1.5 virial radii ($R_{\rm vir} \sim R_{200m}$). We examine the \ion{H}{1} kinematics and find that the majority of absorption lies within $\pm$ 250 km s$^{-1}$ of the galaxy systemic velocity. We examine \ion{H}{1} covering fractions over a range of impact parameters to infer a characteristic size of the CGM, $R^{14}_{\rm CGM}$, as a function of galaxy mass. $R^{14}_{\rm CGM}$ is the impact parameter at which the probability of observing an absorber with \NHI $>$ 10$^{14}$ cm$^{-2}$ is $>$ 50\%. In this framework, the radial extent of the CGM of $M_{\star} > 10^{9.9} M_{\odot}$ galaxies is $R^{14}_{\rm CGM} = 346^{+57}_{-53}$ kpc or $R^{14}_{\rm CGM} \simeq 1.2R_{\rm vir}$. Intermediate-mass galaxies with $10^{9.2} < M_{\star}/M_{\odot} < 10^{9.9}$ have an extent of $R^{14}_{\rm CGM} = 353^{+64}_{-50}$ kpc or $R^{14}_{\rm CGM} \simeq 2.4R_{\rm vir}$. Low-mass galaxies, $M_{\star} < 10^{9.2} M_{\odot}$, show a smaller physical scale $R^{14}_{\rm CGM} = 177_{-65}^{+70}$ kpc and extend to $R^{14}_{\rm CGM} \simeq 1.6R_{\rm vir}$. Our analysis suggests that using $R_{\rm vir}$ as a proxy for the characteristic radius of the CGM likely underestimates its extent., Comment: 30 pages, 14 figures. Submitted to ApJ, Comments welcome

Photodissociation regions (PDRs) are parts of the ISM consisting of predominantly neutral gas, located at the interface between H II regions and molecular clouds. The physical conditions within these regions show variations on very short spatial scales, and therefore PDRs constitute ideal laboratories for investigating the properties and evolution of dust grains. We have mapped IC 63 at high resolution from the UV to the NIR (275 nm to 1.6 $\mu$m), using the Hubble Space Telescope WFC3. Using a Bayesian SED fitting tool, we simultaneously derive a set of stellar ($T_\text{eff}$, $\log(g)$, distance) and extinction ($A_V$, $R_V$) parameters for 520 background stars. We present maps of $A_V$ and $R_V$ with a resolution of 25 arcsec based on these results. The extinction properties vary across the PDR, with values for $A_V$ between 0.5 and 1.4 mag, and a decreasing trend in $R_V$, going from 3.7 at the front of the nebula to values as low as 2.5 further in. This provides evidence for evolution of the dust optical properties. We fit two modified blackbodies to the MIR and FIR SED, obtained by combining the $A_V$ map with data from Spitzer and Herschel. We derive effective temperatures (30 K and 227 K) and the ratio of opacities at 160 $\mu$m to V band $\kappa_{160} / \kappa_V$ ($7.0 \times 10^{-4}$ and $2.9 \times 10^{-9}$) for the two dust populations. Similar fits to individual pixels show spatial variations of $\kappa_{160} / \kappa_{V}$. The analysis of our HST data, combined with these Spitzer and Herschel data, provides the first panchromatic view of dust within a PDR., Comment: 28 pages, 20 figures, accepted for publication in The Astrophysical Journal

The comparison of chemical abundances in the neutral gas of galaxies to photospheric abundances of old and young stars, ionized gas abundances, and abundances in galactic halos can trace the chemical enrichment of the universe through cosmic times. In particular, our understanding of chemical enrichment through spectroscopic observations of damped Lyman alpha systems (DLAs) relies on corrections for depletion of metals from the gas to the dust phase. These corrections must be determined in the nearby universe, where both gas-phase abundances and photospheric abundances of young stars recently formed out of the interstellar medium can be measured. Multi-object high-resolution (R>50,000) ultraviolet (970-2400 A) and optical (300-600 nm) spectroscopy toward massive stars in local volume galaxies (D < 15 Mpc) covering a wide range of metallicities (a few % solar to solar) and morphological types will provide the abundance and depletion measurements needed to obtain a detailed and comprehensive characterization of the lifecycle of metals in neutral gas and dust in galaxies, thereby observationally addressing important questions about chemical enrichment and galaxy evolution., Comment: 7 pages, 3 figures, ASTRO2020 white paper

Absorption line spectroscopy is a powerful way of measuring properties of stars and the interstellar medium. Absorption spectra are often analyzed manually, an approach that limits reproducibility and which cannot practically be applied to modern datasets consisting of thousands or even millions of spectra. Simultaneous probabilistic modeling of absorption features and continuum shape is a promising approach for automating this analysis. Existing implementations of this approach use numerical methods such as Markov chain Monte Carlo (MCMC) to marginalize over the continuum parameters. Numerical marginalization over large numbers of continuum parameters is too slow to be convenient for exploratory analysis, can increase the dimensionality of an inference problem beyond the capacity of simple MCMC samplers, and is in general impractical for the analysis of large datasets. When continua are parameterized as linear functions such as polynomials or splines, it is possible to reduce continuum parameter marginalization to an integral over a multivariate normal distribution, which has a known closed form. In addition to speeding up probabilistic modeling, analytic marginalization makes it trivial to marginalize over continuum parameterizations and to combine continuum description marginalization with optimization for absorption line parameters. These new possibilities allow automatic, probabilistically justified continuum placement in analyses of large spectroscopic datasets. We compare the accuracy to within which absorption line parameters can be recovered using different continuum placement methods and find that marginalization is in many cases an improvement over other methods. We implement analytic marginalization over linear continuum parameters in the open-source package amlc., Comment: Submitted to AAS journals

Metal Evolution, Transport, and Abundance in the LMC (METAL) is a large Cycle 24 program on the Hubble Space Telescope aimed at measuring dust extinction properties and interstellar depletions in the Large Magellanic Cloud (LMC) at half-solar metallicity. The 101-orbit program is comprised of COS and STIS spectroscopy toward 33 LMC massive stars between 1150 A and 3180 A, and parallel WFC3 imaging in 7 NUV-NIR filters. The fraction of silicon in the gas-phase (depletion) obtained from the spectroscopy decreases with increasing hydrogen column density. Depletion patterns for silicon differ between the Milky Way, LMC, and SMC, with the silicon depletion level offsetting almost exactly the metallicity differences, leading to constant gas-phase abundances in those galaxies for a given hydrogen column density. The silicon depletion correlates linearly with the absolute-to-selective extinction, R$_V$, indicating a link between gas depletion and dust grain size. Extinction maps are derived from the resolved stellar photometry in the parallel imaging, which can be compared to FIR images from Herschel and Spitzer to estimate the emissivity of dust at LMC metallicity. The full METAL sample of depletions, UV extinction curves, and extinction maps will inform the abundance, size, composition, and optical properties of dust grains in the LMC, comprehensively improve our understanding of dust properties, and the accuracy with which dust-based gas masses, star formation rates and histories in nearby and high-redshift galaxies are estimated. This overview paper describes the goals, design, data reduction, and initial results of the METAL survey., Comment: 28 pages, 15 figures

Quasar outflows carry mass, momentum and energy into the surrounding environment, and have long been considered a potential key factor in regulating the growth of supermassive black holes and the evolution of their host galaxies. A crucial parameter for understanding the origin of these outflows and measuring their influence on their host galaxies is the distance (R) between the outflow gas and the galaxy center. While R has been measured in a number of individual galaxies, its distribution remains unknown. Here we report the distributions of R and the kinetic luminosities of quasars outflows, using the statistical properties of broad absorption line variability in a sample of 915 quasars from the Sloan Digital Sky Surveys. The mean and standard deviation of the distribution of R are 10^{1.4+/-0.5} parsecs. The typical outflow distance in this sample is tens of parsec, which is beyond the theoretically predicted location (0.01 ~ 0.1 parsecs) where the accretion disc line-driven wind is launched, but is smaller than the scales of most outflows that are derived using the excited state absorption lines. The typical value of the mass-flow rate is of tens to a hundred solar masses per year, or several times the accretion rate. The typical kinetic-to-bolometric luminosity ratio is a few per cent, indicating that outflows are energetic enough to influence the evolution of their host galaxies., Comment: Accepted by Nature Astronomy, 8 pages and 2 figures(main text), 8 supplementary figures

We analyse 2,015 mid-infrared (MIR) spectra of galaxies observed with Spitzer's Infrared Spectrograph, including objects with growing super-massive black holes and objects where most of the infrared emission originates from newly formed stars. We determine if and how accreting super-massive black holes at the centre of galaxies -- known as active galactic nuclei (AGN) -- heat and ionize their host galaxies' dust and molecular gas. We use four MIR diagnostics to estimate the contribution of the AGN to the total MIR emission. We refer to galaxies whose AGN contribute more than 50 per cent of the total MIR emission as AGN-dominated. We compare the relative strengths of PAH emission features and find that PAH grains in AGN-dominated sources have a wider range of sizes and fractional ionizations than PAH grains in non-AGN dominated sources. We measure rotational transitions of H_2 and estimate H_2 excitation temperatures and masses for individual targets, H_2 excitation temperatures for spectra stacked by their AGN contribution to the MIR, and the H_2 excitation temperature distributions via a hierarchical Bayesian model. We find an average 200 K difference between the excitation temperatures of the H_2 S(5) and H_2 S(7) pure rotational molecular hydrogen transition pair in AGN-dominated versus non-AGN dominated galaxies. Our findings suggest that AGN impact the interstellar medium of their host galaxies., Comment: Published in MNRAS, 20 pages, 16 figures, link to full tables and stacked spectra (http://pages.jh.edu/~elambri1/MIR_AGN_MolecularGas/)

In this work, we derive a spatially resolved map of the line-of-sight velocity of the interstellar medium and use it, along with a second map of line-of-sight velocity from Paper I of this series, to determine the nature of gaseous spiral structure in the Milky Way. This map is derived from measurements of the 1.527 $\mu$m diffuse interstellar band (DIB) in stellar spectra from the APOGEE survey and covers the nearest 4-5 kpc of the Northern Galactic plane. We cross-check this new DIB-based line-of-sight velocity map with the map derived in Paper I and find that they agree. We then compare these maps with line-of-sight velocity maps derived from simulations of quasi-stationary density wave spiral structure and dynamic, or material, spiral structure in a Milky Way-like galaxy. While none of the maps derived from these simulations is an exact match to the measured velocity field of the Milky Way, the measurements are more consistent with simulations of dynamic spiral structure. In the dynamic spiral structure simulation that best matches the measurements, the Perseus spiral arm is being disrupted., Comment: Submitted to AAS journals

We use deep $HST$ WFC3/IR imaging to study the Initial Mass Function (IMF) of the ultra faint dwarf galaxy Coma Berenices (Com Ber). Our observations reach the lowest stellar mass ever probed in a resolved galaxy, with 50\% completeness at $\sim 0.17$ M$_{\odot}$. Unresolved background galaxies however limit our purity below $\sim 0.23$ M$_{\odot}$. If modeled with a single power law, we find that the IMF slope is $-1.45^{+0.29}_{-0.3}$ (68\% credible intervals), compared to a Milky Way value of $-2.3$. For a broken power law, we obtain a low-mass slope of $-1.18_{-0.33}^{+0.49}$, a high-mass slope of $-1.88_{-0.49}^{+0.43}$ and a break mass of $0.57_{-0.08}^{+0.12}$ M$_{\odot}$, compared to $-1.3$, $-2.3$ and 0.5 M$_{\odot}$ for a Kroupa IMF. For a log-normal IMF model we obtain values of $0.33_{-0.16}^{+0.15}$ M$_{\odot}$ for the location parameter and of $0.68_{-0.12}^{+0.17}$ for $\sigma$ (0.22 M$_{\odot}$ and 0.57 for the Chabrier system IMF). All three parametrizations produce similar agreement with the data. Our results agree with previous analysis of shallower optical HST data. However analysis of similar optical data of other dwarfs finds IMFs significantly more bottom-light than in the Milky Way. These results suggest two, non mutually exclusive, possibilities: that the discrepancy of the dwarf galaxies IMF with respect to the Milky Way is, at least partly, an artifact of using a single power law model, and that there is real variance in the IMF at low masses between the currently studied nearby dwarfs, with Com Ber being similar to the Milky Way, but other dwarfs differing significantly., Comment: ApJ accepted

Using deep HST/ACS observations, we demonstrate that the sub-solar stellar initial mass function (IMF) of 6 ultra-faint dwarf Milky Way Satellites (UFDs) is more bottom light than the IMF of the Milky Way disk. Our data have a lower mass limit of about 0.45 M$_{\odot}$, while the upper limit is $\sim 0.8$ M$_\odot$, set by the turn-off mass of these old, metal poor systems. If formulated as a single power law, we obtain a shallower IMF slope than the "Salpeter" value of $-2.3$, ranging from $-1.01$ for Leo IV, to $-1.87$ for Bo\"otes I. The significance of such deviations depends on the galaxy and is typically 95\% or more. When modeled as a log-normal, the IMF fit results in a larger peak mass than in the Milky Way disk, however a Milky Way disk value for the characteristic system mass ($\sim0.22$ M$_{\odot}$) is excluded only at 68\% significance, and only for some UFDs in the sample. We find that the IMF slope correlates well with the galaxy mean metallicity and, to a lesser degree, with the velocity dispersion and the total mass. The strength of the observed correlations is limited by shot noise in the number of observed stars, but future space-based missions like JWST and WFIRST will both enhance the number of dwarf Milky Way Satellites that can be studied in such detail, and the observation depth for individual galaxies., Comment: Accepted for publication in ApJ on 18 Jan 2018

We have developed a method for deriving the distribution of the Milky Way's interstellar medium as a function of longitude, latitude, distance and line-of-sight velocity. This method takes as input maps of reddening as a function of longitude, latitude, and distance and maps of line emission as a function of longitude, latitude, and line-of-sight velocity. We have applied this method to datasets covering much of the Galactic plane. The output of this method correctly reproduces the line-of-sight velocities of high-mass star forming regions with known distances from Reid et al. (2014) and qualitatively agrees with results from the Milky Way kinematics literature. These maps will be useful for measuring flows of gas around the Milky Way's spiral arms and into and out of giant molecular clouds., Comment: 11 pages, 5 figures, accepted for publication in The Astronomical Journal

We present the Bayesian Extinction And Stellar Tool (BEAST), a probabilistic approach to modeling the dust extinguished photometric spectral energy distribution of an individual star while accounting for observational uncertainties common to large resolved star surveys. Given a set of photometric measurements and an observational uncertainty model, the BEAST infers the physical properties of the stellar source using stellar evolution and atmosphere models and constrains the line of sight extinction using a newly developed mixture model that encompasses the full range of dust extinction curves seen in the Local Group. The BEAST is specifically formulated for use with large multi-band surveys of resolved stellar populations. Our approach accounts for measurement uncertainties and any covariance between them due to stellar crowding (both systematic biases and uncertainties in the bias) and absolute flux calibration, thereby incorporating the full information content of the measurement. We illustrate the accuracy and precision possible with the BEAST using data from the Panchromatic Hubble Andromeda Treasury. While the BEAST has been developed for this survey, it can be easily applied to similar existing and planned resolved star surveys., Comment: 20 pages, 19 figures, ApJ, in press

We study the stellar mass Tully-Fisher relation (TFR, stellar mass versus rotation velocity) for a morphologically blind selection of emission line galaxies in the field at redshifts 0.1 $<$ z $<$ 0.375. Kinematics ($\sigma_g$, V$_{rot}$) are measured from emission lines in Keck/DEIMOS spectra and quantitative morphology is measured from V- and I-band Hubble images. We find a transition stellar mass in the TFR, $\log$ M$_*$ = 9.5 M$_{\odot}$. Above this mass, nearly all galaxies are rotation-dominated, on average more morphologically disk-like according to quantitative morphology, and lie on a relatively tight TFR. Below this mass, the TFR has significant scatter to low rotation velocity and galaxies can either be rotation-dominated disks on the TFR or asymmetric or compact galaxies which scatter off. We refer to this transition mass as the "mass of disk formation", M$_{\mathrm{df}}$ because above it all star-forming galaxies form disks (except for a small number of major mergers and highly star-forming systems), whereas below it a galaxy may or may not form a disk., Comment: Accepted to MNRAS. 14 pages, 7 figures

We present a new method for deriving the stellar birth function (SBF) of resolved stellar populations. The SBF (stars born per unit mass, time, and metallicity) is the combination of the initial mass function (IMF), the star-formation history (SFH), and the metallicity distribution function (MDF). The framework of our analysis is that of Poisson Point Processes (PPPs), a class of statistical models suitable when dealing with points (stars) in a multidimensional space (the measurement space of multiple photometric bands). The theory of PPPs easily accommodates the modeling of measurement errors as well as that of incompleteness. Compared to most of the tools used to study resolved stellar populations, our method avoids binning stars in the color-magnitude diagram and uses the entirety of the information (i.e., the whole likelihood function) for each data point; the proper combination of the individual likelihoods allows the computation of the posterior probability for the global population parameters. This includes unknowns such as the IMF slope and combination of SFH and MDF, which are rarely solved for simultaneously in the literature, however entangled and correlated they might be. Our method also allows proper inclusion of nuisance parameters, such as distance and extinction distributions. The aim of this paper, is to assess the validity of this new approach under a range of assumptions, using only simulated data. Forthcoming work will show applications to real data. Although it has a broad scope of possible applications, we have developed this method to study multi-band HST observations of the Milky Way Bulge. Therefore we will focus on simulations with characteristics similar to those of the Galactic Bulge., Comment: Accepted for publication in ApJ

We present a study of the composition of gas and dust in the Large and Small Magellanic Clouds (LMC and SMC, together -- the MCs) as measured by UV absorption spectroscopy. We have measured P II and Fe II along 85 sightlines toward the MCs using archival FUSE observations. For 16 of those sightlines, we have measured Si II, Cr II, and Zn II from new HST COS observations. We have combined these measurements with H I and H$_2$ column densities and reference stellar abundances from the literature to derive gas-phase abundances, depletions, and gas-to-dust ratios (GDRs). 80 of our 84 P measurements and 13 of our 16 Zn measurements are depleted by more than 0.1 decades, suggesting that P and Zn abundances are not accurate metallicity indicators at and above the metallicity of the SMC. The maximum P and Zn depletions are the same in the MW, LMC, and SMC. Si, Cr, and Fe are systematically less depleted in the SMC than in the MW or LMC. The minimum Si depletion in the SMC is consistent with zero. Our depletion-derived GDRs broadly agree with GDRs from the literature. The GDR varies from location to location within a galaxy by a factor of up to 2 in the LMC and up to 5 in the SMC. This variation is evidence of dust destruction and/or growth in the diffuse neutral phase of the interstellar medium., Comment: Submitted to ApJ

The nature, composition, abundance, and size distribution of dust in galaxies is determined by the rate at which it is created in the different stellar sources and destroyed by interstellar shocks. Because of their extensive wavelength coverage, proximity, and nearly face-on geometry, the Magellanic Clouds (MCs) provide a unique opportunity to study these processes in great detail. In this paper we use the complete sample of supernova remnants (SNRs) in the MCs to calculate the lifetime and destruction efficiencies of silicate and carbon dust in these galaxies. We find dust lifetimes of 22 +- 13 Myr (30 +- 17 Myr) for silicate (carbon) grains in the LMC, and 54 +- 32 Myr (72 +- 43 Myr) for silicate (carbon) grains in the SMC. The significantly shorter lifetimes in the MCs, as compared to the Milky Way, are explained as the combined effect of their lower total dust mass, and the fact that the dust-destroying isolated SNe in the MCs seem to be preferentially occurring in regions with higher than average dust-to-gas (D2G) mass ratios. We also calculate the supernova rate and the current star formation rate in the MCs, and use them to derive maximum dust injection rates by asymptotic giant branch stars and core collapse supernovae. We find that the injection rates are an order of magnitude lower than the dust destruction rates by the SNRs. This supports the conclusion that, unless the dust destruction rates have been considerably overestimated, most of the dust must be reconstituted from surviving grains in dense molecular clouds. More generally, we also discuss the dependence of the dust destruction rate on the local D2G mass ratio, the ambient gas density and metallicity, as well as the application of our results to other galaxies and dust evolution models., Comment: 15 pages, 8 figures, 5 tables, accepted to ApJ

The dust properties in the Large and Small Magellanic Clouds are studied using the HERITAGE Herschel Key Project photometric data in five bands from 100 to 500 micron. Three simple models of dust emission were fit to the observations: a single temperature blackbody modified by a power- law emissivity (SMBB), a single temperature blackbody modified by a broken power-law emissivity (BEMBB), and two blackbodies with different temperatures, both modified by the same power-law emissivity (TTMBB). Using these models we investigate the origin of the submm excess; defined as the submillimeter (submm) emission above that expected from SMBB models fit to observations < 200 micron. We find that the BEMBB model produces the lowest fit residuals with pixel-averaged 500 micron submm excesses of 27% and 43% for the LMC and SMC, respectively. Adopting gas masses from previous works, the gas-to-dust ratios calculated from our the fitting results shows that the TTMBB fits require significantly more dust than are available even if all the metals present in the interstellar medium (ISM) were condensed into dust. This indicates that the submm excess is more likely to be due to emissivity variations than a second population of colder dust. We derive integrated dust masses of (7.3 +/- 1.7) x 10^5 and (8.3 +/- 2.1) times 10^4 M(sun) for the LMC and SMC, respectively. We find significant correlations between the submm excess and other dust properties; further work is needed to determine the relative contributions of fitting noise and ISM physics to the correlations., Comment: 20 pages, 9 figures, accepted for publication in the ApJ