Source author record

Xuelei Chen

Xuelei Chen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

87works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

87 published item(s)

preprint2026arXiv

The Tianlai-WIYN North Celestial Cap Redshift Survey

We present the results of a small, low redshift spectroscopic survey of galaxies within 3 degrees of the North Celestial Pole (NCP) selected using V-band photometry obtained from the North Celestial Cap Survey (NCCS) (Gorbikov & Brosch 2014). The purpose of the current survey is to create a redshift space template for 21 cm emission from neutral hydrogen with which to correlate radio line intensity observations by the Tianlai dish and cylinder interferometers. A total of 898 redshifts were obtained from the 2102 extended objects in the NCCS with m_V < 19 in the survey area. After accounting for extinction, the survey geometry and selection effects, the number density and clustering pattern of galaxies in the redshift catalog are consistent with other low redshift surveys. We were also able to identify 11 galaxy cluster candidates from this redshift catalog.

preprint2025arXiv

Introduction to the Chinese Space Station Survey Telescope (CSST)

The Chinese Space Station Survey Telescope (CSST) is an upcoming Stage-IV sky survey telescope, distinguished by its large field of view (FoV), high image quality, and multi-band observation capabilities. It can simultaneously conduct precise measurements of the Universe by performing multi-color photometric imaging and slitless spectroscopic surveys. The CSST is equipped with five scientific instruments, i.e. Multi-band Imaging and Slitless Spectroscopy Survey Camera (SC), Multi-Channel Imager (MCI), Integral Field Spectrograph (IFS), Cool Planet Imaging Coronagraph (CPI-C), and THz Spectrometer (TS). Using these instruments, CSST is expected to make significant contributions and discoveries across various astronomical fields, including cosmology, galaxies and active galactic nuclei (AGN), the Milky Way and nearby galaxies, stars, exoplanets, Solar System objects, astrometry, and transients and variable sources. This review aims to provide a comprehensive overview of the CSST instruments, observational capabilities, data products, and scientific potential.

preprint2022arXiv

Anisotropies of Cosmic Optical and Near-IR Background from China Space Station Telescope (CSST)

Anisotropies of the cosmic optical background (COB) and cosmic near-IR background (CNIRB) are capable of addressing some of the key questions in cosmology and astrophysics. In this work, we measure and analyze the angular power spectra of the simulated COB and CNIRB in the ultra-deep field of the China Space Station Telescope (CSST-UDF). The CSST-UDF covers about 9 square degrees, with magnitude limits ~28.3, 28.2, 27.6, 26.7 AB mag for point sources with 5-sigma detection in the r (0.620 um), i (0.760 um), z (0.915 um), and y (0.965 um) bands, respectively. According to the design parameters and scanning pattern of the CSST, we generate mock data, merge images and mask the bright sources in the four bands. We obtain four angular power spectra from l=200 to 2,000,000 (from arcsecond to degree), and fit them with a multi-component model including intrahalo light (IHL) using the Markov chain Monte Carlo (MCMC) method. We find that the signal-to-noise ratio (SNR) of the IHL is larger than 8 over the range of angular scales that are useful for astrophysical studies (l~10,000-400,000). Comparing to previous works, the constraints on the model parameters are improved by factors of 3~4 in this study, which indicates that the CSST-UDF survey can be a powerful probe on the cosmic optical and near-IR backgrounds.

preprint2022arXiv

Constraining Brans-Dicke cosmology with the CSST galaxy clustering spectroscopic survey

The Brans-Dicke (BD) theory is the simplest Scalar-Tensor theory of gravity, which can be considered as a candidate of modified Einstein's theory of general relativity. In this work, we forecast the constraints on BD theory in the CSST galaxy clustering spectroscopic survey with a magnitude limit $\sim 23$ AB mag for point-source 5$σ$ detection. We generate mock data based on the zCOSMOS catalog and consider the observational and instrumental effects of the CSST spectroscopic survey. We predicate galaxy power spectra in the BD theory from $z=0$ to 1.5, and the galaxy bias and other systematical parameters are also included. The Markov Chain Monte Carlo (MCMC) technique is employed to find the best-fits and probability distributions of the cosmological and systematical parameters. A Brans-Dicke parameter $ζ$ is introduced, which satisfies $ζ=\ln \left(1+\frac{1}ω\right)$. We find that the CSST spectroscopic galaxy clustering survey can give $|ζ|<10^{-2}$, or equivalently $|ω|>\mathcal{O}(10^2)$ and $|\dot{G}/G|<10^{-13}$, under the assumption $ζ= 0$. These constraints are almost at the same order of magnitude compared to the joint constraints using the current cosmic microwave background (CMB), baryon acoustic oscillations (BAO), and Type Ia supernova (SN Ia) data, indicating that the CSST galaxy clustering spectroscopic survey would be powerful to constrain the BD theory and other modified gravity theories.

preprint2022arXiv

Effects of Small-Scale Absorption Systems on Neutral Islands during the Late Epoch of Reionization

The reionization process is expected to be prolonged by the small-scale absorbers (SSAs) of ionizing photons, which have been seen as Lyman-limit systems in quasar absorption line observations. We use a set of semi-numerical simulations to investigate the effects of absorption systems on the reionization process, especially their impacts on the neutral islands during the late Epoch of Reionization (EoR). Three models are studied, i.e. the extreme case of no-SSA model with a high level of ionizing background, the moderate-SSA model with a relatively high level of ionizing background, and the dense-SSA model with a low level of ionizing background. We find that while the characteristic scale of neutral regions decreases during the early and middle stages of reionization, it stays nearly unchanged at about 10 comoving Mpc during the late stage for the no-SSA and moderate-SSA models. However, in the case of weak ionizing background in the dense-SSA model, the characteristic island scale shows obvious evolution, as large islands break into many small ones that are slowly ionized. The evolutionary behavior of neutral islands during the late EoR thus provides a novel way to constrain the abundance of SSAs. We discuss the 21-cm observation with the upcoming Square Kilometre Array. The different models can be distinguished by either the 21-cm imaging or the 21-cm power spectrum measurements.

preprint2022arXiv

Extracting Photometric Redshift from Galaxy Flux and Image Data using Neural Networks in the CSST Survey

The accuracy of galaxy photometric redshift (photo-$z$) can significantly affect the analysis of weak gravitational lensing measurements, especially for future high-precision surveys. In this work, we try to extract photo-$z$ information from both galaxy flux and image data expected to be obtained by China Space Station Telescope (CSST) using neural networks. We generate mock galaxy images based on the observational images from the Advanced Camera for Surveys of Hubble Space Telescope (HST-ACS) and COSMOS catalogs, considering the CSST instrumental effects. Galaxy flux data are then measured directly from these images by aperture photometry. The Multi-Layer Perceptron (MLP) and Convolutional Neural Network (CNN) are constructed to predict photo-$z$ from fluxes and images, respectively. We also propose to use an efficient hybrid network, which combines MLP and CNN, by employing transfer learning techniques to investigate the improvement of the result with both flux and image data included. We find that the photo-$z$ accuracy and outlier fraction can achieve $σ_{\rm NMAD} = 0.023$ and $η= 1.43\%$ for the MLP using flux data only, and $σ_{\rm NMAD} = 0.025$ and $η= 1.21\%$ for the CNN using image data only. The result can be further improved in high efficiency as $σ_{\rm NMAD} = 0.020$ and $η= 0.90\%$ for the hybrid transfer network. These approaches result in similar galaxy median and mean redshifts ~0.8 and 0.9, respectively, for the redshift range from 0 to 4. This indicates that our networks can effectively and properly extract photo-$z$ information from the CSST galaxy flux and image data.

preprint2022arXiv

Forecast of Neutrino Cosmology from the CSST Photometric Galaxy Clustering and Cosmic Shear Surveys

China Space Station Telescope (CSST) is a forthcoming powerful Stage IV space-based optical survey equipment. It is expected to explore a number of important cosmological problems in extremely high precision. In this work, we focus on investigating the constraints on neutrino mass and other cosmological parameters under the model of cold dark matter with a constant equation of state of dark energy ($w$CDM), using the mock data from the CSST photometric galaxy clustering and cosmic shear surveys (i.e. 3$\times$2pt). The systematics from galaxy bias, photometric redshift uncertainties, intrinsic alignment, shear calibration, baryonic feedback, non-linear, and instrumental effects are also included in the analysis. We generate the mock data based on the COSMOS catalog considering the instrumental and observational effects of the CSST, and make use of the Markov Chain Monte Carlo (MCMC) method to perform the constraints. Comparing to the results from current similar measurements, we find that CSST 3$\times$2pt surveys can improve the constraints on the cosmological parameters by one order of magnitude at least. We can obtain an upper limit for the sum of neutrino mass $Σm_ν \lesssim 0.36$ (0.56) eV at 68\% (95\%) confidence level, and $Σm_ν \lesssim 0.23$ (0.29) eV at 68\% (95\%) confidence level if ignore the baryonic effect, which is comparable to the {\it Planck} results and much better than the current photometric surveys. This indicates that the CSST photometric surveys can provide stringent constraints on the neutrino mass and other cosmological parameters, and the results also can be further improved by including data from other kinds of CSST cosmological surveys.

preprint2022arXiv

Lunar Orbit Measurement of Cosmic Dawn 21 cm Global Spectrum

A redshifted 21 cm line absorption signature is commonly expected from the cosmic dawn era, when the first stars and galaxies formed. The detailed traits of this signal can provide important insight on the cosmic history. However, high precision measurement of this signal is hampered by the ionosphere refraction and absorption, as well as radio frequency interference (RFI). A space observation can solve the problem of the ionosphere, and the Moon can shield the RFI from the Earth. In this paper, we present simulations of the global spectrum measurement in the 30 -- 120 MHz frequency band on the lunar orbit, from the proposed Discovering the Sky at the Longest wavelength (DSL) project. In particular, we consider how the measured signal varies as the satellite moves along the orbit, take into account the blockage of different parts of the sky by the Moon and the antenna response. We estimate the sensitivity for such a 21 cm global spectrum experiment. An RMS noise level of $\le 0.05$ K is expected at 75 MHz after 10 orbits ($\sim$ 1 day) observation, for a frequency channel width of 0.4 MHz. We also study the influence of a frequency-dependent beam, which may generate complex spectral structures in the spectrum. Estimates of the uncertainties in the foreground and 21 cm model parameters are obtained.

preprint2022arXiv

The Electromagnetic Characteristics of the Tianlai Cylindrical Pathfinder Array

A great challenge for 21 cm intensity mapping experiments is the strong foreground radiation which is orders of magnitude brighter than the 21cm signal. Removal of the foreground takes advantage of the fact that its frequency spectrum is smooth while the redshifted 21cm signal spectrum is stochastic. However, a complication is the non-smoothness of the instrument response. This paper describes the electromagnetic simulation of the Tianlai cylinder array, a pathfinder for 21 cm intensity mapping experiments. Due to the vast scales involved, a direct simulation requires large amount of computing resources. We have made the simulation practical by using a combination of methods: first simulate a single feed, then an array of feed units, finally with the feed array and a cylindrical reflector together, to obtain the response for a single cylinder. We studied its radiation pattern, bandpass response and the effects of mutual coupling between feed units, and compared the results with observation. Many features seen in the measurement result are well reproduced in the simulation, especially the oscillatory features which are associated with the standing waves on the reflector. The mutual coupling between feed units is quantified with S-parameters, which decrease as the distance between the two feeds increases. Based on the simulated S-parameters, we estimate the correlated noise which has been seen in the visibility data, the results show very good agreement with the data in both magnitude and frequency structures. These results provide useful insights on the problem of 21cm signal extraction for real instruments.

