If neutrinos are Dirac particles the existence of light right-handed neutrinos $\nu_{R}$ is implied. Those would contribute to the effective number of relativistic neutrino species $N_{{\rm eff}}$ in the early Universe. With pure standard model interactions, the contribution is negligibly small. In the presence of new interactions, however, the contribution could be significantly enhanced. We consider the most general effective four-fermion interactions for neutrinos (scalar, pseudo-scalar, vector, axial-vector and tensor), and compute the contribution of right-handed neutrinos to $N_{{\rm eff}}$. Taking the Planck 2018 measurement of $N_{{\rm eff}}$, strong constraints on the effective four-fermion coupling are obtained, corresponding to interaction strengths of $10^{-5}\sim10^{-3}$ in units of the Fermi constant. This translates in new physics scales of up to 43 TeV and higher. Future experiments such as CMB-S4 can probe or exclude the existence of effective 4-neutrino operators for Dirac neutrinos. Ways to avoid this conclusion are discussed.

In this work, we explore the cosmological consequences of the latest Type Ia supernova (SN Ia) data-set, Pantheon, by adopting the $wCDM$ model. The Pantheon data-set is the largest SN Ia samples till now, which contains 1048 supernovae on the redshift range $0 < z < 2.3$. Here we take into account three kinds of SN Ia statistics techniques, including: 1. magnitude statistics (MS), which is the traditional SN Ia statistics technique; 2. flux statistics (FS), which bases on the flux-averaging (FA) method; 3. improved flux statistics (IFS), which combines the advantages of MS and FS. It should be mentioned that, The IFS technique need to scan the $(z_{cut},\Delta z)$ parameters plane, where $z_{cut}$ and $\Delta z$ are redshift cut-off and redshift interval of FA, respectively. The results are shown as follows. (1) Using SN data-set only, the best FA recipe for IFS is $(z_{cut},\Delta z)=(0.1,0.08)$; (2) Comparing to the old SN data-set, JLA, adopting Pantheon data-set can reduce the $2\sigma$ error bars of equation of state $w$ by 38\%, 47\% and 53\% for MS, FS and IFS, respectively; (3) FS gives closer results to other observations, such as Baryon acoustic oscillations and Cosmic microwave background; (4) Compared with FS and IFS, MS more favors a Universe that will end in a "big rip".

The inconsistent Hubble constant values derived from cosmic microwave background (CMB) observations and from local distance-ladder measurements may suggest new physics beyond the standard $\Lambda$CDM paradigm. It has been found in earlier works that, at least phenomenologically, non-standard recombination histories can reduce the $\gtrsim 4\sigma$ Hubble tension to $\sim 2\sigma$. Following this path, we vary physical and phenomenological parameters in RECFAST, the standard code to compute ionization history of the universe, to explore possible physics beyond standard recombination. We find that the CMB constraint on the Hubble constant is sensitive to the Hydrogen ionization energy and $2s \rightarrow 1s$ two-photon decay rate, both of which are atomic constants, and is insensitive to other details of recombination. Thus, the Hubble tension is very robust against perturbations of recombination history, unless exotic physics modifies the atomic constants during the recombination epoch.

