M51 and NGC 5195 is an interacting system that can be explored in great details with ground-based telescopes. The H II regions in M51 were observed using the 2.16 m telescope of the National Astronomical Observatories of the Chinese Academy of Sciences and the 6.5 m Multiple Mirror Telescope with spatial resolution of less than $\sim100$ pc. We obtain a total of 113 spectra across the galaxy and combine the literature data of Croxall et al. to derive a series of physical properties, including the gas-phase extinction, stellar population age, star formation rate (SFR) surface density, and oxygen abundance. The spatial distributions and radial profiles of these properties are investigated in order to study the characteristics of M51 and the clues to the formation and evolution of this galaxy. M51 presents a mild radial extinction gradient. The lower gas-phase extinction in the north spiral arms compared to the south arms are possibly caused by the past encounters with the companion galaxy of NGC 5195. A number of H II regions have the stellar age between 50 and 500 Myr, consistent with the recent interaction history by simulations in the literatures. The SFR surface density presents a mild radial gradient, which is ubiquitous in spiral galaxies. There is a negative metallicity gradient of $-0.08$ dex $R_{e}^{-1}$ in the disk region, which is also commonly found in many spiral galaxies. It is supported by the "inside-out" scenario of galaxy formation. We find a positive abundance gradient of 0.26 dex $R_{e}^{-1}$ in the inner region. There are possible reasons causing the positive gradient, including the freezing of the chemical enrichment due to the star-forming quenching in the bulge and the gas infall and dilution due to the pseudobulge growth and/or galactic interaction.

Cold dark matter haloes consist of a relatively smooth dark matter component as well as a system of bound subhaloes. It is the prevailing practice to include all halo mass, including mass in subhaloes, in studies of halo density profiles. However, often in observational studies satellites are treated as having their own distinct dark matter density profiles in addition to the profile of the host. This difference makes comparisons between theoretical and observed results difficult. In this work we investigate density profiles of the smooth components of host haloes by excluding mass contained within subhaloes. We find that the density profiles of the smooth halo component (without subhaloes) differs substantially from the conventional halo density profile. Smooth profiles decline more rapidly at large radii and are not well characterised by the standard NFW profile. We also find that concentrations derived from smooth density profiles exhibit less scatter at fixed mass and a weaker mass dependence than standard concentrations. Both smooth and standard halo profiles can be described by a generalised Einasto profile, an Einasto profile with a modified central slope, with smaller residuals than either an NFW or Einasto profile. These results hold for both Milky Way-mass and cluster-mass haloes. This new characterisation of smooth halo profiles can be useful for many analyses, such as lensing and dark matter annihilation, in which the smooth and clumpy components of a halo should be accounted for separately.

We explore the impact of elastic, anisotropic, velocity-dependent dark matter (DM) self-interactions on the host halo and subhalos of Milky Way (MW)--mass systems. We consider a generic self-interacting dark matter (SIDM) model parameterized by the masses of a light mediator and the DM particle. The ratio of these masses, $w$, sets the velocity scale above which momentum transfer due to DM self-interactions becomes inefficient. We perform high-resolution zoom-in simulations of an MW-mass halo for values of $w$ that span scenarios in which self-interactions either between the host and its subhalos or only within subhalos efficiently transfer momentum, and we study the effects of self-interactions on the host halo and on the abundance, radial distribution, orbital dynamics, and density profiles of subhalos in each case. The abundance and properties of surviving subhalos are consistent with being determined primarily by subhalo--host halo interactions. In particular, subhalos on radial orbits in models with larger values of the cross section at the host halo velocity scale are more susceptible to tidal disruption owing to mass loss from ram pressure stripping caused by self-interactions with the host. This mechanism suppresses the abundance of surviving subhalos relative to collisionless DM simulations, with stronger suppression for larger values of $w$. Thus, probes of subhalo abundance around MW-mass hosts can be used to place upper limits on the self-interaction cross section at velocity scales of $\sim 200\ \rm{km\ s}^{-1}$, and combining these measurements with the orbital properties and internal dynamics of subhalos may break degeneracies among velocity-dependent SIDM models.

We report the discovery of two ultra-faint stellar systems found in early data from the DECam Local Volume Exploration survey (DELVE). The first system, Centaurus I (DELVE J1238-4054), is identified as a resolved overdensity of old and metal-poor stars with a heliocentric distance of ${\rm D}_{\odot} = 116.3_{-0.6}^{+0.6}$ kpc, a half-light radius of $r_h = 2.3_{-0.3}^{+0.4}$ arcmin, an age of $\tau > 12.85$ Gyr, a metallicity of $Z = 0.0002_{-0.0002}^{+0.0001}$, and an absolute magnitude of $M_V = -5.55_{-0.11}^{+0.11}$ mag. This characterization is consistent with the population of ultra-faint satellites, and confirmation of this system would make Centaurus I one of the brightest recently discovered ultra-faint dwarf galaxies. Centaurus I is detected in Gaia DR2 with a clear and distinct proper motion signal, confirming that it is a real association of stars distinct from the Milky Way foreground; this is further supported by the clustering of blue horizontal branch stars near the centroid of the system. The second system, DELVE 1 (DELVE J1630-0058), is identified as a resolved overdensity of stars with a heliocentric distance of ${\rm D}_{\odot} = 19.0_{-0.6}^{+0.5} kpc$, a half-light radius of $r_h = 0.97_{-0.17}^{+0.24}$ arcmin, an age of $\tau = 12.5_{-0.7}^{+1.0}$ Gyr, a metallicity of $Z = 0.0005_{-0.0001}^{+0.0002}$, and an absolute magnitude of $M_V = -0.2_{-0.6}^{+0.8}$ mag, consistent with the known population of faint halo star clusters. Given the low number of probable member stars at magnitudes accessible with Gaia DR2, a proper motion signal for DELVE 1 is only marginally detected. We compare the spatial position and proper motion of both Centaurus I and DELVE 1 with simulations of the accreted satellite population of the Large Magellanic Cloud (LMC) and find that neither is likely to be associated with the LMC.

The concentration parameter is a key characteristic of a dark matter halo that conveniently connects the halo's present-day structure with its assembly history. Using 'Dark Sky', a suite of cosmological $N$-body simulations, we investigate how halo concentration evolves with time and emerges from the mass assembly history. We also explore the origin of the scatter in the relation between concentration and assembly history. We show that the evolution of halo concentration has two primary modes: (1) smooth increase due to pseudo-evolution; and (2) intense responses to physical merger events. Merger events induce lasting and substantial changes in halo structures, and we observe a universal response in the concentration parameter. We argue that merger events are a major contributor to the uncertainty in halo concentration at fixed halo mass and formation time. In fact, even haloes that are typically classified as having quiescent formation histories experience multiple minor mergers. These minor mergers drive small deviations from pseudo-evolution, which cause fluctuations in the concentration parameters and result in effectively irreducible scatter in the relation between concentration and assembly history. Hence, caution should be taken when using present-day halo concentration parameter as a proxy for the halo assembly history, especially if the recent merger history is unknown.

The 21~cm intensity mapping experiments promise to measure the baryon acoustic oscillations (BAO) through the mapping of H{\small I} gas at the post-reionization epoch. Its success depends on the assumption that H{\small I} gas traces the matter density distribution {\it linearly} on large scales. However, both nonlinear halo clustering and nonlinear effects modulating H{\small I} gas in halos may spoil this. We employ three approaches to generate the mock H{\small I} density from a large-scale N-body simulation at low redshifts, and demonstrate that the assumption of H{\small I} linearity is valid at the scale corresponding to the first BAO peak, but breaks down at $k \gtrsim 0.1\,h\, {\rm Mpc}^{-1}$. The nonlinear effects of halo clustering and H{\small I} content modulation counteract each other at small scales, and their competition results in a model-dependent "sweet-spot" redshift near $z$=1 where the H{\small I} bias is scale-independent down to small scales. We also find that the linear H{\small I} bias scales approximately linearly with redshift for $z\le 3$.

The population of Milky Way (MW) satellites contains the faintest known galaxies and thus provides essential insight into galaxy formation and dark matter microphysics. Here we combine a model of the galaxy--halo connection with newly derived observational selection functions based on searches for satellites in photometric surveys over nearly the entire high Galactic latitude sky. In particular, we use cosmological zoom-in simulations of MW-like halos that include realistic Large Magellanic Cloud (LMC) analogs to fit the position-dependent MW satellite luminosity function. We report decisive evidence for the statistical impact of the LMC on the MW satellite population due to an estimated $6\pm 2$ observed LMC-associated satellites, consistent with the number of LMC satellites inferred from Gaia proper-motion measurements, confirming the predictions of cold dark matter models for the existence of satellites within satellite halos. Moreover, we infer that the LMC fell into the MW within the last $2\ \rm{Gyr}$ at high confidence. Based on our detailed full-sky modeling, we find that the faintest observed satellites inhabit halos with peak virial masses below $3.2\times 10^{8}\ M_{\rm{\odot}}$ at $95\%$ confidence, and we place the first robust constraints on the fraction of halos that host galaxies in this regime. We predict that the faintest detectable satellites occupy halos with peak virial masses above $10^{6}\ M_{\rm{\odot}}$, highlighting the potential for powerful galaxy formation and dark matter constraints from future dwarf galaxy searches.

Conventional wisdom was that thermal relics from the epoch of reionization (EOR) would vanish swiftly and hence the usual intergalactic medium (IGM) temperature-density relation would be recovered rapidly. Thus the Ly$\alpha$ forest is one of the primary cosmological probes at the post-reionization epoch. Recently, however, it was shown that the imprint of cosmic reionization can survive to lower redshifts ($z \sim 2$) than previously thought. Given the high sensitivities of upcoming Ly$\alpha$ forest surveys, this effect will be a novel broadband systematic that must be tackled for cosmological application. From the astrophysical point of view, however, the imprint of inhomogeneous reionization can shed light on the EOR and cosmic dawn. We utilize a hybrid method --- which includes two different simulation codes capable of handling the huge dynamical range --- to show the impact of patchy reionization on the Ly$\alpha$ forest and its dependence on different astrophysical scenarios. We found small, but statistically significant, deviations in the 1D Ly$\alpha$ power spectrum that range from a tenth of per cent at $z = 2$ to a few per cent at $z = 4$. The deviations in the 3D Ly$\alpha$ power spectrum are considerably large and range from a few per cent at $z = 2$ up to tens of per cent at $z = 4$. By exploiting different $k$-dependence of power spectrum among various astrophysical scenarios, the effect of patchy reionization on the Ly$\alpha$ forest power spectrum can open a new window into the cosmic reionization and possibly even the cosmic dawn.

The bubble size distribution of ionized hydrogen regions, which could be derived from the tomographic imaging data of the redshifted 21~cm signal from the epoch of reionization, probes the information about the morphology of \HII\ bubbles during the reionization. However, 21~cm imaging is observationally very challenging even for the upcoming large radio interferometers. Given that these interferometers promise to measure the 21~cm power spectrum accurately, we propose a new method, which is based on the artificial neural networks (ANN), to reconstruct the \HII\ bubble size distribution from the 21~cm power spectrum. We demonstrate that the reconstruction from the 21~cm power spectrum can be almost as accurate as directly measured from the imaging data with the fractional error $\lesssim 10\%$, even with thermal noise at the sensitivity level of the Square Kilometre Array. Nevertheless, systematic errors might arise from approximations made in reionization simulations used for training the ANN. This paper, as the first in a series, exemplifies the possibility of recovering from the 21~cm power spectrum with ANN additional statistics of cosmic reionization that could not otherwise be inferred from the power spectrum analysis directly in the conventional methods.

We present measurements of the redshift-dependent clustering of a DESI-like luminous red galaxy (LRG) sample selected from the Legacy Survey imaging dataset, and use the halo occupation distribution (HOD) framework to fit the clustering signal. The LRG sample contains 2.7 million objects over the redshift range of $0.4 < z < 0.9$ over 5655 sq. degrees. We have developed new photometric redshift (photo-$z$) estimates using the Legacy Survey DECam and WISE photometry, with $\sigma_{\mathrm{NMAD}} = 0.02$ precision for LRGs. We compute the projected correlation function using new methods that maximize signal-to-noise while incorporating redshift uncertainties. We present a novel algorithm for dividing irregular survey geometries into equal-area patches for jackknife resampling. For a 5-parameter HOD model fit using the MultiDark halo catalog, we find that there is little evolution in HOD parameters except at the highest-redshifts. The inferred large-scale structure bias is largely consistent with constant clustering amplitude over time. In an appendix, we explore limitations of MCMC fitting using stochastic likelihood estimates resulting from applying HOD methods to N-body catalogs, and present a new technique for finding best-fit parameters in this situation. Accompanying this paper we have released the PRLS (Photometric Redshifts for the Legacy Surveys) catalog of photo-$z$'s obtained by applying the methods used in this work to the full Legacy Survey Data Release 8 dataset. This catalog provides accurate photometric redshifts for objects with $z < 21$ over more than 16,000 square degrees of sky.

Cosmic reionization was driven by the imbalance between early sources and sinks of ionizing radiation, both of which were dominated by small-scale structure and are thus usually treated in cosmological reionization simulations by subgrid modelling. The recombination rate of intergalactic hydrogen is customarily boosted by a subgrid clumping factor, ${\left<n^2\right>/\left<n\right>^2}$, which corrects for unresolved fluctuations in gas density ${n}$ on scales below the grid-spacing of coarse-grained simulations. We investigate in detail the impact of this inhomogeneous subgrid clumping on reionization and its observables, as follows: (1) Previous attempts generally underestimated the clumping factor because of insufficient mass resolution. We perform a high-resolution $N$-body simulation that resolves haloes down to the pre-reionization Jeans mass to derive the time-dependent, spatially-varying local clumping factor and a fitting formula for its correlation with local overdensity. (2) We then perform a large-scale $N$-body and radiative transfer simulation that accounts for this inhomogeneous subgrid clumping by applying this clumping factor-overdensity correlation. Boosting recombination significantly slows the expansion of ionized regions, which delays completion of reionization and suppresses 21 cm power spectra on large scales in the later stages of reionization. (3) We also consider a simplified prescription in which the globally-averaged, time-evolving clumping factor from the same high-resolution $N$-body simulation is applied uniformly to all cells in the reionization simulation, instead. Observables computed with this model agree fairly well with those from the inhomogeneous clumping model, e.g. predicting 21 cm power spectra to within 20% error, suggesting it may be a useful approximation.

