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.

We study the alignments of satellite galaxies, and their anisotropic distribution, with respect to location and orientation of their host central galaxy in MassiveBlack-II and IllustrisTNG simulations. We find that: the shape of the satellite system in halos of mass ($> 10^{13}h^{-1}M_{\odot}$) is well aligned with the shape of the central galaxy at $z=0.06$ with the mean alignment between the major axes being $\sim \Delta \theta = 12^{\circ}$ when compared to a uniform random distribution; that satellite galaxies tend to be anisotropically distributed along the major axis of the central galaxy with a stronger alignment in halos of higher mass or luminosity; and that the satellite distribution is more anisotropic for central galaxies with lower star formation rate, which are spheroidal, and for red central galaxies.Radially we find that satellites tend to be distributed along the major axis of the shape of the stellar component of central galaxies at smaller scales and the dark matter component on larger scales. We find that the dependence of satellite anisotropy on central galaxy properties and the radial distance is similar in both the simulations with a larger amplitude in MassiveBlack-II. The orientation of satellite galaxies tends to point toward the location of the central galaxy at small scales and this correlation decreases with increasing distance, and the amplitude of satellite alignment is higher in high mass halos. However, the projected ellipticities do not exhibit a scale-dependent radial alignment, as has been seen in some observational measurements.

The Wide Field Infrared Survey Telescope (WFIRST), with its field of view and high sensitivity will make surveys of cosmological large-scale structure possible at high redshifts. We investigate the possibility of detecting Baryon Acoustic Oscillations (BAO) at redshifts $z>7.6$ for use as a standard ruler. We use data from the hydrodynamic simulation \textsc{BlueTides} in conjunction with the gigaparsec-scale Outer Rim simulation and a model for patchy reionization to create mock WFIRST High Latitude Survey grism data for \lya\ emission line selected galaxies at redshifts $z=7.4$ to $z=10$, covering 2280 square degrees. We measure the monopoles of galaxies in the mock catalogues and fit the BAO features. We find that for a line flux of $L = 7\times 10^{-17} \ {\rm erg/s/cm}^{2}$, the $5 \sigma$ detection limit for the current design, the BAO feature is partially detectable (measured in three out of four survey quadrants analysed independently). The resulting root mean square error on the angular diameter distance to $z=7.7$ is 7.9$\%$. If we improve the detection sensitivity by a factor of two (i.e. $L = 3.5\times 10^{-17} \ {\rm erg/s/cm}^{2}$), the distance error reduces to $1.4\%$. We caution that many more factors are yet to be modelled, including dust obscuration, the damping wing due to the intergalactic medium, and low redshift interlopers. If these issues do not strongly affect the results, or different observational techniques (such as use of multiple lines) can mitigate them, WFIRST or similar instruments may be able to constrain the angular diameter distance to the high redshift Universe.

Halos and galaxies acquire their angular momentum during the collapse of surrounding large-scale structure. This process imprints alignments between galaxy spins and nearby filaments and sheets. Low mass halos grow by accretion onto filaments, aligning their spins with the filaments, whereas high mass halos grow by mergers along filaments, generating spins perpendicular to the filament. We search for this alignment signal using filaments identified with the "Cosmic Web Reconstruction" algorithm applied to the Sloan Digital Sky Survey Main Galaxy Sample and galaxy spins from the MaNGA integral-field unit survey. MaNGA produces a map of the galaxy's rotational velocity, allowing direct measurement of the galaxy's spin direction, or unit angular momentum vector projected onto the sky. We find no evidence for alignment between galaxy spins and filament directions. We do find hints of a mass-dependent alignment signal, which is in 2-3$\sigma$ tension with the mass-dependent alignment signal in the MassiveBlack-II and Illustris hydrodynamical simulations. However, the tension vanishes when galaxy spin is measured using the H$\alpha$ emission line velocity rather than stellar velocity. Finally, in simulations we find that the mass-dependent transition from aligned to anti-aligned dark matter halo spins is not necessarily present in stellar spins: we find a stellar spin transition in Illustris but not in MassiveBlack-II, highlighting the sensitivity of spin-filament alignments to feedback prescriptions and subgrid physics.

