We predict the number of black holes with stellar companions that are potentially detectable with Gaia astrometry over the course of its five-year mission. Our model estimates that nearly 2$\times$10$^5$ astrometric binaries hosting black holes and stellar companions brighter than Gaia's detection threshold, G $\sim$ 20, should be discovered with 5$\sigma$ sensitivity. Among these detectable binaries, systems with longer orbital periods are favored, and black hole and stellar companion masses in the range M$_{BH} \sim$ 6-10 M$_\odot$ and M$_* \sim$ 1-2 M$_\odot$, respectively, are expected to dominate.

We explore the possibility of planet formation in the carbon-rich protoplanetary disks of carbon-enhanced metal-poor (CEMP) stars, possible relics of the early Universe. The chemically anomalous abundance patterns ([C/Fe] $\geq$ 0.7) in this subset of low-mass stars suggest pollution by primordial core-collapsing supernovae (SNe) ejecta that are particularly rich in carbon dust grains. By comparing the dust-settling timescale in the protoplanetary disks of CEMP stars to the expected disk lifetime (assuming dissipation via photoevaporation), we determine the maximum distance $r_{max}$ from the host CEMP star at which carbon-rich planetesimal formation is possible, as a function of the host star's [C/H] abundance. We then use our linear relation between $r_{max}$ and [C/H], along with the theoretical mass-radius relation derived for a solid, pure carbon planet, to characterize potential planetary transits across host CEMP stars. Given that the related transits are detectable with current and upcoming space-based transit surveys, we suggest initiating an observational program to search for carbon planets around CEMP stars in hopes of shedding light on the question of how early planetary systems may have formed after the Big Bang.

We show that the cumulative CO emission from galaxies throughout cosmic history distorts the spectrum of the cosmic microwave background (CMB) at a level that is well above the detection limit of future instruments, such as the Primordial Inflation Explorer (PIXIE). The modeled CO signal has a prominent bump in the frequency interval 100-200 GHz, with a characteristic peak intensity of ~ 2$\times$10$^{-23}$ W m$^{-2}$ Hz$^{-1}$ sr$^{-1}$. Most of the CO foreground originates from modest redshifts, z ~ 2-5, and needs to be efficiently removed for more subtle distortions from the earlier universe to be detected.

Using the most recent measurements of the ultraviolet (UV) luminosity functions (LFs) and dust estimates of early galaxies, we derive updated dust-corrected star-formation rate functions (SFRFs) at z~4-8, which we model to predict the evolution to higher redshifts, z>8. We employ abundance matching techniques to calibrate a relation between galaxy star formation rate (SFR) and host halo mass M$_{h}$ by mapping the shape of the observed SFRFs at z~4-8 to that of the halo mass function. The resulting scaling law remains roughly constant over this redshift range. We apply the average SFR-M$_{h}$ relation to reproduce the observed SFR functions at 4 <= z <= 8 and also derive the expected UV LFs at higher redshifts. At z~9 and z~10 these model LFs are in excellent agreement with current observed estimates. Our predicted number densities and UV LFs at z>10 indicate that JWST will be able to detect galaxies out to z~15 with an extensive treasury sized program. We also derive the redshift evolution of the star formation rate density and associated reionization history by galaxies. Models which integrate down to the current HUDF12/XDF detection limit (M$_{UV}$ ~ -17.7 mag) result in a SFRD that declines as (1+z)$^{-10.4\pm0.3}$ at high redshift and fail to reproduce the observed CMB electron scattering optical depth, $\tau$ ~ 0.066, to within 1$\sigma$. On the other hand, we find that the inclusion of galaxies with SFRs well below the current detection limit (M$_{UV}$ < -5.7 mag) leads to a fully reionized universe by z ~ 6.5 and an optical depth of $\tau$ ~ 0.054, consistent with the recently derived Planck value at the 1$\sigma$ level.

We present a novel approach to estimating the intensity mapping signal of any CO rotational line emitted during the Epoch of Reionization (EoR). Our approach is based on large velocity gradient (LVG) modeling, a radiative transfer modeling technique that generates the full CO spectral line energy distribution (SLED) for a specified gas kinetic temperature, volume density, velocity gradient, molecular abundance, and column density. These parameters, which drive the physics of CO transitions and ultimately dictate the shape and amplitude of the CO SLED, can be linked to the global properties of the host galaxy, mainly the star formation rate (SFR) and the SFR surface density. By further employing an empirically derived SFR-M relation for high redshift galaxies, we can express the LVG parameters, and thus the specific intensity of any CO rotational transition, as functions of the host halo mass M and redshift z. Integrating over the range of halo masses expected to host CO-luminous galaxies, we predict a mean CO(1-0) brightness temperature ranging from ~0.6 {\mu}K at z= 6 to ~0.03 {\mu} at z= 10 with brightness temperature fluctuations of $\Delta_{CO}^2$ ~ 0.1 and 0.005 {\mu}K respectively, at k = 0.1 Mpc$^{-1}$. In this model, the CO emission signal remains strong for higher rotational levels at z = 6, with <T$_{CO}$>~ 0.3 and 0.05 {\mu}K for the CO J = 6->5 and CO J = 10->9 transitions respectively. Including the effects of CO photodissociation in these molecular clouds, especially at low metallicities, results in the overall reduction in the amplitude of the CO signal, with the low- and high-J lines weakening by 2-20% and 10-45%, respectively, over the redshift range 4 < z < 10.

