The discovery of gravitational waves, which are ripples of space-time itself, opened a new window to test general relativity, because it predicts that there are only plus and cross polarizations for gravitational waves. For alternative theories of gravity, there may be up to six polarizations. The measurement of the polarization is one of the major scientific goals for future gravitational wave detectors. To evaluate the capability of the detector, we need to use the frequency dependent response functions averaged over the source direction and polarization angle. We derive the full analytical formulas of the averaged response functions for all six possible polarizations and present their asymptotic behaviors based on these analytical formulas. Compared with the numerical simulation, the full analytical formulas are more efficient and valid for any equal-arm interferometric gravitational wave detector without optical cavities in the arms and for a time-delay-interferometry Michelson combination.

The possibility that in the mass range around $10^{-12}\ M_\odot$ most of dark matter constitutes of primordial black holes (PBHs) is a very interesting topic. To produce PBHs with this mass, the primordial scalar power spectrum needs to be enhanced to the order of 0.01 at the scale $k\sim 10^{12}\ \text{Mpc}^{-1}$. The enhanced power spectrum also produces large secondary gravitational waves at the mHz band. A phenomenological delta function power spectrum is usually used to discuss the production of PBHs and secondary gravitational waves. Based on G and k inflations, we propose a new mechanism to enhance the power spectrum at small scales by introducing a non-canonical kinetic term $[1-2G(\phi)]X$ with the function $G(\phi)$ having a peak. Away from the peak, $G(\phi)$ is negligible and we recover the usual slow-roll inflation which is constrained by the cosmic microwave background anisotrpy observations. Around the peak, the slow-roll inflation transiently turns to ultra slow-roll inflation. The enhancement of the power spectrum can be obtained with generic potentials, and there is no need to fine tune the parameters in $G(\phi)$. The energy spectrum $\Omega_{GW}(f)$ of secondary gravitational waves have the characteristic power law behaviour $\Omega_{GW}(f)\sim f^{n}$ and is testable by pulsar timing array and space based gravitational wave detectors.

Interesting chemically peculiar field stars may reflect their stellar evolution history and their possible origin in a different environment from where they are found now, which is one of the most important research fields in Galactic archaeology. To explore this further, we have used the CN-CH bands around 4000 A to identify N-rich metal-poor field stars in LAMOST DR3. Here we expand our N-rich metal-poor field star sample to ~100 stars in LAMOST DR5, where 53 of them are newly found in this work. We investigate light elements of the common stars between our sample and APOGEE DR14. While Mg, Al, and Si abundances generally agree with the hypothesis that N-rich metal-poor field stars come from enriched populations in globular clusters, it is still inconclusive for C, N, and O. After integrating the orbits of our N-rich field stars and a control sample of normal metal-poor field stars, we find that N-rich field stars have different orbital parameter distributions compared to the control sample, specifically, apocentric distances, maximum vertical amplitude (Zmax), orbital energy, and z direction angular momentum (Lz). The orbital parameters of N-rich field stars indicate that most of them are inner-halo stars. The kinematics of N-rich field stars support their possible GC origin. The spatial and velocity distributions of our bona fide N-rich field star sample are important observational evidence to constrain simulations of the origin of these interesting objects.

Standard stellar evolution model predicts a severe depletion of lithium (Li) abundance during the first dredge up process (FDU). Yet a small fraction of giant stars are still found to preserve a considerable amount of Li in their atmospheres after the FDU. Those giants are usually identified as Li-rich by a widely used criterion, A(Li) > 1.5 dex. A large number of works dedicated to searching for investigating this minority of the giant family, and the amount of Li-rich giants, has been largely expanded on, especially in the era of big data. In this paper, we present a catalog of Li-rich giants found from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) survey with Li abundances derived from a template-matching method developed for LAMOST low-resolution spectra.