preprint2021arXiv

AlgoSCR: An algorithm for Solar Contamination Removal from radio interferometric data

Hydrogen intensity mapping is a new field in astronomy that promises to make three-dimensional maps of the matter distribution of the Universe using the redshifted $21\,\textrm{cm}$ line of neutral hydrogen gas (HI). Several ongoing and upcoming radio interferometers, such as Tianlai, CHIME, HERA, HIRAX, etc. are using this technique. These instruments are designed to map large swaths of the sky by drift scanning over periods of many months. One of the challenges of the observations is that the daytime data is contaminated by strong radio signals from the Sun. In the case of Tianlai, this results in almost half of the measured data being unusable. We try to address this issue by developing an algorithm for solar contamination removal (AlgoSCR) from the radio data. The algorithm is based on an eigenvalue analysis of the visibility matrix, and hence is applicable only to interferometers. We apply AlgoSCR to simulated visibilities, as well as real daytime data from the Tianlai dish array. The algorithm can remove most of the solar contamination without seriously affecting other sky signals and thus makes the data usable for certain applications.

preprint2021arXiv

Spectroscopic and Photometric Redshift Estimation by Neural Networks For the China Space Station Optical Survey (CSS-OS)

The estimation of spectroscopic and photometric redshifts (spec-z and photo-z) is crucial for future cosmological surveys. It can directly affect several powerful measurements of the Universe, e.g. weak lensing and galaxy clustering. In this work, we explore the accuracies of spec-z and photo-z that can be obtained in the China Space Station Optical Surveys (CSS-OS), which is a next-generation space survey, using neural networks. The 1-dimensional Convolutional Neural Networks (1-d CNN) and Multi-Layer Perceptron (MLP, one of the simplest forms of Artificial Neural Network) are employed to derive the spec-z and photo-z, respectively. The mock spectral and photometric data used for training and testing the networks are generated based on the COSMOS catalog. The networks have been trained with noisy data by creating Gaussian random realizations to reduce the statistical effects, resulting in similar redshift accuracy for both high-SNR (signal to noise ratio) and low-SNR data. The probability distribution functions (PDFs) of the predicted redshifts are also derived via Gaussian random realizations of the testing data, and then the best-fit redshifts and 1-sigma errors also can be obtained. We find that our networks can provide excellent redshift estimates with accuracies ~0.001 and 0.01 on spec-z and photo-z, respectively. Compared to existing photo-z codes, our MLP has similar accuracy but is more efficient in the training process. The fractions of catastrophic redshifts or outliers can be dramatically suppressed comparing to the ordinary template-fitting method. This indicates that the neural network method is feasible and powerful for spec-z and photo-z estimations in future cosmological surveys.

preprint2021arXiv

Underwater Image Enhancement based on Deep Learning and Image Formation Model

Underwater robots play an important role in oceanic geological exploration, resource exploitation, ecological research, and other fields. However, the visual perception of underwater robots is affected by various environmental factors. The main challenge now is that images captured by underwater robots are color-distorted. The hue of underwater images tends to be close to green and blue. In addition, the contrast is low and the details are fuzzy. In this paper, a new underwater image enhancement algorithm based on deep learning and image formation model is proposed. Experimental results show that the advantages of the proposed method are that it eliminates the influence of underwater environmental factors, enriches the color, enhances details, achieves higher scores in PSNR and SSIM metrics, and helps feature key-point point matching get better results. Another significant advantage is that its computation speed is much faster than other methods.

preprint2020arXiv

A Radio Burst Detection Method Based on the Hough Transform

We present a simple and fast method for incoherent dedispersion and fast radio burst (FRB) detection based on the Hough transform, which is widely used for feature extraction in image analysis. The Hough transform maps a point in the time-frequency data to a straight line in the parameter space, and points on the same dispersed $f^{-2}$ curve to a bundle of lines all crossing at the same point, thus the curve is transformed to a single point in the parameter space, enabling an easier way for the detection of radio burst. By choosing an appropriate truncation threshold, in a reasonably radio quiet environment, i.e. with radio frequency interferences (RFIs) present but not dominant, the computing speed of the method is very fast. Using simulation data of different noise levels, we studied how the detected peak varies with different truncation thresholds. We also tested the method with some real pulsar and FRB data.

preprint2020arXiv

Cosmological constraints from line intensity mapping with interlopers

Understanding the formation and evolution of the Universe is crucial for cosmological studies, and the line intensity mapping provides a powerful tool for this kind of study. We propose to make use of multipole moments of redshift-space line intensity power spectrum to constrain the cosmological and astrophysical parameters, such as the equation of state of dark energy, massive neutrinos, primordial non-Gaussianity, and star formation rate density. As an example, we generate mock data of multipole power spectra for H-alpha 6563AA, [OIII] 5007AA and [OII] 3727AA measured by SPHEREx experiment at z=1 considering contaminations from interloper lines, and use Markov Chain Monte Carlo (MCMC) method to constrain the parameters in the model. We find a good fitting result of the parameters compared to their fiducial values, which means that the multipole power spectrum can effectively distinguish signal and interloper lines, and break the degeneracies between parameters, such as line mean intensity and bias. We also explore the cross power spectrum with CSST (Chinese Space Station Telescope) spectroscopic galaxy survey in the constraints. Since more accurate fitting results can be obtained by including measurements of the emission lines at higher redshifts out to z=3 at least and cross-correlations between emission lines can be involved, the line intensity mapping is expected to offer excellent results in future cosmological and astrophysical studies.

preprint2020arXiv

Forecast for FAST: from Galaxies Survey to Intensity Mapping

The Five-Hundred-Meter Aperture Spherical Radio Telescope(FAST) is the largest single-dish radio telescope in the world. In this paper, we make forecast on the FAST HI large scale structure survey by mock observations. We consider a drift scan survey with the L-band 19 beam receiver, which may be commensal with the pulsar search and Galactic HI survey. We also consider surveys at lower frequency, either using the current single feed wide band receiver, or a future multi-beam phased array feed (PAF) in the UHF band. We estimate the number density of detected HI galaxies and the measurement error in positions, the precision of the surveys are evaluated using both Fisher matrix and simulated observations. The measurement error in the HI galaxy power spectrum is estimated, and we find that the error is relatively large even at moderate redshifts, as the number of positively detected galaxies drops drastically with increasing redshift. However, good cosmological measurement could be obtained with the intensity mapping technique where the large scale HI distribution is measured without resolving individual galaxies. The figure of merit (FoM) for the dark energy equation of state with different observation times are estimated, we find that with the existing L-band multi-beam receiver, a good measurement of low redshift large scale structure can be obtained, which complements the existing optical surveys. With a PAF in the UHF band, the constraint can be much stronger, reaching the level of a dark energy task force (DETF) stage IV experiment.

preprint2020arXiv

On the magnetospheric origin of Repeating Fast Radio Bursts

A bright radio burst was newly discovered in SGR 1935+2154, which exhibit some FRB-like temporal- and frequency-properties, suggesting a neutron star (NS)/magnetar magnetospheric origin of FRBs. We propose an explanation of the temporal- and frequency-properties of sub-pulses of repeating FRBs based on the generic geometry within the framework of charged-bunching coherent curvature radiation in the magnetosphere of an NS. The sub-pulses in a radio burst come from bunches of charged particles moving along different magnetic field lines. Their radiation beam sweep across the line of sight at slightly different time, and those radiating at the more curved part tend to be seen earlier and at higher frequency. However, by considering bunches generated at slightly different times, we find there is also a small probability that the emission from the less curved part be seen earlier. We simulate the time--frequency structures by deriving various forms of the electric acceleration field in the magnetosphere. Such structure of sub-pulses is a natural consequence of coherent curvature radiation from an NS magnetosphere with suddenly and violently triggered sparks. We apply this model to explain the time--frequency structure within specific dipolar configuration by invoking the transient pulsar-like sparking from the inner gap of a slowly rotating NS, and have also developed in more generic configurations.

preprint2020arXiv

Testing interacting dark matter and dark energy model with cosmological data

We investigate the model of dark matter-dark energy (DM-DE) interaction with coupling strength proportional to the multiplication of dark sector densities with different power indices $Q = γρ_{\rm c}^α ρ_{\rm d}^β$. We first investigate the modification of the cosmic expansion history, and then further develop the formalism to take into account the cosmological perturbations and dark matter temperature evolution. We then use the latest observational cosmology data, including cosmic microwave background (CMB) data, baryon acoustic oscillations (BAO) data, redshift-space distortion (RSD) data and Type Ia supernovae (SNe) data to constrain the model parameters. We find in the phantom region, a positive $α$ is preferred by the data above $2\, σ$ statistic significance. If we choose the power indices to be integers or half-integers for {\it plausible} physics of particle interaction, the allowed values within $1\, σ$ confidence regions are $α= 0.5$ and $β= 0, 0.5, 1$. The inclusion of BAO and RSD data from large-scale structure and SNe data improves the constraints significantly. Our model predicts lower values of $f(z) σ_8(z)$ at $z<1$ comparing to $Λ$CDM model, which alleviates the tension of $Λ$CDM with various RSD data from optical galaxy surveys. Overall, the DM-DE interaction model is consistent with the current observational data, especially providing a better fit to the RSD data.

preprint2020arXiv

The cosmic atomic hydrogen mass density as a function of mass and galaxy hierarchy from spectral stacking

We use spectral stacking to measure the contribution of galaxies of different masses and in different hierarchies to the cosmic atomic hydrogen (HI) mass density in the local Universe. Our sample includes 1793 galaxies at $z < 0.11$ observed with the Westerbork Synthesis Radio Telescope, for which Sloan Digital Sky Survey spectroscopy and hierarchy information are also available. We find a cosmic HI mass density of $Ω_{\rm HI} = (3.99 \pm 0.54)\times 10^{-4} h_{70}^{-1}$ at $\langle z\rangle = 0.065$. For the central and satellite galaxies, we obtain $Ω_{\rm HI}$ of $(3.51 \pm 0.49)\times 10^{-4} h_{70}^{-1}$ and $(0.90 \pm 0.16)\times 10^{-4} h_{70}^{-1}$, respectively. We show that galaxies above and below stellar masses of $\sim$10$^{9.3}$ M$_{\odot}$ contribute in roughly equal measure to the global value of $Ω_{\rm HI}$. While consistent with estimates based on targeted HI surveys, our results are in tension with previous theoretical work. We show that these differences are, at least partly, due to the empirical recipe used to set the partition between atomic and molecular hydrogen in semi-analytical models. Moreover, comparing our measurements with the cosmological semi-analytic models of galaxy formation {\sc Shark} and GALFORM reveals gradual stripping of gas via ram pressure works better to fully reproduce the properties of satellite galaxies in our sample, than strangulation. Our findings highlight the power of this approach in constraining theoretical models, and confirm the non-negligible contribution of massive galaxies to the HI mass budget of the local Universe.

preprint2019arXiv

Cosmological constraints on ultra-light axion fields

