With the increasing numbers of large stellar survey projects, the quality and quantity of excellent tracers to study the Milky Way is rapidly growing, one of which is the classical Cepheids. Classical Cepheids are high precision standard candles with very low typical uncertainties ($<$ 3%) available via the mid-infrared period-luminosity relation. About 3500 classical Cepheids identified from OGLE, ASAS-SN, Gaia, WISE and ZTF survey data have been analyzed in this work, and their spatial distributions show a clear signature of Galactic warp. Two kinematical methods are adopted to measure the Galactic rotation curve in the Galactocentric distance range of $4\lesssim R_{\rm GC} \lesssim 19$ kpc. Gently declining rotation curves are derived by both the proper motion (PM) method and 3-dimensional velocity vector (3DV) method. The largest sample of classical Cepheids with most accurate 6D phase-space coordinates available to date are modeled in the 3DV method, and the resulting rotation curve is found to decline at the relatively smaller gradient of ($-1.33\pm0.1$) ${\rm km\,s^{-1}\,kpc^{-1}}$. Comparing to results from the PM method, a higher rotation velocity (($232.5\pm0.83$) ${\rm km\,s^{-1}}$) is derived at the position of Sun in the 3DV method. The virial mass and local dark matter density are estimated from the 3DV method which is the more reliable method, $M_{\rm vir} = (0.822\pm0.052)\times 10^{12}\,M_\odot$ and $\rho_{\rm DM,\odot} = 0.33\pm0.03$ GeV ${\rm cm^{-3}}$, respectively.