All stellar mass black holes have hitherto been identified by X-rays emitted by gas that is accreting onto the black hole from a companion star. These systems are all binaries with black holes below 30 M$_{\odot}$$^{1-4}$. Theory predicts, however, that X-ray emitting systems form a minority of the total population of star-black hole binaries$^{5,6}$. When the black hole is not accreting gas, it can be found through radial velocity measurements of the motion of the companion star. Here we report radial velocity measurements of a Galactic star, LB-1, which is a B-type star, taken over two years. We find that the motion of the B-star and an accompanying H$\alpha$ emission line require the presence of a dark companion with a mass of $68^{+11}_{-13}$ M$_{\odot}$, which can only be a black hole. The long orbital period of 78.9 days shows that this is a wide binary system. The gravitational wave experiments have detected similarly massive black holes$^{7,8}$, but forming such massive ones in a high-metallicity environment would be extremely challenging to current stellar evolution theories$^{9-11}$.

Space-based gravitational wave detectors cannot keep rigid structures and precise arm length equality, so the precise equality of detector arms which is required in a ground-based interferometer to cancel the overwhelming laser noise is impossible. The time-delay interferometry method is applied to unequal arm lengths to cancel the laser frequency noise. We give analytical formulas of the averaged response functions for tensor, vector, breathing and longitudinal polarizations in different TDI combinations, and obtain their asymptotic behaviors. At low frequencies, $f\ll f_*$, the averaged response functions of all TDI combinations increase as $f^2$ for all six polarizations. The one exception is the averaged response functions of $\zeta$ for all six polarizations increase as $f^4$ in the equilateral-triangle case. At high frequencies, $f\gg f_*$, the averaged response functions of all TDI combinations for the tensor and breathing modes fall off as $1/f^2$, the averaged response functions of all TDI combinations for the vector mode fall off as $\ln(f)/f^2$ , and the averaged response functions of all TDI combinations for the longitudinal mode fall as $1/f$. We also give LISA and TianQin sensitivity curves in different TDI combinations for tensor, vector, breathing and longitudinal polarizations.

Stellar rotation plays a key role in stellar activity. The rotation period could be detected through light curve variations caused by starspots. Kepler provides two types of light curves, one is the Pre-search Data Conditioning (PDC) light curves, the other is the Simple Aperture Photometer (SAP) light curves. Compared with the PDC light curves, the SAP light curves keep the long-term trend, relatively suitable for searches of long period signals. However, SAP data are inflicted by some artefacts such as quarterly rolls and instrumental errors, making it difficult to find the physical periods in the SAP light curves. We explore a systematic approach based on the light curve pre-processing, period detection and candidate selection. We also develop a simulated light curve test to estimate our detection limits for the SAP-like LCs. After applying our method to the raw SAP light curves, we found more than 1000 main-sequence stars with the period longer than 30 days, 165 are newly discovered. Considering the potential flaw of the SAP, we also inspect the newly found objects with photometry methods, and most of our periodical signals are confirmed.

We report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 background events, which corresponds to 5.6$\sigma$ statistical significance. This is the first detection of photons with $E>100$ TeV from an astrophysical source.

The discovery of Li-rich giant has cast a new challenge for the standard stellar evolution models, and to resolve this issue, the number of this type object has been rapidly increased because of the development of worldwide surveys these days. Taking advantage of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope survey, 44 newly Li-rich giants are reported, which are confirmed with high-resolution observations. Based on the high-resolution and high signal-to-noise spectra, we derived the atmospheric parameters and elemental abundances with the spectral synthesis method. We performed a detailed analysis of their evolutionary stages, infrared excess, projected rotational velocity (v sin i), and the stellar population. We find that (1) The Li-rich giants concentrate at the evolutionary status of the red giant branch bump, red clump, and asymptotic giant branch; (2) Three of them are fast rotators and none exhibit infrared excess. Our results imply that the origins of Li enrichment are most likely to be associated with the extra mixing in the stellar interior, and the external sources maybe only make a minor contribution. Moreover, various Li-rich episodes take place at different evolutionary stages.