The tomographic Alcock-Paczynski (AP) method utilizes the redshift evolution of the AP distortion to place constraints on cosmological parameters. It has proved to be a robust method that can separate the AP signature from the redshift space distortion (RSD) effect, and deliver powerful cosmological constraints using the $\lesssim 40h^{-1}\ \rm Mpc$ clustering region. In previous works, the tomographic AP method was performed via the anisotropic 2-point correlation function statistic. In this work we consider the feasibility of conducting the analysis in the Fourier domain and examine the pros and cons of this approach. We use the integrated galaxy power spectrum (PS) as a function of direction, $\hat P_{\Delta k}(\mu)$, to quantify the magnitude of anisotropy in the large-scale structure clustering, and use its redshift variation to do the AP test. The method is tested on the large, high resolution Big-MultiDark Planck (BigMD) simulation at redshifts $z=0-1$, using the underlying true cosmology $\Omega_m=0.3071,\ w=-1$. Testing the redshift evolution of $\hat P_{\Delta k}(\mu)$ in the true cosmology and cosmologies deviating from the truth with $\delta \Omega_m=0.1,\ \delta w=0.3$, we find that the redshift evolution of the AP distortion overwhelms the effects created by the RSD by a factor of $\sim1.7-3.6$. We test the method in the range of $k\in(0.2,1.8)\ h\ \rm Mpc^{-1}$, and find that it works well throughout the entire regime. We tune the halo mass within the range $2\times 10^{13}$ to $10^{14}\ M_{\odot}$, and find that the change of halo bias results in $\lesssim 5 \%$ change in $\hat P_{\Delta k}(\mu)$, which is less significant compared with the cosmological effect. Our work shows that it is feasible to conduct the tomographic AP analysis in the Fourier space.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

Precision measurements of charged cosmic rays have recently been carried out by space-born (e.g. AMS-02), or ground experiments (e.g. HESS). These measured data are important for the studies of astro-physical phenomena, including supernova remnants, cosmic ray propagation, solar physics and dark matter. Those scenarios usually contain a number of free parameters that need to be adjusted by observed data. Some techniques, such as Markov Chain Monte Carlo and MultiNest, are developed in order to solve the above problem. However, it is usually required a computing farm to apply those tools. In this paper, a genetic algorithm for finding the optimum parameters for cosmic ray injection and propagation is presented. We find that this algorithm gives us the same best fit results as the Markov Chain Monte Carlo but consuming less computing power by nearly 2 orders of magnitudes.

The tomographic Alcock-Paczynski (AP) method can result in tight cosmological constraints by using small and intermediate clustering scales of the large scale structure (LSS) of the galaxy distribution. By focusing on the redshift dependence, the AP distortion can be distinguished from the distortions produced by the redshift space distortions (RSD). In this work, we combine the tomographic AP method with other recent observational datasets of SNIa+BAO+CMB+$H_0$ to reconstruct the dark energy equation-of-state $w$ in a non-parametric form. The result favors a dynamical DE at $z\lesssim1$, and shows a mild deviation ($\lesssim2\sigma$) from $w=-1$ at $z=0.5-0.7$. We find the addition of the AP method improves the low redshift ($z\lesssim0.7$) constraint by $\sim50\%$.

The tomographic Alcock-Paczynski (AP) test is a robust large-scale structure (LSS) measurement that receives little contamination from the redshift space distortion (RSD). It has placed tight cosmological constraints by using small and intermediate clustering scales of the LSS data. However, previous works have neglected the cross-correlation among different redshift bins, which could cause the statistical uncertainty being underestimated by $\sim$20\%. In this work, we further improve this method by including this multi-redshifts full correlation. We apply it to the SDSS DR12 galaxies sample and find out that, for $\Lambda$CDM, the combination of AP with the Planck+BAO dataset slightly reduces (within 1-$\sigma$) $\Omega_m$ to $0.304\pm0.007$ (68.3\% CL). This then leads to a larger $H_0$ and also mildly affects $\Omega_b h^2$, $n_s$ and the derived parameters $z_*$, $r_*$, $z_{re}$ but not $\tau$, $A_s$ and $\sigma_8$. For the flat $w$CDM model, our measurement gives $\Omega_m=0.301\pm 0.010$ and $w=-1.090\pm 0.047$, where the additional AP measurement reduces the error budget by $\sim 25\%$. When including more parameters into the analysis, the AP method also improves the constraints on $\Omega_k$, $\sum m_\mu$, $N_{\rm eff}$ by $20-30\%$. Early universe parameters such as $dn_s/d{\rm ln}k$ and $r$, however, are unaffected. Assuming the dark energy equation of state $w=w_0+w_a \frac{z}{1+z}$, the Planck+BAO+SNIa+$H_0$+AP datasets prefer a dynamical dark energy at $\approx1.5 \sigma$ CL. Finally, we forecast the cosmological constraints expected from the DESI galaxy survey and find that combining AP with CMB+BAO method would improve the $w_0$-$w_a$ constraint by a factor of $\sim 10$.