We demonstrate how the total luminosity in satellite galaxies is a powerful probe of dark matter halos around central galaxies. The method cross-correlates central galaxies in spectroscopic galaxy samples with fainter galaxies detected in photometric surveys. After background subtraction, the excess galaxies around the central galaxies represent faint satellite galaxies within the dark matter halo. Using abundance matching models, we show that the the total galaxy luminosity, L_sat, scales linearly with host halo mass, making L_sat an excellent proxy for M_h. L_sat is also sensitive to the formation time of the halo, as younger halos have more substructure at fixed M_h. We demonstrate that probes of galaxy large-scale environment can break this degeneracy. Although this is an indirect probe of the halo, it can yield a high-S/N measurement for galaxies expected to occupy halos at $<10^{12}$ M_sol, where other methods suffer from larger errors. In this paper we focus on observational and theoretical systematics in the L_sat method. We test the robustness of our method of finding central galaxies and our methods of estimating the number of background galaxies. We implement this method on central galaxies in SDSS data, with satellites identified in faint imaging from the DESI Legacy Imaging Surveys. We find excellent agreement between our theoretical predictions and the observational measurements. Finally, we compare our L_sat measurements to weak lensing estimates of M_h for red and blue subsamples. In the stellar mass range where the measurements overlap, we find consistent results, where red galaxies live in larger halos. However, the L_sat approach allows us to probe significantly lower mass galaxies. At these masses, the L_sat values are equivalent. This example shows the potential of L_sat as a probe of dark halos. (Abridged)

We introduce the Southern Stellar Stream Spectroscopy Survey (${S}^5$), an on-going program to map the kinematics and chemistry of stellar streams in the Southern Hemisphere. The initial focus of ${S}^5$ has been spectroscopic observations of recently identified streams within the footprint of the Dark Energy Survey (DES), with the eventual goal of surveying streams across the entire southern sky. Stellar streams are composed of material that has been tidally striped from dwarf galaxies and globular clusters and hence are excellent dynamical probes of the gravitational potential of the Milky Way, as well as providing a detailed snapshot of its accretion history. Observing with the 3.9-m Anglo-Australian Telescope's 2-degree-Field fibre positioner and AAOmega spectrograph, and combining the precise photometry of DES DR1 with the superb proper motions from $Gaia$ DR2, allows us to conduct an efficient spectroscopic survey to map these stellar streams. So far ${S}^5$ has mapped 9 DES streams and 3 streams outside of DES; the former are the first spectroscopic observations of these recently discovered streams. In addition to the stream survey, we use spare fibres to undertake a Milky Way halo survey and a low-redshift galaxy survey. This paper presents an overview of the ${S}^5$ program, describing the scientific motivation for the survey, target selection, observation strategy, data reduction and survey validation. Finally, we describe early science results on stellar streams and Milky Way halo stars drawn from the survey. Updates on ${S}^5$, including future public data release, can be found at \url{http://s5collab.github.io}.

We develop and implement two new methods for constraining the scatter in the relationship between galaxies and dark matter halos. These new techniques are sensitive to the scatter at low halo masses, making them complementary to previous constraints that are dependent on clustering amplitudes or rich galaxy groups, both of which are only sensitive to more massive halos. In both of our methods, we use a galaxy group finder to locate central galaxies in the SDSS main galaxy sample. Our first technique uses the small-scale cross-correlation of central galaxies with all lower mass galaxies. This quantity is sensitive to the satellite fraction of low-mass galaxies, which is in turn driven by the scatter between halos and galaxies. The second technique uses the kurtosis of the distribution of line-of-sight velocities between central galaxies and neighboring galaxies. This quantity is sensitive to the distribution of halo masses that contain the central galaxies at fixed stellar mass. Theoretical models are constructed using peak halo circular velocity, $V_{\rm peak}$, as our property to connect galaxies to halos. The cross-correlation technique yields a constraint of $\sigma[ M_\ast|V_{\rm peak}]=0.27\pm 0.05$ dex, corresponding to a scatter in $\log M_\ast$ at fixed $M_h$ of $\sigma[ M_\ast|M_h]=0.38\pm 0.06$ dex at $M_h=10^{11.8}$ Msun. The kurtosis technique yields $\sigma[ M_\ast|V_{\rm peak}]=0.30\pm0.03$, corresponding to $\sigma[ M_\ast|M_h]=0.34\pm 0.04$ at $M_h=10^{12.2}$ Msun. The values of $\sigma[ M_\ast|M_h]$ are significantly larger than the constraints at higher masses, in agreement with the results of hydrodynamic simulations. This increase is only partly due to the scatter between $V_{\rm peak}$ and $M_h$, and it represents an increase of nearly a factor of two relative to the values inferred from clustering and group studies at high masses.

Most models for the connection between galaxies and their haloes ignore the possibility that galaxy properties may be correlated with halo properties other than mass, a phenomenon known as galaxy assembly bias. Yet, it is known that such correlations can lead to systematic errors in the interpretation of survey data. At present, the degree to which galaxy assembly bias may be present in the real Universe, and the best strategies for constraining it remain uncertain. We study the ability of several observables to constrain galaxy assembly bias from redshift survey data using the decorated halo occupation distribution (dHOD), an empirical model of the galaxy--halo connection that incorporates assembly bias. We cover an expansive set of observables, including the projected two-point correlation function $w_{\mathrm{p}}(r_{\mathrm{p}})$, the galaxy--galaxy lensing signal $\Delta \Sigma(r_{\mathrm{p}})$, the void probability function $\mathrm{VPF}(r)$, the distributions of counts-in-cylinders $P(N_{\mathrm{CIC}})$, and counts-in-annuli $P(N_{\mathrm{CIA}})$, and the distribution of the ratio of counts in cylinders of different sizes $P(N_2/N_5)$. We find that despite the frequent use of the combination $w_{\mathrm{p}}(r_{\mathrm{p}})+\Delta \Sigma(r_{\mathrm{p}})$ in interpreting galaxy data, the count statistics, $P(N_{\mathrm{CIC}})$ and $P(N_{\mathrm{CIA}})$, are generally more efficient in constraining galaxy assembly bias when combined with $w_{\mathrm{p}}(r_{\mathrm{p}})$. Constraints based upon $w_{\mathrm{p}}(r_{\mathrm{p}})$ and $\Delta \Sigma(r_{\mathrm{p}})$ share common degeneracy directions in the parameter space, while combinations of $w_{\mathrm{p}}(r_{\mathrm{p}})$ with the count statistics are more complementary. Therefore, we strongly suggest that count statistics should be used to complement the canonical observables in future studies of the galaxy--halo connection.

