Accurate background estimation is essential for spectral and temporal analysis in astrophysics. In this work, we construct the in-orbit background model for the High-Energy Telescope (HE) of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). Based on the two-year blank sky observations of Insight-HXMT/HE, we first investigate the basic properties of the background and find that both the background spectral shape and intensity have long-term evolution at different geographical sites. The entire earth globe is then divided into small grids, each with a typical area of 5x5 square degrees in geographical coordinate system. For each grid, an empirical function is used to describe the long-term evolution of each channel of the background spectrum; the intensity of the background can be variable and a modification factor is introduced to account for this variability by measuring the contemporary flux of the blind detector. For a given pointing observation, the background model is accomplished by integrating over the grids that are passed by the track of the satellite in each orbit. Such a background model is tested with both the blank sky observations and campaigns for observations of a series of celestial sources. The results show an average systematic error of 1.5% for the background energy spectrum (26-100 keV) under a typical exposure of 8 ks, and <3% for background light curve estimation (30-150 keV). Therefore, the background model introduced in this paper is included in the Insight-HXMT software as a standard part specialized for both spectral and temporal analyses.

Observations recently suggest $4.4\sigma$ discrepancy between Hubble constant measured locally and that inferred from cosmic microwave background in standard cosmology. Either standard cosmology or local measurement might expect something new. We investigate the possibility that the value of Hubble constants might be affected by observers' motional status in the local measurement. Using adapted coordinate for geodesic observers in FLRW space-time and constraints inferred from observation of cosmic shear, we find that the motional status of reference frame could contribute to about $1.1 \pm 0.3\%$ discrepancy of these Hubble constants. Namely, the difference is $\Delta H_0= 0.81\pm0.22$ km s$^{-1}$ Mpc$^{-1}$, if the locally measured Hubble constant is $ 74$ km s$^{-1}$ Mpc$^{-1}$. This effect seems not negligible as an uncertainty for local measurement of Hubble constant

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR~J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR~J1935+2154. The burst showed two hard peaks with a separation of $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.57$ s, fully consistent with the dispersion delay of FRB~200428. We thus identify the non-thermal X-ray burst is associated with FRB~200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB~200428 share the same physical origin in an explosive event from SGR~J1935+2154.

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

We study the evolution of the temporal properties of MAXI 1820+070 during the 2018 outburst in its hard state from MJD 58190 to 58289 with Insight-HXMT in a broad energy band 1-150 keV. We find different behaviors of the hardness ratio, the fractional rms and time lag before and after MJD 58257, suggesting a transition occurred at around this point. The observed time lags between the soft photons in the 1-5 keV band and the hard photons in higher energy bands, up to 150 keV, are frequency-dependent: the time lags in the low-frequency range, 2-10 mHz, are both soft and hard lags with a timescale of dozens of seconds but without a clear trend along the outburst; the time lags in the high-frequency range, 1-10 Hz, are only hard lags with a timescale of tens of milliseconds; first increase until around MJD 58257 and decrease after this date. The high-frequency time lags are significantly correlated to the photon index derived from the fit to the quasi-simultaneous NICER spectrum in the 1-10 keV band. This result is qualitatively consistent with a model in which the high-frequency time lags are produced by Comptonization in a jet.

With more than 150 blank sky observations at high Galactic latitude, we make a systematic study to the background of the Low Energy Telescope (LE) of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). Both the on-ground simulation and the in-orbit observation indicate that the background spectrum mainly has two components. One is the particle background that dominates above 7 keV and its spectral shape is consistent in every geographical locations. Another is the diffuse X-ray background that dominates below 7 keV and has a stable spectrum less dependent of the sky region. The particle background spectral shape can be obtained from the blind detector data of all the blank sky observations, and the particle background intensity can be measured by the blind detector at 10-12.5 keV. The diffuse X-ray background in the high Galactic latitude can also be obtained from the blank sky spectra after subtracting the particle background. Based on these characteristics, we develop the background model for both the spectrum and the light curve. The systematic error for the background spectrum is investigated with different exposures (T_exp). For the spectrum with T_exp=1 ks, the average systematic errors in 1-7 keV and 1-10 keV are 4.2% and 3.7%, respectively. We also perform the systematic error analyses of the background light curves with different energy bands and time bins. The results show that the systematic errors for the light curves with different time bins are <8% in 1-10 keV.

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches.

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB.

The Medium Energy X-ray Telescope (ME) is one of the main payloads of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). The background of Insight-HXMT/ME is mainly caused by the environmental charged particles and the background intensity is modulated remarkably by the geomagnetic field, as well as the geographical location. At the same geographical location, the background spectral shape is stable but the intensity varies with the level of the environmental charged particles. In this paper, we develop a model to estimate the ME background based on the ME database that is established with the two-year blank sky observations of the high Galactic latitude. In this model, the entire geographical area covered by Insight-HXMT is divided into grids of $5^{\circ}\times5^{\circ}$ in geographical coordinate system. For each grid, the background spectral shape can be obtained from the background database and the intensity can be corrected by the contemporary count rate of the blind FOV detectors. Thus the background spectrum can be obtained by accumulating the background of all the grids passed by Insight-HXMT during the effective observational time. The model test with the blank sky observations shows that the systematic error of the background estimation in $8.9-44.0$ keV is $\sim1.3\%$ for a pointing observation with an average exposure $\sim5.5$ ks. We also find that the systematic error is anti-correlated with the exposure, which indicates the systematic error is partly contributed by the statistical error of count rate measured by the blind FOV detectors.

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate.

Chang'E 4 is the first mission to the far side of the Moon and consists of a lander, a rover, and a relay spacecraft. Lander and rover were launched at 18:23 UTC on December 7, 2018 and landed in the von K\'arm\'an crater at 02:26 UTC on January 3, 2019. Here we describe the Lunar Lander Neutron \& Dosimetry experiment (LND) which is part of the Chang'E 4 Lander scientific payload. Its chief scientific goal is to obtain first active dosimetric measurements on the surface of the Moon. LND also provides observations of fast neutrons which are a result of the interaction of high-energy particle radiation with the lunar regolith and of their thermalized counterpart, thermal neutrons, which are a sensitive indicator of subsurface water content.

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.

LS 5039 is a well-known $\gamma$-ray binary system which consists of an unknown compact object and a massive companion O star. It shows rather stable emissions at high energies over years and hence serves as an ideal laboratory to investigate the emission mechanism for such peculiar systems which emit prominent $\gamma$-rays. To this end, we take the orbital phase resolved energy spectrum as observed by \fermi over 10 years. We divide the orbit into four orbital phases, each with an orbital phase range of 0.25, centered at 0.00, 0.25, 0.50 and 0.75 respectively, where the phase 0.0 is the periastron and phase 0.5 is the apastron. The phases around 0.25 and 0.75 are symmetric and hence are supposed to have identical local acceleration environment. The spectral analysis shows that, the \fermi spectra are largely different from these two symmetric orbital phases: the emission from orbital phase 0.25 turns out to be significantly stronger than that from 0.75. This result does not fit a scenario that $\gamma$-rays are Doppler boosted emission from bow shock tails if LS 5039 has a shock configuration similar to PSR B1259-63, and indicates that the inverse Compton scatterings between the shock accelerated plasma and the stellar particle environment is the underline procedure. We also find that the previous report for a disappearance of the orbital modulation at 3--20 GeV is due to the similar spectral turn-over energies of the different orbital phases. The spectral properties of periastron and apastron regions are addressed in the context of the measurements in phase regions around 0.25 and 0.75.

We study the 2018 outburst of Aql X-1 via the monitor of all sky X-ray image (MAXI) data. We show that the outburst starting in February 2018 is a member of short-low class in the frame of outburst duration and the peak count rate although the outburst morphology is slightly different from the other fast-rise-exponential-decay (FRED) type outbursts with a milder rising stage. We study the partial accretion in the weak propeller stage of Aql X-1 via the MAXI data of the 2018 outburst. We report on the spectral analysis of 3 observations of Aquila X-1 obtained by Insight - hard X-ray modulation telescope (Insight-HXMT) during the late decay stage of the 2018 outburst. We discuss that the data taken by Insight-HXMT is just after the transition to the weak propeller stage. Our analysis shows the necessity of a comptonization component to take into account the existence of an electron cloud resulting photons partly up-scattered.

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

