We evaluate the impact of imaging systematics on the clustering of luminous red galaxies (LRG), emission-line galaxies (ELG) and quasars (QSO) targeted for the upcoming Dark Energy Spectroscopic Instrument (DESI) survey. Using Data Release 7 of the DECam Legacy Survey, we study the effects of astrophysical foregrounds, stellar contamination, differences between north galactic cap and south galactic cap measurements, and variations in imaging depth, stellar density, galactic extinction, seeing, airmass, sky brightness, and exposure time before presenting survey masks and weights to mitigate these effects. With our sanitized samples in hand, we conduct a preliminary analysis of the clustering amplitude and evolution of the DESI main targets. From measurements of the angular correlation functions, we determine power law fits $r_0 = 7.78 \pm 0.26$ $h^{-1}$Mpc, $\gamma = 1.98 \pm 0.02$ for LRGs and $r_0 = 5.45 \pm 0.1$ $h^{-1}$Mpc, $\gamma = 1.54 \pm 0.01$ for ELGs. Additionally, from the angular power spectra, we measure the linear biases and model the scale dependent biases in the weakly nonlinear regime. Both sets of clustering measurements show good agreement with survey requirements for LRGs and ELGs, attesting that these samples will enable DESI to achieve precise cosmological constraints. We also present clustering as a function of magnitude, use cross-correlations with external spectroscopy to infer $dN/dz$ and measure clustering as a function of luminosity, and probe higher order clustering statistics through counts-in-cells moments.

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.

Spectroscopic phase curves provide unique access to the three-dimensional properties of transiting exoplanet atmospheres. However, a modeling framework must be developed to deliver accurate inferences of atmospheric properties for these complex data sets. Here, we develop an approach to retrieve temperature structures and molecular abundances from phase curve spectra at any orbital phase. In the context of a representative hot Jupiter with a large day-night temperature contrast, we examine the biases in typical one-dimensional (1D) retrievals as a function of orbital phase/geometry, compared to two-dimensional (2D) models that appropriately capture the disk-integrated phase geometry. We guide our intuition by applying our new framework on a simulated HST+Spitzer phase curve data set in which the ``truth'' is known, followed by an application to the spectroscopic phase curve of the canonical hot Jupiter, WASP-43b. We also demonstrate the retrieval framework on simulated JWST phase curve observations. We apply our new geometric framework to a joint-fit of all spectroscopic phases, assuming longitudinal molecular abundance homogeneity, resulting in an a factor of 2 improvement in abundances precision when compared to individual phase constraints. With a 1D retrieval model on simulated HST+Spitzer data, we find strongly biased molecular abundances for CH$_4$ and CO$_2$ at most orbital phases. With 2D, the day and night profiles retrieved from WASP-43b remain consistent throughout the orbit. JWST retrievals show that a 2D model is strongly favored at all orbital phases. Based on our new 2D retrieval implementation, we provide recommendations on when 1D models are appropriate and when more complex phase geometries involving multiple TP profiles are required to obtain an unbiased view of tidally locked planetary atmospheres.

Supermassive black holes are commonly found in the center of galaxies and evolve with their hosts. The supermassive binary black holes (SMBBH) are thus expected to exist in close galaxy pairs, however, none has been unequivocally detected. The square kilometre array (SKA) is a multi-purpose radio telescope with a collecting area approaching 1 million square metres, with great potential for detecting nanohertz gravitational waves (GWs). In this paper, we quantify the GW detectability by SKA for a realistic SMBBH population using pulsar timing array (PTA) technique and quantify its impact on revealing SMBBH evolution with redshift for the first time. With only $\sim20$ pulsars, much smaller a requirement than in previous work, the SKA PTA is expected to obtain detection within about 5 years of operation and to achieve a detection rate of more than 100 SMBBHs/yr after about 10 years. Although beyond the scope of this paper, we must acknowledge that the presence of persistent red noise will reduce the number of expected detections here. It is thus imperative to understand and mitigate red noise in the PTA data. The GW signatures from a few well-known SMBBH candidates, such as OJ 287, 3C 66B, NGC 5548 and Ark 120, will be detected given the currently best-known parameters of each system. Within 30 years of operation, about 60 individual SMBBHs detection with $z<0.05$ and more than $10^4$ with $z<1$ are expected. The detection rate drops precipitately beyond $z=1$. The substantial number of expected detections and their discernible evolution with redshift by SKA PTA will make SKA a significant tool for studying SMBBHs.