This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `Outer Rim' run. It covers 440 deg^2 of sky area to a redshift of z=3 and is complete to a magnitude depth of 28 in the r-band. Each galaxy is characterized by a multitude of properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-driven empirical approaches with semi-analytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data.

Due to ionosphere absorption and the interference by natural and artificial radio emissions, ground observation of the sky at the decameter or longer is very difficult. This unexplored part of electromagnetic spectrum has the potential of great discoveries, notably in the study of cosmic dark ages and dawn, but also in heliophysics and space weather, planets, cosmic ray and neutrinos, pulsar and interstellar medium, extragalactic radio sources, and even SETI. At a forum organized by the International Space Science Institute-Beijing (ISSI-BJ), we discussed the prospect of opening up this window for astronomical observations by using a constellation of small or micro-satellites. We discussed the past experiments and the current ones such as the low frequency payload on Chang'e-4 mission lander, relay satellite and the Longjiang satellite, and also the future DSL mission, which is a linear array on lunar orbit which can make synthesized map of the whole sky as well as measure the global spectrum. We also discuss the synergy with other experiments, including ground global experiments such as EDGES, SARAS, SCI-HI and High-z, PRIZM/Albatros, ground imaging facillities such as LOFAR and MWA, and space experiments such as SUNRISE, DARE/DAPPER and PRATUSH. We also discussed some technical aspects of the DSL concept.

On small scales there have been a number of claims of discrepancies between the standard Cold Dark Matter (CDM) model and observations. The 'missing satellites problem' infamously describes the overabundance of subhalos from CDM simulations compared to the number of satellites observed in the Milky Way. A variety of solutions to this discrepancy have been proposed; however, the impact of the specific properties of the Milky Way halo relative to the typical halo of its mass have yet to be explored. Motivated by recent studies that identified ways in which the Milky Way is atypical (e.g., Licquia et al. 2015), we investigate how the properties of dark matter halos with mass comparable to our Galaxy's --- including concentration, spin, shape, and scale factor of the last major merger --- correlate with the subhalo abundance. Using zoom-in simulations of Milky Way-like halos, we build two models of subhalo abundance as functions of host halo properties and conclude that the Milky Way should be expected to have 22%-44% fewer subhalos with low maximum rotation velocities ($V_{\rm max}^{\rm sat} \sim 10$kms$^{-1}$) at the 95% confidence level and up to 72% fewer than average subhalos with high rotation velocities ($V_{\rm max}^{\rm sat} \gtrsim 30$kms$^{-1}$, comparable to the Magellanic Clouds) than would be expected for a typical halo of the Milky Way's mass. Concentration is the most informative single parameter for predicting subhalo abundance. Our results imply that models tuned to explain the missing satellites problem assuming typical subhalo abundances for our Galaxy will be over-correcting.

Hemispherical power asymmetry has emerged as a new challenge to cosmology in early universe. While the cosmic microwave background (CMB) measurements indicated the asymmetry amplitude $A \simeq 0.07$ at the CMB scale $k_{\rm CMB}\simeq 0.0045\,{\rm Mpc}^{-1}$, the high-redshift quasar observations found no significant deviation from statistical isotropy. This conflict can be reconciled in some scale-dependent asymmetry models. We put forward a new parameterization of scale-dependent asymmetric power spectrum, inspired by a multi-speed inflation model. The 21-cm power spectrum from the epoch of reionization can be used to constrain the scale-dependent hemispherical asymmetry. We demonstrate that an optimum, multi-frequency observation by the Square Kilometre Array (SKA) Phase 2 can impose a constraint on the amplitude of the power asymmetry anomaly at the level of $\Delta A \simeq 0.2$ at $0.056 \lesssim k_{\rm 21cm} \lesssim 0.15 \,{\rm Mpc}^{-1}$. This limit may be further improved by an order of magnitude as $\Delta A \simeq 0.01$ with a cosmic variance limited experiment such as the Omniscope.

The three-point correlation function (3PCF) of the 21cm brightness temperature from the Epoch of Reionization (EoR) probes complementary information to the commonly studied two-point correlation function (2PCF) about the morphology of ionized regions. We investigate the 21cm 2PCF and 3PCF in configuration space using semi-numerical simulations and test whether they can be described by the local quadratic bias model. We find that fits of bias model predictions for the 2PCF and 3PCF deviate from our measurements by $\sim 20\%$ at scales above the typical size of ionized regions ($\simeq 30$ Mpc) and at early times with global neutral fractions of $\langle x_{\rm HI} \rangle \gtrsim 0.7$. At later times and smaller scales these deviations increase strongly, indicating a break down of the bias model. The 2PCF and 3PCF fits of the linear bias parameter agree at the $10\%$ level for different EoR model configurations. This agreement holds, when adding redshift space distortions to the simulations. The relation between spatial fluctuations in the matter density and the 21cm signal, as predicted by the bias model, is consistent with direct measurements of this relation in simulations for large smoothing scales ($\gtrsim 30$ Mpc). From this latter test we conclude that negative amplitudes of the 21cm 3PCF result from negative bias parameters, which describe the anti-correlation between the matter over-densities and the 21cm signal during the EoR. However, a more detailed interpretation of the bias parameters may require a description of non-local contributions to the bias model.

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

We discuss how astrophysical observations with the Maunakea Spectroscopic Explorer (MSE), a high-multiplexity (about 4300 fibers), wide field-of-view (1.5 square degree), large telescope aperture (11.25 m) facility, can probe the particle nature of dark matter. MSE will conduct a suite of surveys that will provide critical input for determinations of the mass function, phase-space distribution, and internal density profiles of dark matter halos across all mass scales. N-body and hydrodynamical simulations of cold, warm, fuzzy and self-interacting dark matter suggest that non-trivial dynamics in the dark sector could have left an imprint on structure formation. Analysed within these frameworks, the extensive and unprecedented datasets produced by MSE will be used to search for deviations away from cold and collisionless dark matter model. MSE will provide an improved estimate of the local density of dark matter, critical for direct detection experiments, and will improve estimates of the J-factor for indirect searches through self-annihilation or decay into Standard Model particles. MSE will determine the impact of low mass substructures on the dynamics of Milky Way stellar streams in velocity space, and will allow for estimates of the density profiles of the dark matter halos of Milky Way dwarf galaxies using more than an order of magnitude more tracers. In the low redshift Universe, MSE will provide critical redshifts to pin down the luminosity functions of vast numbers of satellite systems, and MSE will be an essential component of future strong lensing measurements to constrain the halo mass function. Across nearly all mass scales, the improvements offered by MSE, in comparison to other facilities, are such that the relevant analyses are limited by systematics rather than statistics.

The early phases of galaxy formation constitute one of the most exciting frontiers in astrophysics. It is during this era that the first luminous sources reionize the intergalactic medium - the moment when structure formation affects every baryon in the Universe. Here we argue that we will obtain a complete picture of this era by combining observations of galaxies with direct measurements of the reionization process: the former will provide a detailed understanding of bright sources, while the latter will constrain the (substantial) faint source population. We further describe how optimizing the comparison of these two measurements requires near-infrared galaxy surveys covering large volumes and retaining redshift information and also improvements in 21-cm analysis, moving those experiments into the imaging regime.

The epoch of reionization, when photons from early galaxies ionized the intergalactic medium about a billion years after the Big Bang, is the last major phase transition in the Universe's history. Measuring the characteristics of the transition is important for understanding early galaxies and the cosmic web and for modeling dwarf galaxies in the later Universe. But such measurements require probes of the intergalactic medium itself. Here we describe how the 21-cm line of neutral hydrogen provides a powerful probe of the reionization process and therefore important constraints on both the galaxies and intergalactic absorbers at that time. While existing experiments will make precise statistical measurements over the next decade, we argue that improved 21-cm analysis techniques - allowing imaging of the neutral gas itself - as well as improved theoretical models, are crucial for testing our understanding of this important era.

The Dark Ages are the period between the last scattering of the cosmic microwave background and the appearance of the first luminous sources, spanning approximately 1100 < z < 30. The only known way to measure fluctuations in this era is through the 21-cm line of neutral hydrogen. Such observations have enormous potential for cosmology, because they span a large volume while the fluctuations remain linear even on small scales. Observations of 21-cm fluctuations during this era can therefore constrain fundamental aspects of our Universe, including inflation and any exotic physics of dark matter. While the observational challenges to these low-frequency 21-cm observations are enormous, especially from the terrestrial environment, they represent an important goal for cosmology.

The "cosmic dawn" refers to the period of the Universe's history when stars and black holes first formed and began heating and ionizing hydrogen in the intergalactic medium (IGM). Though exceedingly difficult to detect directly, the first stars and black holes can be constrained indirectly through measurements of the cosmic 21-cm background, which traces the ionization state and temperature of intergalactic hydrogen gas. In this white paper, we focus on the science case for such observations, in particular those targeting redshifts z $\gtrsim$ 10 when the IGM is expected to be mostly neutral. 21-cm observations provide a unique window into this epoch and are thus critical to advancing first star and black hole science in the next decade.

In addition to being a probe of Cosmic Dawn and Epoch of Reionization astrophysics, the 21cm line at $z>6$ is also a powerful way to constrain cosmology. Its power derives from several unique capabilities. First, the 21cm line is sensitive to energy injections into the intergalactic medium at high redshifts. It also increases the number of measurable modes compared to existing cosmological probes by orders of magnitude. Many of these modes are on smaller scales than are accessible via the CMB, and moreover have the advantage of being firmly in the linear regime (making them easy to model theoretically). Finally, the 21cm line provides access to redshifts prior to the formation of luminous objects. Together, these features of 21cm cosmology at $z>6$ provide multiple pathways toward precise cosmological constraints. These include the "marginalizing out" of astrophysical effects, the utilization of redshift space distortions, the breaking of CMB degeneracies, the identification of signatures of relative velocities between baryons and dark matter, and the discovery of unexpected signs of physics beyond the $\Lambda$CDM paradigm at high redshifts.

Cosmic dawn and the Epoch of Reionization (EoR) are among the least explored observational eras in cosmology: a time at which the first galaxies and supermassive black holes formed and reionized the cold, neutral Universe of the post-recombination era. With current instruments, only a handful of the brightest galaxies and quasars from that time are detectable as individual objects, due to their extreme distances. Fortunately, a multitude of multi-wavelength intensity mapping measurements, ranging from the redshifted 21 cm background in the radio to the unresolved X-ray background, contain a plethora of synergistic information about this elusive era. The coming decade will likely see direct detections of inhomogenous reionization with CMB and 21 cm observations, and a slew of other probes covering overlapping areas and complementary physical processes will provide crucial additional information and cross-validation. To maximize scientific discovery and return on investment, coordinated survey planning and joint data analysis should be a high priority, closely coupled to computational models and theoretical predictions.

Over the past decade, empirical constraints on the galaxy-dark matter halo connection have significantly advanced our understanding of galaxy evolution. Past techniques have focused on connections between halo properties and galaxy stellar mass and/or star formation rates. Empirical techniques in the next decade will link halo assembly histories with galaxies' circumgalactic media, supermassive black holes, morphologies, kinematics, sizes, colors, metallicities, and transient rates. Uncovering these links will resolve many critical uncertainties in galaxy formation and will enable much higher-fidelity mock catalogs essential for interpreting observations. Achieving these results will require broader and deeper spectroscopic coverage of galaxies and their circumgalactic media; survey teams will also need to meet several criteria (cross-comparisons, public access, and covariance matrices) to facilitate combining data across different surveys. Acting on these recommendations will continue enabling dramatic progress in both empirical modeling and galaxy evolution for the next decade.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

We develop a comprehensive and flexible model for the connection between satellite galaxies and dark matter subhalos in dark matter-only zoom-in simulations of Milky Way (MW)--mass host halos. We systematically identify the physical and numerical uncertainties in the galaxy--halo connection and simulations underlying our method, including (i) the influence of host halo properties; (ii) the relationship between satellite luminosities and subhalo properties, including the effects of reionization; (iii) the relationship between satellite and subhalo locations; (iv) the relationship between satellite sizes and subhalo properties, including the effects of tidal stripping; (v) satellite and subhalo disruption due to baryonic effects; and (vi) artificial subhalo disruption and orphan satellites. To illustrate our approach, we fit this model to the luminosity distribution of both classical MW satellites and those discovered in the Sloan Digital Sky Survey by performing realistic mock observations that depend on the luminosity, size, and distance of our predicted satellites, and we infer the total satellite population that will be probed by upcoming surveys. We argue that galaxy size and surface brightness modeling will play a key role in interpreting current and future observations, as the expected number of observable satellites depends sensitively on their surface brightness distribution.

We report on the analysis of 34 years of photometric observations of the pulsating helium atmosphere white dwarf GD358. The complete data set includes archival data from 1982-2006, and 1195.2 hours of new observations from 2007- 2016. From this data set, we extract 15 frequencies representing g-mode pulsation modes, adding 4 modes to the 11 modes known previously. We present evidence that these 15 modes are ell = 1 modes, 13 of which belong to a consecutive sequence in radial overtone k. We perform a detailed asteroseismic analysis using models that include parameterized, complex carbon and oxygen core composition profiles to fit the periods. Recent spectroscopic analyses place GD358 near the red edge of the DBV instability strip, at 24,000 plus or minus 500 K and a log g of 7.8 plus or minus 0.08 dex. The surface gravity translates to a mass range of 0.455 to 0.540 solar masses. Our best fit model has a temperature of 23,650 K and a mass of 0.5706 solar masses. That is slightly more massive than suggested by most the recent spectroscopy. We find a pure helium layer mass of 10^-5.50, consistent with the result of previous studies and the outward diffusion of helium over time.

Secondary halo bias, commonly known as 'assembly bias,' is the dependence of halo clustering on a halo property other than mass. This prediction of the Lambda-Cold Dark Matter cosmology is essential to modelling the galaxy distribution to high precision and interpreting clustering measurements. As the name suggests, different manifestations of secondary halo bias have been thought to originate from halo assembly histories. We show conclusively that this is incorrect for cluster-size haloes. We present an up-to-date summary of secondary halo biases of high-mass haloes due to various halo properties including concentration, spin, several proxies of assembly history, and subhalo properties. While concentration, spin, and the abundance and radial distribution of subhaloes exhibit significant secondary biases, properties that directly quantify halo assembly history do not. In fact, the entire assembly histories of haloes in pairs are nearly identical to those of isolated haloes. In general, a global correlation between two halo properties does not predict whether or not these two properties exhibit similar secondary biases. For example, assembly history and concentration (or subhalo abundance) are correlated for both paired and isolated haloes, but follow slightly different conditional distributions in these two cases. This results in a secondary halo bias due to concentration (or subhalo abundance), despite the lack of assembly bias in the strict sense for cluster-size haloes. Due to this complexity, caution must be exercised in using any one halo property as a proxy to study the secondary bias due to another property.

We identify subhalos in dark matter-only (DMO) zoom-in simulations that are likely to be disrupted due to baryonic effects by using a random forest classifier trained on two hydrodynamic simulations of Milky Way (MW)-mass host halos from the Latte suite of the Feedback in Realistic Environments (FIRE) project. We train our classifier using five properties of each disrupted and surviving subhalo: pericentric distance and scale factor at first pericentric passage after accretion, and scale factor, virial mass, and maximum circular velocity at accretion. Our five-property classifier identifies disrupted subhalos in the FIRE simulations with an $85\%$ out-of-bag classification score. We predict surviving subhalo populations in DMO simulations of the FIRE host halos, finding excellent agreement with the hydrodynamic results; in particular, our classifier outperforms DMO zoom-in simulations that include the gravitational potential of the central galactic disk in each hydrodynamic simulation, indicating that it captures both the dynamical effects of a central disk and additional baryonic physics. We also predict surviving subhalo populations for a suite of DMO zoom-in simulations of MW-mass host halos, finding that baryons impact each system consistently and that the predicted amount of subhalo disruption is larger than the host-to-host scatter among the subhalo populations. Although the small size and specific baryonic physics prescription of our training set limits the generality of our results, our work suggests that machine-learning classification algorithms trained on hydrodynamic zoom-in simulations can efficiently predict realistic subhalo populations.

Using the N-body simulations of the AEMULUS Project, we construct an emulator for the non-linear clustering of galaxies in real and redshift space. We construct our model of galaxy bias using the halo occupation framework, accounting for possible velocity bias. The model includes 15 parameters, including both cosmological and galaxy bias parameters. We demonstrate that our emulator achieves ~ 1% precision at the scales of interest, 0.1<r<10 h^{-1} Mpc, and recovers the true cosmology when tested against independent simulations. Our primary parameters of interest are related to the growth rate of structure, f, and its degenerate combination fsigma_8. Using this emulator, we show that the constraining power on these parameters monotonically increases as smaller scales are included in the analysis, all the way down to 0.1 h^{-1} Mpc. For a BOSS-like survey, the constraints on fsigma_8 from r<30 h^{-1} Mpc scales alone are more than a factor of two tighter than those from the fiducial BOSS analysis of redshift-space clustering using perturbation theory at larger scales. The combination of real- and redshift-space clustering allows us to break the degeneracy between f and sigma_8, yielding a 9% constraint on f alone for a BOSS-like analysis. The current AEMULUS simulations limit this model to surveys of massive galaxies. Future simulations will allow this framework to be extended to all galaxy target types, including emission-line galaxies.

The rapidly growing statistical precision of galaxy surveys has lead to a need for ever-more precise predictions of the observables used to constrain cosmological and galaxy formation models. The primary avenue through which such predictions will be obtained is suites of numerical simulations. These simulations must span the relevant model parameter spaces, be large enough to obtain the precision demanded by upcoming data, and be thoroughly validated in order to ensure accuracy. In this paper we present one such suite of simulations, forming the basis for the AEMULUS Project, a collaboration devoted to precision emulation of galaxy survey observables. We have run a set of 75 (1.05 h^-1 Gpc)^3 simulations with mass resolution and force softening of 3.51\times 10^10 (Omega_m / 0.3) ~ h^-1 M_sun and 20 ~ h^-1 kpc respectively in 47 different wCDM cosmologies spanning the range of parameter space allowed by the combination of recent Cosmic Microwave Background, Baryon Acoustic Oscillation and Type Ia Supernovae results. We present convergence tests of several observables including spherical overdensity halo mass functions, galaxy projected correlation functions, galaxy clustering in redshift space, and matter and halo correlation functions and power spectra. We show that these statistics are converged to 1% (2%) for halos with more than 500 (200) particles respectively and scales of r>200 ~ h^-1 kpc in real space or k ~ 3 h Mpc^-1 in harmonic space for z\le 1. We find that the dominant source of uncertainty comes from varying the particle loading of the simulations. This leads to large systematic errors for statistics using halos with fewer than 200 particles and scales smaller than k ~ 4 h^-1 Mpc. We provide the halo catalogs and snapshots detailed in this work to the community at https://AemulusProject.github.io.

Existing models for the dependence of the halo mass function on cosmological parameters will become a limiting source of systematic uncertainty for cluster cosmology in the near future. We present a halo mass function emulator and demonstrate improved accuracy relative to state-of-the-art analytic models. In this work, mass is defined using an overdensity criteria of 200 relative to the mean background density. Our emulator is constructed from the AEMULUS simulations, a suite of 40 N-body simulations with snapshots from z=3 to z=0. These simulations cover the flat wCDM parameter space allowed by recent Cosmic Microwave Background, Baryon Acoustic Oscillation and Type Ia Supernovae results, varying the parameters w, Omega_m, Omega_b, sigma_8, N_{eff}, n_s, and H_0. We validate our emulator using five realizations of seven different cosmologies, for a total of 35 test simulations. These test simulations were not used in constructing the emulator, and were run with fully independent initial conditions. We use our test simulations to characterize the modeling uncertainty of the emulator, and introduce a novel way of marginalizing over the associated systematic uncertainty. We confirm non-universality in our halo mass function emulator as a function of both cosmological parameters and redshift. Our emulator achieves better than 1% precision over much of the relevant parameter space, and we demonstrate that the systematic uncertainty in our emulator will remain a negligible source of error for cluster abundance studies through at least the LSST Year 1 data set.

By using the Hectospec 6.5 m Multiple Mirror Telescope (MMT) and the 2.16 m telescope of National Astronomical Observatories, Chinese Academy of Sciences (NAOC), we obtained 188 high signal-to-noise ratio (S/N) spectra of HII regions in the nearby galaxy M101, which are the largest spectroscopic sample of HII regions for this galaxy so far. These spectra cover a wide range of regions on M101, which enables us to analyze two dimensional distributions of its physical properties. The physical parameters are derived from emission lines or stellar continuum, including stellar population age, electron temperature, oxygen abundance and etc. The oxygen abundances are derived using two empirical methods based on O3N2 and R$_{23}$ indicators, as well as the direct Te method when OIII$\lambda4363$ is available. By applying the harmonic decomposition analysis to the velocity field, we obtained line-of-sight rotation velocity of 71 km s$^{-1}$ and a position angle of 36 degree. The stellar age profile shows an old stellar population in galaxy center and a relative young stellar population in outer regions, suggesting an old bulge and a young disk. Oxygen abundance profile exhibits a clear break at $\sim$18 kpc, with a gradient of $-$0.0364 dex kpc$^{-1}$ in the inner region and $-$0.00686 dex kpc$^{-1}$ in the outer region. Our results agree with the "inside-out" disk growth scenario of M101.

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

This paper presents a spectroscopic investigation of 11 HII regions in the nearby galaxy NGC 2403. The HII regions are observed with a long-slit spectrograph mounted on the 2.16 m telescope at XingLong station of National Astronomical Observatories of China. For each of the HII regions, spectra are extracted at different nebular radii along the slit-coverage. Oxygen abundances are empirically estimated from the strong-line indices R23, N2O2, O3N2, and N2 for each spectrophotometric unit, with both observation- and model-based calibrations adopted into the derivation. Radial profiles of these diversely estimated abundances are drawn for each nebula. In the results, the oxygen abundances separately estimated with the prescriptions on the basis of observations and models, albeit from the same spectral index, systematically deviate from each other; at the same time, the spectral indices R23 and N2O2 are distributed with flat profiles, whereas N2 and O3N2 exhibit apparent gradients with the nebular radius. Because our study naturally samples various ionization levels which inherently decline at larger radii within individual HII regions, the radial distributions indicate not only the robustness of R23 and N2O2 against ionization variations but also the sensitivity of N2 and O3N2 to the ionization parameter. The results in this paper provide observational corroboration of the theoretical prediction about the deviation in the empirical abundance diagnostics. Our future work is planned to investigate metal-poor HII regions with measurable T_e, in an attempt to recalibrate the strong-line indices and consequently disclose the cause of the discrepancies between the empirical oxygen abundances.

Determining the positions of halo centres in large-scale structure surveys is crucial for many cosmological studies. A common assumption is that halo centres correspond to the location of their brightest member galaxies. In this paper, we study the dynamics of brightest galaxies with respect to other halo members in the Sloan Digital Sky Survey DR7. Specifically, we look at the line-of-sight velocity and spatial offsets between brightest galaxies and their neighbours. We compare those to detailed mock catalogues, constructed from high-resolution, dark-matter-only $N$-body simulations, in which it is assumed that satellite galaxies trace dark matter subhaloes. This allows us to place constraints on the fraction $f_{\rm BNC}$ of haloes in which the brightest galaxy is not the central. Compared to previous studies we explicitly take into account the unrelaxed state of the host haloes, velocity offsets of halo cores and correlations between $f_{\rm BNC}$ and the satellite occupation. We find that $f_{\rm BNC}$ strongly decreases with the luminosity of the brightest galaxy and increases with the mass of the host halo. Overall, in the halo mass range $10^{13} - 10^{14.5} h^{-1} M_\odot$ we find $f_{\rm BNC} \sim 30\%$, in good agreement with a previous study by Skibba et al. We discuss the implications of these findings for studies inferring the galaxy--halo connection from satellite kinematics, models of the conditional luminosity function and galaxy formation in general.

We report here the results of searching for inelastic scattering of dark matter (initial and final state dark matter particles differ by a small mass splitting) with nucleon with the first 79.6-day of PandaX-II data (Run 9). We set the upper limits for the spin independent WIMP-nucleon scattering cross section up to a mass splitting of 300 keV/c$^2$ at two benchmark dark matter masses of 1 and 10 TeV/c$^2$.

We report new searches for the solar axions and galactic axion-like dark matter particles, using the first low-background data from PandaX-II experiment at China Jinping Underground Laboratory, corresponding to a total exposure of about $2.7\times 10^4$ kg$\cdot$day. No solar axion or galactic axion-like dark matter particle candidate has been identified. The upper limit on the axion-electron coupling ($g_{Ae}$) from the solar flux is found to be about $4.35 \times 10^{-12}$ in mass range from $10^{-5}$ to 1 keV/$c^2$ with 90\% confidence level, similar to the recent LUX result. We also report a new best limit from the $^{57}$Fe de-excitation. On the other hand, the upper limit from the galactic axions is on the order of $10^{-13}$ in the mass range from 1 keV/$c^2$ to 10 keV/$c^2$ with 90\% confidence level, slightly improved compared with the LUX.

We present the first stable release of Halotools (v0.2), a community-driven Python package designed to build and test models of the galaxy-halo connection. Halotools provides a modular platform for creating mock universes of galaxies starting from a catalog of dark matter halos obtained from a cosmological simulation. The package supports many of the common forms used to describe galaxy-halo models: the halo occupation distribution (HOD), the conditional luminosity function (CLF), abundance matching, and alternatives to these models that include effects such as environmental quenching or variable galaxy assembly bias. Satellite galaxies can be modeled to live in subhalos, or to follow custom number density profiles within their halos, including spatial and/or velocity bias with respect to the dark matter profile. The package has an optimized toolkit to make mock observations on a synthetic galaxy population, including galaxy clustering, galaxy-galaxy lensing, galaxy group identification, RSD multipoles, void statistics, pairwise velocities and others, allowing direct comparison to observations. Halotools is object-oriented, enabling complex models to be built from a set of simple, interchangeable components, including those of your own creation. Halotools has an automated testing suite and is exhaustively documented on http://halotools.readthedocs.io, which includes quickstart guides, source code notes and a large collection of tutorials. The documentation is effectively an online textbook on how to build and study empirical models of galaxy formation with Python.

We report a new search of weakly interacting massive particles (WIMPs) using the combined low background data sets in 2016 and 2017 from the PandaX-II experiment in China. The latest data set contains a new exposure of 77.1 live day, with the background reduced to a level of 0.8$\times10^{-3}$ evt/kg/day, improved by a factor of 2.5 in comparison to the previous run in 2016. No excess events were found above the expected background. With a total exposure of 5.4$\times10^4$ kg day, the most stringent upper limit on spin-independent WIMP-nucleon cross section was set for a WIMP with mass larger than 100 GeV/c$^2$, with the lowest exclusion at 8.6$\times10^{-47}$ cm$^2$ at 40 GeV/c$^2$.

Dark matter halo clustering depends not only on halo mass, but also on other properties such as concentration and shape. This phenomenon is known broadly as assembly bias. We explore the dependence of assembly bias on halo definition, parametrized by spherical overdensity parameter, $\Delta$. We summarize the strength of concentration-, shape-, and spin-dependent halo clustering as a function of halo mass and halo definition. Concentration-dependent clustering depends strongly on mass at all $\Delta$. For conventional halo definitions ($\Delta \sim 200\mathrm{m}-600\mathrm{m}$), concentration-dependent clustering at low mass is driven by a population of haloes that is altered through interactions with neighbouring haloes. Concentration-dependent clustering can be greatly reduced through a mass-dependent halo definition with $\Delta \sim 20\mathrm{m}-40\mathrm{m}$ for haloes with $M_{200\mathrm{m}} \lesssim 10^{12}\, h^{-1}\mathrm{M}_{\odot}$. Smaller $\Delta$ implies larger radii and mitigates assembly bias at low mass by subsuming altered, so-called backsplash haloes into now larger host haloes. At higher masses ($M_{200\mathrm{m}} \gtrsim 10^{13}\, h^{-1}\mathrm{M}_{\odot}$) larger overdensities, $\Delta \gtrsim 600\mathrm{m}$, are necessary. Shape- and spin-dependent clustering are significant for all halo definitions that we explore and exhibit a relatively weaker mass dependence. Generally, both the strength and the sense of assembly bias depend on halo definition, varying significantly even among common definitions. We identify no halo definition that mitigates all manifestations of assembly bias. A halo definition that mitigates assembly bias based on one halo property (e.g., concentration) must be mass dependent. The halo definitions that best mitigate concentration-dependent halo clustering do not coincide with the expected average splashback radii at fixed halo mass.

We present the survey strategy and early results of the "Satellites Around Galactic Analogs" (SAGA) Survey. The SAGA Survey's goal is to measure the distribution of satellite galaxies around 100 systems analogous to the Milky Way down to the luminosity of the Leo I dwarf galaxy ($ M_r < -12.3 $). We define a Milky Way analog based on $K$-band luminosity and local environment. Here, we present satellite luminosity functions for 8 Milky Way analog galaxies between 20 to 40 Mpc. These systems have nearly complete spectroscopic coverage of candidate satellites within the projected host virial radius down to $ r_o < 20.75 $ using low redshift $gri$ color criteria. We have discovered a total of 25 new satellite galaxies: 14 new satellite galaxies meet our formal criteria around our complete host systems, plus 11 additional satellites in either incompletely surveyed hosts or below our formal magnitude limit. Combined with 13 previously known satellites, there are a total of 27 satellites around 8 complete Milky Way analog hosts. We find a wide distribution in the number of satellites per host, from 1 to 9, in the luminosity range for which there are five Milky Way satellites. Standard abundance matching extrapolated from higher luminosities predicts less scatter between hosts and a steeper luminosity function slope than observed. We find that the majority of satellites (26 of 27) are star-forming. These early results indicate that the Milky Way has a different satellite population than typical in our sample, potentially changing the physical interpretation of measurements based only on the Milky Way's satellite galaxies.

Dwarf galaxies are known to have remarkably low star formation efficiency due to strong feedback. Adopting the dwarf galaxies of the Milky Way as a laboratory, we explore a flexible semi-analytic galaxy formation model to understand how the feedback processes shape the satellite galaxies of the Milky Way. Using Markov-Chain Monte-Carlo, we exhaustively search a large parameter space of the model and rigorously show that the general wisdom of strong outflows as the primary feedback mechanism cannot simultaneously explain the stellar mass function and the mass--metallicity relation of the Milky Way satellites. An extended model that assumes that a fraction of baryons is prevented from collapsing into low-mass halos in the first place can be accurately constrained to simultaneously reproduce those observations. The inference suggests that two different physical mechanisms are needed to explain the two different data sets. In particular, moderate outflows with weak halo mass dependence are needed to explain the mass--metallicity relation, and prevention of baryons falling into shallow gravitational potentials of low-mass halos (e.g. "pre-heating") is needed to explain the low stellar mass fraction for a given subhalo mass.

The Giant Radio Array for Neutrino Detection (GRAND) aims at detecting ultra-high energy extraterrestrial neutrinos via the extensive air showers induced by the decay of tau leptons created in the interaction of neutrinos under the Earth's surface. Consisting of an array of $\sim10^5$ radio antennas deployed over $\sim 2\times10^5\,\rm {km}^2$, GRAND plans to reach, for the first time, an all-flavor sensitivity of $\sim1.5\times10^{-10} \,\rm GeV\, cm^{-2} \,s^{-1}\, sr^{-1}$ above $5\times10^{17}$ eV and a sub-degree angular resolution, beyond the reach of other planned detectors. We describe here preliminary designs and simulation results, plans for the ongoing, staged approach to the construction of GRAND, and the rich research program made possible by GRAND's design sensitivity and angular resolution.

Empirical models of galaxy formation require assumptions about the correlations between galaxy and halo properties. These may be calibrated against observations or inferred from physical models such as hydrodynamical simulations. In this Letter, we use the EAGLE simulation to investigate the correlation of galaxy size with halo properties. We motivate this analysis by noting that the common assumption of angular momentum partition between baryons and dark matter in rotationally supported galaxies overpredicts both the spread in the stellar mass-size relation and the anticorrelation of size and velocity residuals, indicating a problem with the galaxy-halo connection it implies. We find the EAGLE galaxy population to perform significantly better on both statistics, and trace this success to the weakness of the correlations of galaxy size with halo mass, concentration and spin at fixed stellar mass. Using these correlations in empirical models will enable fine-grained aspects of galaxy scalings to be matched.

We measure how the properties of star-forming central galaxies correlate with large-scale environment, $\delta$, measured on $10$Mpc/h scales. We use group catalogs to isolate a robust sample of central galaxies with high purity and completeness. The properties we investigate are star formation rate (SFR), exponential disk scale length $R_{\rm exp}$, and Sersic index of the light profile, $n$. We find that, at all stellar masses, there is an inverse correlation between SFR and $\delta$, meaning that above-average star forming centrals live in underdense regions. For $n$ and $R_{\rm exp}$, there is no correlation with $\delta$ at $M_{\rm star}\lesssim 10^{10.5}$ $M_\odot$, but at higher masses there are positive correlations; a weak correlation with $R_{\rm exp}$ and a strong correlation with $n$. These data are evidence of assembly bias within the star-forming population. The results for SFR are consistent with a model in which SFR correlates with present-day halo accretion rate, $\dot{M}_h$. In this model, galaxies are assigned to halos using the abundance matching ansatz, which maps galaxy stellar mass onto halo mass. At fixed halo mass, SFR is assigned to galaxies using the same approach,but $\dot{M}_h$ is used to map onto SFR. The best-fit model requires some scatter in the $\dot{M}_h$-SFR relation. The $R_{\rm exp}$ and $n$ measurements are consistent with a model in which these quantities are correlated with the spin parameter of the halo, $\lambda$. Halo spin does not correlate with $\delta$ at low halo masses, but for higher mass halos, high-spin halos live in higher density environments at fixed $M_h$. Put together with the earlier installments of this series, these data demonstrate that quenching processes have limited correlation with halo formation history, but the growth of active galaxies, as well as other detailed properties, are influenced by the details of halo assembly.

Galaxy clustering on small scales is significantly under-predicted by sub-halo abundance matching (SHAM) models that populate (sub-)haloes with galaxies based on peak halo mass, $M_{\rm peak}$. SHAM models based on the peak maximum circular velocity, $V_{\rm peak}$, have had much better success. The primary reason $M_{\rm peak}$ based models fail is the relatively low abundance of satellite galaxies produced in these models compared to those based on $V_{\rm peak}$. Despite success in predicting clustering, a simple $V_{\rm peak}$ based SHAM model results in predictions for galaxy growth that are at odds with observations. We evaluate three possible remedies that could "save" mass-based SHAM: (1) SHAM models require a significant population of "orphan" galaxies as a result of artificial disruption/merging of sub-haloes in modern high resolution dark matter simulations; (2) satellites must grow significantly after their accretion; and (3) stellar mass is significantly affected by halo assembly history. No solution is entirely satisfactory. However, regardless of the particulars, we show that popular SHAM models based on $M_{\rm peak}$ cannot be complete physical models as presented. Either $V_{\rm peak}$ truly is a better predictor of stellar mass at $z\sim 0$ and it remains to be seen how the correlation between stellar mass and $V_{\rm peak}$ comes about, or SHAM models are missing vital component(s) that significantly affect galaxy clustering.

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

Blazars are an extreme subclass of active galactic nuclei. Their rapid variability, luminous brightness, superluminal motion, and high and variable polarization are probably due to a beaming effect. However, this beaming factor (or Doppler factor) is very difficult to measure. Currently, a good way to estimate it is to use the time scale of their radio flares. In this $Letter$, we use multiwavelength data and Doppler factors reported in the literatures for a sample of 86 flaring blazars detected by Fermi to compute their intrinsic multiwavelength data and intrinsic spectral energy distributions, and investigate the correlations among observed and intrinsic data. Quite interestingly, intrinsic data show a positive correlation between luminosity and peak frequency, in contrast with the behavior of observed data, and a tighter correlation between $\gamma$-ray luminosity and the lower energy ones. For flaring blazars detected by Fermi, we conclude that (1) Observed emissions are strongly beamed; (2) The anti-correlation between luminosity and peak frequency from the observed data is an apparent result, the correlation between intrinsic data being positive; and (3) Intrinsic $\gamma$-ray luminosity is strongly correlated with other intrinsic luminosities.

New constraints are presented on the spin-dependent WIMP-nucleon interaction from the PandaX-II experiment, using a data set corresponding to a total exposure of 3.3$\times10^4$ kg-days. Assuming a standard axial-vector spin-dependent WIMP interaction with $^{129}$Xe and $^{131}$Xe nuclei, the most stringent upper limits on WIMP-neutron cross sections for WIMPs with masses above 10 GeV/c$^{2}$ are set in all dark matter direct detection experiments. The minimum upper limit of $4.1\times 10^{-41}$ cm$^2$ at 90\% confidence level is obtained for a WIMP mass of 40 GeV/c$^{2}$. This represents more than a factor of two improvement on the best available limits at this and higher masses. These improved cross-section limits provide more stringent constraints on the effective WIMP-proton and WIMP-neutron couplings.

Using group catalogs from the SDSS DR7, we attempt to measure galactic conformity in the local universe. We measure the quenched fraction of neighbor galaxies around isolated primary galaxies, dividing the isolated sample into star-forming and quiescent objects. We restrict our measurements to scales $>1$ Mpc to probe the correlations between the formation histories of distinct halos. Over the stellar mass range $10^{9.7} \le M_\ast/M_\odot \le 10^{10.9}$, we find minimal statistical evidence for conformity. We further compare these data to predictions of the halo age-matching model, in which the oldest galaxies are associated with the oldest halos at fixed $M_\ast$. For models with strong correlations between halo and stellar age, the conformity signal is too large to be consistent with the data. For weaker implementations of age-matching, galactic conformity is not a sensitive diagnostic of halo assembly bias, and would not produce a detectable signal in SDSS data. We reproduce the results of Kauffmann et al 2013, in which the star formation rates of neighbor galaxies are significantly reduced around primary galaxies when the primaries are themselves low star formers. However, we find this result is mainly driven by contamination in the isolation criterion, when using our group catalog to remove the small fraction of satellite galaxies in the sample, the conformity signal largely goes away. Lastly, we show that small conformity signals, i.e., 2-5% differences in the quenched fractions of neighbor galaxies, can be produced by mechanisms other than halo assembly bias. For example, if passive galaxies occupy more massive halos than star forming galaxies of the same stellar mass, a conformity signal that is consistent with recent measurements from PRIMUS (Berti et al 2016) can be produced.

Empirical methods for connecting galaxies to their dark matter halos have become essential for interpreting measurements of the spatial statistics of galaxies. In this work, we present a novel approach for parameterizing the degree of concentration dependence in the abundance matching method. This new parameterization provides a smooth interpolation between two commonly used matching proxies: the peak halo mass and the peak halo maximal circular velocity. This parameterization controls the amount of dependence of galaxy luminosity on halo concentration at a fixed halo mass. Effectively this interpolation scheme enables abundance matching models to have adjustable assembly bias in the resulting galaxy catalogs. With the new 400 Mpc/h DarkSky Simulation, whose larger volume provides lower sample variance, we further show that low-redshift two-point clustering and satellite fraction measurements from SDSS can already provide a joint constraint on this concentration dependence and the scatter within the abundance matching framework.

We report a new ultra-faint stellar system found in Dark Energy Camera data from the first observing run of the Magellanic Satellites Survey (MagLiteS). MagLiteS J0644-5953 (Pictor II or Pic II) is a low surface brightness ({\mu} = 28.5 mag arcsec$^{-2}$ within its half-light radius) resolved overdensity of old and metal-poor stars located at a heliocentric distance of 45 kpc. The physical size (r$_{1/2}$ = 46 pc) and low luminosity (Mv = -3.2 mag) of this satellite are consistent with the locus of spectroscopically confirmed ultra-faint galaxies. MagLiteS J0644-5953 (Pic II) is located 11.3 kpc from the Large Magellanic Cloud (LMC), and comparisons with simulation results in the literature suggest that this satellite was likely accreted with the LMC. The close proximity of MagLiteS J0644-5953 (Pic II) to the LMC also makes it the most likely ultra-faint galaxy candidate to still be gravitationally bound to the LMC.

We test whether halo age and galaxy age are correlated at fixed halo and galaxy mass. The formation histories, and thus ages, of dark matter halos correlate with their large-scale density $\rho$, an effect known as assembly bias. We test whether this correlation extends to galaxies by measuring the dependence of galaxy stellar age on $\rho$. To clarify the comparison between theory and observation, and to remove the strong environmental effects on satellites, we use galaxy group catalogs to identify central galaxies and measure their quenched fraction, $f_Q$, as a function of large-scale environment. Models that match halo age to central galaxy age predict a strong positive correlation between $f_Q$ and $\rho$. However, we show that the amplitude of this effect depends on the definition of halo age: assembly bias is significantly reduced when removing the effects of splashback halos---those halos that are central but have passed through a larger halo or experienced strong tidal encounters. Defining age using halo mass at its peak value rather than current mass removes these effects. In SDSS data, at M$_{\rm gal}\gtrsim 10^{10.0}$ M_sol/h$^2$, there is a $\sim 5\%$ increase in $f_Q$ from low to high densities, which is in agreement with predictions of dark matter halos using peak halo mass. At lower stellar mass there is little to no correlation of $f_Q$ with $\rho$. For these galaxies, age-matching is inconsistent with the data across the wide range the halo formation metrics that we tested. This implies that halo formation history has a small but statistically significant impact on quenching of star formation at high masses, while the quenching process in low-mass central galaxies is uncorrelated with halo formation history.

Many properties of the Milky Way's dark matter halo, including its mass assembly history, concentration, and subhalo population, remain poorly constrained. We explore the connection between these properties of the Milky Way and its satellite galaxy population, especially the implication of the presence of the Magellanic Clouds for the properties of the Milky Way halo. Using a suite of high-resolution $N$-body simulations of Milky Way-mass halos with a fixed final Mvir ~ 10^{12.1}Msun, we find that the presence of Magellanic Cloud-like satellites strongly correlates with the assembly history, concentration, and subhalo population of the host halo, such that Milky Way-mass systems with Magellanic Clouds have lower concentration, more rapid recent accretion, and more massive subhalos than typical halos of the same mass. Using a flexible semi-analytic galaxy formation model that is tuned to reproduce the stellar mass function of the classical dwarf galaxies of the Milky Way with Markov-Chain Monte-Carlo, we show that adopting host halos with different mass-assembly histories and concentrations can lead to different best-fit models for galaxy-formation physics, especially for the strength of feedback. These biases arise because the presence of the Magellanic Clouds boosts the overall population of high-mass subhalos, thus requiring a different stellar-mass-to-halo-mass ratio to match the data. These biases also lead to significant differences in the mass--metallicity relation, the kinematics of low-mass satellites, the number counts of small satellites associated with the Magellanic Clouds, and the stellar mass of Milky Way itself. Observations of these galaxy properties can thus provide useful constraints on the properties of the Milky Way halo.

Merger trees are routinely used to follow the growth and merging history of dark matter haloes and subhaloes in simulations of cosmic structure formation. Srisawat et al. (2013) compared a wide range of merger-tree-building codes. Here we test the influence of output strategies and mass resolution on tree-building. We find that, somewhat surprisingly, building the tree from more snapshots does not generally produce more complete trees; instead, it tends to short- en them. Significant improvements are seen for patching schemes which attempt to bridge over occasional dropouts in the underlying halo catalogues or schemes which combine the halo-finding and tree-building steps seamlessly. The adopted output strategy does not affec- t the average number of branches (bushiness) of the resultant merger trees. However, mass resolution has an influence on both main branch length and the bushiness. As the resolution increases, a halo with the same mass can be traced back further in time and will encounter more small progenitors during its evolutionary history. Given these results, we recommend that, for simulations intended as precursors for galaxy formation models where of order 100 or more snapshots are analysed, the tree-building routine should be integrated with the halo finder, or at the very least be able to patch over multiple adjacent snapshots.

Xinglong Observatory of the National Astronomical Observatories, Chinese Academy of Sciences (NAOC), is one of the major optical observatories in China, which hosts nine optical telescopes including the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and the 2.16 m reflector. Scientific research from these telescopes is focused on stars, galaxies, and exoplanets using multicolor photometry and spectroscopic observations. Therefore, it is important to provide the observing conditions of the site, in detail, to the astronomers for an efficient use of these facilities. In this article, we present the characterization of observing conditions at Xinglong Observatory based on the monitoring of meteorology, seeing and sky brightness during the period from 2007 to 2014. Results suggest that Xinglong Observatory is still a good site for astronomical observations. Our analysis of the observing conditions at Xinglong Observatory can be used as a reference to the observers on targets selection, observing strategy, and telescope operation.

We study the mass spectrum of destroyed dwarfs that contribute to the accreted stellar mass of Milky Way (MW) mass M_vir ~ 10^12.1 M_sun) halos using a suite of 45 zoom-in, dissipationless simulations. Empirical models are employed to relate (peak) subhalo mass to dwarf stellar mass, and we use constraints from z=0 observations and hydrodynamical simulations to estimate the metallicity distribution of the accreted stellar material. The dominant contributors to the accreted stellar mass are relatively massive dwarfs with M_star ~ 10^8-10^10 M_sun. Halos with more quiescent accretion histories tend to have lower mass progenitors (10^8-10^9 M_sun), and lower overall accreted stellar masses. Ultra-faint mass (M_star < 10^5 M_sun) dwarfs contribute a negligible amount (<< 1%) to the accreted stellar mass and, despite having low average metallicities, supply a small fraction (~2-5 %) of the very metal-poor stars with [Fe/H] < -2. Dwarfs with masses 10^5 < M_star/M_sun < 10^8 provide a substantial amount of the very metal-poor stellar material (~40-80 %), and even relatively metal-rich dwarfs with M_star > 10^8 M_sun can contribute a considerable fraction (~20-60 %) of metal-poor stars if their metallicity distributions have significant metal-poor tails. Finally, we find that the generic assumption of a quiescent assembly history for the MW halo seems to be in tension with the mass spectrum of its surviving dwarfs. We suggest that the MW could be a "transient fossil"; a quiescent halo with a recent accretion event(s) that disguises the preceding formation history of the halo.

We study the abundance of substructure in the matter density near galaxies using ALMA Science Verification observations of the strong lensing system SDP.81. We present a method to measure the abundance of subhalos around galaxies using interferometric observations of gravitational lenses. Using simulated ALMA observations, we explore the effects of various systematics, including antenna phase errors and source priors, and show how such errors may be measured or marginalized. We apply our formalism to ALMA observations of SDP.81. We find evidence for the presence of a $M=10^{8.96\pm 0.12} M_{\odot}$ subhalo near one of the images, with a significance of $6.9\sigma$ in a joint fit to data from bands 6 and 7; the effect of the subhalo is also detected in both bands individually. We also derive constraints on the abundance of dark matter subhalos down to $M\sim 2\times 10^7 M_{\odot}$, pushing down to the mass regime of the smallest detected satellites in the Local Group, where there are significant discrepancies between the observed population of luminous galaxies and predicted dark matter subhalos. We find hints of additional substructure, warranting further study using the full SDP.81 dataset (including, for example, the spectroscopic imaging of the lensed carbon monoxide emission). We compare the results of this search to the predictions of $\Lambda$CDM halos, and find that given current uncertainties in the host halo properties of SDP.81, our measurements of substructure are consistent with theoretical expectations. Observations of larger samples of gravitational lenses with ALMA should be able to improve the constraints on the abundance of galactic substructure.

Cosmological transverse momentum fields, whose directions are perpendicular to Fourier wave vectors, induce temperature anisotropies in the cosmic microwave background via the kinetic Sunyaev-Zeldovich (kSZ) effect. The transverse momentum power spectrum contains the four-point function of density and velocity fields, $\langle\delta\delta v v\rangle$. In the post-reionization epoch, nonlinear effects dominate in the power spectrum. We use perturbation theory and cosmological $N$-body simulations to calculate this nonlinearity. We derive the next-to-leading order expression for the power spectrum with a particular emphasis on the connected term that has been ignored in the literature. While the contribution from the connected term on small scales ($k>0.1\,h\,\rm{Mpc}^{-1}$) is subdominant relative to the unconnected term, we find that its contribution to the kSZ power spectrum at $\ell = 3000$ at $z<6$ can be as large as ten percent of the unconnected term, which would reduce the allowed contribution from the reionization epoch ($z>6$) by twenty percent. The power spectrum of transverse momentum on large scales is expected to scale as $k^2$ as a consequence of momentum conservation. We show that both the leading and the next-to-leading order terms satisfy this scaling. In particular, we find that both of the unconnected and connected terms are necessary to reproduce $k^2$.

We report the discovery of eight new ultra-faint dwarf galaxy candidates in the second year of optical imaging data from the Dark Energy Survey (DES). Six of these candidates are detected at high confidence, while two lower-confidence candidates are identified in regions of non-uniform survey coverage. The new stellar systems are found by three independent automated search techniques and are identified as overdensities of stars, consistent with the isochrone and luminosity function of an old and metal-poor simple stellar population. The new systems are faint (Mv > -4.7 mag) and span a range of physical sizes (17 pc < $r_{1/2}$ < 181 pc) and heliocentric distances (25 kpc < D < 214 kpc). All of the new systems have central surface brightnesses consistent with known ultra-faint dwarf galaxies (\mu < 27.5 mag arcsec$^{-2}$). Roughly half of the DES candidates are more distant, less luminous, and/or have lower surface brightnesses than previously known Milky Way satellite galaxies. Most of the candidates are found in the southern part of the DES footprint close to the Magellanic Clouds. We find that the DES data alone exclude (p < 0.001) a spatially isotropic distribution of Milky Way satellites and that the observed distribution can be well, though not uniquely, described by an association between several of the DES satellites and the Magellanic system. Our model predicts that the full sky may hold ~100 ultra-faint galaxies with physical properties comparable to the DES satellites and that 20-30% of these would be spatially associated with the Magellanic Clouds.

We propose a common terminology for use in describing both temporal merger trees and spatial structure trees for dark-matter halos. We specify a unified data format in HDF5 and provide example I/O routines in C, FORTRAN and PYTHON.

We report the results of a weakly-interacting massive particle (WIMP) dark matter search using the full 80.1\;live-day exposure of the first stage of the PandaX experiment (PandaX-I) located in the China Jin-Ping Underground Laboratory. The PandaX-I detector has been optimized for detecting low-mass WIMPs, achieving a photon detection efficiency of 9.6\%. With a fiducial liquid xenon target mass of 54.0\,kg, no significant excess event were found above the expected background. A profile likelihood analysis confirms our earlier finding that the PandaX-I data disfavor all positive low-mass WIMP signals reported in the literature under standard assumptions. A stringent bound on the low mass WIMP is set at WIMP mass below 10\,GeV/c$^2$, demonstrating that liquid xenon detectors can be competitive for low-mass WIMP searches.

Hierarchical structure formation implies that the number of subhalos within a dark matter halo depends not only on halo mass, but also on the formation history of the halo. This dependence on the formation history, which is highly correlated with halo concentration, can account for the super-Poissonian scatter in subhalo occupation at a fixed halo mass that has been previously measured in simulations. Here we propose a model to predict the subhalo abundance function for individual host halos, that incorporates both halo mass and concentration. We combine results of cosmological simulations with a new suite of zoom-in simulations of Milky Way-mass halos to calibrate our model. We show the model can successfully reproduce the mean and the scatter of subhalo occupation in these simulations. The implications of this correlation between subhalo abundance and halo concentration are further investigated. We also discuss cases in which inferences about halo properties can be affected if this correlation between subhalo abundance and halo concentration is ignored; in these cases our model would give a more accurate inference. We propose that with future deep surveys, satellite occupation in the low-mass regime can be used to verify the existence of halo assembly bias.

Merging haloes with similar masses (i.e., major mergers) pose significant challenges for halo finders. We compare five halo finding algorithms' (AHF, HBT, Rockstar, SubFind, and VELOCIraptor) recovery of halo properties for both isolated and cosmological major mergers. We find that halo positions and velocities are often robust, but mass biases exist for every technique. The algorithms also show strong disagreement in the prevalence and duration of major mergers, especially at high redshifts (z>1). This raises significant uncertainties for theoretical models that require major mergers for, e.g., galaxy morphology changes, size changes, or black hole growth, as well as for finding Bullet Cluster analogues. All finders not using temporal information also show host halo and subhalo relationship swaps over successive timesteps, requiring careful merger tree construction to avoid problematic mass accretion histories. We suggest that future algorithms should combine phase-space and temporal information to avoid the issues presented.

The Large sky Area Multi-Object Spectroscopic Telescope (LAMOST) General Survey is a spectroscopic survey that will eventually cover approximately half of the celestial sphere and collect 10 million spectra of stars, galaxies and QSOs. Objects both in the pilot survey and the first year general survey are included in the LAMOST First Data Release (DR1). The pilot survey started in October 2011 and ended in June 2012, and the data have been released to the public as the LAMOST Pilot Data Release in August 2012. The general survey started in September 2012, and completed its first year of operation in June 2013. The LAMOST DR1 includes a total of 1202 plates containing 2,955,336 spectra, of which 1,790,879 spectra have observed signal-to-noise S/N >10. All data with S/N>2 are formally released as LAMOST DR1 under the LAMOST data policy. This data release contains a total of 2,204,696 spectra, of which 1,944,329 are stellar spectra, 12,082 are galaxy spectra and 5,017 are quasars. The DR1 includes not only spectra, but also three stellar catalogues with measured parameters: AFGK-type stars with high quality spectra (1,061,918 entries), A-type stars (100,073 entries), and M stars (121,522 entries). This paper introduces the survey design, the observational and instrumental limitations, data reduction and analysis, and some caveats. Description of the FITS structure of spectral files and parameter catalogues is also provided.

The linear perturbation theory of inhomogeneous reionization (LPTR) has been developed as an analytical tool for predicting the global ionized fraction and large-scale power spectrum of ionized density fluctuations during reionization. In the original formulation of the LPTR, the ionization balance and radiative transfer equations are linearized and solved in Fourier space. However, the LPTR's approximation to the full solution of the radiative transfer equation is not straightforward to interpret, since the latter is most intuitively conceptualized in position space. To bridge the gap between the LPTR and the language of numerical radiative transfer, we present a new, equivalent, position-space formulation of the LPTR that clarifies the approximations it makes and facilitates its interpretation. We offer a comparison between the LPTR and the excursion-set model of reionization (ESMR), and demonstrate the built-in capability of the LPTR to explore a wide range of reionization scenarios, and to go beyond the ESMR in exploring scenarios involving X-rays.

The study of the Universe on ultra-large scales is one of the major science cases for the Square Kilometre Array (SKA). The SKA will be able to probe a vast volume of the cosmos, thus representing a unique instrument, amongst next-generation cosmological experiments, for scrutinising the Universe's properties on the largest cosmic scales. Probing cosmic structures on extremely large scales will have many advantages. For instance, the growth of perturbations is well understood for those modes, since it falls fully within the linear regime. Also, such scales are unaffected by the poorly understood feedback of baryonic physics. On ultra-large cosmic scales, two key effects become significant: primordial non-Gaussianity and relativistic corrections to cosmological observables. Moreover, if late-time acceleration is driven not by dark energy but by modifications to general relativity, then such modifications should become apparent near and above the horizon scale. As a result, the SKA is forecast to deliver transformational constraints on non-Gaussianity and to probe gravity on super-horizon scales for the first time.

There can arise ubiquitous ultra-light scalar fields in the Universe, such as the pseudo-Goldstone bosons from the spontaneous breaking of an approximate symmetry, which can make a partial contribution to the dark matter and affect the large scale structure of the Universe. While the properties of those ultra-light dark matter are heavily model dependent and can vary in a wide range, we develop a model-independent analysis to forecast the constraints on their mass and abundance using futuristic but realistic 21 cm observables as well as CMB fluctuations, including CMB lensing measurements. Avoiding the highly nonlinear regime, the 21 cm emission line spectra are most sensitive to the ultra-light dark matter with mass m ~10^{-26} eV for which the precision attainable on mass and abundance bounds can be of order of a few percent.

A halo merger tree forms the essential backbone of a semi-analytic model for galaxy formation and evolution. Recent studies have pointed out that extracting merger trees from numerical simulations of structure formation is non-trivial; different tree building algorithms can give differing merger histories. These differences should be carefully understood before merger trees are used as input for models of galaxy formation. We investigate the impact of different halo merger trees on a semi-analytic model. We find that the z=0 galaxy properties in our model show differences between trees when using a common parameter set. The star formation history of the Universe and the properties of satellite galaxies can show marked differences between trees with different construction methods. Independently calibrating the semi-analytic model for each tree can reduce the discrepancies between the z=0 global galaxy properties, at the cost of increasing the differences in the evolutionary histories of galaxies. Furthermore, the underlying physics implied can vary, resulting in key quantities such as the supernova feedback efficiency differing by factors of 2. Such a change alters the regimes where star formation is primarily suppressed by supernovae. Therefore, halo merger trees extracted from a common halo catalogue using different, but reliable, algorithms can result in a difference in the semi-analytic model. Given the uncertainties in galaxy formation physics, however, these differences may not necessarily be viewed as significant.

We report on the first dark-matter (DM) search results from PandaX-I, a low threshold dual-phase xenon experiment operating at the China Jinping Underground Laboratory. In the 37-kg liquid xenon target with 17.4 live-days of exposure, no DM particle candidate event was found. This result sets a stringent limit for low-mass DM particles and disfavors the interpretation of previously-reported positive experimental results. The minimum upper limit, $3.7\times10^{-44}$\,cm$^2$, for the spin-independent isoscalar DM-particle-nucleon scattering cross section is obtained at a DM-particle mass of 49\,GeV/c$^2$ at 90\% confidence level.

Variations in the attenuation law have a significant impact on observed spectral energy distributions for galaxies. As one important observational property for galaxies at ultraviolet and infrared wavelength bands, the correlation between infrared-to-ultraviolet luminosity ratio and ultraviolet color index (or ultraviolet spectral slope), i.e., the IRX-UV relation (or IRX-beta relation), offered a widely used recipe for correcting dust attenuation in galaxies, but the usability appears to be in doubt now because of considerable dispersion in this relation found by many studies. In this paper, on the basis of spectral synthesis modeling and spatially resolved measurements of four nearby spiral galaxies, we provide an interpretation of the deviation in the IRX-UV relation with variations in the attenuation law. From both theoretical and observational viewpoints, two components in the attenuation curve, the linear background and the 2175 Angstrom bump, are suggested to be the parameters in addition to the stellar population age (addressed in the first paper of this series) in the IRX-UV function; different features in the attenuation curve are diagnosed for the galaxies in our sample. Nevertheless, it is often difficult to ascertain the attenuation law for galaxies in actual observations. Possible reasons for preventing the successful detection of the parameters in the attenuation curve are also discussed in this paper, including the degeneracy of the linear background and the 2175 Angstrom bump in observational channels, the requirement for young and dust-rich systems to study, and the difficulty in accurate estimates of dust attenuations at different wavelength bands.

Measurements of the HI 21-cm power spectra from the reionization epoch will be influenced by the evolution of the signal along the line-of-sight direction of any observed volume. We use numerical as well as semi-numerical simulations of reionization in a cubic volume of 607 Mpc across to study this so-called light cone effect on the HI 21-cm power spectrum. We find that the light cone effect has the largest impact at two different stages of reionization: one when reionization is $\sim 20\%$ and other when it is $\sim 80\%$ completed. We find a factor of $\sim 4$ amplification of the power spectrum at the largest scale available in our simulations. We do not find any significant anisotropy in the 21-cm power spectrum due to the light cone effect. We argue that for the power spectrum to become anisotropic, the light cone effect would have to make the ionized bubbles significantly elongated or compressed along the line-of-sight, which would require extreme reionization scenarios. We also calculate the two-point correlation functions parallel and perpendicular to the line-of-sight and find them to differ. Finally, we calculate an optimum frequency bandwidth below which the light cone effect can be neglected when extracting power spectra from observations. We find that if one is willing to accept a $10 \%$ error due to the light cone effect, the optimum frequency bandwidth for $k= 0.056 \, \rm{Mpc}^{-1}$ is $\sim 7.5$ MHz. For $k = 0.15$ and $0.41 \, \rm{Mpc}^{-1}$ the optimum bandwidth is $\sim 11$ and $\sim 16$ MHz respectively.

Merger tree codes are routinely used to follow the growth and merger of dark matter haloes in simulations of cosmic structure formation. Whereas in Srisawat et. al. we compared the trees built using a wide variety of such codes here we study the influence of the underlying halo catalogue upon the resulting trees. We observe that the specifics of halo finding itself greatly influences the constructed merger trees. We find that the choices made to define the halo mass are of prime importance. For instance, amongst many potential options different finders select self-bound objects or spherical regions of defined overdensity, decide whether or not to include substructures within the mass returned and vary in their initial particle selection. The impact of these decisions is seen in tree length (the period of time a particularly halo can be traced back through the simulation), branching ratio (essentially the merger rate of \subhalos) and mass evolution. We therefore conclude that the choice of the underlying halo finder is more relevant to the process of building merger trees than the tree builder itself. We also report on some built-in features of specific merger tree codes that (sometimes) help to improve the quality of the merger trees produced.

Several direct detection experiments, including recently CDMS-II, have reported signals consistent with 5 to 10 GeV dark matter (DM) that appear to be in tension with null results from XENON and LUX experiments; these indicate a careful review of the theoretical basis, including the galactic DM velocity distribution function (VDF). We establish a VDF parameter space from DM-only cosmological simulations and illustrate that seemingly contradictory experimental results can be made consistent within this parameter space. Future experimental limits should be reported after they are marginalized over a range of VDF parameters.

We present the largest-volume (425 Mpc/h=607 Mpc on a side) full radiative transfer simulation of cosmic reionization to date. We show that there is significant additional power in density fluctuations at very large scales. We systematically investigate the effects this additional power has on the progress, duration and features of reionization, as well as on selected reionization observables. We find that comoving simulation volume of ~100 Mpc/h per side is sufficient for deriving a convergent mean reionization history, but that the reionization patchiness is significantly underestimated. We use jackknife splitting to quantify the convergence of reionization properties with simulation volume for both mean-density and variable-density sub-regions. We find that sub-volumes of ~100 Mpc/h per side or larger yield convergent reionization histories, except for the earliest times, but smaller volumes of ~50 Mpc/h or less are not well converged at any redshift. Reionization history milestones show significant scatter between the sub-volumes, of Delta z=0.6-1 for ~50 Mpc/h volumes, decreasing to Delta z=0.3-0.5 for ~100 Mpc/h volumes, and $\Delta z$~0.1 for ~200 Mpc/h volumes. If we only consider mean-density sub-regions the scatter decreases, but remains at Delta z~0.1-0.2 for the different size sub-volumes. Consequently, many potential reionization observables like 21-cm rms, 21-cm PDF skewness and kurtosis all show good convergence for volumes of ~200 Mpc/h, but retain considerable scatter for smaller volumes. In contrast, the three-dimensional 21-cm power spectra at large scales (k<0.25 h/Mpc) do not fully converge for any sub-volume size. These additional large-scale fluctuations significantly enhance the 21-cm fluctuations, which should improve the prospects of detection considerably, given the lower foregrounds and greater interferometer sensitivity at higher frequencies. (abridged)

Measuring the small primordial nonGaussianity (PNG) predicted by cosmic inflation theories may help diagnose them. The detectability of PNG by its imprint on the 21cm power spectrum from the epoch of reionization is reassessed here in terms of $f_{NL}$, the local nonlinearity parameter. We find that an optimum, multi-frequency observation by SKA can achieve $\Delta f_{NL} \sim 3$ (comparable to recent Planck CMB limits), while a cosmic-variance-limited array of this size like Omniscope can even detect $\Delta f_{NL} \sim 0.2$. This substantially revises the methods and results of previous work.

Merger trees follow the growth and merger of dark-matter haloes over cosmic history. As well as giving important insights into the growth of cosmic structure in their own right, they provide an essential backbone to semi-analytic models of galaxy formation. This paper is the first in a series to arise from the SUSSING MERGER TREES Workshop in which ten different tree-building algorithms were applied to the same set of halo catalogues and their results compared. Although many of these codes were similar in nature, all algorithms produced distinct results. Our main conclusions are that a useful merger-tree code should possess the following features: (i) the use of particle IDs to match haloes between snapshots; (ii) the ability to skip at least one, and preferably more, snapshots in order to recover subhaloes that are temporarily lost during merging; (iii) the ability to cope with (and ideally smooth out) large, temporary flucuations in halo mass. Finally, to enable different groups to communicate effectively, we defined a common terminology that we used when discussing merger trees and we encourage others to adopt the same language. We also specified a minimal output format to record the results.

One of the most promising ways to study the epoch of reionization (EoR) is through radio observations of the redshifted 21-cm line emission from neutral hydrogen. These observations are complicated by the fact that the mapping of redshifts to line-of-sight positions is distorted by the peculiar velocities of the gas. Such distortions can be a source of error if they are not properly understood, but they also encode information about cosmology and astrophysics. We study the effects of redshift space distortions on the power spectrum of 21-cm radiation from the EoR using large scale $N$-body and radiative transfer simulations. We quantify the anisotropy introduced in the 21-cm power spectrum by redshift space distortions and show how it evolves as reionization progresses and how it relates to the underlying physics. We go on to study the effects of redshift space distortions on LOFAR observations, taking instrument noise and foreground subtraction into account. We find that LOFAR should be able to directly observe the power spectrum anisotropy due to redshift space distortions at spatial scales around $k \sim 0.1$ Mpc$^{-1}$ after $\gtrsim$ 1000 hours of integration time. At larger scales, sample errors become a limiting factor, while at smaller scales detector noise and foregrounds make the extraction of the signal problematic. Finally, we show how the astrophysical information contained in the evolution of the anisotropy of the 21-cm power spectrum can be extracted from LOFAR observations, and how it can be used to distinguish between different reionization scenarios.

(ABRIDGED)The rise of cosmic structure depends upon the statistical distribution of initial density fluctuations generated by inflation. While the simplest models predict an almost perfectly Gaussian distribution, more-general models predict a level of primordial non-Gaussianity (PNG) that observations might yet be sensitive enough to detect. Recent Planck Collaboration measurements of the CMB temperature anisotropy bispectrum significantly tighten the observational limits, but they are still far from the PNG level predicted by the simplest models of inflation. Probing levels below CMB sensitivities will require other methods, such as searching for the statistical imprint of PNG on galactic halo clustering. During the epoch of reionization (EoR), the first stars and galaxies released radiation into the intergalactic medium (IGM) that created ionized patches whose large-scale geometry and evolution reflected the underlying abundance and large-scale clustering of the star-forming galaxies. This statistical connection between ionized patches in the IGM and galactic halos suggests that observing reionization may be another way to constrain PNG. We employ the linear perturbation theory of reionization and semi-analytic models based on the excursion-set formalism to model the effects of PNG on the EoR. We quantify the effects of PNG on the large-scale structure of reionization by deriving the ionized density bias, i.e. ratio of ionized atomic to total matter overdensities in Fourier space, at small wavenumber. Just as previous studies found that PNG creates a scale-dependent signature in the halo bias, so, too, we find a scale-dependent signature in the ionized density bias. Our results, which differ significantly from previous attempts in the literature to characterize this PNG signature, will be applied elsewhere to predict its observable consequences, e.g. in the cosmic 21cm background.

Multi-band photometric images from ultraviolet and optical to infrared are collected to derive spatially resolved properties of a nearby Scd type galaxy M 101. With evolutionary stellar population synthesis models, two-dimensional distributions and radial profiles of age, metallicity, dust attenuation, and star formation timescale in the form of the Sandage star formation history are obtained. When fitting with the models, we use the IRX-$A_\mathrm{FUV}$ relation, found to depend on a second parameter of birth rate b (ratio of present and past-averaged star formation rate), to constrain the dust attenuation. There are obvious parameter gradients in the disk of M101, which supports the theory of an "inside-out disk" growth scenario. Two distinct disc regions with different gradients of age and color are discovered, similar to another late-type galaxy NGC 628. The metallicity gradient of the stellar content is flatter than that of H {\sc ii} regions. The stellar disk is optically thicker inside than outside and the global dust attenuation of this galaxy is lower, compared with galaxies of similar and earlier morphological type. We highlight that a variational star formation timescale describes the real star formation history of a galaxy. The timescale increases steadily from the center to the outskirt. We also confirm that the bulge in this galaxy is a disk-like pseudobulge, whose evolution is likely to be induced by some secular processes of the small bar with relatively young age, rich metal, and much dust.

We discuss the properties of subhalos in cluster-size halos, using a high-resolution statistical sample: the Rhapsody simulations introduced in Wu et al. (2012). We demonstrate that the criteria applied to select subhalos have significant impact on the inferred properties of the sample, including the scatter in the number of subhalos, the correlation between the subhalo number and formation time, and the shape of subhalos' spatial distribution and velocity structure. We find that the number of subhalos, when selected using the peak maximum circular velocity in their histories (a property expected to be closely related to the galaxy luminosity), is uncorrelated with the formation time of the main halo. This is in contrast to the previously reported correlation from studies where subhalos are selected by the current maximum circular velocity; we show that this difference is a result of the tidal stripping of the subhalos. We also find that the dominance of the main halo and the subhalo mass fraction are strongly correlated with halo concentration and formation history. These correlations are important to take into account when interpreting results from cluster samples selected with different criteria. Our sample also includes a fossil cluster, which is presented separately and placed in the context of the rest of the sample.

The 21cm background from the epoch of reionization is a promising cosmological probe: line-of-sight velocity fluctuations distort redshift, so brightness fluctuations in Fourier space depend upon angle, which linear theory shows can separate cosmological from astrophysical information. Nonlinear fluctuations in ionization, density and velocity change this, however. The validity and accuracy of the separation scheme are tested here for the first time, by detailed reionization simulations. The scheme works reasonably well early in reionization (< 40% ionized), but not late (> 80% ionized).

We present the first results from the Rhapsody cluster re-simulation project: a sample of 96 "zoom-in" simulations of dark matter halos of 10^14.8 +- 0.05 Msun/h, selected from a 1 (Gpc/h)^3 volume. This simulation suite is the first to resolve this many halos with ~5x10^6 particles per halo in the cluster-mass regime, allowing us to statistically characterize the distribution of and correlation between halo properties at fixed mass. We focus on the properties of the main halos and how they are affected by formation history, which we track back to z=12, over five decades in mass. We give particular attention to the impact of the formation history on the density profiles of the halos. We find that the deviations from the Navarro-Frenk-White (NFW) model and the Einasto model depend on formation time. Late-forming halos tend to have considerable deviations from both models, partly due to the presence of massive subhalos, while early-forming halos deviate less but still significantly from the NFW model and are better described by the Einasto model. We find that the halo shapes depend only moderately on formation time. Departure from spherical symmetry impacts the density profiles through the anisotropic distribution of massive subhalos. Further evidence of the impact of subhalos is provided by analyzing the phase-space structure. A detailed analysis of the properties of the subhalo population in Rhapsody is presented in a companion paper.

We examine the Velocity Distribution Function (VDF) in dark matter halos from Milky Way to cluster mass scales. We identify an empirical model for the VDF with a wider peak and a steeper tail than a Maxwell--Boltzmann distribution, and discuss physical explanations. We quantify sources of scatter in the VDF of cosmological halos and their implication for direct detection of dark matter. Given modern simulations and observations, we find that the most significant uncertainty in the VDF of the Milky Way arises from the unknown radial position of the solar system relative to the dark matter halo scale radius.

Large-scale reionization simulations are described which combine the results of cosmological N-body simulations that model the evolving density and velocity fields and identify the galactic halo sources, with ray-tracing radiative transfer calculations which model the nonequilibrium ionization of the intergalactic medium. These simulations have been used to predict some of the signature effects of reionization on cosmic radiation backgrounds, including the CMB, near-IR, and redshifted 21cm backgrounds. We summarize some of our recent progress in this work, and address the question of whether observations of such signature effects can be used to distinguish the relative contributions of galaxies of different masses to reionization.

The correlation between infrared-to-ultraviolet luminosity ratio and ultraviolet color, i.e. the IRX-UV relation, was regarded as a prevalent recipe for correcting extragalactic dust attenuation. Considerable dispersion in this relation discovered for normal galaxies, however, complicates its usability. In order to investigate the cause of the dispersion, in this paper, we select five nearby spiral galaxies, and perform spatially resolved studies on each of the galaxies, with a combination of ultraviolet and infrared imaging data. We measure all positions within each galaxy and divide the extracted regions into young and evolved stellar populations. By means of this approach, we attempt to discover separate effects of dust attenuation and stellar population age on the IRX-UV relation for individual galaxies. In this work, in addition to dust attenuation, stellar population age is interpreted to be another parameter in the IRX-UV function, and the diversity of star formation histories is suggested to disperse the age effects. At the same time, strong evidence shows the necessity of more parameters in the interpretation of observational data, such as variations in attenuation/extinction law. Fractional contributions of different components to the integrated luminosities of the galaxies suggest that the integrated measurements of galaxies which comprise different populations would weaken the effect of the age parameter on IRX-UV diagrams. The dependance of the IRX-UV relation on luminosity and radial distance in galaxies presents weak trends, which offers an implication of selective effects. The two-dimensional maps of the UV color and the infrared-to-ultraviolet ratio are displayed and show a disparity in the spatial distributions between the two parameters in galaxies, which offers a spatial interpretation of the scatter in the IRX-UV relation.

The intergalactic medium was reionized before redshift z~6, most likely by starlight which escaped from early galaxies. The very first stars formed when hydrogen molecules (H2) cooled gas inside the smallest galaxies, minihalos of mass between 10^5 and 10^8 solar masses. Although the very first stars began forming inside these minihalos before redshift z~40, their contribution has, to date, been ignored in large-scale simulations of this cosmic reionization. Here we report results from the first reionization simulations to include these first stars and the radiative feedback that limited their formation, in a volume large enough to follow the crucial spatial variations that influenced the process and its observability. We show that, while minihalo stars stopped far short of fully ionizing the universe, reionization began much earlier with minihalo sources than without, and was greatly extended, which boosts the intergalactic electron-scattering optical depth and the large-angle polarization fluctuations of the cosmic microwave background significantly. Although within current WMAP uncertainties, this boost should be readily detectable by Planck. If reionization ended as late as z_ov<~7, as suggested by other observations, Planck will thereby see the signature of the first stars at high redshift, currently undetectable by other probes.

Observations of redshifted 21-cm radiation from neutral hydrogen during the epoch of reionization (EoR) are considered to constitute the most promising tool to probe that epoch. One of the major goals of the first generation of low frequency radio telescopes is to measure the 3D 21-cm power spectrum. However, the 21-cm signal could evolve substantially along the line of sight (LOS) direction of an observed 3D volume, since the received signal from different planes transverse to the LOS originated from different look-back times and could therefore be statistically different. Using numerical simulations we investigate this so-called light cone effect on the spherically averaged 3D 21-cm power spectrum. For this version of the power spectrum, we find that the effect mostly `averages out' and observe a smaller change in the power spectrum compared to the amount of evolution in the mean 21-cm signal and its rms variations along the LOS direction. Nevertheless, changes up to 50% at large scales are possible. In general the power is enhanced/suppressed at large/small scales when the effect is included. The cross-over mode below/above which the power is enhanced/suppressed moves toward larger scales as reionization proceeds. When considering the 3D power spectrum we find it to be anisotropic at the late stages of reionization and on large scales. The effect is dominated by the evolution of the ionized fraction of hydrogen during reionization and including peculiar velocities hardly changes these conclusions. We present simple analytical models which explain qualitatively all the features we see in the simulations.

