We measure the effective opacity ($\tau_{eff}$) of the Intergalactic Medium (IGM) from the composite spectra of 281 Lyman-Break Galaxies (LBGs) in the redshift range $2 \lesssim z \lesssim 3$. Our spectra are taken from the COSMOS Lyman-Alpha Mapping And Tomographic Observations (CLAMATO) survey derived from the Low Resolution Imaging Spectrometer (LRIS) on the W.M. Keck I telescope. We generate composite spectra in two redshift intervals and fit them with spectral energy distribution (SED) models composed of simple stellar populations. Extrapolating these SED models into the Ly$\alpha$ forest, we measure the effective Ly$\alpha$ opacity ($\tau_{eff}$) in the $2.02 \leq z \leq 2.44$ range. At $z = 2.22$, we estimate $\tau_{eff} = 0.159 \pm 0.001$ from a power-law fit to the data. These measurements are consistent with estimates from quasar analyses at $z<2.5$ indicating that the systematic errors associated with normalizing quasar continua are not substantial. We provide a Gaussian Processes model of our results and previous $\tau_{eff}$ measurements that describes the steep redshift evolution in $\tau_{eff}$ from $z = 1.5 - 4$.

Identifying the mechanisms driving the escape of Lyman Continuum (LyC) photons is crucial to find Lyman Continuum Emitter (LCE) candidates. To understand the physical properties involved in the leakage of LyC photons, we investigate the connection between the HI covering fraction, HI velocity width, the Lyman alpha (LyA) properties and escape of LyC photons in a sample of 22 star-forming galaxies including 13 LCEs. We fit the stellar continua, dust attenuation, and absorption lines between 920 and 1300 A to extract the HI covering fractions and dust attenuation. Additionally, we measure the HI velocity widths of the optically thick Lyman series and derive the LyA equivalent widths (EW), escape fractions (fesc), peak velocities and fluxes at the minimum of the LyA profiles. Overall, we highlight strong correlations between the presence of low HI covering fractions and (1) low LyA peak velocities; (2) more flux at the profile minimum; and (3) larger EW(LyA), fesc(LyA), and fesc(LyC). Hence, low column density channels are crucial ISM ingredients for the leakage of LyC and LyA photons. Additionally, galaxies with narrower HI absorption velocity widths have higher LyA equivalent widths, larger LyA escape fractions, and lower LyA peak velocity separations. This suggests that these galaxies have low HI column density. Finally, we find that dust regulates the amount of LyA and LyC radiation that actually escapes the ISM. Overall, the ISM porosity is one origin of strong LyA emission and enables the escape of ionizing photons in low-z leakers. However, this is not enough to explain the largest fesc(LyC) observed, which indicates that the most extreme LCEs are likely density-bounded along all lines of sight to the observer. Overall, the neutral gas porosity constrains a lower limit to the escape fraction of LyC and LyA photons, providing a key estimator of the leakage of ionizing photons.

During the epoch of reionisation, neutral gas in the early Universe was ionized by hard ultraviolet radiation emitted by young stars in the first galaxies. To do so, ionizing ultraviolet photons must escape from the host galaxy. We present Hubble Space Telescope observations of the gravitationally lensed galaxy PSZ1-ARC G311.6602-18.4624, revealing bright, multiply-imaged ionizing photon escape from a compact star-forming region through a narrow channel in an optically thick gas. The gravitational lensing magnification shows how ionizing photons escape this galaxy, contributing to the re-ionization of the Universe. The multiple sight lines to the source probe absorption by intergalactic neutral hydrogen on scales of no more than a few hundred, perhaps even less than ten, parsec.

In the past decade, our understanding of how stars and galaxies formed during the first 5 billion years after the Big Bang has been revolutionized by observations that leverage gravitational lensing by intervening masses, which act as natural cosmic telescopes to magnify background sources. Previous studies have harnessed this effect to probe the distant universe at ultraviolet, optical, infrared and millimeter wavelengths. However, strong lensing studies of young, star-forming galaxies have never extended into X-ray wavelengths, which uniquely trace high-energy phenomena. Here we report an X-ray detection of star formation in a highly magnified, strongly lensed galaxy. This lensed galaxy, seen during the first third of the history of the Universe, is a low--mass, low--metallicity starburst with elevated X-ray emission, and is a likely analog to the first generation of galaxies. Our measurements yield insight into the role that X-ray emission from stellar populations in the first generation of galaxies may play in re-ionizing the Universe. This observation paves the way for future strong lensing-assisted X-ray studies of distant galaxies reaching orders of magnitude below the detection limits of current deep fields, and previews the depths that will be attainable with future X-ray observatories.