Ultra-light axions (ULAs) with mass less than 10^-20 eV have interesting behaviors that may contribute to either dark energy or dark matter at different epochs of the Universe. Its properties can be explored by cosmological observations, such as expansion history of the Universe, cosmic large-scale structure, cosmic microwave background, etc. In this work, we study the ULAs with a mass around 10^-33 eV, which means the ULA field still rolls slowly at present with the equation of state w=-1 as dark energy. In order to investigate the mass and other properties of this kind of ULA field, we adopt the measurements of Type Ia supernova (SN Ia), baryon acoustic oscillation (BAO), and Hubble parameter H(z). The Markov Chain Monte Carlo (MCMC) technique is employed to perform the constraints on the parameters. Finally, by exploring four cases of the model, we find that the mass of this ULA field is about 3x10^-33 eV if assuming the initial axion field phi_i=M_pl. We also investigate a general case by assuming phi_i< M_pl and find that the fitting results of phi_i/M_pl are consistent with or close to 1 for the datasets we use.

preprint2019arXiv

FRB 171019: An event of binary neutron star merger?

The fast radio burst, FRB 171019, was relatively bright when discovered first by ASKAP, but was identified as a repeater with three faint bursts detected later by GBT and CHIME. These observations lead to the discussion of whether the first bright burst shares the same mechanism with the following repeating bursts. A model of binary neutron star merger is proposed for FRB 171019, in which the first bright burst occurred during the merger event, while the subsequent repeating bursts are starquake-induced, and generally fainter, as the energy release rate for the starquakes can hardly exceed that of the catastrophic merger event. This scenario is consistent with the observation that no burst detected is as bright as the first one.

preprint2018arXiv

Cosmology with Phase 1 of the Square Kilometre Array; Red Book 2018: Technical specifications and performance forecasts

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

preprint2016arXiv

Axion decay and anisotropy of near-IR extragalactic background light

The extragalactic background light (EBL) is comprised of the cumulative radiation from all galaxies and active galactic nuclei over the cosmic history. In addition to point sources, EBL also contains information from diffuse sources of radiation. The angular power spectra of the near-infrared intensities could contain additional signals and a complete understanding of the nature of the IR background is still lacking in the literature. Here we explore the constraints that can be placed on particle decays, especially candidate dark matter models involving axions that trace dark matter halos of galaxies. Axions with a mass around a few eV will decay via two photons with wavelengths in the near-IR band, and will leave a signature in the IR background intensity power spectrum. Using recent power spectra measurements from the Hubble Space Telescope (HST) and Cosmic Infrared Background Experiment (CIBER), we find that the 0.6 to 1.6 micron power spectra can be explained by axions with masses around 4 eV. The total axion abundance Omega_a~0.05, and it is comparable to the baryon density of the Universe. The suggested mean axion mass and abundance are not ruled out by existing cosmological observations. Interestingly, the axion model with a mass distribution is preferred by the data, which cannot be explained by the standard quantum chromodynamics (QCD) theory and needs further discussion.

preprint2016arXiv

Measurement of Neutrino Masses from Relative Velocities

We present a new technique to measure neutrino masses using their flow field relative to dark matter. Present day streaming motions of neutrinos relative to dark matter and baryons are several hundred km/s, comparable with their thermal velocity dispersion. This results in a unique dipole anisotropic distortion of the matter-neutrino cross power spectrum, which is observable through the dipole distortion in the cross correlation of different galaxy populations. Such a dipole vanishes if not for this relative velocity and so it is a clean signature for neutrino mass. We estimate the size of this effect and find that current and future galaxy surveys may be sensitive to these signature distortions.

preprint2016arXiv

Precise measurements of inflationary features with 21 cm observations

Future observations of 21~cm emission using HI intensity mapping will enable us to probe the large scale structure of the Universe over very large survey volumes within a reasonable observation time. We demonstrate that the three-dimensional information contained in such surveys will be an extremely powerful tool in searching for features that were imprinted in the primordial power spectrum and bispectrum during inflation. Here we focus on the "resonant" and "step" inflation models, and forecast the potential of upcoming 21~cm experiments to detect these inflationary features in the observable power- and bispectrum. We find that the full scale Tianlai experiment and the Square Kilometre Array (SKA) have the potential to improve on the sensitivity of current Cosmic Microwave Background (CMB) experiments by several orders of magnitude.

preprint2016arXiv

Primordial density and BAO reconstruction

We present a new method to reconstruct the primordial (linear) density field using the estimated nonlinear displacement field. The divergence of the displacement field gives the reconstructed density field. We solve the nonlinear displacement field in the 1D cosmology and show the reconstruction results. The new reconstruction algorithm recovers a lot of linear modes and reduces the nonlinear damping scale significantly. The successful 1D reconstruction results imply the new algorithm should also be a promising technique in the 3D case.

preprint2016arXiv

Probing Neutrino Hierarchy and Chirality via Wakes

The relic neutrinos are expected to acquire a bulk relative velocity with respect to the dark matter at low redshifts, and neutrino wakes are expected to develop downstream of the dark matter halos. We propose a method of measuring the neutrino mass based on this mechanism. This neutrino wake will cause a dipole distortion of the galaxy-galaxy lensing pattern. This effect could be detected by combining upcoming lensing surveys with a low redshift galaxy survey or a 21 cm intensity mapping survey, which can map the neutrino flow field. The data obtained with LSST and Euclid should enable us to make a positive detection if the three neutrino masses are quasidegenerate with each neutrino mass of $\sim$0.1 eV, and a future high precision 21 cm lensing survey would allow the normal hierarchy and inverted hierarchy cases to be distinguished, and even the right-handed Dirac neutrinos may be detectable.

preprint2016arXiv

Sky reconstruction for the Tianlai cylinder array

In this paper, we apply our sky map reconstruction method for transit type interferometers to the Tianlai cylinder array. The method is based on the spherical harmonic decomposition, and can be applied to cylindrical array as well as dish arrays and we can compute the instrument response, synthesised beam, transfer function and the noise power spectrum. We consider cylinder arrays with feed spacing larger than half wavelength, and as expected, we find that the arrays with regular spacing have grating lobes which produce spurious images in the reconstructed maps. We show that this problem can be overcome, using arrays with different feed spacing on each cylinder. We present the reconstructed maps, and study the performance in terms of noise power spectrum, transfer function and beams for both regular and irregular feed spacing configurations.

preprint2016arXiv

Sky reconstruction from transit visibilities: PAON-4 and Tianlai Dish Array

The spherical harmonics $m$-mode decomposition is a powerful sky map reconstruction method suitable for radio interferometers operating in transit mode. It can be applied to various configurations, including dish arrays and cylinders. We describe the computation of the instrument response function, the point spread function (PSF), transfer function, the noise covariance matrix and noise power spectrum. The analysis in this paper is focused on dish arrays operating in transit mode. We show that arrays with regular spacing have more pronounced side lobes as well as structures in their noise power spectrum, compared to arrays with irregular spacing, specially in the north-south direction. A good knowledge of the noise power spectrum $C^{\mathrm{noise}}(\ell)$ is essential for intensity mapping experiments as non uniform $C^{\mathrm{noise}}(\ell)$ is a potential problem for the measurement of the HI power spectrum. Different configurations have been studied to optimise the PAON-4 and Tianlai dish array layouts. We present their expected performance and their sensitivities to the 21-cm emission of the Milky Way and local extragalactic HI clumps

preprint2016arXiv

Two-dimensional Topology of Cosmological Reionization

We study the two-dimensional topology of the 21-cm differential brightness temperature for two hydrodynamic radiative transfer simulations and two semi-numerical models. In each model, we calculate the two dimensional genus curve for the early, middle and late epochs of reionization. It is found that the genus curve depends strongly on the ionized fraction of hydrogen in each model. The genus curves are significantly different for different reionization scenarios even when the ionized faction is the same. We find that the two-dimensional topology analysis method is a useful tool to constrain the reionization models. Our method can be applied to the future observations such as those of the Square Kilometer Array.

preprint2015arXiv

Cosmology from the EoR/Cosmic Dawn with the SKA

The SKA will build upon early detections of the EoR by precursor instruments, such as MWA, PAPER, and LOFAR, and planned instruments, such as HERA, to make the first high signal-to-noise measurements of fluctuations in the 21 cm brightness temperature from both reionization and the cosmic dawn. This will allow both imaging and statistical maps of the 21cm signal at redshifts z = 6 - 27 and constrain the underlying cosmology and evolution of the density field. This era includes nearly 60% of the (in principle) observable volume of the Universe and many more linear modes than the CMB, presenting an opportunity for SKA to usher in a new level of precision cosmology. This optimistic picture is complicated by the need to understand and remove the effect of astrophysics, so that systematics rather than statistics will limit constraints. This chapter describes the cosmological, as opposed to astrophysical, information available to SKA. Key areas for discussion include: cosmological parameters constraints using 21cm fluctuations as a tracer of the density field; lensing of the 21cm signal, constraints on heating via exotic physics such as decaying or annihilating dark matter; impact of fundamental physics such as non-Gaussianity or warm dark matter on the source population; and constraints on the bulk flows arising from the decoupling of baryons and photons at z = 1000. The chapter explores the path to separating cosmology from astrophysics, for example via velocity space distortions and separation in redshift. We discuss new opportunities for extracting cosmology made possible by the sensitivity of SKA Phase 1 and explores the advances achievable with SKA2.

preprint2015arXiv

Dense magnetized plasma associated with a fast radio burst

Fast Radio Bursts are bright, unresolved, non-repeating, broadband, millisecond flashes, found primarily at high Galactic latitudes, with dispersion measures much larger than expected for a Galactic source. The inferred all-sky burst rate is comparable to the core-collapse supernova rate out to redshift 0.5. If the observed dispersion measures are assumed to be dominated by the intergalactic medium, the sources are at cosmological distances with redshifts of 0.2 to 1. These parameters are consistent with a wide range of source models. One fast radio burst showed circular polarization [21(7)%] of the radio emission, but no linear polarization was detected, and hence no Faraday rotation measure could be determined. Here we report the examination of archival data revealing Faraday rotation in a newly detected burst - FRB 110523. It has radio flux at least 0.6 Jy and dispersion measure 623.30(5) pc cm$^{-3}$. Using Galactic contribution 45 pc cm$^{-3}$ and a model of intergalactic electron density, we place the source at a maximum redshift of 0.5. The burst has rotation measure -186.1(1.4) rad m$^{-2}$, much higher than expected for this line of sight through the Milky Way and the intergalactic medium, indicating magnetization in the vicinity of the source itself or within a host galaxy. The pulse was scattered by two distinct plasma screens during propagation, which requires either a dense nebula associated with the source or a location within the central region of its host galaxy. Keeping in mind that there may be more than one type of fast radio burst source, the detection in this instance of source-local magnetization and scattering favours models involving young stellar populations such as magnetars over models involving the mergers of older neutron stars, which are more likely to be located in low density regions of the host galaxy.

preprint2015arXiv

Probing First Galaxies and Their Impact on the Intergalactic Medium through the 21-cm Observation of the Cosmic Dawn with the SKA

We present an overview of the theory of high-redshift star and X-ray source formation, and how they affect the 21-cm background. Primary focus is given to Lyman alpha pumping and X-ray heating mechanisms at cosmic dawn, opening a new observational window for high-redshift astrophysics by generating sizable fluctuations in the 21-cm background. We describe observational prospects for power spectrum analysis and 3D tomography (imaging) of the signature of these early astrophysical sources by SKA1-LOW and SKA2.

preprint2015arXiv

Quiescent luminous red galaxies (LRGs) as cosmic chronometers: on the significance of mass and environmental dependence

