We present a fast implementation of the next-to-leading order (1-loop) redshift-space galaxy power spectrum by using FFTlog-based methods. [V. Desjacques, D. Jeong, and F. Schmidt, JCAP 1812 (12), 035] have shown that the 1-loop galaxy power spectrum in redshift space can be computed with 28 independent loop integrals with 22 bias parameters. Analytical calculation of the angular part of the loop integrals leaves the radial part in the form of a spherical Bessel transformation that is ready to be integrated numerically by using the FFTLog transformation. We find that the original 28 loop integrals can be solved with a total of 85 unique FFTLog transformations, yet leading to a few orders of magnitude speed up over traditional multi-dimensional integration. The code used in this work is publicly available at https://github.com/JosephTomlinson/GeneralBiasPk

The galaxy catalogs generated from low-resolution emission line surveys often contain both foreground and background interlopers due to line misidentification, which can bias the cosmological parameter estimation. In this paper, we present a method for correcting the interloper bias by using the joint-analysis of auto- and cross-power spectra of the main and the interloper samples. In particular, we can measure the interloper fractions from the cross-correlation between the interlopers and survey galaxies, because the true cross-correlation must be negligibly small. The estimated interloper fractions, in turn, remove the interloper bias in the cosmological parameter estimation. For example, in the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) low-redshift ($z<0.5$) [O II] $\lambda3727${\AA} emitters contaminate high-redshift ($1.9<z<3.5$) Lyman-$\alpha$ line emitters. We demonstrate that the joint-analysis method yields a high signal-to-noise ratio measurement of the interloper fractions while only marginally increasing the uncertainties in the cosmological parameters relative to the case without interlopers. We also show the same is true for the high-latitude spectroscopic survey of Wide-Field Infrared Survey Telescope (WFIRST) mission where contamination occurs between the Balmer-$\alpha$ line emitters at lower redshifts ($1.1<z<1.9$) and Oxygen ([O III] $\lambda5007${\AA}) line emitters at higher redshifts ($1.7<z<2.8$).

We explore the possible spectrum of binary mergers of sub-solar mass black holes formed out of dark matter particles interacting via a dark electromagnetism. We estimate the properties of these dark black holes by assuming that their formation process is parallel to Population-III star formation; except that dark molecular cooling can yield smaller opacity limit. We estimate the binary coalescence rates for the Advanced LIGO and Einstein telescope, and find that scenarios compatible with all current constraints could produce dark black holes at rates high enough for detection by Advanced LIGO.

