We examine the counts-in-cells probability distribution functions that describe dark matter halos in the Dark Energy Universe Simulations (DEUS). We describe the measurements between redshifts $z=0$ to $z=4$ on both linear and non-linear scales. The best-fits of the gravitational quasi-equilibrium distribution (GQED), the negative binomial distribution (NBD), the Poisson-Lognormal distribution (PLN), and the Poisson-Lognormal distribution with a bias parameter (PLNB) are compared to simulations. The fits agree reasonably consistently over a range of redshifts and scales. To distinguish quintessence (RPCDM) and phantom ($w$CDM) dark energy from $\Lambda$ dark energy, we present a new method that compares the model parameters of the counts-in-cells probability distribution functions. We find that the mean and variance of the halo counts-in-cells on $2-25h^{-1}$Mpc scales between redshifts $0.65<z<4$ show significant percentage differences for different dark energy cosmologies. On $15-25h^{-1}$Mpc scales, the $g$ parameter in NBD, $\omega$ parameter in PLN, $b$ and $C_b$ parameters in PLNB show larger percentage differences for different dark energy cosmologies than on smaller scales. On $2-6h^{-1}$Mpc scales, kurtosis and the $b$ parameter in the GQED show larger percentage differences for different dark energy cosmologies than on larger scales. For cosmologies explored in the DEUS simulations, the percentage differences between these statistics for the RPCDM and $w$CDM dark energy cosmologies relative to $\Lambda$CDM generally increases with redshift from a few percent to significantly larger percentages at $z=4$. Applying our method to simulations and galaxy surveys can provide a useful way to distinguish among dark energy models and cosmologies in general.

Similar to the compactness parameter ($\beta=M/R$), the gravitational binding energy (GBE) is also a characteristic parameter which can reflect the internal structure of a neutron star and thus can be used to expressing the universal relations. Scaling by the stellar mass, this investigation demonstrates a perfect universal relation between the GBE and the moment of inertia, where both of the normal neutron stars and the quark stars satisfy the same universal relation. Moreover, a fine empirical relation between the GBE and the tidal deformability is proposed, where the difference of the relations can be used to distinguish wether a pulsar is a normal neutron star or a quark star if the stellar mass and the tidal deformability can be observed or estimated rather accurately. These universal relations provide a potential way to estimate the GBE if the stellar mass and the moment of inertia/the tidal deformability are precisely measured.

Determining the Equation of State (EOS) of dense neutron-rich nuclear matter is a shared goal of both nuclear physics and astrophysics. Except possible phase transitions, the density dependence of nuclear symmetry \esym is the most uncertain part of the EOS of neutron-rich nucleonic matter especially at supra-saturation densities. Much progresses have been made in recent years in predicting the symmetry energy and understanding why it is still very uncertain using various microscopic nuclear many-body theories and phenomenological models. Simultaneously, significant progresses have also been made in probing the symmetry energy in both terrestrial nuclear laboratories and astrophysical observatories. In light of the GW170817 event as well as ongoing or planned nuclear experiments and astrophysical observations probing the EOS of dense neutron-rich matter, we review recent progresses and identify new challenges to the best knowledge we have on several selected topics critical for understanding astrophysical effects of the nuclear symmetry energy.

Within a minimum model for neutron stars consisting of nucleons, electrons and muons at $\beta$-equilibrium using about a dozen Equation of States (EOSs) from microscopic nuclear many-body theories and 40,000 EOSs randomly generated using an explicitly isospin-dependent parametric EOS model for high-density neutron-rich nucleonic matter within its currently known uncertainty range, we study correlations among the f-mode frequency, its damping time and the tidal deformability as well as the compactness of neutron stars. Except for quark stars, both the f-mode frequency and damping time of canonical neutron stars are found to scale with the tidal deformability independent of the EOSs used. Applying the constraint on the tidal deformability of canonical neutron stars $\Lambda_{1.4}=190^{+390}_{-120}$ extracted by the LIGO+VIRGO Collaborations from their improved analyses of the GW170817 event, the f-mode frequency and its damping time of canonical neutron stars are limited to 1.67 kHz - 2.18 kHz and 0.155 s - 0.255 s, respectively, providing a useful guidance for the ongoing search for gravitational waves from the f-mode oscillations of isolated neutron stars. Moreover, assuming either or both the f-mode frequency and its damping time will be measured precisely in future observations with advanced gravitational wave detectors, we discuss how information about the mass and/or radius as well as the still rather elusive nuclear symmetry energies at supra-saturation densities may be extracted.