Massive luminous red galaxies (LRGs) are believed to be evolving passively and can be used as cosmic chronometers to estimate the Hubble constant. However, different LRGs may be located in different environments. The environmental effects may limit the use of the LRGs as cosmic chronometers. We aim to investigate the environmental and mass dependence of the formation of 'quiescent' LRGs selected from the Sloan Digital Sky Survey Date Release 8 and to pave the way for using LRGs as cosmic chronometers. Using the population synthesis software STARLIGHT, we derive the stellar populations in each LRG through the full spectrum fitting and obtain the mean age distribution and the mean star formation history (SFH) of those LRGs. We find that there is no apparent dependence of the mean age and the SFH of quiescent LRGs on their environment, while the ages of those quiescent LRGs depend weakly on their mass. We compare the SFHs of the SDSS LRGs with those obtained from a semi-analytical galaxy formation model, and find that they are roughly consistent with each other if we consider the errors in the STARLIGHT-derived ages. We find that a small fraction of later star formation in LRGs leads to a systematical overestimation (~28 %) of the Hubble constant by the differential age method, and the systematical errors in the STARLIGHT-derived ages may lead to an underestimation (~ 16 %) of the Hubble constant. However, these errors can be corrected by a detailed study of the mean SFH of those LRGs and by calibrating the STARLIGHT-derived ages to those obtained independently by other methods. The environmental effects do not play significant role in the age estimates of quiescent LRGs, and the quiescent LRGs as a population can be used securely as cosmic chronometers.

preprint2015arXiv

Quintessence and phantom emerging from the split-complex field and the split-quaternion field

Motivated by the mathematic theory of split-complex numbers (or hyperbolic numbers, also perplex numbers) and the split-quaternion numbers (or coquaternion numbers), we define the notion of split-complex scalar field and the split-quaternion scalar field. Then we explore the cosmic evolution of these scalar fields in the background of spatially flat Friedmann-Robertson-Walker Universe. We find that both the quintessence field and the phantom field could naturally emerge in these scalar fields. Introducing the metric of field space, these theories fall into a subclass of the multi-field theories which have been extensively studied in inflationary cosmology.

preprint2015arXiv

The consistency test on the cosmic evolution

We propose a new and robust method to test the consistency of the cosmic evolution given by a cosmological model. It is realized by comparing the combined quantity r_d^CMB/D_V^SN, which is derived from the comoving sound horizon r_d from cosmic microwave background (CMB) measurements and the effective distance D_V derived from low-redshift Type-Ia supernovae (SNe Ia) data, with direct and independent r_d/D_V obtained by baryon acoustic oscillation (BAO) measurements at median redshifts. We apply this test method for the Lambda-CDM and wCDM models, and investigate the consistency of the derived value of r_d/D_V from Planck 2015 and the SN Ia data sets of Union2.1 and JLA (z<1.5), and the r_d/D_V directly given by BAO data from six-degree-field galaxy survey (6dFGS), Sloan Digital Sky Survey Data Release 7 Main Galaxy Survey (SDSS-DR7 MGS), DR11 of SDSS-III, WiggleZ and Ly-alpha forecast surveys from Baryon Oscillation Spectroscopic Data (BOSS) DR-11 over 0.1<z<2.36. We find that r_d^CMB/D_V^SN for both non-flat Lambda-CDM and flat wCDM models with Union2.1 and JLA data are well consistent with the BAO and CMB measurements within 1-sigma CL. Future surveys will further tight up the constraints significantly, and provide stronger test on the consistency.

preprint2015arXiv

Topology of neutral hydrogen distribution with the Square Kilometer Array

Morphology of the complex HI gas distribution can be quantified by statistics like the Minkowski functionals, and can provide a way to statistically study the large scale structure in the HI maps both at low redshifts, and during the epoch of reionization (EoR). At low redshifts, the 21cm emission traces the underlying matter distribution. Topology of the HI gas distribution, as measured by the genus, could be used as a "standard ruler". This enables the determination of distance-redshift relation and also the discrimination of various models of dark energy and of modified gravity. The topological analysis is also sensitive to certain primordial non-Gaussian features. Compared with two-point statistics, the topological statistics are more robust against the nonlinear gravitational evolution, bias, and redshift-space distortion. The HI intensity map observation naturally avoids the sparse sampling distortion, which is an important systematic in optical galaxy survey. The large cosmic volume accessible to SKA would provide unprecedented accuracy using such a measurement... [abridged]

preprint2014arXiv

Cosmological parameter fittings with the BICEP2 data

Combining the latest Planck, Wilkinson Microwave Anisotropy Probe (WMAP), and baryon acoustic oscillation (BAO) data, we exploit the recent cosmic microwave background (CMB) B-mode power spectra data released by the BICEP2 collaboration to constrain the cosmological parameters of the $Λ$CDM model, esp. the primordial power spectra parameters of the scalar and the tensor modes, $n_s,α_s, r, n_t$. We obtained constraints on the parameters for a lensed $Λ$CDM model using the Markov Chain Monte Carlo (MCMC) technique, the marginalized $68\%$ bounds are $r = 0.1043\pm_{0.0914}^{0.0307}$, $n_s = 0.9617\pm_{0.0061}^{ 0.0061}$, $α_s = -0.0175\pm_{0.0097}^{0.0105}$, $n_t = 0.5198\pm_{0.4579}^{0.4515}$. We found that a blue tilt for $n_t$ is favored slightly, but it is still well consistent with flat or even red tilt. Our $r$ value is slightly smaller than the one obtained by the BICEP team, as we allow $n_t$ as a free parameter without imposing the single-field slow roll inflation consistency relation. If we impose this relation, $r=0.2130\pm_{0.0609}^{0.0446}$. For most other parameters, the best fit values and measurement errors are not altered much by the introduction of the BICEP2 data.

preprint2014arXiv

Forecasts on the Dark Energy and Primordial Non-Gaussianity Observations with the Tianlai Cylinder Array

The Tianlai experiment is dedicated to the observation of large scale structures (LSS) by the 21 cm intensity mapping technique. In this paper we make forecasts on its capability at observing or constraining the dark energy parameters and the primordial non-Gaussianity. From the LSS data one can use the baryon acoustic oscillation (BAO) and the growth rate derived from the redshift space distortion (RSD) to measure the dark energy density and equation of state. The primordial non-Gaussianity can be constrained either by looking for scale-dependent bias in the power spectrum, or by using the bispectrum. Here we consider three cases: the Tianlai cylinder array pathfinder which is currently being built, an upgrade of the pathfinder array with more receiver units, and the full-scale Tianlai cylinder array. Using the full-scale Tianlai experiment, we expect $σ_{w_0} \sim 0.082$ and $σ_{w_a} \sim 0.21$ from the BAO and RSD measurements, $σ_{\rm f_{NL}}^{\rm local} \sim 14$ from the power spectrum measurements with scale-dependent bias, and $σ_{\rm f_{NL}}^{\rm local} \sim 22$ and $σ_{\rm f_{NL}}^{\rm equil} \sim 157$ from the bispectrum measurements.

preprint2014arXiv

Numerical Simulation of Radio Signal from Extended Air Showers

The burst of radio emission by the extensive air shower provides a promising alternative for detecting ultra-high energy cosmic rays.We have developed an independent numerical program to simulate these radio signals. Our code is based on a microscopic treatment, with both the geosynchrotron radiation and charge excess effect included. Here we make a first presentation of our basic program and its results. The time signal for different polarizations are computed, we find that the pulses take on a bipolar pattern, the spectrum is suppressed towards the lower frequencies.We investigate how the shower at different heights in atmosphere contribute to the total signal, and examine the signal strength and distribution at sites of different elevations. We also study the signal from showers of different inclination angles and azimuth directions. In all these cases we find the charge excess effect important.

preprint2014arXiv

TeV cosmic-ray proton and helium spectra in the myriad model II

Recent observations show that the cosmic ray nuclei spectra start to harden above 100 GeV, in contradiction with the conventional steady-state cosmic ray model. We had suggested that this anomaly is due to the propagation effect of cosmic rays released from local young cosmic ray sources, the total flux of the cosmic ray should be computed with the myriad model, where contribution from sources in local catalog is added to the background. However, while the hardening could be elegantly explained in this model, the model parameters obtained from the fit skew toward a region with fast diffusion and low supernova rate in the Galaxy, in tension with other observations. In this paper, we further explore this model in order to set up a concordant picture. Two possible improvements related to the cosmic ray sources have been considered. Firstly, instead of the usual axisymmetric disk model, we considered a spiral model of source distribution. Secondly, for the nearby and young sources which are paramount to the hardening, we allow for an energy-dependent escape time. We find that major improvement comes from the energy-dependent escape time of the local sources, and with both modifications, not only the cosmic ray proton and helium anomalies are solved, but also the parameters attain reasonable range values compatible with other analysis.

preprint2014arXiv

The brief era of direct collapse black hole formation

It has been proposed that the first, intermediate-mass ($\approx 10^{5-6}~M_\odot$) black holes might form through direct collapse of unpolluted gas in atomic-cooling halos exposed to a strong Lyman-Werner (LW) or near-infrared (NIR) radiation. As these systems are expected to be Compton-thick, photons above 13.6 eV are largely absorbed and re-processed into lower energy bands. It follows that direct collapse black holes (DCBHs) are very bright in the LW/NIR bands, typically outshining small high-redshift galaxies by more than 10 times. Once the first DCBHs form, they then trigger a runaway process of further DCBH formation, producing a sudden rise in their cosmic mass density. The universe enters the "DCBH era" at $z \approx 20$ when a large fraction of atomic-cooling halos are experiencing DCBH formation. By combining the clustering properties of the radiation sources with Monte Carlo simulations we show that in this scenario the DCBH mass density rises from $\sim 5$~$M_\odot$ Mpc$^{-3}$ at $z\sim 30$ to the peak value $\sim5\times10^5 M_\odot$ Mpc$^{-3}$ at $z \sim 14$ in our fiducial model. However, the abundance of \textit{active} (accreting) DCBHs drops after $z \sim 14$, as gas in the potential formation sites (unpolluted halos with virial temperature slightly above $10^4$~K) is photoevaporated. This effect almost completely suppresses DCBH formation after $z\sim 13$. The DCBH formation era lasts only $\approx 150$ Myr, but it might crucially provide the seeds of the supermassive black holes (SMBHs) powering $z\sim6$ quasars.

preprint2013arXiv

An analytical model of the large neutral regions during the late stage of reionization

In this paper we investigate the nature and distribution of large neutral regions during the late epoch of reionization. In the "bubble model" of reionization, the mass distribution of large ionized regions ("bubbles") during the early stage of reionization is obtained by using the excursion set model, where the ionization of a region corresponds to the first up-crossing of a barrier by random trajectories. We generalize this idea, and develop a method to predict the distribution of large scale neutral regions during the late stage of reionization, taking into account the ionizing background after the percolation of HII regions. The large scale neutral regions which we call "neutral islands" are not individual galaxies or minihalos, but larger regions where fewer galaxies formed and hence ionized later, and they are identified in the excursion set model with the first down-crossings of the island barrier. Assuming that the consumption rate of ionizing background photons is proportional to the surface area of the neutral islands, we obtained the size distribution of the neutral islands. We also take the "bubbles-in-island" effect into account by considering the conditional probability of up-crossing a bubble barrier after down-crossing the island barrier. We find that this effect is very important. An additional barrier is set to avoid islands being percolated through. We find that there is a characteristic scale for the neutral islands, while the small islands are rapidly swallowed up by the ionizing background, this characteristic scale does not change much as the reionization proceeds.

