RV variable stars are important in astrophysics. The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) spectroscopic survey has provided ~ 6.5 million stellar spectra in its Data Release 4 (DR4). During the survey, ~ 4.7 million unique sources were targeted and ~ 1 million stars observed repeatedly. The probabilities of stars being RV variables are estimated by comparing the observed radial velocity variations with the simulated ones. We build a catalog of 80,702 RV variable candidates with probability greater than 0.60 by analyzing the duplicate-observed multi-epoch sources covered by the LAMOST DR4. Simulations and cross-identifications show that the purity of the catalog is higher than 80%. The catalog consists of 77% binary systems and 7% pulsating stars as well as 16% pollution by single stars. 3,138 RV variables are classified through cross-identifications with published results in literatures. By using the 3,138 sources common to both LAMOST and a collection of published RV variable catalogs we are able to analyze LAMOST's RV variable detection rate. The efficiency of the method adopted in this work relies not only on the sampling frequency of observations but also periods and amplitudes of RV variables. With the progress of LAMOST, Gaia and other surveys, more and more RV variables would will be confirmed and classified. This catalog is valuable for other large-scale surveys, especially for RV variable searches. The catalog will be released according to the LAMOST Data Policy via http://dr4.lamost.org.

Solar system materials are variably depleted in moderately volatile elements (MVEs) relative to the proto-solar composition. To address the origin of this MVE depletion, we conducted a systematic study of high-precision K isotopic composition on 16 carbonaceous chondrites (CCs) of types CM1-2, CO3, CV3, CR2, CK4-5 and CH3 and 28 ordinary chondrites (OCs) covering petrological types 3 to 6 and chemical groups H, L, and LL. We observed significant overall K isotope (delta41K) variations (-1.54 to 0.70 permil). The K isotope compositions of CCs are largely higher than the Bulk Silicate Earth (BSE) value, whereas OCs show typically lower values than BSE. Neither CCs nor OCs show resolvable correlations between K isotopes and chemical groups, petrological types, shock levels, exposure ages, fall or find occurrence, or terrestrial weathering. The lack of a clear trend between K isotopes and K content indicates that the K isotope fractionations were decoupled from the relative elemental K depletions. The range of K isotope variations in the CCs is consistent with a four-component (chondrule, refractory inclusion, matrix and water) mixing model that is able to explain the bulk elemental and isotopic compositions of the main CC groups, but requires a fractionation in K isotopic compositions in chondrules. We propose that the major control of the isotopic compositions of group averages is condensation or vaporization in nebular environments that is preserved in the compositional variation of chondrules. Parent-body processes (aqueous alteration, thermal metamorphism, and metasomatism) can mobilize K and affect the K isotopes in individual samples. In the case of the OCs, the full range of K isotopic variations can only be explained by the combined effects of the size and relative abundances of chondrules, parent-body aqueous and thermal alteration.

Enstatite chondrites and aubrites are meteorites that show the closest similarities to the Earth in many isotope systems that undergo mass-independent and mass-dependent isotope fractionations. Due to the analytical challenges to obtain high-precision K isotopic compositions in the past, potential differences in K isotopic compositions between enstatite meteorites and the Earth remained uncertain. We report the first high-precision K isotopic compositions of eight enstatite chondrites and four aubrites and find that there is a significant variation of K isotopic compositions among enstatite meteorites (from -2.34 permil to -0.18 permil). However, K isotopic compositions of nearly all enstatite meteorites scatter around the Bulk Silicate Earth (BSE) value. The average K isotopic composition of the eight enstatite chondrites (-0.47 +/- 0.57 permil) is indistinguishable from the BSE value (-0.48 +/- 0.03 permil), thus further corroborating the isotopic similarity between Earth' building blocks and enstatite meteorite precursors. We found no correlation of K isotopic compositions with the chemical groups, petrological types, shock degrees, and terrestrial weathering conditions; however, the variation of K isotopes among enstatite meteorite can be attributed to the parent body processing. Our sample of the main group aubrite MIL 13004 is exceptional and has an extremely light K isotopic composition (delta 41K= -2.34 +/- 0.12 permil). We attribute this unique K isotopic feature to the presence of abundant djerfisherite inclusions in our sample because this K-bearing sulfide mineral is predicted to be enriched in 39K during equilibrium exchange with silicates.

Monitoring chromospheric and photospheric indexes of magnetic activity can provide valuable information, especially the interaction between different parts of the atmosphere and their response to magnetic fields. We extract chromospheric indexes, S and Rhk+, for 59,816 stars from LAMOST spectra in the LAMOST-Kepler program, and photospheric index, Reff, for 5575 stars from Kepler light curves. The log Reff shows positive correlation with log Rhk+. We estimate the power-law indexes between Reff and Rhk+ for F-, G-, and K-type stars, respectively. We also confirm the dependence of both chromospheric and photospheric activity on stellar rotation. Ca II H and K emissions and photospheric variations generally decrease with increasing rotation periods for stars with rotation periods exceeding a few days. The power-law indexes in exponential decay regimes show different characteristics in the two activity-rotation relations. The updated largest sample including the activity proxies and reported rotation periods provides more information to understand the magnetic activity for cool stars.

We present a catalog of 5,290 RR Lyrae stars (RRLs) with metallicities estimated from spectra of the LAMOST Experiment for Galactic Understanding and Exploration (LEGUE) and the Sloan Extension for Galactic Understanding and Exploration (SEGUE) surveys. Nearly 70 per cent of them (3,642 objects) also have systemic radial velocities measured. Given the pulsating nature of RRLs, metallicity estimates are based on spectra of individual exposures, by matching them with the synthetic templates. The systemic radial velocities are measured by fitting the observed velocity as a function of phase assuming an empirical pulsating velocity template curve. Various tests show that our analyses yield metallicities with a typical precision of 0.20\,dex and systemic radial velocities with uncertainties ranging from 5 to 21\,km\,s$^{-1}$ (depending on the number of radial velocity measurements available for a given star). Based on the well calibrated near-infrared $PM_{W1}Z$ or $PM_{K_{\rm s}}Z$, and $M_{V}$-[Fe/H] relations, precise distances are derived for these RRLs. Finally, we include Gaia DR2 proper motions in our catalog. The catalog should be very useful for various Galactic studies, especially of the Galactic halo.

We present a large, homogeneous catalogue of molecular clouds within 4 kpc from the Sun at low Galactic latitudes ($|b|$ $<$ 10\degr) with unprecedented accurate distance determinations. Based on the three-dimensional dust reddening map and estimates of colour excesses and distances of over 32 million stars presented in Chen et al, we have identified 567 dust/molecular clouds with a hierarchical structure identification method and obtained their distance estimates by a dust model fitting algorithm. The typical distance uncertainty is less than 5 per cent. As far as we know, this is the first large catalogue of molecular clouds in the Galactic plane with distances derived in a direct manner. The clouds are seen to lie along the Sagittarius, Local and Perseus Arms. In addition to the known structures, we propose the existence of a possible {\it spur}, with a pitch angle of about 34\degr, connecting the Local and the Sagittarius Arms in the fourth quadrant. We have also derived the physical properties of those molecular clouds. The distribution of cloud properties in different parameter spaces agrees grossly with the previous results. Our cloud sample is an ideal starting point to study the concentration of dust and gas in the solar vicinity and their star formation activities.

According to \emph{Wind} observations between June 2004 and May 2019, this Letter investigates the proton and alpha particle temperatures in the space of ($\theta_d$, $V_d/V_A$) for the first time, where $\theta_d$ and $V_d$ are the radial angle and magnitude of alpha$-$proton differential flow vector respectively, $V_A$ is the local Alfv\'en speed. Results show that the temperatures significantly depend on $\theta_d$ as well as $V_d/V_A$. In case of low proton parallel beta ($\beta_{p{\parallel}} < 1$), it is found that the proton perpendicular temperature is clearly enhanced when $\theta_d$ is small ($\lesssim 45^\circ$) and $V_d/V_A \gtrsim 0.5$. On the contrary, the perpendicular temperature of alpha particles is considerably enhanced when $\theta_d$ is large ($\gtrsim 90^\circ$) or $V_d/V_A$ is sufficiently small. The maximum of proton parallel temperature takes place at $\theta_d \sim 90^\circ$ accompanied by higher $\beta_{p{\parallel}}$ and by larger turbulence amplitude of magnetic fluctuations in inertial range. This study should present strong evidence for cyclotron resonance heating of protons and alpha particles in the solar wind. Other mechanisms including Landau resonance and stochastic heating are also proposed, which tend to have different ($\theta_d$, $V_d/V_A$) spaces than cyclotron resonance heating.