Motivated by the recent discovery of the near-ubiquity of Ly$\alpha$ emission around $z \gtrsim 3$ QSOs, we performed a systematic study of QSO circumgalactic Ly$\alpha$ emission at $z\approx2$, utilizing the unique capability of the Keck Cosmic Web Imager (KCWI) -- a new wide-field, blue sensitive integral-field spectrograph (IFU). In this paper, we present KCWI observations on a sample of 16 ultraluminous Type-I QSOs at $z=2.1-2.3$ with ionizing luminosities of $L_{\rm{\nu_{\rm{LL}}}}=10^{31.1-32.3}$ erg s$^{-1}$ Hz$^{-1}$. We found that 14 out of 16 QSOs are associated with Ly$\alpha$ nebulae with projected linear-sizes larger than 50 physical kpc (pkpc). Among them, four nebulae have enormous Ly$\alpha$ emission with the Ly$\alpha$ surface brightness $SB_{\rm{Ly\alpha}} >10^{-17}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ on the $>100$ kpc scale, extending beyond the field of view of KCWI. Our KCWI observations reveal that most $z\approx2$ QSO nebulae have a more irregular morphology compared to those at $z\gtrsim3$. In turn, we measure that the circularly-averaged surface brightness (SB) at $z\approx2$ is 0.4 dex fainter than the redshift-corrected, median SB at $z\gtrsim 3$. The Ly$\alpha$ SB profile (SB$_{\rm{Ly\alpha}}$) of QSOs at $z\approx 2$ can be described by a power law of SB$_{\rm{Ly\alpha},z\approx2.3}= 3.7\times10^{-17}\times(r/40)^{-1.8}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$, with the slope similar to that at $z\gtrsim3$. The observed lower redshift-corrected, circularly-averaged SB may be mainly due to the lower covering factor of cool gas clouds in massive halos at $z\approx2$.

We compare the molecular and ionized gas kinematics of two strongly lensed galaxies at $z\sim1$ based on observations from ALMA and MUSE. We derive the CO and [OII] rotation curves and dispersion profiles of these two galaxies. We find a difference between the observed molecular and ionized gas rotation curves for one of the galaxies, the Cosmic Snake, for which we obtain a spatial resolution of few hundred parsecs along the major axis. The rotation curve of the molecular gas is steeper than the rotation curve of the ionized gas. In the second galaxy, A521, the molecular and ionized gas rotation curves are consistent, but the spatial resolution is only of few kpc on the major axis. Using simulations, we investigate the effect of the thickness of the gas disk and effective radius on the observed rotation curves and find that a more extended and thicker disk smooths the curve. We also find that the presence of a strongly inclined thick disk (>1 kpc) can smooth the rotation curve because it degrades the spatial resolution along the line of sight. By building a model using a stellar disk and two gas disks, we reproduce the rotation curves of the Cosmic Snake with a molecular gas disk that is more massive and more radially and vertically concentrated than the ionized gas disk. Finally, we also obtain an intrinsic velocity dispersion in the Cosmic Snake of 18.5+-7 km/s and 19.5+-6 km/s for the molecular and ionized gas, respectively, which is consistent with a molecular disk with a smaller and thinner disk. For A521, the intrinsic velocity dispersion values are 11+-8 km/s and 54+-11 km/s, with a higher value for the ionized gas. This could indicate that the ionized gas disk is thicker and more turbulent in this galaxy. These results highlight the different spatial distribution of the molecular and ionized gas disks at $z\sim1$ and suggest the presence of thick ionized gas disks at this epoch.

In the last few years, prominent high-ionization nebular emission lines (i.e., OIII], CIII], CIV, HeII) have been observed in the deep UV spectra of z~5-7 galaxies, indicating that extreme radiation fields characterize reionization-era systems. These lines have been linked to the leakage of Lyman continuum photons both theoretically and observationally. Consequently, high-ionization UV emission lines present our best probe to detect and characterize the most distant galaxies that we will observe in the coming years, and are key to understanding the sources of reionization, yet the physics governing their production is poorly understood. Here we present recent high-resolution Hubble Space Telescope spectra of two nearby extreme UV emission-line galaxies, J104457 and J141851. We report the first observations of intense nebular HeII and double-peaked, resonantly-scattered CIV emission, a combination that suggests these galaxies both produce and transmit a significant number of very high-energy ionizing photons (E>47.89 eV) through relatively low column densities of high-ionization gas. This suggests that, in addition to photons at the H-ionizing edge, the very hard ionizing photons that escape from these galaxies may provide a secondary source of ionization that is currently unconstrained observationally. Simultaneous radiative transfer models of LyA and CIV are needed to understand how ionizing radiation is transmitted through both low- and high-ionization gas. Future rest-frame FUV observations of galaxies within the epoch of reionization using the JWST or ELTs will allow us to constrain the escape of helium-ionizing photons and provide an estimate for their contribution to the reionization budget.

