A first-order quantum chromodynamics (QCD) phase transition (PT) may take place in the protocompact star (PCS) produced by a core-collapse supernova (CCSN). In this work, we study the consequences of such a PT in a non-rotating CCSN with axisymmetric hydrodynamic simulations. We find that the PT leads to the collapse of the PCS and results in a loud burst of gravitational waves (GWs). The amplitude of this GW burst is $\sim30$ times larger than the post-bounce GW signal normally found for non-rotating CCSN. It shows a broad peak at high frequencies ($\sim2500-4000$ Hz) in the spectrum, has a duration of $\lesssim5 {\rm ms}$, and carries $\sim3$ orders of magnitude more energy than the other episodes. Also, the peak frequency of the PCS oscillation increases dramatically after the PT-induced collapse. In addition to a second neutrino burst, the GW signal, if detected by the ground-based GW detectors, is decisive evidence of the first-order QCD PT inside CCSNe and provides key information about the structure and dynamics of the PCS.

Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

We perform two-dimensional hydrodynamic simulations for the thermonuclear explosion of Chandrasekhar-mass white dwarfs with dark matter (DM) cores in Newtonian gravity. We include a 19-isotope nuclear reaction network and make use of the pure turbulent deflagration model as the explosion mechanism in our simulations. Our numerical results show that the general properties of the explosion depend quite sensitively on the mass of the DM core M$_{{\rm DM}}$: a larger M$_{{\rm DM}}$ generally leads to a weaker explosion and a lower mass of synthesized iron-peaked elements. In particular, the total mass of $^{56}$Ni produced can drop from about 0.3 to 0.03 $M_{\odot}$ as M$_{{\rm DM}}$ increases from 0.01 to 0.03 $M_{\odot}$. We have also constructed the bolometric light curves obtained from our simulations and found that our results match well with the observational data of sub-luminous Type-Ia supernovae.

Recently observed pulsars with masses $\sim 1.1 ~M_{\odot}$ challenge the conventional neutron star (NS) formation path by core-collapse supernova (CCSN). Using spherically symmetric hydrodynamics simulations, we follow the collapse of a massive white dwarf (WD) core triggered by electron capture, until the formation of a proto-NS (PNS). For initial WD models with the same central density, we study the effects of a static, compact dark matter (DM) admixed core on the collapse and bounce dynamics and mass of the PNS, with DM mass $\sim 0.01 ~M_{\odot}$. We show that increasing the admixed DM mass generally leads to slower collapse and smaller PNS mass, down to about 1.0 $M_{\odot}$. Our results suggest that the accretion-induced collapse of dark matter admixed white dwarfs can produce low-mass neutron stars, such as the observed low-mass pulsar J0453+1559, which cannot be obtained by conventional NS formation path by CCSN.

The clustering amplitude of 7143 clusters from the Sloan Digital Sky Survey (SDSS) is found to increase with cluster mass, closely agreeing with the Gaussian random field hypothesis for structure formation. The amplitude of the observed cluster correlation exceeds the predictions from pure cold dark matter (CDM) simulation by $\simeq 6\%$ for the standard Planck-based values of the cosmological parameters. We show that this excess can be naturally accounted for by free streaming of light neutrinos, which opposes gravitational growth, so clusters formed at fixed mass are fewer and hence more biased than for a pure CDM density field. An enhancement of the cluster bias by 7\% matches the observations, corresponding to a total neutrino mass, $m_{\nu} = 0.119 \pm 0.034$ eV at 67\% confidence level, for the standard relic neutrino density. If ongoing laboratory experiments favor a normal neutrino mass hierarchy then we may infer a somewhat larger total mass than the minimum oscillation based value, $\sum m_{\nu} \simeq 0.056eV$, with 90\% confidence. Much higher precision can be achieved by applying our method to a larger sample of more distant clusters with weak lensing derived masses.

We present axisymmetric hydrodynamical simulations of accretion-induced collapse (AIC) of dark matter (DM) admixed rotating white dwarfs (WD) and their burst gravitational-wave (GW) signals. For initial WD models with the same central baryon density, the admixed DM is found to delay the plunge and bounce phases of AIC, and decrease the central density and mass of the proto-neutron star (PNS) produced. The bounce time, central density and PNS mass generally depend on two parameters, the admixed DM mass $M_\mathrm{DM}$ and the ratio between the rotational kinetic and gravitational energies of the inner core at bounce $\beta_\mathrm{ic,b}$. The emitted GWs have generic waveform shapes and the variation of their amplitudes $h_+$ show a degeneracy on $\beta_\mathrm{ic,b}$ and $M_\mathrm{DM}$. We found that the ratios between the GW amplitude peaks around bounce allow breaking the degeneracy and extraction of both $\beta_\mathrm{ic,b}$ and $M_\mathrm{DM}$. Even within the uncertainties of nuclear matter equation of state, a DM core can be inferred if its mass is greater than 0.03 $M_{\odot}$. We also discuss possible DM effects on the GW signals emitted by PNS g-mode oscillations. GW may boost the possibility for the detection of AIC, as well as open a new window in the indirect detection of DM.