The peculiar velocity of the intergalactic gas responsible for the cosmic 21cm background from the epoch of reionization and beyond introduces an anisotropy in the three-dimensional power spectrum of brightness temperature fluctuations. Measurement of this anisotropy by future 21cm surveys is a promising tool for separating cosmology from 21cm astrophysics. However, previous attempts to model the signal have often neglected peculiar velocity or only approximated it crudely. This paper re-examines the effects of peculiar velocity on the 21cm signal in detail, improving upon past treatment and addressing several issues for the first time. (1) We show that properly accounting for finite optical depth eliminates the unphysical divergence of 21cm brightness temperature in overdense regions of the IGM found by previous work that employed the usual optically-thin approximation. (2) The approximation made previously to circumvent the diverging brightness temperature problem by capping velocity gradient can misestimate the power spectrum on all scales. (3) The observed power spectrum in redshift-space remains finite even in the optically-thin approximation if one properly accounts for the redshift-space distortion. However, results that take full account of finite optical depth show that this approximation is only accurate in the limit of high spin temperature. (4) The linear theory for redshift-space distortion results in ~30% error in the observationally relevant wavenumber range, at the 50% ionized epoch. (5) We describe and test two numerical schemes to calculate the 21cm signal from reionization simulations to incorporate peculiar velocity effects in the optically-thin approximation accurately. One is particle-based, the other grid-based, and while the former is most accurate, we demonstrate that the latter is computationally more efficient and can achieve sufficient accuracy. [Abridged]