We give the procedure to reconstruct the extended inflationary potentials for attractors and use the $\alpha$ attractor and the constant-roll model as examples to show how to reconstruct the class of extended inflationary potentials in the strong coupling limit. We also derive the strong coupling condition which the coupling constant $\xi$ satisfies. The class of extended inflationary potentials has the same attractor, and the reconstructed extended inflationary potentials are consistent with the observational constraints. The existence of the attractors from extended inflationary models further challenges the model discrimination just by the observables $n_s$ and $r$.

About one percent of giants\upcite{Brown1989} are detected to have anomalously high lithium (Li) abundances in their atmospheres, conflicting directly with the prediction of the standard stellar evolution models\upcite{Iben1967}, and making the production and evolution of Li more intriguing, not only in the sense of the Big Bang nucleosynthesis\upcite{Cyburt2016,Spite1982} or the Galactic medium\upcite{Tajitsu2015}, but also the evolution of stars. Decades of efforts have been put into explaining why such outliers exist\upcite{Sackmann1999, Denissenkov2004, Charbonnel2010}, yet the origins of Li-rich giants are still being debated. Here we report the discovery of the most Li-rich giant known to date, with a super-high Li abundance of 4.51. This rare phenomenon was snapshotted together with another short-term event that the star is experiencing its luminosity bump on the red giant branch. Such high Li abundance indicates that the star might be at the very beginning of its Li-rich phase, which provides a great opportunity to investigate the origin and evolution of Li in the Galaxy. A detailed nuclear simulation is presented with up-to-date reaction rates to recreate the Li enriching process in this star. Our results provide tight constraints on both observational and theoretical points of view, suggesting that low-mass giants can produce Li inside themselves to a super high level via $^{7}$Be transportation during the red giant phase.

For the constant-roll tachyon inflation, we derive the analytical expressions for the scalar and tensor power spectra, the scalar and tensor spectral tilts and the tensor to scalar ratio up to the first order by using the method of Bessel function approximation. The derived $n_s-r$ results for the constant-roll inflation are also compared with the observations, we find that only one constant-roll inflation is consistent with the observations and observations constrain the constant-roll inflation to be slow-roll inflation. The tachyon potential is also reconstructed for the constant-roll inflation which is consistent with the observations.

The lithium abundances in a few percent of giants exceed the value predicted by the standard stellar evolution models, and the mechanisms of Li enhancement are still under debate. The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) survey has obtained over six million spectra in the past five years, and thus provides a great opportunity to search these rare objects and to more clearly understand the mechanisms of Li enhancement. Based on the high-resolution spectrum we obtained the stellar parameters ($T_\mathrm{eff}$, $\log g$, [Fe/H]), and determined the elemental abundances of Li, C, N, $\alpha$, Fe-peak, r-process, s-process elements, and the projected rotational velocity. For a better understanding of the effect of mixing processes, we also derived the $^{12}\rm{C}$ to $^{13}\rm{C}$ ratio, and constrained the evolutionary status of TYC\,3251-581-1 based on the BaSTI stellar isochrones. The super Li-rich giant TYC\,3251-581-1 has $\rm{A(Li)} = 3.51$, the average abundance of two lithium lines at $\lambda = 6708$ \AA\ and 6104 \AA\ based on the non-local thermodynamic equilibrium (NLTE) analysis. The atmospheric parameters show that our target locates on the luminosity function bump. The low carbon isotopic ratio ($^{12}\rm{C}/^{13}\rm{C} = 9.0 $), a slow rotational velocity $v\sin i = 2.2\ \rm{km\,s}^{-1}$, and no sign of IR excess suggest that additional mixing after first dredge up (FDU) should occur to bring internal synthesized Li to the surface. The low carbon ($[\rm{C}/\rm{Fe}] \sim -0.34$ ) and enhanced nitrogen ($[\rm{N}/\rm{Fe}] \sim 0.33$) are also consistent with the sign of mixing. Given the evolutionary stage of TYC\,3251-581-1 with the relatively low $^{12}\rm{C}/^{13}\rm{C}$, the internal production which replenishes Li in the outer layer is the most likely origin of Li enhancement for this star.

