The gravitational instabilities are important to the evolution of the disks and the planet formation in the disks. We calculate the evolution of the disks which form from the collapse of the molecular cloud cores. By changing the properties of the cloud cores and the hydrodynamical viscosity parameters, we explore their effects on the properties of the gravitational instabilities. We find that the disk is unstable when the angular velocity of the molecular cloud core is larger than a critical value. The time duration of the instability increases as the angular velocity of the core increases. The increase of the hydrodynamical viscosity parameter hardly affects the stability of the disk, but decreases the time duration of the critical state of the gravitational instability in the disk. The instability of the disks can happen at very early time of evolution of the disk, which is consistent with the observations.

We calculate the cross-correlation function $\langle (\Delta T/T)(\mathbf{v}\cdot \mathbf{n}/\sigma_{v}) \rangle$ between the kinetic Sunyaev-Zeldovich (kSZ) effect and the reconstructed peculiar velocity field using linear perturbation theory, to constrain the optical depth $\tau$ and peculiar velocity bias of central galaxies with Planck data. We vary the optical depth $\tau$ and the velocity bias function $b_{v}(k)=1+b(k/k_{0})^{n}$, and fit the model to the data, with and without varying the calibration parameter $y_{0}$ that controls the vertical shift of the correlation function. By constructing a likelihood function and constraining $\tau$, $b$ and $n$ parameters, we find that the quadratic power-law model of velocity bias $b_{v}(k)=1+b(k/k_{0})^{2}$ provides the best-fit to the data. The best-fit values are $\tau=(1.18 \pm 0.24) \times 10^{-4}$, $b=-0.84^{+0.16}_{-0.20}$ and $y_{0}=(12.39^{+3.65}_{-3.66})\times 10^{-9}$ ($68\%$ confidence level). The probability of $b>0$ is only $3.12 \times 10^{-8}$ for the parameter $b$, which clearly suggests a detection of scale-dependent velocity bias. The fitting results indicate that the large-scale ($k \leq 0.1\,h\,{\rm Mpc}^{-1}$) velocity bias is unity, while on small scales the bias tends to become negative. The value of $\tau$ is consistent with the stellar mass--halo mass and optical depth relation proposed in the previous literatures, and the negative velocity bias on small scales is consistent with the peak background-split theory. Our method provides a direct tool to study the gaseous and kinematic properties of galaxies.

In this paper, we develop a full statistical method for the pairwise velocity estimator previously proposed, and apply Cosmicflows-2 catalogue to this method to constrain cosmology. We first calculate the covariance matrix for line-of-sight velocities for a given catalogue, and then simulate the mock full-sky surveys from it, and then calculate the variance for the pairwise velocity field. By applying the $8315$ independent galaxy samples and compressed $5224$ group samples from Cosmicflows-2 catalogue to this statistical method, we find that the joint constraint on $\Omega^{0.6}_{\rm m}h$ and $\sigma_{8}$ is completely consistent with the WMAP 9-year and Planck 2015 best-fitting cosmology. Currently, there is no evidence for the modified gravity models or any dynamic dark energy models from this practice, and the error-bars need to be reduced in order to provide any concrete evidence against/to support $\Lambda$CDM cosmology.

Based on our previous work, we deduce a general formula for pressure of degenerate and relativistic electrons,Pe, which is suitable for superhigh magnetic fields, discuss the quantization of Landau levels of electrons, and consider the quantum electrodynam-ic(QED) effects on the equations of states (EOSs) for different matter systems. The main conclusions are as follows:Pe is related to the magnetic field B, matter density ?, and electron fraction Ye ; the stronger the magnetic field, the higher the electron pressure becomes; the high electron pressure could be caused by high Fermi energy of electrons in a superhigh magnetic field; compared with a common radio pulsar, a magnetar could be a more compact oblate spheroid-like deformed neutron star due to the anisotropic total pressure; and an increase in the maximum mass of a magnetar is expected because of the positive contribution of the magnetic field energy to the EOS of the star.