The hot massive luminous blue variables (LBVs) represent an important evolutionary phase of massive stars. Here, we report the discovery of a new LBV -- LAMOST J0037+4016 in the distant outskirt of the Andromeda galaxy. It is located in the south-western corner (a possible faint spiral arm) of M31 with an unexpectedly large projection distance of $\sim$ 22 kpc from the center. The optical light curve shows a 1.2 mag variation in $V$ band and its outburst and quiescence phases both last over several years. The observed spectra indicate an A-type supergiant at epoch close to the outburst phase and a hot B-type supergiant with weak [Fe II] emission lines at epoch of much dimmer brightness. The near-infrared color-color diagram further shows it follows the distribution of Galactic and M31 LBVs rather than B[e] supergiants. All the existing data strongly show that LAMOST J0037+4016 is an LBV. By spectral energy distribution fitting, we find it has a luminosity ($4.42 \pm 1.64$)$\times 10^5$ $L_{\odot}$ and an initial mass $\sim 30$ $M_{\odot}$, indicating its nature of less luminosity class of LBV.

In this paper,we report a detailed analysis of recurrent jets originated from a location with emerging, canceling and converging negative magnetic field at the east edge of NOAA active region AR11166 from 2011 March 09 to 10. The event presented several interesting features. First, a satellite sunspot appeared and collided with a pre-existing opposite polarity magnetic field and caused a recurrent solar jet event. Second, the evolution of the jets showed blowout-like nature and standard characteristics. Third, the satellite sunspot exhibited a motion toward southeast of AR11166 and merged with the emerging flux near the opposite polarity sunspot penumbra, which afterward, due to flux convergence and cancellation episodes, caused recurrent jets. Fourth, three of the blowout jets associated with coronal mass ejections (CMEs), were observed from field of view of the Solar Terrestrial Relations Observatory. Fifth, almost all the blowout jet eruptions were accompanied with flares or with more intense brightening in the jet base region, while almost standard jets did not manifest such obvious feature during eruptions. The most important, the blowout jets were inclined to faster and larger scale than the standard jets. The standard jets instead were inclined to relative longer-lasting. The obvious shearing and twisting motions of the magnetic field may be interpreted as due to the shearing and twisting motions for a blowout jet eruption. From the statistical results, about 30% blowout jets directly developed into CMEs. It suggests that the blowout jets and CMEs should have a tight relationship.

The spatial-dependent propagation model has been successfully used to explain diverse observational phenomena, including the spectral hardening of cosmic-ray nuclei above $200$ GV, the large-scale dipole anisotropy and the diffusive gamma distribution. In this work, we further apply the spatial-dependent propagation model to both electrons and positrons. To account for the excess of positrons above $10$ GeV, an additional local source is introduced. And we also consider a more realistic spiral distribution of background sources. We find that due to the gradual hardening above $10$ GeV, the hardening of electron spectrum above tens of GeV can be explained in the SDP model and both positron and electron spectra less than TeV energies could be naturally described. The spatial-dependent propagation with spiral-distributed sources could conforms with the total electron spectrum in the whole energy. Meanwhile compared with the conventional model, the spatial-dependent propagation with spiral-distributed sources could produce larger background positron flux, so that the multiplier of background positron flux is $1.42$, which is much smaller than the required value by the conventional model. Thus the shortage of background positron flux could be solved. Furthermore we compute the anisotropy of electron under spatial-dependent propagation model, which is well below the observational limit of Fermi-LAT experiment.

With the photometric data from the SDSS survey, the spectroscopic data from the SDSS/SEGUE and the LAMOST surveys, and the astrometric data from the Gaia DR2, we have identified 71 highly-probable member stars of the GD-1 cold stellar stream spread along almost its entire length (i.e. from 126 to 203 degree in Right Ascension). With the accurate spectroscopic (i.e. metallicity and line-of-sight velocity) and astrometric (i.e. proper motions) information, the position-velocity diagrams, i.e. $\phi_{1}$-$\mu_{\alpha}$, $\phi_{1}$-$\mu_{\delta}$ and $\phi_{1}$-${v_{\rm gsr}}$, of the GD-1 stream are well mapped. The stream has an average metallicity [Fe/H] $= -1.95 \pm 0.23$. The rich information of member stars of the stream now available allow one not only to model its origin, but also to place strong constraints on the mass distribution and the gravitational potential of the Milky Way.

We have investigated the distributions of stellar azimuthal and radial velocity components $V_{\Phi}$ and $V_{R}$ in the vertical position-velocity plane $Z$-$V_{Z}$ across the Galactic disc of $6.34 \lesssim R \lesssim 12.34$\,kpc and $|\Phi| \lesssim 7.5^{\circ}$ using a Gaia and Gaia-LAMOST sample of stars. As found in previous works, the distributions exhibit significant spiral patterns. The $V_{R}$ distributions also show clear quadrupole patterns, which are the consequence of the well-known tilt of the velocity ellipsoid. The observed spiral and quadrupole patterns in the phase space plane vary strongly with radial and azimuthal positions. The phase spirals of $V_{\Phi}$ become more and more relaxed as $R$ increases. The spiral patterns of $V_{\Phi}$ and $V_{R}$ and the quadrupole patterns of $V_{R}$ are strongest at $-2^{\circ} < \Phi < 2^{\circ}$ but negligible at $4^{\circ} < \Phi < 6^{\circ}$ and $-6^{\circ} < \Phi < -4^{\circ}$. Our results suggest an external origin of the phase spirals. In this scenario, the intruder, most likely the previously well-known Sagittarius dwarf galaxy, passed through the Galactic plane in the direction towards either Galactic center or anti-center. The azimuthal variations of the phase spirals also help us constrain the passage duration of the intruder. A detailed model is required to reproduce the observed radial and azimuthal variations of the phase spirals of $V_{\Phi}$ and $V_{R}$.

We present a catalog of stellar age and mass estimates for a sample of 640\,986 red giant branch (RGB) stars of the Galactic disk from the LAMOST Galactic Spectroscopic Survey (DR4). The RGB stars are distinguished from the red clump stars utilizing period spacing derived from the spectra with a machine learning method based on kernel principal component analysis (KPCA). Cross-validation suggests our method is capable of distinguishing RC from RGB stars with only 2 per cent contamination rate for stars with signal-to-noise ratio (SNR) higher than 50. The age and mass of these RGB stars are determined from their LAMOST spectra with KPCA method by taking the LAMOST - $Kepler$ giant stars having asteroseismic parameters and the LAMOST-TGAS sub-giant stars based on isochrones as training sets. Examinations suggest that the age and mass estimates of our RGB sample stars with SNR $>$ 30 have a median error of 30 per cent and 10 per cent, respectively. Stellar ages are found to exhibit positive vertical and negative radial gradients across the disk, and the age structure of the disk is strongly flared across the whole disk of $6<R<13$\,kpc. The data set demonstrates good correlations among stellar age, [Fe/H] and [$\alpha$/Fe]. There are two separate sequences in the [Fe/H] -- [$\alpha$/Fe] plane: a high--$\alpha$ sequence with stars older than $\sim$\,8\,Gyr and a low--$\alpha$ sequence composed of stars with ages covering the whole range of possible ages of stars. We also examine relations between age and kinematic parameters derived from the Gaia DR2 parallax and proper motions. Both the median value and dispersion of the orbital eccentricity are found to increase with age. The vertical angular momentum is found to fairly smoothly decrease with age from 2 to 12\,Gyr, with a rate of about $-$50\,kpc\,km\,s$^{-1}$\,Gyr$^{-1}$. A full table of the catalog is public available online.

We present new three-dimensional (3D) interstellar dust reddening maps of the Galactic plane in three colours, E(G-Ks), E(Bp-Rp) and E(H-Ks). The maps have a spatial angular resolution of 6 arcmin and covers over 7000 deg$^2$ of the Galactic plane for Galactic longitude 0 deg $<$ $l$ $<$ 360 deg and latitude $|b|$ $<$ $10$ deg. The maps are constructed from robust parallax estimates from the Gaia Data Release 2 (Gaia DR2) combined with the high-quality optical photometry from the Gaia DR2 and the infrared photometry from the 2MASS and WISE surveys. We estimate the colour excesses, E(G-Ks), E(Bp-Rp) and E(H-Ks), of over 56 million stars with the machine learning algorithm Random Forest regression, using a training data set constructed from the large-scale spectroscopic surveys LAMOST, SEGUE and APOGEE. The results reveal the large-scale dust distribution in the Galactic disk, showing a number of features consistent with the earlier studies. The Galactic dust disk is clearly warped and show complex structures possibly spatially associated with the Sagittarius, Local and Perseus arms. We also provide the empirical extinction coefficients for the Gaia photometry that can be used to convert the colour excesses presented here to the line-of-sight extinction values in the Gaia photometric bands.