preprint2013arXiv

Can Large Extra Dimensions Solve the Proton Radius Puzzle?

The proton charge radius extracted from the recent muonic hydrogen spectroscopy [Antognini et al. 2013; Pohl et al. 2010] differs from the CODATA 2010 recommended value [Mohr et al. 2012] by more than 4% or $4.4 σ$. This discrepancy, dubbed as the "Proton Radius Puzzle", is a big challenge to the Standard Model of particle physics, and has triggered a number of works on the quantum electrodynamic calculations recently. The proton radius puzzle may indicate the presence of an extra correction which enlarges the 2S-2P energy gap in muonic hydrogen. Here we explore the possibility of large extra dimensions which could modify the Newtonian gravity at small scales and lower the 2S state energy while leaving the 2P state nearly unchanged. We find that such effect could be produced by four or more large extra dimensions which are allowed by the current constraints from low energy physics.

preprint2013arXiv

Constraints on the Brans-Dicke gravity theory with the Planck data

Based on the new cosmic CMB temperature data from the Planck satellite, the 9 year polarization data from the WMAP, the BAO distance ratio data from the SDSS and 6dF surveys, we place a new constraint on the Brans-Dicke theory. We adopt a parametrization $ζ=\ln(1+1/ω})$, where the general relativity (GR) limit corresponds to $ζ= 0$. We find no evidence of deviation from general relativity. At 95% probability, $-0.00246 < ζ< 0.00567$, correspondingly, the region $-407.0 < ω<175.87$ is excluded. If we restrict ourselves to the $ζ>0$ (i.e. $ω>0$) case, then the 95% probability interval is $ζ<0.00549$, corresponding to $ω> 181.65$. We can also translate this result to a constraint on the variation of gravitational constant, and find the variation rate today as $\dot{G}=-1.42^{+2.48}_{-2.27} \times 10^{-13} $yr$^{-1} $ ($1σ$ error bar), the integrated change since the epoch of recombination is $δG/G = 0.0104^{+0.0186}_{-0.0067} $ ($1σ$ error bar). These limits on the variation of gravitational constant are comparable with the precision of solar system experiments.

preprint2013arXiv

Infrared background signatures of the first black holes

Angular fluctuations of the Near InfraRed Background (NIRB) intensity are observed up to scales $\simlt 1^{\ensuremath{^{\circ}}}$. Their interpretation is challenging as even after removing the contribution from detected sources, the residual signal is $>10$ times higher than expected from distant galaxies below the detection limit and first stars. We propose here a novel interpretation in which early, intermediate mass, accreting direct collapse black holes (DCBH), which are too faint to be detected individually in current surveys, could explain the observed fluctuations. We find that a population of highly obscured ($N_{\rm H}\simgt 10^{25} \rm cm^{-2}$) DCBHs formed in metal-free halos with virial temperature $10^4$ K at $z\simgt 12$, can explain the observed level $\approx 10^{-3}$ (nW m$^{-2}$ sr$^{-1})^2$ of the 3.6 and 4.5 $μ$m fluctuations on scales $>100''$. The signal on smaller scales is instead produced by undetected galaxies at low and intermediate redshifts. Albeit Compton-thick, at scales $θ> 100''$ DCBHs produce a CXB (0.5-2 keV)-NIRB ($4.5 \rm μm$) cross-correlation signal of $\simeq 10^{-11}$ erg s$^{-1}$ cm$^{-2}$ nW m$^{-2}$ sr$^{-1}$ slightly dependent on the specific value of the absorbing gas column ($N_{\rm H} \approx 10^{25} \rm cm^{-2}$) adopted and in agreement with the recent measurements by \cite{2012arXiv1210.5302C}. At smaller scales the cross-correlation is dominated by the emission of high-mass X-ray binaries (HMXB) hosted by the same low-$z$, undetected galaxies accounting for small scale NIRB fluctuations. These results outline the great potential of the NIRB as a tool to investigate the nature of the first galaxies and black holes.

preprint2013arXiv

On the Full Spectrum Fitting of Luminous Red Galaxies by Using ULySS and STARLIGHT

In this paper, we select a sample of quiescent luminous red galaxies (LRGs) from Sloan Digital Sky Survey Data Release 7 (SDSS DR7) with signal-to-noise ratio (S/N) substantially high to study the consistency of the full spectrum fitting method by using different packages, mainly, ULySS and STARLIGHT. The spectrum of each galaxy in the sample is fitted by the full spectrum fitting packages ULySS and STARLIGHT, respectively. We find: (1) for spectra with higher S/Ns, the ages of stellar populations obtained from ULySS are slightly older than that from STARLIGHT, and metallicities derived from ULySS are slightly richer than that from STARLIGHT. In general, both packages can give roughly consistent fitting results. (2) for low S/N spectra, it is possible that the fitting by ULySS can be trapped some local minimum parametric regions during execution and thus may give unreliable results, while STARLIGHT can still give reliable results. Based on the fitting results of LRGs, we further analyse their star formation history (SFH) and the relation between their age and velocity dispersion, and find that those highly agrees with conclusions from others' previous work.

preprint2013arXiv

Relic Abundance of Asymmetric Dark Matter in Quintessence

We investigate the relic abundance of asymmetric Dark Matter particles in quintessence model with a kination phase. The analytic calculation of the asymmetric Dark Matter in the standard cosmological scenario is extended to the nonstandard cosmological scenario where we specifically discuss the quintessence model with a kination phase. We found that the enhancement of Hubble rate changes the relic density of particles and anti--particles. We use the present day Dark Matter abundance to constrain the Hubble rate in quintessence model with a kination phase for asymmetric Dark Matter.

preprint2013arXiv

The contribution of high redshift galaxies to the Near-Infrared Background

Several independent measurements have confirmed the existence of fluctuations ($δF_{\rm obs}\approx 0.1 \rm nW/m^{2}/sr$ at $3.6 \rm μm$) up to degree angular scales in the source-subtracted Near InfraRed Background (NIRB) whose origin is unknown. By combining high resolution cosmological N-body/hydrodynamical simulations with an analytical model, and by matching galaxy Luminosity Functions (LFs) and the constraints on reionization simultaneously, we predict the NIRB absolute flux and fluctuation amplitude produced by high-$z$ ($z > 5$) galaxies (some of which harboring Pop III stars, shown to provide a negligible contribution). This strategy also allows us to make an empirical determination of the evolution of ionizing photon escape fraction: we find $f_{\rm esc} = 1$ at $z \ge 11$, decreasing to $\approx 0.05$ at $z = 5$. In the wavelength range $1.0-4.5 \rm μm$, the predicted cumulative flux is $F =0.2-0.04 \rm nW/m^2/sr$. However, we find that the radiation from high-$z$ galaxies (including those undetected by current surveys) is insufficient to explain the amplitude of the observed fluctuations: at $l=2000$, the fluctuation level due to $z > 5$ galaxies is $δF = 0.01-0.002 \rm nW/m^2/sr$, with a relative wavelength-independent amplitude $δF/F = 4$%. The source of the missing power remains unknown. This might indicate that an unknown component/foreground, with a clustering signal very similar to that of high-$z$ galaxies, dominates the source-subtracted NIRB fluctuation signal.

preprint2013arXiv

The contribution of ultracompact dark matter minihalos to the isotropic radio background

The ultracompact minihalos could be formed during the earlier epoch of the universe. The dark matter annihilation within them is very strong due to the steep density profile, $ρ\sim r^{-2.25}$. The high energy electrons and positrons from the dark matter annihilation can inverse Compton scatter (ICS) with the background photons, such as CMB photons, to acquire higher energy. On the other hand, the synchrotron radiation can also be produced when they meet the magnetic field. In this paper, we study the signals from the UCMHs due to the dark matter annihilation for the radio, X-ray and $γ$-ray band. We found that for the radio emission the UCMHs can provide one kind of source for the radio excess observed by ARCADE 2. But the X-ray signals due to the ICS effect or the $γ$-ray signals mainly due to the prompt emission from dark matter would exceed the present observations, such as Fermi, COMPTEL and CHANDRA. We found that the strongest limits on the fraction of UCMHs come from the X-ray observations and the constraints from the radio data are the weakest.

preprint2012arXiv

Constraints on the dark matter annihilation scenario of Fermi 130 GeV $γ$-ray line emission by continuous gamma-rays, Milky Way halo, galaxy clusters and dwarf galaxies observations

It was recently reported that there may exist monochromatic $γ$-ray emission at $\sim 130$ GeV from the Galactic center in the Fermi Large Area Telescope data, which might be related with dark matter (DM) annihilation. In this work we carry out a comprehensive check of consistency of the results with the DM annihilation scenario, using the 3.7 yrs Fermi observation of the inner Galaxy, Galactic halo, clusters of galaxies and dwarf galaxies. The results found are as follows. 1) Very strong constraints on the DM annihilation into continuous $γ$-rays from the Galactic center are set, which are as stringent as the "natural" scale assuming thermal freeze-out of DM. Such limit sets strong constraint on the DM models to explain the line emission. 2) No line emission from the Galactic halo is found in the Fermi data, and the constraints on line emission is marginally consistent with the DM annihilation interpretation of the $\sim 130$ GeV line emission from the inner Galaxy. 3) No line emission from galaxy clusters and dwarf galaxies is detected, although possible concentration of photons from clusters in 120-140 GeV is revealed. The constraints from clusters and dwarf galaxies are weak and consistent with the DM annihilation scenario to explain the $\sim 130$ GeV line emission.

preprint2012arXiv

Cosmic Tides

We apply CMB lensing techniques to large scale structure and solve for the 3-D cosmic tidal field. We use small scale filamentary structures to solve for the large scale tidal shear and gravitational potential. By comparing this to the redshift space density field, one can measure the gravitational growth factor on large scales without cosmic variance. This potentially enables accurate measurements of neutrino masses and reconstruction of radial modes lost in 21 cm intensity mapping, which are essential for CMB and other cross correlations. We relate the tidal fields to the squeezed limit bispectrum, and present initial results from simulations and data from the SDSS.

preprint2012arXiv

Observational constraints on cosmic neutrinos and dark energy revisited

Using several cosmological observations, i.e. the cosmic microwave background anisotropies (WMAP), the weak gravitational lensing (CFHTLS), the measurements of baryon acoustic oscillations (SDSS+WiggleZ), the most recent observational Hubble parameter data, the Union2.1 compilation of type Ia supernovae, and the HST prior, we impose constraints on the sum of neutrino masses ($\mnu$), the effective number of neutrino species ($\neff$) and dark energy equation of state ($w$), individually and collectively. We find that a tight upper limit on $\mnu$ can be extracted from the full data combination, if $\neff$ and $w$ are fixed. However this upper bound is severely weakened if $\neff$ and $w$ are allowed to vary. This result naturally raises questions on the robustness of previous strict upper bounds on $\mnu$, ever reported in the literature. The best-fit values from our most generalized constraint read $\mnu=0.556^{+0.231}_{-0.288}\rm eV$, $\neff=3.839\pm0.452$, and $w=-1.058\pm0.088$ at 68% confidence level, which shows a firm lower limit on total neutrino mass, favors an extra light degree of freedom, and supports the cosmological constant model. The current weak lensing data are already helpful in constraining cosmological model parameters for fixed $w$. The dataset of Hubble parameter gains numerous advantages over supernovae when $w=-1$, particularly its illuminating power in constraining $\neff$. As long as $w$ is included as a free parameter, it is still the standardizable candles of type Ia supernovae that play the most dominant role in the parameter constraints.