The AGN-host co-evolution issue is investigated here by focusing on the evolution of the [\ion{O}{3}]$\lambda5007$ emission-line profile. In order to simultaneously measure both [\ion{O}{3}] line profile and circumnuclear stellar population in individual spectrum, a large sample of narrow emission-line galaxies is selected from the MPA/JHU SDSS DR7 catalog. By requiring that 1) the [\ion{O}{3}] line signal-to-noise ratio is larger than 30, 2) the [\ion{O}{3}] line width is larger than the instrumental resolution by a factor of 2, our sample finally contains 2,333 Seyfert galaxies/LINERs (AGNs), 793 transition galaxies, and 190 starforming galaxies. In additional to the commonly used profile parameters (i.e., line centroid, relative velocity shift and velocity dispersion), two dimensionless shape parameters, skewness and kurtosis, are used to quantify the line shape deviation from a pure Gaussian function. We show that the transition galaxies are systematically associated with narrower line widths and weaker [\ion{O}{3}] broad wings than the AGNs, which implies that the kinematics of the emission-line gas is different in the two kinds of objects. By combining the measured host properties and line shape parameters, we find that the AGNs with stronger blue asymmetries tend to be associated with younger stellar populations. However, the similar trend is not identified in the transition galaxies. The failure is likely resulted from a selection effect in which the transition galaxies are systematically associated with younger stellar populations than the AGNs. The evolutionary significance revealed here suggests that both NLR kinematics and outflow feedback in AGNs co-evolve with their host galaxies.