We present a detailed timing study of the brightest persistent X-ray source Sco X-1 using the data collected by the Hard X-ray Modulation Telescope ($Insight$-HXMT) from July 2017 to August 2018. A complete $Z$-track hardness-intensity diagram (HID) is obtained. The normal branch oscillations (NBOs) at $\sim$ 6 Hz in the lower part of the normal branch (NB) and the flare branch oscillations (FBOs) at $\sim$ 16 Hz in the beginning part of the flaring branch (FB) are found in observations with the Low Energy X-ray Telescope (LE) and the Medium Energy X-ray Telescope (ME) of $Insight$-HXMT, while the horizontal branch oscillations (HBOs) at $\sim$ 40 Hz and the kilohertz quasi-periodic oscillations (kHz QPOs) at $\sim$ 800 Hz are found simultaneously up to 60 keV for the first time on the horizontal branch (HB) by the High Energy X-ray Telescope (HE) and ME. We find that for all types of the observed QPOs, the centroid frequencies are independent of energy, while the root mean square (rms) increases with energy; the centroid frequencies of both the HBOs and kHz QPOs increase along the $Z$-track from the top to the bottom of the HB; and the NBOs show soft phase lags increasing with energy. A continuous QPO transition from the FB to NB in $\sim$ 200 s are also detected. Our results indicate that the non-thermal emission is the origin of all types of QPOs, the innermost region of the accretion disk is non-thermal in nature, and the corona is nonhomogeneous geometrically.

The evidences for the influence of thermonuclear (type-I) X-ray bursts upon the surrounding environments in neutron star low-mass X-ray binaries (LMXB) were detected previously via spectral and timing analyses. Benefitting from a broad energy coverage of Insight-HXMT, we analyze one photospheric radius expansion (PRE) burst, and find an emission excess at soft X-rays. Our spectral analysis shows that, such an excess is not likely relevant to the disk reflection induced by the burst emission and can be attributed to an enhanced pre-burst/persistent emission. We find that the burst and enhanced persistent emissions sum up to exceed Eddington luminosity by $\sim$ 40 percentages. We speculate that the enhanced emission is from a region beyond the PRE radius, or through the Comptonization of the corona.

We test the possible dipole anisotropy of a Finslerian cosmological model and other three dipole-modulated cosmological models, i.e., the dipole-modulated $\rm{\Lambda}$CDM, $w$CDM and Chevallier--Polarski--Linder (CPL) model by using the recently released Pantheon sample of SNe Ia. The Markov chain Monte Carlo (MCMC) method is used to explore the whole parameter space. We find that the dipole anisotropy is very weak in all cosmological models used. Although the dipole amplitudes of four cosmological models are consistent with zero within $1\sigma$ uncertainty, the dipole directions are close to the axial direction to the plane of SDSS subsample among Pantheon. It may imply that the weak dipole anisotropy in the Pantheon sample originates from the inhomogeneous distribution of the SDSS subsample. More homogeneous distribution of SNe Ia is necessary to constrain the cosmic anisotropy.

Corona cooling was detected previously from stacking a series of short type-I bursts occurred during the low/had state of atoll outburst. Type-I bursts are hence regarded as sharp probe to our better understanding on the basic property of the corona. The launch of the first Chinese X-ray satellite Insight-HXMT has large detection area at hard X-rays which provide almost unique chance to move further in this research field. We report the first detection of the corona cooling by Insight-HXMT from single short type-I burst showing up during {\bf flare} of 4U 1636-536. This type-I X-ray burst has a duration of $\sim$13 seconds and hard X-ray shortage is detected with significance 6.2~$\sigma$ in 40-70 keV. A cross-correlation analysis between the lightcurves of soft and hard X-ray band, shows that the corona shortage lag the burst emission by 1.6 $\pm$1.2~s. These results are consistent with those derived previously from stacking a large amount of bursts detected by RXTE/PCA within a series of {\bf flares} of 4U 1636-536. Moreover, the broad bandwidth of Insight-HXMT allows as well for the first time to infer the burst influence upon the continuum spectrum via performing the spectral fitting of the burst, which ends up with the finding that hard X-ray shortage appears at around 40 keV in the continuum spectrum. These results suggest that the evolution of the corona along with the outburst{\bf /flare} of NS XRB may be traced via looking into a series of embedded type-I bursts by using Insight-HXMT.

The Hard X-ray Modulation Telescope (HXMT) named Insight is China's first X-ray astronomical satellite. The Low Energy X-ray Telescope (LE) is one of its main payloads onboard. The detectors of LE adopt swept charge device CCD236 with L-shaped transfer electrodes. Charges in detection area are read out continuously along specific paths, which leads to a time response distribution of photons readout time. We designed a long exposure readout mode to measure the time response distribution. In this mode, CCD236 firstly performs exposure without readout, then all charges generated in preceding exposure phase are read out completely. Through analysis of the photons readout time in this mode, we obtained the probability distribution of photons readout time.

The Insight-Hard X-ray Modulation Telescope (Insight-HXMT) is a broad band X-ray and gamma-ray (1-3000 keV) astronomy satellite. The High Energy X-ray telescope (HE) is one of its three main telescopes. The main detector plane of HE is composed of 18 NaI(Tl)/CsI(Na) phoswich detectors, where NaI(Tl) serves as primary detector to measure ~ 20-250 keV photons incident from the field of view (FOV) defined by the collimators, and CsI(Na) is used as an active shield detector to NaI(Tl) by pulse shape discrimination. CsI(Na) is also used as an omnidirectional gamma-ray monitor. The HE collimators have a diverse FOV: 1.1{\deg}x 5.7{\deg} (15 units), 5.7{\deg}x 5.7{\deg} (2 units) and blocked (1 unit), thus the combined FOV of HE is about 5.7{\deg}x 5.7{\deg}. Each HE detector has a diameter of 190 mm, resulting in the total geometrical area of about 5100 cm_2. The energy resolution is ~15% at 60 keV. The timing accuracy is better than 10 {\mu}s and dead-time for each detector is less than 10 {\mu}s. HE is devoted to observe the spectra and temporal variability of X-ray sources in the 20-250 keV band either by pointing observations for known sources or scanning observations to unveil new sources, and to monitor the gamma-ray sky in 0.2-3 MeV. This paper presents the design and performance of the HE instruments. Results of the on-ground calibration experiments are also reported.

We present a study of timing properties of the accreting pulsar 2S 1417-624 observed during its 2018 outburst, based on Swift/BAT, Fermi/GBM, Insight-HXMT and NICER observations. We report a dramatic change of the pulse profiles with luminosity. The morphology of the profile in the range 0.2-10.0keV switches from double to triple peaks at $\sim2.5$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$ and from triple to quadruple peaks at $\sim7$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$. The profile at high energies (25-100keV) shows significant evolutions as well. We explain this phenomenon according to existing theoretical models. We argue that the first change is related to the transition from the sub to the super-critical accretion regime, while the second to the transition of the accretion disc from the gas-dominated to the radiation pressure-dominated state. Considering the spin-up as well due to the accretion torque, this interpretation allows to estimate the magnetic field self-consistently at $\sim7\times 10^{12}$G.

The long-term evolution of the centroid energy of the CRSF in Her X-1 is still a mystery. We report a new measurement from a campaign between {\sl Insight}-HXMT and {\sl NuSTAR} performed in February 2018. Generally, the two satellites show well consistent results of timing and spectral properties. The joint spectral analysis confirms that the previously observed long decay phase has ended, and that the line energy instead keeps constant around 37.5 keV after flux correction.

We report on analysis of observations of the bright transient X-ray pulsar \src obtained during its 2017-2018 giant outburst with Insight-HXMT, \emph{NuSTAR}, and \textit{Swift} observatories. We focus on the discovery of a sharp state transition of the timing and spectral properties of the source at super-Eddington accretion rates, which we associate with the transition of the accretion disk to a radiation pressure dominated (RPD) state, the first ever directly observed for magnetized neutron star. This transition occurs at slightly higher luminosity compared to already reported transition of the source from sub- to super-critical accretion regime associate with onset of an accretion column. We argue that this scenario can only be realized for comparatively weakly magnetized neutron star, not dissimilar to other ultra-luminous X-ray pulsars (ULPs), which accrete at similar rates. Further evidence for this conclusion is provided by the non-detection of the transition to the propeller state in quiescence which strongly implies compact magnetosphere and thus rules out magnetar-like fields.

In this work, we report the in-orbit demonstration of X-ray pulsar navigation with Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which was launched on Jun. 15th, 2017. The new pulsar navigation method 'Significance Enhancement of Pulse-profile with Orbit-dynamics' (SEPO) is adopted to determine the orbit with observations of only one pulsar. In this test, the Crab pulsar is chosen and observed by Insight-HXMT from Aug. 31th to Sept. 5th in 2017. Using the 5-day-long observation data, the orbit of Insight-HXMT is determined successfully with the three telescopes onboard - High Energy X-ray Telescope (HE), Medium Energy X-ray Telescope (ME) and Low Energy X-ray Telescope (LE) - respectively. Combining all the data, the position and velocity of the Insight-HXMT are pinpointed to within 10 km (3 sigma) and 10 m/s (3 sigma), respectively.