In this work, we explore the evolution of the dark energy equation of state {\omega} by using Chevalliear-PolarskiLinder (CPL) parametrization and the binned parametrizations. For binned parametrizations, we adopt three methods to choose the redshift interval: I. Ensure that "{\Delta}z = const", where {\Delta}z is the width of each bin; II.Ensure that "n{\Delta}z = const", where n is the number of SNIa in each bin; III. Treat redshift discontinuity points as models parameters, i.e. "free {\Delta}z". For observational data, we adopt JLA type Ia supernova (SNIa) samples, SDSS DR12 data, and Planck 2015 distance priors. In particular, for JLA SNIa samples, we consider three statistic techniques: I. Magnitude statistics, which is the traditional method; II. Flux statistics, which reduces the systematic uncertainties of SNIa; III. Improve flux statistics, which can reduce the systematic uncertainties and give tighter constrains at the same time. The results are as follows: (1) For all the cases, {\omega} = -1 is always satisfied at 1{\sigma} confidence regions; It means that {\Lambda}CDM is still favored by current observations. (2) For magnitude statistics, "free {\Delta}z" model will give the smallest error bars; this conclusion does not hold true for flux statistics and improved flux statistic. (3) The improved flux statistic yields a largest present fractional density of matter {\Omega}m; in addition, this technique will give a largest current deceleration parameter q0 , which reveals a universe with a slowest cosmic acceleration.

We make a comparison for thirteen dark energy (DE) models by using current cosmological observations, including type Ia supernova, baryon acoustic oscillations, and cosmic microwave background. To perform a systematic and comprehensive analysis, we consider three statistics methods of SNIa, including magnitude statistic (MS), flux statistic (FS), and improved flux statistic (IFS), as well as two kinds of BAO data. In addition, Akaike information criteria (AIC) and Bayesian information criteria (BIC) are used to assess the worth of each model. We find that: (1) The thirteen models can be divided into four grades by performing cosmology-fits. The cosmological constant model, which is most favored by current observations, belongs to grade one; $\alpha$DE, constant $w$ and generalized Chaplygin gas models belong to grade two; Chevalliear-Polarski-Linder (CPL) parametrization, Wang parametrization, doubly coupled massive gravity, new generalized Chaplygin gas and holographic DE models belong to grade three; agegraphic DE, Dvali-Gabadadze-Porrati, Vacuum metamorphosis and Ricci DE models, which are excluded by current observations, belong to grade four. (2) For parameter estimation, adopting IFS yields the biggest $\Omega_m$ and the smallest $h$ for all the models. In contrast, using different BAO data does not cause significant effects. (3) IFS has the strongest constraint ability on various DE models. For examples, adopting IFS yields the smallest value of $\Delta$AIC for all the models; in addition, making use of this technique yields the biggest figure of merit for CPL and Wang parametrizations.