We have investigated the metallicity distributions of mono-age stellar populations across the disc of $6 \lesssim R \lesssim 12$\,kpc and $|Z| \lesssim 2$\,kpc using samples selected from the main-sequence turn-off and sub-giant (MSTO-SG) stars targeted by the LAMOST Galactic Spectroscopic surveys. Both the mean values and the profiles of the distributions exhibit significant variations with age and position. We confirm that the oldest ($>11$\,Gyr) stars have nearly flat radial [Fe/H] gradients at all heights above the disc but show negative vertical [Fe/H] gradients. For stars younger than 11\,Gyr, the radial [Fe/H] gradients flatten with $|Z|$, while the vertical [Fe/H] gradients flatten with $R$. Stars of 4--6\,Gyr exhibit steeper negative radial [Fe/H] gradients than those of either younger or older ages. Values of [$\alpha$/Fe] of mono-age stellar populations also show significant radial and vertical gradients, with patterns varying with age. The [Fe/H] distribution profiles of old ($>8$\,Gyr) stars vary little with $R$, while those of younger stars exhibit strong radial variations, probably a consequence of significant radial migration. The [$\alpha$/Fe] radial distribution profiles show opposite patterns of variations with age compared to those of [Fe/H]. We have also explored the impacts of stellar mixing by epicycle motions (blurring) on the [Fe/H] and [$\alpha$/Fe] distributions, and found that blurring mainly change the widths of the distribution profiles. Our results suggest that the disc may have experienced a complex assemblage history, in which both the "inside-out" and "upside-down" formation processes may have played an important role.

This study expands the coverage and improves the homogeneity of the distribution of MILES template stars in the parameter space, as well as extends the wavelength coverage of the template spectra to the far red beyond the Ca\,{\sc{ii}} triplet. To achieve this we have carried out a major observational campaign using the OMR long-slit spectrograph mounted on the NAOC 2.16\,m telescope and the YFOSC long-slit spectrograph mounted on the YNAO 2.4\,m telescope. The original sample is based on the MILES library, supplemented by 918 stars selected from PASTEL database. In total, 822 OMR and 1,324 YFOSC spectra, covering respectively the wavelength ranges $\lambda\lambda$3800--5180 and $\lambda\lambda$5150--9000, have been collected and reduced. The spectra have a mean resolution FWHM (full-width at half-maximum) of $\sim 3.3$\,{\AA} and are wavelength- and flux-calibrated to an accuracy of $\sim 20$\,km\,s$^{-1}$ and $\sim 5$\,per\,cent, respectively. The spectra are further corrected for systematic errors in the wavelength calibration to an accuracy of $\sim 4$\,km\,s$^{-1}$ by cross-correlating with the theoretical spectra. Almost all the spectra have an average signal to noise ratio (SNR) better than 100 per pixel. Combined with the MILES spectra, there are now 1,731, 1,542, 1,324 and 1,273 stars with spectra covering respectively $\lambda\lambda$3800--5180, $\lambda\lambda$3800--7500, $\lambda\lambda$5150--9000 and $\lambda\lambda$3800--9000. This paper describes our template star selection, the observation and data reduction, and presents the reduced spectra collected hitherto.

We present a detailed determination and analysis of 3D stellar mass distribution of the Galactic disk for mono-age populations using a sample of 0.93 million main-sequence turn-off and subgiant stars from the LAMOST Galactic Surveys. Our results show (1) all stellar populations younger than 10\,Gyr exhibit strong disk flaring, which is accompanied with a dumpy vertical density profile that is best described by a $sech^n$ function with index depending on both radius and age; (2) Asymmetries and wave-like oscillations are presented in both the radial and vertical direction, with strength varying with stellar populations; (3) As a contribution by the Local spiral arm, the mid-plane stellar mass density at solar radius but 400--800\,pc (3--6$^\circ$) away from the Sun in the azimuthal direction has a value of $0.0594\pm0.0008$\,$M_\odot$/pc$^3$, which is 0.0164\,$M_\odot$/pc$^3$ higher than previous estimates at the solar neighborhood. The result causes doubts on the current estimate of local dark matter density; (4) The radial distribution of surface mass density yields a disk scale length evolving from $\sim$4\,kpc for the young to $\sim$2\,kpc for the old populations. The overall population exhibits a disk scale length of $2.48\pm0.05$\,kpc, and a total stellar mass of $3.6(\pm0.1)\times10^{10}$\,$M_\odot$ assuming $R_{\odot}=8.0$\,kpc, and the value becomes $4.1(\pm0.1)\times10^{10}$\,$M_\odot$ if $R_{\odot}=8.3$\,kpc; (5) The disk has a peak star formation rate ({\rm SFR}) changing from 6--8\,Gyr at the inner to 4--6\,Gyr ago at the outer part, indicating an inside-out assemblage history. The 0--1\,Gyr population yields a recent disk total {\rm SFR} of $1.96\pm0.12$\,$M_\odot$/yr.

Using a sample of about 123,000 stars with accurate 3D velocity measurements from the LAMOST-TGAS data, we confirm the kinematic signature of the Galactic warp found by Sch\"onrich & Dehnen recently. The data reveal a clear trend of increasing mean vertical velocity $\overline{V_{z}}$ as a function of absolute vertical angular momentum $L_{z}$ and azimuthal velocity $V_{\phi}$ for guiding center radius $R_{g}$ between 6.0 and 10.5 kpc. The trend is consistent with a large-scale Galactic warp. Similar to Sch\"onrich & Dehnen, we also find a wave-like pattern of $\overline{V_{z}}$ versus $L_{z}$ with an amplitude of $\sim 0.9$ km s$^{-1}$ on a scale of $\sim 2.0$ kpc, which could arise from bending waves or a winding warp. Finally, we confirm a prominent, localized peak in $\overline{V_z}$ near $L_z \sim 2150$ kpc km s$^{-1}$ (corresponding to $R_{g} \sim 9$ kpc and $V_{\phi} \sim 255$ km s$^{-1}$). The additional line-of-sight velocity information from LAMOST reveals that stars in this feature have a large, inward radial velocity of $V_{R} \sim -13.33 \pm 0.59$ km s$^{-1}$ and a small radial velocity dispersion $\sigma_{R} \sim 25.27 \pm 0.89$ km s$^{-1}$, suggesting that a stellar stream gives rise to this feature.

We present a new catalogue of 18 080 radial velocity standard stars selected from the APOGEE data. These RV standard stars are observed at least three times and have a median stability ($3\sigma_{\rm RV}$) around 240 m s$^{-1}$ over a time baseline longer than 200 days. They are largely distributed in the northern sky and could be extended to the southern sky by the future APOGEE-2 survey. Most of the stars are red giants ($J - K_{\rm s} \ge 0.5$) owing to the APOGEE target selection criteria. Only about ten per cent of them are main-sequence stars. The $H$ band magnitude range of the stars is 7-12.5 mag with the faint limit much fainter than the magnitudes of previous RV standard stars. As an application, we show the new set of standard stars to determine the radial velocity zero points of the RAVE, the LAMOST {and the Gaia-RVS} Galactic spectroscopic surveys.

We present the observations of a double-decker filament to study its formation, triggering, and eruption physics. It is observed that the double-decker filament was formed by splitting of an original single filament. During the splitting process, intermittent bright point bursts are observed in the filament channel, which resulted in the generation of the upper filament branch. The eruption of the newly formed double-decker filament was possibly triggered by two recurrent two-sided loop jets in the filament channel and the continuous mass unloading from the upper filament body. The interaction between the first jet and the filament directly resulted in the unstable of the lower branch and the fast rising phase of the upper branch. The second jet occurred at the same site about three hours after the first one, which further disturbed and accelerated the rising of the lower filament branch. It is interesting that the rising lower branch overtook the upper one, and then the two branches probably merged into one filament. Finally, the whole filament erupted violently and caused a large-scale coronal mass ejection, leaving behind a pair of flare ribbons and two dimming regions on the both sides of the filament channel. We think that the intermittent bursts may directly result in the rearrangement of the filament magnetic field and therefore the formation of the double-decker filament, then the recurrent jets further caused the fully eruption of the entire filament system. The study provides convincing evidence for supporting the scenario that a double-decker filament can be formed by splitting a single filament into two branches.