The ionizing photons produced by massive stars are key actors in galaxy evolution. Ionizing photon production and escape is poorly understood. Improved space-based, spatially-resolved, multiplexed spectroscopic capabilities covering observed wavelengths of 1000 to 3000 Angstrom, in concert with spectroscopy from the ELTs and JWST, would lead to definitive answers as to how ionizing photons are produced and leaked, what populations of galaxies are responsible for ionizing photon leakage, what determines whether escape is possible, and how ionizing galaxy populations evolve over cosmic time.

We infer the properties of massive star populations using the far-ultraviolet stellar continua of 61 star-forming galaxies: 42 at low-z observed with HST and 19 at z~2 from the Megasaura sample. We fit each stellar continuum with a linear combination of up to 50 single age and single metallicity Starburst99 models. From these fits, we derive light-weighted ages and metallicities, which agree with stellar wind and photospheric spectral features, and infer the spectral shapes and strengths of the ionizing continua. Inferred light-weighted stellar metallicities span 0.05-1.5 Z$_\odot$ and are similar to the measured nebular metallicities. We quantify the ionizing continua using the ratio of the ionizing flux at 900\AA\ to the non-ionizing flux at 1500\AA\ and demonstrate the evolution of this ratio with stellar age and metallicity using theoretical single burst models. These single burst models only match the inferred ionizing continua of half of the sample, while the other half are described by a mixture of stellar ages. Mixed age populations produce stronger and harder ionizing spectra than continuous star formation histories, but, contrary to previous studies that assume constant star formation, have similar stellar and nebular metallicities. Stellar population age and metallicity affect the far-UV continua in different and distinguishable ways; assuming a constant star formation history diminishes the diagnostic power. Finally, we provide simple prescriptions to determine the ionizing photon production efficiency ($\xi_{ion}$) from the stellar population properties. $\xi_{ion}$ has a range of log($\xi_{ion})=24.4-25.7$ Hz erg$^{-1}$ that depends on stellar age, metallicity, star formation history, and contributions from binary star evolution. These stellar population properties must be observationally determined to determine the number of ionizing photons generated by massive stars.

Diagnosing the physical and chemical conditions within star-forming galaxies (SFGs) is of paramount importance to understanding key components of galaxy formation and evolution: star-formation, gas enrichment, outflows, and accretion. Well established optical emission-line diagnostics used to discern such properties (i.e., metal content, density, strength/shape of ionizing radiation) will be observationally inaccessible for the earliest galaxies, emphasizing the need for robust, reliable interstellar medium (ISM) diagnostics at ultraviolet (UV) wavelengths. Calibrating these UV diagnostics requires a comprehensive comparison of the UV and optical emission lines in nearby SFGs. Optical integral field unit (IFU) surveys have revealed the inhomogeneous nature of the ISM in SFGs, which leads to non-systematic biases in the interpretation of unresolved sources. Spatial variations are especially important to consider at UV wavelengths, where the strongest emission features originate from only the highest excitation regions of the nebula and are challenging to distinguish from competing high-ionization sources (e.g., shocks, AGN, etc.). Since surveys collecting large-scale optical integral field unit (IFU) spectroscopy are already underway, this white paper calls for an IFU or multi-object far-UV (FUV) spectroscopic instrument with high sensitivity, high spatial resolution, and large field of view (FoV). Given the impact of large-scale optical IFU surveys over the past decade, this white paper emphasizes the scientific need for a comparable foundation of spatially-resolved far-UV spectroscopy survey of nearby galaxies that will lay the foundation of diagnostics critical to the interpretation of the distant universe.

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

The gas-phase metallicity of low-mass galaxies increases with increasing stellar mass ($M_\ast$) and is nearly constant for high-mass galaxies. Theory suggests that this tight mass-metallicity relationship is shaped by galactic outflows removing metal-enriched gas from galaxies. Here, we observationally model the outflow metallicities of the warm outflowing phase from a sample of seven local star-forming galaxies with stellar masses between 10$^{7}$-10$^{11}$ M$_\odot$. We estimate the outflow metallicities using four weak rest-frame ultraviolet absorption lines, the observed stellar continua, and photoionization models. The outflow metallicity is flat with $M_\ast$, with a median metallicity of $1.0\pm0.6$ Z$_\odot$. The observed outflows are metal-enriched: low and high-mass galaxies have outflow metallicities 10-50 and 2.6 times larger than their ISM metallicities, respectively. The observed outflows are mainly composed of entrained ISM gas with at most 22% of the metals directly coming from recent supernovae enrichment. The metal outflow rate shallowly increases with $M_\ast$, as $M_\ast^{0.2 \pm 0.1}$, because the mass outflow rate shallow increases with $M_\ast$. Finally, we normalize the metal outflow rate by the rate at which star formation retains metals to calculate the metal-loading factor. The metal-loading factor inversely scales with $M_\ast$. The normalization and scaling of the metal-loading factor agree with analytic expressions that reproduce observed mass-metallicity relations. Galactic outflows fundamentally shape the observed mass-metallicity relationship.