In a previous study, we formulated a framework of the entropy-based equilibrium statistical mechanics for self-gravitating systems. This theory is based on the Boltzmann-Gibbs entropy and includes the generalized virial equations as additional constraints. With the truncated distribution function to the lowest order, we derived a set of second-order equations for the equilibrium states of the system. In this work, the numerical solutions of these equations are investigated. It is found that there are three types of solutions for these equations. Both the isothermal and divergent solutions are thermally unstable and have unconfined density profiles with infinite mass, energy and spatial extent. The convergent solutions, however, seem to be reasonable. Although the results cannot match the simulation data well, because of the truncations of the distribution function and its moment equations, these lowest-order convergent solutions show that the density profiles of the system are confined, the velocity dispersions are variable functions of the radius, and the velocity distributions are also anisotropic in different directions. The convergent solutions also indicate that the statistical equilibrium of self-gravitating systems is by no means the thermodynamic equilibrium. These solutions are just the lowest-order approximation, but they have already manifested the qualitative success of our theory. We expect that higher-order solutions of our statistical-mechanical theory will give much better agreement with the simulation results concerning dark matter haloes.

The long-standing puzzle surrounding the statistical mechanics of self-gravitating systems has not yet been solved successfully. We formulate a systematic theoretical framework of entropy-based statistical mechanics for spherically symmetric collisionless self-gravitating systems. We use an approach that is very different from that of the conventional statistical mechanics of short-range interaction systems. We demonstrate that the equilibrium states of self-gravitating systems consist of both mechanical and statistical equilibria, with the former characterized by a series of velocity-moment equations and the latter by statistical equilibrium equations, which should be derived from the entropy principle. The velocity-moment equations of all orders are derived from the steady-state collisionless Boltzmann equation. We point out that the ergodicity is invalid for the whole self-gravitating systems, but it can be re-established locally. Based on the local ergodicity, using Fermi-Dirac-like statistics, with the nondegenerate condition and the spatial independence of the local microstates, we rederive the Boltzmann-Gibbs entropy. This is consistent with the validity of the collisionless Boltzmann equation, and should be the correct entropy form for collisionless self-gravitating systems. Apart from the usual constraints of mass and energy conservation, we demonstrate that the series of moment or virialization equations must be included as additional constraints on the entropy functional when performing the variational calculus; this is an extension to the original prescription by White & Narayan. Any possible velocity distribution can be produced by the statistical-mechanical approach that we have developed with the extended Boltzmann-Gibbs/White-Narayan statistics. Finally, we discuss the questions of negative specific heat and ensemble inequivalence for self-gravitating systems.

The statistical mechanics of self-gravitating systems has not been well understood, and still remains an open question so far. In a previous study by Kang & He, we showed that the fluid approximation may give a clue to further investigate this problem. In fact, there are indeed many dynamical similarities between self-gravitating and fluid systems. Based on a fluid-like entropy, that work explained successfully the outer density profiles of dark matter halos, but there left some drawbacks with the calculation concerning extremizing process of the entropy. In the current paper, with the improved extremizing calculation -- including an additional differential constraint of dynamical equilibrium and without any other assumptions, we confirm that statistical-mechanical methods can give a density profile with finite mass and finite energy. Moreover, this density profile is also consistent with the observational surface brightness of the elliptical galaxy NGC 3379. In our methods, the density profile is derived from the equation of state, which is obtained from entropy principle but does not correspond to the maximum entropy of the system. Finally, we suggest an alternative entropy form, a hybrid of Boltzmann-Gibbs and Tsallis entropy, whose global maximum may give rise to this equation of state.

In this Letter, we investigate the stability of the statistical equilibrium of spherically symmetric collisionless self-gravitating systems. By calculating the second variation of the entropy, we find that perturbations of the relevant physical quantities should be classified as long- and short-range perturbations, which correspond to the long- and short-range relaxation mechanisms, respectively. We show that the statistical equilibrium states of self-gravitating systems are neither maximum nor minimum, but complex saddle-point entropy states, and hence differ greatly from the case of ideal gas. Violent relaxation should be divided into two phases. The first phase is the entropy-production phase, while the second phase is the entropy-decreasing phase. We speculate that the second-phase violent relaxation may just be the long-wave Landau damping, which would work together with short-range relaxations to keep the system equilibrated around the saddle-point entropy states.