We present extensive spectroscopic observations of supernova remnant (SNR) S147 collected with the Large sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST). The spectra were carefully sky-subtracted taking into account the complex filamentary structure of S147. We have utilized all available LAMOST spectra toward S147, including sky and stellar spectra. By measuring the prominent optical emission lines including H{\alpha}, [NII]{\lambda}6584, and [SII]{\lambda}{\lambda}6717,6731, we present maps of radial velocity and line intensity ratio covering the whole nebula of S147 with unprecedented detail. The maps spatially correlated well with the complex filamentary structure of S147. For the central 2 deg of S147, the radial velocity varies from -100 to 100 km/s and peaks between ~ 0 and 10 km/s. The intensity ratios of H{\alpha}/[SII]{\lambda}{\lambda}6717,6731, [SII]{\lambda}6717/{\lambda}6731 and H{\alpha}/[NII]{\lambda}6584 peak at about 0.77, 1.35 and 1.48, respectively, with a scatter of 0.17, 0.19 and 0.37, respectively. The intensity ratios are consistent with the literature values. However, the range of variations of line intensity ratios estimated here and representative of the whole nebula, are larger than previously estimated.

Stellar systems composed of single, double, triple or high-order systems are rightfully regarded as the fundamental building blocks of the Milky Way. Binary stars play an important role in formation and evolution of the Galaxy. Through comparing the radial velocity variations from multi-epoch observations, we analyze the binary fraction of dwarf stars observed with the LAMOST. Effects of different model assumptions such as orbital period distributions on the estimate of binary fractions, are investigated. The results based on log-normal distribution of orbital periods reproduce the previous complete analyses better than the power-law distribution. We find that the binary fraction increases with $T_{\rm eff}$ and decreases with [Fe/H]. We first investigate the relation between $\alpha$-elements and binary fraction in such a large sample as the LAMOST. The old stars with high [$\alpha$/Fe] dominate higher binary fraction than young stars with low [$\alpha$/Fe]. At the same mass, former forming stars possess a higher binary fraction than newly forming ones, which may be related with the evolution of the Galaxy.

We present a detailed analysis of the selection function of the LAMOST Spectroscopic Survey of the Galactic Anti-centre (LSS-GAC). LSS-GAC was designed to obtain low resolution optical spectra for a sample of more than 3 million stars in the Galactic anti-centre. The second release of value-added catalogues of the LSS-GAC (LSS-GAC DR2) contains stellar parameters, including radial velocity, atmospheric parameters, elemental abundances and absolute magnitudes deduced from 1.8 million spectra of 1.4 million unique stars targeted by the LSS-GAC between 2011 and 2014. For many studies using this database, such as those investigating the chemodynamical structure of the Milky Way, a detailed understanding of the selection function of the survey is indispensable. In this paper, we describe how the selection function of the LSS-GAC can be evaluated to sufficient detail and provide selection function corrections for all spectroscopic measurements with reliable parameters released in LSS-GAC DR2. The results, to be released as new entries in the LSS-GAC value-added catalogues, can be used to correct the selection effects of the catalogue for scientific studies of various purposes.

The sky brightness is a critical parameter for estimating the coronal observation conditions for solar observatory. As part of a site-survey project in Western China, we measured the sky brightness continuously at the Lijiang Observatory in Yunnan province in 2011. A sky brightness monitor (SBM) was adopted to measure the sky brightness in a region extending from 4.5 to 7.0 apparent solar radii based on the experience of the Daniel K. Inouye Solar Tele- scope (DKIST) site survey. Every month, the data were collected manually for at least one week. We collected statistics of the sky brightness at four bandpasses located at 450, 530, 890, and 940 nm. The results indicate that aerosol scattering is of great importance for the diurnal variation of the sky brightness. For most of the year, the sky brightness remains under 20 millionths per airmass before local Noon. On average, the sky brightness is less than 20 millionths accounts for 40.41 % of the total observing time in a clear day. The best observation time is from 9:00 to 13:00 (Beijing time). Lijiang Observatory is therefore suitable for coronagraphs investigating the structures and dynamics of the corona.

We analyze an asteroseismic binary system: KIC 7107778, a non-eclipsing, unresolved target, with solar-like oscillations in both components. We used \emph{Kepler} short cadence time series spanning nearly two years to obtain the power spectrum. Oscillation mode parameters were determined using Bayesian inference and a nested sampling Monte Carlo algorithm with the \texttt{DIAMONDS} package. The power profiles of the two components fully overlap, indicating their close similarity. We modelled the two stars with \texttt{MESA} and calculated oscillation frequencies with \texttt{GYRE}. Stellar fundamental parameters (mass, radius and age) were estimated by grid modelling with atmospheric parameters and the oscillation frequencies of $l=0,2$ modes as constraints. Most $l=1$ mixed modes were identified with models searched using a bisection method. Stellar parameters for the two sub-giant stars are: $M_A=1.42\pm0.06$ $M_{\astrosun}$, $M_B=1.39\pm0.03$ $M_{\astrosun}$, $R_A=2.93\pm0.05$ $R_{\astrosun}$, $R_B=2.76\pm0.04$ $R_{\astrosun}$, $t_A=3.32\pm0.54$ Gyr and $t_B=3.51\pm0.33$ Gyr. The mass difference of the system is $\sim$1\%. The results confirm their simultaneous birth and evolution, as is expected from binary formation. KIC 7107778 comprises almost identical twins, and is the first asteroseismic sub-giant binary to be detected.

The Xuyi Schmidt Telescope Photometric Survey of the Galactic Anti-center (XSTPS-GAC) is a photometric sky survey that covers nearly 6 000 deg^2 towards Galactic anti-center in g r i bands. Half of its survey field locates on the Galactic Anti-center disk, which makes XSTPS-GAC highly suitable for searching new open clusters in the GAC region. In this paper, we report new open cluster candidates discovered in this survey, as well as properties of these open cluster candidates, such as age, distance and reddening, derived by isochrone fitting in the color-magnitude diagram (CMD). These open cluster candidates are stellar density peaks detected in the star density maps by applying the method from Koposov et al. (2008). Each candidate is inspected on its true color image composed from XSTPS-GAC three band images. Then its CMD is checked, in order to identify whether the central region stars have a clear isochrone-like trend differing from the background stars. The parameters derived from isochrone fitting for these candidates are mainly based on three band photometry of XSTPS-GAC. Meanwhile, when these new candidates are able to be seen clearly on 2MASS, their parameters are also derived based on 2MASS (J-H, J) CMD. Finally, there are 320 known open clusters rediscovered and 24 new open cluster candidates discovered in this work. Further more, the parameters of these new candidates, as well as another 11 known recovered open clusters, are properly determined for the first time.

Obtaining accurate and precise masses and ages for large numbers of giant stars is of great importance for unraveling the assemblage history of the Galaxy. In this paper, we estimate masses and ages of 6940 red giant branch (RGB) stars with asteroseismic parameters deduced from \emph{Kepler} photometry and stellar atmospheric parameters derived from LAMOST spectra. The typical uncertainties of mass is a few per cent, and that of age is $\sim$\,20 per cent. The sample stars reveal two separate sequences in the age -- [$\alpha$/Fe] relation -- a high--$\alpha$ sequence with stars older than $\sim$\,8\,Gyr and a low--$\alpha$ sequence composed of stars with ages ranging from younger than 1\,Gyr to older than 11\,Gyr. We further investigate the feasibility of deducing ages and masses directly from LAMOST spectra with a machine learning method based on kernel based principal component analysis, taking a sub-sample of these RGB stars as a training data set. We demonstrate that ages thus derived achieve an accuracy of $\sim$\,24 per cent. We also explored the feasibility of estimating ages and masses based on the spectroscopically measured carbon and nitrogen abundances. The results are quite satisfactory and significantly improved compared to the previous studies.

Main sequence turn-off (MSTO) stars have advantages as indicators of Galactic evolution since their ages could be robustly estimated from atmospheric parameters. Hundreds of thousands of MSTO stars have been selected from the LAMOST Galactic sur- vey to study the evolution of the Galaxy, and it is vital to derive accurate stellar parameters. In this work, we select 150 MSTO star candidates from the MSTO stars sample of Xiang that have asteroseismic parameters and determine accurate stellar parameters for these stars combing the asteroseismic parameters deduced from the Kepler photometry and atmospheric parameters deduced from the LAMOST spectra.With this sample, we examine the age deter- mination as well as the contamination rate of the MSTO stars sample. A comparison of age between this work and Xiang shows a mean difference of 0.53 Gyr (7%) and a dispersion of 2.71 Gyr (28%). The results show that 79 of the candidates are MSTO stars, while the others are contaminations from either main sequence or sub-giant stars. The contamination rate for the oldest stars is much higher than that for the younger stars. The main cause for the high contamination rate is found to be the relatively large systematic bias in the LAMOST surface gravity estimates.