We have obtained the first complete ultraviolet (UV) spectrum of a strong Lyman continuum(LyC) emitter at low redshift -- the compact, low-metallicity, star-forming galaxy J1154+2443 -- with a Lyman continuum escape fraction of 46% discovered recently. The Space Telescope Imaging Spectrograph spectrum shows strong Lya and CIII] 1909 emission, as well as OIII] 1666. Our observations show that strong LyC emitters can have UV emission lines with a high equivalent width (e.g. EW(CIII])$=11.7 \pm 2.9 \AA$ rest-frame), although their equivalent widths should be reduced due to the loss of ionizing photons. The intrinsic ionizing photon production efficiency of J1154+2443 is high, $\log(\xi_{\rm ion}^0)=25.56$ erg$^{-1}$ Hz, comparable to that of other recently discovered $z \sim 0.3-0.4$ LyC emitters. Combining our measurements and earlier determinations from the literature, we find a trend of increasing $\xi_{\rm ion}^0$ with increasing CIII] 1909 equivalent width, which can be understood by a combination of decreasing stellar population age and metallicity. Simple ionization and density-bounded photoionization models can explain the main observational features including the UV spectrum of J1154+2443.

The processes allowing the escape of ionizing photons from galaxies into the intergalactic medium are poorly known. To understand how Lyman continuum (LyC) photons escape galaxies, we constrain the HI covering fractions and column densities using ultraviolet HI and metal absorption lines of 18 star-forming galaxies which have Lyman series observations. Nine of these galaxies are confirmed LyC emitters. We fit the stellar continuum, dust attenuation, metal, and HI properties to consistently determine the UV attenuation, as well as the column densities and covering factors of neutral hydrogen and metals. We use synthetic interstellar absorption lines to explore the systematics of our measurements. Then we apply our method to the observed UV spectra of low-redshift and z-2 galaxies. The observed HI lines are found to be saturated in all galaxies. An indirect approach using OI column densities and the observed O/H abundances yields HI column densities of 18.6 to 20 cm-2. These columns are too high to allow the escape of ionizing photons. We find that the known LyC leakers have HI covering fractions less than unity. Ionizing photons escape through optically thin holes/channels in a clumpy interstellar medium. Our simulations confirm that the HI covering fractions are accurately recovered. The SiII and HI covering fractions scale linearly, in agreement with observations from stacked Lyman break galaxy spectra at z-3. Thus, with an empirical correction, the SiII absorption lines can also be used to determine the HI coverage. Finally, we show that a consistent fitting of dust attenuation, continuum and absorption lines is required to properly infer the covering fraction of neutral gas and subsequently to infer the escape fraction of ionizing radiation. These measurements can estimate the LyC escape fraction, as we demonstrate in a companion paper.

The fraction of ionizing photons that escape high-redshift galaxies sensitively determines whether galaxies reionized the early universe. However, this escape fraction cannot be measured from high-redshift galaxies because the opacity of the intergalactic medium is large at high redshifts. Without methods to indirectly measure the escape fraction of high-redshift galaxies, it is unlikely that we will know what reionized the universe. Here, we analyze the far-ultraviolet (UV) H I (Lyman series) and low-ionization metal absorption lines of nine low-redshift, confirmed Lyman continuum emitting galaxies. We use the H I covering fractions, column densities, and dust attenuations measured in a companion paper to predict the escape fraction of ionizing photons. We find good agreement between the predicted and observed Lyman continuum escape fractions (within $1.4\sigma$) using both the H I and ISM absorption lines. The ionizing photons escape through holes in the H I, but we show that dust attenuation reduces the fraction of photons that escape galaxies. This means that the average high-redshift galaxy likely emits more ionizing photons than low-redshift galaxies. Two other indirect methods accurately predict the escape fractions: the Ly$\alpha$ escape fraction and the optical [O III]/[O II] flux ratio. We use these indirect methods to predict the escape fraction of a sample of 21 galaxies with rest-frame UV spectra but without Lyman continuum observations. Many of these galaxies have low escape fractions ($f_{\rm esc} \le 1$\%), but 11 have escape fractions $>1$\%. The methods presented here will measure the escape fractions of high-redshift galaxies, enabling future telescopes to determine whether star-forming galaxies reionized the early universe.