We investigate the redshift evolution of the intrinsic alignments (IA) of galaxies in the \texttt{MassiveBlackII} (MBII) simulation. We select galaxy samples above fixed subhalo mass cuts ($M_h>10^{11,12,13}~M_{\odot}/h$) at $z=0.6$ and trace their progenitors to $z=3$ along their merger trees. Dark matter components of $z=0.6$ galaxies are more spherical than their progenitors while stellar matter components tend to be less spherical than their progenitors. The distribution of the galaxy-subhalo misalignment angle peaks at $\sim10~\mathrm{deg}$ with a mild increase with time. The evolution of the ellipticity-direction~(ED) correlation amplitude $\omega(r)$ of galaxies (which quantifies the tendency of galaxies to preferentially point towards surrounding matter overdensities) is governed by the evolution in the alignment of underlying dark matter~(DM) subhaloes to the matter density of field, as well as the alignment between galaxies and their DM subhaloes. At scales $\sim1~\mathrm{cMpc}/h$, the alignment between DM subhaloes and matter overdensity gets suppressed with time, whereas the alignment between galaxies and DM subhaloes is enhanced. These competing tendencies lead to a complex redshift evolution of $\omega(r)$ for galaxies at $\sim1~\mathrm{cMpc}/h$. At scales $>1~\mathrm{cMpc}/h$, alignment between DM subhaloes and matter overdensity does not evolve significantly; the evolution of the galaxy-subhalo misalignment therefore leads to an increase in $\omega(r)$ for galaxies by a factor of $\sim4$ from $z=3$ to $0.6$ at scales $>1~\mathrm{cMpc}/h$. The balance between competing physical effects is scale dependant, leading to different conclusions at much smaller scales($\sim0.1~\mathrm{Mpc}/h$).

Many theoretical models predict that quasar driven outflows account for the observed quenching of star formation in massive galaxies. There is growing observational evidence for quasar-launched massive outflows at z>6, while the details of outflow-host galaxy interaction remain obscure. In this paper, we study the feedback around the highest redshift quasar in the BlueTides simulation, the largest volume cosmological hydrodynamic simulation so far carried out. We present predictions for gas outflows around the brightest $z = 7.54$ quasar which hosts the most massive black hole in the simulation volume, which has grown to black hole mass $6.7\times 10^{8}{\rm M}_\odot$ consistent with the current record holder for high-z quasars. We introduce a method to identify and trace the gas outflowing from the halo. We find that the total mass of the outflow gas is about $3.6\times 10^{9}{\rm M}_\odot$, constituting 6\% of the total gas in the halo. The outflow gas contains a cold, dense molecular component with mass about $2.6\times 10^{8}{\rm M}_\odot$, that originates from the inner region of the halo, within a few kpc of the central black hole. The velocities of the outflow gas reach thousands of km/s, within which the molecular component has mass averaged outward radial velocity of $1300$ km/s, consistent with observations. The averaged outflow rate is about $200-300 {\rm M}_\odot/yr$, with the outflowing gas mainly in a hotter ($T \sim 10^7$ K) and lower density state than the average of the host halo gas.

