In this work, we considered 2 schemes (a high-rigidity break in primary source injections and a high-rigidity break in diffusion coefficient) to reproduce the newly released AMS-02 nuclei spectra (He, C, N, O, Li, Be, and B) when the rigidity larger than 50 GV. The fitting results show that current data set favors a high-rigidity break at $\sim 325 \mathrm{GV}$ in diffusion coefficient rather than a break at $\sim 365 \mathrm{GV}$ in primary source injections. Meanwhile, the fitted values of the factors to rescale the cosmic-ray (CR) flux of secondary species/components after propagation show us that the secondary flux are underestimated in current propagation model. It implies that we might locate in a slow diffusion zone, in which the CRs propagate with a small value of diffusion coefficient compared with the averaged value in the galaxy. Another hint from the fitting results show that extra secondary CR nuclei injection may be needed in current data set. All these new hints should be paid more attention in future research.

In Phys. Rev. D 98, 103023 (2018), a novel scenario was proposed to probe the interactions between dark matter (DM) particles and electrons, via hydrogen-atmosphere pulsating white dwarfs (DAVs) in globular clusters. The estimation showed that the scenario could hopefully test the parameter space: $5 \mathrm{GeV} \le m_{\chi} \le 10^{4} \mathrm{GeV}$ and $\sigma_{\chi,e} \ge 10^{-40} \mathrm{cm}^{2}$, where $m_{\chi}$ is the DM particle's mass and $\sigma_{\chi,e}$ is the elastic scattering cross section between DM and electron. In this comment, we have determined the exact lower limit of the testable DM particle mass $\sim 1.38 - 1.58 \mathrm{GeV}$, which depends on $\sigma_{\chi,e}$. This gives us a credible lower limit of the testable DM particle mass in above scenario, and provide a clear upper limit of the DM particle mass which we should consider in future research.

An increasing number of white dwarf stars show atmospheric chemical composition polluted by heavy elements accreted from debris disk material. The existence of such debris disks strongly suggests the presence of one or more planet(s) whose gravitational interaction with rocky planetesimals is responsible for their disruption by tidal effect. The ZZ Ceti pulsator and polluted DAZ white dwarf GD 133 is a good candidate for searching for such a potential planet. We started in 2011 a photometric follow-up of its pulsations. As a result of this work in progress, we used the data gathered from 2011 to 2015 to make an asteroseismological analysis of GD 133, providing the star parameters from a best fit model with $M$/$M_{\odot}$ = 0.630 $\pm$ 0.002, $T_{\rm eff}$ = 12400 K $\pm$ 70 K, log($M_{\rm He}/M$) = -2.00 $\pm$ 0.02, log($M_{\rm H}/M$) = -4.50 $\pm$ 0.02 and determining a rotation period of $\approx$ 7 days.

Based on the precise nuclei data released by AMS-02, we study the spectra hardening of both the primary (proton, helium, carbon, oxygen, and the primary component of nitrogen) and the secondary (anti-proton, lithium, beryllium, boron and the secondary component of nitrogen) cosmic ray (CR) nuclei. With the diffusion-reacceleration model, we consider two schemes to reproduce the hardening in the spectra: (i) A high-rigidity break in primary source injection; (ii) A high-rigidity break in diffusion coefficient. The global fitting results show that both schemes could reproduce the spectra hardening in current status. More precise multi-TV data (especially the data of secondary CR species) is needed if one wants to distinguish these two schemes. In our global fitting, each of the nuclei species is allocated an independent solar modulation potential and a re-scale factor (which accounts for the isotopic abundance for primary nuclei species and uncertainties of production cross section or inhomogeneity of CR sources and propagation for secondary nuclei species). The fitting values of these two parameter classes show us some hints on some new directions in CR physics. All the fitted re-scale factors of primary nuclei species have values that systematically smaller than 1.0, while that of secondary nuclei species are systematically larger than 1.0. Moreover, both the re-scale factor and solar modulation potential of beryllium have values which are obviously different from other species. This might indicate that beryllium has the specificity not only on its propagation in the heliosphere, but also on its production cross section. All these new results should be seriously studied in the future.

