In the coming decade, a new generation of telescopes, including JWST and WFIRST, will probe the period of the formation of first galaxies and quasars, and open up the last frontier for structure formation. Recent simulations as well as observations have suggested that these galaxies are strongly clustered (with large scale bias $\gtrsim6$), and therefore have significant cosmic variance. In this work, we use \texttt{BlueTides}, the largest volume cosmological simulation of galaxy formation, to directly estimate the cosmic variance for current and upcoming surveys. Given its resolution and volume, \texttt{BlueTides} can probe the bias and cosmic variance of $z>7$ galaxies between magnitude $M_{UV}\sim-16$ to $M_{UV}\sim-22$ over survey areas $\sim0.1\ \mathrm{arcmin}^2$ to $\sim 10~\mathrm{deg}^2$. Within this regime, the cosmic variance decreases with survey area/ volume as a power law with exponents between $\sim-0.25$ to $\sim-0.45$. For the planned $10~\mathrm{deg}^2$ field of WFIRST, the cosmic variance is between $3\%$ to $10\%$. Upcoming JWST medium/ deep surveys with areas up to $A\sim100\ \mathrm{arcmin}^2$ will have cosmic variance ranging from $\sim 20-50\%$. Lensed surveys have the highest cosmic variance $\gtrsim 40\%$; the cosmic variance of $M_{UV}\lesssim-16$ galaxies is $\lesssim100\%$ up to $z\sim11$. At higher redshifts such as $z\sim12~(14)$, effective volumes of $\gtrsim(8~\mathrm{Mpc}/h)^3$ ($\gtrsim(12\ \mathrm{Mpc}/h)^3$) are required to limit the cosmic variance to within $100\%$. Finally, we find that cosmic variance is larger than Poisson variance and forms the dominant component of the overall uncertainty in all current and upcoming surveys. We present our calculations in the form of simple fitting functions and an online cosmic variance calculator (CV_AT_COSMIC_DAWN) which we publicly release.

In this study we quantify the properties of the gas and dark matter around active galactic nuclei (AGN) in simulated galaxy groups and clusters and analyze the effect of AGN feedback on the surrounding intra-cluster (group) medium. Our results suggest downsizing of AGN luminosity with host halo mass, supporting the results obtained from clustering studies of AGN. By examining the temperature and density distribution of the gas in the vicinity of AGN we show that due to feedback from the central engine, the gas gets displaced from the centre of the group/cluster resulting in a reduction of the density but an enhancement of temperature. We show that these effects are pronounced at both high and low redshifts and propose new observables to study the effect of feedback in higher redshift galaxies. We also show that the average stellar mass is decreased in halos in the presence of AGN feedback confirming claims from previous studies. Our work for the first time uses a fully cosmological-hydrodynamic simulation to evaluate the global effects of AGN feedback on their host dark matter halos as well as galaxies at scales of galaxy groups and clusters.

We examine recent high-precision measurements of small-scale quasar clustering (at $z\sim0.5-2$ on scales of $\sim25~\mathrm{kpc/h}$) from the SDSS in the context of the MassiveBlackII (MBII) cosmological hydrodynamic simulation and conditional luminosity function (CLF) modeling. At these high luminosities ($g < 20.85$ quasars), the MBII simulation volume ($100~\mathrm{cMpc}/h$ comoving boxsize) has only 3 quasar pairs at distances of $1-4$ Mpc. The black-hole masses for the pairs range between $M_{bh}\sim1-3\times 10^{9}~M_{\odot}/h$ and the quasar hosts are haloes of $M_h\sim1-3\times10^{14}~M_{\odot}/h$. Such pairs show signs of recent major mergers in the MBII simulation. By modeling the central and satellite AGN CLFs as log-normal and Schechter distributions respectively (as seen in MBII AGNs), we arrive at CLF models which fit the simulation predictions and observed luminosity function and the small-scale clustering measured for the SDSS sample. The small-scale clustering of our mock quasars is well-explained by central-satellite quasar pairs that reside in $M_h>10^{14}~M_{\odot}/h$ dark matter haloes. For these pairs, satellite quasar luminosity is similar to that of central quasars. Our CLF models imply a relatively steep increase in the maximum satellite luminosity, $L^*_{\mathrm{sat}}$, in haloes of $M_h>10^{14}~M_{\odot}/h$ with associated larger values of $L^*_{\mathrm{sat}}$ at higher redshift. This leads to increase in the satellite fraction that manifests itself in an enhanced clustering signal at $\lesssim$ 1 Mpc/h. For the ongoing eBOSS-CORE sample, we predict $\sim 200-500$ quasar pairs at $z\sim1.5$ (with $M_h \gtrsim10^{13}~M_{\odot}/h$ and $M_{bh} \gtrsim10^{8}~M_{\odot}/h$) at $\sim25~\mathrm{kpc}$ scales. Such a sample would be $\gtrsim10$ times larger than current pair samples.