The most distant known quasar recently discovered by Ba\~nados et al. (2018) is at $z=7.5$ (690 Myr after the Big Bang), at the dawn of galaxy formation. We explore the host galaxy of the brightest quasar in the large volume cosmological hydrodynamic simulation BlueTides, which in Phase II has reached these redshifts. The brightest quasar in BlueTides has a luminosity of a $\sim$ few $10^{13} L_{\odot}$ and a black hole mass of $6.4 \times 10^{8} M_{\odot}$ at $z \sim 7.5$, comparable to the observed quasar (the only one in this large volume). The quasar resides in a rare halo of mass $M_{H} \sim 10^{12} M_{\odot}$ and has a host galaxy of stellar mass of $4 \times 10^{10}M_{\odot}$ with an ongoing (intrinsic) star formation rate of $\sim 80 M_{\odot} yr^{-1}$. The corresponding intrinsic UV magnitude of the galaxy is $-23.1$, which is roughly $2.7$ magnitudes fainter than the quasar's magnitude of $-25.9$. We find that the galaxy is highly metal enriched with a mean metallicity equal to the solar value. We derive quasar and galaxy spectral energy distribution (SED) in the mid and near infrared JWST bands. We predict a significant amount of dust attenuation in the rest-frame UV corresponding to $A_{1500} \sim 1.7$ giving an UV based SFR of $\sim 14 M_{\odot} yr^{-1}$. We present mock JWST images of the galaxy with and without central point source, in different MIRI and NIRCam filters. The host galaxy is detectable in NIRCam filters, but it is extremely compact ($R_{E}=0.35$ kpc). It will require JWST's exquisite sensitivity and resolution to separate the galaxy from the central point source. Finally within the FOV of the quasar in BlueTides there are two more sources that would be detectable by 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.

The use of high-quality simulated sky catalogs is essential for the success of cosmological surveys. The catalogs have diverse applications, such as investigating signatures of fundamental physics in cosmological observables, understanding the effect of systematic uncertainties on measured signals and testing mitigation strategies for reducing these uncertainties, aiding analysis pipeline development and testing, and survey strategy optimization. The list of applications is growing with improvements in the quality of the catalogs and the details that they can provide. Given the importance of simulated catalogs, it is critical to provide rigorous validation protocols that enable both catalog providers and users to assess the quality of the catalogs in a straightforward and comprehensive way. For this purpose, we have developed the DESCQA framework for the Large Synoptic Survey Telescope Dark Energy Science Collaboration as well as for the broader community. The goal of DESCQA is to enable the inspection, validation, and comparison of an inhomogeneous set of synthetic catalogs via the provision of a common interface within an automated framework. In this paper, we present the design concept and first implementation of DESCQA. In order to establish and demonstrate its full functionality we use a set of interim catalogs and validation tests. We highlight several important aspects, both technical and scientific, that require thoughtful consideration when designing a validation framework, including validation metrics and how these metrics impose requirements on the synthetic sky catalogs.

Supermassive blackholes with masses of a billion solar masses or more are known to exist up to $z=7$. However, the present-day environments of the descendants of first quasars is not well understood and it is not known if they live in massive galaxy clusters or more isolated galaxies at $z=0$. We use a dark matter-only realization (BTMassTracer) of the BlueTides cosmological hydrodynamic simulation to study the halo properties of the descendants of the most massive black holes at $z=8$. We find that the descendants of the quasars with most massive black holes are not amongst the most massive halos. They reside in halos of with group-like ($\sim 10^{14}M_{\odot}$) masses, while the most massive halos in the simulations are rich clusters with masses $\sim 10^{15} M_{\odot}$. The distribution of halo masses at low redshift is similar to that of the descendants of least massive black holes, for a similar range of halo masses at $z=8$, which indicates that they are likely to exist in similar environments. By tracing back to the $z = 8$ progenitors of the most massive (cluster sized) halos at $z=0$; we find that their most likely black hole mass is less than $10^7 M_{\odot}$; they are clearly not amongst the most massive black holes. We also provide estimates for the likelihood of finding a high redshift quasar hosting a black hole with masses above $10^{7} M_{\odot}$ for a given halo mass at $z=0$. For halos above $10^{15} M_{\odot}$, there is only $20 \%$ probability that their $z=8$ progenitors hosted a black hole with mass above $10^{7} M_{\odot}$.

A strong correlation has been measured between the observed centripetal accelerations in galaxies and the accelerations implied by the baryonic components of galaxies. This empirical radial acceleration relation must be accounted for in any viable model of galaxy formation. We measure and compare the radial accelerations contributed by baryons and by dark matter in disk galaxies in the MassiveBlack-II hydrodynamic galaxy formation simulation. The sample of 1594 galaxies spans three orders of magnitude in luminosity and four in surface brightness, comparable to the observed sample from the Spitzer Photometry & Accurate Rotation Curves (SPARC) dataset used by McGaugh et al. (2016). We find that radial accelerations contributed by baryonic matter only and by total matter are highly correlated, with only small scatter around their mean or median relation, despite the wide ranges of galaxy luminosity and surface brightness. We further find that the radial acceleration relation in this simulation differs from that of the SPARC sample, and can be described by a simple power law in the acceleration range we are probing.