In this work, we propose a novel scenario to probe the interactions between dark matter (DM) particles and electrons, via hydrogen-atmosphere pulsating white dwarfs (DAVs) in globular clusters. In this special configuration, the DM particles, which are predominantly captured by multiple scattering with the electrons in a DAV, would annihilate by pairs and provide extra energy source to the DAV. This mechanism slows down the natural cooling evolution which can be presented by the period variation rates of pulsation modes. The differences between the secular rates predicted by the precise asteroseismology and the secular rates obtained from observation can reveal the DM-electron interactions. An important observable has been proposed and corresponding estimations have been made. According to the estimation, if this scenario could be implemented in the near future, the potential sensitivity on $m_{\chi}$ (DM particle's mass) and $\sigma_{\chi,e}$ (elastic scattering cross section between DM and electron) could be hopefully extended to a region $5 \mathrm{GeV} \lesssim m_{\chi} \lesssim 10^{4} \mathrm{GeV}$ and $\sigma_{\chi,e} \gtrsim 10^{-40} \mathrm{cm}^{2}$. Combining with indirect DM detection results, this could give us a cross check on the existence of such leptonphilic DM particles to some extent.

In this work, we have studied the variable stars in the young open cluster NGC 1893 based on a multi-year photometric survey covering a sky area around the cluster up to $31' \times 31'$ wide. More than 23\,000 images in the $V$ band taken from January 2008 to February 2017 with different telescopes, complemented with 90 images in the $B$ band in 2014 and 2017, were reduced, and light curves were derived in $V$ for 5653 stars. By analyzing these light curves, we detected 147 variable stars (85 of them being new discoveries), including 110 periodic variables, 15 eclipsing binaries and 22 non-periodic variables. Proper motions, radial velocities, color-magnitude and two-color diagrams were used to identify the cluster membership of these variable stars, resulting in 84 members. Periodic variable members were then classified into different variability types, mainly according to their magnitudes and to their periods of variability, as well as to their positions in the Hertzsprung-Russell diagram for the early-type stars. As a result, among main-sequence periodic variable members, we identified five $\beta$ Cep candidates, seven slowly pulsating B-type candidates, and thirteen fast-rotating pulsating B-type (FaRPB) candidates (one of which is a confirmed classical Be star). While most of the FaRPB stars display properties similar to the ones discovered in NGC 3766 by Mowlavi et al. (2013), five of them have periods below 0.1~d, contrary to expectations. Additional observations, including spectroscopic, are called for to further characterize these stars. We also find a binary candidate harboring a $\delta$-Scuti candidate.

Bi-site time-series photometric and high-resolution spectroscopic observations were made for the double-mode high-amplitude $\delta$ Scuti star VX\ Hya. The fundamental frequency $f_{0}=4.4763\ \rm{c\ days^{-1}}$, the first overtone $f_{1}=5.7897\ \rm{c\ days^{-1}}$ and 23 harmonics and linear combinations of $f_{0}$ and $f_{1}$ are detected by pulsation analysis. From the spectroscopic data, we get $\rm{[Fe/H] = -0.2\pm0.1\ dex}$. The period change rate of the fundamental mode is obtained by using the Fourier-phase diagram method, which gives the value of $(1/P_{0})(dP_{0}/dt)=(1.81\pm0.09) \times 10^{-7}\ \rm{yr}^{-1}$. With these results from the observations, we perform theoretical explorations with the stellar evolution code MESA, and constrain the models by fitting $f_{0}$, $f_{1}$, and $(1/P_{0})(dP_{0}/dt)$ within $3\sigma$ deviations. The results show that the period change of VX\ Hya could be ascribed to the evolutionary effect. The stellar parameters of VX\ Hya could be derived as: the mass $2.385\pm0.025\ M_{\odot}$, the luminosity $\log(L/L_{\odot})=1.93\pm0.02$ and the age $(4.43\pm0.13)\times 10^8$ years. VX\ Hya is found to locate at the post-main-sequence stage with a helium core and a hydrogen-burning shell on the H${-}$R diagram.

Recently, DAMPE has released its first results on the high-energy cosmic-ray electrons and positrons (CREs) from about $25$ GeV to $4.6$ TeV, which directly detect a break at $\sim 1$ TeV. This result gives us an excellent opportunity to study the source of the CREs excess. In this work, we used the data fo proton and helium flux (from AMS-02 and CREAM), $\bar{\mathrm{p}}/\mathrm{p}$ ratio (from AMS-02), positron flux (from AMS-02) and CREs flux (from DAMPE without the peak signal point at $\sim 1.4$ TeV) to do global fitting simultaneously, which can account the influence from the propagation model, the nuclei and electron primary source injection and the secondary lepton production precisely. For extra source to interpret the excess in lepton spectrum, we consider two separate scenarios (pulsar and dark matter annihilation via leptonic channels) to construct the bump ($\gtrsim 100$ GeV) and the break at $\sim 1$ TeV. The result shows: (i) in pulsar scenario, the spectral index of the injection should be $\nu_{\mathrm{psr}} \sim 0.65$ and the cut-off should be $R_{c} \sim 650$ GV; (ii) in dark matter scenario, the dark matter particle's mass is $m_{\chi} \sim 1208$ GeV and the cross section is $\langle \sigma v \rangle \sim 1.48 \times 10^{-23} \mathrm{cm}^{3} \mathrm{s}^{-1}$. Moreover, in the dark matter scenario, the $\tau \bar{\tau}$ annihilation channel is highly suppressed, and a DM model is built to satisfy the fitting results.

The variability of young stellar objects is mostly driven by star-disk interactions. In long-term photometric monitoring of the accreting T Tauri star GI Tau, we detect extinction events with typical depths of $\Delta V \sim 2.5$ mag that last for days-to-months and often appear to occur stochastically. In 2014 - 2015, extinctions that repeated with a quasi-period of 21 days over several months is the first empirical evidence of slow warps predicted from MHD simulations to form at a few stellar radii away from the central star. The reddening is consistent with $R_V=3.85\pm0.5$ and, along with an absence of diffuse interstellar bands, indicates that some dust processing has occurred in the disk. The 2015 -- 2016 multi-band lightcurve includes variations in spot coverage, extinction, and accretion, each of which results in different traces in color-magnitude diagrams. This lightcurve is initially dominated by a month-long extinction event and return to the unocculted brightness. The subsequent light-curve then features spot modulation with a 7.03 day period, punctuated by brief, randomly-spaced extinction events. The accretion rate measured from $U$-band photometry ranges from $1.3\times10^{-8}$ to $1.1\times10^{-10}$ M$_\odot$ yr$^{-1}$ (excluding the highest and lowest 5% of high- and low- accretion rate outliers), with an average of $4.7 \times 10^{-9}$ M$_\odot$ yr$^{-1}$. A total of 50% of the mass is accreted during bursts of $>12.8\times10^{-9}$ M$_\odot$ yr${^{-1}}$, which indicates limitations on analyses of disk evolution using single-epoch accretion rates.

We propose a novel and feasible method to detect dark matter (DM) electron interaction via pulsating white dwarfs (WDs) in the central region of globular clusters. Annihilation of the DM particles captured by those WDs can provide additional energy source along the natural cooling evolution of WDs and the cooling print can be well offered by precise asteroseismology. The measurement of the long time scale physical quantity -- the rates of period variation of pulsation modes -- could be used to constrain the cross section between DM particles and electrons ($\sigma_{e}$), when DM particle mass $m_{\chi} \gtrsim 5 \mathrm{GeV}$. We construct estimations to prove that this original method is feasible and can be implemented in the challenging time-series photometry in the near future.

We analyze the photometric data and spectroscopic data that collect on the $\delta$ Scuti star AE UMa. The fundamental and the first overtone frequencies are confirmed as $f_{0} = 11.62560$\,c\,d$^{-1}$ and $f_{1} = 15.03124$\,c\,d$^{-1}$, respectively, from the frequency content by analyzing of the 40 nights light curve spanning from 2009 to 2012. Additionally, another 37 frequencies are identified as either the harmonics or the linear combinations of the fundamental and the first overtone frequencies, among which 25 are newly detected. The rate of period change of the fundamental mode is determined as $(1/P_{0})(dP_{0}/dt) = 5.4(\pm 1.9) \times 10^{-9}$\,yr$^{-1}$ as revealed from the $O-C$ diagram based on the 84 newly determined times of maximum light combined with those derived from the literature. The spectroscopic data suggests that AE UMa is a population\,I $\delta$ Scuti star. With these physical properties, we perform theoretical explorations based on the stellar evolution code MESA on this target, considering that the variation of pulsation period is caused by secular evolutionary effects. We finally constraint the AE UMa with the physical parameters as: the mass of $1.805 \pm 0.055\ M_{\odot}$, the radius of $1.647 \pm 0.032 \times 10^{11}$\,cm, the luminosity of $1.381 \pm 0.048~(\log L/L_{\odot})$ and the age of $1.055 \pm 0.095 \times 10^{9}$\,yr. AE UMa can be the (Pop. I) $\delta$ Scuti star that locates just after the second turn-off of its evolutional track leaving the main sequence, a star in the phase of the Hertzsprung Gap with a helium core and a hydrogen-burning shell.

O-C diagram is a useful technique to analyse the period changes of a pulsator by using the maximum (or minimum) value points which have been obtained from the historical data. But if an object is a double-mode or multi-mode pulsator, the extreme value points are the results of all the modes other than just the fundamental mode. We discussed these situations and give out some criteria to judge whether the O-C diagram can be used in these situations.