We introduce Project MEGaSaURA: The Magellan Evolution of Galaxies Spectroscopic and Ultraviolet Reference Atlas. MEGaSaURA comprises medium-resolution, rest-frame ultraviolet spectroscopy of N=15 bright gravitationally lensed galaxies at redshifts of 1.68$<$z$<$3.6, obtained with the MagE spectrograph on the Magellan telescopes. The spectra cover the observed-frame wavelength range $3200 < \lambda_o < 8280$ \AA ; the average spectral resolving power is R=3300. The median spectrum has a signal-to-noise ratio of $SNR=21$ per resolution element at 5000 \AA . As such, the MEGaSaURA spectra have superior signal-to-noise-ratio and wavelength coverage compared to what COS/HST provides for starburst galaxies in the local universe. This paper describes the sample, the observations, and the data reduction. We compare the measured redshifts for the stars, the ionized gas as traced by nebular lines, and the neutral gas as traced by absorption lines; we find the expected bulk outflow of the neutral gas, and no systemic offset between the redshifts measured from nebular lines and the redshifts measured from the stellar continuum. We provide the MEGaSaURA spectra to the astronomical community through a data release.

Galactic outflows regulate the amount of gas galaxies convert into stars. However, it is difficult to measure the mass outflows remove because they span a large range of temperatures and phases. Here, we study the rest-frame ultraviolet spectrum of a lensed galaxy at z~2.9 with prominent interstellar absorption lines from O I, tracing neutral gas, up to O VI, tracing transitional phase gas. The O VI profile mimics weak low-ionization profiles at low velocities, and strong saturated profiles at high velocities. These trends indicate that O VI gas is co-spatial with the low-ionization gas. Further, at velocities blueward of -200 km/s the column density of the low-ionization outflow rapidly drops while the O VI column density rises, suggesting that O VI is created as the low-ionization gas is destroyed. Photoionization models do not reproduce the observed O VI, but adequately match the low-ionization gas, indicating that the phases have different formation mechanisms. Photoionized outflows are more massive than O VI outflows for most of the observed velocities, although the O VI mass outflow rate exceeds the photoionized outflow at velocities above the galaxy's escape velocity. Therefore, most gas capable of escaping the galaxy is in a hot outflow phase. We suggest that the O VI absorption is a temporary by-product of conduction transferring mass from the photoionized phase to an unobserved hot wind, and discuss how this mass-loading impacts the observed circum-galactic medium.

We stack the rest-frame ultraviolet spectra of N=14 highly magnified gravitationally lensed galaxies at redshifts 1.6<z<3.6. The resulting new composite spans $900< \lambda_{rest} < 3000$ \AA, with a peak signal-to-noise ratio of 103 per spectral resolution element ($\sim$100 km/s). It is the highest signal-to-noise ratio, highest spectral resolution composite spectrum of $z\sim2$--3 galaxies yet published. The composite reveals numerous weak nebular emission lines and stellar photospheric absorption lines that can serve as new physical diagnostics, particularly at high redshift with the James Webb Space Telescope (JWST). We report equivalent widths to aid in proposing for and interpreting JWST spectra. We examine the velocity profiles of strong absorption features in the composite, and in a matched composite of $z\sim0$ COS/HST galaxy spectra. We find remarkable similarity in the velocity profiles at $z\sim 0$ and $z\sim2$, suggesting that similar physical processes control the outflows across cosmic time. While the maximum outflow velocity depends strongly on ionization potential, the absorption-weighted mean velocity does not. As such, the bulk of the high-ionization absorption traces the low-ionization gas, with an additional blueshifted absorption tail extending to at least $-2000$ km/s . We interpret this tail as arising from the stellar wind and photospheres of massive stars. Starburst99 models are able to replicate this high-velocity absorption tail. However, these theoretical models poorly reproduce several of the photospheric absorption features, indicating that improvements are needed to match observational constraints on the massive stellar content of star-forming galaxies at $z \sim 2$. We publicly release our composite spectra.

To reionize the early universe, high-energy photons must escape the galaxies that produce them. It has been suggested that stellar feedback drives galactic outflows out of star-forming regions, creating low density channels through which ionizing photons escape into the inter-galactic medium. We compare the galactic outflow properties of confirmed Lyman continuum (LyC) leaking galaxies to a control sample of nearby star-forming galaxies to explore whether the outflows from leakers are extreme as compared to the control sample. We use data from the Cosmic Origins Spectrograph on the Hubble Space Telescope to measure the equivalent widths and velocities of Si II and Si III absorption lines, tracing neutral and ionized galactic outflows. We find that the Si II and Si III equivalent widths of the LyC leakers reside on the low-end of the trend established by the control sample. The leakers' velocities are not statistically different than the control sample, but their absorption line profiles have a different asymmetry: their central velocities are closer to their maximum velocities. The outflow kinematics and equivalent widths are consistent with the scaling relations between outflow properties and host galaxy properties -- most notably metallicity -- defined by the control sample. Additionally, we use the Ly\alpha\ profiles to show that the Si II equivalent width scales with the Ly\alpha\ peak velocity separation. We determine that the low equivalent widths of the leakers are likely driven by low metallicities and low H I column densities, consistent with a density-bounded ionization region, although we cannot rule out significant variations in covering fraction. While we do not find that the LyC leakers have extreme outflow velocities, the low maximum-to-central velocity ratios demonstrate the importance of the acceleration and density profiles for LyC and Ly\alpha\ escape. [abridged]