Flare research is becoming a burgeoning realm in the study of stellar activity due to the launch of {\it Kepler} in 2009. {\it Kepler} provides data with two time resolutions, i.e., the long-cadence (LC) data with a time resolution of 30 minutes and the short-cadence(SC) data with a time resolution of 1 minute, both of which can be used to study stellar flares. In this paper, we search flares in light curves with both LC data and SC data, and compare them in aspects of the true-flare rate, the flare energy, the flare amplitude, and the flare duration. It is found that LC data systematically underestimated the energies of flares by 25\%, and underestimated the amplitudes of flares by 60\% compared with SC flares. The duration are systematically overestimated by 50\% compared with SC flares. However, the above percentages are poorly constrained and there is a lot of scatter. About 60\% SC flares have not been detected by LC data. We investigate the limitation of LC data, and suggest that although LC data cannot reflect the detailed profiles of flares, they also can capture the basic properties of stellar flares.

We discuss the constant-roll inflation with constant $\epsilon_2$ and constant $\bar\eta$. By using the method of Bessel function approximation, the analytical expressions for the scalar and tensor power spectra, the scalar and tensor spectral tilts, and the tensor to scalar ratio are derived up to the first order of $\epsilon_1$. The model with constant $\epsilon_2$ is ruled out by the observations at the $3\sigma$ confidence level, and the model with constant $\bar\eta$ is consistent with the observations at the $1\sigma$ confidence level. The potential for the model with constant $\bar\eta$ is also obtained from the Hamilton-Jacobi equation. Although the observations constrain the constant-roll inflation to be slow-roll inflation, the $n_s-r$ results from the constant-roll inflation are not the same as those from the slow-roll inflation even when $\bar\eta\sim 0.01$.

We have studied the variability of \object{S5 0716+714} at radio 15 GHz and $\gamma$-ray band using three different methods. A possible periodicity of $P_{15 GHz}=266.0\pm11.5$ and $P_{\gamma}=344.0\pm16.4$ days are obtained for radio 15 GHz and $\gamma$-ray light curves, respectively. The variability may be related to the intrinsically emission mechanism. The difference between the variability timescales of radio 15 GHz and $\gamma$-ray may be due to that the emission of radio 15 GHz is produced via the synchrotron process, while the $\gamma$-ray is produced by both the SSC and EC processes.

The hybrid Tibet AS array was successfully constructed in 2014. It has 4500 m$^{2}$ underground water Cherenkov pools used as the muon detector (MD) and 789 scintillator detectors covering 36900 m$^{2}$ as the surface array. At 100 TeV, cosmic-ray background events can be rejected by approximately 99.99\%, according to the full Monte Carlo (MC) simulation for $\gamma$-ray observations. In order to use the muon detector efficiently, we propose to extend the surface array area to 72900 m$^{2}$ by adding 120 scintillator detectors around the current array to increase the effective detection area. A new prototype scintillator detector is developed via optimizing the detector geometry and its optical surface, by selecting the reflective material and adopting dynode readout. This detector can meet our physics requirements with a positional non-uniformity of the output charge within 10\% (with reference to the center of the scintillator), time resolution FWHM of $\sim$2.2 ns, and dynamic range from 1 to 500 minimum ionization particles.

By using the relations between the slow-roll parameters and the power spectrum for the single field slow-roll inflation, we derive the scalar spectral tilt $n_s$ and the tensor to scalar ratio $r$ for the constant slow-roll inflation and obtain the constraint on the slow-roll parameter $\eta$ from the Planck 2015 results. The inflationary potential for the constant slow-roll inflation is then reconstructed in the framework of both general relativity and scalar-tensor theory of gravity, and compared with the recently reconstructed E model potential. In the strong coupling limit, we show that the $\eta$ attractor is reached.