The evolution of an accreting white dwarf (WD) with strong magnetic field toward a type Ia supernova (SN Ia) may differ from the classical single-degenerate (SD) channel. In this paper, we perform binary population synthesis (BPS) simulations for the SD channel with a main-sequence (MS) companion, including the strongly magnetized WD accretion. Under a reasonable assumption that the fraction of such systems is $\sim 15$%, the resulting delay time distribution (DTD) roughly follows the $t^{-1}$ power-law distribution. Within the (WD/MS) SD channel, the contribution from the highly magnetized WD is estimated to be comparable to that from the classical, non-magnetized WD channel. The contribution of the SD channel toward SNe Ia can be at least $\sim 30$% among the whole SN Ia population. We suggest that the SNe Ia resulting from the highly-magnetized WD systems would not share observational properties expected for the classical SD channel; for every (potentially peculiar) SN observationally associated with the SD channel, we expect a comparable number of the `hidden' SD population in the normal class.

Redbacks (RBs) and black widows (BWs) are two peculiar classes of eclipsing millisecond pulsars (MSPs). The accretion-induced collapse (AIC) of an oxygen/neon/magnesium composition white dwarf to a neutron star has been suggested as one possible formation pathway for those two classes of MSPs. However, it is difficult to produce all known MSPs with the traditional AIC scenario. In this study by using the MESA stellar evolution code, we investigate the detailed pre-AIC evolution of magnetized white dwarf binaries with the magnetic confinement model where the high magnetic field strength of the white dwarf can confine the accreted matter in the polar caps. We find that the initial donor mass and orbital periods in our model can be lower than that of previous traditional AIC models. We also present post-AIC evolution models to form RBs and BWs with and without the spin down luminosity evaporation of MSPs. Under the magnetic confinement model and evaporative winds (with corresponding angular momentum loss from the surface of the donor star), the companion masses and orbital periods of all known RBs can be covered and a number of binaries can evolve to become BWs.

With the increasing number of observed magnetic white dwarfs (WDs), the role of magnetic field of the WD in both single and binary evolutions should draw more attentions. In this study, we investigate the WD/main-sequence star binary evolution with the Modules for Experiments in Stellar Astrophysics (MESA code), by considering WDs with non-, intermediate and high magnetic field strength. We mainly focus on how the strong magnetic field of the WD (in a polar-like system) affects the binary evolution towards type Ia supernovae (SNe Ia). The accreted matter goes along the magnetic field lines and falls down onto polar caps, and it can be confined by the strong magnetic field of the WD, so that the enhanced isotropic pole-mass transfer rate can let the WD grow in mass even with a low mass donor with the low Roche-lobe overflow mass transfer rate. The results under the magnetic confinement model show that both initial parameter space for SNe Ia and characteristics of the donors after SNe Ia are quite distinguishable from those found in pervious SNe Ia progenitor models. The predicted natures of the donors are compatible with the non-detection of a companion in several SN remnants and nearby SNe.

We report on the discovery of five carbon-enhanced metal-poor (CEMP) stars in the metallicity range of $-3.3<$ [Fe/H] $<-2.4$. These stars were selected from the LAMOST DR3 low-resolution (R$ \sim$ 2,000) spectroscopic database as metal-poor candidates and followed-up with high-resolution spectroscopy (R$ \sim$110,000) with the LICK/APF. Stellar parameters and individual abundances for 25 chemical elements (from Li to Eu) are presented for the first time. These stars exhibit chemical abundance patterns that are similar to those reported in other literature studies of very and extremely metal-poor stars. One of our targets, J2114$-$0616, shows high enhancement in carbon ([C/Fe]=1.37), nitrogen ([N/Fe]= 1.88), barium ([Ba/Fe]=1.00), and europium ([Eu/Fe]=0.84). Such chemical abundance pattern suggests that J2114$-$0616 can be classified as CEMP-r/s star. In addition, the star J1054+0528 can be classified as a CEMP-rI star, with [Eu/Fe]=0.44 and [Ba/Fe]=$-$0.52. The other stars in our sample show no enhancements in neutron-capture elements and can be classified as CEMP-no stars. We also performed a kinematic and dynamical analysis of the sample stars based on Gaia DR2 data. The kinematic parameters, orbits, and binding energy of these stars, show that J2114$-$0616 is member of the outer halo population, while the remaining stars belong to the inner halo population but with an accreted origin. Collectively, these results add important constraints on the origin and evolution of CEMP stars as well as on their possible formation scenarios.

There are important but unresolved processes in the standard formation scenarios of double compact star binaries (DCBs; BH-BH, BH-NS, NS-NS systems), such as mass transfer and the common envelope (CE) phase. We analyze the effects of different assumptions on key physical processes and binary initial conditions on massive star binary evolution with binary population synthesis (BPS), including a survey of proposed prescriptions for the mass transfer ($\rm q_{\rm cr}$) and the binding energy parameter ($\lambda$) in the CE phase. We find that $\rm q_{\rm cr}$ clearly affects the properties of NS-NS systems while $\lambda$ has influence on the mass distributions of BH-BH systems. The merger rates of DCBs are increased by efficient CE ejection, which in our prescription is related to the binding energy parameter including all the possible budgets to the energy content. It has been suggested that the difference in the properties of GW150914 and GW151226 may reflect different metallicity. We reproduce their properties with our BPS calculations and find that the property of BH-BH systems at low metallicity is sensitive to $\lambda$; the efficient CE ejection leads to a top-heavy mass distribution both for the primary and secondary BHs, which is favored to explain the nature of GW150914. The efficient CE ejection also leads to enhancement of both the BH-BH and NS-NS merger rates to the level consistent with the observational constraints from the detected gravitational wave sources including GW170817.

Most of known RR Lyraes are type ab RR Lyraes (RRLab), and they are the excellent tool to map the Milky Way and its substructures. We find that 1148 RRLab stars determined by Drake et al.(2013) have been observed by spectroscopic surveys of SDSS and LAMOST. We derived radial velocity dispersion, circular velocity and mass profile from 860 halo tracers in our paper I. Here, we present the stellar densities and radial velocity distributions of thick disk and halo of the Milky Way. The 288 RRLab stars located in the thick disk have the mean metallicity of [Fe/H]$=-1.02$. Three thick disk tracers have the radial velocity lower than 215 km $\rm s^{-1}$. With 860 halo tracers which have a mean metallicity of [Fe/H]$=-1.33$, we find a double power-law of $n(r) \propto r^{-2.8}$ and $n(r) \propto r^{-4.8}$ with a break distance of 21 kpc to express the halo stellar density profile. The radial velocity dispersion at 50 kpc is around 78 km $\rm s^{-1}$.

We present a sample of 1148 ab-type RR Lyrae (RRLab) variables identified from Catalina Surveys Data Release 1, combined with SDSS DR8 and LAMOST DR4 spectral data. We firstly use a large sample of 860 Galactic halo RRLab stars and derive the circular velocity distributions for the stellar halo. With the precise distances and carefully determined radial velocities (the center-of-mass radial velocities) by considering the pulsation of the RRLab stars in our sample, we can obtain a reliable and comparable stellar halo circular velocity curve. We take two different prescriptions for the velocity anisotropy parameter {\beta} in the Jeans equation to study the circular velocity curve and mass profile. We test two different solar peculiar motions in our calculation. Our best result with the adopted solar peculiar motion 1 of (U, V, W) = (11.1, 12, 7.2) km/s is that the enclosed mass of the Milky Way within 50 kpc is (3.75 +/- 1.33) *10^11Msun based on \beta = 0 and the circular velocity 180 +/- 31.92 (km/s) at 50 kpc. This result is consistent with dynamical model results, and it is also comparable to the previous similar works.

Binary population synthesis (BPS) study provides a comprehensive way to understand evolutions of binaries and their end products. Close white dwarf (WD) binaries have crucial characteristics in examining in uence of yetunresolved physical parameters on the binary evolution. In this paper, we perform Monte Carlo BPS simulations, investigating the population of WD/main sequence (WD/MS) binaries and double WD binaries, with a publicly available binary star evolution code under 37 different assumptions on key physical processes and binary initial conditions. We considered different combinations of the binding energy parameter (lambda_g:considering gravitational energy only, lambda_b: considering both gravitational energy and internal energy, and lambda_e:considering gravitational energy, internal energy, and entropy of the envelope, the values of them derived with the MESA code), CE effciency, critical mass ratio, initial primary mass function and metallicity. We find that a larger number of post-CE WD/MS binaries in tight orbits are formed when the binding energy parameters is set by lambda_e than the cases adopting the other prescriptions. We also find effects of the other input parameters on orbital period and mass distributions of post-CE WD/MS binaries as well. Containing at least one CO WD, the double WD system evolved from WD/MS binaries may explode as type Ia supernovae (SNe Ia) by merging. In this work, we also investigate a frequency of two WD mergers and compare it to the SNe Ia rate.

CAL 87 is one of the best known supersoft X-ray sources. However, the measured masses, orbital period and orbital period evolution of CAL 87 cannot be addressed by the standard thermal-timescale mass-transfer model for supersoft X-ray sources. In this work we explore the orbital evolution of CAL 87 with both analytic and numerical methods. We demonstrate that the characteristics mentioned above can be naturally accounted for by the excited-wind-driven mass-transfer model.

Accretion-induced collapse of massive white dwarfs (WDs) has been proposed to be an important channel to form binary millisecond pulsars (MSPs). Recent investigations on thermal timescale mass transfer in WD binaries demonstrate that the resultant MSPs are likely to have relatively wide orbit periods ($\gtrsim 10$ days). Here we calculate the evolution of WD binaries taking into account the excited wind from the companion star induced by X-ray irradiation of the accreting WD, which may drive rapid mass transfer even when the companion star is less massive than the WD. This scenario can naturally explain the formation of the strong-field neutron star in the low-mass X-ray binary 4U 1822$-$37. After AIC the mass transfer resumes when the companion star refills its Roche lobe, and the neutron star is recycled due to mass accretion. A large fraction of the binaries will evolve to become binary MSPs with a He WD companion, with the orbital periods distributed between $\gtrsim 0.1$ day and $\lesssim 30$ days, while some of them may follow the cataclysmic variable-like evolution towards very short orbits. If we instead assume that the newborn neutron star appears as an MSP and part of its rotational energy is used to ablate its companion star, the binaries may also evolve to be the redback-like systems.

In the single degenerate scenario for the progenitors of type Ia supernovae (SNe Ia), a white dwarf rapidly accretes hydrogen- or helium-rich material from its companion star, and appears as a supersoft X-ray source. This picture has been challenged by the properties of the supersoft X-ray sources with very low-mass companions and the observations of several nearby SNe Ia. It has been pointed out that, the X-ray radiation or the wind from the accreting white dwarf can excite wind or strip mass from the companion star, thus significantly influence the mass transfer processes. In this paper we perform detailed calculations of the wind-driven evolution of white dwarf binaries. We present the parameter space for the possible SN Ia progenitors, and for the surviving companions after the SNe. The results show that the ex-companion stars of SNe Ia have characteristics more compatible with the observations, compared to those in the traditional single degenerate scenario.