Light but massive cosmological neutrinos do not cluster significantly on small scales, due to their high thermal velocities. With finite masses, cosmological neutrinos become part of the total matter field and contribute to its smoothing. Structure formation in the presence of massive neutrinos is therefore impeded compared to that in the standard $\Lambda$CDM cosmology with massless neutrinos. Neutrinos' masses also distort the anisotropy power spectrum of cosmic microwave background (CMB). Furthermore, a finite chemical potential $\mu$ for cosmological neutrinos, still allowed by current data, would have a non-negligible impact on CMB and structure formation. We consistently evaluate effects of neutrino masses and chemical potentials on the matter power spectrum by use of a neutrino-involved N-body simulation, with cosmological parameters obtained from a Markov-Chian Moonte-Carlo (MCMC) refitting of CMB data. Our results show that while a finite averaged neutrino mass $m_\nu$ tends to suppress the matter power spectrum in a range of wave numbers, the neutrino degeneracy parameters ${\xi_i \equiv \mu_i /T}$ ($i=$1, 2, 3) enhance the latter, leading to a large parameter degeneracy between $m_\nu$ and $\xi_i$. We provide an empirical formula for the effects on the matter power spectrum in a selected range of wave numbers induced by $m_\nu$ and $\eta \equiv \sqrt{\sum_i \xi^2_i}$. Observing a strong correlation between $m_\nu$ and $\eta$, we propose a single redshift-independent parameter $m_\nu - \frac{4}{3}\eta^2$ to characterize the neutrino effects on the matter power spectrum.

Fuzzy Dark Matter (FDM), motivated by string theory, has recently become a hot candidate for dark matter. The rest mass of FDM is believed to be $\sim 10^{-22}$eV and the corresponding de-Broglie wave length is $\sim 1$kpc. Therefore, the quantum effect of FDM plays an important role in structure formation. In order to study the cosmological structure formation in FDM model, several simulation techniques have been introduced. We review the current status and challenges in the cosmological simulation for the FDM model in this paper.

With recent Lyman-alpha forest data from BOSS and XQ-100, some studies suggested that the lower mass limit on the fuzzy dark matter (FDM) particles is lifted up to $10^{-21}\,\mathrm{eV}$. However, such a limit was obtained by $\Lambda$CDM simulations with the FDM initial condition and the quantum pressure (QP) was not taken into account which could have generated non-trivial effects in large scales structures. We investigate the QP effects in cosmological simulations systematically, and find that the QP leads to further suppression of the matter power spectrum at small scales, as well as the halo mass function in the low mass end. We estimate the suppressing effect of QP in the 1D flux power spectrum of Lyman-alpha forest and compare it with data from BOSS and XQ-100. The rough uncertainties of thermal gas properties in the flux power spectrum model calculation were discussed. We conclude that more systematic studies, especially with QP taken into account, are necessary to constrain FDM particle mass using Lyman-alpha forest.

Percolation analysis has long been used to quantify the connectivity of the cosmic web. Most of the previous work is based on density fields on grids. By smoothing into fields, we lose information about galaxy properties like shape or luminosity. Lack of mathematical model also limits our understanding of percolation analysis. In order to overcome these difficulties, we have studied percolation analysis based on discrete points. Using a Friends-of-Friends (FoF) algorithm, we generate the S-bb relation, between the fractional mass of the largest connected group (S) and the FoF linking length (bb). We propose a new model, the Probability Cloud Cluster Expansion Theory (PCCET) to relate the S-bb relation with correlation functions. We show that the S-bb relation reflects a combination of all orders of correlation functions. Using N-body simulation, we find that the S-bb relation is robust against redshift distortion and incompleteness in observation. From the Bolshoi simulation, with Halo Abundance Matching (HAM), we have generated a mock galaxy catalogue. Good matching of the projected two-point correlation function with observation is confirmed. However, comparing the mock catalogue with the latest galaxy catalogue from SDSS DR12, we have found significant differences in their S-bb relations. This indicates that the mock galaxy catalogue cannot accurately retain higher order correlation functions than the two-point correlation function, which reveals the limit of HAM method. As a new measurement, S-bb relation is applicable to a wide range of data types, fast to compute, robust against redshift distortion and incompleteness, and it contains information of all orders of correlation function.

The Ultra-Light Axion (ULA) is a dark matter candidate with mass $\mathcal{O}(10^{-22})$eV and de Broglie wavelength of order kpc. Such an axion, also called the Fuzzy Dark Matter (FDM), thermalizes via the gravitational force and forms a Bose-Einstein condensate. Recent studies suggested that the quantum pressure from the FDM can significantly affect the structure formation in small scales, thus alleviating the so-called "small-scale crisis". In this paper, we develop a new technique to discretize the quantum pressure and illustrate the interactions among FDM particles in the $N$-body simulation, which accurately simulates the formation of the dark-matter halo and its inner structure in the region outside the softening length. In a self-gravitationally-bound virialized halo, we find a constant density, solitonic core, which is consistent with the theoretical prediction. The existence of the solitonic core reveals the non-linear effect of quantum pressure and impacts the structure formation in the FDM model.

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

The angular momentum distribution in dark matter haloes and galaxies is a key ingredient in understanding their formation. Especially, the internal distribution of angular momenta is closely related to the formation of disk galaxies. In this article, we use haloes identified from a high-resolution simulation, the Bolshoi simulation, to study the spatial distribution of specific angular momenta, $j(r,\theta)$. We show that by stacking haloes with similar masses to increase the signal-to-noise ratio, the profile can be fitted as a simple function, $j(r,\theta)=j_s \sin^2(\theta/\theta_s) (r/r_s)^2/(1+r/r_s)^4 $, with three free parameters, $j_s, r_s$, and $\theta_s$. Specifically, $j_s$ correlates with the halo mass $M_\mathrm{vir}$ as $j_s\propto M_\mathrm{vir}^{2/3}$, $r_s$ has a weak dependence on the halo mass as $r_s \propto M_\mathrm{vir}^{0.040}$, and $\theta_s$ is independent of $M_\mathrm{vir}$. This profile agrees with that from a rigid shell model, though its origin is unclear. Our universal specific angular momentum profile $j(r,\theta)$ is useful in modelling haloes' angular momenta. Furthermore, by using an empirical stellar mass - halo mass relation, we can infer the averaged angular momentum distribution of a dark matter halo. The specific angular momentum - stellar mass relation within a halo computed from our profile is shown to share a similar shape as that from the observed disk galaxies.