We construct halo occupation distribution (HOD) models of high redshift ($z \gtrsim 7.5$) galaxies with $M_{*}>10^8~M_{\odot}/h$ using the BlueTides hydrodynamic simulation suite, with a particular emphasis on modelling the small scale / 1-halo clustering ($0.01\lesssim r \lesssim 1~ h^{-1}\rm{Mpc}$). Similar to low redshift studies, we find that the central and satellite mean HODs ($\left<N_{\mathrm{cen}}\right>$ and $\left<N_{\mathrm{sat}}\right>$) can be modeled by a smoothed step function and a power law respectively. The number density of satellite galaxies is however significantly suppressed compared to low redshift (satellite fractions drop from $\sim 50 \%$ at $z=0$ to $\lesssim 10 \%$ at $z=7.5$). The mean number of satellites, $\left<N_{\mathrm{sat}}\right> < 1$ for halo masses below $3 \times 10^{11} M_{\odot}/h$ (a rare halo at these redshifts). For the radial number density profiles, satellites with $10^8 \lesssim M^* \lesssim 10^{9} M_{\odot}/h$ in halos with $M_H \gtrsim 3 \times10^{11} M_{\odot}/h$ are consistent with NFW (with concentrations $c_{\mathrm{sat}} \sim 10-40$). Within halos of mass $M_H\lesssim 3 \times 10^{11} M_{\odot}/h$ satellites exhibit a power law profile with slope -3. Because these halos dominate the small scale clustering, the resulting 1-halo term is steeper than predicted using standard NFW profiles. Using this power-law profile for satellites, we can successfully reproduce the small-scale clustering exhibited by BlueTides galaxies using HOD modelling. We predict the highest probability of detecting satellites at $z>7.5$ is around centrals of $M^*\sim 3 \times 10^{10} M_{\odot}/h$ (with $M^{*}\gtrsim$ a few $10^{7} M_{\odot}/h$ ). This should be achievable with the James Webb Space Telescope (JWST).

SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$\mu$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$\sigma$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop.

In this paper, we study the filamentary structures and the galaxy alignment along filaments at redshift $z=0.06$ in the MassiveBlack-II simulation, a state-of-the-art, high-resolution hydrodynamical cosmological simulation which includes stellar and AGN feedback in a volume of (100 Mpc$/h$)$^3$. The filaments are constructed using the subspace constrained mean shift (SCMS; Ozertem & Erdogmus (2011) and Chen et al. (2015a)). First, we show that reconstructed filaments using galaxies and reconstructed filaments using dark matter particles are similar to each other; over $50\%$ of the points on the galaxy filaments have a corresponding point on the dark matter filaments within distance $0.13$ Mpc$/h$ (and vice versa) and this distance is even smaller at high-density regions. Second, we observe the alignment of the major principal axis of a galaxy with respect to the orientation of its nearest filament and detect a $2.5$ Mpc$/h$ critical radius for filament's influence on the alignment when the subhalo mass of this galaxy is between $10^9M_\odot/h$ and $10^{12}M_\odot/h$. Moreover, we find the alignment signal to increase significantly with the subhalo mass. Third, when a galaxy is close to filaments (less than $0.25$ Mpc$/h$), the galaxy alignment toward the nearest galaxy group depends on the galaxy subhalo mass. Finally, we find that galaxies close to filaments or groups tend to be rounder than those away from filaments or groups.

At the extreme densities and temperatures typical of the central engine of GRBs, the accreting torus is cooled mainly by advection and by neutrino emission. The latter process is dominated by electron and positron capture onto nucleons ($\beta$ reactions). We calculate the reaction rates and the nuclear composition of matter, assuming that the torus consists of helium, eletron-positron pairs, free neutrons and protons. After determining the equation of state and solving for the disk structure for a given initial accretion rate, we subsequently follow its time evolution. We find that, for accretion rates of the order of $10 M_{\odot}$/s, likely typical for the early stages of the accretion event, the disk becomes unstable, giving rise to variable energy output. This instability may play an important role for producing internal shocks.

We report evidence for an Fe K-alpha fluorescence line feature and disk reflection in the Very High, High, and Low State X-ray spectra of the galactic microquasar XTE J1748-288 during its June~1998 outburst. Spectral analyses are made on data gathered throughout the outburst by the Rossi X-ray Timing Explorer Proportional Counter Array. Gaussian line, relativistic disk emission line, and ionized disk reflection models are fit to the data. In the Very High State the line profile appears strongly redshifted, consistent with disk emission from the innermost stable orbits around a maximally rotating Kerr black hole. In the High State the line profile is less redshifted and increasingly prominent. The Low State line profile is very strong (~0.5 keV equivalent width) and centered at 6.7 +/- 0.10 keV; disk line emission model fits indicate that the inner edge of the disk fluctuates between ~20 and ~100 gravitational radii in this state. The disk reflection fraction is traced through the outburst; reflection from an ionized disk is preferred in the VHS and HS, and reflection from a relatively neutral disk is preferred in the LS. We discuss the implications of our findings for the binary system dynamics and accretion flow geometry in XTE J1748-288.