Using high temporal and high spatial resolution observations taken by the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory, we present the detailed observational analysis of a high quality quasi-periodic fast- propagating (QFP) magnetosonic wave that was associated with the eruption of a magnetic flux rope and a GOES C5.0 flare. For the first time, we find that the QFP wave lasted during the entire flare lifetime rather than only the rising phase of the accompanying flare as reported in previous studies. In addition, the propagation of the different parts of the wave train showed different kinematics and morphologies. For the southern (northern) part, the speed, duration, intensity variation are about 875 +/- 29 (1485 +/- 233) km/s, 45 (60) minutes, and 4% (2%), and the pronounced periods of them are 106 +/- 12 and 160 +/- 18 (75 +/- 10 and 120 +/- 16) seconds, respectively. It is interesting that the northern part of the wave train showed obvious refraction effect when they pass through a region of strong magnetic field. Periodicity analysis result indicates that all the periods of the QFP wave can be found in the period spectrum of the accompanying flare, suggesting their common physical origin. We propose that the quasi-periodic nonlinear magnetohydrodynamics process in the magnetic reconnection that produces the accompanying flare should be important for exciting of QFP wave, and the different magnetic distribution along different paths can account for the different speeds and morphology evolution of the wave fronts.

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.

We present observations of a small-scale Extreme-ultraviolet (EUV) wave that was associated with a mini-filament eruption and a GOES B1.9 micro-fare in the quiet Sun region. The initiation of the event was due to the photospheric magnetic emergence and cancellation in the eruption source region, which first caused the ejection of a small plasma ejecta, then the ejecta impacted on a nearby mini-filament and thereby led to the filament's eruption and the associated fare. During the filament eruption, an EUV wave at a speed of 182 { 317 km/s was formed ahead of an expanding coronal loop, which propagated faster than the expanding loop and showed obvious deceleration and refection during the propagation. In addition, the EUV wave further resulted in the transverse oscillation of a remote filament whose period and damping time are 15 and 60 minutes, respectively. Based on the observational results, we propose that the small-scale EUV wave should be a fast-mode magnetosonic wave that was driven by the the expanding coronal loop. Moreover, with the application of filament seismology, it is estimated that the radial magnetic field strength is about 7 Gauss. The observations also suggest that small-scale EUV waves associated with miniature solar eruptions share similar driving mechanism and observational characteristics with their large-scale counterparts.

We present the observations of a blowout jet that experienced two distinct ejection stages. The first stage started from the emergence of a small positive magnetic polarity, which cancelled with the nearby negative magnetic field and caused the rising of a mini-filament and its confining loops. This further resulted in a small jet due to the magnetic reconnection between the rising confining loops and the overlying open field. The second ejection stage was mainly due to the successive removal of the confining field by the reconnection. Thus that the filament erupted and the erupting cool filament material directly combined with the hot jet originated form the reconnection region and therefore formed the cool and hot components of the blowout jet. During the two ejection stages, cool H{\alpha} jets are also observed cospatial with their coronal counterparts, but their appearance times are earlier than the hot coronal jets a few minutes. Therefore, the hot coronal jets are possibly caused by the heating of the cool H{\alpha} jets, or the rising of the reconnection height from chromosphere to the corona. The scenario that magnetic reconnection occurred between the confining loops and the overlying open loops are supported by many observational facts, including the bright patches on the both sides of the mini-filament, hot plasma blobs along the jet body, and periodic metric radio type III bursts at the very beginnings of the two stages. The evolution and characteristics of these features manifest the detailed non-linear process in the magnetic reconnection.

We present a three dimensional (3D) extinction analysis in the region toward the supernova remnant (SNR) S147 (G180.0-1.7) using multi-band photometric data from the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anticentre (XSTPS-GAC), 2MASS and WISE. We isolate a previously unrecognised dust structure likely to be associated with SNR S147. The structure, which we term as "S147 dust cloud", is estimated to have a distance $d$ = 1.22 $\pm$ 0.21 kpc, consistent with the conjecture that S147 is associated with pulsar PSR J0538 + 2817. The cloud includes several dense clumps of relatively high extinction that locate on the radio shell of S147 and coincide spatially with the CO and gamma-ray emission features. We conclude that the usage of CO measurements to trace the SNR associated MCs is unavoidably limited by the detection threshold, dust depletion, and the difficulty of distance estimates in the outer Galaxy. 3D dust extinction mapping may provide a better way to identify and study SNR-MC interactions.

We report the discovery of two new unbound hypervelocity stars (HVSs) from the LAMOST spectroscopic surveys. They are respectively a B2V type star of ~ 7 M$_{\rm \odot}$ with a Galactic rest-frame radial velocity of 502 km/s at a Galactocentric radius of ~ 21 kpc and a B7V type star of ~ 4 M$_{\rm \odot}$ with a Galactic rest-frame radial velocity of 408 km/s at a Galactocentric radius of ~ 30 kpc. The origins of the two HVSs are not clear given their currently poorly measured proper motions. However, the future data releases of Gaia should provide proper motion measurements accurate enough to solve this problem. The ongoing LAMOST spectroscopic surveys are expected to yield more HVSs to form a statistical sample, providing vital constraint on understanding the nature of HVSs and their ejection mechanisms.

We present estimates of stellar age and mass for 0.93 million Galactic disk main sequence turn-off and sub-giant stars from the LAMOST Galactic Spectroscopic Surveys. The ages and masses are determined by matching with stellar isochrones using Bayesian algorithm, utilizing effective temperature $T_{\rm eff}$, absolute magnitude ${\rm M}_V$, metallicity [Fe/H] and $\alpha$-element to iron abundance ratio [$\alpha$/Fe] deduced from the LAMOST spectra. Extensive examinations suggest the age and mass estimates are robust. The overall sample stars have a median error of 34 per cent for the age estimates, and half of the stars older than 2\,Gyr have age uncertainties of only 20--30 per cent. Median error for the mass estimates of the whole sample stars is $\sim8$ per cent. The huge dataset demonstrates good correlations among stellar age, [Fe/H] ([$\alpha$/H]) and [$\alpha$/Fe]. Particularly, double sequence features are revealed in the both the age--[$\alpha$/Fe] and age--[Fe/H]([$\alpha$/H]) spaces. In the [Fe/H]--[$\alpha$/Fe] space, stars of 8--10\,Gyr exhibit both the thin and thick disk sequences, while younger (older) stars show only the thin (thick) disk sequence, indicating that the thin disk became prominent 8--10\,Gyr ago, while the thick disk formed earlier and almost quenched 8\,Gyr ago. Stellar ages exhibit positive vertical and negative radial gradients across the disk, and the outer disk of $R\gtrsim$\,9\,kpc exhibits a strong flare in stellar age distribution.

We present observational analysis of two successive two-sided loop jets observed by the ground-based New Vacuum Solar Telescope (NVST) and the space-borne Solar Dynamics Observatory ( SDO). The two successive two-sided loop jets manifested similar evolution process and both were associated with the interaction of two small-scale adjacent filamentary threads, magnetic emerging and cancellation processes at the jet's source region. High temporal and high spatial resolution observations reveal that the two adjacent ends of the two filamentary threads are rooted in opposite magnetic polarities within the source region. The two threads approached to each other, and then an obvious brightening patch is observed at the interaction position. Subsequently, a pair of hot plasma ejections are observed heading to opposite directions along the paths of the two filamentary threads, and with a typical speed of two-sided loop jets of the order 150 km/s. Close to the end of the second jet, we report the formation of a bright hot loop structure at the source region, which suggests the formation of new loops during the interaction. Based on the observational results, we propose that the observed two-sided loop jets are caused by the magnetic reconnection between the two adjacent filamentary threads, largely different from the previous scenario that a two-sided loop jet is generated by magnetic reconnection between an emerging bipole and the overlying horizontal magnetic fields.

High precision and long-lasting \emph{Kepler} data enabled us to estimate stellar properties with asteroseismology as an accurate tool. We performed asteroseismic analysis on six solar-like stars observed by the \emph{Kepler} mission: KIC 6064910, KIC 6766513, KIC 7107778, KIC 10079226, KIC 10147635 and KIC 12069127. The extraction of seismic information includes two parts. First, we obtained two global asteroseismic parameters, mean large separation $\Delta\nu$ and frequency of maximum power $\nu_{\rm{max}}$, with autocorrelation function and collapsed autocorrelation function. Second, we extracted individual oscillation modes $\nu_{nl}$ with low-$l$ degree using a least-squares fit. Stellar grid models were built with Yale Rotating Stellar Evolution Code (YREC) to analyze stellar properties. They covered the range of $M=0.8\sim1.8$ $M_{\odot}$ with a step of 0.02 $M_{\odot}$ and $\rm{[Fe/H]}=-0.3\sim0.4$ $\rm{dex}$ with a step of 0.1 $\rm{dex}$. We used a Bayesian approach to estimate stellar fundamental parameters of the six stars, under the constraints of asteroseismic parameters ($\Delta\nu$, $\nu_{\rm{max}}$) and non-asteroseismic parameters ($T_{\rm{eff}}$, $\rm{[Fe/H]}$). We discover that the six targets include five sub-giant stars with $1.2\sim1.5$ $M_{\odot}$ and one main-sequence star with $1.08M_{\odot}$, and with ages in the range of $3\sim5$ $\rm{Gyr}$.