Cosmological N-body simulations have revealed many empirical relationships of dark matter halos, yet the physical origin of these halo properties still remains unclear. On the other hand, the attempts to establish the statistical mechanics for self-gravitating systems have encountered many formal difficulties, and little progress has been made for about fifty years. The aim of this work is to strengthen the validity of the statistical-mechanical approach we have proposed previously to explain the dark matter halo properties. By introducing an effective pressure instead of the radial pressure to construct the specific entropy, we use the entropy principle and proceed in a similar way as previously to obtain an entropy stationary equation. An equation of state for equilibrated dark halos is derived from this entropy stationary equation, by which the dark halo density profiles with finite mass can be obtained. We also derive the anisotropy parameter and pseudo-phase-space density profile. All these predictions agree well with numerical simulations in the outer regions of dark halos. Our work provides further support to the idea that statistical mechanics for self-gravitating systems is a viable tool for investigation.

(abbreviated) The statistical mechanics of self-gravitating systems is a long-held puzzle. In this work, we employ a phenomenological entropy form of ideal gas, first proposed by White & Narayan, to revisit this issue. By calculating the first-order variation of the entropy, subject to the mass- and energy-conservation constraints, we obtain an entropy stationary equation. Incorporated with the Jeans equation, and by specifying some functional form for the anisotropy parameter beta, we numerically solve the two equations, and demonstrate that the velocity anisotropy parameter plays an important role to attain a density profile that is finite in mass, energy, and spatial extent. If incorporated again with some empirical density profile from simulations, our theoretical predictions of the anisotropy parameter, and the radial pseudo-phase-space density in the outer non-gravitationally degenerate region of the dark matter halo, agree rather well with the simulation data, and the predictions are also acceptable in the middle weak-degenerate region of the dark halo. The second-order variational calculus reveals the seemingly paradoxical but actually complementary consequence that the equilibrium state of self-gravitating systems is the global minimum entropy state for the whole system under long-range violent relaxation, but simultaneously the local maximum entropy state for every and any small part of the system under short-range two-body relaxation and Landau damping. This minimum-maximum entropy duality means that the standard second law of thermodynamics needs to be re-expressed or generalized for self-gravitating systems. We believe that our findings, especially the complementary second law of thermodynamics, may provide crucial clues to the development of the statistical physics of self-gravitating systems as well as other long-range interaction systems.

N-body simulations of dark matter halos show that the density is cusped near the center of the halo. The density profile behaves as $r^{-\gamma}$ in the inner parts, where $\gamma \simeq 1$ for the NFW model and $\gamma \simeq 1.5$ for the Moore's model, but in the outer parts, both models agree with each other in the asymptotic behavior of the density profile. The simulations also show the information about anisotropy parameter $\beta(r)$ of velocity distribution. $\beta\approx 0$ in the inner part and $\beta\approx 0.5$ (radially anisotropic) in the outer part of the halo. We provide some distribution functions $F(E,L)$ with the constant anisotropy parameter $\beta$ for the two spherical models of dark matter halos: a new generalized NFW model and a generalized Moore model. There are two parameters $\alpha$ and $\epsilon$ for those two generalized models to determine the asymptotic behavior of the density profile. In this paper, we concentrate on the situation of $\beta(r)=1/2$ from the viewpoint of the simulation.

N-body simulations of dark matter halos show that the density profiles of halos behave as $\rho(r)\propto r^{-\alpha(r)}$, where the density logarithmic slope $\alpha \simeq 1\sim1.5$ in the center and $\alpha \simeq 3\sim 4$ in the outer parts of halos. However, some observations are not in agreement with simulations in the very central region of halos. The simulations also show that velocity dispersion anisotropy parameter $\beta\approx 0$ in the inner part of the halo and the so called "pseudo phase-space density" $\rho/\sigma^3$ behaves as a power-law in radius $r$. With these results in mind, we study the distribution function and the pseudo phase-space density $\rho/\sigma^3$ of the center of dark matter halos and find that they are closely-related.