We present the 2-point function from Fast and Accurate Spherical Bessel Transformation (2-FAST) algorithm for a fast and accurate computation of integrals involving one or two spherical Bessel functions. These types of integrals occur when projecting the galaxy power spectrum $P(k)$ onto the configuration space, $\xi_\ell^\nu(r)$, or spherical harmonic space, $C_\ell(\chi,\chi')$. First, we employ the FFTlog transformation of the power spectrum to divide the calculation into $P(k)$-dependent coefficients and $P(k)$-independent integrations of basis functions multiplied by spherical Bessel functions. We find analytical expressions for the latter integrals in terms of special functions, for which recursion provides a fast and accurate evaluation. The algorithm, therefore, circumvents direct integration of highly oscillating spherical Bessel functions.

We compare the rest-frame ultraviolet and rest-frame optical morphologies of 2 < z < 3 star-forming galaxies in the GOODS-S field using Hubble Space Telescope WFC3 and ACS images from the CANDELS, GOODS, and ERS programs. We show that the distribution of sizes and concentrations for 1.90 < z < 2.35 galaxies selected via their rest-frame optical emission-lines are statistically indistinguishable from those of Lyman-alpha emitting systems found at z ~ 2.1 and z ~ 3.1. We also show that the z > 2 star-forming systems of all sizes and masses become smaller and more compact as one shifts the observing window from the UV to the optical. We argue that this offset is due to inside-out galaxy formation over the first ~ 2 Gyr of cosmic time.

We compare the physical and morphological properties of z ~ 2 Lyman-alpha emitting galaxies (LAEs) identified in the HETDEX Pilot Survey and narrow band studies with those of z ~ 2 optical emission line galaxies (oELGs) identified via HST WFC3 infrared grism spectroscopy. Both sets of galaxies extend over the same range in stellar mass (7.5 < logM < 10.5), size (0.5 < R < 3.0 kpc), and star-formation rate (~1 < SFR < 100). Remarkably, a comparison of the most commonly used physical and morphological parameters -- stellar mass, half-light radius, UV slope, star formation rate, ellipticity, nearest neighbor distance, star formation surface density, specific star formation rate, [O III] luminosity, and [O III] equivalent width -- reveals no statistically significant differences between the populations. This suggests that the processes and conditions which regulate the escape of Ly-alpha from a z ~ 2 star-forming galaxy do not depend on these quantities. In particular, the lack of dependence on the UV slope suggests that Ly-alpha emission is not being significantly modulated by diffuse dust in the interstellar medium. We develop a simple model of Ly-alpha emission that connects LAEs to all high-redshift star forming galaxies where the escape of Ly-alpha depends on the sightline through the galaxy. Using this model, we find that mean solid angle for Ly-alpha escape is 2.4+/-0.8 steradians; this value is consistent with those calculated from other studies.

We explore the evolution of the Stellar Mass-Star Formation Rate-Metallicity Relation using a set of 256 COSMOS and GOODS galaxies in the redshift range 1.90 < z < 2.35. We present the galaxies' rest-frame optical emission-line fluxes derived from IR-grism spectroscopy with the Hubble Space Telescope and combine these data with star formation rates and stellar masses obtained from deep, multi-wavelength (rest-frame UV to IR) photometry. We then compare these measurements to those for a local sample of galaxies carefully matched in stellar mass (7.5 < log(M*/Msol) < 10.5) and star formation rate (-0.5 < log(SFR) < 2.5 in Msol yr^-1). We find that the distribution of z ~ 2.1 galaxies in stellar mass-SFR-metallicity space is clearly different from that derived for our sample of similarly bright (L_H\b{eta} > 3 . 10^40 ergs s^-1) local galaxies, and this offset cannot be explained by simple systematic offsets in the derived quantities. At stellar masses above ~10^9 Msol and star formation rates above ~10 Msol yr^-1, the z ~ 2.1 galaxies have higher oxygen abundances than their local counterparts, while the opposite is true for lower-mass, lower-SFR systems.

We derive the mean wavelength dependence of stellar attenuation in a sample of 239 high redshift (1.90 < z < 2.35) galaxies selected via Hubble Space Telescope (HST) WFC3 IR grism observations of their rest-frame optical emission lines. Our analysis indicates that the average reddening law follows a form similar to that derived by Calzetti et al. for local starburst galaxies. However, over the mass range 7.2 < log M/Msolar < 10.2, the slope of the attenuation law in the UV is shallower than that seen locally, and the UV slope steepens as the mass increases. These trends are in qualitative agreement with Kriek & Conroy, who found that the wavelength dependence of attenuation varies with galaxy spectral type. However, we find no evidence of an extinction "bump" at 2175 A in any of the three stellar mass bins, or in the sample as a whole. We quantify the relation between the attenuation curve and stellar mass and discuss its implications.

We use near-IR grism spectroscopy from the Hubble Space Telescope to examine the strength of [Ne~III] 3869 relative to H-beta, [O~II] 3727 and [O~III] 5007 in 236 low mass (7.5 < log (Mstar/Msolar) < 10.5) star-forming galaxies in the redshift range 1.90 < z < 2.35. By stacking the data by stellar mass, we show that the [Ne~III]/[O~II] ratios of the z ~ 2 universe are marginally higher than those seen in a comparable set of local SDSS galaxies, and that [Ne~III]/[O~III] is enhanced by ~0.2 dex. We consider the possible explanations for this ~4-sigma result, including higher oxygen depletion out of the gas-phase, denser H~II regions, higher production of Ne22 via Wolf-Rayet stars, and the existence of a larger population of X-ray obscured AGN at z ~ 2 compared to z ~ 0. None of these simple scenarios, alone, are favored to explain the observed line ratios. We conclude by suggesting several avenues of future observations to further explore the mystery of enhanced [Ne~III] emission.

We compare the H-beta line strengths of 1.90 < z < 2.35 star-forming galaxies observed with the near-IR grism of the Hubble Space Telescope with ground-based measurements of Ly-alpha from the HETDEX Pilot Survey and narrow-band imaging. By examining the line ratios of 73 galaxies, we show that most star-forming systems at this epoch have a Ly-alpha escape fraction below ~6%. We confirm this result by using stellar reddening to estimate the effective logarithmic extinction of the H-beta emission line (c_Hbeta = 0.5) and measuring both the H-beta and Ly-alpha luminosity functions in a ~ 100,000 cubic Mpc volume of space. We show that in our redshift window, the volumetric Ly-alpha escape fraction is at most 4.4+/-2.1(1.2)%, with an additional systematic ~25% uncertainty associated with our estimate of extinction. Finally, we demonstrate that the bulk of the epoch's star-forming galaxies have Ly-alpha emission line optical depths that are significantly greater than that for the underlying UV continuum. In our predominantly [O~III] 5007-selected sample of galaxies, resonant scattering must be important for the escape of Ly-alpha photons.

We use Hubble Space Telescope near-IR grism spectroscopy to examine the H-beta line strengths of 260 star-forming galaxies in the redshift range 1.90 < z < 2.35. We show that at these epochs, the H-beta star formation rate (SFR) is a factor of ~1.8 higher than what would be expected from the systems' rest-frame UV flux density, suggesting a shift in the standard conversion between these quantities and star formation rate. We demonstrate that at least part of this shift can be attributed to metallicity, as H-beta is more greatly enhanced in systems with lower oxygen abundance. This offset must be considered when measuring the star formation rate history of the universe. We also show that the relation between stellar and nebular extinction in our z ~ 2 sample is consistent with that observed in the local universe.

We use broadband photometry extending from the rest-frame UV to the near-IR to fit the individual spectral energy distributions (SEDs) of 63 bright (L(Ly-alpha) > 10^43 ergs/s) Ly-alpha emitting galaxies (LAEs) in the redshift range 1.9 < z < 3.6. We find that these LAEs are quite heterogeneous, with stellar masses that span over three orders of magnitude, from 7.5 < log M < 10.5. Moreover, although most LAEs have small amounts of extinction, some high-mass objects have stellar reddenings as large as E(B-V) ~0.4. Interestingly, in dusty objects the optical depths for Ly-alpha and the UV continuum are always similar, indicating that Ly-alpha photons are not undergoing many scatters before escaping their galaxy. In contrast, the ratio of optical depths in low-reddening systems can vary widely, illustrating the diverse nature of the systems. Finally, we show that in the star formation rate (SFR)-log mass diagram, our LAEs fall above the "main-sequence" defined by z ~ 3 continuum selected star-forming galaxies. In this respect, they are similar to sub-mm-selected galaxies, although most LAEs have much lower mass.