We present timing solutions and analyses of 11 pulsars discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). These pulsars were discovered using an ultra-wide bandwidth receiver in drift-scan observations made during the commissioning phase of FAST, and were then confirmed and timed using the 64-m Parkes Radio Telescope. Each pulsar has been observed over a span of at least one year. Highlighted discoveries include PSR J0344-0901, which displays mode-changing behaviour and may belong to the class of so-called `swooshing' pulsars (alongside PSRs B0919+06 and B1859+07); PSR J0803-0942, whose emission is almost completely linearly polarised; and PSRs J1900-0134 and J1945+1211, whose well defined polarisation angle curves place stringent constraints on their emission geometry. We further discuss the detectability of these pulsars by earlier surveys, and highlight lessons learned from our work in carrying out confirmation and monitoring observations of pulsars discovered by a highly sensitive telescope, many of which may be applicable to next-generation pulsar surveys. This paper marks one of the first major releases of FAST-discovered pulsars, and paves the way for future discoveries anticipated from the Commensal Radio Astronomy FAST Survey (CRAFTS).

We revisit techniques for performing cosmological simulations with both baryons and cold dark matter when each fluid has different initial conditions, as is the case at the end of the radiation era. Most simulations do not reproduce the linear prediction for the difference between the cold dark matter and baryon perturbations. We show that this is due to the common use of offset regular grids when setting up the particle initial conditions. The correct behaviour can be obtained without any loss of simulation resolution by using a Lagrangian glass for the baryon particles. We further show that the difference between cold dark matter and baryons may affect predictions for the Lyman-alpha forest flux power spectrum at the 5% level, potentially impacting current cosmological constraints.

The paper examines the early growth of supermassive black holes (SMBHs) in cosmological hydrodynamic simulations with different BH seeding scenarios. Employing the constrained Gaussian realization, we reconstruct the initial conditions in the large-volume BlueTides simulation and run them to $z=6$ to cross-validate that the method reproduces the first quasars and their environments. Our constrained simulations in a volume of $(15\, h^{-1}{\rm Mpc})^3$ successfully recover the evolution of large-scale structure and the stellar and BH masses in the vicinity of a $\sim10^{12}\, M_{\odot}$ halo which we identified in BlueTides at $z\sim7$ hosting a $\sim10^9\, M_{\odot}$ SMBH. Among our constrained simulations, only the ones with a low-tidal field and high-density peak in the initial conditions induce the fastest BH growth required to explain the $z>6$ quasars. We run two sets of simulations with different BH seed masses of $5\times10^3$, $5\times10^4$, and $5\times10^5\, h^{-1}M_{\odot}$, (a) with the same ratio of halo to BH seed mass and (b) with the same halo threshold mass. At $z=6$, all the SMBHs converge in mass to $\sim10^9\, M_{\odot}$ except for the one with the smallest seed in (b) undergoing critical BH growth and reaching $10^8$ -- $10^9\, M_{\odot}$, albeit with most of the growth in (b) delayed compared to set (a). The finding of eight BH mergers in the small-seed scenario (four with masses $10^4$ -- $10^6\, M_{\odot}$ at $z>12$), six in the intermediate-seed scenario, and zero in the large-seed scenario suggests that the vast BHs in the small-seed scenario merge frequently during the early phases of the growth of SMBHs. The increased BH merger rate for the low-mass BH seed and halo threshold scenario provides an exciting prospect for discriminating BH formation mechanisms with the advent of multi-messenger astrophysics and next-generation gravitational wave facilities.

We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the corresponding fluence $>0.2$~Jy~ms. Based on the observed dispersion measure of 1812~cm$^{-3}$~pc, we infer a redshift of $\sim 1.9$. From this, we estimate the peak luminosity and isotropic energy to be $\lesssim 2\times10^{43}$~erg~s$^{-1}$ and $\lesssim 2\times10^{40}$~erg, respectively. With only one FRB from the survey detected so far, our constraints on the event rate are limited. We derive a 95\% confidence lower limit for the event rate of 900 FRBs per day for FRBs with fluences $>0.025$~Jy~ms. We performed follow-up observations of the source with FAST for four hours and have not found a repeated burst. We discuss the implications of this discovery for our understanding of the physical mechanisms of FRBs.

We present \textit{Spitzer} full-orbit thermal phase curves of the hot Jupiter Qatar-1b, a planet with the same equilibrium temperature---and intermediate surface gravity and orbital period---as the well-studied planets HD 209458b and WASP-43b. We measure secondary eclipse of $0.21 \pm 0.02 \%$ at $3.6~\mu$m and $0.30 \pm 0.02 \%$ at $4.5~\mu$m, corresponding to dayside brightness temperatures of $1542^{+32}_{-31}$~K and $1557^{+35}_{-36}$~K, respectively, consistent with a vertically isothermal dayside. The respective nightside brightness temperatures are $1117^{+76}_{-71}$~K and $1167^{+69}_{-74}$~K, in line with a trend that hot Jupiters all have similar nightside temperatures. We infer a Bond albedo of $0.12_{-0.16}^{+0.14}$ and a moderate day-night heat recirculation efficiency, similar to HD 209458b. General circulation models for HD 209458b and WASP-43b predict that their bright-spots should be shifted east of the substellar point by tens of degrees, and these predictions were previously confirmed with \textit{Spitzer} full-orbit phase curve observations. The phase curves of Qatar-1b are likewise expected to exhibit eastward offsets. Instead, the observed phase curves are consistent with no offset: $11^{\circ}\pm 7^{\circ}$ at $3.6~\mu$m and $-4^{\circ}\pm 7^{\circ}$ at $4.5~\mu$m. The discrepancy in circulation patterns between these three otherwise similar planets points to the importance of secondary parameters like rotation rate and surface gravity, and the presence or absence of clouds, in determining atmospheric conditions on hot Jupiters.

Pulsar timing arrays (PTAs) have been proven as one of the most potential approaches to detect low-frequency gravitational waves in the near future. Although PTAs have not captured any GW signals at present, there have been a large number of related theoretical researches and some meaningful detection limits. In this paper, we focus on the nanohertz gravitational waves (GWs) from individual supermassive binary black holes (SMBBHs). Given a specific pulsar PSR J1909$-$3744, the corresponding GW$-$induced timing residuals in PTAs with Gaussian white noise can be simulated. Then we present the classification of simulated PTA data and parameter estimation for potential GW sources using machine learning based on neural networks. As a classifier, the convolutional neural network (CNN) shows a great performance when the signal to noise ratio $\geq$7 (the corresponding dimensionless amplitude of GW is $h\geq 2\times 10^{-15}$). On the other hand, we apply recurrent neural network (RNN) and Bayesian neural networks (BNN) to chirp mass ($\mathcal{M}$) estimation, which are able to provide the mean relative error and uncertainties of $\mathcal{M}$. That is crucial to astrophysical observation. In our case, the mean relative error of chirp mass estimation is less than $15\%$ with PTA sensitivity. Although these results are achieved for simulated PTA data, it's important for realizing intelligent processing on PTA data analysis.

Nebular emission lines associated with galactic HII regions carry information about both physical properties of the ionised gas and the source of ionising photons as well as providing the opportunity of measuring accurate redshifts and thus distances once a cosmological model is assumed. While nebular line emission has been extensively studied at lower redshift there are currently only few constraints within the epoch of reionisation (EoR, $z>6$), chiefly due to the lack of sensitive near-IR spectrographs. However, this will soon change with the arrival of the Webb Telescope providing sensitive near-IR spectroscopy covering the rest-frame UV and optical emission of galaxies in the EoR. In anticipation of Webb we combine the large cosmological hydrodynamical simulation Bluetides with photoionisation modelling to predict the nebular emission line properties of galaxies at $z=8\to 13$. We find good agreement with the, albeit limited, existing direct and indirect observational constraints on equivalent widths though poorer agreement with luminosity function constraints.

The Habitable Exoplanet Observatory, or HabEx, has been designed to be the Great Observatory of the 2030s. For the first time in human history, technologies have matured sufficiently to enable an affordable space-based telescope mission capable of discovering and characterizing Earthlike planets orbiting nearby bright sunlike stars in order to search for signs of habitability and biosignatures. Such a mission can also be equipped with instrumentation that will enable broad and exciting general astrophysics and planetary science not possible from current or planned facilities. HabEx is a space telescope with unique imaging and multi-object spectroscopic capabilities at wavelengths ranging from ultraviolet (UV) to near-IR. These capabilities allow for a broad suite of compelling science that cuts across the entire NASA astrophysics portfolio. HabEx has three primary science goals: (1) Seek out nearby worlds and explore their habitability; (2) Map out nearby planetary systems and understand the diversity of the worlds they contain; (3) Enable new explorations of astrophysical systems from our own solar system to external galaxies by extending our reach in the UV through near-IR. This Great Observatory science will be selected through a competed GO program, and will account for about 50% of the HabEx primary mission. The preferred HabEx architecture is a 4m, monolithic, off-axis telescope that is diffraction-limited at 0.4 microns and is in an L2 orbit. HabEx employs two starlight suppression systems: a coronagraph and a starshade, each with their own dedicated instrument.

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.

We perform a search for gravitational waves (GWs) from several supermassive binary black hole (SMBBH) candidates (NGC 5548, Mrk 231, OJ 287, PG 1302-102, NGC 4151, Ark 120 and 3C 66B) in long-term timing observations of the pulsar PSR J1909$-$3744 obtained using the Parkes radio telescope. No statistically significant signals were found. We constrain the chirp masses of those SMBBH candidates and find the chirp mass of NGC 5548 and 3C 66B to be less than $2.4 \times 10^9\,\rm M_{\odot}$ and $2.5 \times 10^9\,\rm M_{\odot}$ (with 95% confidence), respectively. Our upper limits remain a factor of 3 to 370 above the likely chirp masses for these candidates as estimated from other approaches. The observations processed here provide upper limits on the GW strain amplitude that improve upon the results from the first Parkes Pulsar Timing Array data release by a factor of 2 to 7. We investigate how information about the orbital parameters can help improve the search sensitivity for individual SMBBH systems. Finally, we show that these limits are insensitive to uncertainties in the Solar System ephemeris model.

High-$z$ AGNs hosted in gas rich galaxies are expected to grow through significantly obscured accretion phases. This may limit or bias their observability. In this work, we use \textsc{BlueTides}, a large volume cosmological simulation of galaxy formation to examine quasar obscuration for the highest-redshift ($z \geq 7$) supermassive black holes residing in the center of galaxies. We find that for the bright quasars, most of the high column density gas ($>90\%$) resides in the innermost regions of the host galaxy, (typically within $< 10$ ckpc), while the gas in the outskirts is a minor contributor to the $N_\mathrm H$. The brightest quasars can have large angular variations in galactic obscuration, over 2 orders of magnitude, where the lines of sight with the lowest obscuration are those formed via strong gas outflows driven by AGN feedback. We find that for the overall AGN population, the mean $N_\mathrm H$ is generally larger for high luminosity and BH mass, while the $N_\mathrm H$ distribution is significantly broadened, developing a low $N_\mathrm H $ wing due to the angular variations driven by the AGN outflows/feedback. The obscured fraction P($N_{\rm H} > 10^{23} {\rm cm}^{-2}$) typically range from 0.6 to 1.0 for increasing $L_{X}$ (with $L_X > 10^{43} \rm{ergs/s}$), with no clear trend of redshift evolution. With respect to the galaxy host property, we find a linear relation between $N_{\rm H}$, $M_*$ and $M_{\rm H_2}$ with $\log N_{\rm H} = (0.24 \pm 0.03) \log M_{*} + (20.7 \pm 0.3)$ and $\log N_{\rm H} = (0.47 \pm 0.03) \log M_{\rm H_2} + (18.4 \pm 0.3)$. The dust optical depth in the UV band $\tau_{\mathrm UV}$ has tight positive correlation with $N_{\rm H}$. Our dust extincted UVLF is about 1.5 dex lower than the intrinsic UVLF, implying that more than 99\% of the $z \sim 7$ AGNs are heavily dust extincted and therefore would be missed by the UV band observation.

The characterization of rocky, Earth-like planets is an important goal for future large ground- and space-based telescopes. In support of developing an efficient observational strategy, we have applied Bayesian statistical inference to interpret the albedo spectrum of cloudy true-Earth analogs that include a diverse spread in their atmospheric water vapor mixing ratios. We focus on detecting water-bearing worlds by characterizing their atmospheric water vapor content via the strong 0.94$\,\mu$m H$_2$O absorption feature, with several observational configurations. Water vapor is an essential signpost when assessing planetary habitability, and determining its presence is important in vetting whether planets are suitable for hosting life. We find that R=140 spectroscopy of the absorption feature combined with a same-phase green optical photometric point at $0.525-0.575\,\mu$m is capable of distinguishing worlds with less than $0.1\times$ Earth-like water vapor levels from worlds with $1\times$ Earth-like levels or greater at a signal-to-noise ratio of 5 or better with $2\sigma$ confidence. This configuration can differentiate between $0.01\times$ and $0.1\times$ Earth-like levels when the signal-to-noise ratio is 10 or better at the same confidence. However, strong constraints on the water vapor mixing ratio remained elusive with this configuration even at signal-to-noise of 15. We find that adding the same-phase optical photometric point does not significantly help characterize the H$_2$O mixing ratio, but does enable an upper limit on atmospheric ozone levels. Finally, we find that a 0.94$\,\mu$m photometric point, instead of spectroscopy, combined with the green-optical point, fails to produce meaningful information about atmospheric water content.

This paper introduces the Hidden Valley simulations, a set of trillion-particle N-body simulations in gigaparsec volumes aimed at intensity mapping science. We present details of the simulations and their convergence, then specialize to the study of 21-cm fluctuations between redshifts 2 and 6. Neutral hydrogen is assigned to halos using three prescriptions, and we investigate the clustering in real and redshift-space at the 2-point level. In common with earlier work we find the bias of HI increases from near 2 at z = 2 to 4 at z = 6, becoming more scale dependent at high z. The level of scale-dependence and decorrelation with the matter field are as predicted by perturbation theory. Due to the low mass of the hosting halos, the impact of fingers of god is small on the range relevant for proposed 21-cm instruments. We show that baryon acoustic oscillations and redshift-space distortions could be well measured by such instruments. Taking advantage of the large simulation volume, we assess the impact of fluctuations in the ultraviolet background, which change HI clustering primarily at large scales.

The Quijote simulations are a set of 43100 full N-body simulations spanning more than 7000 cosmological models in the $\{\Omega_{\rm m}, \Omega_{\rm b}, h, n_s, \sigma_8, M_\nu, w \}$ hyperplane. At a single redshift the simulations contain more than 8.5 trillions of particles over a combined volume of 43100 $(h^{-1}{\rm Gpc})^3$. Billions of dark matter halos and cosmic voids have been identified in the simulations, whose runs required more than 35 million core hours. The Quijote simulations have been designed for two main purposes: 1) to quantify the information content on cosmological observables, and 2) to provide enough data to train machine learning algorithms. In this paper we describe the simulations and show a few of their applications. We also release the Petabyte of data generated, comprising hundreds of thousands of simulation snapshots at multiple redshifts, halo and void catalogs, together with millions of summary statistics such as power spectra, bispectra, correlation functions, marked power spectra, and estimated probability density functions.

Generating mocks for future sky surveys requires large volumes and high resolutions, which is computationally expensive even for fast simulations. In this work we try to develop numerical schemes to calibrate various halo and matter statistics in fast low resolution simulations compared to high resolution N-body and hydrodynamic simulations. For the halos, we improve the initial condition accuracy and develop a halo finder "relaxed-FOF", where we allow different linking length for different halo mass and velocity dispersions. We show that our relaxed-FoF halo finder improves the common statistics, such as halo bias, halo mass function, halo auto power spectrum in real space and in redshift space, cross correlation coefficient with the reference halo catalog, and halo-matter cross power spectrum. We also incorporate the potential gradient descent (PGD) method into fast simulations to improve the matter distribution at nonlinear scale. By building a lightcone output, we show that the PGD method significantly improves the weak lensing convergence tomographic power spectrum. With these improvements FastPM is comparable to the high resolution full N-body simulation of the same mass resolution, with two orders of magnitude fewer time steps. These techniques can be used to improve the halo and matter statistics of FastPM simulations for mock catalogs of future surveys such as DESI and LSST.

Matter evolved under influence of gravity from minuscule density fluctuations. Non-perturbative structure formed hierarchically over all scales, and developed non-Gaussian features in the Universe, known as the Cosmic Web. To fully understand the structure formation of the Universe is one of the holy grails of modern astrophysics. Astrophysicists survey large volumes of the Universe and employ a large ensemble of computer simulations to compare with the observed data in order to extract the full information of our own Universe. However, to evolve trillions of galaxies over billions of years even with the simplest physics is a daunting task. We build a deep neural network, the Deep Density Displacement Model (hereafter D$^3$M), to predict the non-linear structure formation of the Universe from simple linear perturbation theory. Our extensive analysis, demonstrates that D$^3$M outperforms the second order perturbation theory (hereafter 2LPT), the commonly used fast approximate simulation method, in point-wise comparison, 2-point correlation, and 3-point correlation. We also show that D$^3$M is able to accurately extrapolate far beyond its training data, and predict structure formation for significantly different cosmological parameters. Our study proves, for the first time, that deep learning is a practical and accurate alternative to approximate simulations of the gravitational structure formation of the Universe.

Lijiang 2.4-meter Telescope(LJT), the largest common-purpose optical telescope in China, has been applied to the world-wide astronomers since 2008. It is located at Gaomeigu site, Lijiang Observatory(LJO), the southwest of China. The site has very good observational conditions. Since 10-year operation, several instruments have been equipped on the LJT. Astronomers can perform both photometric and spectral observations. The main scientific goals of LJT include photometric and spectral evolution of supernova, reverberation mapping of active galactic nucleus, physical properties of binary star and near-earth object(comet and asteroid), identification of exoplanet, and all kinds of transients. Until now, the masses of 41 high accretion rate black holes have been measured, and more than 168 supernova have been identified by the LJT. More than 190 papers related to the LJT have been published. In this paper, the general observation condition of the Gaomeigu site is introduced at first. Then, the LJT structure is described in detail, including the optical, mechanical, motion and control system. The specification of all the instruments, and some detailed parameters of the YFOSC is also presented. Finally, some important scientific results and future expectations are summarized.

We present the UNIT $N$-body cosmological simulations project, designed to provide precise predictions for nonlinear statistics of the galaxy distribution. We focus on characterizing statistics relevant to emission line and luminous red galaxies in the current and upcoming generation of galaxy surveys. We use a suite of precise particle mesh simulations (FastPM) as well as with full $N$-body calculations with a mass resolution of $\sim 1.2\times10^9\,h^{-1}$M$_{\odot}$ to investigate the recently suggested technique of Angulo & Pontzen 2016 to suppress the variance of cosmological simulations We study redshift space distortions, cosmic voids, higher order statistics from $z=2$ down to $z=0$. We find that both two- and three-point statistics are unbiased. Over the scales of interest for baryon acoustic oscillations and redshift-space distortions, we find that the variance is greatly reduced in the two-point statistics and in the cross correlation between halos and cosmic voids, but is not reduced significantly for the three-point statistics. We demonstrate that the accuracy of the two-point correlation function for a galaxy survey with effective volume of 20 ($h^{-1}$Gpc)$^3$ is improved by about a factor of 40, indicating that two pairs of simulations with a volume of 1 ($h^{-1}$Gpc)$^3$ lead to the equivalent variance of $\sim$150 such simulations. The $N$-body simulations presented here thus provide an effective survey volume of about seven times the effective survey volume of DESI or Euclid. The data from this project, including dark matter fields, halo catalogues, and their clustering statistics, are publicly available at http://www.unitsims.org.

Based on the Chinese historical sunspots drawings, a data set consisting of the scanned images and all their digitized parameters from 1925 to 2015 have been constructed. In this paper, we briefly describe the developmental history of sunspots drawings in China. This paper describes the preliminary processing processes that strat from the initial data (inputing to the scanning equipment) to the parameters extraction, and finally summarizes the general features of this dataset. It is the first systematic project in Chinese solar-physics community that the historical observation of sunspots drawings were digitized. Our data set fills in an almost ninety years historical gap, which span 60 degrees from east to west and 50 degrees from north to south and have no continuous and detailed digital sunspot observation information. As a complementary to other sunspots observation in the world, our dataset provided abundant information to the long term solar cycles solar activity research.

During the last decades, rapid progress has been made in measurements of the rest-frame ultraviolet (UV) luminosity function (LF) for high-redshift galaxies ($z \geq 6$). The faint-end of the galaxy LF at these redshifts provides powerful constraints on different dark matter models that suppress small-scale structure formation. In this work we perform full hydrodynamical cosmological simulations of galaxy formation using an alternative DM model composed of extremely light bosonic particles ($m \sim 10^{-22}$ eV), also known as fuzzy dark matter (FDM), and examine the predictions for the galaxy stellar mass function and luminosity function at $z \geq 6$ for a range of FDM masses. We find that for FDM models with bosonic mass $m = 5\times10^{-22}$ eV, the number density of galaxies with stellar mass $\rm M_* \sim 10^7 M_{\odot}$ is suppressed by $\sim 40\%$ at z = 9, $\sim 20\%$ at z = 5, and the UV LFs within magnitude range of -16 < $M_{\rm UV}$ < -14 is suppressed by $\sim 60\%$ at $z = 9$, $\sim 20\%$ at $z = 5$ comparing to the CDM counterpart simulation. Comparing our predictions with current measurements of the faint-end LFs ($-18 \leqslant M_{\rm UV} \leqslant -14$), we find that FDM models with $m_{22} < 5\times10^{-22}$ are ruled out at $3\sigma$ confidence level. We expect that future LF measurements by James Webb Space Telescope (JWST), which will extend down to $M_{\rm UV} \sim -13$ for $z \lesssim 10$, with a survey volume that is comparable to the Hubble Ultra Deep Field (HUDF) would have the capability to constrain FDM models to $m\; \gtrsim 10^{-21}$ eV.

Provided that sufficient resources are deployed, we can look forward to an extraordinary future in which we will characterize potentially habitable planets. Until now, we have had to base interpretations of observations on habitability hypotheses that have remained untested. To test these theories observationally, we propose a statistical comparative planetology approach to questions of planetary habitability. The key objective of this approach will be to make quick and cheap measurements of critical planetary characteristics on a large sample of exoplanets, exploiting statistical marginalization to answer broad habitability questions. This relaxes the requirement of obtaining multiple types of data for a given planet, as it allows us to test a given hypothesis from only one type of measurement using the power of an ensemble. This approach contrasts with a "systems science" approach, where a few planets would be extensively studied with many types of measurements. A systems science approach is associated with a number of difficulties which may limit overall scientific return, including: the limited spectral coverage and noise of instruments, the diversity of exoplanets, and the extensive list of potential false negatives and false positives. A statistical approach could also be complementary to a systems science framework by providing context to interpret extensive measurements on planets of particular interest. We strongly recommend future missions with a focus on exoplanet characterization, and with the capability to study large numbers of planets in a homogenous way, rather than exclusively small, intense studies directed at a small sample of planets.

Understanding the sources of lunar water is crucial for studying the history of lunar evolution, and also the solar wind interaction with the Moon and other airless bodies. Recent observations revealed lunar hydration is very likely a surficial dynamic process driven by solar wind. Solar wind is shielded over a period of 3-5 days as the Moon passes through the Earth's magnetosphere, during which a significant loss of hydration is expected from previous works.Here we study lunar hydration inside the magnetosphere using orbital spectral data, which unexpectedly found that the polar surficial OH/H2O abundance remains at the same level when in the solar wind and in the magnetosphere. We suggest that particles from the magnetosphere (Earth wind, naturally different from solar wind) contribute to lunar hydration. From lunar orbital plasma observations, we find the existence of optimal energy ranges, other than 1 keV as previously thought, for surface hydration formation. These optimal energy ranges deduced from space observations may provide strong implications for laboratory experiments simulating lunar hydration processes.

Optimal analyses using the 2-point functions of large-scale structure probes require accurate covariance matrices. A covariance matrix of the 2-point function comprises the disconnected part and the connected part. While the connected covariance only becomes important on small scales, the disconnected covariance is dominant on large scales, where the survey window has a significant impact. In this work, we develop an analytical method to compute the disconnected covariance, accounting for the window effect. Derived under the flat-sky approximation, our formalism is applicable to wide surveys by swapping in the curved-sky window functions. Our method works for both the power spectrum and the correlation function, and applies to the covariances of various probes including the multipoles and the wedges of 3D clustering, the angular and the projected statistics of clustering and shear, as well as the cross covariances between different probes. We verify the analytic covariance against the sample covariance from the galaxy mock simulations in two test cases: (1) the power spectrum multipole covariance, and (2) the joint covariance of the projected correlation function and the correlation function multipoles. Our method achieve good agreement with the mocks, while at a negligible computational cost. Unlike mocks, our analytic covariance is free of sampling noise, which often leads to numerical problems and the need to inflate the errors. In addition, our method can use the best-fit power spectrum as input, in contrast to the standard procedure of using a fiducial model that may deviate significantly from the truth. We also show that a naive diagonal power spectrum covariance underestimates the signal-to-noise ratio compared to our analytic covariance. The code that accompanies this paper is available at https://github.com/eelregit/covdisc.

We present an efficient and accurate method for simulating massive neutrinos in cosmological structure formation simulations, together with an easy to use public implementation. Our method builds on our earlier implementation of the linear response approximation (LRA) for neutrinos, coupled with an N-body code for cold dark matter particles. The LRA's good behaviour at early times and in the linear regime is preserved, while better following the non-linear clustering of neutrinos on small scales. Massive neutrinos are split into initially "fast" and "slow" components. The fast component is followed analytically with the LRA all the way to redshift zero. The slow component is evolved with the LRA only down to a switch-on redshift $z_\nu = 1$, below which it is followed with the particle method, in order to fully account for its non-linear evolution. The slow neutrino particles are initialized at $z = 99$ in order to have accurate positions and velocities at the switch-on time, but are not used to compute the potential until $z \leq 1$, thus avoiding the worst effect of particle shot noise. We show that our hybrid method matches (and for small neutrino masses, exceeds) the accuracy of neutrino particle simulations with substantially lower particle load requirements.

The large scale structure (LSS) of the universe is generated by the linear density gaussian modes, which are evolved into the observed nonlinear LSS. The posterior surface of the modes is convex in the linear regime, leading to a unique global maximum (MAP), but this is no longer guaranteed in the nonlinear regime. In this paper we investigate the nature of posterior surface using the recently developed MAP reconstruction method, with a simplified but realistic N-body simulation as the forward model. The reconstruction method uses optimization with analytic gradients from back-propagation through the simulation. For low noise cases we recover the initial conditions well into the nonlinear regime ($k\sim 1$ h/Mpc) nearly perfectly. We show that the large scale modes can be recovered more precisely than the linear expectation, which we argue is a consequence of nonlinear mode coupling. For noise levels achievable with current and planned LSS surveys the reconstruction cannot recover very small scales due to noise. We see some evidence of non-convexity, specially for smaller scales where the non-injective nature of the mappings: several very different initial conditions leading to the same near perfect final data reconstruction. We investigate the nature of these phenomena further using a 1-d toy gravity model, where many well separated local maximas are found to have identical data likelihood but differ in the prior. We also show that in 1-d the prior favors some solutions over the true solution, though no clear evidence of these in 3-d. Our main conclusion is that on very small scales and for a very low noise the posterior surface is multi-modal and the global maximum may be unreachable with standard methods, while for realistic noise levels in the context of the current and next generation LSS surveys MAP optimization method is likely to be nearly optimal.

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).

MARST (Multi-channel Antarctic Solar Telescope) aims to observe the Sun in multiple wavebands in Antarctica and will be China's first solar telescope in Antarctica. The telescope has two tubes, corresponding to Photosphere observation which uses 11 filters and Chromosphere observation in H{\alpha} waveband. The two tubes need to be coordinated to observe at the same time. The telescope will observe the Sun for a long time, so a self-guiding module is needed to improve sun tracking. Besides, performing solar specific flat-field exposure is necessary for analyzing. EPICS is introduced to control each hardware and an autonomous observation system based on RTS2 is designed under such demands. EPICS application modules are implemented for each device: telescope mount & focuser, filter wheel, H{\alpha} filter, dome with webcams, Andor CCD and PI CCD. We also integrate EPICS modules into RTS2 framework with an XML format configuration. To control these applications autonomously we have developed a RTS2 executor module where two plan classes are instantiated to control two sets of filters and CCDs, and to ensure only one could control the mount at the same time. Different types of observation plans are designed to describe different series of processes with different priorities. To improve sun tracking, we calculate the centroid of each image to get the offset, then apply the correction to the telescope during observation process. For frontend users, a GUI based on PyQt5 and QML is implemented and connects to rts2-httpd and rts2-proxy modules so that users can control devices, check images and get logs.

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.

Future Baryon Acoustic Oscillation surveys aim at observing galaxy clustering over a wide range of redshift and galaxy populations at great precision, reaching tenths of a percent, in order to detect any deviation of dark energy from the $\LCDM$ model. We utilize a set of paired quasi-\Nb\, FastPM simulations that were designed to mitigate the sample variance effect on the BAO feature and evaluated the BAO systematics as precisely as $\sim 0.01\%$. We report anisotropic BAO scale shifts before and after density field reconstruction in the presence of redshift-space distortions over a wide range of redshift, galaxy/halo biases, and shot noise levels. We test different reconstruction schemes and different smoothing filter scales, and introduce physically-motivated BAO fitting models. For the first time, we derive a Galilean-invariant infrared resummed model for halos in real and redshift space. We test these models from the perspective of robust BAO measurements and non-BAO information such as growth rate and nonlinear bias. We find that pre-reconstruction BAO scale has moderate fitting-model dependence at the level of $0.1\%-0.2\%$ for matter while the dependence is substantially reduced to less than $0.07\%$ for halos. We find that post-reconstruction BAO shifts are generally reduced to below $0.1\%$ in the presence of galaxy/halo bias and show much smaller fitting model dependence. Different reconstruction conventions can potentially make a much larger difference on the line-of-sight BAO scale, upto $0.3\%$. Meanwhile, the precision (error) of the BAO measurements is quite consistent regardless of the choice of the fitting model or reconstruction convention.

We present a method to reconstruct the initial conditions of the universe using observed galaxy positions and luminosities under the assumption that the luminosities can be calibrated with weak lensing to give the mean halo mass. Our method relies on following the gradients of forward model and since the standard way to identify halos is non-differentiable and results in a discrete sample of objects, we propose a framework to model the halo position and mass field starting from the non-linear matter field using Neural Networks. We evaluate the performance of our model with multiple metrics. Our model is more than $95\%$ correlated with the halo-mass fields up to $k\sim 0.7 {\rm h/Mpc}$ and significantly reduces the stochasticity over the Poisson shot noise. We develop a data likelihood model that takes our modeling error and intrinsic scatter in the halo mass-light relation into account and show that a displaced log-normal model is a good approximation to it. We optimize over the corresponding loss function to reconstruct the initial density field and develop an annealing procedure to speed up and improve the convergence. We apply the method to halo number densities of $\bar{n} = 2.5\times 10^{-4} -10^{-3}({\rm h/Mpc})^3$, typical of current and future redshift surveys, and recover a Gaussian initial density field, mapping all the higher order information in the data into the power spectrum. We show that our reconstruction improves over the standard reconstruction. For baryonic acoustic oscillations (BAO) the gains are relatively modest because BAO is dominated by large scales where standard reconstruction suffices. We improve upon it by $\sim 15-20\%$ in terms of error on BAO peak as estimated by Fisher analysis at $z=0$. We expect larger gains will be achieved when applying this method to the broadband linear power spectrum reconstruction on smaller scales.

W Ursa Majoris (W UMa)-type contact binary systems (CBs) are useful statistical distance indicators because of their large numbers. Here, we establish (orbital) period-luminosity relations (PLRs) in 12 optical-to-mid-infrared bands (GBV RIJHKsW1W2W3W4) based on 183 nearby W UMa-type CBs with accurate Tycho-Gaia parallaxes. The 1{\sigma} dispersion of the PLRs decreases from optical to near- and mid-infrared wavelengths. The minimum scatter, 0.16 mag, implies that W UMa-type CBs can be used to recover distances to 7% precision. Applying our newly determined PLRs to 19 open clusters containing W UMa-type CBs demonstrates that the PLR and open cluster CB distance scales are mutually consistent to within 1%. Adopting our PLRs as secondary distance indicators, we compiled a catalog of 55,603 CBs candidates, of which 80% have distance estimates based on a combination of optical, near-, and mid-infrared photometry. Using Fourier decomposition, 27,318 high-probability W UMa-type CBs were selected. The resulting 8% distance accuracy implies that our sam- ple encompasses the largest number of objects with accurate distances within a local volume with a radius of 3 kpc available to date. The distribution of W UMa-type CBs in the Galaxy suggests that in different environments, the CB luminosity function may be different: larger numbers of brighter (longer-period) W UMa-type CBs are found in younger environments.

Fast N-body PM simulations with a small number of time steps such as FastPM or COLA have been remarkably successful in modeling the galaxy statistics, but their lack of small scale force resolution and long time steps cannot give accurate halo matter profiles or matter power spectrum. High resolution N-body simulations can improve on this, but lack baryonic effects, which can only be properly included in hydro simulations. Here we present a scheme to calibrate the fast simulations to mimic the precision of the hydrodynamic simulations or high resolution N-body simulations. The scheme is based on a gradient descent of either effective gravitational potential, which mimics the short range force, or of effective enthalpy, which mimics gas hydrodynamics and feedback. The scheme is fast and differentiable, and can be incorporated as a post-processing step into any simulation. It gives very good results for the matter power spectrum for several of the baryonic feedback and dark matter simulations, and also gives improved dark matter halo profiles. The scheme is even able to find the large subhalos, and increase the correlation coefficient between the fast simulations and the high resolution N-body or hydro simulations. It can also be used to add baryonic effects to the high resolution N-body simulations. While the method has free parameters that can be calibrated on various simulations, they can also be viewed as astrophysical nuisance parameters describing baryonic effects that can be marginalized over during the data analysis. In this view these parameters can be viewed as an efficient parametrization of baryonic effects.

X-ray-dim isolated neutron stars (XDINSs), also known as the Magnificent Seven, exhibits a Planck-like soft X-ray spectrum. In the optical/ultraviolet(UV) band, there is an excess of radiation compared to the extrapolation from the X-ray spectrum. A model of bremsstrahlung emission from a nonuniform plasma atmosphere is proposed in the regime of a strangeon star to explain the optical/UV excess and its spectral deviation. The nonuniform atmosphere is formed by the accretion of interstellar matter which moves along the magnetic field lines to polar caps, and these particles may spread out of the pole regions. The modeled electron temperatures are $\sim100-200$ eV with radiation radii $R_{\rm opt}^{\infty}\sim5-14$ km. The spectra of five sources (RX J0720.4--3125, RX J0806.4--4123, RX J1308.6+2127, RX J1605.3+3249, RX J1856.5--3754) from optical/UV to X-ray bands could be well fitted by the radiative model, and exhibit gaussian absorption lines at $\sim 100-500$ eV as would be expected. Furthermore, the surroundings (i.e., fallback disks or dusty belts) of XDINSs could be tested by future observations with the Large Optical Telescope (LOT).

We develop a new method to constraint primordial non-Gaussianities of the local kind using unclustered tracers of the Large Scale Structure. We show that in the limit of low noise, zero bias tracers yield large improvement over standard methods, mostly due to vanishing sampling variance. We propose a simple technique to construct such a tracer, using environmental information obtained from the original sample, and validate our method with N-body simulations. Our results indicate that $\sigma_{f_{\mathrm{NL}}^{\mathrm{loc}}}\simeq1$ can be reached using only information on a single tracer of sufficiently high number density.

Space-based high contrast imaging mission concepts for studying rocky exoplanets in reflected light are currently under community study. We develop an inverse modeling framework to estimate the science return of such missions given different instrument design considerations. By combining an exoplanet albedo model, an instrument noise model, and an ensemble Markov chain Monte Carlo sampler, we explore retrievals of atmospheric and planetary properties for Earth twins as a function of signal-to-noise ratio (SNR) and resolution ($R$). Our forward model includes Rayleigh scattering, single-layer water clouds with patchy coverage, and pressure-dependent absorption due to water vapor, oxygen, and ozone. We simulate data at $R = 70$ and $R = 140$ from 0.4-1.0 $\mu$m with SNR $ = 5, 10, 15, 20$ at 550 nm (i.e., for HabEx/LUVOIR-type instruments). At these same SNR, we simulate data for WFIRST paired with a starshade, which includes two photometric points between 0.48-0.6 $\mu$m and $R = 50$ spectroscopy from 0.6-0.97 $\mu$m. Given our noise model for WFIRST-type detectors, we find that weak detections of water vapor, ozone, and oxygen can be achieved with observations with at least $R = 70$ / SNR$\ = 15$, or $R = 140$ / SNR$\ = 10$ for improved detections. Meaningful constraints are only achieved with $R = 140$ / SNR$\ = 20$ data. The WFIRST data offer limited diagnostic information, needing at least SNR = 20 to weakly detect gases. Most scenarios place limits on planetary radius, but cannot constrain surface gravity and, thus, planetary mass.

One of the main unsolved problems of cosmology is how to maximize the extraction of information from nonlinear data. If the data are nonlinear the usual approach is to employ a sequence of statistics (N-point statistics, counting statistics of clusters, density peaks or voids etc.), along with the corresponding covariance matrices. However, this approach is computationally prohibitive and has not been shown to be exhaustive in terms of information content. Here we instead develop a Bayesian approach, expanding the likelihood around the maximum posterior of linear modes, which we solve for using optimization methods. By integrating out the modes using perturbative expansion of the likelihood we construct an initial power spectrum estimator, which for a fixed forward model contains all the cosmological information if the initial modes are gaussian distributed. We develop a method to construct the window and covariance matrix such that the estimator is explicitly unbiased and nearly optimal. We then generalize the method to include the forward model parameters, including cosmological and nuisance parameters, and primordial non-gaussianity. We apply the method in the simplified context of nonlinear structure formation, using either simplified 2-LPT dynamics or N-body simulations as the nonlinear mapping between linear and nonlinear density, and 2-LPT dynamics in the optimization steps used to reconstruct the initial density modes. We demonstrate that the method gives an unbiased estimator of the initial power spectrum, providing among other a near optimal reconstruction of linear baryonic acoustic oscillations.

Observations of the 21 cm line from neutral hydrogen indicate that an Epoch of Heating (EoH) might have preceded the later Epoch of Reionization (EoR). Here we study the effects on the ionization state and the thermal history of the Intergalactic Medium (IGM) during the EoH induced by different assumptions on ionizing sources in the high redshift Universe: (i) stars, (ii) X-ray binaries (XRBs), (iii) thermal bremsstrahlung of the hot Interstellar Medium (ISM), and (iv) accreting nuclear black holes (BHs). To this aim, we post-process outputs from the ($100 h^{-1}$ cMpc)$^3$ hydrodynamical simulation MassiveBlack-II with the cosmological 3D radiative transfer code CRASH, which follows the propagation of UV and X-ray photons, computing the thermal and ionization state of hydrogen and helium through the EoH. We find that stars determine the fully ionized morphology of the IGM, while the spectrally hard XRBs pave way for efficient subsequent heating and ionization by the spectrally softer ISM. With the seeding prescription in MassiveBlack-II, BHs do not contribute significantly to either ionization or heating. With only stars, most of the IGM remains in a cold state (with a median $T=11$ K at $z=10$), however, the presence of more energetic sources raises the temperature of regions around the brightest and more clustered sources above that of the CMB, opening the possibility to observing the 21 cm signal in emission.

The scaling relations between the mass of supermassive black holes ($M_{\bullet}$) and host galaxy properties (stellar mass, $M_{\star}$, and velocity dispersion, $\sigma$), provide a link between the growth of black holes (BHs) and that of their hosts. Here we investigate if and how the BH-galaxy relations are established in the high-$z$ universe using \textsc{BlueTides}, a high-resolution large volume cosmological hydrodynamic simulation. We find the $M_{\bullet}-M_{\star}$ and $M_{\bullet}-\sigma$ relations at $z=8$: $\log_{10}(M_{\bullet}) = 8.25 + 1.10 \ \log_{10}(M_{\star}/10^{11}M_{\odot})$ and $\log_{10}(M_{\bullet}) = 8.35 + 5.31 \ \log_{10}(\sigma/200kms^{-1})$ at $z=8$, both fully consistent with the local measurements. The slope of the $M_{\bullet}-\sigma$ relation is slightly steeper for high star formation rate and $M_{\star}$ galaxies while it remains unchanged as a function of Eddington accretion rate onto the BH. The intrinsic scatter in $M_{\bullet}-\sigma$ relation in all cases ($\epsilon \sim 0.4$) is larger at these redshifts than inferred from observations and larger than in $M_{\bullet}-M_{\star}$ relation ($\epsilon \sim 0.14$). We find the gas-to-stellar ratio $f=M_{\rm gas}/M_{\star}$ in the host (which can be very high at these redshifts) to have the most significant impact setting the intrinsic scatter of $M_{\bullet}-\sigma$. The scatter is significantly reduced when galaxies with high gas fractions ($\epsilon = 0.28$ as $f<10$) are excluded (making the sample more comparable to low-$z$ galaxies); these systems have the largest star formation rates and black hole accretion rates, indicating that these fast-growing systems are still moving toward the relation at these high redshifts. Examining the evolution (from $z=10$ to 8) of high mass black holes in $M_{\bullet}-\sigma$ plane confirms this trend.

We study the clustering of the highest-z galaxies (from ~ $0.1$ to a few tens Mpc scales) using the BLUETIDES simulation and compare it to current observational constraints from Hubble legacy and Hyper Suprime Cam (HSC) fields (at $z=6-7.2$). With a box length of $400$ $Mpc/h$ on each side and $0.7$ trillion particles, BLUETIDES is the largest high resolution cosmological hydrodynamic simulation to date ideally suited for studies of high-z galaxies. We find that galaxies with magnitude $m_{UV}<27.7$ have a bias ($b_g$) of $8.1\pm 1.2$ at $z=8$, and typical halo masses $M_H \gtrsim 6\times10^{10} M_{\odot}$. Given the redshift evolution between $z=8$ to $z=10$ ($b_g\propto(1+z)^{1.6}$), our inferred values of the bias and halo masses are consistent with measured angular clustering at $z \sim 6.8$ from these brighter samples. The bias of fainter galaxies (in the Hubble legacy field at $H_{160} \lesssim29.5$) is $5.9\pm0.9$ at $z=8$ corresponding to halo masses $M_H \gtrsim 10^{10} M_{\odot}$. We investigate directly the 1-halo term inthe clustering and show that it dominates on scales $r \lesssim 0.1$ Mpc/$h$ ($\Theta \lesssim 3"$) with non-linear effect at transition scales between the 1-halo and 2-halo term affecting scales 0.1 $\lesssim r \lesssim $ 20 Mpc/$h$ ($3"\lesssim \Theta \lesssim 90"$). Current clustering measurements probe down to the scales in the transition between 1-halo to 2-halo regime where non-linear effects are important. The amplitude of the 1-halo term implies that occupation numbers for satellites in \texttt{BLUETIDES} are somewhat higher than standard HODs adopted in these analyses (which predict amplitudes in the 1-halo regime suppressed by a factor 2-3). That possibly implies a higher number of galaxies detected by JWST (at small scales and even fainter magnitudes) observing these fields.

The hybrid Tibet AS array was successfully constructed in 2014. It has 4500 m$^{2}$ underground water Cherenkov pools used as the muon detector (MD) and 789 scintillator detectors covering 36900 m$^{2}$ as the surface array. At 100 TeV, cosmic-ray background events can be rejected by approximately 99.99\%, according to the full Monte Carlo (MC) simulation for $\gamma$-ray observations. In order to use the muon detector efficiently, we propose to extend the surface array area to 72900 m$^{2}$ by adding 120 scintillator detectors around the current array to increase the effective detection area. A new prototype scintillator detector is developed via optimizing the detector geometry and its optical surface, by selecting the reflective material and adopting dynode readout. This detector can meet our physics requirements with a positional non-uniformity of the output charge within 10\% (with reference to the center of the scintillator), time resolution FWHM of $\sim$2.2 ns, and dynamic range from 1 to 500 minimum ionization particles.

We present nbodykit, an open-source, massively parallel Python toolkit for analyzing large-scale structure (LSS) data. Using Python bindings of the Message Passing Interface (MPI), we provide parallel implementations of many commonly used algorithms in LSS. nbodykit is both an interactive and scalable piece of scientific software, performing well in a supercomputing environment while still taking advantage of the interactive tools provided by the Python ecosystem. Existing functionality includes estimators of the power spectrum, 2 and 3-point correlation functions, a Friends-of-Friends grouping algorithm, mock catalog creation via the halo occupation distribution technique, and approximate N-body simulations via the FastPM scheme. The package also provides a set of distributed data containers, insulated from the algorithms themselves, that enable nbodykit to provide a unified treatment of both simulation and observational data sets. nbodykit can be easily deployed in a high performance computing environment, overcoming some of the traditional difficulties of using Python on supercomputers. We provide performance benchmarks illustrating the scalability of the software. The modular, component-based approach of nbodykit allows researchers to easily build complex applications using its tools. The package is extensively documented at http://nbodykit.readthedocs.io, which also includes an interactive set of example recipes for new users to explore. As open-source software, we hope nbodykit provides a common framework for the community to use and develop in confronting the analysis challenges of future LSS surveys.

Motivated by recent observational constraints on dust reprocessed emission in star forming galaxies at $z\sim 6$ and above we use the very-large cosmological hydrodynamical simulation \bluetides\ to explore predictions for the amount of dust obscured star formation in the early Universe ($z>8$). \bluetides\ matches current observational constraints on both the UV luminosity function and galaxy stellar mass function and predicts that approximately $90\%$ of the star formation in high-mass ($M_{*}>10^{10}\,{\rm M_{\odot}}$) galaxies at $z=8$ is already obscured by dust. The relationship between dust attenuation and stellar mass predicted by \bluetides\ is consistent with that observed at lower redshift. However, observations of several individual objects at $z>6$ are discrepant with the predictions, though it is possible their uncertainties may have been underestimated. We find that the predicted surface density of $z\ge 8$ sub-mm sources is below that accessible to current {\em Herschel}, SCUBA-2, and ALMA sub-mm surveys. However, as ALMA continues to accrue additional surface area the population of $z>8$ dust-obscured galaxies may become accessible in the near future.

We introduce FastPM, a highly-scalable approximated particle mesh N-body solver, which implements the particle mesh (PM) scheme enforcing correct linear displacement (1LPT) evolution via modified kick and drift factors. Employing a 2-dimensional domain decomposing scheme, FastPM scales extremely well with a very large number of CPUs. In contrast to COmoving-LAgrangian (COLA) approach, we do not require to split the force or track separately the 2LPT solution, reducing the code complexity and memory requirements. We compare FastPM with different number of steps ($N_s$) and force resolution factor ($B$) against 3 benchmarks: halo mass function from Friends of Friends halo finder, halo and dark matter power spectrum, and cross correlation coefficient (or stochasticity), relative to a high resolution TreePM simulation. We show that the modified time stepping scheme reduces the halo stochasticity when compared to COLA with the same number of steps and force resolution. While increasing $N_s$ and $B$ improves the transfer function and cross correlation coefficient, for many applications FastPM achieves sufficient accuracy at low $N_s$ and $B$. For example, $N_s=10$ and $B=2$ simulation provides a substantial saving (a factor of 10) of computing time relative to $N_s=40$, $B=3$ simulation, yet the halo benchmarks are very similar at $z=0$. We find that for abundance matched halos the stochasticity remains low even for $N_s=5$. FastPM compares well against less expensive schemes, being only 7 (4) times more expensive than 2LPT initial condition generator for $N_s=10$ ($N_s=5$). Some of the applications where FastPM can be useful are generating a large number of mocks, producing non-linear statistics where one varies a large number of nuisance or cosmological parameters, or serving as part of an initial conditions solver.

Reliable extraction of cosmological information from clustering measurements of galaxy surveys requires estimation of the error covariance matrices of observables. The accuracy of covariance matrices is limited by our ability to generate sufficiently large number of independent mock catalogs that can describe the physics of galaxy clustering across a wide range of scales. Furthermore, galaxy mock catalogs are required to study systematics in galaxy surveys and to test analysis tools. In this investigation, we present a fast and accurate approach for generation of mock catalogs for the upcoming galaxy surveys. Our method relies on low-resolution approximate gravity solvers to simulate the large scale dark matter field, which we then populate with halos according to a flexible nonlinear and stochastic bias model. In particular, we extend the \textsc{patchy} code with an efficient particle mesh algorithm to simulate the dark matter field (the \textsc{FastPM} code), and with a robust MCMC method relying on the \textsc{emcee} code for constraining the parameters of the bias model. Using the halos in the BigMultiDark high-resolution $N$-body simulation as a reference catalog, we demonstrate that our technique can model the bivariate probability distribution function (counts-in-cells), power spectrum, and bispectrum of halos in the reference catalog. Specifically, we show that the new ingredients permit us to reach percentage accuracy in the power spectrum up to $k\sim 0.4\; \,h\,{\rm Mpc}^{-1}$ (within 5\% up to $k\sim 0.6\; \,h\,{\rm Mpc}^{-1}$) with accurate bispectra improving previous results based on Lagrangian perturbation theory.

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}$.

We describe a simple and fast algorithm for identifying friends-of-friends features and prove its correctness. The algorithm avoids unnecessary expensive neighbor queries, uses minimal memory overhead, and rejects slowdown in high over-density regions. We define our algorithm formally based on pair enumeration, a problem that has been heavily studied in fast 2-point correlation codes and our reference implementation employs a dual KD-tree correlation function code. We construct features in a hierarchical tree structure, and use a splay operation to reduce the average cost of identifying the root of a feature from $O[\log L]$ to $O[1]$ ($L$ is the size of a feature) without additional memory costs. This reduces the overall time complexity of merging trees from $O[L\log L]$ to $O[L]$, reducing the number of operations per splay by orders of magnitude. We next introduce a pruning operation that skips merge operations between two fully self-connected KD-tree nodes. This improves the robustness of the algorithm, reducing the number of merge operations in high density peaks from $O[\delta^2]$ to $O[\delta]$. We show that for cosmological data set the algorithm eliminates more than half of merge operations for typically used linking lengths $b \sim 0.2$ (relative to mean separation). Furthermore, our algorithm is extremely simple and easy to implement on top of an existing pair enumeration code, reusing the optimization effort that has been invested in fast correlation function codes.

We employ the very large cosmological hydrodynamical simulation BLUETIDES to investigate the predicted properties of the galaxy population during the epoch of reionisation ($z>8$). BLUETIDES has a resolution and volume ($(400/h\approx 577)^{3}\,{\rm cMpc^3}$) providing a population of galaxies which is well matched to depth and area of current observational surveys targeting the high-redshift Universe. At $z=8$ BLUETIDES includes almost 160,000 galaxies with stellar masses $>10^{8}\,{\rm M_{\odot}}$. The population of galaxies predicted by BLUETIDES closely matches observational constraints on both the galaxy stellar mass function and far-UV ($150\,{\rm nm}$) luminosity function. Galaxies in BLUETIDES are characterised by rapidly increasing star formation histories. Specific star formation rates decrease with redshift though remain largely insensitive to stellar mass. As a result of the enhanced surface density of metals more massive galaxies are predicted to have higher dust attenuation resulting in a significant steepening of the observed far-UV luminosity function at high luminosities. The contribution of active SMBHs to the UV luminosities of galaxies with stellar masses $10^{9-10}\,{\rm M_{\odot}}$ is around $3\%$ on average. Approximately $25\%$ of galaxies with $M_{*}\approx 10^{10}\,{\rm M_{\odot}}$ are predicted to have active SMBH which contribute $>10\%$ of the total UV luminosity.

Previous measurements of heat redistribution efficiency (the ability to transport energy from a planet's highly irradiated dayside to its eternally dark nightside) show considerable variation between exoplanets. Theoretical models predict a positive correlation between heat redistribution efficiency and temperature for tidally locked planets; however, recent HST WASP-43b spectroscopic phase curve results are inconsistent with current predictions. Using the Spitzer Space Telescope, we obtained a total of three phase curve observations of WASP-43b (P=0.813 days) at 3.6 and 4.5 microns. The first 3.6 micron visit exhibits spurious nightside emission that requires invoking unphysical conditions in our cloud-free atmospheric retrievals. The two other visits exhibit strong day-night contrasts that are consistent with the HST data. To reconcile the departure from theoretical predictions, WASP-43b would need to have a high-altitude, nightside cloud/haze layer blocking its thermal emission. Clouds/hazes could be produced within the planet's cool, nearly retrograde mid-latitude flows before dispersing across its nightside at high altitudes. Since mid-latitude flows only materialize in fast-rotating ($\lesssim1$ day) planets, this may explain an observed trend connecting measured day-night contrast with planet rotation rate that matches all current Spitzer phase curve results. Combining independent planetary emission measurements from multiple phases, we obtain a precise dayside hemisphere H2O abundance ($2.5\times 10^{-5} - 1.1\times 10^{-4}$ at 1$\sigma$ confidence) and, assuming chemical equilibrium and a scaled solar abundance pattern, we derive a corresponding metallicity estimate that is consistent with being solar (0.4 -- 1.7). Using the retrieved global CO+CO2 abundance under the same assumptions, we estimate a comparable metallicity of 0.3 - 1.7$\times$ solar. (Abridged)

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

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.

The determination of atmospheric structure and molecular abundances of planetary atmospheres via spectroscopy involves direct comparisons between models and data. While varying in sophistication, most model-spectra comparisons fundamentally assume "1D" model physics. However, knowledge from general circulation models and of solar system planets suggests that planetary atmospheres are inherently "3D" in their structure and composition. We explore the potential biases resulting from standard "1D" assumptions within a Bayesian atmospheric retrieval framework. Specifically, we show how the assumption of a single 1-dimensional thermal profile can bias our interpretation of the thermal emission spectrum of a hot Jupiter atmosphere that is composed of two thermal profiles. We retrieve upon spectra of unresolved model planets as observed with a combination of $HST$ WFC3+$Spitzer$ IRAC as well as $JWST$ under varying differences in the two thermal profiles. For WFC3+IRAC, there is a significantly biased estimate of CH$_4$ abundance using a 1D model when the contrast is 80%. For $JWST$, two thermal profiles are required to adequately interpret the data and estimate the abundances when contrast is greater than 40%. We also apply this preliminary concept to the recent WFC3+IRAC phase curve data of the hot Jupiter WASP-43b. We see similar behavior as present in our simulated data: while the H$_2$O abundance determination is robust, CH$_4$ is artificially well-constrained to incorrect values under the 1D assumption. Our work demonstrates the need to evaluate model assumptions in order to extract meaningful constraints from atmospheric spectra and motivates exploration of optimal observational setups.

The growth of the most massive black holes in the early universe, consistent with the detection of highly luminous quasars at $z> 6$ implies sustained, critical accretion of material to grow and power them. Given a black hole seed scenario, it is still uncertain which conditions in the early Universe allow the fastest black hole growth. Large scale hydrodynamical cosmological simulations of structure formation allow us to explore the conditions conducive to the growth of the earliest supermassive black holes. We use the cosmological hydrodynamic simulation BlueTides, which incorporates a variety of baryon physics in a (400 Mpc/h)^3 volume with 0.7 trillion particles to follow the earliest phases of black hole critical growth. At z=8 the most massive black holes (a handful) approach masses of 10^8 Msun with the most massive (with M_BH = 4 x 10^8 Msun ) being found in an extremely compact spheroid-dominated host galaxy. Examining the large-scale environment of hosts, we find that the initial tidal field is more important than overdensity in setting the conditions for early BH growth. In regions of low tidal fields gas accretes 'cold' onto the black hole and falls along thin, radial filaments straight into the center forming the most compact galaxies and most massive black holes at earliest times. Regions of high tidal fields instead induce larger coherent angular momenta and influence the formation of the first population of massive compact disks. The extreme early growth depends on the early interplay of high gas densities and the tidal field that shapes the mode of accretion. Mergers play a minor role in the formation of the first generation, rare massive BHs.

We use the BlueTides simulation to predict the properties of the high-$z$ galaxy and active galactic nuclei (AGN) populations for the planned 2200deg$^2$ Wide-Field Infrared Survey Telescope's (WFIRST)-AFTA High Latitude Survey (HLS). BlueTides is a cosmological hydrodynamic simulation, which incorporates a variety of baryon physics in a $(400h^{-1} \mathrm{Mpc})^3$ volume evolved to $z=8$ with 0.7 trillion particles. The galaxy luminosity functions in the simulation show good agreement with all the current observational constraints (up to $z=11$) and predicts an enhanced number of UV bright galaxies. At the proposed depth of the HLS ($m < 26.75$), BlueTides predicts $10^6$ galaxies at $z=8$ with a few up to $z\sim 15$ due to the enhanced bright end of the galaxy luminosity function. At $z=8$, galaxies in the mock HLS have specific star formation rates of $\sim 10 {\rm Gyr}^{-1}$ and ages of $\sim 80 {\rm Myr}$ (both evolving linearly with redshift) and a non-evolving mass-metallicity relation. BlueTides also predicts $\sim 10^4$ AGN in WFIRST HLS from $z=8$ out to $z\sim 14$. These AGN host black holes of $M\sim 10^6-10^8 M_\odot$ accreting close to their Eddington luminosity. Galaxies and AGN have host halo masses of $M_{halo}\sim 10^{11-12} M_\odot$ and a linear bias $b\approx 13-20$. Given the expected galaxy space densities, their high bias and large volume probed we speculate that it may be feasible for WFIRST HLS detect the Baryon Acoustic Oscillation peak in the galaxy power spectrum out to $z=8-9$.

We use the large cosmological hydro-dynamic simulation BlueTides to predict the photometric properties of galaxies during the epoch of reionisation ($z=8-15$). These properties include the rest-frame UV to near-IR broadband spectral energy distributions, the Lyman continuum photon production, the UV star formation rate calibration, and intrinsic UV continuum slope. In particular we focus on exploring the effect of various modelling assumptions, including the assumed choice of stellar population synthesis model, initial mass function, and the escape fraction of Lyman continuum photons, upon these quantities. We find that these modelling assumptions can have a dramatic effect on photometric properties leading to consequences for the accurate determination of physical properties from observations. For example, at $z=8$ we predict that nebular emission can account for up-to $50\%$ of the rest-frame $R$-band luminosity, while the choice of stellar population synthesis model can change the Lyman continuum production rate up to a factor of $\times 2$.

We explore the relationship between the nonlinear matter power spectrum and the various Lagrangian and Standard Perturbation Theories (LPT and SPT). We first look at it in the context of one dimensional (1-d) dynamics, where 1LPT is exact at the perturbative level and one can exactly resum the SPT series into the 1LPT power spectrum. Shell crossings lead to non-perturbative effects, and the PT ignorance can be quantified in terms of their ratio, which is also the transfer function squared in the absence of stochasticity. At the order of PT we work, this parametrization is equivalent to the results of effective field theory (EFT), and can thus be expanded in terms of the same parameters. We find that its radius of convergence is larger than the SPT loop expansion. The same EFT parametrization applies to all SPT loop terms and, if stochasticity can be ignored, to all N-point correlators. In 3-d, the LPT structure is considerably more complicated, and we find that LPT models with parametrization motivated by the EFT exhibit running with $k$ and that SPT is generally a better choice. Since these transfer function expansions contain free parameters that change with cosmological model their usefulness for broadband power is unclear. For this reason we test the predictions of these models on baryonic acoustic oscillations (BAO) and other primordial oscillations, including string monodromy models, for which we ran a series of simulations with and without oscillations. Most models are successful in predicting oscillations beyond their corresponding PT versions, confirming the basic validity of the model.

Using deep Hubble and Spitzer observations Oesch et al. (2016) have identified a bright ($M_{\rm UV}\approx -22$) star forming galaxy candidate at $z \approx 11$. The presence of GN-$z11$ implies a number density $\sim 10^{-6}\,{\rm Mpc^{-3}}$, roughly an order of magnitude higher than the expected value based on extrapolations from lower redshift. Using the unprecedented volume and high resolution of the BlueTides cosmological hydrodynamical simulation, we study the population of luminous rare objects at $z > 10$. The luminosity function in BlueTides implies an enhanced number of massive galaxies, consistent with the observation of GN-$z11$. We find about 30 galaxies at $M_{\rm UV}\approx -22$ at $z = 11$ in the BlueTides volume, including a few objects about 1.5 magnitudes brighter. The probability of observing GN-$z11$ in the volume probed by Oesch et al. (2016) is $\sim 13$ per cent. The predicted properties of the rare bright galaxies at $z = 11$ in BlueTides closely match those inferred from the observations of GN-$z11$. BlueTides predicts a negligible contribution from faint AGN in the observed SED. The enormous increase in volume surveyed by WFIRST will provide observations of $\sim1000$ galaxies with $M_{\rm UV} < -22$ beyond $z = 11$ out to $z = 13.5$.

The Lyman Continuum photon production efficiency ($\xi_{\rm ion}$) is a critical ingredient for inferring the number of photons available to reionise the intergalactic medium. To estimate the theoretical production efficiency in the high-redshift Universe we couple the BlueTides cosmological hydrodynamical simulation with a range of stellar population synthesis models. We find Lyman Continuum photon production efficiencies of $\log_{10}(\xi_{\rm ion}/{\rm erg^{-1}\, Hz})\approx 25.1-25.5$ depending on the choice of stellar population synthesis model. These results are broadly consistent with recent observational constraints at high-redshift though favour a model incorporating the effects of binary evolution

As galaxy surveys begin to measure the imprint of baryonic acoustic oscillations (BAO) on large-scale structure at the sub-percent level, reconstruction techniques that reduce the contamination from nonlinear clustering become increasingly important. Inverting the nonlinear continuity equation, we propose an Eulerian growth-shift reconstruction algorithm that does not require the displacement of any objects, which is needed for the standard Lagrangian BAO reconstruction algorithm. In real-space DM-only simulations the algorithm yields 95% of the BAO signal-to-noise obtained from standard reconstruction. The reconstructed power spectrum is obtained by adding specific simple 3- and 4-point statistics to the pre-reconstruction power spectrum, making it very transparent how additional BAO information from higher-point statistics is included in the power spectrum through the reconstruction process. Analytical models of the reconstructed density for the two algorithms agree at second order. Based on similar modeling efforts, we introduce four additional reconstruction algorithms and discuss their performance.

We describe mock data-sets generated to simulate the high-redshift quasar sample in Data Release 11 (DR11) of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). The mock spectra contain Ly{\alpha} forest correlations useful for studying the 3D correlation function including Baryon Acoustic Oscillations (BAO). They also include astrophysical effects such as quasar continuum diversity and high-density absorbers, instrumental effects such as noise and spectral resolution, as well as imperfections introduced by the SDSS pipeline treatment of the raw data. The Ly{\alpha} forest BAO analysis of the BOSS collaboration, described in Delubac et al. 2014, has used these mock data-sets to develop and cross-check analysis procedures prior to performing the BAO analysis on real data, and for continued systematic cross checks. Tests presented here show that the simulations reproduce sufficiently well important characteristics of real spectra. These mock data-sets will be made available together with the data at the time of the Data Release 11.

In gamma ray astronomy, the energy range from sub-100GeV to TeV is crucial due to where there is a gap between space experiments and ground-based ones. In addition, observations in this energy range are expected to provide more details about the high energy emission from GRBs,and thus to understand EBL. Based on the observation results and the related knowledge, scientists may be able to unveil the mysteries of galaxy formation and the evolution of early universe. One of the principal issues for next generation Imaging Atmospheric Cherenkov Telescopes (IACT) is to achieve larger field of view (FoV). In this work, we report a refractive water convex lens as light collector to test the feasibility of a new generation of IACT, and some preliminary test results on the optical properties (the focal length, spot size, transmittance, etc.) of a 0.9 m diameter water lens, the photodetectors and DAQ system of a prototype are presented and discussed.

We examine the luminosity function (LF) of [OII] emission-line galaxies in the high-resolution cosmological simulation MassiveBlack-II (MBII). From the spectral energy distribution of each galaxy, we select a sub-sample of star-forming galaxies at $0.06 \le z \le 3.0$ using the [OII] emission line luminosity L([OII]). We confirm that the specific star formation rate matches that in the GAMA survey. We show that the [OII] LF at z=1.0 from the MBII shows a good agreement with the LFs from several surveys below L([OII])=$10^{43.0}$ erg/s while the low redshifts ($z \le 0.3$) show an excess in the prediction of bright [OII] galaxies, but still displaying a good match with observations below L([OII])=$10^{41.6}$ erg/s. Based on the validity in reproducing the properties of [OII] galaxies at low redshift ($z \le 1$), we forecast the evolution of the [OII] LF at high redshift ($z \le 3$), which can be tested by upcoming surveys such as the HETDEX and DESI. The slopes of the LFs at bright and faint ends range from -3 to -2 showing minima at z=2. The slope of the bright end evolves approximately as 1/(z+1) at z=2 while the faint end evolves as ~3/(z+1) at $0.6 \le z \le 2$. In addition, a similar analysis is applied for the evolution of [OIII] LFs, which is to be explored in the forthcoming survey WFIRST-AFTA. Finally, we show that the auto-correlation function of [OII] and [OIII] emitting galaxies shows a rapid evolution from z=2 to 1.

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 introduce the BlueTides simulation and report initial results for the luminosity functions of the first galaxies and AGN, and their contribution to reionization. BlueTides was run on the BlueWaters cluster at NCSA from $z=99$ to $z=8.0$ and includes 2$\times$7040$^3$ particles in a $400$Mpc/h per side box, making it the largest hydrodynamic simulation ever performed at high redshift. BlueTides includes a pressure-entropy formulation of smoothed particle hydrodynamics, gas cooling, star formation (including molecular hydrogen), black hole growth and models for stellar and AGN feedback processes. The star formation rate density in the simulation is a good match to current observational data at $z\sim 8-10$. We find good agreement between observations and the predicted galaxy luminosity function in the currently observable range $-18\le M_{\mathrm UV} \le -22.5$ with some dust extinction required to match the abundance of brighter objects. BlueTides implements a patchy reionization model that produces a fluctuating UV background. BlueTides predicts number counts for galaxies fainter than current observational limits which are consistent with extrapolating the faint end slope of the luminosity function with a power law index $\alpha\sim -1.8$ at $z\sim 8$ and redshift dependence of $\alpha\sim (1+z)^{-0.4}$. The AGN population has a luminosity function well fit by a power law with a slope $\alpha\sim -2.4$ that compares favourably with the deepest CANDELS-Goods fields. We investigate how these luminosity functions affect the progress of reionization, and find that a high Lyman-$\alpha$ escape fraction ($f_\mathrm{esc} \sim 0.5$) is required if galaxies dominate the ionising photon budget during reionization. Smaller galaxy escape fractions imply a large contribution from faint AGN (down to $M_\mathrm{UV}=-12$) which results in a rapid reionization, disfavoured by current observations.

An on-going effort in the characterization of exoplanetary systems is the accurate determination of host star properties. This effort extends to the relatively bright host stars of planets discovered with the radial velocity method. The Transit Ephemeris Refinement and Monitoring Survey (TERMS) is aiding in these efforts as part of its observational campaign for exoplanet host stars. One of the first known systems is that of 70 Virginis, which harbors a jovian planet in an eccentric orbit. Here we present a complete characterization of this system with a compilation of TERMS photometry, spectroscopy, and interferometry. We provide fundamental properties of the host star through direct interferometric measurements of the radius (1.5\% uncertainty) and through spectroscopic analysis. We combined 59 new Keck HIRES radial velocity measurements with the 169 previously published from the ELODIE, Hamilton, and HIRES spectrographs, to calculate a refined orbital solution and construct a transit ephemeris for the planet. These newly determined system characteristics are used to describe the Habitable Zone of the system with a discussion of possible additional planets and related stability simulations. Finally, we present 19 years of precision robotic photometry that constrain stellar activity and rule out central planetary transits for a Jupiter-radius planet at the 5$\sigma$ level, with reduced significance down to an impact parameter of $b = 0.95$.

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.

Exoplanetary systems closest to the Sun, with the brightest host stars, provide the most favorable opportunities for characterization studies of the host star and their planet(s). The Transit Ephemeris Refinement and Monitoring Survey uses both new radial velocity measurements and photometry in order to greatly improve planetary orbit uncertainties and the fundamental properties of the star, in this case HD 130322. The only companion, HD 130322b, orbits in a relatively circular orbit, e = 0.029 every ~10.7 days. Radial velocity measurements from multiple sources, including 12 unpublished from the Keck I telescope, over the course of ~14 years have reduced the uncertainty in the transit midpoint to ~2 hours. The transit probability for the b-companion is 4.7%, where M_p sin i = 1.15 M_J and a = 0.0925 AU. In this paper, we compile photometric data from the T11 0.8m Automated Photoelectric Telescope at Fairborn Observatory taken over ~14 years, including the constrained transit window, which results in a dispositive null result for both full transit exclusion of HD 130322b to a depth of 0.017 mag and grazing transit exclusion to a depth of ~0.001 mag. Our analysis of the starspot activity via the photometric data reveals a highly accurate stellar rotation period: 26.53 +/-0.70 days. In addition, the brightness of the host with respect to the comparison stars is anti-correlated with the Ca II H and K indices, typical for a young solar-type star.

We present an update to seven stars with long-period planets or planetary candidates using new and archival radial velocities from Keck-HIRES and literature velocities from other telescopes. Our updated analysis better constrains orbital parameters for these planets, four of which are known multi-planet systems. HD 24040 b and HD 183263 c are super-Jupiters with circular orbits and periods longer than 8 yr. We present a previously unseen linear trend in the residuals of HD 66428 indicative on an additional planetary companion. We confirm that GJ 849 is a multi-planet system and find a good orbital solution for the c component: it is a $1 M_{\rm Jup}$ planet in a 15 yr orbit (the longest known for a planet orbiting an M dwarf). We update the HD 74156 double-planet system. We also announce the detection of HD 145934 b, a $2 M_{\rm Jup}$ planet in a 7.5 yr orbit around a giant star. Two of our stars, HD 187123 and HD 217107, at present host the only known examples of systems comprising a hot Jupiter and a planet with a well constrained period $> 5$ yr, and with no evidence of giant planets in between. Our enlargement and improvement of long-period planet parameters will aid future analysis of origins, diversity, and evolution of planetary systems.

We use the high-resolution simulation MassiveBlackII to examine scaling relations between black hole mass (MBH) and host galaxy properties (sigma, M*, and LV), finding good agreement with observational data, especially at the high-mass end. The simulations have less intrinsic scatter than observations, and the MBH-LV correlation has the largest scatter, suggesting it may the the least fundamental of the three relations. We find Gaussian scatter about all three relations, except among the highest mass galaxies, which host more massive black holes. Below z~2 the slopes for the full population remain roughly z-independent, and only steepen by 50% by z~4. The normalization of the sigma, LV relations evolve by 0.3, 0.43 dex, while the MBH correlation does not evolve to at least z~2. Testing for selection biases, we find samples selected by MBH or M* have steeper slopes than randomly selected samples. If unaccounted for, such a selection function would find faster evolution than inferred from a randomly selected sample, as objects at the highend of the relation tend to evolve more rapidly. We find a potential bias among high-LBH subsamples (tending to reside in higher mass galaxies), but these bright-AGN exhibit no intrinsic bias relative to fainter ones in equivalent-mass hosts, nor is there a significant difference between active- and inactive-samples. Finally we characterize the evolution of individual black holes along the scaling planes. Below the local relation, black holes grow faster than their host (72% of black holes 0.3 dex below the mean relation have a MBH-M* trajectory steeper than the local relation), while those above have shallower trajectories (only 14% are steeper than local). Black holes tend to grow faster than their hosts until surpassing the local relation, at which point their growth is suppressed while their hosts continue to grow, returning them to the mean relation.

The Exoplanet Orbit Database (EOD) compiles orbital, transit, host star, and other parameters of robustly detected exoplanets reported in the peer-reviewed literature. The EOD can be navigated through the Exoplanet Data Explorer (EDE) Plotter and Table, available on the World Wide Web at exoplanets.org. The EOD contains data for 1492 confirmed exoplanets as of July 2014. The EOD descends from a table in Butler et al. (2002) and the Catalog of Nearby Exoplanets (Butler et al. 2006), and the first complete documentation for the EOD and the EDE was presented in Wright et al. (2011). In this work, we describe our work since then. We have expanded the scope of the EOD to include secondary eclipse parameters, asymmetric uncertainties, and expanded the EDE to include the sample of over 3000 Kepler Objects of Interest (KOIs), and other real planets without good orbital parameters (such as many of those detected by microlensing and imaging). Users can download the latest version of the entire EOD as a single comma separated value file from the front page of exoplanets.org.

The abundances of elements in stars are a critical clue to their origins. Observed star-to-star variations in logarithmic abundance within an open cluster are typically only $\sim 0.01-0.05$ over many elements, significantly smaller than the variation of $\sim 0.06-0.3$ seen in the interstellar medium from which the stars form. It is unknown why clusters are so homogenous, and whether homogeneity should also prevail in regions of lower star formation efficiency that do not produce bound clusters. Here we report adaptive mesh simulations using passively-advected scalars in order to trace the mixing of chemical elements as star-forming clouds form and collapse. We show that turbulent mixing during cloud assembly naturally produces a stellar abundance scatter at least ~5 times smaller than that in the gas, sufficient to fully explain the observed chemical homogeneity of stars. Moreover, mixing occurs very early, so that regions with efficiencies $\varepsilon \sim 10\%$ are nearly as well-mixed as those with $\varepsilon\sim 50\%$. This implies that even regions that do not form bound clusters are likely to be well-mixed, and enhances the prospects for using chemical tagging to reconstruct dissolved star clusters via their unique chemical signatures.

Using the suite of high-resolution zoom re-simulations of individual haloes by Martig et al., and the large-scale simulation \emph{MassiveBlack-II}, we examine the differences in measured galaxy properties from techniques with various aperture definitions of where galaxies end. We perform techniques popular in the literature and present a new technique of our own, where the aperture radius is based on the baryonic mass profiles of simulated (sub)haloes. For the average Milky-Way-mass system, we find the two most popular techniques in the literature return differences of order 30 per cent for stellar mass, a factor of 3 for gas mass, 40 per cent for star formation rate, and factors of several for gas accretion and ejection rates. Individual cases can show variations greater than this, with the severity dependent on the concentration of a given system. The average difference in integrated properties for a more general galaxy population are not as striking, but are still significant for stellar and gas mass, especially for optical-limit apertures. The large differences that can occur are problematic for comparing results from various publications. We stress the importance of both defining and justifying a technique choice and discourage using popular apertures that use an exact fraction of the virial radius, due to the unignorable variation in galaxy-to-(sub)halo size. Finally, we note that technique choice does not greatly affect simulated galaxies from lying within the scatter of observed scaling relations, but it can alter the derived best-fit slope for the Kennicutt-Schmidt relation.

A direct WIMP (Weakly Interacting Massive Particle) detector with a neutron veto system is designed to better reject neutrons. An experimental configuration is studied in the present paper: a WIMP detectors with CsI(Na) target is placed inside a reactor neutrino detector. The neutrino detector is used as a neutron veto device. The neutron background for the experimental design has been estimated using the Geant4 simulation. The results show that the neutron background can decrease to O(0.01) events per year per tonne of CsI(Na). We calculate the sensitivity to spin-independent WIMP-nucleon elastic scattering. An exposure of one tonne $\times$ year could reach a cross-section of about 3$\times$$10^{-11}$ pb.

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.

Though there are now many hundreds of confirmed exoplanets known, the binarity of exoplanet host stars is not well understood. This is particularly true of host stars which harbor a giant planet in a highly eccentric orbit since these are more likely to have had a dramatic dynamical history which transferred angular momentum to the planet. Here we present observations of four exoplanet host stars which utilize the excellent resolving power of the Differential Speckle Survey Instrument (DSSI) on the Gemini North telescope. Two of the stars are giants and two are dwarfs. Each star is host to a giant planet with an orbital eccentricity > 0.5 and whose radial velocity data contain a trend in the residuals to the Keplerian orbit fit. These observations rule out stellar companions 4-8 magnitudes fainter than the host star at passbands of 692nm and 880nm. The resolution and field-of-view of the instrument result in exclusion radii of 0.05-1.4 arcsecs which excludes stellar companions within several AU of the host star in most cases. We further provide new radial velocities for the HD 4203 system which confirm that the linear trend previously observed in the residuals is due to an additional planet. These results place dynamical constraints on the source of the planet's eccentricities, constraints on additional planetary companions, and informs the known distribution of multiplicity amongst exoplanet host stars.

(Abridged for arXiv)We investigate the properties of halos, galaxies and blackholes to z=0 in the high resolution hydrodynamical simulation MassiveBlack-II (MBII) which evolves a LCDM cosmology in a comoving volume Vbox=100(Mpc/h)^3. MBII is the highest resolution simulation of this size which includes a self-consistent model for star formation, black hole accretion and associated feedback. We provide a simulation browser web application which enables interactive search and tagging of halos, subhalos and their properties and publicly release our galaxy catalogs. Our analysis of the halo mass function (MF) in MBII reveals that baryons have strong effects, with changes in the halo abundance of 20-35% below the knee of the MF (Mhalo < 10^13.2 Msun/h at z=0) when compared to fits based on dark matter only simulations. We provide a fitting function for the halo MF out to redshift z=11 and discuss how the onset of non-universality in the MF limits the accuracy of our fit. We study the halo occupation distribution and clustering of galaxies, in particular the evolution and scale dependence of stochasticity and bias finding reasonable agreement with observational data. The shape of the cosmic spectral energy distribution predicted by MBII is consistent with observations, but lower in amplitude. The Galaxy Stellar Mass Function (GSMF) function is broadly consistent with observations at z>=2. At z<2, the population of passive low mass (for M*<10^9 Msun) galaxies in MBII makes the GSMF too steep compared to observations whereas at the high mass end (M*>10^11 Msun) galaxies hosting bright AGNs make significant contributions to the GSMF. The quasar bolometric luminosity function is also largely consistent with observations. We note however that more efficient AGN feedback (beyond simple thermal coupling used here) is likely necessary for the largest, rarest objects/clusters at low redshifts.

We use zoom-in techniques to re-simulate three high-redshift (z > 5.5) halos which host 10^9 solar mass blackholes from the ~ Gpc volume, MassiveBlack cosmological hydrodynamic simulation. We examine a number of factors potentially affecting supermassive blackhole growth at high redshift in cosmological simulations. These include numerical resolution, feedback prescriptions and formulation of smoothed particle hydrodynamics. We find that varying the size of the region over which feedback energy is deposited directly, either for fixed number of neighbours or fixed volume makes very little difference to the accretion history of blackholes. Changing mass resolution by factors of up to 64 also does not change the blackhole growth history significantly. We find that switching from the density-entropy formulation to the pressure-entropy formulation of smoothed particle hydrodynamics slightly increases the accretion rate onto blackholes. In general numerical details appear to have small effects on the main fueling mechanism for blackholes at these high redshifts. We examine the fashion by which this occurs, finding that the insensitivity to simulation technique seems to be a hallmark of the cold flow feeding picture of these high-z supermassive blackholes. We show that the gas that participates in critical accretion phases, in these massive objects at z > 6~7 is in all cases colder, denser, and forms more coherent streams than the average gas in the halo. This is also mostly the case when the blackhole accretion is feedback regulated (z < 6), however the distinction is less prominent. For our resimulated halos, cold flows appear to be a viable mechanism for forming the most massive blackholes in the early universe, occurring naturally in LambdaCDM models of structure formation. Not requiring fine tuning of numerical parameters, they seem to be physically inevitable in these objects.

By combining optical and near-IR observations from the Hubble Space Telescope with NIR photometry from the Spitzer Space Telescope it is possible to measure the rest-frame UV-optical colours of galaxies at z=4-8. The UV-optical spectral energy distribution of star formation dominated galaxies is the result of several different factors. These include the joint distribution of stellar masses, ages, and metallicities, and the subsequent reprocessing by dust and gas in the ISM. Using a large cosmological hydrodynamical simulation we investigate the predicted spectral energy distributions of galaxies at high-redshift with a particular emphasis on assessing the potential contribution of nebular emission. We find that the average pure stellar UV-optical colour correlates with both luminosity and redshift such that galaxies at lower-redshift and higher-luminosity are typically redder. Assuming the escape fraction of ionising photons is close to zero, the effect of nebular emission is to redden the UV-optical 1500-V_w colour by, on average, 0.4 mag at z=8 declining to 0.25 mag at z=4. Young and low-metallicity stellar populations, which typically have bluer pure stellar UV-optical colours, produce larger ionising luminosities and are thus more strongly affected by the reddening effects of nebular emission. This causes the distribution of 1500-V_w colours to narrow and the trends with luminosity and redshift to weaken. The strong effect of nebular emission leaves observed-frame colours critically sensitive to the source redshift. For example, increasing the redshift by 0.1 can result in observed frame colours changing by up to ~0.6. These predictions reinforce the need to include nebular emission when modelling the spectral energy distributions of galaxies at high-redshift and also highlight the difficultly in interpreting the observed colours of individual galaxies without precise redshifts.

We present the direct imaging detection of a faint tertiary companion to the single-lined spectroscopic binary HD 8375 AB. Initially noticed as an 53 m/s/yr Doppler acceleration by Bowler et al. 2010, we have obtained high-contrast adaptive optics observations at Keck using NIRC2 that spatially resolve HD 8375 C from its host(s). Astrometric measurements demonstrate that the companion shares a common proper-motion. We detect orbital motion in a clockwise direction. Multiband relative photometry measurements are consistent with a spectral-type of M1V. Our combined Doppler and imaging observations place a lower-limit of m>0.297Msun on its dynamical mass. We also provide a refined orbit for the inner pair using recent RV measurements obtained with HIRES. HD 8375 is one of many triple-star systems that are apparently missing in the solar neighborhood.

The observed UV continuum slope of star forming galaxies is strongly affected by the presence of dust. Its observation is then a potentially valuable diagnostic of dust attenuation, particularly at high-redshift where other diagnostics are currently inaccesible. Interpreting the observed UV continuum slope in the context of dust attenuation is often achieved assuming the empirically calibrated Meurer et al. (1999) relation. Implicit in this relation is the assumption of an intrinsic UV continuum slope ($\beta=-2.23$). However, results from numerical simulations suggest that the intrinsic UV continuum slopes of high-redshift star forming galaxies are bluer than this, and moreover vary with redshift. Using values of the intrinsic slope predicted by numerical models of galaxy formation combined with a Calzetti et al. (2000) reddening law we infer UV attenuations ($A_{1500}$) $0.35-0.5\,{\rm mag}$ ($A_{V}$: $0.14-0.2\,{\rm mag}$ assuming Calzetti et al. 2000) greater than simply assuming the Meurer relation. This has significant implications for the inferred amount of dust attenuation at very-high ($z\approx 7$) redshift given current observational constraints on $\beta$, combined with the Meurer relation, suggest dust attenuation to be virtually zero in all but the most luminous systems.

We investigate the evolution of the galaxy stellar mass function at high-redshift ($z\ge 5$) using a pair of large cosmological hydrodynamical simulations: {\em MassiveBlack} and {\em MassiveBlack-II}. By combining these simulations we can study the properties of galaxies with stellar masses greater than $10^{8}\,{\rm M_{\odot}}\,h^{-1}$ and (co-moving) number densities of $\log_{10}(\phi\, [{\rm Mpc^{-3}\,dex^{-1}}\,h^{3}])>-8$. Observational determinations of the galaxy stellar mass function at very-high redshift typically assume a relation between the observed UV luminosity and stellar mass-to-light ratio which is applied to high-redshift samples in order to estimate stellar masses. This relation can also be measured from the simulations. We do this, finding two significant differences with the usual observational assumption: it evolves strongly with redshift and has a different shape. Using this relation to make a consistent comparison between galaxy stellar mass functions we find that at $z=6$ and above the simulation predictions are in good agreement with observed data over the whole mass range. Without using the correct UV luminosity and stellar mass-to-light ratio, the discrepancy would be up to two orders of magnitude for large galaxies $>10^{10}\,{\rm M_{\odot}}\,h^{-1}$. At $z=5$, however the stellar mass function for low mass $<10^{9}\,{\rm M_{\odot}}\,h^{-1}$ galaxies is overpredicted by factors of a few, consistent with the behaviour of the UV luminosity function, and perhaps a sign that feedback in the simulation is not efficient enough for these galaxies.

We use radiative transfer to study the growth of ionized regions around the brightest, z=8 quasars in a large cosmological hydrodynamic simulation that includes black hole growth and feedback (the MassiveBlack simulation). We find that in the presence of the quasar s the comoving HII bubble radii reach 10 Mpc/h after 20 My while with the stellar component alone the HII bubbles are smaller by at least an order of magnitude. Our calculations show that several features are not captured within an analytical growth model of Stromgren spheres. The X-ray photons from hard quasar spectra drive a smooth transition from fully neutral to partially neutral in the ionization front. However the transition from partially neutral to fully ionized is significantly more complex. We measure the distance to the edge of bubbles as a function of angle and use the standard deviation of these distances as a diagnostic of the isotropy of ionized regions. We find that the overlapping of nearby ionized regions from clustered halos not only increases the anisotropy, but also is the main mechanism which allows the outer radius to grow. We therefore predict that quasar ionized bubbles at this early stage in the reionization process should be both significantly larger and more irregularly shaped than bubbles around star forming galaxies. Before the star formation rate increases and the Universe fully reionizes, quasar bubbles will form the most striking and recognizable features in 21cm maps.

We investigate the formation and properties of galaxies hosting z~6 quasars, in the gigaparsec scale cosmological hydrodynamical simulation: MassiveBlack, which includes a self-consistent model for star formation, black hole accretion and associated feedback. We show that the MassiveBlack reproduces current estimates of the galaxy stellar mass function z=5, 6. We find that quasar hosts in the simulation are compact gas rich systems with high star formations rates of SFR ~ 100-1000 Msun/yr consistent with observed properties of Sloan quasar hosts in the redshift range 5.5 < z < 6.5. We show that the star-forming gas in these galaxies predominantly originates from high density cold streams which are able to penetrate the halo and grow the galaxy at the center. MassiveBlack predicts a deviation from the local Mbh-sigma and Mbh-Mstar relation implying that black holes are relatively more massive for a given stellar host at these redshifts.

The increasing size of cosmological simulations has led to the need for new visualization techniques. We focus on Smoothed Particle Hydrodynamical (SPH) simulations run with the GADGET code and describe methods for visually accessing the entire simulation at full resolution. The simulation snapshots are rastered and processed on supercomputers into images that are ready to be accessed through a web interface (GigaPan). This allows any scientist with a web-browser to interactively explore simulation datasets in both in spatial and temporal dimensions, datasets which in their native format can be hundreds of terabytes in size or more. We present two examples, the first a static terapixel image of the MassiveBlack simulation, a P-GADGET SPH simulation with 65 billion particles, and the second an interactively zoomable animation of a different simulation with more than one thousand frames, each a gigapixel in size. Both are available for public access through the GigaPan web interface. We also make our imaging software publicly available.

Observations of the most distant bright quasars imply that billion solar mass supermassive black holes (SMBH) have to be assembled within the first eight hundred million years. Under our standard galaxy formation scenario such fast growth implies large gas densities providing sustained accretion at critical or supercritical rates onto an initial black hole seed. It has been a long standing question whether and how such high black hole accretion rates can be achieved and sustained at the centers of early galaxies. Here we use our new cosmological hydrodynamic simulation (MassiveBlack) covering a volume (0.75 \Gpc)^3 appropriate for studying the rare first quasars to show that steady high density cold gas flows responsible for assembling the first galaxies produce the high gas densities that lead to sustained critical accretion rates and hence rapid growth commensurate with the existence of ~10^9 solar mass black holes as early as z~7. We find that under these conditions quasar feedback is not effective at stopping the cold gas from penetrating the central regions and hence cannot quench the accretion until the host galaxy reaches M_halo > 10^{12} solar masses. This cold-flow driven scenario for the formation of quasars implies that they should be ubiquitous in galaxies in the early universe and that major (proto)galaxy mergers are not a requirement for efficient fuel supply and growth, particularly for the earliest SMBHs.

We present a database of well determined orbital parameters of exoplanets. This database comprises spectroscopic orbital elements measured for 427 planets orbiting 363 stars from radial velocity and transit measurements as reported in the literature. We have also compiled fundamental transit parameters, stellar parameters, and the method used for the planets discovery. This Exoplanet Orbit Database includes all planets with robust, well measured orbital parameters reported in peer-reviewed articles. The database is available in a searchable, filterable, and sortable form on the Web at http://exoplanets.org through the Exoplanets Data Explorer Table, and the data can be plotted and explored through the Exoplanets Data Explorer Plotter. We use the Data Explorer to generate publication-ready plots giving three examples of the signatures of exoplanet migration and dynamical evolution: We illustrate the character of the apparent correlation between mass and period in exoplanet orbits, the selection different biases between radial velocity and transit surveys, and that the multiplanet systems show a distinct semi-major axis distribution from apparently singleton systems.

We have analysed four ASCA observations (1994--1995, 1996--1997) and three XMM-Newton observations (2005) of this source, in all of which the source is in high/soft state. We modeled the continuum spectra with relativistic disk model kerrbb, estimated the spin of the central black hole, and constrained the spectral hardening factor f_col and the distance. If kerrbb model applies, for normally used value of f_col, the distance cannot be very small, and f_col changes with observations.

Understanding the properties of the hot corona is important for studying the accretion disks in black hole X-ray binary systems. Using the Monte-Carlo technique to simulate the inverse Compton scattering between photons emitted from the cold disk and electrons in the hot corona, we have produced two table models in the $XSPEC$ format for the spherical corona case and the disk-like (slab) corona case. All parameters in our table models are physical properties of the system and can be derived from data fitting directly. Applying the models to broad-band spectra of the black hole candidate XTE J2012+381 observed with BeppoSAX, we estimated the size of the corona and the inner radius of the disk. The size of the corona in this system is several tens of gravitational radius, and the substantial increase of the inner disk radius during the transit from hard-state to soft-state is not found.

Cyg X-1 was observed in the high state at the conjunction orbital phase (0) with HETG/Chandra. Strong and asymmetric absorption lines of highly ionized species were detected, such as Fe XXV, Fe XXIV, Fe XXIII, Si XIV, S XVI, Ne X, and etc. In the high state the profile of the absorption lines are composed of an extended red wing and a less extended blue wing. The red wings of higher ionized species are more extended than that of lower ionized species. The detection of these lines provides a way to probe the properties of the flow around the companion and the black hole in Cyg X-1 during the high state. A broad emission feature around 6.5 keV was significantly detected from the both spectra of HETG/Chandra and PCA/RXTE. This feature appears to be symmetric and can be fitted with a Gaussian function rather than the Laor disk line model of fluorescent Fe K$_\alpha$ line from an accretion disk. The implications of these results on the structure of the accretion flow of Cyg X-1 in the high state are discussed.

With excellent angular resolution, good energy resolution and broad energy band, the Chandra ACIS is the best instrument for studying the X-ray halos around some galactic X-ray point sources caused by the dust scattering of X-rays in the interstellar medium. However, the direct images of bright sources obtained with ACIS usually suffer from severe pile-up. Making use of the fact that an isotropic image could be reconstructed from its projection into any direction, we can reconstruct the images of the X-ray halos from the data obtained with the HETGS and/or in CC mode. These data have no or less serious pile-up and enable us to take full advantage of the excellent angular resolution of Chandra. With the reconstructed high resolution images, we can probe the X-ray halos as close as 1$''$ to their associated point sources. Applying this method to Cygnus X-1 observed with Chandra HETGS in CC mode, we derived an energy dependent radial halo flux distribution and concluded that, in a circular region (2$'$ in radius) centered at the point source: (1) relative to the total intensity, the fractional halo intensity (FHI) is about 15% at $\sim$1 keV and drops to about 5% at $\sim$6 keV; (2) about 50% of the halo photons are within the region of a radius less than 40$''$; and (3) the spectrum of the point source is slightly distorted by the halo contamination.

For an accretion disk around a black hole, the strong relativistic effects affect every aspect of the radiation from the disk, including its spectrum, light-curve, and image. This work investigates in detail how the images of a thin disk around a black hole will be distorted, and what the observer will see from different viewing angles and in different energy bands.

We present results from observations of persistent black hole candidates with the High Energy Transmission Gratings aboard the Chandra X-ray Observatory. The sources include LMC X-1, LMC X-3, GRS 1758-258, and Cyg X-1. Along with the published results on 1E1740.7-2942, we have completed a high-resolution spectroscopic survey of such systems. The observed X-ray spectra of LMC X-1 and LMC X-3 show no prominent discrete features, while absorption edges (Mg K and Si K) are detected in the spectrum of GRS 1758-258. The edges are likely to be of interstellar origin. In most cases, the X-ray continuum can be described well by models that are often adopted in low-resolution studies of black hole candidates: a multi-temperature disk spectrum plus a Comptonization component. However, the relative contribution of the two components varies greatly among different sources. For instance, only the disk component is present for LMC X-1 and GRS 1758-258, while the Comptonized component is required for other sources. We discuss general issues related to obtaining disk parameters from modeling X-ray continuum.

We observed Cyg X-1 with {\em RXTE} contiguously over its 5.6-day binary orbit. The source was found to be in the hard state throughout the observation. Many intensity dips were detected in the X-ray light curves. We found that the dips fell into two distinct categories based on their spectral properties. One type exhibits strong energy-dependent attenuation of X-ray emission at low energies during a dip, which is characteristic of photoelectric absorption, but the other type shows nearly energy-independent attenuation. While the first type of dips are likely caused by density enhancement in an inhomogeneous wind of the companion star, as previous studies have indicated, the second type might be due to partial obscuration of an extended X-ray emitting region by optically thick ``clumps'' in the accretion flow. It is also possible that the latter are caused by a momentary decrease in the X-ray luminosity of the source, due, for instance, to a decrease in the mass accretion rate, or by Thomson scattering in highly ionized ``clumps''. We discuss the implications of these scenarios.

Recent monitoring of Cyg X-1 with {\em RXTE} revealed a period of intense flaring, which started in October of 2000 and lasted until March of 2001. The source exhibited some quite unusual behaviors during this period. The soft X-ray flux of the source went up and down three times on a timescale of about one month, as discovered by the ASM aboard RXTE, before finally returning to the normal level (of the hard state). The observed spectral and temporal X-ray properties of Cyg X-1 are mostly intermediate between the canonical hard and soft states. This is known previously for strong X-ray flares, however, we show that the source did enter a period that resembles, in many ways, a sustained soft state during the last of the three flares. We make detailed comparisons between this flare and the 1996 state transition, in terms of the observed X-ray properties, such as flux--hardness correlation, X-ray spectrum, and power density spectrum. We point out the similarities and differences, and discuss possible implications of the results on our understanding of the phenomena of flares and state transitions associated with Cyg X-1.

GX 339-4 was regularly monitored with RXTE during a period (in 1999) when its X-ray flux decreased significantly (from 4.2$\times 10^{-10}$ erg cm$^{-2} s^{-1}$ to 7.6$\times 10^{-12}$ erg cm$^{-2}$s$^{-1}$ in the 3--20 keV band), as the source settled into the ``off state''. Our spectral analysis revealed the presence of a prominent iron K$_{\alpha}$ line in the observed spectrum of the source for all observations. The line shows an interesting evolution: it is centered at $\sim$6.4 keV when the measured flux is above 5$\times 10^{-11}$ erg cm$^{-2} s^{-1}$, but is shifted to $\sim$6.7 keV at lower fluxes. The equivalent width of the line appears to increase significantly toward lower fluxes, although it is likely to be sensitive to calibration uncertainties. While the fluorescent emission of neutral or mildly ionized iron atoms in the accretion disk can perhaps account for the 6.4 keV line, as is often invoked for black hole candidates, it seems difficult to understand the 6.7 keV line with this mechanism, because the disk should be less ionized at lower fluxes (unless its density changes drastically). On the other hand, the 6.7 keV line might be due to recombination cascade of hydrogen or helium like iron ions in an optically thin, highly ionized plasma. We discuss the results in the context of proposed accretion models.

Temporal and spectral properties of X-ray rapid variability of Cyg X-1 are studied by an approach of correlation analysis in the time domain on different time scales. The correlation coefficients between the total intensity in 2-60 keV and the hardness ratio of 13-60 keV to 2-6 keV band on the time scale of about 1 ms are always negative in all states. For soft states, the correlation coefficients are positive on all the time scales from about 0.01 s to 100 s, which is significantly different with that for transition and low states. Temporal structures in high energy band are narrower than that in low energy band in quite a few cases. The delay of high energy photons relative to low energy ones in the X-ray variations has also been revealed by the correlation analysis. The implication of observed temporal and spectral characteristics to the production region and mechanism of Cyg X-1 X-ray variations is discussed.

X-ray shots of Cyg X-1 in different energy bands and spectral states have been studied with PCA/RXTE observations. The detailed shot structure is obtained by superposing many shots with one millisecond time bin through aligning their peaks with an improved algorithm. In general, the shots are composed of a slow rise and fast decay. The shot structures in the different states are different. The duration of shot in the high state is shorter than that in the low and transition states. The shot profile in the high energy band is more asymmetric and narrower than that in the low energy band. The average hardness of shot is lower than that of steady emission in the transition and low states but higher than that in the high state. The time lags between the shots in higher and lower energy bands have been found in the different states. In transition states, the time lag is the largest among the different states of Cyg X-1, and it is the smallest in the low state. The implications of the observed shot features for shot models are discussed.