We study the X-ray emission of baryon fluid in the universe using the WIGEON cosmological hydrodynamic simulations. It has been revealed that cosmic baryon fluid in the nonlinear regime behaves like Burgers turbulence, i.e. the fluid field consists of shocks. Like turbulence in incompressible fluid, the Burgers turbulence plays an important role in converting the kinetic energy of the fluid to thermal energy and heats the gas. We show that the simulation sample of the $\Lambda$CDM model without adding extra heating sources can fit well the observed distributions of X-ray luminosity versus temperature ($L_{\rm x}$ vs. $T$) of galaxy groups and is also consistent with the distributions of X-ray luminosity versus velocity dispersion ($L_{\rm x}$ vs. $\sigma$). Because the baryonic gas is multiphase, the $L_{\rm x}-T$ and $L_{\rm x}-\sigma$ distributions are significantly scattered. If we describe the relationships by power laws $L_{\rm x}\propto T^{\alpha_{LT}}$ and $L_{\rm x}\propto \sigma^{\alpha_{LV}}$, we find $\alpha_{LT}>2.5$ and $\alpha_{LV}>2.1$. The X-ray background in the soft $0.5-2$ keV band emitted by the baryonic gas in the temperature range $10^5<T<10^7$ K has also been calculated. We show that of the total background, (1) no more than 2% comes from the region with temperature less than $10^{6.5}$ K, and (2) no more than 7% is from the region of dark matter with mass density $\rho_{\rm dm}<50 \bar{\rho}_{\rm dm}$. The region of $\rho_{\rm dm}>50\bar{\rho}_{\rm dm}$ is generally clustered and discretely distributed. Therefore, almost all of the soft X-ray background comes from clustered sources, and the contribution from truly diffuse gas is probably negligible. This point agrees with current X-ray observations.

We investigate the statistical properties of cosmic baryon fluid in the nonlinear regime, which is crucial for understanding the large-scale structure formation of the universe. With the hydrodynamic simulation sample of the Universe in the cold dark matter model with a cosmological constant, we show that the intermittency of the velocity field of cosmic baryon fluid at redshift z=0 in the scale range from the Jeans length to about 16 Mpc/h can be extremely well described by She-Leveque's universal scaling formula. The baryon fluid also possesses the following features: (1) for volume weight statistics, the dissipative structures are dominated by sheets, and (2) the relation between the intensities of fluctuations is hierarchical. These results imply that the evolution of highly evolved cosmic baryon fluid is similar to a fully developed turbulence.

The dynamic evolution of the baryonic intergalactic medium (IGM) caused by the underlying dark matter gravity is governed by the Navier-Stokes equations in which many cooling and heating processes are involved. However, it has long been recognized that the growth mode dynamics of cosmic matter clustering can be sketched by a random force driven Burgers' equation if cooling and heating are ignored. Just how well the dynamics of the IGM can be described as a Burgers fluid has not been fully investigated probably because cooling and heating are essential for a detailed understanding of the IGM. Using IGM samples produced by a cosmological hydrodynamic simulation in which heating and cooling processes are properly accounted for, we show that the IGM velocity field in the nonlinear regime shows the features of a Burgers fluid, that is, when the Reynolds number is high, the velocity field consists of an ensemble of shocks. Consequently, (1) the IGM velocity $v$ is generally smaller than that of dark matter; (2) for the smoothed field, the IGM velocity shows tight correlation with dark matter given by $v \simeq s v_{dm}$, with $s<1$, such that the lower the redshift, the smaller $s$; (3) the velocity PDFs are asymmetric between acceleration and deceleration events; (4) the PDF of velocity difference $\Delta v=v(x+r)-v(x)$ satisfies the scaling relation for a Burgers fluid, i.e., $P(\Delta v)=(1 r^y)F(\Delta v/r^y)$. We find the scaling function and parameters for the IGM which are applicable to the entire scale range of the samples (0.26 - 8 h$^{-1}$ Mpc). These properties show that the similarity mapping between the IGM and dark matter is violated on scales much larger than the Jeans length of the IGM.