We study the shapes and intrinsic alignments of disks and elliptical galaxies in the MassiveBlack-II (MBII) and Illustris cosmological hydrodynamic simulations, with volumes of $(100h^{-1}Mpc)^{3}$ and $(75h^{-1}Mpc)^{3}$ respectively. We find that simulated disk galaxies are more oblate in shape and more misaligned with the shape of their host dark matter subhalo when compared with ellipticals. The disk major axis is found to be oriented towards the location of nearby elliptical galaxies. We also find that the disks are thinner in MBII and misalignments with dark matter halo orientations are smaller in both disks and ellipticals when compared with Illustris. As a result, the intrinsic alignment correlation functions at fixed mass have a higher amplitude in MBII than in Illustris. Finally, at scales above $\sim 0.1$\hmpc, the intrinsic alignment two-point correlation functions for disk galaxies in both simulations are consistent with a null detection, unlike those for ellipticals. Despite significant differences in the treatments of hydrodynamics and baryonic physics in the simulations, we find that the $w_{\delta +}$ correlation function scales similarly with transverse separation. However, the less massive galaxies show different scale dependence in the ED correlation. This result indicates that, while hydrodynamic simulations are a promising tool to study intrinsic alignments, further study is needed to understand the impact of differences in the implementations of hydrodynamics and baryonic feedback.

We explore the effects of specific assumptions in the subgrid models of star formation and stellar and AGN feedback on intrinsic alignments of galaxies in cosmological simulations of "MassiveBlack-II" family. Using smaller volume simulations, we explored the parameter space of the subgrid star formation and feedback model and found remarkable robustness of the observable statistical measures to the details of subgrid physics. The one observational probe most sensitive to modeling details is the distribution of misalignment angles. We hypothesize that the amount of angular momentum carried away by the galactic wind is the primary physical quantity that controls the orientation of the stellar distribution. Our results are also consistent with a similar study by the EAGLE simulation team.

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.

Whether among the myriad tiny proto-galaxies there exists a population with similarities to present day galaxies is an open question. We show, using BlueTides, the first hydrodynamic simulation large enough to resolve the relevant scales, that the first massive galaxies to form are %in fact predicted to have extensive rotationally-supported disks. Although their morphology resembles in some ways Milky-way types seen at much lower redshifts, these high-redshift galaxies are smaller, denser, and richer in gas than their low redshift counterparts. From a kinematic analysis of a statistical sample of 216 galaxies at redshift $z=8-10$ we have found that disk galaxies make up 70\% of the population of galaxies with stellar mass $10^{10} M_\odot$ or greater. Cold Dark Matter cosmology therefore makes specific predictions for the population of large galaxies 500 million years after the Big Bang. We argue that wide-field satellite telescopes (e.g. WFIRST) will in the near future discover these first massive disk galaxies. The simplicity of their structure and formation history should make possible new tests of cosmology.

We compare the shapes and intrinsic alignments of galaxies in the MassiveBlack-II cosmological hydrodynamic simulation (MBII) to those in a dark matter-only (DMO) simulation performed with the same volume (100$h^{-1}$Mpc)$^{3}$, cosmological parameters, and initial conditions. Understanding the impact of baryonic physics on galaxy shapes and alignments and their relation to the dark matter distribution should prove useful to map the intrinsic alignments of galaxies from hydrodynamic to dark matter-only simulations. We find that dark matter subhalos are typically rounder in MBII, and the shapes of stellar matter in low mass galaxies are more misaligned with the shapes of the dark matter of the corresponding subhalos in the DMO simulation. At $z=0.06$, the fractional difference in the mean misalignment angle between MBII and DMO simulations varies from $\sim 28 \% - 12 \%$ in the mass range $10^{10.8} - 6.0 \times 10^{14} h^{-1}M_{\odot}$. We study the dark matter halo shapes and alignments as a function of radius, and find that while galaxies in MBII are more aligned with the inner parts of their dark matter subhalos, there is no radial trend in their alignments with the corresponding subhalo in the DMO simulation. This result highlights the importance of baryonic physics in determining the alignment of the galaxy with respect to the inner parts of the halo. Finally, we compare the ellipticity-direction (ED) correlation for galaxies to that for dark matter halos, finding that it is suppressed on all scales by stellar-dark matter misalignment. In the projected shape-density correlation ($w_{\delta+}$), which includes ellipticity weighting, this effect is partially canceled by the higher mean ellipticities of the stellar component, but differences of order $30-40\%$ remain on scales $> 1$ Mpc over a range of subhalo masses, with scale-dependent effects below $1$ Mpc.