We derive the general formulae for the the scalar and tensor spectral tilts to the second order for the inflationary models with non-minimally derivative coupling without taking the high friction limit. The non-minimally kinetic coupling to Einstein tensor brings the energy scale in the inflationary models down to be sub-Planckian. In the high friction limit, the Lyth bound is modified with an extra suppression factor, so that the field excursion of the inflaton is sub-Planckian. The inflationary models with non-minimally derivative coupling are more consistent with observations in the high friction limit. In particular, with the help of the non-minimally derivative coupling, the quartic power law potential is consistent with the observational constraint at 95\% CL.

We present spectra of the extreme polar AR Ursae Majoris (AR UMa) which display a clear Al I absorption doublet, alongside spectra taken less than a year earlier in which that feature is not present. Re-examination of earlier SDSS spectra indicates that the Al I absorption doublet was also present $\approx$8 years before our first non-detection. We conclude that this absorbing material is unlikely to be on the surface of either the white dwarf (WD) or the donor star. We suggest that this Al I absorption feature arises in circumstellar material, perhaps produced by the evaporation of asteroids as they approach the hot WD. The presence of any remaining reservoir of rocky material in AR UMa might help to constrain the prior evolution of this unusual binary system. We also apply spectral decomposition to find the stellar parameters of the M dwarf companion, and attempt to dynamically measure the mass of the WD in AR UMa by considering both the radial velocity curves of the H$_\beta$ emission line and the Na I absorption line. Thereby we infer a mass range for the WD in AR UMa of 0.91 $M_{\odot}$ $<$ $M_{\mathrm{WD}}$ $<$ 1.24 $M_{\odot}$.

The observational data on the anisotropy of the cosmic microwave background constraints the scalar spectral tilt $n_s$ and the tensor to scalar ratio $r$ which depend on the first and second derivatives of the inflaton potential. The information can be used to reconstruct the inflaton potential in the polynomial form up to some orders. However, for some classes of potentials, $n_s$ and $r$ behave as $n_s(N)$ and $r(N)$ universally in terms of the number of e-folds $N$. The universal behaviour of $n_s(N)$ can be used to reconstruct a class of inflaton potentials. By parametrizing one of the parameters $n_s(N)$, $\epsilon(N)$ and $\phi(N)$, and fitting the parameters in the models to the observational data, we obtain the constraints on the parameters and reconstruct the classes of the inflationary models which include the chaotic inflation, T-model, hilltop inflation, s-dual inflation, natural inflation and $R^2$ inflation.

We present the detection of Cepheids in the barred spiral galaxy NGC 1313, using the Wide Field and Planetary Camera 2 on the Hubble Space Telescpoe. Twenty B (F450W) and V (F555W) epochs of observations spanning over three weeks were obtained, on which the profile-fitting photometry of all stars in the monitored field was performed using the package HSTphot. A sample of 26 variable stars have been identified to be Cepheids, with periods between 3 and 14 days. Based on the derived period-luminosity relations in B and V bands, we obtain an extinction-corrected distance modulus of mu = 28.32 +- 0.08 (random) +- 0.06 (systematic), employing the Large Magellanic Cloud as the distance zero point calibrator. The above moduli correspond to a distance of 4.61 +- 0.17 (random) +- 0.13 (systematic) Mpc, consistent with previous measurements reported in the literature with uncertainties. In addition, the reddening to NGC 1313 is found to be small.