By employing four typical equation of states (EOSs) of the inner crust matters, the constraint on properties of low-mass neutron stars are theoretically investigated. Based on the well known fact that there is a big gap between the neutron stars and white dwarfs in the mass-radius sequence of compact stars, here we explicitly point out the forbidden region for the central density and stellar radius to form a compact stars, that is, there is no compact star in the nature having central density in the region from about $10^{12}~kg/m^3$ to $10^{17}~kg/m^3$, and there is also no compact star having a radius in the region from about 400 km to 2000 km. It is found that for a low-mass neutron star with central densities at subnuclear densities, the stellar mass is mainly determined by the softness of the EOS at the higher density region ($ > 10^{16}~\textrm{kg.m}^{-3}$) while the stellar radius is dependent on the softness of the EOS at lower density region ($ < 10^{16}~\textrm{kg.m}^{-3}$). For a stable neutron star with stellar mass near the minimum mass, the Keplerian period is several hundred millisecond, the moment of inertia is in an order of $10^{37} kg\cdot m^{2}$, and the surface gravitational redshift is in an order of $10^{-4}$.

Dark matter is believed to be a major component of our universe. In this paper we propose a new mechanism based on dark matter inspired super-Chandrasekhar mass white dwarf to explain the recent observation of super luminous type Ia supernovae explosions. Our calculation shows when a white dwarf accretes enough dark matter, due to the Pauli exclusive principle between fermionic dark matter particles, the mass of corresponding dark white dwarf (which means the white dwarf mixed with dark matters) can significantly exceeds the Chandrasekhar limit. Moreover, we investigate some physical observable quantities, such as the redshift and moment of inertia of the dark white dwarf and found that these quantities are sensitive to the dark matter particle's distributions and thus can be potentially used to probe the relevant information of dark matter particles in the future.

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

A number of observed phenomena associated with individual neutron star systems or neutron star populations find explanations in models in which the neutron star crust plays an important role. We review recent work examining the sensitivity to the slope of the symmetry energy $L$ of such models, and constraints extracted on $L$ from confronting them with observations. We focus on six sets of observations and proposed explanations: (i) The cooling rate of the neutron star in Cassiopeia A, confronting cooling models which include enhanced cooling in the nuclear pasta regions of the inner crust, (ii) the upper limit of the observed periods of young X-ray pulsars, confronting models of magnetic field decay in the crust caused by the high resistivity of the nuclear pasta layer, (iii) glitches from the Vela pulsar, confronting the paradigm that they arise due to a sudden re-coupling of the crustal neutron superfluid to the crustal lattice after a period during which they were decoupled due to vortex pinning, (iv) The frequencies of quasi-periodic oscillations in the X-ray tail of light curves from giant flares from soft gamma-ray repeaters, confronting models of torsional crust oscillations, (v) the upper limit on the frequency to which millisecond pulsars can be spun-up due to accretion from a binary companion, confronting models of the r-mode instability arising above a threshold frequency determined in part by the viscous dissipation timescale at the crust-core boundary, and (vi) the observations of precursor electromagnetic flares a few seconds before short gamma-ray bursts, confronting a model of crust shattering caused by resonant excitation of a crustal oscillation mode by the tidal gravitational field of a companion neutron star just before merger.

Recent astrophysics observation reveals the existence of some super luminous type Ia supernovae. One natural explanation of such a peculiar phenomenon is to require the progenitor of such a supernova to be a highly super-Chandrasekhar mass white dwarf. Along this line, in this paper, we propose a possible mechanism to explain this phenomenon based on a charged white dwarf. In particular, by choosing suitable new variables and a representative charge distribution, an analytic solution is obtained. The stability issue is also discussed, remarkably, it turns out that the charged white dwarf configuration can be dynamically stable. Moreover, we investigate the general relativistic effects and it is shown that the general relativistic effects can be negligible when the mass of the charged white dwarf is below about $3M_\odot$.