The intrinsic alignment of galaxies with the large-scale density field is an important astrophysical contaminant in upcoming weak lensing surveys. We present detailed measurements of the galaxy intrinsic alignments and associated ellipticity-direction (ED) and projected shape ($w_{g+}$) correlation functions for galaxies in the cosmological hydrodynamic MassiveBlack-II (MB-II) simulation. We carefully assess the effects on galaxy shapes, misalignment of the stellar component with the dark matter shape and two-point statistics of iterative weighted (by mass and luminosity) definitions of the (reduced and unreduced) inertia tensor. We find that iterative procedures must be adopted for a reliable measurement of the reduced tensor but that luminosity versus mass weighting has only negligible effects. Both ED and $w_{g+}$ correlations increase in amplitude with subhalo mass (in the range of $10^{10} - 6.0\times 10^{14}h^{-1}M_{\odot}$), with a weak redshift dependence (from $z=1$ to $z=0.06$) at fixed mass. At $z \sim 0.3$, we predict a $w_{g+}$ that is in reasonable agreement with SDSS LRG measurements and that decreases in amplitude by a factor of $\sim 5$--18 for galaxies in the LSST survey. We also compared the intrinsic alignments of centrals and satellites, with clear detection of satellite radial alignments within their host halos. Finally, we show that $w_{g+}$ (using subhalos as tracers of density) and $w_{\delta+}$ (using dark matter density) predictions from the simulations agree with that of non-linear alignment models (NLA) at scales where the 2-halo term dominates in the correlations (and tabulate associated NLA fitting parameters). The 1-halo term induces a scale dependent bias at small scales which is not modeled in the NLA model.

The intrinsic alignment of galaxy shapes with the large-scale density field is a contaminant to weak lensing measurements, as well as being an interesting signature of galaxy formation and evolution (albeit one that is difficult to predict theoretically). Here we investigate the shapes and relative orientations of the stars and dark matter of halos and subhalos (central and satellite) extracted from 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{h^{-1}\mathrm{Mpc}})^3$. We consider redshift evolution from $z=1$ to $0.06$ and mass evolution within the range of subhalo masses, $10^{10} -6.0 \times 10^{14.0}{h^{-1}M_{\odot}}$. The shapes of the dark matter distributions are generally more round than the shapes defined by stellar matter. The projected root-mean-square (RMS) ellipticity per component for stellar matter is measured to be $e_\text{rms} = 0.28$ at $z=0.3$ for $M_{subhalo}> 10^{12.0}{h^{-1}M_{\odot}}$, which compares favourably with observational measurements. We find that the shapes of stellar and dark matter are more round for less massive subhalos and at lower redshifts. By directly measuring the relative orientation of the stellar matter and dark matter of subgroups, we find that, on average, the misalignment between the two components is larger for less massive subhalos. The mean misalignment angle varies from $\sim 30^{\circ}-10^{\circ}$ for $M \sim 10^{10} - 10^{14} {h^{-1}M_{\odot}}$ and shows a weak dependence on redshift. We also compare the misalignment angles in central and satellite subhalos at fixed subhalo mass, and find that centrals are more misaligned than satellites. We present fitting formulae for the shapes of dark and stellar matter in subhalos and also the probability distributions of misalignment angles.