preprint2012arXiv

Reionization in the Warm Dark Matter Model

Compared with the cold dark matter (CDM) model, in the warm dark matter (WDM) model formation of small scale structure is suppressed. It is often thought that this would delay the reionization of the intergalactic medium (IGM), as the star formation rate during the epoch of reionization (EOR) would be lowered. However, during the later stage of the EOR, a large portion of the ionizing photons are consumed by recombination inside the minihalos, where the gas has higher density and recombination rates than the gas in the IGM. The suppression of small scale structure would therefore reduce the recombination rate, and could potentially shorten the reionization process. This effect is investigated here by using the analytical "bubble model" of reionization. We find that in some cases, though the initiation of the EOR is delayed in the WDM model, its completion could be even earlier than the CDM case, but the effect is generally small. We obtain limits on the WDM particles mass for different reionization redshifts.

preprint2012arXiv

Semi-Numerical Simulation of Reionization with Semi-Analytical Modeling of Galaxy Formation

In a semi-numerical model of reionization, the evolution of ionization fraction is simulated approximately by the ionizing photon to baryon ratio criterion. In this paper we incorporate a semi-analytical model of galaxy formation based on the Millennium II N-body simulation into the semi-numerical modeling of reionization. The semi-analytical model is used to predict the production of ionizing photons, then we use the semi-numerical method to model the reionization process. Such an approach allows more detailed modeling of the reionization, and also connects observations of galaxies at low and high redshifts to the reionization history. The galaxy formation model we use was designed to match the low-$z$ observations, and it also fits the high redshift luminosity function reasonably well, but its prediction on the star formation falls below the observed value, and we find that it also underpredicts the stellar ionizing photon production rate, hence the reionization can not be completed at $z \sim 6$ without taking into account some other potential sources of ionization photons. We also considered simple modifications of the model with more top heavy initial mass functions (IMF), with which the reionization can occur at earlier epochs. The incorporation of the semi-analytical model may also affect the topology of the HI regions during the EoR, and the neutral regions produced by our simulations with the semi-analytical model appeared less poriferous than the simple halo-based models.

preprint2012arXiv

The Age-Redshift Relation For Luminous Red Galaxies Obtained From the Full Spectrum Fitting and Its Cosmological Implications

The relative age of galaxies at different redshifts can be used to infer the Hubble parameter and put constraints on cosmological models. We select 23,883 quiescent luminous red galaxies (LRGs) from the SDSS DR7 and divide them into four sub-samples according to their velocity dispersions and each sub-sample is further divided into 12 redshift bins. The spectra of the LRGs in each redshift and velocity bin are co-added in order to obtain a combined spectrum with relatively high $S/N$. Adopting the GalexEV/SteLib model, we estimate the mean ages of the LRGs from these combined spectra by the full-spectrum fitting method. We check the reliability of the estimated age by using Monte-Carlo simulations and find that the estimates are robust and reliable. Assuming that the LRGs in each sub-sample and each redshift bin were on average formed at the same time, the Hubble parameter at the present time $H_0$ is estimated from the age--redshift relation obtained for each sub-sample, which is compatible with the $H_0$ value measured by other methods. We demonstrate that a systematic bias (up to $\sim 20%$) may be introduced to the $H_0$ estimation because of recent star formation in the LRGs due to the later major mergers at $z\la 0.4$, but this bias may be negligible for those sub-samples with large velocity dispersions. Using the age--redshift relations obtained from the sub-sample with the largest velocity dispersion or the two sub-samples with high velocity dispersions, we find $H_0= 65^{+7}_{-3}\kmsmpc$ or $H_0= 74^{+5}_{-4}\kmsmpc$ by assuming a spatially flat $Λ$CDM cosmology. With upcoming surveys, such as the Baryon Oscillation Spectroscopic Survey (BOSS), even larger samples of quiescent massive LRGs may be obtained, and thus the Hubble parameter can be measured with high accuracy through the age--redshift relation.

preprint2012arXiv

The Tianlai project: a 21cm cosmology experiment

In my talk at the 2nd Galileo-Xu Meeting, I presented several different topics in 21cm cosmology for which I have done research. These includes the 21cm signature of the first stars[1,2], the 21cm signal from the IGM and minihalos[3], effect of dark matter annihila- tions on 21cm signal[4], the 21cm forest by ionized/neutral region[5], and the 21cm forest by minihalo and earliest galaxies[6,7]. In this conference proceeding I shall not repeat these discussions, but instead focus on the last part of my talk, i.e. the Tianlai project, an experiment effort on low redshift 21cm intensity mapping observation for dark energy measurements.

preprint2011arXiv

Black Holes in the Universe: Generalized Lemaitre-Tolman-Bondi Solutions

We present new exact solutions {which presumably describe} black holes in the background of a spatially flat, pressureless dark matter (DM)-, or dark matter plus dark energy (DM+DE)-, or quintom-dominated universe. These solutions generalize Lemaitre-Tolman-Bondi metrics. For a DM- or (DM+DE)-dominated universe, the area of the black hole apparent horizon (AH) decreases with the expansion of the universe while that of the cosmic AH increases. However, for a quintom-dominated universe, the black hole AH first shrinks and then expands, while the cosmic AH first expands and then shrinks. A (DM+DE)-dominated universe containing a black hole will evolve to the Schwarzschild-de Sitter solution with both AHs approaching constant size. In a quintom-dominated universe, the black hole and cosmic AHs will coincide at a certain time, after which the singularity becomes naked, violating Cosmic Censorship.

preprint2011arXiv

Constraint on dark matter annihilation with dark star formation using Fermi extragalactic diffuse gamma-ray background data

It has been proposed that during the formation of the first generation stars there might be a "dark star" phase in which the power of the star comes from dark matter annihilation. The adiabatic contraction process to form the dark star would result in a highly concentrated density profile of the host halo at the same time, which may give enhanced indirect detection signals of dark matter. In this work we investigate the extragalactic $γ$-ray background from dark matter annihilation with such a dark star formation scenario, and employ the isotropic $γ$-ray data from Fermi-LAT to constrain the model parameters of dark matter. The results suffer from large uncertainties of both the formation rate of the first generation stars and the subsequent evolution effects of the host halos of the dark stars. We find, in the most optimistic case for $γ$-ray production via dark matter annihilation, the expected extragalactic $γ$-ray flux will be enhanced by 1-2 orders of magnitude. In such a case, the annihilation cross section of the supersymmetric dark matter can be constrained to the thermal production level, and the leptonic dark matter model which is proposed to explain the positron/electron excesses can be well excluded. Conversely, if the positron/electron excesses are of a dark matter annihilation origin, then the early Universe environment is such that no dark star can form.

preprint2011arXiv

Constraints on ultracompact minihalos from extragalactic γ-ray background

It has been proposed that ultracompact minihalos (UCMHs) might be formed in earlier epoch. If dark matter consists of Weakly Interacting Massive Particles (WIMPs), UCMHs can be treated as the γ-ray sources due to dark matter annihilation within them. In this paper, we investigate the contributions of UCMHs formed during three phase transi- tions (i.e., electroweak symmetry breaking, QCD confinement and e+ e- annihilation) to the extragalactic γ-ray background. Moreover, we use the Fermi-LAT observation data of the extragalactic γ-ray background to get the constraints on the current abundance of UCMHs produced during these phase transitions. We also compare these results with those obtained from Cosmic Microwave Background (CMB) observations and find that the constraints from the Fermi-LAT are more stringent than those from CMB

preprint2011arXiv

Cosmological models with Lagrange Multiplier Field

We first consider the Einstein-aether theory with a gravitational coupling and a Lagrange multiplier field, and then consider the non-minimally coupled quintessence field theory with Lagrange multiplier field. We study the influence of the Lagrange multiplier field on these models. We show that the energy density evolution of the Einstein-aether field and the quintessence field are significantly modified. The energy density of the Einstein-aether is nearly a constant during the entire history of the Universe. The energy density of the quintessence field can also be kept nearly constant in the matter dominated Universe, or even exhibit a phantom-like behavior for some models. This suggests a possible dynamical origin of the cosmological constant or dark energy. Further more, for the canonical quintessence in the absence of gravitational coupling, we find that the quintessence scalar field can play the role of cold dark matter with the introduction of a Lagrange multiplier field. We conclude that the Lagrange multiplier field could play a very interesting and important role in the construction of cosmological models.

preprint2011arXiv

New Constraints on Primordial Minihalo Abundance Using Cosmic Microwave Background Observations

It was proposed that the massive compact halo objects (MCHOs) would be produced during an earlier epoch of cosmology if the density perturbations are between $3 \times 10^{-4}$ and 0.3. Then these objects can accrete dark matter particles onto them due to their high density. If the dark matter is in the form of the weakly interacting massive particles, the MCHOs can have a significant effect on the evolution of cosmology due to the dark matter annihilation within them. Using the WMAP-7 years data, we investigated the constraints on the current abundance of MCHOs ($f_{MCHOs}$) formed during the $e^{+} e^{-}$ annihilation phase transition. We have found that the $2σ$ constraint is $f_{MCHOs} \lesssim 10^{-4}$ for dark matter masses in the range between 1 GeV and 1 TeV.

preprint2011arXiv

Perturbation Theory of the Cosmological Log-Density Field

The matter density field exhibits a nearly lognormal probability density distribution (PDF) after entering into the nonlinear regime. Recently, it has been shown that the shape of the power spectrum of a logarithmically transformed density field is very close to the linear density power spectrum, motivating an analytic study of it. In this paper, we develop cosmological perturbation theory for the power spectrum of this field. Our formalism is developed in the context of renormalized perturbation theory, which helps to regulate the convergence behavior of the perturbation series, and of the Taylor- series expansion we use of the logarithmic mapping. This approach allows us to handle the critical issue of density smoothing in a straightforward way. We also compare our perturbative results with simulation measurements.

preprint2011arXiv

Relic Abundance of Asymmetric Dark Matter

We investigate the relic abundance of asymmetric Dark Matter particles that were in thermal equilibrium in the early universe. The standard analytic calculation of the symmetric Dark Matter is generalized to the asymmetric case. We calculate the asymmetry required to explain the observed Dark Matter relic abundance as a function of the annihilation cross section. We show that introducing an asymmetry always reduces the indirect detection signal from WIMP annihilation, although it has a larger annihilation cross section than symmetric Dark Matter. This opens new possibilities for the construction of realistic models of MeV Dark Matter.

preprint2011arXiv

The Abundance of New Kind of Dark Matter Structures

A new kind of dark matter structures, ultracompact minihalos (UCMHs) was proposed recently. They would be formed during the radiation dominated epoch if the large density perturbations are existent. Moreover, if the dark matter is made up of weakly interacting massive particles, the UCMHs can have effect on cosmological evolution because of the high density and dark matter annihilation within them. In this paper, one new parameter is introduced to consider the contributions of UCMHs due to the dark matter annihilation to the evolution of cosmology, and we use the current and future CMB observations to obtain the constraint on the new parameter and then the abundance of UCMHs. The final results are applicable for a wider range of dark matter parameters

preprint2011arXiv

The OH line contamination of 21 cm intensity fluctuation measurements for z=1~4