The multi-TeV $\gamma$-rays from the Galactic Center (GC) have a cutoff at tens of TeV, whereas the diffuse emission has no such cutoff, which is regarded as an indication of PeV proton acceleration by the HESS experiment. It is important to understand the inconsistency and study the possibility that PeV cosmic-ray acceleration could account for the apparently contradictory point and diffuse $\gamma$-ray spectra. In this work, we propose that the cosmic rays are accelerated up to $>$PeV in GC. The interaction between cosmic rays and molecular clouds is responsible for the multi-TeV $\gamma$-ray emissions from both the point source and diffuse sources today. Enhanced by the small volume filling factor (VFF) of the clumpy structure, the absorption of the $\gamma$-rays leads to a sharp cutoff spectrum at tens of TeV produced in the GC. Away from galactic center, the VFF grows and the absorption enhancement becomes negligible. As a result, the spectra of $\gamma$-ray emissions for both point source and diffuse sources can be successfully reproduced under such self-consistent picture. In addition, a "surviving-tail" at $\sim$100 TeV is expected from the point source, which can be observed by future projects CTA and LHAASO. Neutrinos are simultaneously produced during proton-proton (PP) collision. With 5-10 years observations, the KM3Net experiment will be able to detect the PeV source according to our calculation.

Recent investigations have shown that the extended main-sequence turnoffs (eMSTOs) are a common feature of intermediate-age star clusters in the Magellanic Clouds. The eMSTOs are also found in the color-magnitude diagram (CMD) of young-age star clusters. The origin of the eMSTOs is still an open question. Moreover, asteroseismology shows that the value of the overshooting parameter $\delta_{\rm ov}$ of the convective core is not fixed for the stars with an approximatelly equal mass. Thus the MSTO of star clusters may be affected by the overshooting of the convective core (OVCC). We calculated the effects of the OVCC with different $\delta_{\rm ov}$ on the MSTO of young- and intermediate-age star clusters. \textbf{If $\delta_{\rm ov}$ varies between stars in a cluster,} the observed eMSTOs of young- and intermediate-age star clusters can be explained well by the effects. The equivalent age spreads of MSTO caused by the OVCC are related to the age of star clusters and are in good agreement with observed results of many clusters. Moreover, the observed eMSTOs of NGC 1856 are reproduced by the coeval populations with different $\delta_{\rm ov}$. The eMSTOs of star clusters may be relevant to the effects of the OVCC. The effects of the OVCC \textbf{are similar to that of rotation in some respects. But the effects cannot result in a significant split of main sequence of young star clusters at $m_{U}\lesssim 21$.} The presence of a rapid rotation can make the split of main sequence of young star clusters more significant.

We present the second release of value-added catalogues of the LAMOST Spectroscopic Survey of the Galactic Anticentre (LSS-GAC DR2). The catalogues present values of radial velocity $V_{\rm r}$, atmospheric parameters --- effective temperature $T_{\rm eff}$, surface gravity log$g$, metallicity [Fe/H], $\alpha$-element to iron (metal) abundance ratio [$\alpha$/Fe] ([$\alpha$/M]), elemental abundances [C/H] and [N/H], and absolute magnitudes ${\rm M}_V$ and ${\rm M}_{K_{\rm s}}$ deduced from 1.8 million spectra of 1.4 million unique stars targeted by the LSS-GAC since September 2011 until June 2014. The catalogues also give values of interstellar reddening, distance and orbital parameters determined with a variety of techniques, as well as proper motions and multi-band photometry from the far-UV to the mid-IR collected from the literature and various surveys. Accuracies of radial velocities reach 5kms$^{-1}$ for late-type stars, and those of distance estimates range between 10 -- 30 per cent, depending on the spectral signal-to-noise ratios. Precisions of [Fe/H], [C/H] and [N/H] estimates reach 0.1dex, and those of [$\alpha$/Fe] and [$\alpha$/M] reach 0.05dex. The large number of stars, the contiguous sky coverage, the simple yet non-trivial target selection function and the robust estimates of stellar radial velocities and atmospheric parameters, distances and elemental abundances, make the catalogues a valuable data set to study the structure and evolution of the Galaxy, especially the solar-neighbourhood and the outer disk.

We intent to use stellar models with and without alpha-enhancement, as well as asteroseismic analysis, to study two alpha-enhanced stars, KIC 7976303 and KIC 8694723. For the alpha-enhanced models, we adopt [alpha/Fe] = 0.2, and 0.4. For the evolved sub-giant KIC 7976303 with mixed-modes, alpha-enhanced models fit the observations better than the models without alpha-enhancement, and point to a star with M = 1.20 $\pm$ 0.02 M_odot, t = 4.88 $\pm$ 0.08 Gyr, R = 2.04 $\pm$ 0.01 R_odot, and L = 5.07 $\pm$ 0.05 L_odot. For the post turnoff star KIC 8694723, we find that the models fit the observations well in all three cases ([alpha/Fe] = 0.0, 0.2, 0.4). According to the observed metal abundances, only alpha- enhanced models are used to estimate the stellar parameters, which are M = 1.13 $\pm$ 0.06 M_odot, t = 5.47 $\pm$ 0.43 Gyr, R = 1.55 $\pm$ 0.13 R_odot, and L = 3.18 $\pm$ 0.26 L_odot. Our alpha-enhanced models indicate KIC 7976303 has a larger mass and younger age than in previous works based on standard models. Furthermore, the differences of estimated mass and age between the three cases are $\sim$ 0.1 M_odot and $\sim$ 0.5 - 1.3 Gyr. These results suggest that we include alpha-enhancement in the modeling of alpha-enhanced stars, such as members of GCs (Globular Clusters), metal-poor stars in the disc and in the halo.

Photometric data from the Xuyi Schmidt Telescope Photometric Survey of the Galactic Anticentre (XSTPS-GAC) and the Sloan Digital Sky Survey (SDSS) are used to derive the global structure parameters of the smooth components of the Milky Way. The data, which cover nearly 11,000 deg$^2$ sky area and the full range of Galactic latitude, allow us to construct a globally representative Galactic model. The number density distribution of Galactic halo stars is fitted with an oblate spheroid that decays by power law. The best-fit yields an axis ratio and a power law index $\kappa=0.65$ and $p=2.79$, respectively. The $r$-band differential star counts of three dwarf samples are then fitted with a Galactic model. The best-fit model yielded by a Markov Chain Monte Carlo analysis has thin and thick disk scale heights and lengths of $H_{1}=$ 322\,pc and $L_{1}=$2343\,pc, $H_{2}=$794\,pc and $L_{2}=$3638\,pc, a local thick-to-thin disk density ratio of $f_2=$11\,per\,cent, and a local density ratio of the oblate halo to the thin disk of $f_h=$0.16\,per\,cent. The measured star count distribution, which is in good agreement with the above model for most of the sky area, shows a number of statistically significant large scale overdensities, including some of the previously known substructures, such as the Virgo overdensity and the so-called "north near structure", and a new feature between 150\degr $< l < $ 240\degr~and $-1$5\degr $< b < $ $-$5\degr, at an estimated distance between 1.0 and 1.5\,kpc. The Galactic North-South asymmetry in the anticentre is even stronger than previously thought.

Recent precise measurements of cosmic ray spectral revealed an anomalous hardening at ~200 GV for nuclei from PAMELA, CREAM, ATIC, AMS02 experiments and at tens of GeV for primary electron derived from AMS02 experiment. Particularly, the latest observation of pbar/p ratio by AMS02 demonstrated a flat distribution, which further validated the spectrum anomalies of secondary particles. All those new phenomena indicated that the conventional propagation model of cosmic rays meet challenge. In this work, the spatial-dependent propagation coefficient D(r,z,\rho) is employed by tracing the source distribution under the physical picture of two-halo model in DRAGON package. Under such scenario, the model calculation will result in a two-component spectral for primary nuclei and electron. Simultaneously, due to the smaller rigidity dependence of D(r,z,\rho) in galactic disk, the ratio of secondary-to-primary will be inevitablly flatter than the calculation in the conventional propagation model. As a result, we can reproduce the spectral hardening of proton, electron and the flat ratio of pbar/p and B/C by only adopting the spatial-dependent propagation coefficient D(r,z,\rho) in galactic disk.

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.