A thorough understanding of many astrophysical phenomena associated with compact objects requires reliable knowledge about both the equation of state (EOS)of super-dense nuclear matter and the theory of strong-field gravity simultaneously because of the EOS-gravity degeneracy. Currently, variations of the neutron star (NS) mass-radius correlation from using alternative gravity theories are much larger than those from changing the NS matter EOS within known constraints. At least two independent observables are required to break the EOS-gravity degeneracy. Using model EOSs for hybrid stars and a Yukawa-type non-Newtonian gravity, we investigate both the mass-radius correlation and pulsating frequencies of NSs. While the maximum mass of NSs increases, the frequencies of the $f$, $p_1$, $p_2$, and $w_I$ pulsating modes are found to decrease with the increasing strength of the Yukawa-type non-Newtonian gravity, providing a useful reference for future determination simultaneously of both the strong-field gravity and the supranuclear EOS by combining data of x-ray and gravitational wave emissions of neutron stars.

The modeling of many neutron star observables incorporates the microphysics of both the stellar crust and core, which is tied intimately to the properties of the nuclear matter equation of state (EoS). We explore the predictions of such models over the range of experimentally constrained nuclear matter parameters, focusing on the slope of the symmetry energy at nuclear saturation density $L$. We use a consistent model of the composition and EoS of neutron star crust and core matter to model the binding energy of pulsar B of the double pulsar system J0737-3039, the frequencies of torsional oscillations of the neutron star crust and the instability region for r-modes in the neutron star core damped by electron-electron viscosity at the crust-core interface. By confronting these models with observations, we illustrate the potential of astrophysical observables to offer constraints on poorly known nuclear matter parameters complementary to terrestrial experiments, and demonstrate that our models consistently predict $L<70$ MeV.

The recently discovery of a massive neutron star (PSR J1614-2230 of $1.97\pm0.04M_{\odot}$) rules out the soft equation of states (EOSs) such as those included hyperons or kaon condensates at high densities, while the nuclear theory or the terrestrial laboratory data prefer a soft EOS. Here we propose one possible mechanism to allow that the observed massive neutron star can be supported by a soft EOS, that is, if the the gravitational constant $G$ varies at super strong field, a soft EOS can support the massive neutron stars.

Using a simple model of a neutron star with a perfectly rigid crust constructed and a set of crust and core equations of state that span the range of nuclear experimental uncertainty in the density dependence of the symmetry energy from 25 MeV (soft EOS) to 115 MeV (stiff EOS), we calculate the instability window for the onset of the Chandrasekhar-Friedmann-Schutz (CFS) instability in r-mode oscillations for canonical neutron stars (1.4 M_{\odot}) and massive neutron stars (2.0 M_{\odot}). In these models the crust-core transition density, and thus crustal thickness, is calculated consistently with the core equation of state (EOS). For the canonical neutron star, the lower bound of the r-mode instability window is reduced in frequency by \approx150 Hz from the softest to the stiffest symmetry energy used, independent of mass and temperature. The instability window also drops by \approx 100 Hz independent of EOS when the mass is raised from 1.4 M_{\odot} to 2.0 M_{\odot}. Where temperature estimates are available, the observed neutron stars in low mass X-ray binaries (LMXBs) have frequencies below the instability window for the 1.4 M_{\odot} models, while some LMXBs fall within the instability window for 2.0 M_{\odot} stars if the symmetry energy is relatively stiff, indicating that a softer symmetry energy is more consistent with observations within this model. Thus we conclude that smaller values of L help stabilize neutron stars against runaway r-mode oscillations. The critical temperature, below which no star can reach the instability window without exceeding its Kepler frequency, varies by nearly an order of magnitude from soft to stiff symmetry energies. When the crust thickness and core EOS are treated consistently, a thicker crust corresponds to a lower critical temperature, the opposite result to previous studies in which the transition density was independent of the core EOS.