We show that, in the halo model of large-scale structure formation, the difference between the Fourier and the DWT (discrete wavelet transform) power spectra provides a statistical measurement of the halos. This statistical quantity is free from parameters related to the shape of the mass profile and the identification scheme of halos. That is, the statistical measurement is invariant in the sense that models with reasonably defined and selected parameters of the halo models should yield the same difference of the Fourier and DWT spectra. This feature is useful to extract ensemble averaged properties of halos, which cannot be obtained with the identification of individual halo. To demonstrate this point, we show with WIGEON hydrodynamical simulation samples that the spectrum difference provides a quantitative measurement of the discrepancy of the distribution of baryonic gas from that of the underlying dark matter field within halos. We also show that the mass density profile of halos in physical space can be reconstructed with this statistical measurement. This profile essentially is the average over an ensemble of halos, including well virialized halos as well as halos with significant internal substructures. Moreover, this reconstruction is sensitive to the tail of the mass density profile. We showed that the profile with $1/r^3$ tail gives very different result from that of $1/r^2$. Other possible applications of this method are discussed as well.

(Abridged) The nonlinear evolution of a system consisting of baryons and dark matter is generally characterized by strong shocks and discontinuities. The baryons slow down significantly at postshock areas of gravitational strong shocks, which can occur in high overdense as well as low overdense regions. Consequently, the baryon fraction would be nonuniform on large scales. We studied these phenomena with simulation samples produced by the WENO hybrid cosmological hydrodynamic/N-body code. We find that the baryon fraction in high mass density regions is lower on average than the cosmic baryon fraction, and many baryons accumulate in the regions with moderate mass density to form a high baryon fraction phase (HBFP). In dense regions with rho>100, which are the possible hosts for galaxy clusters, the baryon fraction can be lower than the cosmic baryon fraction by about 10%--20% at z ~ 0. Our simulation samples show that about 3% of the cosmic baryon budget was hidden in the HBFP at redshift z=3, while this percentage increases to about 14% at the present day. The gas in the HBFP cannot be detected either by Ly-alpha forests of QSO absorption spectra or by soft X-ray background. That is, the HBFP would be missed in the baryon budget given by current observations.

We investigate the 21 cm emission from the epoch between reionization z_r and Gunn-Peterson transparency z_{gp}. According to the lognormal model of the thermal history around reionization, hydrogen clouds in z_r > z > z_{gp} are hot and a predominant part of baryonic gas is ionized, but still opaque to Ly-alpha photons. Therefore, 21 cm emission is a distinctive characteristic of this epoch. We show that the 21 cm emission comes from both uncollapsed and collapsing hydrogen clouds. The spatial distribution of the brightness temperature excess $\delta T_b$ is highly non-Gaussian. It consists of spikes with high $\delta T_b$ and a low $\delta T_b$ area between the spikes. The field has the following statistical features: (1) the one-point distributions of $\delta T_b$ are described approximately by power-law tailed probability distribution functions; (2) the nth-order moment of $\delta T_b$ is increasing much faster with n than that of a Gaussian field, but slower than that of a lognormal field; (3) the scale-scale correlation of the $\delta T_b$ field is significant for all scales larger than the Jeans length of the gas. These features would be useful for distinguishing the 21 cm emission of the early clustering from the noise of foreground contamination.

(Abridged)The temperature (T) and entropy (S) fields of baryonic gas/IGM are analyzed using simulation samples by a cosmological hydro/N-body code. We show that in the nonlinear regime the dynamical similarity between the IGM and dark matter will be broken by strong shocks. The heating and entropy production by the shocks breaks the IGM into multiple phases. The mean entropy, or the cosmological entropy floor, is found to be more than 100 h^(-1/3) keV cm^2 in all regions when z < 1. At redshift z ~ 2 - 3, high-entropy gas (S > 50 h^(-1/3) keV cm^2) mostly resides in areas on scales larger than 1 h^(-1) Mpc and with density $\rho_{dm}>10^2$. Therefore, gravitational shocks are an effective preheating mechanism of the IGM, and probably enough to provide the entropy excess of clusters and groups if the epoch of the gas falling in cluster cores is not earlier than z ~ 2 - 3. On the other hand, at redshifts z < 4, there is always a more than 90% volume of the low dark matter mass density regions filled by the IGM with T < 10^4.5 K. Therefore, the multiphased character and non-Gaussianity of the IGM field would explain the high-temperature and high-entropy gas observed in groups and clusters.