The prospect of detecting the first galaxies by observing their impact on the intergalactic medium as they reionized it during the first billion years leads us to ask whether such indirect observations are capable of diagnosing which types of galaxies were most responsible for reionization. We attempt to answer this by considering a set of large-scale radiative transfer simulations of reionization in sufficiently large volumes to make statistically meaningful predictions of observable signatures, while also directly resolving all atomically-cooling halos down to 10^8 M_solar. We focus here on predictions of the 21-cm background, to see if upcoming observations are capable of distinguishing a universe ionized primarily by high-mass halos from one in which both high-mass and low-mass halos are responsible, and to see how these results depend upon the uncertain source efficiencies. We find that 21-cm fluctuation power spectra observed by the first generation EoR/21-cm radio interferometer arrays should be able to distinguish the case of reionization by high-mass halos alone from that by both high- and low-mass halos, together. Some reionization scenarios yield very similar power spectra and rms evolution and thus can only be discriminated by their different mean reionization history and 21-cm PDF distributions. We find that the skewness of the 21-cm PDF distribution smoothed over LOFAR-like window shows a clear feature correlated with the rise of the rms due to patchiness. Measurements of the mean photoionization rates are sensitive to the average density of the regions being studied and therefore could be strongly skewed in certain cases. (abridged)