We study the empirical relation between an astronomical object's angular momentum $J$ and mass $M$, $J=\beta M^\alpha$, the $J-M$ relation, using N-body simulations. In particular, we investigate the time evolution of the $J-M$ relation to study how the initial power spectrum and cosmological model affect this relation, and to test two popular models of its origin - mechanical equilibrium and tidal torque theory. We find that in the $\Lambda$CDM model, $\alpha$ starts with a value of $\sim 1.5$ at high redshift $z$, increases monotonically, and finally reaches $5/3$ near $z=0$, whereas $\beta$ evolves linearly with time in the beginning, reaches a maximum and decreases, and stabilizes finally. A three-regime scheme is proposed to understand this newly observed picture. We show that the tidal torque theory accounts for this time evolution behaviour in the linear regime, whereas $\alpha=5/3$ comes from the virial equilibrium of haloes. The $J-M$ relation in the linear regime contains the information of the power spectrum and cosmological model. The $J-M$ relations for haloes in different environments and with different merging histories are also investigated to study the effects of a halo's non-linear evolution. An updated and more complete understanding of the $J-M$ relation is thus obtained.

We construct a large, redshift complete sample of distant galaxy clusters by correlating Sloan Digital Sky Survey (SDSS) Data Release 12 (DR12) redshifts with clusters identified with the red-sequence Matched-filter Probabilistic Percolation (redMaPPer) algorithm. Our spectroscopic completeness is 97% for ~ 7,000 clusters within the redMaPPer selection limit, $z \leqslant$ 0.325, so that our cluster correlation functions are much more precise than earlier work and not suppressed by photometric redshifts. We derive an accurate power-law mass-richness relation from the observed abundance with respect to the mass function from Millennium XXL (MXXL) simulation, adjusted to the Planck weighted cosmology. The number density of clusters is found to decline by 20% over the range 0.1 $< z <$ 0.3, in good agreement with the evolution predicted by MXXL. Our projected three-dimensional correlation function scales with richness, $\lambda$, rising from $r_0=$ 14 h$^{-1}$ Mpc at $\lambda\simeq$ 25, to $r_0=$ 22 h$^{-1}$ Mpc at $\lambda\simeq$ 60, with a gradient that matches MXXL when applying our mass-richness relation, whereas the observed amplitude of the correlation function at $\left<z\right>=$ 0.24 exceeds the MXXL prediction by 20% at the $\simeq$ 2.5$\sigma$ level. This tension cannot be blamed on spurious, randomly located clusters as this would reduce the correlation amplitude. Full consistency between the correlation function and the abundances is achievable for the pre-Planck values of $\sigma_8=$ 0.9, $\Omega_m=$ 0.25, and h = 0.73, matching the improved distance ladder estimate of the Hubble constant.

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

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

The standard cold dark matter (CDM) model predicts too many and too dense small structures. We consider an alternative model that the dark matter undergoes two-body decays with cosmological lifetime $\tau$ into only one type of massive daughters with non-relativistic recoil velocity $V_k$. This decaying dark matter model (DDM) can suppress the structure formation below its free-streaming scale at time scale comparable to $\tau$. Comparing with warm dark matter (WDM), DDM can better reduce the small structures while being consistent with high redshfit observations. We study the cosmological structure formation in DDM by performing self-consistent N-body simulations and point out that cosmological simulations are necessary to understand the DDM structures especially on non-linear scales. We propose empirical fitting functions for the DDM suppression of the mass function and the mass-concentration relation, which depend on the decay parameters lifetime $\tau$ and recoil velocity $V_k$, and redshift. The fitting functions lead to accurate reconstruction of the the non-linear power transfer function of DDM to CDM in the framework of halo model. Using these results, we set constraints on the DDM parameter space by demanding that DDM does not induce larger suppression than the Lyman-$\alpha$ constrained WDM models. We further generalize and constrain the DDM models to initial conditions with non-trivial mother fractions and show that the halo model predictions are still valid after considering a global decayed fraction. Finally, we point out that the DDM is unlikely to resolve the disagreement on cluster numbers between the Planck primary CMB prediction and the Sunyaev-Zeldovich (SZ) effect number count for $\tau \sim H_{0}^{-1}$.

A two-dimensional hydrodynamics code for Type Ia supernovae (SNIa) simulations is presented. The code includes a fifth-order shock-capturing scheme WENO, detailed nuclear reaction network, flame-capturing scheme and sub-grid turbulence. For post-processing we have developed a tracer particle scheme to record the thermodynamical history of the fluid elements. We also present a one-dimensional radiative transfer code for computing observational signals. The code solves the Lagrangian hydrodynamics and moment-integrated radiative transfer equations. A local ionization scheme and composition dependent opacity are included. Various verification tests are presented, including standard benchmark tests in one and two dimensions. SNIa models using the pure turbulent deflagration model and the delayed-detonation transition model are studied. The results are consistent with those in the literature. We compute the detailed chemical evolution using the tracer particles' histories, and we construct corresponding bolometric light curves from the hydrodynamics results. We also use a Graphics Processing Unit (GPU) to speed up the computation of some highly repetitive subroutines. We achieve an acceleration of 50 times for some subroutines and a factor of 6 in the global run time.

The annihilation of dark matter (DM) particles in the Milky Way can contribute to the diffuse gamma-ray background (DGRB). Due to the presence of substructures, this emission will appear anisotropic in a predictable way. We generate full-sky maps of the gamma-ray emission in galactic substructures from results of the high-resolution Via Lactea II N-body simulation of the Milky Way DM halo. We calculate the anisotropy pattern, taking into account different radial profiles of the DM distribution in substructures, cosmic variance, and the detection threshold, and compare it to the anisotropy in the DGRB observed by the Fermi Large Area Telescope (LAT). By comparing the upper limits on the DM self-annihilation cross-section, <$\sigma v$>, implied by the anisotropy to the intensity of the DGRB and detected sources in the LAT 2-yr Point Source Catalog, we find that galactic substructure cannot contribute to the anisotropies in the DGRB without strongly violating these observations. Our results challenge the perception that small-scale anisotropies in the DGRB can be used as a probe of DM annihilation in galactic subhaloes.