The large-scale structure of the Universe can be mapped with unresolved intensity fluctuations of the 21 cm line. The power spectrum of the intensity fluctuations has been proposed as a probe of the baryon acoustic oscillations at low to moderate redshifts with interferometric experiments now under consideration. We discuss the contamination to the low-redshift 21 cm intensity power spectrum generated by the 18 cm OH line since the intensity fluctuations of the OH line generated at a slightly higher redshift contribute to the intensity fluctuations observed in an experiment. We assume the OH megamaser luminosity is correlated with the star formation rate, and use the simulation to estimate the OH signal and the spatial anisotropies. We also use a semi-analytic simulation to predict the 21 cm power spectrum. At z=1 to 3, we find that the OH contamination could reach 0.1% to 1% of the 21 cm rms fluctuations at the scale of the first peak of the baryon acoustic oscillation. When z>3 the OH signal declines quickly, so that the contamination on the 21-cm becomes negligible at high redshifts.

preprint2011arXiv

Topology of large scale structure as test of modified gravity

The genus of the iso-density contours is a robust measure of the topology of large scale structure, and it is relatively insensitive to nonlinear gravitational evolution, galaxy bias and redshift-space distortion. We show that the growth of density fluctuations is scale-dependent even in the linear regime in some modified gravity theories, which opens a new possibility of testing the theories observationally. We propose to use the genus of the iso-density contours, an intrinsic measure of the topology of large scale structure, as a statistic to be used in such tests. In Einstein's general theory of relativity, density fluctuations are growing at the same rate on all scales in the linear regime, and the genus per comoving volume is almost conserved as structures are growing homologously, so we expect that the genus-smoothing scale relation is basically time-independent. However, in some modified gravity models where structures grow with different rates on different scales, the genus-smoothing scale relation should change over time. This can be used to test the gravity models with large scale structure observations. We studied the case of the f(R) theory, DGP braneworld theory as well as the parameterized post-Friedmann (PPF) models. We also forecast how the modified gravity models can be constrained with optical/IR or redshifted 21cm radio surveys in the near future.

preprint2010arXiv

Cosmic microwave background with Brans-Dicke gravity II: constraints with the WMAP and SDSS data

Using the covariant formalism developed in a companion paper (paper I), we derive observational constraint on the Brans-Dicke model in a flat FLRW universe with cosmological constant and cold dark matter. The CMB observations we use include the WMAP five year data, the ACBAR 2007 data, the CBI polarization data, and the BOOMERanG 2003 flight data. For the large scale structure we use the matter power spectrum data measured with the LRG survey of the SDSS DR4. We parameterize the Brans-Dicke parameter $ω$ with a new parameter $ζ=\ln(1/ω+1)$, and use the Markov-Chain Monte Carlo method to explore the parameter space. We find that using CMB data alone, one could place some constraint on positive $ζ$ or $ω$, but negative $ζ$ or $ω$ could not be constrained effectively. However, with additional large scale structure data, one could break the degeneracy at $ζ<0$. The $2σ$ (95.5%) bound on $ζ$ is $-0.00837<ζ<0.01018$ (corresponding to $ω< -120.0$ or $ω> 97.8$). We also obtained constraint on $\dot{G}/G$, the rate of change of $G$ at present, as $-1.75 \times 10^{-12} \yr^{-1}<\dot{G}/G < 1.05 \times 10^{-12}\yr^{-1}$, and $δG/G$, the total variation of $G$ since the epoch of recombination, as $ -0.083 < δ{G}/G < 0.095$ at $2σ$ confidence level.

preprint2010arXiv

Cosmic microwave background with Brans-Dicke Gravity: I. Covariant Formulation

In the covariant cosmological perturbation theory, a 1+3 decomposition ensures that all variables in the frame-independent equations are covariant, gauge-invariant and have clear physical interpretations. We develop this formalism in the case of Brans-Dicke gravity, and apply this method to the calculation of cosmic microwave background (CMB) anisotropy and large scale structures (LSS). We modify the publicly available Boltzmann code CAMB to calculate numerically the evolution of the background and adiabatic perturbations, and obtain the temperature and polarization spectra of the Brans-Dicke theory for both scalar and tensor mode, the tensor mode result for the Brans-Dicke gravity are obtained numerically for the first time. We first present our theoretical formalism in detail, then explicitly describe the techniques used in modifying the CAMB code. These techniques are also very useful to other gravity models. Next we compare the CMB and LSS spectra in Brans-Dicke theory with those in the standard general relativity theory. At last, we investigate the ISW effect and the CMB lensing effect in the Brans-Dicke theory. Constraints on Brans-Dicke model with current observational data is presented in a companion paper (paper II).

preprint2010arXiv

Cosmology with Photometric Surveys of Type Ia Supernovae

We discuss the extent to which photometric measurements alone can be used to identify Type Ia supernovae (SNIa) and to determine redshift and other parameters of interest for cosmological studies. We fit the light curve data of the type expected from a survey such as the one planned with Large Synoptic Survey Telescope (LSST) and also to remove the contamination from the core-collapse supernovae to SNIa samples. We generate 1000 SNIa mock flux data for each of the LSST filters based on existing design parameters, then use a Markov Chain Monte-Carlo (MCMC) analysis to fit for the redshift, apparent magnitude, stretch factor and the phase of the SNIa. We find that the model fitting works adequately well when the true SNe redshift is below 0.5, while at $z < 0.2$ the accuracy of the photometric data is almost comparable with spectroscopic measurements of the same sample. We discuss the contamination of Type Ib/c (SNIb/c) and Type II supernova (SNII) on the SNIa data set. We find it is easy to distinguish the SNII through the large $χ^2$ mismatch when fitting to photometric data with Ia light curves. This is not the case for SNIb/c. We implement a statistical method based on the Bayesian estimation in order to statistically reduce the contamination from SNIb/c for cosmological parameter measurements from the whole SNe sample. The proposed statistical method also evaluate the fraction of the SNIa in the total SNe data set, which provides a valuable guide to establish the degree of contamination.

preprint2010arXiv

Couplings between holographic dark energy and dark matter

We consider the interaction between dark matter and dark energy in the framework of holographic dark energy, and propose a natural and physically plausible form of interaction, in which the interacting term is proportional to the product of the powers of the dark matter and dark energy densities. We investigate the cosmic evolution in such models. The impact of the coupling on the dark matter and dark energy components may be asymmetric. While the dark energy decouples from the dark matter at late time, just as other components of the cosmic fluid become decoupled as the universe expands, interestingly, the dark matter may actually become coupled to the dark energy at late time. We shall name such a phenomenon as "incoupling". We use the latest type Ia supernovae data from the SCP team, baryon acoustics oscillation data from SDSS and 2dF surveys, and the position of the first peak of the CMB angular power spectrum to constrain the model. We find that the interaction term which is proportional to the the first power product of the dark energy and dark matter densities gives excellent fit to the current data.

preprint2010arXiv

Dark matter annihilation and non-thermal Sunyaev-Zel'dovich effect: II. dwarf spheroidal galaxy