We investigate the constraint on the abundance of primordial black holes (PBHs) and the spectral index $n$ of primeval density perturbations given by the ionizing photon background at the epoch of reionization. Within the standard inflationary cosmogony, we show that the spectral index $n$ of the power-law power spectrum of primeval density perturbations should be $n<$1.27. Since the universe is still optical thick at the reionization redshift $z\sim 6$ - 8, this constraint is independent of the unknown parameter of reheating temperature of the inflation. The ionizing photon background from the PBHs can be well approximated by a power law spectrum $J(\nu)\propto{\nu}^3$, which is greatly different from those given by models of massive stars and quasars.

Pure luminosity evolution models for galaxies provide an unacceptable fit to the redshifts and colors of faint galaxies. In this paper we demonstrate, using HST morphological number counts derived both from the $I_{814}$-band of WFPC2 in the Medium Deep Survey (MDS) and the Hubble Deep Field (HDF) and from the $H_{1.6}$-band of NICMOS, and ground-based spectroscopic data of the Hawaii Deep Field and the Canada-France Redshift Survey, that number evolution is necessary for galaxies, regardless of whether the cosmic geometry is flat, open, or $\Lambda$-dominated. Furthermore, we show that the number evolution is small at redshifts of z<1, but large at z>1, and that this conclusion is valid for all the three cosmological models under consideration. If the universe is open or $\Lambda$-dominated, the models, which are subject to the constraint of the conservation of the comoving mass density of galaxies, naturally predict a population of star-forming galaxies with the redshift distribution peaking at $z=2\sim3$, which seems to be consistent with the recent findings from Lyman-break photometric selection techniques. If the cosmological model is flat, however, the conservation of the comoving mass density is invalid. Hence, in order to account for the steep slope of B-band number counts at faint magnitudes in the flat universe, such a star-forming galaxy population has to be introduced ad hoc into the modelling alongside the merger assumption.

The explanation for the existence of an excess population of faint blue galaxies (FBGs) has been a mystery for nearly two decades, and remains one of the grand astronomical issues to date. Existing models cannot explain all of the observational data such as galaxy number counts in the optical and infrared passbands and the redshift distributions of galaxies. Here, by modelling the morphological number counts derived from the Hubble Space Telescope, as well as the number counts in optical and infrared passbands, and the redshift and color distributions of galaxies obtained from ground-based observations, we show that the `FBG problem' cannot be resolved if elliptical galaxies are assumed to have formed in an instantaneous burst of star formation at high redshift with no subsequent star formation events, which is just the conventional scenario for formation and evolution of ellipticals. There exist great discrepancies between the observed $B-K$ color distribution and the predicted distribution for ellipticals by such a pure luminosity evolution (PLE) model in the context of the conventional scenario. Neither can the mild evolution (i.e., the star formation events have lasted for a longer time than those of the instantaneous burst and passive evolution since the formation of galaxies) for ellipticals be accepted in the context of PLE assumption. The introduction of dust extinction also cannot save the PLE models. This conclusion holds for each of the three cosmological models under consideration: flat, open and $\Lambda$-dominated. Hence, our investigation suggests that PLE assumption for elliptical galaxies is questionable, and number evolution may be essential for ellipticals.

In this paper, we explore the plausible luminosity evolution of early-type galaxies in different cosmological models by constructing a set of pure luminosity evolution (PLE) models via the choices of the star formation rate (SFR) parameters and formation redshift $z_f$ of galaxies, with the observational constraints derived from the Hubble Space Telescope (HST) morphological number counts for elliptical and S0 galaxies of the Medium Deep Survey (MDS) and the Hubble Deep Field (HDF). We find that the number counts of early-type galaxies can be explained by the pure luminosity evolution models, without invoking exotic scenarios such as merging or introducing an additional population. But the evolution should be nearly passive, with a high $z_f$ assumed. The conclusion is valid in all of the three cosmological models we adopted in this paper. We also present the redshift distributions for three bins of observed magnitudes in F814w pass-band, to show at which redshift are the objects that dominate the counts at a given magnitude. The predictions of the redshift distribution of $22.5<b_j<24.0$ are also presented for comparison with future data.