One of the main goals of modern cosmic microwave background (CMB) missions is to measure the tensor-to-scalar ratio $r$ accurately to constrain inflation models. Due to ignorance about the reionization history $X_{e}(z)$, this analysis is usually done by assuming an instantaneous reionization $X_{e}(z)$ which, however, can bias the best-fit value of $r$. Moreover, due to the strong mixing of B-mode and E-mode polarizations in cut-sky measurements, multiplying the sky coverage fraction $f_{sky}$ by the full-sky likelihood would not give satisfactory results. In this work, we forecast constraints on $r$ for the Planck mission taking into account the general reionization scenario and cut-sky effects. Our results show that by applying an N-point interpolation analysis to the reionization history, the bias induced by the assumption of instantaneous reionization is removed and the value of $r$ is constrained within $5\%$ error level, if the true value of $r$ is greater than about 0.1 .

We study the equilibrium structures of white dwarfs with dark matter cores formed by non-self-annihilating dark matter DM particles with mass ranging from 1 GeV to 100 GeV, which are assumed to form an ideal degenerate Fermi gas inside the stars. For DM particles of mass 10 GeV and 100 GeV, we find that stable stellar models exist only if the mass of the DM core inside the star is less than O(10^-3) Msun and O(10^-6) Msun, respectively. The global properties of these stars, and in particular the corresponding Chandrasekhar mass limits, are essentially the same as those of traditional white dwarf models without DM. Nevertheless, in the 10 GeV case, the gravitational attraction of the DM core is strong enough to squeeze the normal matter in the core region to densities above neutron drip, far above those in traditional white dwarfs. For DM with particle mass 1 GeV, the DM core inside the star can be as massive as around 0.1 Msun and affects the global structure of the star significantly. In this case, the radius of a stellar model with DM can be about two times smaller than that of a traditional white dwarf. Furthermore, the Chandrasekhar mass limit can also be decreased by as much as 40%. Our results may have implications on to what extent type Ia supernovae can be regarded as standard candles - a key assumption in the discovery of dark energy.

With a simple power-law approximation of high-redshift ($\gtrsim3.5$) star formation history, i.e., $\dot{\rho}_*(z)\propto [(1+z)/4.5]^{-\alpha}$, we investigate the reionization of intergalactic medium (IGM) and the consequent Thomson scattering optical depth for cosmic microwave background (CMB) photons. A constraint on the evolution index $\alpha$ is derived from the CMB optical depth measured by the {\it Wilkinson Microwave Anisotropy Probe} (WMAP) experiment, which reads $\alpha\approx2.18\lg{\mathscr{N}_{\gamma}}-3.89$, where the free parameter $\mathscr{N}_\gamma$ is the number of the escaped ionizing ultraviolet photons per baryon. Moreover, the redshift for full reionization, $z_f$, can also be expressed as a function of $\alpha$ as well as $\mathscr{N}_{\gamma}$. By further taking into account the implication of the Gunn-Peterson trough observations to quasars for the full reionization redshift, i.e., $6\lesssim z_f \lesssim7$, we obtain $0.3\lesssim\alpha\lesssim1.3$ and $80\lesssim\mathscr{N}_{\gamma}\lesssim230$. For a typical number of $\sim4000$ of ionizing photons released per baryon of normal stars, the fraction of these photons escaping from the stars, $f_{\rm esc}$, can be constrained to within the range of $(2.0-5.8)%$.

If a component of cosmological dark matter is made up of massive particles - such as sterile neutrinos - that decay with cosmological lifetime to emit photons, the reionization history of the universe would be affected, and cosmic microwave background anisotropies can be used to constrain such a decaying particle model of dark matter. The optical depth depends rather sensitively on the decaying dark matter particle mass m_{dm}, lifetime tau_{dm}, and the mass fraction of cold dark matter f that they account for in this model. Assuming that there are no other sources of reionization and using the WMAP 7-year data, we find that 250 eV < m_{dm} < 1 MeV, whereas 2.23*10^3 yr < tau_{dm} < 1.23*10^18 yr. The best fit values for m_{dm} and tau_{dm}/f are 17.3 keV and 2.03*10^16 yr respectively.

In [Leung et al., Phys. Rev. D 84, 107301 (2011)], we presented our results on using a general relativistic two-fluid formalism to study the hydrostatic equilibrium configuration of an admixture of degenerate dark matter and normal nuclear matter. In this work, we present more analysis to complement our previous findings. We study the radial oscillation modes of these compact stars in detail. We find that these stars in general have two classes of oscillation modes. For a given total mass of the star, the first class of modes is insensitive to the dark-matter particle mass. They also reduce properly to the oscillation modes of the corresponding ordinary neutron star, with the same total mass, when the mass fraction of dark matter tends to zero. On the other hand, the second class of modes is due mainly to the dark-matter fluid. In the small dark-matter mass fraction limit, these modes are characterized purely by the oscillations of dark matter, while the normal matter is essentially at rest. In the intermediate regime where the mass fractions of the two fluids are comparable, the normal matter oscillates with the dark matter due to their coupling through gravity. In contrast to the first class of modes, the frequencies of these dark-matter dominated modes depend sensitively on the mass of dark-matter particles.