We calculate the CMB temperature distortion due to the energetic electrons and positrons produced by dark matter annihilation (Sunyaev-Zel'dovich effect), in dwarf spheroidal galaxies (dSphs). In the calculation we have included two important effects which were previously ignored. First we show that the electron-positron pairs with energy less than GeV, which were neglected in previous calculation, could contribute a significant fraction of the total signal. Secondly we also consider the full effects of diffusion loss, which could significantly reduce the density of electron-positron pairs at the center of cuspy halos. For neutralinos, we confirm that detecting such kind of SZ effect is beyond the capability of the current or even the next generation experiments. In the case of light dark matter (LDM) the signal is much larger, but even in this case it is only marginally detectable with the next generation of experiment such as ALMA. We conclude that similar to the case of galaxy clusters, in the dwarf galaxies the $SZ_2DM}$ effect is not a strong probe of DM annihilations.

preprint2010arXiv

Distribution of satellite galaxies in high redshift groups

We use galaxy groups at redshifts between 0.4 and 1.0 selected from the Great Observatories Origins Deep Survey (GOODS) to study the color-morphological properties of satellite galaxies, and investigate possible alignment between the distribution of the satellites and the orientation of their central galaxy. We confirm the bimodal color and morphological type distribution for satellite galaxies at this redshift range: the red and blue classes corresponds to the early and late morphological types respectively, and the early-type satellites are on average brighter than the late-type ones. Furthermore, there is a {\it morphological conformity} between the central and satellite galaxies: the fraction of early-type satellites in groups with an early-type central is higher than those with a late-type central galaxy. This effect is stronger at smaller separations from the central galaxy. We find a marginally significant signal of alignment between the major axis of the early-type central galaxy and its satellite system, while for the late-type centrals no significant alignment signal is found. We discuss the alignment signal in the context of shape evolution of groups.

preprint2010arXiv

Leptonic dark matter annihilation in the evolving universe: constraints and implications

The cosmic electron and positron excesses have been explained as possible dark matter (DM) annihilation products. In this work we investigate the possible effects of such a DM annihilation scenario during the evolution history of the Universe. We first calculate the extragalactic $γ$-ray background (EGRB), which is produced through the final state radiation of DM annihilation to charged leptons and the inverse Compton scattering between electrons/positrons and the cosmic microwave background. The DM halo profile and the minimal halo mass, which are not yet well determined from the current N-body simulations, are constrained by the EGRB data from EGRET and Fermi telescopes. Then we discuss the impact of such leptonic DM models on cosmic evolution, such as the reionization and heating of intergalactic medium, neutral Hydrogen 21 cm signal and suppression of structure formation. We show that the impact on the Hydrogen 21 cm signal might show interesting signatures of DM annihilation, but the influence on star formation is not remarkable. Future observations of the 21 cm signals could be used to place new constraints on the properties of DM.

preprint2010arXiv

Searching for the earliest galaxies in the 21 cm forest

We use a model developed by Xu et al. (2010) to compute the 21 cm line absorption signatures imprinted by star-forming dwarf galaxies (DGs) and starless minihalos (MHs). The method, based on a statistical comparison of the equivalent width (W_ν) distribution and flux correlation function, allows us to derive a simple selection criteria for candidate DGs at very high (z >= 8) redshift. We find that ~ 18% of the total number of DGs along a line of sight to a target radio source (GRB or quasar) can be identified by the condition W_ν< 0; these objects correspond to the high-mass tail of the DG distribution at high redshift, and are embedded in large HII regions. The criterion W_ν> 0.37 kHz instead selects ~ 11% of MHs. Selected candidate DGs could later be re-observed in the near-IR by the JWST with high efficiency, thus providing a direct probe of the most likely reionization sources.

preprint2010arXiv

The earliest galaxies seen in 21 cm line absorption

We investigate the 21 cm absorption lines produced by non-linear structures during the early stage of reionization, i.e. the starless minihalos and the dwarf galaxies. After a detailed modelling of their properties, with particular attention to the coupling physics, we determine their 21 cm absorption line profiles. The infalling gas velocity around minihalos/dwarf galaxies strongly affects the line shape, and with the low spin temperatures outside the virial radii of the systems, gives rise to horn-like line profiles. The optical depth of a dwarf galaxy is reduced for lines of sight penetrating through its HII region, and especially, a large HII region created by a dwarf galaxy with higher stellar mass and/or a top-heavy initial mass function results in an optical depth trough rather than an absorption line. We compute synthetic spectra of 21 cm forest for both high redshift quasars and radio afterglows of gamma ray bursts (GRBs). Even with the planned SKA, radio afterglows of most if not all GRBs would still be too dim to be the background sources for high resolution (1 kHz) observations, but absorption lines can be easily detected towards a high-z quasar. Broadband observation against GRB afterglows can also be used to reveal the evolving 21 cm signal from both minihalos and dwarf galaxies if there was no X-ray background or it was extremely weak, but it becomes difficult if an early X-ray background existed. Hence the 21 cm absorption could be a powerful probe of the presence/intensity of the X-ray background and the thermal history of the early universe.

preprint2010arXiv

The Herschel-SPIRE Legacy Survey (HSLS): the scientific goals of a shallow and wide submillimeter imaging survey with SPIRE

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

preprint2009arXiv

Cosmological Implications of 5-dimensional Brans-Dicke Theory

The five dimensional Brans-Dicke theory naturally provides two scalar fields by the Killing reduction mechanism. These two scalar fields could account for the accelerated expansion of the universe. We test this model and constrain its parameter by using the type Ia supernova (SN Ia) data. We find that the best fit value of the 5-dimensional Brans-Dicke coupling contant is $ω= -1.9$. This result is also consistent with other observations such as the baryon acoustic oscillation (BAO).

preprint2009arXiv

Dark matter annihilation and non-thermal Sunyaev-Zel'dovich effect: I. galaxy cluster

In this work we calculate the Sunyaev-Zel'dovich (SZ) effect due to the $e^+e^-$ from dark matter (DM) annihilation in galaxy clusters. Two candidates of DM particle, (1) the weakly-interacting massive particle (WIMP) and (2) the light dark matter (LDM) are investigated. For each case, we also consider several DM profiles with and without central cusp. We generally find smaller signals than previously reported. Moreover, the diffusion of electrons and positrons in the galaxy clusters, which was generally thought to be negligible, is considered and found to have significant effect on the central electron/positron distribution for DM profile with large spatial gradient. We find that the SZ effect from WIMP is almost always non-observable, even for the highly cuspy DM profile, and using the next generation SZ interferometer such as ALMA. Although the signal of the LDM is much larger than that of the WIMP, the final SZ effect is still very small due to the smoothing effect of diffusion. Only for the configuration with large central cusp and extremely small diffusion effect, the LDM induced SZ effect might have a bit chance of being detected.

preprint2009arXiv

Features of holographic dark energy under the combined cosmological constraints

We investigate the observational signatures of the holographic dark energy models in this paper, including both the original model and a model with an interaction term between the dark energy and dark matter. We first delineate the dynamical behavior of such models, especially whether they would have "Big Rip" for different parameters, then we use several recent observations, including 182 high-quality type Ia supernovae data observed with the Hubble Space Telescope, the SNLS and ESSENCE surveys, 42 latest Chandra X-ray cluster gas mass fraction, 27 high-redshift gamma-ray burst samples, the baryon acoustic oscillation measurement from the Sloan Digital Sky Survey, and the CMB shift parameter from WMAP three years result to give more reliable and tighter constraints on the holographic dark energy models. The results of our constraints for the holographic dark energy model without interaction is $c=0.748^{+0.108}_{-0.009}$, $Ω_{\mathrm{m0}}=0.276^{+0.017}_{-0.016}$, and for model with interaction ($c=0.692^{+0.135}_{-0.107}$, $Ω_{\mathrm{m0}}=0.281^{+0.017}_{-0.017}$ ,$α=-0.006 ^{+0.021}_{-0.024}$, where $α$ is an interacting parameter). As these models have more parameters than the $Λ$CDM model, we use the Bayesian evidence as a model selection criterion to make comparison. We found that the holographic dark energy models are mildly favored by the observations compared with the $\mathrm{ΛCDM}$ model.

preprint2009arXiv

Primordial Non-Gaussianity from LAMOST Surveys

The primordial non-Gaussianity (PNG) in matter density perturbation is a very powerful probe of the physics of the very early Universe. The local PNG can induce a distinct scale-dependent bias on the large scale structure distribution of galaxies and quasars, which could be used for constraining it. We study the detection limits on PNG from the surveys of the LAMOST telescope. The cases of the main galaxy survey, the luminous red galaxy (LRG) survey, and the quasar survey of different magnitude limits are considered. We find that the MAIN1 sample (i.e. the main galaxy survey with one magnitude deeper than the SDSS main galaxy survey, or r<18.8) could only provide very weak constraint on PNG. For the MAIN2 sample (r<19.8) and the LRG survey, the 2σ(95.5%) limit on the PNG parameter f_{NL} are |f_{NL}|<145 and |f_{NL}|<114 respectively, comparable to the current limit from cosmic microwave background (CMB) data. The quasar survey could provide much more stringent constraint, and we find that the 2σlimit for |f_{NL}| is between 50 and 103, depending on the magnitude limit of the survey. With Planck-like priors on cosmological parameters, the quasar survey with g<21.65 would improve the constraints to |f_{NL}|<43 (2σ). We also discuss the possibility of further tightening the constraint by using the relative bias method proposed by Seljak(2008).

preprint2009arXiv

The 21 cm Forest as a Probe of the Reionization and the Temperature of the Intergalactic Medium

Using high redshift radio sources as background, the 21cm forest observations probe the neutral hydrogen absorption signatures of early structures along the lines of sight. Directly sensitive to the spin temperature, it complements the 21cm tomography observations, and provides information on the temperature as well as the ionization state of the intergalactic medium (IGM). We use a radiative transfer simulation to investigate the 21cm forest signals during the epoch of reionization. We first confirmed that the optical depth and equivalent width (EW) are good representations of the ionization and thermal state of the IGM. The features selected by their relative optical depth are excellent tracers of ionization fields, while those selected by absolute optical depth are very sensitive to the IGM temperature, so the IGM temperature information could potentially be extracted, thus breaking a degeneracy in 21cm tomographic observation. From the distributions of EWs and the number evolution of absorbers and leakers with different EWs, we see clearly the cosmological evolution of ionization state of the IGM. The number density of potentially observable features decreases rapidly with increasing gas temperature. The sensitivity of the proposed EW statistic to the IGM temperature makes it a unique and potentially powerful probe of reionization. Missing small-scale structures, such as small filaments and minihalos that are unresolved in our current simulation, and lack of an accurate calculation of the IGM temperature, however, likely have rendered the presented signals quantitatively inaccurate. Finally we discuss the requirements of the background radio sources for such observations, and find that signals with equivalent widths larger than 1kHz are hopeful to be detected.(Abridged)

preprint2009arXiv

The contribution of the IGM and minihalos to the 21 cm signal of reionization

We study the statistical properties of the cosmological 21 cm signal from both the intergalactic medium (IGM) and minihalos, using a reionization simulation that includes a self--consistent treatment of minihalo photoevaporation. We consider two models for minihalo formation and three typical thermal states of the IGM -- heating purely by ionization, heating from both ionizing and $\rm Lyα$ photons, and a maximal "strong heating" model. We find that the signal from the IGM is almost always dominant over that from minihalos. In our calculation, the differential brightness temperature, $δT_b,$ of minihalos is never larger than 2 mK. Although there are indeed some differences in the signals from the minihalos and from the IGM, even with the planned generation of radio telescopes it will be unfeasible to detect them. However, minihalos significantly affect the ionization state of the IGM and the corresponding 21 cm flux.

preprint2007arXiv

Detecting First Star Lyman-$α$ Spheres through Gravitational Telescopes

Lyman-$α$ spheres, i.e. regions around the first stars which are illuminated by Lyman-$α$ photons and show 21cm absorption feature against the CMB, are smoking guns at the dawn of the reionization epoch. Though overwhelming radio foreground makes their detections extremely difficult, we pointed out that, strong gravitational lensing can significantly improve their observational feasibility. Since Lyman-$α$ spheres have ~10" sizes, comparable to the caustic size of galaxy clusters, individual images of each strongly lensed Lyman-$α$ sphere often merge together and form single structures in the 21cm sky with irregular shapes. Using high-resolution N-body LCDM simulations, we found that the lensing probability to have magnification bigger than 10 is ~10^{-5}. This results in $\ga 10^6$ strongly lensed Lyman-$α$ spheres across the sky, which should be the primary targets for first detections of Lyman-$α$ spheres. Although the required total radio array collecting area for their detection is large (~100 km^2), the design of long fixed cylindrical reflectors can significantly reduce the total cost of such array to the level of the square kilometer array (SKA) and makes the detection of these very first objects feasible.

preprint2007arXiv

The 21cm Signature of the First Stars

We predict the 21-cm signature of the first metal-free stars. The soft X-rays emitted by these stars penetrate the atomic medium around their host halos, generating Lyman alpha photons that couple the spin and kinetic temperatures. These creates a region we call the Lyman alpha sphere, visible in 21-cm against the CMB, which is much larger than the HII region produced by the same star. The spin and kinetic temperatures are strongly coupled before the X-rays can substantially heat the medium, implying that a strong 21-cm absorption signal from the adiabatically cooled gas in Hubble expansion around the star is expected when the medium has not been heated previously. A central region of emission from the gas heated by the soft X-rays is also present although with a weaker signal than the absorption. The Lyman alpha sphere is a universal signature that should be observed around any first star illuminating its vicinity for the first time. The 21-cm radial profile of the Lyman alpha sphere can be calculated as a function of the luminosity, spectrum and age of the star. For a star of a few hundred solar masses and zero metallicity (as expected for the first stars), the physical radius of the Lyman alpha sphere can reach tens of kiloparsecs. The first metal-free stars should be strongly clustered because of high cosmic biasing; this implies that the regions producing a 21-cm absorption signal may contain more than one star and will generally be irregular and not spherical, because of the complex distribution of the gas. We discuss the feasiblity of detecting these Lyman alpha spheres, which would be present at redshifts $z\sim 30$ in the Cold Dark Matter model. Their observation would represent a direct proof of the detection of a first star.

preprint2004arXiv

Particle decays during the cosmic dark ages

We consider particle decays during the cosmic dark ages with two aims: (1) to explain the high optical depth reported by WMAP, and (2) to provide new constraints to the parameter space for decaying particles. We delineate the decay channels in which most of the decay energy ionizes and heats the IGM gas (and thus affects the CMB), and those in which most of the energy is carried away -- e.g. photons with energies 100 keV < E < 1 TeV -- and thus appears as a contribution to diffuse x-ray and gamma-ray backgrounds. The new constraints to the decay-particle parameters from the CMB power spectrum thus complement those from the cosmic X-ray and gamma-ray backgrounds. Although decaying particles can indeed produce an optical depth consistent with that reported by WMAP, in so doing they produce new fluctuations in the CMB temperature/polarization power spectra. For decay lifetimes less than the age of the Universe, the induced power spectra generally violate current constraints, while the power spectra are usually consistent if the lifetime is longer than the age of the Universe.

preprint1997arXiv

On the instability of the one-texture universe

The one-texture universe, introduced by Davis in 1987, is a homogeneous mapping of a scalar field with an $S^3$ vacuum into a closed universe. It has long been known to mathematicians that such solutions, although static, are unstable. We show by explicit construction that there is only one unstable mode, corresponding to collapse of the texture towards a single point, in the case where gravitational backreaction is neglected. We discuss the instability timescale in both static and expanding space-times; in the latter case it is of order of the present age of the universe, suggesting that, though unstable, the one-texture universe could survive to the present. The cosmic microwave background constrains the initial magnitude of this unstable perturbation to be less than of order 10^{-3}.