The ratios $r_{01}$ and $r_{10}$ of small to large separations of KIC 2837475 primarily exhibit an increase behavior in the observed frequency range. The calculations indicate that only the models with overshooting parameter $\delta_{\rm ov}$ between approximately 1.2 and 1.6 can reproduce the observed ratios $r_{01}$ and $r_{10}$ of KIC 2837475. The ratios $r_{01}$ and $r_{10}$ of the frequency separations of p-modes with inner turning points that are located in the overshooting region of convective core can exhibit an increase behavior. The frequencies of the modes that can reach the overshooting region decrease with the increase in $\delta_{\rm ov}$. Thus the ratio distributions are more sensitive to $\delta_{\rm ov}$ than to other parameters. The increase behavior of the KIC 2837475 ratios results from a direct effect of the overshooting of convective core. The characteristic of the ratios provides a strict constraint on stellar models. Observational constraints point to a star with $M=1.490\pm0.018$ $M_{\odot}$, $R=1.67\pm0.01$ $R_{\odot}$, age $=2.8\pm0.4$ Gyr, and $1.2\lesssim$ $\delta_{\rm ov}$ $\lesssim1.6$ for KIC 2837475.

Asteroseismology provides a powerful way to constrain stellar parameters. Solar-like oscillations have been observed on subgiant stars with the \emph{Kepler\/} mission. The continuous and high-precision time series enables us to carry out a detailed asteroseismic study for these stars. We carry out data processing of two subgiants of spectral type G: KIC 6442183 and KIC 11137075 observed with the \emph{Kepler} mission, and perform seismic analysis for the two evolved stars. We estimate the values of global asteroseismic parameters: $\Delta\nu=64.9\pm 0.2 $ $\mu$Hz and $\nu_{\rm max}=1225 \pm 17$ $\mu$Hz for KIC 6442183, $\Delta\nu=65.5\pm 0.2 $ $\mu$Hz and $\nu_{\rm max}=1171 \pm 8$ $\mu$Hz for KIC 11137075, respectively. In addition, we extract the individual mode frequencies of the two stars. We compare stellar models and observations, including mode frequencies and mode inertias. The mode inertias of mixed modes, which are sensitive to the stellar interior, are used to constrain stellar models. We define a quantity $d\nu_{\rm m-p}$ that measures the difference between the mixed modes and the expected pure pressure modes, which is related to the inertia ratio of mixed modes to radial modes. Asteroseismic together with spectroscopic constraints provide the estimations of the stellar parameters: $M = 1.04_{-0.04}^{+0.01} M_{\odot}$, $R = 1.66_{-0.02}^{+0.01} R_{\odot}$ and $t=8.65_{-0.06}^{+1.12}$ Gyr for KIC 6442183, and $M = 1.00_{-0.01}^{+0.01} M_{\odot}$, $R = 1.63_{-0.01}^{+0.01} R_{\odot}$ and $t=10.36_{-0.20}^{+0.01}$ Gyr for KIC 11137075. Either mode inertias or $d\nu_{\rm m-p}$ could be used to constrain stellar models.

In this paper we develop two models for the steady states and evolution of two dimensional isothermal self gravitating and rotating incompressible gas which are based on the hydrodynamic equations for stratified fluid. The first model is for the steady states of the gas while the second addresses the time evolution of the gas subject to some constraints. These models reduce the initial five partial differential equations that govern this system to two for the steady state model and to three for the time dependent model. Analytical and numerical solutions of the model equations are used to study the structure of the resulting steady and time dependent states of the fluid with some possible astrophysical applications.

The standard model of cosmic ray propagation has been very successful in explaining all kinds of the Galactic cosmic ray spectra. However, high precision measurement recently revealed the appreciable discrepancy between data and model expectation, from spectrum observations of $\gamma$-rays, $e^+/e^-$ and probably the $B/C$ ratio starting from $\sim$10 GeV energy. In this work, we propose that the fresh cosmic rays, which are supplied by the young accelerators and detained by local magnetic field, can contribute additional secondary particles interacting with local materials. As this early cosmic ray has a hard spectrum, the model calculation results in a two-component $\gamma$-ray spectrum, which agree very well with the observation. Simultaneously, the expected neutrino number from the galactic plane could contribute $\sim60\%$ of IceCube observation neutrino number below a few hundreds of TeV. The same pp-collision process can account for a significant amount of the positron excesses. Under this model, it is expected that the excesses in $\overline p/p$ and $B/C$ ratio will show up when energy is above $\sim$10 GeV. We look forward that the model will be tested in the near future by new observations from AMS02, IceCube, AS$\gamma$, HAWC and future experiments such as LHASSO, HiSCORE and CTA.

The origin of the knee in cosmic ray spectrum remains to be an unsolved fundamental problem. There are various kinds of models which predict different break positions and the compositions of the knee. In this work, we suggest to use diffuse $\gamma$-rays and neutrinos as probes to test these models. Based on several typical types of the composition models, the diffuse $\gamma$-ray and neutrino spectra are calculated, which show distinctive cutoff behaviours at energies from tens of TeV to multi-PeV. The expected flux will be observable by the newly upgraded Tibet-AS$\gamma$+MD (muon detector) experiment as well as more sensitive future projects, such as LHAASO and HiSCORE. By comparing the neutrino spectrum with the recent observations by IceCube experiment, we find that the diffuse neutrinos from interactions between the cosmic rays and the interstellar medium may not be responsible to the majority of the IceCube events. Future measurements of the neutrinos may be able to identify the Galactic diffuse component and further shed light on the problem of the knee of cosmic rays.

The Sky Brightness Monitor (SBM) is an important instrument to measure the brightness level for the sky condition, which is a critical parameter for judging a site for solar coronal observations. In this paper we present an automatic method for the processing of SBM data in large quantity, which can separate the regions of the Sun and the nearby sky as well as recognize the regions of the supporting arms in the field of view. These processes are implemented on the data acquired by more than one SBM instruments during our site survey project in western China. An analysis applying the result from our processes has been done for the assessment of the scattered-light levels by the instrument. Those results are considerably significant for further investigations and studies, notably to derive a series of the other important atmospheric parameters such as extinctions, aerosol content and precipitable water vapor content for candidate sites. Our processes also provide a possible way for full-disk solar telescopes to track the Sun without an extra guiding system.

The aim of this paper is to determinate the fundamental parameters of six exoplanet host (EH) stars and their planets. While techniques for detecting exoplanets yield properties of the planet only as a function of the properties of the host star, hence, we must accurately determine parameters of EH stars at first. For this reason, we constructed a grid of stellar models including diffusion and rotation-induced extra-mixing with given ranges of input parameters (i.e. mass, metallicity, and initial rotation rate). In addition to the commonly used observational constraints such as the effective temperature T_{eff}, luminosity L and metallicity [Fe/H], we added two observational constraints, the lithium abundance log N (Li) and the rotational period P_{rot}. These two additional observed parameters can make further constrains on the model due to their correlations with mass, age and other stellar properties. Hence, our estimations of fundamental parameters for these EH stars and their planets are with higher precision than previous works. Therefore, the combination of rotational period and lithium help us to obtain more accurate parameters for stars, leading to an improvement of the knowledge of the physical state about the EH stars and their planets.

Winking (oscillating) filaments have been observed for many years. However, observations of successive winking filaments in one event have not been reported yet. In this paper, we present the observations of a chain of winking filaments and a subsequent jet that are observed right after the X2.1 are in AR11283. The event also produced an Extreme-ultraviolet (EUV) wave that has two components: upward dome-like wave (850 km/s) and lateral surface wave (554 km/s) which was very weak (or invisible) in imaging observations. By analyzing the temporal and spatial relationships between the oscillating filaments and the EUV waves, we propose that all the winking filaments and the jet were triggered by the weak (or invisible) lateral surface EUV wave. The oscillation of the filaments last for two or three cycles, and their periods, Doppler velocity amplitudes, and damping times are 11 - 22 minutes, 6 - 14 km/s, and 25 - 60 minutes, respectively. We further estimate the radial component magnetic field and the maximum kinetic energy of the filaments, and they are 5 - 10 Gauss and ~ 10^19 J, respectively. The estimated maximum kinetic energy is comparable to the minimum energy of ordinary EUV waves, suggesting that EUV waves can efficiently launch filament oscillations on their path. Based on our analysis results, we conclude that the EUV wave is a good agent for triggering and connecting successive but separated solar activities in the solar atmosphere, and it is also important for producing solar sympathetic eruptions.

The oscillations of the solar-like star HD 49933 have been observed thoroughly by CoRot. Two dozens of frequency shifts, which are closely related with the change in magnetic activity, have been measured. To explore the effects of the magnetic activity on the frequency shifts, we calculate frequency shifts for the radial and $l = 1$ p-modes of HD 49933 with the general variational method, which evaluates the shifts using a spatial integral of the product of a kernel and some sources. The theoretical frequency shifts well reproduce the observation. The magnitudes and positions of the sources are determined according to the $\chi^2$ criterion. We predict the source that contributes to both $l = 0$ and $l = 1$ modes is located at $0.48 - 0.62$Mm below the stellar surface. In addition, based on the assumption that $A_{0}$ is proportional to the change in the MgII activity index $\Delta{i}_{MgII}$, we obtained that the change of MgII index between minimum and maximum of HD 49933 cycle period is about 0.665. The magnitude of the frequency shifts compared to the Sun already told us that HD 49933 is much more active than the Sun, which is further confirmed in this paper. Furthermore, our calculation on the frequency shifts of $l = 1$ modes indicates the variation of turbulent velocity in the stellar convective zone may be an important source for the $l = 1$ shifts.