If massive sterile neutrinos exist, their decays into photons and/or electron-positron pairs may give rise to observable consequences. We consider the possibility that MeV sterile neutrino decays lead to the diffuse positron annihilation line in the Milky Way center, and we thus obtain bounds on the sterile neutrino decay rate $\Gamma_e \ge 10^{-28}$ s$^{-1}$ from relevant astrophysical/cosmological data. Also, we expect a soft gamma flux of $1.2 \times 10^{-4}-9.7 \times 10^{-4}$ ph cm$^{-2}$ s$^{-1}$ from the Milky Way center which shows up as a small MeV bump in the background photon spectrum. Furthermore, we estimate the flux of active neutrinos produced by sterile neutrino decays to be $0.02-0.1$ cm$^{-2}$ s$^{-1}$ passing through the earth.

We study the hydrostatic equilibrium configuration of an admixture of degenerate dark matter and normal nuclear matter by using a general relativistic two-fluid formalism. We consider non- self-annihilating dark matter particles of mass ~ 1 GeV. The mass-radius relations and moments of inertia of these dark-matter admixed neutron stars are investigated and the stability of these stars is demonstrated by performing a radial perturbation analysis. We find a new class of compact stars which consists of a small normal matter core with radius of a few km embedded in a ten-kilometer- sized dark matter halo. These stellar objects may be observed as extraordinarily small neutron stars that are incompatible with realistic nuclear matter models.

We show that the existence of a degenerate halo of sterile neutrinos with rest mass of 17.4 keV near the Galactic Center can account for both the excess 8.7 keV emission observed by the $Suzaku$ mission and the power needed ($10^{40}$ erg s$^{-1}$) to maintain the high temperature of the hot gas (8 keV) near the Galactic Center. The required decay rate and mixing angle of the sterile neutrinos are $\Gamma \ge 5 \times 10^{-20}$ s$^{-1}$ and $\sin^22 \theta \sim 10^{-3}-10^{-4}$ respectively. These values are consistent with a low reheating temperature in the inflation model, and suggest the exciting possibility that the sterile - active neutrino oscillation can be visible in near future experiments.