The recently discovered neutron star transient Swift J0243.6+6124 has been monitored by {\it the Hard X-ray Modulation Telescope} ({\it Insight-\rm HXMT). Based on the obtained data, we investigate the broadband spectrum of the source throughout the outburst. We estimate the broadband flux of the source and search for possible cyclotron line in the broadband spectrum. No evidence of line-like features is, however, found up to $\rm 150~keV$. In the absence of any cyclotron line in its energy spectrum, we estimate the magnetic field of the source based on the observed spin evolution of the neutron star by applying two accretion torque models. In both cases, we get consistent results with $B\rm \sim 10^{13}~G$, $D\rm \sim 6~kpc$ and peak luminosity of $\rm >10^{39}~erg~s^{-1}$ which makes the source the first Galactic ultraluminous X-ray source hosting a neutron star.

We employ the hemisphere comparison (HC) method and the dipole fitting (DF) method to investigate the cosmic anisotropy in the recently released Pantheon sample of type Ia supernovae (SNe Ia) and five combinations among Pantheon. For the HC method, we find the maximum anisotropy level in the full Pantheon sample is $\mathrm{AL}_{max}=0.361\pm0.070$ and corresponding direction $(l,b)=({123.05^{\circ}}^{+11.25^{\circ}}_{-4.22^{\circ}}, {4.78^{\circ}}^{+1.80^{\circ}}_{-8.36^{\circ}})$. A robust check shows the statistical significance of maximum anisotropy level is about $2.1\sigma$. We also find that the Low-$z$ and SNLS subsamples have decisive impact on the overall anisotropy while other three subsamples have little impact. Moreover, the anisotropy level map significantly rely on the inhomogeneous distribution of SNe Ia in the sky. For the DF method, we find the dipole anisotropy in the Pantheon sample is very weak. The dipole magnitude is constrained to be less than $1.16\times10^{-3}$ at $95\%$ confidence level. However, the dipole direction is well inferred by MCMC method and it points towards $(l,b)=({306.00^{\circ}}^{+82.95^{\circ}}_{-125.01^{\circ}}, {-34.20^{\circ}}^{+16.82^{\circ}}_{-54.93^{\circ}})$. This direction is very close to the axial direction to the plane of SDSS subsample. It may imply that SDSS subsample is the decisive part to the dipole anisotropy in the full Pantheon sample. All these facts imply that the cosmic anisotropy found in Pantheon sample significantly rely on the inhomogeneous distribution of SNe Ia in the sky. More homogeneous distribution of SNe Ia is necessary to search for a more convincing cosmic anisotropy.

Milgrom's modified Newtonian dynamics (MOND) can explain well the mass discrepancy problem in galaxy without invoking dark matter. The MOND theory predicts a universal constant acceleration scale in galaxy, below which the Newtonian dynamics is no longer valid. Recently, Rodrigues et al. analyzed 193 high-quality disk galaxies by using Bayesian inference and they concluded that the probability of existence of a fundamental acceleration scale is essentially 0. In this paper, we use the same method to revisit the same question: Is there a fundamental acceleration scale in galaxies? We fit each galaxy rotation curve of 175 SPARC galaxies with Gaussian priors on galaxy parameters and a variable acceleration scale. After marginalizing over the nuisance parameters, we find that the marginalized posterior distributions of the acceleration scale become broad. The incompatibility between the global best fit and the marginalized posterior distributions of the acceleration scale is greatly reduced. However, there still exists evidence that rejects MOND as a fundamental theory. More accurate observations are needed to exclude or confirm the existence of a fundamental acceleration scale in galaxies.

We study the long-term evolution of the centroid energy of cyclotron lines - often referered to as Cyclotron Resonance Scattering Features (CRSF) - in Her X-1, Vela X-1 and Cen X-3, using survey observations of the Burst Alert Telescope (BAT) onboard Swift. We find a significant decrease of the fundamental CRSF energy in Her X-1 and the first harmonic line energy in Vela X-1, since the launch of Swift in 2004 and until 2010 and 2012, respectively.In both sources the decreases stopped at some time, with a quite stable centroid energy thereafter. Unlike in Her X-1 and Vela X-1, the CRSF energy in Cen X-3 does not show a long-term decrease. It is observed not to change for at least the past 14 years. The long-term variation of the line energy is a direct way to investigate the magnetic field structure in the polar regions of pulsars. Our results may stimulate the development of theoretical models, especially regarding to how the accreted mass accumulates in the accretion mound or how the magnetic field distorts around the polar cap.

PSR B1259-63 is a $\gamma$-ray emitting high mass X-ray binary system, in which the compact object is a millisecond pulsar. The system has an orbital period of 1236.7 d and shows peculiar $\gamma$-ray flares when the neutron star moves out of the stellar disk of the companion star. The $\gamma$-ray flare events were firstly discovered by using Fermi-LAT around the 2010 periastron passage, which was repeated for the 2014 and 2017 periastron passages. We analyze the Fermi-LAT data for all the three periastron passages and found that in each flare the energy spectrum can be represented well by a simple power law. The $\gamma$-ray light curves show that in 2010 and 2014 after each periastron there are two main flares, but in 2017 there are four flares including one precursor about 10 d after the periastron passage. The first main flares in 2010 and 2014 are located at around 35 d after the periastron passage, and the main flare in 2014 is delayed by roughly 1.7 d with respect to that in 2010. In the 2017 flare, the source shows a precursor about 10 d after the periastron passage, but the following two flares become weaker and lag behind those in 2010 by roughly 5 d. The strongest flares in 2017 occurred 58 d and 70 d after the periastron passage. These results challenge the previous models.

We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system.

We propose a framework of Finsler space-time to explain the observed parity asymmetry and the power deficit in the low-$\ell$ ($2\leqslant \ell \leqslant 29$) multipole range of cosmic microwave background (CMB) temperature anisotropies. In the $3+1$ dimensional space-time, the three-dimensional space is described by a Randers-Finsler space, which is spatially irreversible, inducing the parity asymmetry in the CMB angular power spectrum. We estimate the constraints on the two parameters introduced by Finsler space-time via analyzing the low-$\ell$ angular power spectrum in PLANCK 2015 CMB temperature data. We see that the low-$\ell$ power suppression in the CMB temperature anisotropies can also be resolved in this scenario. Our study shows that the two low-$\ell$ anomalies, i.e., parity asymmetry and power deficit, may have a common origin.

We report a possible dipole anisotropy on acceleration scale $g_{\dag}$ with 147 rotationally supported galaxies in local Universe. It is found that a monopole and dipole correction for the radial acceleration relation can better describe the SPARC data set. The monopole term is negligible but the dipole magnitude is significant. It is also found that the dipole correction is mostly induced by the anisotropy on the acceleration scale. The magnitude of $\hat{g}_{\dag}$-dipole reaches up to $0.25\pm0.04$, and its direction is aligned to $(l,b) = (171.30^{\circ}\pm7.18^{\circ}, -15.41^{\circ}\pm4.87^{\circ})$, which is very close to the maximum anisotropy direction from the hemisphere comparison method. Furthermore, robust check shows that the dipole anisotropy couldn't be reproduced by isotropic mock data set. However, it is still premature to claim that the Universe is anisotropic due to the small data samples and uncertainty in the current observations.

We test the cosmological principle by fitting a dipolar modulation of distance modulus and searching for an evolution of this modulation with respect to cosmological redshift. Based on a redshift tomographic method, we divide the Joint Light-curve Analysis compilation of supernovae of type Ia into different redshift bins, and employ a Markov-Chain Monte-Carlo method to infer the anisotropic amplitude and direction in each redshift bin. However, we do not find any significant deviations from the cosmological principle, and the anisotropic amplitude is stringently constrained to be less than a few thousandths at $95\%$ confidence level.

In this work, we study a direction dependent power spectrum in anisotropic Finsler space-time. We use this direction dependent power spectrum to address the low-l power observed in WMAP and PLANCK data. The angular power spectrum of the temperature fluctuations has a lower amplitude in comparison to the $\Lambda$CDM model in the multipole range l = 2-40. Our theoretical model gives a correction to the isotropic angular power spectrum $C^{TT}_l$ due to the breaking of the rotational invariance of the primordial power spectrum. We estimate best-fit model parameters along with the six $\Lambda$CDM cosmological parameters using PLANCK likelihood code in CosmoMC software. We see that this modified angular power spectrum fits the CMB temperature data in the multipole range l=2-10 to a good extent but fails for the whole multipole range l = 2-40.

Finding the electromagnetic (EM) counterpart of binary compact star merger, especially the binary neutron star (BNS) merger, is critically important for gravitational wave (GW) astronomy, cosmology and fundamental physics. On Aug. 17, 2017, Advanced LIGO and \textit{Fermi}/GBM independently triggered the first BNS merger, GW170817, and its high energy EM counterpart, GRB 170817A, respectively, resulting in a global observation campaign covering gamma-ray, X-ray, UV, optical, IR, radio as well as neutrinos. The High Energy X-ray telescope (HE) onboard \textit{Insight}-HXMT (Hard X-ray Modulation Telescope) is the unique high-energy gamma-ray telescope that monitored the entire GW localization area and especially the optical counterpart (SSS17a/AT2017gfo) with very large collection area ($\sim$1000 cm$^2$) and microsecond time resolution in 0.2-5 MeV. In addition, \textit{Insight}-HXMT quickly implemented a Target of Opportunity (ToO) observation to scan the GW localization area for potential X-ray emission from the GW source. Although it did not detect any significant high energy (0.2-5 MeV) radiation from GW170817, its observation helped to confirm the unexpected weak and soft nature of GRB 170817A. Meanwhile, \textit{Insight}-HXMT/HE provides one of the most stringent constraints (~10$^{-7}$ to 10$^{-6}$ erg/cm$^2$/s) for both GRB170817A and any other possible precursor or extended emissions in 0.2-5 MeV, which help us to better understand the properties of EM radiation from this BNS merger. Therefore the observation of \textit{Insight}-HXMT constitutes an important chapter in the full context of multi-wavelength and multi-messenger observation of this historical GW event.

It is well known that the Milgrom's MOND (modified Newtonian dynamics) explains well the mass discrepancy problem in galaxy rotation curves. The MOND predicts a universal acceleration scale below which the Newtonian dynamics is invalid yet. The universal acceleration scale we got from the SPARC dataset is $g_{\dag}=1.02\times10^{-10} \rm m~s^{-2}$. Milgrom suggested that the acceleration scale may be a fingerprint of cosmology on local dynamics and related with the Hubble constant $g_{\dag}\sim cH_0$. In this paper, we use the hemisphere comparison method with the SPARC dataset to investigate the spatial anisotropy on the acceleration scale. We find that the hemisphere of the maximum acceleration scale is in the direction $(l,b) = ({175.5^\circ}^{+6^\circ}_{-10^\circ}, {-6.5^\circ}^{+8^\circ}_{-3^\circ})$ with $g_{\dag,max}=1.10\times10^{-10} \rm m~s^{-2}$, while the hemisphere of the minimum acceleration scale is in the opposite direction $(l,b) = ({355.5^\circ}^{+6^\circ}_{-10^\circ}, {6.5^\circ}^{+3^\circ}_{-8^\circ})$ with $g_{\dag,min}=0.76\times10^{-10} \rm m~s^{-2}$. The maximum anisotropy level reaches up to $0.37\pm0.04$. Robust tests present that such a level of anisotropy can't be reproduced by a statistically isotropic data. In addition, we show that the spatial anisotropy on the acceleration scale has little correlation with the non-uniform distribution of the SPARC data points in sky. We also find that the maximum anisotropy direction is close with other cosmological preferred directions, especially the direction of the "Australia dipole" for the fine structure constant.

The statistical properties of the soft gamma repeater SGR J1550--5418 are investigated carefully. We find that the cumulative distributions of fluence, peak flux and duration can be well fitted by a bent power law, while the cumulative distribution of waiting time follows a simple power law. In particular, the probability density functions of fluctuations of fluence, peak flux, and duration have a sharp peak and fat tails, which can be well fitted by a $q$-Gaussian function. The $q$ values keep approximately steady for different scale intervals, indicating a scale-invariant structure of soft gamma repeaters. Those results support that the origin of soft gamma repeaters is crustquakes of neutron stars with extremely strong magnetic fields.

Silicon Photo-Multipliers (SiPM) are regarded as novel photo-detectors to replace conventional Photo-Multiplier Tubes (PMTs). However, the breakdown voltage dependence on the ambient temperature results in a gain variation of $\sim$3$\% /^{\circ} \mathrm C$. This severely limits the application of this device in experiments with wide range of operating temperature, especially in space missions. An experimental setup was established to investigate the temperature and bias voltage dependence of gain for the Multi-Pixel Photon Counter (MPPC). The gain and breakdown voltage dependence on operating temperature of an MPPC can be approximated by a linear function, which is similar to the behavior of a zener diode. The measured temperature coefficient of the breakdown voltage is $(59.4 \pm 0.4$ mV)$/^{\circ} \mathrm C$. According to this fact, an analog power supply based on two zener diodes and an operational amplifier was designed with a positive temperature coefficient. The measured temperature dependence for the designed power supply is between 63.65 to 64.61~mV/$^{\circ} \mathrm C$ at different output voltages. The designed power supply can bias the MPPC at an over-voltage with a temperature variation of $\sim$ 5~mV$/^{\circ} \mathrm C$. The gain variation of the MPPC biased at over-voltage of 2~V was reduced from 2.8$\% /^{\circ} \mathrm C$ to 0.3$\% /^{\circ} \mathrm C$ when biased the MPPC with the designed power supply for gain control. Detailed design and performance of the analog power supply in the temperature range from -42.7$^{\circ}\mathrm{C}$ to 20.9$^{\circ}\mathrm{C}$ will be discussed in this paper.

Linear polarization have been observed in both the prompt phase and afterglow of some bright gamma-ray bursts (GRBs). Polarization in the prompt phase spans a wide range, and may be as high as $\gtrsim 50\%$. In the afterglow phase, however, it is usually below $10\%$. According to the standard fireball model, GRBs are produced by synchrotron radiation and Compton scattering process in a highly relativistic jet ejected from the central engine. It is widely accepted that prompt emissions occur in the internal shock when shells with different velocities collide with each other, and the magnetic field advected by the jet from the central engine can be ordered in large scale. On the other hand, afterglows are often assumed to occur in the external shock when the jet collides with interstellar medium, and the magnetic field produced by the shock through, e.g. Weibel instability, is possibly random. In this paper, we calculate the polarization properties of the synchrotron self-Compton process from a highly relativistic jet, in which the magnetic field is randomly distributed in the shock plane. We also consider the generalized situation where a uniform magnetic component perpendicular to the shock plane is superposed on the random magnetic component. We show that, the polarization is hardly to be larger than $10\%$ if the seed electrons are isotropic in the jet frame. This may account for the observed upper limit of polarization in the afterglow phase of GRBs. In addition, if the random and uniform magnetic components decay with time in different speeds, then the polarization angle may change $90^{\circ}$ duration the temporal evolution.

It has been observed that photons in the prompt emission of some gamma-ray bursts (GRBs) are highly polarized. The high polarization is used by some authors to give a strict constraint on the Lorentz invariance violation (LIV). If the Lorentz invariance is broken, the polarization vector of a photon may rotate during its propagation. The rotation angle of polarization vector depends on both the photon energy and the distance of source. It is believed that if high polarization is observed, then the relative rotation angle (denoted by $\alpha$) of polarization vector of the highest energy photon with respect to that of the lowest energy photon should be no more than $\pi/2$. Otherwise, the net polarization will be severely suppressed, thus couldn't be as high as what was actually observed. In this paper, we will give a detailed calculation on the evolution of GRB polarization arising from LIV effect duration the propagation. It is shown that the polarization degree rapidly decrease as $\alpha$ increases, and reaches a local minimum at $\alpha\approx \pi$, then increases until $\alpha\approx 3\pi/2$, after that decreases again until $\alpha \approx 2\pi$, etc. The polarization degree as a function of $\alpha$ oscillates with a quasi-period $T\approx \pi$, while the oscillating amplitude gradually decreases to zero. Moreover, we find that a considerable amount (more than $60\%$ of the initial polarization) of polarization degree can be conserved when $\alpha\approx \pi/2$. The polarization observation in a higher and wider energy band, a softer photon spectrum, and a higher redshift GRB is favorable in order to tightly constrain LIV effect.

LS 5039 is so far the best studied $\gamma$-ray binary system at multi-wavelength energies. A time resolved study of its spectral energy distribution (SED) shows that above 1 keV its power output is changing along its binary orbit as well as being a function of energy. To disentangle the energy dependence of the power output as a function of orbital phase, we investigated in detail the orbital light curves as derived with different telescopes at different energy bands. We analysed the data from all existing \textit{INTEGRAL}/IBIS/ISGRI observations of the source and generated the most up-to-date orbital light curves at hard X-ray energies. In the $\gamma$-ray band, we carried out orbital phase-resolved analysis of \textit{Fermi}-LAT data between 30 MeV and 10 GeV in 5 different energy bands. We found that, at $\lesssim$100 MeV and $\gtrsim$1 TeV the peak of the $\gamma$-ray emission is near orbital phase 0.7, while between $\sim$100 MeV and $\sim$1 GeV it moves close to orbital phase 1.0 in an orbital anti-clockwise manner. This result suggests that the transition region in the SED at soft $\gamma$-rays (below a hundred MeV) is related to the orbital phase interval of 0.5--1.0 but not to the one of 0.0--0.5, when the compact object is "behind" its companion. Another interesting result is that between 3 and 20 GeV no orbital modulation is found, although \textit{Fermi}-LAT significantly ($\sim$18$\sigma$) detects LS 5039. This is consistent with the fact that at these energies, the contributions to the overall emission from the inferior conjunction phase region (INFC, orbital phase 0.45 to 0.9) and from the superior conjunction phase region (SUPC, orbital phase 0.9 to 0.45) are equal in strength. At TeV energies the power output is again dominant in the INFC region and the flux peak occurs at phase $\sim$0.7.

The response of a detector to gravitational wave is a function of frequency. When the time a photon moving around in the Fabry-Perot cavities is the same order of the period of a gravitational wave, the phase-difference due to the gravitational wave should be an integral along the path. We present a formula description for detector response to gravitational wave with varied frequencies. The LIGO data for GW150914 and GW 151226 are reexamined in this framework. For GW150924, the traveling time of a photon in the LIGO detector is just a bit larger than a half period of the highest frequency of gravitational wave and the similar result is obtained with LIGO and Virgo collaborations. However, we are not always so luck. In the case of GW151226, the time of a photon traveling in the detector is larger than the period of the highest frequency of gravitational wave and the announced signal cannot match well the template with the initial black hole masses 14.2M$_\odot$ and 7.5M$_\odot$.

Einstein's equivalence principle (EEP) can be tested by the time delay between photons with different energies passing through a gravitational field. As one of the most energetic explosions in the Universe, gamma-ray bursts (GRBs) provide an effective tool to test the accuracy of EEP. In this paper, we use the continuous spectra of 20 short GRBs detected by the Swift/BAT to test the validity of EEP. Taking the duration of GRBs as the upper limit of the time delay induced by EEP violation (assuming that the high energy photons arrive later than the low energy photons), the difference of the parameterized post-Newtonian parameter is constrained with high accuracy. The strictest constraint, $|\gamma(150~{\rm keV})-\gamma(15~{\rm keV})|<5.59\times 10^{-10}$ from GRB 150101B, is about $1\sim 2$ orders of magnitude tighter than previous constraints. Moreover, our result is more statistically significant than previous results because we use the continuous spectra instead of isolated photons.

We use two different methods, i.e., dipole-fitting (DF) and hemisphere comparison (HC), to search for the anisotropic signals hiding in the Union2.1 data set. We find that the directions of maximum matter density derived using these two methods are about $114^{\circ}$ away from each other. We construct four Union2.1-like mock samples to test the statistical significance of these two methods. It is shown that DF method is statistically significant, while HC method is strongly biased by the distribution of data points in the sky. Then we assume that the anisotropy of distance modulus is mainly induced by the anisotropy of matter density, which is modeled to be the dipole form $\Omega_M=\Omega_{M0}(1-\cos\theta)$. We fit our model to Union2.1, and find that the direction of maximum matter density is well consistent with the direction derived using DF method, but it is very different from the direction previously claimed. Monte Carlo simulations show that our method is statistically more significant than HC method, although it is not as significant as DF method. The statistical significance can be further improved if the data points are homogeneously distributed in the sky. Due to the low quality of present supernovae data, however, it is still premature to claim that the Universe has any preferred direction.

The spectra of gamma-ray bursts (GRBs) in a wide energy range can usually be well described by the Band function, which is a two smoothly jointed power laws cutting at a breaking energy. Below the breaking energy, the Band function reduces to a cut-off power law, while above the breaking energy it is a simple power law. However, for some detectors (such as the Swift-BAT) whose working energy is well below or just near the breaking energy, the observed spectra can be fitted to cut-off power law with enough precision. Besides, since the energy band of Swift-BAT is very narrow, the spectra of most GRBs can be fitted well even using a simple power law. In this paper, with the most up-to-date sample of Swift-BAT GRBs, we study the effect of different spectral models on the empirical luminosity correlations, and further investigate the effect on the reconstruction of GRB Hubble diagram. We mainly focus on two luminosity correlations, i.e., the Amati relation and Yonetoku relation. We calculate these two luminosity correlations on both the case that the GRB spectra are modeled by Band function and cut-off power law. It is found that both luminosity correlations only moderately depend on the choice of GRB spectra. Monte Carlo simulations show that Amati relation is insensitive to the high-energy power-law index of the Band function. As a result, the GRB Hubble diagram calibrated using luminosity correlations is almost independent on the GRB spectra.

In the centenary year of Einstein's General Theory of Relativity, this paper reviews the current status of gravitational wave astronomy across a spectrum which stretches from attohertz to kilohertz frequencies. Sect. 1 of this paper reviews the historical development of gravitational wave astronomy from Einstein's first prediction to our current understanding the spectrum. It is shown that detection of signals in the audio frequency spectrum can be expected very soon, and that a north-south pair of next generation detectors would provide large scientific benefits. Sect. 2 reviews the theory of gravitational waves and the principles of detection using laser interferometry. The state of the art Advanced LIGO detectors are then described. These detectors have a high chance of detecting the first events in the near future. Sect. 3 reviews the KAGRA detector currently under development in Japan, which will be the first laser interferometer detector to use cryogenic test masses. Sect. 4 of this paper reviews gravitational wave detection in the nanohertz frequency band using the technique of pulsar timing. Sect. 5 reviews the status of gravitational wave detection in the attohertz frequency band, detectable in the polarisation of the cosmic microwave background, and discusses the prospects for detection of primordial waves from the big bang. The techniques described in sects. 1-5 have already placed significant limits on the strength of gravitational wave sources. Sects. 6 and 7 review ambitious plans for future space based gravitational wave detectors in the millihertz frequency band. Sect. 6 presents a roadmap for development of space based gravitational wave detectors by China while sect. 7 discusses a key enabling technology for space interferometry known as time delay interferometry.

Part of a review paper entitled "Gravitational wave astronomy: the current status.", appeared in " Science China Physics, Mechanics & Astronomy 58.12 (2015): 1-41.

We probe the possible anisotropy of the Universe by using the JLA compilation of type-Ia supernovae. We apply the Markov Chain Monte Carlo (MCMC) method to constrain the amplitude and direction of anisotropy in three cosmological models. For the dipole-modulated $\Lambda$CDM model, the anisotropic amplitude is consistent with zero at $68\%$ C.L., and has an upper bound $A_D<1.98\times10^{-3}$ at $95\%$ C.L. Regardless of much larger uncertainty, we find the dipole direction of JLA is amazingly opposite to that of Union2. Similar results are found for the dipole-modulated $w$CDM and CPL models. Thus, the Universe is still well consistent with the isotropy according to the JLA compilation.

In quantum gravity, a foamy structure of space-time leads to Lorentz invariance violation (LIV). As the most energetic astrophysical processes in the Universe, gamma-ray bursts (GRBs) provide an effective way to probe quantum gravity effects. We use continuous spectra of 20 short GRBs detected by the Swift satellite to give a conservative lower limit of quantum gravity energy scale $M_\textrm{QG} $. Due to the LIV effect, photons with different energy have different velocities. This will lead to the delayed arrival of high energy photons relative to the low energy ones. Based on the fact that the LIV-induced time delay can't be longer than the duration of a GRB, we present the most conservative estimation of the quantum gravity energy scales from 20 short GRBs. The most strict constraint, $M_\textrm{QG}>5.05\times10^{14}$ GeV, is from GRB 140622A.

Gamma-ray bursts (GRBs) are widely proposed as an effective probe to trace the Hubble diagram of the Universe in high redshift range. However, the calibration of GRBs is not as easy as that of type-Ia supernovae (SNe Ia). Most calibrating methods at present make use one or some of the empirical luminosity correlations, e.g., Amati relation. One of the underlying assumptions of these calibrating methods is that the empirical correlation is universal over all redshifts. In this paper, we check to what extent this assumption holds. Assuming that SNe Ia exactly trace the Hubble diagram of the Universe, we re-investigate the Amati relation for low redshift ($z<1.4$) and high redshift ($z>1.4$) GRBs, respectively. It is found that the Amati relation of low-$z$ GRBs differs from that of high-$z$ GRBs at more than $3\sigma$ confidence level. This result is insensitive to cosmological models. We should be cautious when using Amati relation to reconstruct the Hubble diagram of the Universe.

We propose the gravity's rainbow scenario as a possible alternative of the inflation paradigm to account for the flatness and horizon problems. We focus on studying the cosmological scalar perturbations which are seeded by the quantum fluctuations in the very early universe. The scalar power spectrum is expected to be nearly scale-invariant. We estimate the rainbow index $\lambda$ and energy scale $M$ in the gravity's rainbow scenario by analyzing the Planck temperature and WMAP polarization datasets. The constraints on them are given by $\lambda=2.933\pm0.012$ and $\ln (10^5M/M_p)= -0.401^{+0.457}_{-0.451}$ at the $68\%$ confidence level.

Type-I X-ray bursts on the surface of a neutron star are a unique probe to the accretion in X-ray binary systems. However, we know little about the feedback of the burst emission upon accretion. Hard X-ray shortages and enhancements of the persistent emission at soft X-rays have been observed. To put these findings in context with the aim of understanding the possible mechanism underneath, we investigated 68 bursts seen by RXTE from the clocked burster GS 1826--238. We diagnosed jointly the burst influence at both soft and hard X-rays, and found that the observations can be described as the CompTT model with variable normalization, electron temperature and optical depth. Putting these results in a scenario of coronal Compton cooling via the burst emission would lead to a shortage of the cooling power, which may suggest that additional consideration like the influence of the burst on the corona formation should be accounted for as well.

We suggest the universe is Finslerian in the stage of inflation. The Finslerian background spacetime breaks rotational symmetry and induces parity violation. The primordial power spectrum is given for quantum fluctuation of the inflation field. It depends not only on the magnitude of wavenumber but also on the preferred direction. We derive the gravitational field equations in the perturbed Finslerian background spacetime, and obtain a conserved quantity outside the Hubble horizon. The angular correlation coefficients are presented in our anisotropic inflation model. The parity violation feature of Finslerian background spacetime requires that the anisotropic effect only appears in angular correlation coefficients if $l'=l+1$. The numerical results of the angular correlation coefficients are given to describe the anisotropic effect.

Under the steady state condition, the spectrum of electrons is investigated by solving the continuity equation under the complex radiation of both the synchrotron and Compton processes. The resulted GRB spectrum is a broken power law in both the fast and slow cooling phases. On the basis of this electron spectrum, the spectral indices of the Band function in four different phases are presented. In the complex radiation frame, the detail investigation on physical parameters reveals that both the reverse shock photosphere model and the forward shock with strong coupling model can answer the $\alpha \sim -1$ problem. A possible marginal to fast cooling phase transition in GRB 080916C is discussed. The time resolved spectra in different pulses of GRB 100724B, GRB 100826A and GRB 130606B are investigated. We found that the flux is proportional to the peak energy in almost all pulses. The phases for different pulses are determined according to the spectral index revolution. We found the strong correlations between spectral indices and the peak energy in GRB 100826A, which can be explained by the Compton effect in the fast cooling phase. However, the complex scenario predicts a steeper index for the injected electrons, which challenges the acceleration mechanism in GRBs.

The influence of cosmological constant type dark energy in the early universe is investigated. This is accommodated by a new dispersion relation in de Sitter spacetime. We perform a global fitting to explore the cosmological parameters space by using the CosmoMC package with the recently released Planck TT and WMAP Polarization datasets. Using the results from global fitting, we compute a new CMB temperature-temperature spectrum. The obtained TT spectrum has lower power compared with the one based on $\Lambda$CDM model at large scales.

The polarization measurement is an important tool to probe the prompt emission mechanism in gamma-ray bursts (GRBs). The synchrotron photons can be scattered by cold electrons in the outflow via Compton scattering processes. The observed polarization depends on both the photon energy and the viewing angle. With the typical bulk Lorentz factor $\Gamma \sim 200$, photons with energy $E>10$ MeV tend to have smaller polarization than photons with energy $E<1$ MeV. At the right viewing angle, i.e. $\theta \sim \Gamma^{-1}$, the polarization achieves its maximal value, and the polarization angle changes $90^{\circ}$ relative to the initial polarization direction. Thus, the synchrotron radiation plus Compton scattering model can naturally explain the $90^{\circ}$ change of the polarization angle in GRB 100826A.

The polarimetric measurement on the prompt phase of GRB 100826A shows that the polarization angle changes $\sim 90^{\circ}$ between two adjacent time intervals. We will show that this phenomenon can be naturally interpreted in the framework of the magnetic-dominated jet (MDJ) model. The MDJ model suggests that the bulk Lorentz factor of the outflow increases as $\Gamma\propto r^{1/3}$, until reaching a saturated value $\Gamma_{\rm sat}$. Electrons move in the globally ordered magnetic field advected by the jet from the central engine and produce synchrotron photons. The polarized synchrotron photons travel alone the jet direction and then collide with the cold electrons at the front of the jet. After the Compton scattering process, these photons escape from the jet and are detected by the observer locating slightly off-axis. If photons are emitted before the bulk Lorentz factor saturates, the change of polarization angle is a natural result of the acceleration of the outflow.

Type-Ia supernovae (SNe Ia) are often used as the standard candles to probe the anisotropic expansion of the Universe. In this paper, we make a comprehensive comparison between the hemisphere comparison (HC) method and dipole-fitting (DF) method in searching for the cosmological preferred direction using the Union2 dataset, a compilation of 557 well-calibrated SNe Ia. We find that the directions of the faintest SNe Ia derived from these two methods are approximately opposite. Monte Carlo simulations show that the results of the HC method strongly depend on the distribution of the data points in the sky. The coincidence that the HC method and DF method give two completely opposite directions may be due to the extremely nonuniform distribution of the Union2 dataset.

The high polarization observed in the prompt phase of some gamma-ray bursts (GRBs) arouses extensive studies on the emission mechanism. In this paper, we investigate the polarization properties of the synchrotron-self-Compton (SSC) process from a highly relativistic jet. A magnetic-dominated, baryon-loaded jet ejected from the central engine travels with a large Lorentz factor. Shells with slightly different velocities collide with each other and produce shocks. The shocks accelerate electrons to power-law distribution, and at the same time, magnify the magnetic field. Electrons move in the magnetic field and produce synchrotron photons. The synchrotron photons suffer from the Compton scattering (CS) process and then are detected by an observer locating slightly off-axis. We derive analytically the formulae of photon polarization in the SSC process in two magnetic configurations: magnetic field in the shock plane and perpendicular to the shock plane. We show that photons induced by the SSC process can be highly polarized, with the maximum polarization $\Pi \sim 24\%$ in the energy band $[0.5,5]$ MeV. The polarization depends on the viewing angles, peaking in the plane perpendicular to the magnetic field. In the energy band $[0.05,0.5]$ MeV, in which most $\gamma$-ray polarimeters are active, the polarization is about twice of that in the Thomson limit, reaching to $\Pi\sim 20\%$. This implies that the Klein-Nishina effect, which is often neglected in literatures, should be carefully considered.

CosmoMC is a Fortran 95 Markov-Chain Monte-Carlo (MCMC) engine to explore the cosmological parameter space, plus a Python suite for plotting and presenting results (see http://cosmologist.info/cosmomc/). This document describes the installation of the CosmoMC on a Linux system (Ubuntu 14.04.1 LTS 64-bit version). It is written for those who want to use it in their scientific research but without much training on Linux and the program. Besides a step-by-step installation guide, we also give a brief introduction of how to run the program on both a desktop and a cluster. We share our way to generate the plots that are commonly used in the references of cosmology. For more information, one can refer to the CosmoCoffee forum (http://cosmocoffee.info/viewforum.php?f=11) or contact the authors of this document. Questions and comments would be much appreciated.

Recently, an anisotropic cosmological model was proposed. An arbitrary 1-form, which picks out a privileged axis in the universe, was added to the Friedmann-Robertson-Walker line element. The distance-redshift relation was modified such that it is direction dependent. In this paper, we use the Union2 dataset and 59 high-redshift gamma-ray bursts to give constraints on the anisotropy of the universe. The results show that the magnitude of anisotropy is about $D=-0.044\pm0.018$, and the privileged axis points towards the direction $(l_0,b_0)=(306.1^{\circ}\pm 18.7^{\circ},-18.2^{\circ}\pm 11.2^{\circ})$ in the galactic coordinate system. The anisotropy is small and the isotropic cosmological model is an excellent approximation.

Both the Wilkinson Microwave Anisotropy Probe (WMAP) and Planck observations reported the hemispherical asymmetry of the cosmic microwave background (CMB) temperature fluctuation. The hemispherical asymmetry might be stemmed from the primordial statistical anisotropy during the inflationary era of the universe. In this paper, we study possible implications of the primordial power spectra with dipolar anisotropy on the CMB temperature fluctuation and polarizations. We explicitly show that the statistical dipolar anisotropy may induce the off-diagonal (\(\ell'\neq\ell\)) $TT$, $EE$, $BB$, and $TE$ correlations, as well as the diagonal (\(\ell'=\ell\)) $TB$ and $EB$ spectra. In particular, these correlation coefficients are expected to be \(m\)-dependent generically. These signals of statistical anisotropy might be tested by CMB observations in future.

Recent released Planck data and other astronomical observations show that the universe may be anisotropic on large scales. Inspired by this, anisotropic cosmological models have been proposed. We note that the Finsler-Randers spacetime provides an appropriate framework for the anisotropic cosmology. By adding an arbitrary 1-form to the Friedmann-Robertson-Walker (FRW) line element, a privileged axis in the universe is picked out. The distance-redshift relation is modified to be direction-dependent. We wish that the anisotropic cosmological model may be tested crossly by independent observations. Type-Ia supernovae (SNe Ia) calibrated from four different light curve fitters are used to constrain possible anisotropy of the universe. The magnitudes of anisotropy are all between 2% --- 5%, but the systematic uncertainty cannot be excluded. The directions of privileged axis seem to differ from each other. The statistical significance is not high enough to make a convincing conclusion. Nevertheless, the $1\sigma$ contours in the $(l,b)$ plane obtained from four groups of SNe Ia have an overlap, centering at $(l,b)\approx (170^{\circ},0^{\circ})$. Monte Carlo simulation shows that the anisotropy is unlikely to be caused by selection effect.

Synchrotron radiation and Compton scattering are widely accepted as the most likely emission mechanisms of some astrophysical phenomena, such as gamma-ray bursts (GRBs) and active galactic nuclei (AGNs). The measurement on polarization of photons provides a useful tool to distinguish different emission mechanisms and structures of the emission region. Based on the differential cross section of a polarized photon scattered by an unpolarized electron of any initial momentum, we derive analytical formula of polarization for beamed photons scattered by isotropic electrons with a power law distribution. Numerical calculations are carried out in four special cases: electrons at rest, Thomson limit, head-on collision and monochromatic electrons. It is found that the maximum polarization can be as high as $100\%$ for low energy photons, if the electrons are at rest. Although polarization is highly suppressed due to the isotropic electrons, a maximum value of $\sim 10\% \-- 20\%$ can still be achieved. Compton scattering process can be used to explain the polarization of GRB 041219A and GRB 100826A.

The WMAP and Planck observations show that the quadrupole and octopole orientations of the CMB might align with each other. We reveal that the quadrupole--octopole alignment is a natural implication of the primordial power spectrum in an anisotropic spacetime. The primordial power spectrum is presented with a dipolar modulation. We obtain the privileged plane by employing the ``power tensor'' technique. At this plane, there is the maximum correlation between quadrupole and octopole. The probability for the alignment is much larger than what in the isotropic universe. We find that this model would lead to deviations from the statistical isotropy only for low--\(\ell\) multipoles.

Recently, the Planck 2013 results showed possible evidence for a dipolar power modulation of the CMB temperature fluctuations at low--(\ell) multipoles. This anomaly might imply certain deviations from statistical isotropy. To incorporate the Planck's data into standard cosmological model, we propose an inflation model of the very early universe in an anisotropic spacetime. A generalized Friedmann-Robertson-Walker (FRW) metric is presented in the Randers-Finsler spacetime. We obtain the primordial power spectrum of the scalar perturbation with direction dependence, such as the dipolar modulation. This is consistent with the dipolar power modulation of the CMB anisotropy signaled by the Planck observation.

We present a scalar perturbation for the $\Lambda$CDM model, which breaks the isotropic symmetry of the universe. Based on the Union2 data, the least-$\chi^2$ fit of the scalar perturbed $\Lambda$CDM model shows that the universe has a preferred direction $(l,b)=(287^\circ\pm25^\circ,11^\circ\pm22^\circ)$. The magnitude of scalar perturbation is about $-2.3\times10^{-5}$. The scalar perturbation for the $\Lambda$CDM model implies a peculiar velocity, which is perpendicular to the radial direction. We show that the maximum peculiar velocities at redshift $z=0.15$ and $z=0.015$ equal to $73\pm28 \rm km\cdot s^{-1}$ and $1099\pm427 \rm km\cdot s^{-1}$, respectively. They are compatible with the constraints on peculiar velocity given by Planck Collaboration.

Based on the quantum electrodynamics, we present a generic formalism of the polarization for beamed monochromatic photons scattered by electrons in any spectral distribution. The formulae reduce to the components of the Fano matrix when electrons are at rest. We mainly investigate the polarization in three scenarios, i.e., electrons at rest, isotropic electrons with a power law spectrum and thermal electrons. If the incident beam is polarized, the polarization is reduced significantly by isotropic electrons at large viewing angles, and the degree of polarization due to thermal electrons is about one times less than that of electrons in a power law. If the incident bean is unpolarized, soft $\gamma$-rays can lead to about 15% polarization at viewing angles around $\pi/4$. For isotropic electrons, one remarkable feature is that the polarization as a function of the incident photon energy always peaks roughly at 1 MeV, this is valid for both the thermal and power law cases. This feature can be used to distinguish the model of the inverse Compton scattering from that of the synchrotron radiation.

The gravitational lensing of Bullet Clusters and early-type galaxies pose serious challenges on the validity of MOND. Recently, Finslerian MOND, a generalization of MOND in the framework of Finsler gravity, has been proposed to explain the mass discrepancy problem of Bullet Cluster 1E 0657\ 558. In this paper, we check the validity of the Finslerian MOND in describing the strong gravitational lensing of early-type galaxies. The investigation on ten strong lenses of the CASTLES samples shows that there is no strong evidence for the existence of dark matter.

Based on Cohen \& Glashow's very special relativity [A. G. Cohen and S. L. Glashow, Phys. Rev. Lett. {\bf 97} (2006) 021601], we propose an anisotropic modification to the Friedmann-Robertson-Walker (FRW) line element. An arbitrarily oriented 1-form is introduced and the FRW spacetime becomes of the Randers-Finsler type. The 1-form picks out a privileged axis in the universe. Thus, the cosmological redshift as well as the Hubble diagram of the type Ia supernovae (SNe Ia) becomes anisotropic. By directly analyzing the Union2 compilation, we obtain the privileged axis pointing to \((l,b)=({304^\circ}\pm{43^\circ},{-27^\circ}\pm{13^\circ})\) (\(68\%~\rm{C.L.}\)). This privileged axis is close to those obtained by comparing the best-fit Hubble diagrams in pairs of hemispheres. It should be noticed that the result is consistent with isotropy at the \(1\sigma\) level since the anisotropic magnitude is \(D=0.03\pm 0.03\).

The observations on galaxy rotation curves show significant discrepancies from the Newtonian theory. This issue could be explained by the effect of the anisotropy of the spacetime. Conversely, the spacetime anisotropy could also be constrained by the galaxy rotation curves. Finsler geometry is a kind of intrinsically anisotropic geometry. In this paper, we study the effect of the spacetime anisotropy at the galactic scales in the Finsler spacetime. It is found that the Finslerian model has close relations with the Milgrom's MOND. By performing the best-fit procedure to the galaxy rotation curves, we find that the anisotropic effects of the spacetime become significant when the Newtonian acceleration \(GM/r^2\) is smaller than the critical acceleration \(a_0\). Interestingly, the critical acceleration \(a_0\), although varies between different galaxies, is in the order of magnitude \(cH_0/2\pi\sim 10^{-10} \rm{m\,\, s^{-2}}\).

Astronomical observations showed that there may exist a bulk flow with peculiar velocities in the universe, which contradicts with the (\Lambda)CDM model. The bulk flow reveals that the observational universe is anisotropic at large scales. In fact, a more reliable observation on the anisotropy of spacetime comes from the CMB power spectra. The WMAP and Planck satellites both show that there is a hemispherical power asymmetry at large-angular scales. In this paper, we propose a "wind" scenario to the bulk flow (or the anisotropy of spacetime). Under the influence of the "wind", the spacetime metric could become a Finsler structure. By resolving the null geodesic equation, we obtain the modified luminosity distance, which has a dipolar form at the leading order. Thus, the "wind" describes well the bulk flow. In addition, we perform a least-(\chi^2) fit to the data of type Ia supernovae (SNe Ia) in the Union2.1 compilation. The peculiar velocity of the bulk flow has an upper limit (v_{bulk}\lesssim 4000 \rm{km/s}), which is compatible with all the existing observational values.

The data of the Bullet Cluster 1E0657-558 released on November 15, 2006 reveal that the strong and weak gravitational lensing convergence $\kappa$-map has an $8\sigma$ offset from the $\Sigma$-map. The observed $\Sigma$-map is a direct measurement of the surface mass density of the Intracluster medium(ICM) gas. It accounts for 83% of the averaged mass-fraction of the system. This suggests a modified gravity theory at large distances different from Newton's inverse-square gravitational law. In this paper, as a cluster scale generalization of Grumiller's modified gravity model (D. Grumiller, Phys. Rev. Lett. 105, 211303 (2010)), we present a gravity model with a generalized linear Rindler potential in Randers-Finslerian spacetime without invoking any dark matter. The galactic limit of the model is qualitatively consistent with the MOND and Grumiller's. It yields approximately the flatness of the rotational velocity profile at the radial distance of several kpcs and gives the velocity scales for spiral galaxies at which the curves become flattened. Plots of convergence $\kappa$ for a galaxy cluster show that the peak of the gravitational potential has chances to lie on the outskirts of the baryonic mass center. Assuming an isotropic and isothermal ICM gas profile with temperature $T=14.8$ keV (which is the center value given by observations), we obtain a good match between the dynamical mass $M_\textmd{T}$ of the main cluster given by collisionless Boltzmann equation and that given by the King $\beta$-model. We also consider a Randers+dark matter scenario \textbf{and a $\Lambda$-CDM model} with the NFW dark matter distribution profile. We find that a mass ratio $\eta$ between dark matter and baryonic matter about 6 fails to reproduce the observed convergence $\kappa$-map for the isothermal temperature $T$ taking the observational center value.

Gravitational field equations in Randers-Finsler space of approximate Berwald type are investigated. A modified Friedmann equation and a new luminosity distance-redshift relation is proposed. A best-fit to the Type Ia supernovae (SNe) observations yields that the $\Omega_{\Lambda}$ in the $\Lambda$-CDM model is suppressed to almost zero. This fact indicates that the astronomical observations on the Type Ia SNe can be described well without invoking any form of dark energy. The best-fit age of the universe is given. It is in agreement with the age of our galaxy.

It is known that theory of MOND with spherical symmetry cannot account for the convergence $\kappa$-map of Bullet Cluster 1E0657-558. In this paper, we try to set up a Finslerian MOND, a generalization of MOND in Finsler spacetime. We use $Ric=0$ to obtain the gravitational vacuum field equation in a four-dimensional Finsler spacetime. To leading order in the post-Newtonian approximation, we obtain the explicit form of the Finslerian line element. It is simply the Schwarzschild's metric except for the Finslerian rescaling coefficient $f(v)$ of the radial coordinate $r$, i.e. $R=f(v(r))r$. By setting $f(v(r))=(1-\sqrt{a_0r^2/GM})^{-1}$, we obtain the famous MOND in a Finslerian framework. Taking a dipole and a quadrupole term into consideration, we give the convergence $\kappa$ in gravitational lensing astrophysics in our model. Numerical analysis shows that our prediction is to a certain extent in agreement with the observations of Bullet Cluster 1E0657-558. With the theoretical temperature $T$ taking the observed value 14.8 keV, the mass density profile of the main cluster obtained in our model is the same order as that given by the best-fit King $\beta$-model.

The effective potential in the Grumiller's modified gravity [D. Grumiller, Phys. Rev. Lett {\bf 105}, 211303 (2010)] includes the Newtonian potential and a Rindler term. The fitting to the rotation curve data of eight galaxies suggests a universal Rindler acceleration $a\approx 0.30\times 10^{-10}$ m s$^{-2}$. We do a two-parameter fit first, with the mass-to-light ratio ($\Upsilon_*$) and the Rindler acceleration ($a$) as free parameters. It is found that the data of six out of the eight galaxies fit well with the prediction of theory in the range $0\lesssim r \lesssim 40$ kpc, although the theoretical curves show a tendency of arising beyond this range. The Rindler accelerations of the six well-fitted galaxies have the same magnitude, with an average value $\bar{a}\approx0.30\times 10^{-10}$ m s$^{-2}$. Inspired by this fact, we then carry out a one-parameter ($\Upsilon_*$) fit to the six galaxies, with $a$ fixed at $\bar{a}$, and find that the theory can still reproduce the observation. The value of Rindler acceleration we get here is a quarter of that of Milgrom's MOND. For the rest two galaxies, NGC5055 and DDO154, the fitting results are significantly improved if the photometric scale length ($h$) is included as another free parameter.

Recent observations by the Fermi-LAT showed that there are delayed arrivals of GeV photons relative to the onset of MeV photons in some GRBs. In order to avoid a large optical depth, the minimal value of the Lorentz factor has been estimated to be higher than 1000 in some brightest bursts. In this paper, we present a detailed calculation of the time delay between the MeV and GeV photons in the framework of the magnetic-dominated jet model. We find that the time delay strongly depends on the saturated bulk Lorentz factor of the jet. Inspired by this fact, we use this model to calculate the Lorentz factors of the four brightest Fermi bursts. The results indicate that the Lorentz factors are much smaller than that obtained from the "single-zone" scenario. The short GRB 090510 has a minimal Lorentz factor 385, while the three long bursts GRB 080916c, GRB090902b and GRB 090926 have almost the same Lorentz factors, with an average value near 260. Another interesting result is that, for long bursts, GeV photons are emitted after the bulk Lorentz factor saturates. For the short GRB, however, MeV and GeV photons are emitted at the same phase, i.e., either at the expansion phase or at the coast phase.

Possible Lorentz invariance violation (LIV) has been investigated for a long time based on observations of GRBs . These arguments relied on the assumption that photons with different energy are emitted at the same place and time. In this work, we try to take account of the intrinsic time delay $\Delta t_{\rm int}$ between emissions of low and high energy photons by using the magnetic jet model. The possible LIV effects are discussed in a unified scenario both for long and short {\it Fermi}-detected GRBs. This leads to a unique quantum gravity energy scale $M_1c^2 \sim 1.0 \times 10^{20}$ GeV respecting the linear dispersion relation.

We investigate the possible effect of cosmological-constant type dark energy during the inflation period of the early universe. This is accommodated by a new dispersion relation in de Sitter space. The modified inflation model of a minimally-coupled scalar field is still able to yield an observation-compatible scale-invariant primordial spectrum, simultaneously having potential to generate a spectrum with lower power at large scales. A qualitative match to the WMAP 7-year data is presented. We obtain an $\Omega_\Lambda$ of the same order of that in the $\Lambda$-CDM model. Possible relations between the de Sitter scenario and the Doubly Special Relativity(DSR) are also discussed.

In Verlinde's entropic force scenario of gravity, Newton's laws and Einstein equations can be obtained from the first pinciples and general assumptions. However, the equipartition law of energy is invalid at very low temperatures. We show clearly that the threshold of the equipartition law of energy is related with horizon of the universe. Thus, a one-dimension Debye (ODD) model in the direction of radius of the modified entropic force (MEF) maybe suitable in description of the accelerated expanding universe. We present a Friedmann cosmic dynamical model in the ODD-MEF framework. We examine carefully constraints on the ODD-MEF model from the Union2 compilation of the Supernova Cosmology Project (SCP) collaboration, the data from the observation of the large-scale structure (LSS) and the cosmic microwave background (CMB), i.e. SNe Ia+LSS+CMB. The combined numerical analysis gives the best-fit value of the model parameters $\zeta\simeq10^{-9}$ and $\Omega_{m0}=0.224$, with $\chi_{min}^2=591.156$. The corresponding age of the universe agrees with the result of D. Spergel {\it et al.}\cite{Spergel2003} at 95% confidence level. The numerical result also yields an accelerated expanding universe without invoking any kind of dark energy. Taking $\zeta(\equiv 2\pi \omega_D/H_0)$ as a running parameter associated with the structure scale $r$, we obtain a possible unified scenario of the asymptotic flatness of the radial velocity dispersion of spiral galaxies, the accelerated expanding universe and the Pioneer 10/11 anomaly in the entropic force framework of Verlinde.

We use the Chandra measurements of the X-ray gas mass fraction of 26 rich clusters released by Allen et al. to perform constraints on the holographic dark energy model. The constraints are consistent with those from other cosmological tests, especially with the results of a joint analysis of supernovae, cosmic microwave background, and large scale structure data. From this test, the holographic dark energy also tends to behave as a quintom-type dark energy.

We discuss the threshold anomaly in ultra-high energy cosmic ray physics by assuming that the matter world just be perturbation of a de Sitter universe, which is consistent with the recent astronomical observations: about two thirds of the whole energy in the universe is contributed by a small positive cosmological constant. One-particle states are presented explicitly. It is noticed that the dispersion relation of free particles is dependent on the degrees of freedom of angular momentum. This fact can be regarded as the effects of the cosmological constant on kinematics of particles.

We show that a tiny but non-zero positive cosmological constant, which has been strongly suggested by the recent astronomical observations on supernovae and CMBR, may change notably the behaviors of the ultrahigh energy cosmic ray interacting with soft photons. The threshold anomaly of the ultrahigh energy cosmic ray disappears naturally after the effects of the cosmological constant are taken into account.

We present a covariant framework of kinematics in the dS spacetime, which is a natural postulate of recent astronomical observations ($\Lambda>0)$. One-particle states are presented explicitly. It is noticed that the dispersion relation of free particles is dependent on the degrees of freedom of angular momentum and spin. This fact can be referred to as the effects of the cosmological constant on kinematics of particles. The kinematics in dS spacetime is used to investigate the phenomenon of ultra high energy cosmic rays. We emphasize the possibility of solving the threshold anomalies of the interactions between ultra high energy cosmic rays and soft photons in the covariant framework of kinematics in the dS spacetime.

We calculate the vacuum energy density by taking account of different massive scalar fields in AdS spacetime. It is found that the mass spectrum of a scalar field in AdS spacetime is discrete because of a natural boundary condition. The results match well with the observed cosmological constant.