Based on Kepler data, we present the results of a search for white-light flares on 1049 close binaries. We identify 234 flare binaries, on which 6818 flares are detected. We compare the flare-binary fraction in different binary morphologies ("detachedness"). The result shows that the fractions in over-contact and ellipsoidal binaries are approximately 10-20 percent lower than those in detached and semi-detached systems. We calculate the binary flares activity level (AL) of all the flare binaries, and discuss its variations along the orbital period (P_orb) and rotation period (P_rot, calculated for only detached binaries). We find that AL increases with decreasing P_orb or P_rot up to the critical values at P_orb near 3 days or P_rot near 1.5 days, thereafter, the AL starts decreasing no matter how fast the stars rotate. We examine the flaring rate as a function of orbital phase in 2 eclipsing binaries on which a large number of flares are detected. It appears that there is no correlation between flaring rate and orbital phase in these 2 binaries. In contrast, when we examine the function with 203 flares on 20 non-eclipse ellipsoidal binaries, bimodal distribution of amplitude weighted flare numbers shows up at orbital phase 0.25 and 0.75. Such variation could be larger than what is expected from the cross-section modification.

In gamma ray astronomy, the energy range from sub-100GeV to TeV is crucial due to where there is a gap between space experiments and ground-based ones. In addition, observations in this energy range are expected to provide more details about the high energy emission from GRBs,and thus to understand EBL. Based on the observation results and the related knowledge, scientists may be able to unveil the mysteries of galaxy formation and the evolution of early universe. One of the principal issues for next generation Imaging Atmospheric Cherenkov Telescopes (IACT) is to achieve larger field of view (FoV). In this work, we report a refractive water convex lens as light collector to test the feasibility of a new generation of IACT, and some preliminary test results on the optical properties (the focal length, spot size, transmittance, etc.) of a 0.9 m diameter water lens, the photodetectors and DAQ system of a prototype are presented and discussed.

To reconcile the BICEP2 measurement on the tensor-to-scalar ratio $r$ with Planck constraint, a large negative running of scalar spectral index $n_s$ is needed. So the inflationary observable such as $n_s$ should be expanded at least to the second-order slow-roll parameters for single-field inflationary models. The large value of $r$ and the Lyth bound indicate that it is impossible to obtain the sub-Planckian excursion for the inflaton. However, we derive an absolutely minimal bound $\Delta\phi/M_{\rm Pl}>\sqrt{r/2}$ on the inflaton excursion for single-field inflationary models, which can be applied to non-slow-roll inflationary models as well. This bound excludes the possibility of the small-field inflation with $\Delta\phi<0.1 M_{\rm Pl}$ if the BICEP2 result on $r$ stands and it opens the window of sub-Planckian excursion with $\Delta\phi<M_{\rm Pl}$ even if $r$ is as large as $0.1$. To get the sub-Planckian excursion with $\Delta\phi<0.1 M_{\rm Pl}$, our modified bound requires $r<0.02$. Using a fifth-order polynomial potential as an explicit example, we show that it not only agrees with the observational results, but also violates the Lyth bound.

The BICEP2 experiment confirms the existence of primordial gravitational wave with the tensor-to-scalar ratio $r=0$ ruled out at $7\sigma$ level. The consistency of this large value of $r$ with the {\em Planck} data requires a large negative running $n'_s$ of the scalar spectral index. Herein we propose two types of the single field inflation models with simple potentials to study the possibility of the consistency of the models with the BICEP2 and {\em Planck} observations. One type of model suggested herein is realized in the supergravity model building. These models fail to provide the needed $n'_s$ even though both can fit the tensor-to-scalar ratio and spectral index.

The detection of B-mode power spectrum by the BICEP2 collaboration constrains the tensor-to-scalar ratio $r=0.20^{+0.07}_{-0.05}$ for the lensed-$\Lambda$CDM model. The consistency of this big value with the {\em Planck} results requires a large running of the spectral index. The large values of the tensor-to-scalar ratio and the running of the spectral index put a challenge to single field inflation. For the chaotic inflation, the larger the value of the tensor-to-scalar ratio is, the smaller the value of the running of the spectral index is. For the natural inflation, the absolute value of the running of the spectral index has an upper limit.

