Researcher profile

M. -C. Chu

M. -C. Chu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

astro-ph.CO astro-ph gr-qc hep-lat astro-ph.HE astro-ph.GA astro-ph.IM astro-ph.SR nucl-th

Trust snapshot

Quick read

Trust 21 - Emerging

18works

0followers

9topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

preprint2020arXiv

Metamorphosis of dwarf halo density profile under dark matter decays

We study the density profile of a dwarf halo in the decaying dark matter (DDM) cosmology, using a new algorithm that resolves halo density profiles down to the innermost $700$ pc robustly with high efficiency. Following Schwarzschild's orbit-based method, we have also developed a simplified model to calculate the DDM halo density profiles, which agree remarkably well with those from N-body simulations. Both zoom-in simulations and the simplified model reveal that dark matter decays lead to the flattening of central density and overall reduction of density in dwarf halos, and the underlying physics mechanisms are well illustrated by the simplified model. The slowly-rising scaled rotation curves of DDM dwarf halos agree with the observation of local dwarf galaxies. Our results suggest that the DDM holds great potential for resolving the small-scale problems faced by the cold dark matter (CDM) model.

preprint2015arXiv

A new hydrodynamics code for Type Ia Supernovae

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.

preprint2015arXiv

Cosmological Structure Formation in Decaying Dark Matter Models

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 $τ$ 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 $τ$. 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 $τ$ 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-$α$ 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 $τ\sim H_{0}^{-1}$.

preprint2014arXiv

Can galactic dark matter substructure contribute to the cosmic gamma-ray anisotropy?

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, <$σ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.

preprint2014arXiv

Cosmic microwave background constraints on the tensor-to-scalar ratio

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 .

preprint2013arXiv

Dark-matter admixed white dwarfs

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.

preprint2012arXiv

Cosmic Microwave Background constraints of decaying dark matter particle properties

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.

preprint2012arXiv

Equilibrium Structure and Radial Oscillations of Dark Matter Admixed Neutron Stars

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.

preprint2011arXiv

Dark-Matter Admixed Neutron Stars

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.

preprint2011arXiv

Electron-positron Annihilation Lines and Decaying Sterile Neutrinos

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 $Γ_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.

preprint2009arXiv

Is the universe rotating?

Models of a rotating universe have been studied widely since G{ö}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 $Λ$-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 $Λ$CDM universe.

preprint2008arXiv

Decaying Sterile Neutrinos as a Heating Source in the Milky Way Center

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.

preprint2008arXiv

The Existence of Sterile Neutrino Halos in Galactic Centers as an Explanation of the Black Hole mass - Velocity Dispersion Relation

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})=α\log (σ/200 {\rm km s^{-1}})+ β$ with $α\approx 4$ and $β\approx 8$.

preprint2006arXiv

Gravitational Waves from Phase-Transition Induced Collapse of Neutron Stars

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

preprint2006arXiv

Heating the intergalactic medium by radiative decay of neutrinos

Assuming that there exists a species of heavy sterile neutrinos($m_ν>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.

preprint1997arXiv

Temperature dependence of instantons in QCD

We investigate the temperature dependence of the instanton contents of gluon fields, using unquenched lattice QCD and the cooling method. The instanton size parameter deduced from the correlation function decreases from 0.44fm below the phase-transition temperature $T_c$ ($\approx 150$MeV) to 0.33fm at 1.3 $T_c$. The instanton charge distribution is Poissonian above $T_c$, but it deviates from the convoluted Poisson at low temperature. The topological susceptibility decreases rapidly below $T_c$, showing the apparent restoration of the $U(1)_A$ symmetry already at $T \approx T_c$.

preprint1994arXiv

Instanton content of finite temperature QCD matter

We investigate the temperature dependence of the instanton content of gluon fields and their contribution to quark correlation using quenched lattice QCD and the cooling method. We found a suppression of the topological susceptibility at finite temperature, agreeing with the PCAC expectation at low temperature and enhanced suppression for temperatures at and above the deconfinement transition. For temperatures up to about 334 MeV, the topological charge correlation agrees well with a single instanton profile, though the size parameter seems to change across the phase transition. The screening wave functions for the pion and the rho become slightly more compact at higher temperatures. Lattice cooling shows no contribution from the instantons to the screening wave functions at temperatures above the phase transition.

preprint1992arXiv

Lattice Analysis of Two-Point Hadronic Correlators in the QCD Vacuum

Results from the first lattice QCD analysis of vacuum correlators of local hadronic currents using dispersion relations are presented. We have explored the vector, pseudoscalar, axial, and scalar meson channels, and the proton-like and delta-like baryon channels. The lattice results are shown to agree qualitatively with experimental results in channels where experimental data exist, and shed insight into interacting instanton approximations and sum rule calculations in the other channels. 3 figures included (as PostScript files): pi.ps, nd.ps, va.ps : Search for %%% to retrieve the files. Latex file, needs espcrc2.sty (Search for %%%).