Galactic outflows are believed to play an important role in regulating star formation in galaxies, but estimates of the outflowing mass and momentum have historically been based on uncertain assumptions. Here, we measure the mass, momentum, and energy outflow rates of seven nearby star-forming galaxies using ultraviolet absorption lines and observationally motivated estimates for the density, metallicity, and radius of the outflow. Low-mass galaxies generate outflows faster than their escape velocities with mass outflow rates up to twenty times larger than their star formation rates. These outflows from low-mass galaxies also have momenta larger than provided from supernovae alone, indicating that multiple momentum sources drive these outflows. Only 1-20\% of the supernovae energy is converted into kinetic energy, and this fraction decreases with increasing stellar mass such that low-mass galaxies drive more efficient outflows. We find scaling relations between the outflows and the stellar mass of their host galaxies (M$_\ast$) at the 2-3$\sigma$ significance level. The mass-loading factor, or the mass outflow rate divided by the star formation rate, scales as M$_\ast^{-0.4}$ and with the circular velocity as v$_\mathrm{circ}^{-1.6}$. The scaling of the mass-loading factor is similar to recent simulations, but the observations are a factor of five smaller, possibly indicating that there is a substantial amount of unprobed gas in a different ionization phase. The outflow momenta are consistent with a model where star formation drives the outflow while gravity counteracts this acceleration.

We study the properties of 66 galaxies with kinematically misaligned gas and stars from MaNGA survey. The fraction of kinematically misaligned galaxies varies with galaxy physical parameters, i.e. M*, SFR and sSFR. According to their sSFR, we further classify these 66 galaxies into three categories, 10 star-forming, 26 "Green Valley" and 30 quiescent ones. The properties of different types of kinematically misaligned galaxies are different in that the star-forming ones have positive gradient in D4000 and higher gas-phase metallicity, while the green valley/quiescent ones have negative D4000 gradients and lower gas-phase metallicity on average. There is evidence that all types of the kinematically misaligned galaxies tend to live in more isolated environment. Based on all these observational results, we propose a scenario for the formation of star forming galaxies with kinematically misaligned gas and stars - the progenitor accretes misaligned gas from a gas-rich dwarf or cosmic web, the cancellation of angular momentum from gas-gas collisions between the pre-existing gas and the accreted gas largely accelerates gas inflow, leading to fast centrally-concentrated star-formation. The higher metallicity is due to enrichment from this star formation. For the kinematically misaligned green valley and quiescent galaxies, they might be formed through gas-poor progenitors accreting kinematically misaligned gas from satellites which are smaller in mass.

Galaxies grow through both internal and external processes. In about 10% of nearby red galaxies with little star formation, gas and stars are counter-rotating, demonstrating the importance of external gas acquisition in these galaxies. However, systematic studies of such phenomena in blue, star-forming galaxies are rare, leaving uncertain the role of external gas acquisition in driving evolution of blue galaxies. Based on new measurements with integral field spectroscopy of a large representative galaxy sample, we find an appreciable fraction of counter-rotators among blue galaxies (9 out of 489 galaxies). The central regions of blue counter-rotators show younger stellar populations and more intense, ongoing star formation than their outer parts, indicating ongoing growth of the central regions. The result offers observational evidence that the acquisition of external gas in blue galaxies is possible; the interaction with pre-existing gas funnels the gas into nuclear regions (< 1 kpc) to form new stars.