Planck measurements of the cosmic microwave background power spectra find a lower value of the Hubble constant $H_0$ and a higher value of the fractional matter energy density $\Omega_{m0}$ for the concordance $\Lambda$CDM model, and these results are in tension with other measurements. The {\em Planck} group argued that the tension came either from some sources of unknown systematic errors in some astrophysical measurements or the wrong $\Lambda$CDM model applied in fitting the data. We studied the reason for the tension on $H_0$ from different measurements by considering two dynamical dark energy models. We found that there is no tension between different data, the constraint on $H_0$ is almost unchanged for different dark energy models and the tension with the local measurements remains when the error bar on $H_0$ is tightened to be around 1. We argue that the tension on $H_0$ is not caused by the fitting model.

The dynamics of scalar fields as dark energy is well approximated by some general relations between the equation of state parameter $w(z)$ and the fraction energy density $\Omega_\phi$. Based on the approximation, for slowly-rolling scalar fields, we derived the analytical expressions of $w(z)$ which reduce to the popular Chevallier-Polarski-Linder parametrization with explicit degeneracy relation between $w_0$ and $w_a$. The models approximate the dynamics of scalar fields well and help eliminate the degeneracies among $w_a$, $w_0$ and $\Omega_{m0}$. With the explicit degeneracy relations, we test their effects on the constraints of cosmological parameters. We find that: (1) The analytical relations between $w_0$ and $w_a$ for the two models are consistent with observational data; (2) The degeneracies have little effect on $\Omega_{m0}$; (3) The $1\sigma$ error of $w_0$ was reduced about 30% with the degeneracy relations.

We consider a dark energy model with a relation between the equation of state parameter $w$ and the energy density parameter $\Omega_\phi$ derived from thawing scalar field models. Assuming the variation of the fine structure constant is caused by dark energy, we use the observational data of the variation of the fine structure constant to constrain the current value of $w_0$ and $\Omega_{\phi 0}$ for the dark energy model. At the $1\sigma$ level, the observational data excluded some areas around $w_0=-1$, which explains the positive detection of the variation of the fine structure constant at the $1\sigma$ level, but $\Lambda$CDM model is consistent with the data at the $2\sigma$ level.

The constraints on the $\Lambda$CDM model from type Ia supernova (SNe Ia) data alone and BAO data alone are similar, so it is worthwhile to compare their constraints on the property of dark energy. We apply the SNLS3 compilation of 472 SNe Ia data, the baryon acoustic oscillation measurement of distance, the cosmic microwave background radiation data from the WMAP7, and the Hubble parameter data to study the effect of their different combinations on the fittings of cosmological parameters in the modified holographic dark energy model and the Chevallier-Polarski-Linder model. Neither BAO nor WMAP7 data alone give good constraint on the equation of state parameter of dark energy, but both WMAP7 data and BAO data help SNe Ia data break the degeneracies among the model parameters, hence tighten the constraint on the variation of equation of state parameter $w_a$, and WMAP7 data do the job a little better. Although BAO and WMAP7 data provide reasonably good constraints on $\Omega_m$ and $\Omega_k$, they are not able to constrain the dynamics of dark energy. On the other hand, SNe Ia data do not provide good constraints on $\Omega_m$ and $\Omega_k$, but they provide good constraint on the dynamics of dark energy, especially the variation of the equation of state parameter of dark energy, so we need to combine SNe Ia with BAO and WMAP7 data to probe the property of dark energy. The addition of $H(z)$ data helps better constrain the geometry of the Universe $\Omega_k$ and the property of dark energy. For the SNLS SNe Ia data, the nuisance parameters $\alpha$ and $\beta$ are consistent for all different combinations of the above data, and their impacts on the fittings of cosmological parameters are minimal. By fitting the data to different models, $\Lambda$CDM model is still consistent with all the observational data.