SDSS J094857.3+002225 is a very radio-loud narrow-line Seyfert 1 (NLS1) galaxy. Here, we report our discovery of the intranight optical variability (INOV) of this galaxy through the optical monitoring in the B and R bands that covered seven nights in 2009. Violent rapid variability in the optical bands was identified in this RL-NLS1 for the first time, and the amplitudes of the INOV reaches 0.5 mag in both the B and R bands on the timescale of several hours. The detection of the INOV provides a piece of strong evidence supporting the fact that the object carries a relativistic jet with a small viewing angle, which confirms the conclusion drawn from the previous multi-wavelength studies.

We introduce a new physical parameter, the optical variability amplitude, to the well-established Eigenvector 1 space of quasars and test a sample of long-term B-band light curves of 42 PG quasars monitored by Giveon et al. (1999). We find that the optical variability amplitude strongly correlates with the intensity ratio of Fe II to H$\beta$, H$\beta$ width and peak luminosity at 5007$\rm{\AA}$. We briefly discuss the physical meaning of our findings and suggest that the Eddington ratio may be a key factor in determining a quasar's variability.

We present our analysis of the three HII-buried-AGN: SDSS J091053+333008, SDSS J121837+091324, and SDSS J153002-020415, by studying their optical spectra extracted from SDSS. The location in the BPT diagnostic diagrams of the three galaxies indicates that the narrow emission lines are mainly exited from HII regions. However, after the removal of the host galaxy's stellar emission, the emission lines display the typical feature of Narrow-line Seyfert 1-like. All of the three objects have large Eddington ratio, small black hole mass, and low star formation rate. We propose that the three galaxies are at the transit stage from the starburst-dominated phase to AGN-dominated phase.