We report the detection of far-infrared (FIR) CO rotational emission from nearby active galactic nuclei (AGN) and starburst galaxies, as well as several merging systems and Ultra-Luminous Infrared Galaxies (ULIRGs). Using Herschel-PACS, we have detected transitions in the J$_{upp}$ = 14 - 20 range ($\lambda \sim$ 130 - 185 $\mu$m, $\nu \sim$ 1612 - 2300 GHz) with upper limits on (and in two cases, detections of) CO line fluxes up to J$_{upp}$ = 30. The PACS CO data obtained here provide the first well-sampled FIR extragalactic CO SLEDs for this range, and will be an essential reference for future high redshift studies. We find a large range in the overall SLED shape, even amongst galaxies of similar type, demonstrating the uncertainties in relying solely on high-J CO diagnostics to characterize the excitation source of a galaxy. Combining our data with low-J line intensities taken from the literature, we present a CO ratio-ratio diagram and discuss its potential diagnostic value in distinguishing excitation sources and physical properties of the molecular gas. The position of a galaxy on such a diagram is less a signature of its excitation mechanism, than an indicator of the presence (or absence) of warm, dense molecular gas. We then quantitatively analyze the CO emission from a subset of the detected sources with Large Velocity Gradient (LVG) radiative transfer models to fit the CO SLEDs. Using both single-component and two-component LVG models to fit the kinetic temperature, velocity gradient, number density and column density of the gas, we derive the molecular gas mass and the corresponding CO-to-H$_2$ conversion factor, $\alpha_{CO}$, for each respective source. For the ULIRGs we find $\alpha$ values in the canonical range 0.4 - 5 M$_\odot$/(K kms$^{-1}$pc$^2$), while for the other objects, $\alpha$ varies between 0.2 and 14.} Finally, we compare our best-fit LVG model ..

We simulate the effects of gravitational lensing on the source count of high redshift galaxies as projected to be observed by the Hubble Frontier Fields program and the James Webb Space Telescope (JWST) in the near future. Taking the mass density profile of the lensing object to be the singular isothermal sphere (SIS) or the Navarro-Frenk-White (NFW) profile, we model a lens residing at a redshift of z_L = 0.5 and explore the radial dependence of the resulting magnification bias and its variability with the velocity dispersion of the lens, the photometric sensitivity of the instrument, the redshift of the background source population, and the intrinsic maximum absolute magnitude (M_{max}) of the sources. We find that gravitational lensing enhances the number of galaxies with redshifts z >= 13 detected in the angular region \theta_E/2 <= \theta <= 2\theta_E (where \theta_E is the Einstein angle) by a factor of ~ 3 and 1.5 in the HUDF (df/d\nu_0 ~ 9 nJy) and medium-deep JWST surveys (df/d\nu_0 ~ 6 nJy). Furthermore, we find that even in cases where a negative magnification bias reduces the observed number count of background sources, the lensing effect improves the sensitivity of the count to the intrinsic faint-magnitude cut-off of the Schechter luminosity function. In a field centered on a strong lensing cluster, observations of z >= 6 and z >= 13 galaxies with JWST can be used to infer this cut-off magnitude for values as faint as M_{max} ~ -14.4 and -16.1 mag (L_{min} ~ 2.5*10^{26} and 1.2*10^{27} erg s^{-1} Hz^{-1}) respectively, within the range bracketed by existing theoretical models. Gravitational lensing may therefore offer an effective way of constraining the low-luminosity cut-off of high-redshift galaxies.

Walter et al. (20012) have recently identified the J=6-5, 5-4, and 2-1 CO rotational emission lines, and [C_{II}] fine-structure emission line from the star-forming interstellar medium in the high-redshift submillimeter source HDF 850.1, at z = 5.183. We employ large velocity gradient (LVG) modeling to analyze the spectra of this source assuming the [C_{II}] and CO emissions originate from (i) separate unvirialized regions, (ii) separate virialized regions, (iii) uniformly mixed unvirialized region, and (iv) uniformly mixed virialized regions. We present the best fit set of parameters, including for each case the ratio $\alpha$ between the total hydrogen/helium gas mass and the CO(1-0) line luminosity. We also present computations of the ratio of H_{2} mass to [C_{II}] line-luminosity for optically thin conditions, for a range of gas temperatures and densities, for direct conversion of [C_{II}] line-luminosities to "dark-H_{2}" masses. For HDF 850.1 we find that a model in which the CO and C^{+} are uniformly mixed in gas that is shielded from UV radiation, requires a cosmic-ray or X-ray ionization rate of $\zeta \approx$ 10^{-13} s^{-1}, plausibly consistent with the large star-formation rate ($\sim$ 10^{3} M$_{\odot}$ yr^{-1}) observed in this source. Enforcing the cosmological constraint posed by the abundance of dark matter halos in the standard $\Lambda$CDM cosmology and taking into account other possible contributions to the total gas mass, we find that three of these four models are less likely at the 2$\sigma$ level. We conclude that modeling HDF 850.1's ISM as a collection of unvirialized molecular clouds with distinct CO and C^{+} layers, for which $\alpha$ = 0.6 M$_{\odot}$ (K km s^{-1} pc^{2})^{-1} for the CO to H_{2} mass-to-luminosity ratio, (similar to the standard ULIRG value), is most consistent with the $\Lambda$CDM cosmology.