Baryons cycle into galaxies from the inter-galactic medium, are converted into stars, and a fraction of the baryons are ejected out of galaxies by stellar feedback. Here we present new high resolution (3.9"; 68 pc) CO(2-1) and CO(3-2) images that probe these three stages of the baryon cycle in the nearby starburst M 82. We combine these new observations with previous CO(1-0) and [Fe II] images to study the physical conditions within the molecular gas. Using a Bayesian analysis and the radiative transfer code RADEX, we model molecular Hydrogen temperatures and densities, as well as CO column densities. Besides the disc, we concentrate on two regions within the galaxy: an expanding super-bubble and the base of a molecular streamer. Shock diagnostics, kinematics, and optical extinction suggest that the streamer is an inflowing filament, with a molecular gas mass inflow rate of 3.5 M$_\odot$ yr$^{-1}$. We measure the molecular gas mass outflow rate of the expanding super-bubble to be 17 M$_\odot$ yr$^{-1}$, 5 times higher than the inferred inflow rate, and 1.3 times the star formation rate of the galaxy. The high mass outflow rate and large star formation rate will deplete the galaxy of molecular gas within eight million years, unless there are additional sources of molecular gas.

To evaluate the impact of stellar feedback, it is critical to estimate the mass outflow rates of galaxies. Past estimates have been plagued by uncertain assumptions about the outflow geometry, metallicity, and ionization fraction. Here we use Hubble Space Telescope ultraviolet spectroscopic observations of the nearby starburst NGC 6090 to demonstrate that many of these quantities can be constrained by the data. We use the Si~{\sc IV} absorption lines to calculate the scaling of velocity (v), covering fraction (C$_f$), and density with distance from the starburst (r), assuming the Sobolev optical depth and a velocity law of the form: v$~\propto(1 -\mathrm{R}_\mathrm{i}/\mathrm{r} )^\beta$ (where R$_\mathrm{i}$ is the inner outflow radius). We find that the velocity ($\beta$=0.43) is consistent with an outflow driven by an r$^{-2}$ force with the outflow radially accelerated, while the scaling of the covering fraction ($C_f \propto \mathrm{r}^{-0.82}$) suggests that cool clouds in the outflow are in pressure equilibrium with an adiabatically expanding medium. We use the column densities of four weak metal lines and CLOUDY photoionization models to determine the outflow metallicity, the ionization correction, and the initial density of the outflow. Combining these values with the profile fitting, we find R$_\mathrm{i}$ = 63 pc, with most of the mass within 300~pc of the starburst. Finally, we find that the maximum mass outflow rate is 2.3~M$_\odot$ yr$^{-1}$ and the mass loading factor (outflow divided by the star formation rate) is 0.09, a factor of 10 lower than the value calculated using common assumptions for the geometry, metallicity and ionization structure of the outflow.

We study the ionization structure of galactic outflows in 37 nearby, star forming galaxies with the Cosmic Origins Spectrograph on the Hubble Space Telescope. We use the O I, Si II, Si III, and Si IV ultraviolet absorption lines to characterize the different ionization states of outflowing gas. We measure the equivalent widths, line widths, and outflow velocities of the four transitions, and find shallow scaling relations between them and galactic stellar mass and star formation rate. Regardless of the ionization potential, lines of similar strength have similar velocities and line widths, indicating that the four transitions can be modeled as a co-moving phase. The Si equivalent width ratios (e.g. Si IV/Si II) have low dispersion, and little variation with stellar mass; while ratios with O I and Si vary by a factor of 2 for a given stellar mass. Photo-ionization models reproduce these equivalent width ratios, while shock models under predict the relative amount of high ionization gas. The photo-ionization models constrain the ionization parameter (U) between -2.25 < log(U) < -1.5, and require that the outflow metallicities are greater than 0.5 Z$_\odot$. We derive ionization fractions for the transitions, and show that the range of ionization parameters and stellar metallicities leads to a factor of 1.15-10 variation in the ionization fractions. Historically, mass outflow rates are calculated by converting a column density measurement from a single metal ion into a total Hydrogen column density using an ionization fraction, thus mass outflow rates are sensitive to the assumed ionization structure of the outflow.

We investigate the supernova-driven galactic wind of the barred spiral galaxy NGC 7552, using both ground-based optical nebular emission lines and far-ultraviolet absorption lines measured with the Hubble Space Telescope Cosmic Origins Spectrograph. We detect broad (~300 km/s) blueshifted (-40 km/s) optical emission lines associated with the galaxy's kpc-scale star-forming ring. The broad line kinematics and diagnostic line ratios suggest that the H-alpha emission comes from clouds of high density gas entrained in a turbulent outflow. We compare the H-alpha emission line profile to the UV absorption line profile measured along a coincident sight line and find significant differences. The maximum blueshift of the H-alpha-emitting gas is ~290 km/s, whereas the UV line profile extends to blueshifts upwards of 1000 km/s. The mass outflow rate estimated from the UV is roughly nine times greater than that estimated from H-alpha. We argue that the H-alpha emission traces a cluster-scale outflow of dense, low velocity gas at the base of the large-scale wind. We suggest that UV absorption line measurements are therefore more reliable tracers of warm gas in starburst-driven outflows.

We report on a sample of 51 nearby, star-forming galaxies observed with the Cosmic Origin Spectrograph on the Hubble Space Telescope. We calculate Si II kinematics and densities arising from warm gas entrained in galactic outflows. We use multi-wavelength ancillary data to estimate stellar masses (M$_\ast$), star-formation rates (SFR), and morphologies. We derive significant correlations between outflow velocity and SFR$^{\sim 0.1}$, M$_\ast^{\sim 0.1}$ and v$_\text{circ}^{\sim 1/2}$. Some mergers drive outflows faster than these relations prescribe, launching the outflow faster than the escape velocity. Calculations of the mass outflow rate reveal strong scaling with SFR$^{\sim 1/2}$ and M$_\ast^{\sim 1/2}$. Additionally, mass-loading efficiency factors (mass outflow rate divided by SFR) scale approximately as M$_\ast^{-1/2}$. Both the outflow velocity and mass-loading scaling suggest that these outflows are powered by supernovae, with only 0.7% of the total supernovae energy converted into the kinetic energy of the warm outflow. Galaxies lose some gas if log(M$_\ast$/M$_\odot$) < $9.5$, while more massive galaxies retain all of their gas, unless they undergo a merger. This threshold for gas loss can explain the observed shape of the mass-metallicity relation.

We use principal component analysis (PCA) to estimate stellar masses, mean stellar ages, star formation histories (SFHs), dust extinctions and stellar velocity dispersions for ~290,000 galaxies with stellar masses greater than $10^{11}Msun and redshifts in the range 0.4<z<0.7 from the Baryon Oscillation Spectroscopic Survey (BOSS). We find the fraction of galaxies with active star formation first declines with increasing stellar mass, but then flattens above a stellar mass of 10^{11.5}Msun at z~0.6. This is in striking contrast to z~0.1, where the fraction of galaxies with active star formation declines monotonically with stellar mass. At stellar masses of 10^{12}Msun, therefore, the evolution in the fraction of star-forming galaxies from z~0.6 to the present-day reaches a factor of ~10. When we stack the spectra of the most massive, star-forming galaxies at z~0.6, we find that half of their [OIII] emission is produced by AGNs. The black holes in these galaxies are accreting on average at ~0.01 the Eddington rate. To obtain these results, we use the stellar population synthesis models of Bruzual & Charlot (2003) to generate a library of model spectra with a broad range of SFHs, metallicities, dust extinctions and stellar velocity dispersions. The PCA is run on this library to identify its principal components over the rest-frame wavelength range 3700-5500A. We demonstrate that linear combinations of these components can recover information equivalent to traditional spectral indices such as the 4000A break strength and HdA, with greatly improved S/N. This method is able to recover physical parameters such as stellar mass-to-light ratio, mean stellar age, velocity dispersion and dust extinction from the relatively low S/N BOSS spectra. We examine the sensitivity of our stellar mass estimates to the input parameters in our model library and the different stellar population synthesis models.

Primordial Black Holes (PBHs) that form from the collapse of density perturbations are more clustered than the underlying density field. In a previous paper, we showed the constraints that this has on the prospects of PBH dark matter. In this paper we examine another consequence of this clustering: the formation of bound systems of PBHs in the early universe. These would hypothetically be the earliest gravitationally collapsed structures, forming when the universe is still radiation dominated. Depending upon the size and occupation of the clusters, PBH merging occurs before they would have otherwise evaporated due to Hawking evaporation.

In this paper we revisit thermodynamic constraints on primordial black hole (PBH) formation in the early universe. Under the assumption that PBH mass is equal to the cosmological horizon mass, one can use the 2nd Law of Thermodynamics to put a lower limit on the PBH mass. In models of PBH formation, however, PBHs are created at some fraction of the horizon mass. We show that this thermodynamic constraint still holds for sub-horizon PBH formation.

We investigate the spatial clustering properties of primordial black holes (PBHs). With minimal assumptions, we show that PBHs created in the radiation era are highly clustered. Using the peaks theory model of bias, we compute the PBH two-point correlation function and power spectrum. For creation from an initially adiabatic power spectrum of perturbations, the PBH power spectrum contains both isocurvature and adiabatic components. The absence of observed isocurvature fluctuations today constrains the mass range in which PBHs may serve as dark matter. We briefly discuss other consequences of PBH clustering.

The violation of Lorentz invariance (LI) has been invoked in a number of ways to explain issues dealing with ultrahigh energy cosmic ray (UHECR) production and propagation. These treatments, however, have mostly been limited to examples in the proton-neutron system and photon-electron system. In this paper we show how a broader violation of Lorentz invariance would allow for a series of previously forbidden decays to occur, and how that could lead to UHECR primaries being heavy baryonic states or Higgs bosons.

We search for nearby, isolated, accreting, ``stellar-mass'' (3 to $100M_\odot$) black holes. Models suggest a synchrotron spectrum in visible wavelengths and some emission in X-ray wavelengths. Of 3.7 million objects in the Sloan Digital Sky Survey Early Data Release, about 150,000 objects have colors and properties consistent with such a spectrum, and 87 of these objects are X-ray sources from the ROSAT All Sky Survey. Thirty-two of these have been confirmed not to be black-holes using optical spectra. We give the positions and colors of these 55 black-hole candidates, and quantitatively rank them on their likelihood to be black holes. We discuss uncertainties the expected number of sources, and the contribution of blackholes to local dark matter.