Models of a rotating universe have been studied widely since G{\"o}del \cite{1}, who showed an example that is consistent with General Relativity (GR). By now, the possibility of a rotating universe has been discussed comprehensively in the framework of some types of Bianchi's models, such as Type V, VII and IX \cite{2,3}, and different approaches have been proposed to constrain the rotation. Recent discoveries of some non-Gaussian properties of the Cosmic Microwave Background Anisotropies (CMBA) \cite{nG1,nG2,nG3,nG4,nG5,nG6,nG7}, such as the suppression of the quadrupole and the alignment of some multipoles draw attention to some Bianchi models with rotation \cite{bi1,bi2}. However, cosmological data, such as those of the CMBA, strongly prefer a homogeneous and isotropic model. Therefore, it is of interest to discuss the rotation of the universe as a perturbation of the Robertson-Walker metric, to constrain the rotating speed by cosmological data and to discuss whether it could be the origin of the non-Gaussian properties of the CMBA mentioned above. Here, we derive the general form of the metric (up to 2nd-order perturbations) which is compatible with the rotation perturbation in a flat $\Lambda$-CDM universe. By comparing the 2nd-order Sachs-Wolfe effect \cite{4,5,6,7,8} due to rotation with the CMBA data, we constrain the angular speed of the rotation to be less than $10^{-9}$ rad yr$^{-1}$ at the last scattering surface. This provides the first constraint on the shear-free rotation of a $\Lambda$CDM universe.

If sterile neutrinos exist and form halos in galactic centers, they can give rise to observational consequences. In particular, the sterile neutrinos decay radiatively and heat up the gas in the protogalaxy to achieve hydrostatic equilibrium, and they provide the mass to form supermassive blackholes. A natural correlation between the blackhole mass and velocity dispersion thus arises $\log(M_{BH,f}/M_{\odot})=\alpha \log (\sigma /200 {\rm km s^{-1}})+ \beta$ with $\alpha \approx 4$ and $\beta \approx 8$.

Based on the assumptions that a fraction of cluster dark matter is composed of degenerate neutrinos and they are in hydrostatic equilibrium with other matter, we predict a relation between the density profile and temperature of the cluster hot gas. The predicted relation agrees with observational data of 103 clusters.

Recent Chandra and Newton observations indicate that there are two-temperature components ($T \sim$ 8 keV, 0.8 keV) of the diffuse x-rays emitted from deep inside the center of Milky Way. We show that this can be explained by the existence of sterile neutrinos, which decay to emit photons that can be bound-free absorbed by the isothermal hot gas particles in the center of Milky Way. This model can account for the two-temperature components naturally as well as the energy needed to maintain the $\sim$ 8 keV temperature in the hot gas. The predicted sterile neutrino mass is between 16-18 keV.

In the standard framework of cosmology, primordial density fluctuations are assumed to have an isotropic Gaussian distribution. We search for deviations from this assumption in the WMAP data for the low $l$ modes of Cosmic Microwave Background Anisotropies (CMBA), by studying the directions of the z-axis that maximize the $l=m$ modes and the resulting amplitudes of these modes. We find a general alignment of the directions for $l=2$ to 10 modes to within 1/4 of the northern hemisphere. This alignment can be regarded as a generalization of the recently discovered alignment of the $l=2$ and 3 modes - the so-called `Axis of Evil'. Furthermore, we find abnormally high (low) powers in the $l=m=6$, 12 - 17 ($l=m=5$) modes; the probabilities for having the anomalous amplitudes of the $l=m=5$, 6, 17 modes are about 0.1%, 1% and 1% respectively according to the Gaussian conjecture. The alignment and anomalous amplitudes for these low $l$ modes are very robust against foreground contamination or different cleaning strategies, suggesting a cosmological origin and possibly supporting a non-standard inflation.

The constraints on the radion evolution in the Universal Extra Dimension (UED) model from Cosmic Microwave Background (CMB) and Type Ia supernovae (SNe Ia) data are studied. In the UED model, where both the gravity and standard model fields can propagate in the extra dimensions, the evolution of the extra dimensional volume, the radion, induces variation of fundamental constants. We discuss the effects of variation of the relevant constants in the context of UED for CMB power spectrum and SNe Ia data. We then use the three-year WMAP data to constrain the radion evolution at z \sim 1100, and the 2 \sigma constraint on \dot{\rho} / \rho_0 (\rho is a function of the radion, to be defined in the text) is [ -8.8, 6.6] \times 10 ^{-13} yr^-1. The SNe Ia gold sample yields a constraint on \dot{\rho} / \rho_0, for redshift between 0 and 1, to be [-4.7, 14] \times 10^{-13} yr^-1. Furthermore, the constraints from SNe Ia can be interpreted as bounds on the evolution QCD scale parameter, \dot{\Lambda}_{QCD} / \Lambda_{QCD, 0}, [-1.4, 2.8] \times 10^{-11} yr^-1, without reference to the UED model.

In many versions of brane model, the modulus field of extra dimensions, the radion, could have cosmological evolution, which induces variation of the Higgs vacuum expectation value, $<H>$, resulting in cosmological variation of the electron mass $m_e$. The formation of Cosmic Microwave Background (CMB) anisotropies is thus affected, causing changes both in the peaks positions and amplitudes in the CMB power spectra. Using the three-year Wilkinson Microwave Anisotropies Probe (WMAP) CMB data, with the Hubble parameter $H_0$ fixed to be the Hubble Space Telescope (HST) result 72 km s$^{-1}$ Mpc$^{-1}$, we obtain a constraint on $\rho$, the ratio of the value of $<H>$ at CMB recombination to its present value, to be [0.97, 1.02].

In this work we derive the covariant and gauge invariant perturbation equations in general theories of $f(R)$ gravity in the Palatini formalism to linear order and calculate the cosmic microwave background (CMB) and matter power spectra for an extensively discussed model, $f(R) = R + \alpha(-R)^\beta$, which is a possible candidate for the late-time cosmic accelerating expansion found recently. These spectra are discussed and found to be sensitively dependent on the value of $\beta$. We are thus able to make stringent constraints on $\beta$ from cosmological data on CMB and matter power spectra: The three-year Wilkinson Microwave Anisotropy Probe (WMAP) data alone gives a constraint $|\beta| \lesssim \mathcal{O}(10^{{-3}})$ while the joint WMAP, Supernova Lagacy Survey (SNLS) and Sloan Digital Sky Survey (SDSS) data sets tightens this to $\beta \sim \mathcal{O}(10^{{-6}})$, about an order of magnitude more stringent than the constraint from SDSS data alone, which makes this model practically indistinguishable from the standard $\Lambda\mathrm{CDM}$ paradigm.

We use the Cosmic Microwave Background Anisotropies (CMBA) power spectra to constrain the cosmological variation of gravitational constant G. It is found that the sensitivity of CMBA to the variation of G is enhanced when G is required to converge to its present value. The variations of G from the CMB decoupling epoch z ~ 1000 to the present time are modelled by a step function and a linear function of scale factor $a$ respectively, and the corresponding 95% confidence intervals for G/G_0 are [0.95, 1.05] and [0.89, 1.13], G_0 being the present value. The CMBA constraint is unique in the sense that it entails the range of redshift from z \approx 1000 to 0.

We calculate in this article the CMB and matter power spectra of a class of early $f(R)$ cosmologies, which takes the form of $f(R) = R + \lambda_1 H_0^2\text{exp}[R/(\lambda_2 H_0^2)]$. Unlike the late-time $f(R)$ cosmologies such as $f(R) = R + \alpha(-R)^\beta$ ($\beta<0$), the deviation from $\Lambda\text{CDM}$ of this model occurs at a higher redshift (thus the name \emph{Early $f(R)$ Cosmology}), and this important feature leads to rather different ISW effect and CMB spectrum. The matter power spectrum of this model is, at the same time, again very sensitive to the chosen parameters, and LSS observations such as SDSS should constrain the parameter space stringently. We expect that our results are applicable at least qualitatively to other models that produce $f(R)$ modification to GR at earlier times (\emph{e.g.}, redshifts $\mathcal{O}(10) \lesssim z \lesssim \mathcal{O}(1)$) than when dark energy begins to dominate -- such models are strongly disfavored by data on CMB and matter power spectra.

We study hydrodynamic evolution of cosmological background neutrinos. By using a spherically symmetric Newtonian hydrodynamic code, we calculate the time evolution of the density profiles of neutrino matter in cluster and galactic scales. We discuss the possible observational consequences of such evolution and the resulting density profiles of the degenerate neutrino `stars' in galaxies and clusters.

Assuming that there exists a species of heavy sterile neutrinos($m_{\nu}>1$ keV) and their decays can serve as a heating source of the hot gas in galaxy clusters, we study how the observational constraints on cooling flow limit the mass of these sterile neutrinos. We predict a relation among the luminosity, total mass and the redshift of a cluster, and we compare this relation with data from 12 clusters to obtain an estimate of the decay rate of the sterile neutrino.

We present a method for reducing cosmological data to constraints on the amplitudes of modes of the dark energy density as a function of redshift. The modes are chosen so that (1) one of them has constant density and (2) the others are non-zero only if there is time-variation in the dark energy density and (3) the amplitude errors for the time-varying modes are uncorrelated with each other. We apply our method to various combinations of three-year WMAP data, baryon acoustic oscillation data, the 'Gold' supernova data set, and the Supernova Legacy Survey data set. We find no significant evidence for a time-varying dark energy density or for non-zero mean curvature. Although by some measure the limits on four of the time-varying mode amplitudes are quite tight, they are consistent with the expectation that the dark energy density does not vary on timescales shorter than a Hubble time. Since we do not expect detectable time variation in these modes, our results should be viewed as a systematic error test which the data have passed. We discuss a procedure to identify modes with maximal signal-to-noise ratio.

(Abridged) We study the gravitational radiation from gravitational collapses of rapidly rotating neutron stars induced by a phase-transition from normal nuclear matter to a mixed phase of quark and nuclear matter in the core of the stars. The study is based on self-consistent three dimensional hydrodynamic simulations with Newtonian gravity and a high resolution shock capturing scheme, and the quadrupole formula of gravitational radiation. The quark matter of the mixed phase is described by the MIT bag model and the normal nuclear matter is described by an ideal fluid equation of state (EOS). 1. We determined the magnitudes of the emitted gravitational waves for several collapse scenarios. 2. We determined the types and frequencies of the fluid oscillation modes excited by the process. In particular, we find that the gravitational wave signals produced by the collapses are dominated by the fundamental quadrupole and quasi-radial modes of the final equilibrium configurations. In some collapse scenarios, we find that the oscillations are damped out within a few dynamical timescales due to the growth of differential rotations and the formation of strong shock waves. 3. We showed that the spectrum of the gravitational wave signals is sensitive to the EOS, implying that the detection of such gravitational waves could provide useful constraints on the EOS of newly born quark stars. 4. For the range of rotation periods we have studied, we found no sign of the development of nonaxisymmetric dynamical instabilities in the collapse process.

We consider the big bang nucleosynthesis (BBN) in the Brane world scenario, where all matter fields are confined on our 3-brane and the radion of the Brane evolves cosmologically. In the Einstein frame fundamental fermion masses vary and the results of standard BBN (SBBN) are modified. We can thus use the observational primordial element abundances to impose constraints on the possible variations of the radion. The possibility of using the evolving radion to resolve the discrepancies between the Wilkinson Microwave Anisotropy Probe (WMAP) and SBBN values of the baryon-to-photon ratio ($\eta$) is also discussed. The results and constraints presented here are applicable to other models in which fundamental fermion masses vary.

The successful prediction of light element abundances from Big Bang Nucleosynthesis has been a pillar of the standard model of Cosmology. Because many of the relevant reaction rates are sensitive to the values of fundamental constants, such as the fine structure constant and the strong coupling constant, BBN is a useful tool to probe and to put constraints on possible cosmological variations of these constants, which arise naturally from many versions of extra-dimensional theories. In this paper, we study the dependences of fundamental constants on the radion field of the universal extra dimension model, and calculate the effects of such varying constants on BBN. We show that the discrepancy between BBN and the Wilkinson Microwave Anisotropy Probe data on the baryon-to-photon ratio may be reduced if the volume of the extra dimensions was slightly larger - by $\mathcal{O} (10^{-3})$ - at the BBN era compared to its present value, which would result in smaller gauge couplings at BBN by the same factor.

One of the most important and surprising discoveries in cosmology in recent years is the realization that our Universe is dominated by a mysterious dark energy, which leads to an accelerating expansion of space-time. A simple generalization of the standard Friedmann-Robertson-Walker equations based on General Relativity and the Cosmological Principles with the inclusion of a number of closed extra dimensions reproduces the currently observed data on dark energy, without the introduction of any cosmological constant or new particles. In particular, with a few extra dimensions, we obtain the redshift dependence of the deceleration parameter, the dark energy equation of state, as well as the age of the Universe in agreement with data, with essentially no free parameter. This model predicts that the extra dimensions have been compactified throughout the cosmic history, and it therefore suggests that signals from early Universe may give promising signatures of extra dimensions.

We study the Blackwell-Rao (BR) estimator of the probability distribution of the angular power spectrum, P(C_l|d), by applying it to samples of full-sky no-noise CMB maps generated via Gibbs sampling. We find the estimator, given a set of samples, to be very fast and also highly accurate, as determined by tests with simulated data. We also find that the number of samples required for convergence of the BR estimate rises rapidly with increasing l, at least at low l. Our existing sample chains are only long enough to achieve convergence at l less than about 40. In comparison with P(C_l|d) as reported by the WMAP team we find significant differences at these low l values. These differences lead to up to approximately 0.5 sigma shifts in the estimates of parameters in a 7-parameter Lambda CDM model with non-zero dn_s/dln(k). Fixing dn_s/dln(k)= 0 makes these shifts much less significant. Unlike existing analytic approximations, the BR estimator can be straightforwardly extended for the case of power spectra from correlated fields, such as temperature and polarization. We discuss challenges to extending the procedure to higher l and provide some solutions.

We study the cooling behaviour of an isolated strange quark star, using an equation of state derived from perturbative QCD up to second order in strong coupling constant, and we compare it with that of a neutron star. After an initial rapid cooling, a quark star may undergo the QCD phase transition to become a neutron star. We propose several signatures for such a scenario: a large amount of energy can be released due to latent heat, a long duration $\gamma$-ray source, and a second neutrino burst after a supernova explosion.

Cosmological parameters affect observables in physically distinct ways. For example, the baryon density, omega_b, affects the ionization history and also the pressure of the pre-recombination fluid. To investigate the relative importance of different physical effects to the determination of omega_b, and to test the cosmological model, we artificially split omega_b into two `metaparameters': omega_{be} which controls the ionization history and omega_{bp} which plays the role of omega_b for everything else. In our demonstration of the technique we find omega_b = .0229 +/- .0012 (with no parameter splitting), omega_{bp} = .0238 +/- .0021, omega_{be}= .0150 +/- .0034 and omega_{bp}-omega_{be} = .0088 +/- .0039.

We study the scenario that a new born strange quark star cools to the Quantum ChromoDynamics (QCD) phase transition temperature and converts to a neutron star, and we calculate the evolution of temperature and luminosity of the compact star. We argue that the conversion energy released can be of the order $10^{53}$ erg. We also propose that a second neutrino burst will be emitted at the completion of this phase transition.

Recent Strong lensing data and rotation curves of dwarf galaxies indicate that many galactic clusters may have a soft core instead of a central cusp in their density distribution. This result challenges the standard CDM (Cold Dark Matter) theory prediction based on N-body simulations. We find that the observed density profile is consistent with that of a spherical gas of degenerate fermions such as neutrinos in hydrostatic equilibrium. Also, we compare two different models of the dark matter halo and their predictions for the hot gas profile.

The atomic recombination process leads to a softening of the matter equation of state as reflected by a reduced generalized adiabatic index, with accompanying heat release. We study the effects of this recombination softening and reheating of the cosmic plasma on the ionization history, visibility function, Cold Dark Matter (CDM) transfer function, and the Cosmic Microwave Background (CMB) spectra. The resulting modifications of the CMB spectrm is 1/10 of WMAP's current error and is comparable to PLANCK's error. Therefore, this effect should be considered when data with higher accuracy are analysed.

The normal parameters are a non--linear transformation of the cosmological parameters whose likelihood function is very well--approximated by a normal distribution. This transformation serves as an extreme form of data compression allowing for practically instantaneous calculation of the likelihood of any given model, as long as the model is in the parameter space originally considered. The compression makes all the information about cosmological parameter constraints from a given set of experiments available in a useable manner. Here we explicitly define the normal parameters that work for the current CMB data, and give their mean and covariance matrix which best fit the likelihood function calculated by the Monte Carlo Markov Chain method. Along with standard parameter estimation results, we propose that future CMB parameter analyses define normal parameters and quote their mean and covariance matrix.

We study the gravitational wave signals emitted from phase-transition induced collapses of rapidly rotating neutron stars to strange stars by performing 3D numerical simulations. Our preliminary results suggest that the complete conversion of neutron stars to strange stars would occur within a fraction of millisecond. We also find that the gravitational waves generated from the collapse process may be detectable by the advanced LIGO for reasonable source distance. In addition, the study such gravitational wave signals would put useful constraint on the parameters of QCD.

The recent discovery of a Gunn--Peterson (GP) trough in the spectrum of the redshift 6.28 SDSS quasar has raised the tantalizing possibility that we have detected the reionization of the universe. However, a neutral fraction (of hydrogen) as small as 0.1% is sufficient to cause the GP trough, hence its detection alone cannot rule out reionization at a much earlier epoch. The Cosmic Microwave Background (CMB) polarization anisotropy offers an alternative way to explore the dark age of the universe. We show that for most models constrained by the current CMB data and by the discovery of a GP trough (showing that reionization occurred at z > 6.3), MAP can detect the reionization signature in the polarization power spectrum. The expected 1-sigma error on the measurement of the electron optical depth is around 0.03 with a weak dependence on the value of that optical depth. Such a constraint on the optical depth will allow MAP to achieve a 1-sigma error on the amplitude of the primordial power spectrum of 6%. MAP with two years (Planck with one year) of observation can distinguish a model with 50% (6%) partial ionization between redshifts of 6.3 and 20 from a model in which hydrogen was completely neutral at redshifts greater than 6.3. Planck will be able to distinguish between different reionization histories even when they imply the same optical depth to electron scattering for the CMB photons.

As the chiral symmetry is widely recognized as an important driver of the strong interaction dynamics, current strange stars models based on MIT bag models do not obey such symmetry. We investigate properties of bare strange stars using the Cloudy Bag Model, in which a pion cloud coupled to the quark-confining bag is introduced such that chiral symmetry is conserved. We find that in this model the decay of pions is a very efficient cooling way. In fact it can carry out most the thermal energy in a few milliseconds and directly convert them into 100MeV photons via pion decay. This may be a very efficient $\gamma$-ray burst mechanism. Furthermore, the cooling behavior may provide a possible way to distinguish a compact object between a neutron star, MIT strange star and Cloudy Bag strange star in observations.

The energy source has remained to be the great mystery in understanding of the gamma-ray bursts (GRBs) if the events are placed at cosmological distances as indicated by a number of recent observations. The currently popular models include (1)the merger of two neutron stars or a neutron star and a black hole binary and (2)the hypernova scenario of the collapse of a massive member in a close binary. Since a neutron star will inevitably collapse into a black hole if its mass exceeds the limit $M_{max}\approx3M_{\odot}$, releasing a total binding gravitational energy of $\sim10^{54}$ erg, we explore semi-empirically the possibility of attributing the energy source of GRB to the accretion- induced collapse of a neutron star (AICNS) in a massive X-ray binary system consisting of a neutron star and a type O/B companion. This happens because a significant mass flow of $\sim10^{-3}$--$10^{-4}M_{\odot}$ yr$^{-1}$ may be transferred onto the neutron star through the Roche-lobe overflow and primarily during the spiral-in phase when it plunges into the envelope of the companion, which may eventually lead to the AICNS before the neutron star merges with the core of the companion. In this scenario, a ``dirty'' fireball with a moderate amount of beaming is naturally expected because of the nonuniformity of the stellar matter surrounding the explosion inside the companion, and a small fraction ($\sim0.1%$) of the energy is sufficient to create the observed GRBs. In addition, the bulk of the ejecting matter of the companion star with a relatively slow expansion rate may act as the afterglow. Assuming a non-evolutionary model for galaxies, we estimate that the birthrate of the AICNS events is about 2 per day within a volume to redshift $z=1$ for an $\Omega_0=1$ universe, consistent with the reported GRB rate.