Significant progress has been made in recent years in constraining the density dependence of nuclear symmetry energy using terrestrial nuclear laboratory data. Around and below the nuclear matter saturation density, the experimental constraints start to merge in a relatively narrow region. At supra-saturation densities, there are, however, still large uncertainties. After summarizing the latest experimental constraints on the density dependence of nuclear symmetry energy, we highlight a few recent studies examining imprints of nuclear symmetry energy on the binding energy, energy release during hadron-quark phase transitions as well as the $w$-mode frequency and damping time of gravitational wave emission of neutron stars.

Probably No! As an example, using soft EOSs consistent with existing terrestrial nuclear laboratory experiments for hybrid neutron stars containing a quark core described with MIT bag model using reasonable parameters, we show that the recently discovered new holder of neutron star maximum mass PSR J1614-2230 of $1.97\pm0.04M_{\odot}$ can be well described by incorporating a Yukawa gravitational correction that is consistent with existing constraints from neutron-proton and neutron-lead scatterings as well as the spectroscopy of antiproton atoms.

Describing the hyperonic and quark phases of neutron stars with an isospin- and momentum-dependent effective interaction for the baryon octet and the MIT bag model, respectively, and using the Gibbs conditions to construct the mixed phase, we study the energy release due to the hadron-quark phase transition. Moreover, the frequency and damping time of the first axial $w$-mode of gravitational waves are studied for both hyperonic and hybrid stars. We find that the energy release is much more sensitive to the bag constant than the density dependence of the nuclear symmetry energy. Also, the frequency of the $w$-mode is found to be significantly different with or without the hadron-quark phase transition and depends strongly on the value of the bag constant. Effects of the density dependence of the nuclear symmetry energy become, however, important for large values of the bag constant that lead to higher hadron-quark transition densities.

Considering the non-Newtonian gravity proposed in the grand unification theories, we show that the stability and observed global properties of neutron stars can not rule out the super-soft nuclear symmetry energies at supra-saturation densities. The degree of possible violation of the Inverse-Square-Law of gravity in neutron stars is estimated using an Equation of State (EOS) of neutron-rich nuclear matter consistent with the available terrestrial laboratory data.

The eigen-frequencies of the axial w-modes of oscillating neutron stars are studied using the continued fraction method with an Equation of State (EOS) partially constrained by the recent terrestrial nuclear laboratory data. It is shown that the density dependence of the nuclear symmetry energy $E_{sym}(\rho)$ affects significantly both the frequencies and the damping times of these modes. Besides confirming the previously found universal behavior of the mass-scaled eigen-frequencies as functions of the compactness of neutron stars, we explored several alternative universal scaling functions. Moreover, the $w_{II}$-mode is found to exist only for neutron stars having a compactness of $M/R\geq 0.1078$ independent of the EOS used.

Heavy-ion reactions especially those induced by radioactive beams provide useful information about the density dependence of the nuclear symmetry energy, thus the Equation of State of neutron-rich nuclear matter, relevant for many astrophysical studies. The latest developments in constraining the symmetry energy at both sub- and supra-saturation densities from analyses of the isopsin diffusion and the $\pi^-/\pi^+$ ratio in heavy-ion collisions using the IBUU04 transport model are discussed. Astrophysical ramifications of the partially constrained symmetry energy on properties of neutron star crusts, gravitational waves emitted by deformed pulsars and the w-mode oscillations of neutron stars are presented briefly.

The influence of the rotation on the total masses and radii of the neutron stars are calculated by the Hartle's slow rotation formalism, while the equation of state is considered in a relativistic $\sigma-\omega$ model. Comparing with the observation, the calculating result shows that the double neutron star binaries are more like hyperon stars and the neutron stars of X-ray binaries are more like traditional neutron stars. As the changes of the mass and radius to a real neutron star caused by the rotation are very small comparing with the total mass and radius, one can see that Hartle's approximate method is rational to deal with the rotating neutron stars. If three property values: mass, radius and period are observed to the same neutron star, then the EOS of this neutron star could be decided entirely.