Three-dimensional neutral hydrogen mapping using the redshifted 21 cm line has recently emerged as a promising cosmological probe. Within the framework of slow-roll reconstruction, we analyze how well the inflationary potential can be reconstructed by combining data from 21 cm experiments and cosmic microwave background data from the Planck satellite. We consider inflationary models classified according to the amplitude of their tensor component, and show that 21 cm measurements can significantly improve constraints on the slow-roll parameters and determine the shape of the inflationary potential.

We have mapped the central region of the HH 111 protostellar system in 1.33 mm continuum, C18O(J=2-1), 13CO (J=2-1), and SO (N_J=5_6-4_5) emission at ~3" resolution with the Submillimeter Array. There are two sources, VLA 1 (=IRAS 05491+0247) and VLA 2, with the VLA 1 source driving the HH 111 jet. Thermal emission is seen in 1.33 mm continuum tracing the dust in the envelope and the putative disks around the sources. A flattened, torus-like envelope is seen in C18O and 13CO around the VLA 1 source surrounding the dust lane perpendicular to the jet axis, with an inner radius of ~ 400 AU (1"), an outer radius of ~ 3200 AU (8"), and a thickness of ~ 1000 AU (2.5"). It seems to be infalling toward the center with conservation of specific angular momentum rather than with a Keplerian rotation as assumed by Yang et al. 1997. An inner envelope is seen in SO, with a radius of ~ 500 AU (1.3"). The inner part of this inner envelope, which is spatially coincident with the dust lane, seems to have a differential rotation and thus may have formed a rotationally supported disk. The outer part of this inner envelope, however, may have a rotation velocity decreasing toward the center and thus represent a region where an infalling envelope is in transition to a rotationally supported disk. A brief comparison with a collapsing model suggests that the flattened, torus-like envelope seen in C18O and 13CO could result from a collapse of a magnetized rotating toroid.

We examine the accretion properties in a sample of 42 hard (3-60keV) X-ray selected nearby broad-line AGNs. The energy range in the sample is harder than that usually used in the similar previous studies. These AGNs are mainly complied from the RXTE All Sky Survey (XSS), and complemented by the released INTEGRAL AGN catalog. The black hole masses, bolometric luminosities of AGN, and Eddington ratios are derived from their optical spectra in terms of the broad H$\beta$ emission line. The tight correlation between the hard X-ray (3-20keV) and bolometric/line luminosity is well identified in our sample. Also identified is a strong inverse Baldwin relationship of the H$\beta$ emission line. In addition, all these hard X-ray AGNs are biased toward luminous objects with high Eddington ratio (mostly between 0.01 to 0.1) and low column density ($<10^{22} \mathrm{cm^{-2}}$), which is most likely due to the selection effect of the surveys. The hard X-ray luminosity is consequently found to be strongly correlated with the black hole mass. We believe the sample completeness will be improved in the next few years by the ongoing Swift and INTEGRAL missions, and by the next advanced missions, such as NuSTAR, Simbol-X, and NeXT. Finally, the correlation between RFe (=optical FeII/H$\beta$) and disk temperature as assessed by $T\propto (L/L_{\mathrm{Edd}})M_{\mathrm{BH}}^{-1}$ leads us to suggest that the strength of the FeII emission is mainly determined by the shape of the ionizing spectrum.

We use astrophysical data to shed light on fundamental physics by constraining parametrized theoretical cosmological and gravitational models. Gravitational parameters are those constants that parametrize possible departures from Einstein's general theory of relativity. We develop a general framework to describe torsion in the spacetime around the Earth, and show that certain observables of the Gravity Probe B experiment can be computed in this framework. We also search for viable theories of gravity where the Ricci scalar R in the Lagrangian is replaced by an arbitrary function f(R). Making use of the equivalence between such theories and scalar-tensor gravity, we find that models can be made consistent with solar system constraints either by giving the scalar a high mass or by exploiting the so-called Chameleon Effect. Cosmology can successfully describe the evolution of our universe using six or more adjustable cosmological parameters. There is growing interest in using 3-dimensional neutral hydrogen mapping with the redshifted 21 cm line as a cosmological probe. We quantify how the precision with which cosmological parameters can be measured depends on a broad range of assumptions. We present an accurate and robust method for measuring cosmological parameters that exploits the fact that the ionization power spectra are rather smooth functions that can be accurately fit by 7 phenomenological parameters. We find that a future square kilometer array optimized for 21 cm tomography could have great potential, improving the sensitivity to spatial curvature and neutrino masses by up to two orders of magnitude, to Delta-Omega_k ~ 0.0002 and Delta-m_nu ~ 0.007 eV, and giving a 4 sigma detection of the spectral index running predicted by the simplest inflation models.

There is growing interest in using 3-dimensional neutral hydrogen mapping with the redshifted 21 cm line as a cosmological probe, as it has been argued to have a greater long-term potential than the cosmic microwave background. However, its utility depends on many assumptions. To aid experimental planning and design, we quantify how the precision with which cosmological parameters can be measured depends on a broad range of assumptions. We cover assumptions related to modeling of the ionization power spectrum and associated nonlinearity, experimental specifications like array layout and noise, cosmological assumptions about reionization history and inter-galactic medium (IGM) evolution, and assumptions about astrophysical foregrounds. We derive simple analytic approximations for how various assumptions affect the results, and find that ionization power modeling is most important, followed by array layout (crudely, the more compact, the better). We also present an accurate yet robust method for measuring cosmological parameters in practice, separating the physics from the astrophysics by exploiting both gravitationally induced clustering anisotropy and the fact that the ionization power spectra are rather smooth functions that can be accurately fit by 7 phenomenological parameters. For example, a future square kilometer array optimized for 21 cm tomography could improve the sensitivity of the Planck CMB satellite to spatial curvature and neutrino masses by up to two orders of magnitude, to Delta-Omega_k ~ 0.0002 and Delta m_nu ~ 0.007 eV, and give a 4 sigma detection of the spectral index running predicted by the simpliest inflation models.

It is well-entrenched folklore that torsion gravity theories predict observationally negligible torsion in the solar system, since torsion (if it exists) couples only to the intrinsic spin of elementary particles, not to rotational angular momentum. We argue that this assumption has a logical loophole which can and should be tested experimentally. In the spirit of action=reaction, if a rotating mass like a planet can generate torsion, then a gyroscope should also feel torsion. Using symmetry arguments, we show that to lowest order, the torsion field around a uniformly rotating spherical mass is determined by seven dimensionless parameters. These parameters effectively generalize the PPN formalism and provide a concrete framework for further testing GR. We construct a parametrized Lagrangian that includes both standard torsion-free GR and Hayashi- Shirafuji maximal torsion gravity as special cases. We demonstrate that classic solar system tests rule out the latter and constrain two observable parameters. We show that Gravity Probe B (GPB) is an ideal experiment for further constraining torsion theories, and work out the most general torsion-induced precession of its gyroscope in terms of our torsion parameters

We search for viable f(R) theories of gravity, making use of the equivalence between such theories and scalar-tensor gravity. We find that models can be made consistent with solar system constraints either by giving the scalar a high mass or by exploiting the so-called chameleon effect. However, in both cases, it appears likely that any late-time cosmic acceleration will be observationally indistinguishable from acceleration caused by a cosmological constant. We also explore further observational constraints from, e.g., big bang nucleosynthesis and inflation.