The origin of the progenitors of type Ia supernovae (SNe Ia) is still uncertain. The core-degenerate (CD) scenario has been proposed as an alternative way for the production of SNe Ia. In this scenario, SNe Ia are formed at the final stage of common-envelope evolution from a merger of a carbon-oxygen white dwarf (CO WD) with the CO core of an asymptotic giant branch companion. However, the birthrates of SNe Ia from this scenario are still not well determined. In this work, we performed a detailed investigation on the CD scenario based on a binary population synthesis approach. The SN Ia delay times from this scenario are basically in the range of 90Myr-2500Myr, mainly contributing to the observed SNe Ia with short and intermediate delay times although this scenario can also produce some old SNe Ia. Meanwhile, our work indicates that the Galactic birthrates of SNe Ia from this scenario are no more than 20% of total SNe Ia due to more careful treatment of mass transfer. Although the SN Ia birthrates in the present work are lower than those in Ilkov & Soker, the CD scenario cannot be ruled out as a viable mechanism for the formation of SNe Ia. Especially, SNe Ia with circumstellar material from this scenario contribute to 0.7-10% of total SNe Ia, which means that the CD scenario can reproduce the observed birthrates of SNe Ia like PTF 11kx. We also found that SNe Ia happen systemically earlier for a high value of metallicity and their birthrates increase with metallicity.

The observation of $1.97\pm0.04$ solar-mass neutron-like star gives constraint on the equation of state (EOS) of cold, condensed matter. In this paper, the EOS for both pure quark star and hybrid star with a quark core described by quasi-particle model are considered. The parameters of quasi-particle model which affect the mass of both quark star and hybrid star can be constrained by the observation.

In this paper, incorporating the property of the vacuum negative pressure, namely, the bag constant, we presented a new model of the equation of state (EOS) of quark matter at finite chemical potential and zero temperature. By comparing our EOS with Fraga {\it et~al.}'s EOS and SQM1 model, one find that our EOS is softer than Fraga {\it et~al.}'s EOS and SQM1 model. The reason for this difference is analyzed. With these results we investigate the structure of quark star. A comparison between our model of quark star and other models is made. The obtained mass of quark star is $ 1.3 \sim 1.66 M_\odot$ and the radius is $9.5 \sim 14 Km$. One can see that our star's compactness is smaller than that of other two models.

The temperature of hot accretion flows around black holes is sufficiently high for the ignition of nuclear reactions. This is potentially an important nucleosynthesis mechanism in the universe. As the first step in studying this problem, we need to measure physical quantities such as density and temperature of the accretion flow. In usual studies of the hot accretion flow, viscous dissipation is considered to be the only heating mechanism, while the heating caused by nuclear reactions is not considered. In this paper, we investigate whether the energy generation rate of nuclear reaction is important compared to the viscous heating. Our calculation indicates that the former is at most one percent of the latter and thus is not important. The dynamics of accretion flow can be therefore calculated in the usual way, without the need to consider heating due to nuclear reactions.

ASTROD I is a planned interplanetary space mission with multiple goals. The primary aims are: to test General Relativity with an improvement in sensitivity of over 3 orders of magnitude, improving our understanding of gravity and aiding the development of a new quantum gravity theory; to measure key solar system parameters with increased accuracy, advancing solar physics and our knowledge of the solar system and to measure the time rate of change of the gravitational constant with an order of magnitude improvement and the anomalous Pioneer acceleration, thereby probing dark matter and dark energy gravitationally. It is an international project, with major contributions from Europe and China and is envisaged as the first in a series of ASTROD missions. ASTROD I will consist of one spacecraft carrying a telescope, four lasers, two event timers and a clock. Two-way, two-wavelength laser pulse ranging will be used between the spacecraft in a solar orbit and deep space laser stations on Earth, to achieve the ASTROD I goals. A second mission, ASTROD II is envisaged as a three-spacecraft mission which would test General Relativity to one part per billion, enable detection of solar g-modes, measure the solar Lense-Thirring effect to 10 parts per million, and probe gravitational waves at frequencies below the LISA bandwidth. In the third phase (ASTROD III or Super-ASTROD), larger orbits could be implemented to map the outer solar system and to probe primordial gravitational-waves at frequencies below the ASTROD II bandwidth.