The frequency ratios and of HD 49933 exhibit an increase at high frequencies. This behavior also exists in the ratios of other stars, which is considered to result from the low signal-to-noise ratio and the larger line width at the high-frequency end and could not be predicted by stellar models in previous work. Our calculations show that the behavior not only can be reproduced by stellar models, but can be predicted by asymptotic formulas of the ratios. The frequency ratios of the Sun, too, can be reproduced well by the asymptotic formulas. The increased behavior derives from the fact that the gradient of mean molecular weight at the bottom of the radiative region hinders the propagation of p-modes, while the hindrance does not exist in the convective core. This behavior should exist in the ratios of stars with a large convective core. The characteristic of the ratios at high frequencies provides a strict constraint on stellar models and aids in determining the size of the convective core and the extent of overshooting. Observational constraints point to a star with $M=1.28\pm0.01 M_{\odot}$, $R=1.458\pm0.005 R_{\odot}$, $t=1.83\pm0.1$ Gyr, $r_{cc}=0.16\pm0.02 R_{\odot}$, $\alpha=1.85\pm0.05$, and $\delta_{ov}=0.6\pm0.2$ for HD 49933.

We present an observational study of a quasi-periodic fast propagating (QFP) magnetosonic wave on 2012, April 23. The multiple wave trains were observed along an active region open loop system which has a divergence geometry. The wave trains were first observed in 171 A observations at a distance of 150 Mm from the footpoint of the guiding loop system and with a speed of 689 km/s, then they appeared in 193 A observations after their interaction with a perpendicular underlaying loop system on the path, in the meantime, the wave speed decelerated to 343 km/s quickly within a short timescale. The sudden deceleration of the wave trains and their appearance in 193 A observations caused by the interaction are interpreted through geometric effect and the density increase of the guiding loop system, respectively. On the other hand, with Wavelet and Fourier analysis methods we find that the wave trains has a common period of 80 s with the associated flare. In addition, a few low frequencies are also identified in the QFP wave. We propose that the generation of the period of 80 s was caused by the periodic releasing of energy busts through some nonlinear processes in magnetic reconnection or the so-called oscillatory reconnection mechanism, while the low frequencies detected in the QFP wave were possibly the manifestations of the leakage of pressure-driven oscillations from the photosphere or chromosphere, which could be an important source for driving QFP waves in the low corona. Our observational results also indicate that the properties of the guiding magnetic structure such as the distributions of magnetic field and density as well as geometry are crucial for modulating the propagation behaviors of QFP waves. Therefore, QFP waves could be used for remote diagnostics of the local physical properties of the solar corona.

We present observations of the diffraction, refraction, and reflection of a global extreme-ultraviolet (EUV) wave propagating in the solar corona. These intriguing phenomena are observed when the wave interacts with two remote active regions, and they together exhibit the wave property of this EUV wave. When the wave approached AR11465, it became weaker and finally disappeared in the active region, but a few minutes latter a new wavefront appeared behind the active region, and it was not concentric with the incoming wave. In addition, a reflected wave was also observed simultaneously on the wave incoming side. When the wave approached AR11459, it transmitted through the active region directly and without reflection. The formation of the new wavefront and the transmission could be explained with diffraction and refraction effects, respectively. We propose that the different behaviors observed during the interactions may caused by different speed gradients at the boundaries of the two active regions. For the origin of the EUV wave, we find that it formed ahead of a group of expanding loops a few minutes after the start of the loops' expansion, which represents the initiation of the associated coronal mass ejection (CME). Based on these results, we conclude that the EUV wave should be a nonlinear magnetosonic wave or shock driven by the associated CME, which propagated faster than the ambient fast-mode speed and gradually slowed down to an ordinary linear wave. Our observations support the hybrid model that includes both fast wave and slow non-wave components.

It was found that the equatorial velocity of A-type stars undergoes an acceleration in the first third of the main sequence (MS) stage, but the velocity decreases as if the stars were not undergoing any redistribution of angular momentum in the external layers in the last stage of the MS phase. Our calculations show that the acceleration and the decrease of the equatorial velocity can be reproduced by the evolution of the differential rotation zero-age MS model with the angular momentum transport caused by hydrodynamic instabilities during the MS stage. The acceleration results from the fact that the angular momentum stored in the interiors of the stars is transported outwards. In the last stage, the core and the radiative envelope are uncoupling, and the rotation of the envelope is a quasi-solid rotation; the uncoupling and the expansion of the envelope lead to that the decrease of the equatorial velocity approximately follows the slope for the change in the equatorial velocity of the model without any redistribution of angular momentum. When the fractional age 0.3 $\lesssim\mathrm{t/t_{MS}}\lesssim$ 0.5, the equatorial velocity remains almost constant for the stars whose central density increases with age in the early stage of the MS phase, while the velocity decreases with age for the stars whose central density decreases with age in the early stage of the MS phase.

We calculated the populations of core-helium-burning (CHeB) stars and found that the secondary red clump (SRC) stars can form an SRC peak in the distributions of the frequency of maximum seismic amplitude ($\nu_{max}$) and mean large-frequency separation ($\Delta\nu$) of CHeB stars when metallicity $Z \geq$ 0.02. The $\nu_{max}$ and $\Delta\nu$ of CHeB stars are dependent not only on He core mass but on H-shell burning. The SRC peak is composed of the CHeB stars with mass roughly between the critical mass M_{Hef} and M_{Hef}+0.2 while He core mass is between about 0.33 and 0.36 M_{sun}. The location of the SRC peak can be affected by the mixing-length parameter $\alpha$, metallicity $Z$, and overshooting parameter $\delta_{ov}$. A decrease in $\alpha$ or increase in $Z$ or $\delta_{ov}$ leads to a movement of the SRC peak towards a lower frequency. However, the change in $Z$ and $\alpha$ only slightly affects the value of M_{Hef} but the variation in $\delta_{ov}$ can significantly affects the value of M_{Hef}. Thus the SRC peak might aid in determining the value of M_{Hef} and calibrating $\delta_{ov}$. In addition, the effects of convective acceleration of SRC stars and the $\nu_{max}$ of `semi-degenerate' stars decreasing with mass result in the appearance of a shoulder between about 40 and 50 $\mu$hz in the \dnu{} distribution. However, the convective acceleration of stars with M < M_{Hef} leads to the deficit in the $\nu_{max}$ distribution between about 9 and 20 $\mu$hz{}. Moreover, the value of the parameter $b$ of the relation between $\nu_{max}$ and $\Delta\nu$ for the populations with M > M_{Hef} is obviously larger than that for the populations with $M <$ \dmhef{}.

We have calculated the ground state electronic structure of He under pressure from 0 to 1500 GPa using both all-electron full-potential and pseudopotential methods based on the density functional theory (DFT). We find that throughout this pressure range, pseudopotentials yield essentially the same energy-volume curve for all of bcc, fcc, and hcp configurations as does the full-potential method, a strong indication that pseudopotential approximation works well for He both as the common element in some giant planets and as detrimental impurities in fusion reactor materials. The hcp lattice is always the most stable structure and bcc the least stable one. Since the energy preference of hcp over fcc and bcc is within 0.01 eV below 100 GPa and about 0.1 eV at 1500 GPa, on the same order of the error bar in local or semi-local density approximations in DFT, phase transitions can only be discussed with more precise description of electron correlation in Quantum Monte Carlo or DFT-based GW methods.

Double or extended main-sequence turn-offs (DMSTOs) and dual red clump (RC) were observed in intermediate-age clusters, such as in NGC 1846 and 419. the DMSTOs are interpreted as that the cluster has two distinct stellar popula- tions with differences in age of about 200-300 Myr but with the same metallicity. The dual RC is interpreted as a result of a prolonged star formation. Using a stellar population-synthesis method, we calculated the evolutions of binary-star stellar population (BSP). We found that binary interactions and merging can reproduce the dual RC in the color-magnitude diagrams of an intermediate-age cluster, whereas in actuality only a single population exists. Moreover, the binary interactions can lead to an extended MSTO rather than DMSTOs. However, the rest of main sequence, subgiant branch and first giant branch are hardly spread by the binary interactions. Part of the observed dual RC and extended MSTO may be the results of binary interactions and merger.