We present the discovery of a peculiar blue core in the elliptical galaxy IC 225 by using images and spectrum from the Sloan Digital Sky Survey (SDSS). The outer parts of the surface brightness profiles of u-, g-, r-, i- and z-band SDSS images for IC 225 are well fitted with an exponential function. The fitting results show that IC 225 follows the same relations between the magnitude, scale length and central surface brightness for dwarf elliptical galaxies. Its absolute blue magnitude (M_B) is -17.14 mag, all of which suggest that IC 225 is a typical dwarf elliptical galaxy. The g-r color profile indicates a very blue core with a radius of 2 arcseconds, which is also clearly seen in the RGB image made of g-, r- and i-band SDSS images. The SDSS optical spectrum exhibits strong and very narrow nebular emission lines. The metal abundances derived by the standard methods, which are 12+log(O/H) = 8.98, log(N/O) = -0.77 and 12+log(S+/H+) = 6.76, turn out to be significantly higher than that predicted by the well-known luminosity-metallicity relation. After carefully inspecting the central region of IC 225, we find that there are two distinct nuclei, separated by 1.4 arcseconds, the off-nucleated one is even bluer than the nucleus of IC 225. The asymmetric line profiles of higher-order Balmer lines indicate that the emission lines are bluer shifted relative to the absorption lines, suggesting that the line emission arises from the off-center core, whose nature is a metal-rich Hii region. To the best of our knowledge, it is the first high-metallicity Hii region detected in a dwarf elliptical galaxy.

We present the results of 11 elliptical galaxies with strong nebular emission lines during our study of star formation history along the Hubble sequence. After removing the dilution from the underlying old stellar populations by use of stellar population synthesis model, we derive the accurate fluxes of all emission lines for these objects, which are later classified with emission line ratios into one Seyfert 2, six LINERs and four HII galaxies. We also identify one HII galaxy (A1216+04) as a hitherto unknown Wolf-Rayet galaxy from the presence of the Wolf-Rayet broad bump at 4650 \AA. We propose that the star-forming activities in elliptical galaxies are triggered by either galaxy-galaxy interaction or the merging of a small satellite/a massive star cluster, as already suggested by recent numerical simulations.

On the basis of the neutrino emission from the isotropic $^1$S$_0$ neutron superfluid vortexes in neutron star interiors, we propose a rocket model of neutrino jet for the observed pulsar kick.

We show that the baryonic gas content of the halo of our galaxy can be probed by performing a multipole expansion on the distribution of diffuse background gamma-ray emission. While the monopole moment (isotropic background) can be used to constrain the baryonic fraction of the halo gas, the quadrupole to monopole ratio is a sensitive probe of the distribution of gas in the halo, {\it i.e.} the degree of flattening of the gas distribution. The predicted diffuse gamma ray flux is found to be very sensitive to the adopted cosmic ray density distribution throughout the halo. If the cosmic rays are uniform, then the upper bound on the gas fraction is $16.6\%$ regardless of the flattening of the halo gas distribution. However, this bound can be weakened by taking into account the removal of flux in and close to the galactic plane, especially for a oblate ($e <1$) gas distribution. On the other hand, in the more realistic situation that the cosmic rays linearly trace the smoothed halo mass (and halo gas) distribution, then a stringent bound on the baryonic gas fraction in the halo $\eta$, $\eta \ltwid 3\%$, can be placed with existing data, regardless of halo flattening.

Non-Gaussian imprints on the cosmic microwave background radiation (CBR) sky are within the grasp of current experiments. A clear non-Gaussian signature would be point-like sources. We have examined the nature of possible point sources that were tentatively identified in a recent high frequency CBR experiment with half-degree resolution (Cheng et al. 1994: MSAM). The effects of local foreground sources, including cold dust clouds, radio sources and the Sunyaev-Zeldovich effect due to foreground rich clusters are considered, and the effective spectral slopes of these various foreground contaminations are calculated. Radio source emission and the Sunyaev-Zeldovich effect are ruled out as the explanation of the possible MSAM sources. Modeles are examined of extremely cold dust clouds which are located in the solar neighborhood, the interstellar medium, the galactic halo or at cosmological distances. We estimate the cloud mass and dust grain parameters, and in particular the grain size distribution, that are required in order to produce a detectable signal in an MSAM-type experiment. It is shown that cold dust clouds can have an important effect on CBR experiments only if the clouds are nearby, and located within a disc scale height of the solar neighborhood. Cold dust emission remains a possible source of far-infrared signal at the level of the detected CBR fluctuations on degree scales, but may be ruled out if the dust emissivity index satisfies $\alpha = 1.5 \pm 0.5$ .

Non-Gaussian distributions of cosmic microwave background (CMB) anisotropies have been proposed to reconcile the discrepancies between different experiments at half-degree scales (Coulson et al. 1994). Each experiment probes a different part of the sky, furthermore, sky coverage is very small, hence the sample variance of each experiment can be large, especially when the sky signal is non-Gaussian. We model the degree-scale CMB sky as a $\chi_{n}^{2}$ field with $n$-degrees of freedom and show that the sample variance is enhanced over that of a Gaussian distribution by a factor of $ {(n +6)/ n}$. The sample variance for different experiments are calculated, both for Gaussian and non-Gaussian distributions. We also show that if the distribution is highly non-Gaussian $(n \ltwid 4)$ at half-degree scales, then the non-Gaussian signature of the CMB could be detected in the FIRS map, though probably not in the COBE map.

COBE has provided us with a whole-sky map of the CBR anisotropies. However, even if the noise level is negligible when the four year COBE data are available, the cosmic variance will prevent us from obtaining information about the Gaussian nature of the primordial fluctuations. This important issue is addressed here by studying the angular bispectrum of the cosmic microwave background anisotropies. A general form of the angular bispectrum is given and the cosmic variance of the angular bispectrum for Gaussian fluctuations is calculated. The advantage of using the angular bispectrum is that one can choose to use the multipole moments which minimize the cosmic variance term. The non-Gaussian signals in most physically motivated non-Gaussian models are small compared with cosmic variance. Unless the amplitudes are large, the non-Gaussian signals are only detectable in the COBE data in those models where the angular bispectrum is flat or increases with increasing multipole moment.

Information about the physical processes that generate the primordial fluctuations in the early universe can be gained by testing the Gaussian nature of the fluctuations through cosmic microwave background radiation (CBR) temperature anisotropy experiments. One of the crucial aspects of density perturbations that are produced by the standard inflation scenario is that they are Gaussian, whereas seeds produced by topological defects left over from an early cosmic phase transition tend to be non-Gaussian. To carry out this test, sophisticated statistical tools are required. In this paper, we will discuss several such statistical tools, including multivariant skewness and kurtosis, Euler-Poincare characteristics, the three point temperature correlation function, and the Hotelling's $T^{2}$ statistic defined through bispectral estimates of a one dimensional dataset. The effect of noise present in the current data is discussed in detail and the COBE 53 GHz dataset is analyzed. Our analysis shows that, on the large angular scale to which COBE is sensitive, the statistics are probably Gaussian. On the small angular scales, the importance of Hotelling's $T^{2}$ statistic is stressed, and the minimum sample size required to test Gaussianity is estimated. Although the current dataset available from various experiments at half-degree scales is still too small, improvement of the dataset by roughly a factor of two will be enough to test the Gaussianity statistically. On the arcminute scales, we analyze the recent RING data through bispectral analysis, and the result

One of the crucial aspects of density perturbations that are produced by the standard inflation scenario is that they are Gaussian where seeds produced by topological defects tend to be non-Gaussian. The three point correlation function of the temperature anisotropy of the cosmic microwave background radiation (CBR) provides a sensitive test of this aspect of the primordial density field. In this paper, this function is calculated in the general context of various allowed non-Gaussian models. It is shown that by COBE and the forthcoming South Pole and Balloon CBR anisotropy data may be able to test provide a crucial test of Gaussianity.