Type Iax supenovae (SNe Iax) are an important sub-luminous class of SNe Ia. However, their progenitors and explosion mechanism have not been understood yet. It has been suggested that SNe Iax may be produced from weak deflagration explosions of Chandrasekhar-mass white dwarfs (WDs) in binary systems with a helium (He) star donor. In such progenitor systems, the stripped He companion material caused by the ejecta-companion interaction is expected to present some He features in their late-time spectra. However, the detection of He lines in late-time spectra of SNe Iax has not yet been successful, which gives an upper limit on the amount of stripped He mass of < 0.002-0.1 solar masses. In this work, we study the interaction between SN Iax ejecta and a He star companion by performing three-dimensional hydrodynamical simulations with a weak pure deflagration explosion model. We find that about 0.004 solar masses of He material can be stripped off from the companion star by SN explosion, which is very close to (or lower than) the observational upper-limit on the total stripped He mass in SNe Iax. We, therefore, conclude that non-detection of He lines in late-time spectra of SNe Iax could be reasonably explained if they are indeed resulted from weak pure deflagration explosions of Chandrasekhar-mass WDs in progenitor systems with a He star donor.

Binary evolution is indispensable in stellar evolution to understand the formation and evolution of most peculiar and energetic objects, such as binary compact objects, Type Ia supernovae, X-ray binaries, cataclysmic variables, blue stragglers, hot subdwarfs, and central binaries in planetary nebulae. Mass transfer in binary stars can change the evolutionary path and fate of the corresponding objects relative to what is expected from single stellar evolution. What is the critical mass ratio at which unstable mass transfer occurs is an unsolved fundamental problem in binary evolution. To resolve this issue, we construct the thermal equilibrium mass loss model and derive critical mass ratios for both thermal timescale mass transfer and unstable mass transfer, the latter of which occurs when the outer Lagrangian point, L2, is overfilled. Using several 3.2 Msun stellar models as examples, we study the stellar response to thermal equilibrium mass loss and present the thresholds for thermal timescale mass transfer. We study the possible mass transfer channels of binary systems containing a 3.2 Msun donor star, taking into account thermal timescale mass transfer, unstable mass transfer through L2, and dynamical timescale mass transfer. We repeat this simulation for a grid of donor stars with different masses (from 0.1 to 100 Msun with Z = 0.02) and at different evolutionary stages, and present our results. The results show that unstable mass transfer due to the overfilling of the outer Lagrangian point may also play an essential role in the formation of common envelopes for late red giant branch and asymptotic giant branch donors.

Since September 2018, LAMOST starts a new 5-year medium-resolution spectroscopic survey (MRS) using bright/gray nights. We present the scientific goals of LAMOST-MRS and propose a near optimistic strategy of the survey. A complete footprint is also provided. Not only the regular medium-resolution survey, but also a time-domain spectroscopic survey is being conducted since 2018 and will be end in 2023. According to the detailed survey plan, we expect that LAMOST-MRS can observe about 2 million stellar spectra with ~7500 and limiting magnitude of around G=15 mag. Moreover, it will also provide about 200 thousand stars with averagely 60-epoch observations and limiting magnitude of G~14 mag. These high quality spectra will give around 20 elemental abundances, rotational velocities, emission line profiles as well as precise radial velocity with uncertainty less than 1 km/s. With these data, we expect that LAMOST can effectively leverage sciences on stellar physics, e.g. exotic binary stars, detailed observation of many types of variable stars etc., planet host stars, emission nebulae, open clusters, young pre-main-sequence stars etc.

Several dozen hypervelocity star (HVS) candidates have been reported based on the second data release of Gaia (Gaia DR2). However, it has been proven that the radial velocities of some Gaia HVS candidates are not reliable. In this paper, we employ refined astrometric criteria to re-examine Gaia DR2, arriving at a more reliable sample of HVS and high velocity star candidates than those found by previous authors.We develop a method called Binary Escape Probability Analysis to identify some HVS candidates. This method allows us to work with stars having only two epochs of measured radial velocity. These stars were usually discarded in previous similar studies. A scrutiny of our final results sheds light on selection effects present in our studies, which we propose to be the focus of future studies. In total, we find three late-type (2 G-type and 1 K-type) HVS and 21 high velocity star candidates, 3 and 11 of which are new, respectively. Judging by their historical trajectories, which we calculate, all three HVS candidates could not have had Galactic center origins. Further monitoring is required to confirm their status.

This paper reports the discovery of an Algol system KIC 10736223 that just past the rapid mass transfer stage. From the light curve and radial-velocity modelling we find KIC 10736223 to be a detached Algol system with the less-massive secondary nearly filling its Roche lobe. Based on the short-cadence Kepler data, we analyzed intrinsic oscillations of the pulsator and identified six secured independent $\delta$ Scuti-type pulsation modes ($f_{1}$, $f_3$, $f_{9}$, $f_{19}$, $f_{42}$, and $f_{48}$). We compute two grids of theoretical models to reproduce the $\delta$ Scuti freqiencies, and find fitting results of mass-accreting models meet well with those of single-star evolutionary models. The fundamental parameters of the primary star yielded with asteroseismology are $M$ = $1.57^{+0.05}_{-0.09}$ $M_{\odot}$, $Z$ = 0.009 $\pm$ 0.001, $R$ = $1.484^{+0.016}_{-0.028}$ $R_{\odot}$, $\log g$ = $4.291^{+0.004}_{-0.009}$, $T_{\rm eff}$ = $7748^{+230}_{-378}$ K, $L$ = $7.136^{+1.014}_{-1.519}$ $L_{\odot}$. The asteroseismic parameters match well with the dynamical parameters derived from the binary model. Moreover, our asteroseismic results show that the pulsator is an almost unevolved star with an age between 9.46-11.65 Myr for single-star evolutionary models and 2.67-3.14 Myr for mass-accreting models. Thereofore, KIC 10736223 may be an Algol system that has just undergone the rapid mass-transfer process.

Double Degenerate systems (DDs) are supposed to be significant gravitational wave (GW) sources for future space-based gravitational-wave detectors, e.g., Laser Interferometer Space Antenna (LISA). Recently, one type of DDs with Extremely low-mass WD (ELM WD; $\leq 0.30\; M_\odot$) companions has been largely found in the ELM Survey. They have very short orbital periods and are therefore important sources for LISA detection. Besides, due to the thick envelope of ELM WDs compared with massive WDs (e.g. CO WDs), they are much easier to be found by the combination of electromagnetic (EM) and GW observations. In this paper, we first obtain the population of ELM WDs in DDs with considering the detailed evolutionary tracks of ELM WDs, and then analyse the GW radiation of these systems. We found that about $6\times10^3$ sources could be solely detected by LISA, including $\sim2\times10^3$ chirping sources, and $\sim13$ ($\sim107$) more sources are expected to be detected by both LISA and ELM Survey (Gaia).

Kilonovae are optical flashes produced in the aftermath of neutron star-neutron star mergers (NNMs) or neutron star-black hole mergers (NBMs). In this work, we use the Millennium Simulation, combined with a semi-analytic galaxy formation model--GABE (Galaxy Assembly with Binary Evolution) which adopts binary stellar population synthesis models, to explore the cosmic event rate of kilonovae, and the properties of their host galaxies in a cosmological context. We find that model with supernova kick velocity of 0 km/s fits the observation best, in agreement with the exception of some formation channels of binary neutron star. This indicates that NNMs prefer to originate from binary systems with low kick velocities. With V$_{\rm kick}$=0 km/s, the cosmic event rate of NNMs and NBMs at z=0 are 283 Gpc$^{-3}$yr$^{-1}$ and 91 Gpc$^{-3}$yr$^{-1}$, respectively, marginally consistent with the constraint from LVC GWTC-1. For Milky Way-mass galaxies, we predict the NNM rate is $25.7^{+59.6}_{-7.1}$ Myr$^{-1}$, which is also in good agreement with the observed properties of binary neutron stars in the Milky Way. Taking all the NNMs into account in the history of Milky Way-mass galaxies, we find that the averaged r-process elements yield with A>79 in a NNM and NBM event should be 0.01 M$_{\odot}$ to be consistent with observation. We conclude that NGC 4993, the host galaxy of GW170817, is a typical host galaxy for NNMs. However, generally NNMs and NBMs tend to reside in young, blue, star-forming, late-type galaxies, with stellar mass and gaseous metallicity distribution peaking at $10^{10.65}$ M$_{\odot}$ and 8.72-8.85, respectively. By studying kilonovae host galaxies in the cosmological background, it is promising to constrain model details better when we have more events in the forthcoming future. (abridged)

Context: Subdwarf B stars (sdBs) play a crucial role in stellar evolution, asteroseismology, and far-UV radiation of early-type galaxies, and have been intensively studied with observation and theory. It has theoretically been predicted that sdBs with neutron star (NS) companions exist in the Galaxy, but none have been discovered yet. This remains a puzzle in this field. In a previous study (hereafter Paper I), we have studied the formation channels of sdB+NS binaries from main-sequence (MS) stars plus NS binaries by establishing a model grid, but it is still unclear how these binaries consisting of MS stars and NS binaries came to be in the first place. Aims: We systematically study the formation of sdB+NS binaries from their original zero-age main-sequence progenitors. We bridge the gap left by our previous study in this way. We obtain the statistical population properties of sdB+NS binaries and provide some guidance for observational efforts. Methods: We first used Hurley's rapid binary evolution code BSE to evolve 10^7 primordial binaries to the point where the companions of NS+MS, NS+Hertzsprung gap (HG) star, and NS+Giant Branch (GB) star binaries have just filled their Roche lobes. Next, we injected these binaries into the model grid we developed in Paper I to obtain the properties of the sdB+NS populations. We adopted two prescriptions of NS natal kicks. Different values of common-envelope ejection efficiency were chosen to examine the effect of common-envelope evolution on the results. Conclusions: Most sdB+NS binaries are located in the Galactic disk with small RV semi-amplitudes. SdB+NS binaries with large RV semi-amplitudes are expected to be strong GWR sources, some of which could be detected by LISA in the future.

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

The radial velocity (RV) is a basic physical quantity which can be determined through Doppler shift of the spectrum of a star. The precision of RV measurement depends on the resolution of the spectrum we used and the accuracy of wavelength calibration. In this work, radial velocities of LAMOST-II medium resolution (R ~ 7500) spectra are measured for 1,594,956 spectra (each spectrum has two wavebands) through matching with templates. A set of RV standard stars are used to recalibrate the zero point of the measurement, and some reference sets with RVs derived from medium/high-resolution observations are used to evaluate the accuracy of the measurement. Comparing with reference sets, the accuracy of our measurement can get 0.0227 km s/1 with respect to radial velocities standard stars. The intrinsic precision is estimated with the multiple observations of single stars, which can achieve to 1.36 km s/1,1.08 km s/1, 0.91 km s/1 for the spectra at signal-to-noise levels of 10, 20, 50, respectively.

Novae are the observational manifestations of thermonuclear runaways on the surface of accreting white dwarfs (WDs). Although novae are an ubiquitous phenomenon, their properties at low metallicity are not well understood. Using the publicly-available stellar evolution code Modules for Experiments in Stellar Astrophysics (MESA), we model the evolution of accreting carbon-oxygen WDs and consider models which accrete matter with metallicity Z=0.02 or $10^{-4}$. We consider both models without mixing and with matter enriched by CO-elements assuming that mixing occurs in the process of accretion (with mixing fraction 0.25). We present and contrast ignition mass, ejected mass, recurrence period and maximum luminosity of novae for different WD masses and accretion rates for these metallicities and mixing cases. We find that models with Z = 0.02 have ignition masses and recurrence periods smaller than models with low Z, while the ejected mass and maximum luminosity are larger. Retention efficiency during novae outbursts decreases with increasing metallicity. In our implementation, inclusion of mixing at the H/He interface reduces accreted mass, ejected mass and recurrence period as compared to the no-mixing case, while the maximum luminosity becomes larger. Retention efficiency is significantly reduced, becoming negative in most of our models. For ease of use, we provide a tabular summary of our results.

Subdwarf A-type stars (sdAs) are objects that have hydrogen-rich spectra with surface gravity similar to that of hot subdwarf stars but effective temperature below the zero-age horizontal branch (ZAHB). They are considered to be metal-poor main sequence (MS) stars or extremely low-mass white dwarfs (ELM WDs). In this work, using the stellar evolution code, Modules for Experiments in Stellar Astrophysics (MESA), we investigate the sdAs formed both by the evolution of (pre-)ELM WDs in double degenerate systems (DDs) and metal-poor main sequence with the single evolution models. We find that both of the evolutionary tracks of ELM WDs and metal-poor MSs can explain the observation properties of sdAs. However, the proportion between these two populations are uncertain. In this work, we adopt the method of binary population synthesis of both ELM WDs in the disk and metal-poor MS in the halo to obtain the populations of ELM WDs and metal-poor MSs at different stellar population ages and calculate their proportions. We find that the proportion of metal-poor MSs to sdA population for a stellar population of 10 Gyr is $\sim$98.5 percent, which is consistent with the conclusion that most sdAs ($>$95 percent) are metal-poor MSs. And the proportion of ELM WDs (metal-poor MSs) to sdA population increases (decreases) from 0.1\% (99.9\%) to 20\% (80\%) with the stellar population ages from 5 to 13.7 Gyr.

Combing Gaia DR2 with LAMOST DR5, we spectroscopically identified 924 hot subdwarf stars, among which 32 stars exhibit strong double-lined composite spectra. We measured the effective temperature $T_{\rm eff}$, surface gravity $\log\,g$, helium abundance $y=n{\rm He}/n{\rm H}$, and radial velocities of 892 non-composite spectra hot subdwarf stars by fitting LAMOST observations with Tlusty/Synspec non-LTE synthetic spectra. We outlined four different groups in the $T_{\rm eff}-\log\,g$ diagram with our helium abundance classification scheme and two nearly parallel sequences in the $T_{\rm eff}-\log(y)$ diagram. 3D Galactic space motions and orbits of 747 hot subdwarf stars with $(G_{BP}-G_{RP})_{0}<-0.36$ mag were computed using LAMOST radial velocities and Gaia parallaxes and proper motions. Based on the $U-V$ velocity diagram, $J_{z}-$eccentricity diagram, and Galactic orbits, we derived Galactic population classifications and the fractional distributions of the four hot subdwarf helium groups in the halo, thin disk and thick disk. Comparisons with the predictions of binary population synthesis calculations (Han 2008) suggest that He-rich hot subdwarf stars with $\log(y)\ge0$ are from the double helium white dwarfs merger, He-deficient hot subdwarf stars with $-2.2\le\log(y)<-1$ from the common envelope ejection, and He-deficient hot subdwarf stars with $\log(y)<-2.2$ from the stable Roche lobe overflow channels. The relative number of He-rich hot subdwarf stars with $-1\le\log(y)<0$ and $\log(y)\ge0$ in the halo is more than twice the prediction of Zhang et al.(2017), even more than six times in the thin disk, which implies that the mergers of helium white dwarfs with low mass main sequence stars may not be the main formation channel of He-rich hot subdwarf stars with $-1\le\log(y)<0$, specially in younger environments.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

Extremely low-mass white dwarfs (ELM WDs) are helium WDs with a mass less than $\sim$$0.3\rm\;M_\odot$. Most ELM WDs are found in double degenerates (DDs) in the ELM Survey led by Brown and Kilic. These systems are supposed to be significant gravitational-wave sources in the mHz frequency. In this paper, we firstly analyzed the observational characteristics of ELM WDs and found that there are two distinct groups in the ELM WD mass and orbital period plane, indicating two different formation scenarios of such objects, i.e. a stable Roche lobe overflow channel (RL channel) and common envelope ejection channel (CE channel). We then systematically investigated the formation of ELM WDs in DDs by a combination of detailed binary evolution calculation and binary population synthesis. Our study shows that the majority of ELM WDs with mass less than $0.22\rm\;M_\odot$ are formed from the RL channel. The most common progenitor mass in this way is in the range of $1.15-1.45\rm\;M_\odot$ and the resulting ELM WDs have a peak around $0.18\rm\;M_\odot$ when selection effects are taken into account, consistent with observations. The ELM WDs with a mass larger than $0.22\rm\;M_\odot$ are more likely to be from the CE channel and have a peak of ELM WD mass around $0.25\rm\;M_\odot$ which needs to be confirmed by future observations. By assuming a constant star formation rate of 2$\rm\;M_\odot yr^{-1}$ for a Milky Way-like galaxy, the birth rate and local density are $5\times10^{-4}\rm\;yr^{-1}$ and $1500\rm\;kpc^{-3}$, respectively, for DDs with an ELM WD mass less than $0.25\rm\;M_\odot$.

Fitting the multi-wavelength spectral energy distributions (SEDs) of galaxies is a widely used technique to extract information about the physical properties of galaxies. However, a major difficulty lies in the numerous uncertainties regarding almost all ingredients of the SED modeling of galaxies. The Bayesian methods provide a consistent conceptual basis for dealing with the problem of inference with many uncertainties. While the Bayesian parameter estimation method have become quite popular in the field of SED fitting of galaxies, the Bayesian model comparison method, which is based on the same Bayes' rule, is still not widely used in this field. With the application of Bayesian model comparison method in a series of papers, we show that the results obtained with Bayesian model comparison are understandable in the context of stellar/galaxy physics. These results indicate that Bayesian model comparison is a reliable and very powerful method for the SED fitting of galaxies.

We have investigated the outburst properties of low-mass X-ray binary transients (LMXBTs) based on a comprehensive study of the outbursts observed in the past few decades. The outburst rates were estimated based on the X-ray monitoring data from Swift/BAT, RXTE/ASM and MAXI, and previous reports in the literature. We found that almost all LMXBTs with the orbital period below $\sim$12 hr showed only one outburst in these observations. There are systematic difference in the outburst rate between long-period ($P_{\rm orb} \gtrsim$ 12 hr) and short-period ($P_{\rm orb} \lesssim$ 12 hr) systems. We infer that mass transfer rate is responsible for the systematic difference, since the disk instability model (DIM) suggested that the mass transfer rate is a key factor affecting the quiescence time. The difference in outburst rate between long-period and short-period LMXBTs is probably due to the different mass transfer mechanism at different evolutionary stages of the donors. Based on the evolutionary tracks of single stars, we derived the critical orbital period for X-ray binaries that harbor a subgiant donor in various metallicity. The critical orbital period ($P_{\rm orb,crit}=$12.4 hr) is consistent with the above orbital period boundary obtained from the statistics of outburst rates. Furthermore, we found a negative correlation between the outburst rate and the orbital period in the samples for which the luminosity class of the donor star is III/IV. The best-fitting power-law index for the black hole subsamples is roughly consistent with the theoretical prediction for those systems with a donor star evolved off the main sequence.

We present a timing analysis of the eclipsing post-common envelope binary (PCEB) DE CVn. Based on new CCD photometric observations and the published data, we found that the orbital period in DE CVn has a cyclic period oscillation with an amplitude of $28.08$ s and a period of $11.22$ years plus a rapid period decrease at a rate of $\dot{P}=-3.35\times10^{-11}ss^{-1}$. According to the evolutionary theory, secular period decreases in PCEBs arise from angular momentum losses (AMLs) driven by gravitational radiation (GR) and magnetic braking (MB). However, the observed orbital decay is too fast to be produced by AMLs via GR and MB, indicating that there could be other AML mechanism. We suggest that a circumbinary disk around DE CVn may be responsible for the additional AML. The disk mass was derived as a few$\times$$10^{-4}$-$10^{-3}$$M_{\odot}$ , which is in agreement with that inferred from previous studies in the order of magnitude. The cyclic change is most likely result of the gravitational perturbation by a circumbinary object due to the Applegate's mechanism fails to explain such a large period oscillation. The mass of the potential third body is calculated as $M_{3}\sin{i'}=0.011(\pm0.003)M_{\odot}$. Supposing the circumbinary companion and the eclipsing binary is coplanar, its mass would correspond to a giant planet. This hypothetical giant planet is moving in a circular orbit of radius $\sim5.75(\pm2.02)$ AU around its host star.

When modeling and interpreting the spectral energy distributions (SEDs) of galaxies, the simple stellar population (SSP) model, star formation history (SFH) and dust attenuation law (DAL) are three of the most important components. However, each of them carries significant uncertainties which have seriously limited our ability to reliably recover the physical properties of galaxies from the analysis of their SEDs. In this paper, we present a Bayesian framework to deal with these uncertain components simultaneously. Based on the Bayesian evidence, a quantitative implement of the principle of Occam's razor, the method allows a more objective and quantitative discrimination among the different assumptions about these uncertain components. With a Ks-selected sample of 5467 low-redshift (mostly with $z\lesssim 1$) galaxies in the COSMOS/UltraVISTA field and classified into passively evolving galaxies (PEGs) and star-forming galaxies (SFGs) with UVJ diagram, we present a Bayesian discrimination of a set of 16 SSP models from five research groups (BC03 and CB07, M05, GALEV, Yunnan-II, BPASS V2.0), five forms of SFH (Burst, Constant, Exp-dec, Exp-inc, Delayed-$\tau$), and four kinds of DAL (Calzetti law, MW, LMC, SMC). We show that the results obtained with the method are either obvious or understandable in the context of galaxy physics. We conclude that the Bayesian model comparison method, especially that for a sample of galaxies, is very useful for discriminating the different assumptions in the SED modeling of galaxies. The new version of the BayeSED code, which is used in this work, is publicly available at https://bitbucket.org/hanyk/bayesed/.

Wide binaries with hot subdwarf-B (sdB) primaries and main sequence companions are thought to form only through stable Roche lobe overflow (RLOF) of the sdB progenitor near the tip of the red giant branch (RGB). We present the orbital parameters of eleven new long period composite sdB binaries based on spectroscopic observations obtained with the UVES, FEROS and CHIRON spectrographs. Using all wide sdB binaries with known orbital parameters, 23 systems, the observed period distribution is found to match very well with theoretical predictions. A second result is the strong correlation between the orbital period (P) and the mass ratio (q) in the observed wide sdB binaries. In the P-q plane two distinct groups emerge, with the main group (18 systems) showing a strong correlation of lower mass ratios at longer orbital periods. The second group are systems that are thought to be formed from higher mass progenitors. Based on theoretical models, a correlation between the initial mass ratio at the start of RLOF and core mass of the sdB progenitor is found, which defines a mass-ratio range at which RLOF is stable on the RGB.

Using a hybrid binary population synthesis approach, we modelled the formation and evolution of populations of accreting WDs for differing star formation histories. We found that the delay time distribution of SNe Ia in the single degenerate scenario is inconsistent with observations. Additionally, we found that our predicted X-ray and UV emission of populations of accreting WDs are consistent with the X-ray luminosities of early-type galaxies observed by Chandra and the HeII 4686\AA/H\beta line ratio measured in stacked SDSS spectra of passively evolving galaxies. Moreover, we found that the majority of current novae in elliptical-like galaxies have low-mass WDs, long decay times, long recurrence periods and are relatively faint. In contrast, the majority of current novae in spiral-like galaxies have massive WDs, short decay times, short recurrence periods and are relatively bright. Our predicted distribution of mass-loss timescales in an M31-like galaxy is consistent with observations for Andromeda.

Binary population synthesis predicts the existence of subdwarf B stars (sdBs) with neutron star (NS) or black hole (BH) companions. We systematically investigate the formation of sdB+NS binaries from binary evolution and aim to obtain some clues for a search for such systems. We started from a series of MS+NS systems and determined the parameter spaces for producing sdB+NS binaries from the stable Roche-lobe overflow (RLOF) channel and from the common envelope (CE) ejection channel. Various NS accretion efficiencies and NS masses were examined to investigate the effects they have. We show the characteristics of the produced sdB+NS systems, such as the mass of components, orbital period, the semi-amplitude of the radial velocity (K), and the spin of the NS component. In the stable RLOF channel, the orbital period of sdB+NS binaries produced in this way ranges from several days to more than 1000 days and moves toward the short-period (~ hr) side with increasing initial MS mass. the sdB+NS systems that result from CE ejection have very short orbital periods and then high values of K (up to 800km s^-1). Such systems are born in very young populations (younger than 0.3 Gyr) and are potential gravitational wave sources that might be resolved by the Laser Interferometer Space Antenna (LISA) in the future. Gravitational wave radiation may again bring them into contact on a timescale of only ~Myr. As a consequence, they are rare and hard to discover. The pulsar signal is likely a feature of sdB+NS systems caused by stable RLOF, and some NS components in sdB binaries may be millisecond pulsars.

Mixing in the convective core is quite uncertain in core helium burning stars. In order to explore the overshooting mixing beyond the convective core, we incorporated the $k$-$\omega$ proposed by Li (2012, 2017) into MESA, and investigated the overshooting mixing in evolution of sdB models. We found that the development of the convective core can be divided into three stages. When the radiative temperature gradient $\nabla_{\rm rad}$ monotonically decreases from the stellar center, the overshooting mixing presents a behaviour of exponential decay similar with Herwig (2000), and the overshooting distance is to make $\nabla_{\rm rad} \simeq \nabla_{\rm ad}$ at the boundary of the convective core, in agreement with the prediction of the self-driving mechanism Castellani, Giannone \& Renzini (1971a). When the radiative temperature gradient $\nabla_{\rm rad}$ shows a minimum value in the convective core, the convective core may be divided into two convection zones if the minimum value of $\nabla_{\rm rad}$ is smaller than the adiabatic temperature gradient $\nabla_{\rm ad}$. For the single-zone case, the overshooting mixing still shows an exponential decay behaviour, but the overshooting distance is much smaller than in the initial stage. For the double-zone case, the overshooting mixing is similar to that of the single-case beyond the convective core, while it almost stops on both sides of the above convective shell. Our overshooting mixing scheme is similar to the maximal overshoot scheme of Constantino et al. (2015). In the final stage, continuous injection of fresh helium into the convective core by the overshooting mixing happens, which is similar to the "core breathing pulses" (Sweigart \& Demarque 1973; Castellani et al. 1985).

Stellar systems consisting of multiple stars tend to undergo tidal interactions when the separations between the stars are short. While tidal phenomena have been extensively studied, a certain tidal effect exclusive to hierarchical triples (triples in which one component star has a much wider orbit than the others) has hardly received any attention, mainly due to its complexity and consequent resistance to being modelled. This tidal effect is the tidal perturbation of the tertiary by the inner binary, which in turn depletes orbital energy from the inner binary, causing the inner binary separation to shrink. In this paper, we develop a fully numerical simulation of these "tertiary tides" by modifying established tidal models. We also provide general insight as to how close a hierarchical triple needs to be in order for such an effect to take place, and demonstrate that our simulations can effectively retrieve the orbital evolution for such systems. We conclude that tertiary tides are a significant factor in the evolution of close hierarchical triples, and strongly influence at least ~ 1% of all multiple star systems.

Using the stellar evolution code---Modules for Experiments in Stellar Astrophysics (MESA), we investigate the evolution of massive stars with different rotational velocities and metallicities towards Wolf-Rayet stars. In our simulations, the initial rotating velocities are taken as 0, 250, 500 and 650 km s$^{-1}$, and the metallicities equal to 0.02, 0.014, 0.008, 0.006, 0.004 and 0.002. We show our rapid rotation models in the HR diagram compared with the observations. We find that the rotational mixing is less efficient at high metallicity, and these stars become Wolf-Rayet (WR) stars when the helium in their center is ignited. However, rapid rotating massive stars at low metallicity can easily evolve into WR stars due to the rotation resulted in chemically homogeneous evolution. This can explain the origin of single WR stars in galaxy at low metallicity. In our models, the observed SMC WR stars are consistent with the single-star evolution models. However at higher metallicities our single-star evolution models can only explain the luminous, hydrogen-rich WN stars and O stars (which are classified as WR stars previously).

As an eclipsing polar with 3.39 hrs orbital period, MN Hya was going through state change when we observed it during 2009-2016. 10 new mid-eclipse times, along with others obtained from literature, allow us to give a new ephemeris. The residuals of linear fit show that period decreased during the phase of state change. It means angular momentum was lost during this phase. The X-ray observation indicates the mass accretion rate as about $3.6\times10^{-9}M_{\odot}yr^{-1}$. The period decrease gives that at least 60 percent of mass being transfered from secondary was lost, maybe in form of the spherically symmetric stellar wind. In high state, the data shows the intensity of the flickering reduced when system had higher accretion rate, and that flickering sticks out with primary timescale about 2 minutes, which implies the position of the threading point as about 30 radius of the white dwarf above the surface of it. The trend of light curves of the system in high state follows that of low state for a large fraction of phase interval from phase 0 to phase 0.4 since which the cyclotron feature is visible, and the primary intensity hump of light curves near phase 0.7 when the system is in high state did not appear on the curve when it is in low state. Those facts contradict the predictions of the two-pole model.

Although type Ia supernovae (SNe Ia) are very useful in many astrophysical fields, their exact progenitor nature is still unclear. A basic method to distinguish the different progenitor models is to search the signal from the single degenerate (SD) model, e.g. the signal for the existence of a non-degenerate companion before or after supernova explosion. Observationally, some SNe Ia show such signal, while the others do not. Here, we propose a universal model to explain these observations based on the spin-up/spin-down model, in which a white dwarf (WD) will experience a spin-down phase before supernova explosion, and the spin-down timescale is determined by its initial mass, i.e. the more massive the initial WD, the shorter the spin-down timescale and then the more likely the SN Ia to show the SD signature. Therefore, our model predicts that the SNe Ia from hybrid carbon-oxygen-neon WDs are more likely to show the SD signature observationally, as some peculiar SNe Ia showed.

Context. The companions of the exploding carbon-oxygen white dwarfs (CO WDs) for producing type Ia supernovae (SNe Ia) are still not conclusively confirmed. A red-giant (RG) star has been suggested to be the mass donor of the exploding WD, named as the symbiotic channel. However, previous studies on the this channel gave a relatively low rate of SNe Ia. Aims. We aim to systematically investigate the parameter space, Galactic rates and delay time distributions of SNe Ia from the symbiotic channel by employing a revised mass-transfer prescription. Methods. We adopted an integrated mass-transfer prescription to calculate the mass-transfer process from a RG star onto the WD. In this prescription, the mass-transfer rate varies with the local material states. Results. We evolved a large number of WD+RG systems, and found that the parameter space of WD+RG systems for producing SNe Ia is significantly enlarged. This channel could produce SNe Ia with intermediate and old ages, contributing to at most 5% of all SNe Ia in the Galaxy. Our model increases the SN Ia rate from this channel by a factor of 5. We suggest that the symbiotic systems RS Oph and T CrB are strong candidates for the progenitors of SNe Ia.

It has been suggested that accretion-induced collapse (AIC) is a non-negligible path for the formation of the observed neutron stars (NSs). An ONe white dwarf (WD) that accretes material from a He star may experience AIC process and eventually produce intermediate-mass binary pulsars (IMBPs), named as the ONe WD+He star scenario. Note that previous studies can only account for part of the observed IMBPs with short orbital periods. In this work, we investigate the evolution of about 900 ONe WD+He star binaries to explore the distribution of IMBPs. We found that the ONe WD+He star scenario could form IMBPs including pulsars with 5-340 ms spin periods and 0.75-1.38 Msun WD companions, in which the orbital periods range from 0.04 to 900 d. Compared with the 20 observed IMBPs, this scenario can cover the parameters of 13 sources in the final orbital period-WD mass plane and the Corbet diagram, most of which has short orbital periods. We found that the ONe WD+He star scenario can explain almost all the observed IMBPs with short orbital periods. This work can well match the observed parameters of PSR J1802-2124 (one of the two precisely observed IMBPs), providing a possible evolutional path for its formation. We also speculate that the compact companion of HD 49798 (a hydrogen depleted sdO6 star) may be not a NS based on the present work.

In this work, we present a catalog of 2651 carbon stars from the fourth Data Release (DR4) of the Large Sky Area Multi-Object Fiber Spectroscopy Telescope (LAMOST). Using an efficient machine-learning algorithm, we find out these stars from more than seven million spectra. As a by-product, 17 carbon-enhanced metal-poor (CEMP) turnoff star candidates are also reported in this paper, and they are preliminarily identified by their atmospheric parameters. Except for 176 stars that could not be given spectral types, we classify the other 2475 carbon stars into five subtypes including 864 C-H, 226 C-R, 400 C-J, 266 C-N, and 719 barium stars based on a series of spectral features. Furthermore, we divide the C-J stars into three subtypes of CJ( H), C-J(R), C-J(N), and about 90% of them are cool N-type stars as expected from previous literature. Beside spectroscopic classification, we also match these carbon stars to multiple broadband photometries. Using ultraviolet photometry data, we find that 25 carbon stars have FUV detections and they are likely to be in binary systems with compact white dwarf companions.

About 766 $\delta$ Scuti stars were observed by LAMOST by June 16, 2017. Stellar atmospheric parameters of 525 variables were determined, while spectral types were obtained for all samples. In the paper, those spectroscopic data are catalogued. We detect a group of 131 unusual and cool variable stars (UCVs) that are distinguished from the normal $\delta$ Scuti stars (NDSTs). On the H-R diagram and the $\log g-T$ diagram, the UCVs are far beyond the red edge of pulsational instability trip. Their effective temperatures are lower than 6700\,K with periods in the range from 0.09 to 0.22\,days. NDSTs have metallicity close to that of the Sun as expected, while UCVs are slightly metal poor than NDSTs. The two peaks on the distributions of the period and stellar atmospheric parameters are all caused by the existence of UCVs. When those UCVs are excluded, it is discovered that the effective temperature, the gravitational acceleration and the metallicity all are correlated with the pulsating period for NDSTs and their physical properties and evolutionary states are discussed. The detecting of those UCVs indicates that about 25\% of the known $\delta$ Scuti stars may be misclassified. However, if some of them are confirmed to be pulsating stars, they will be a new-type pulsator and their investigations will shed light on theoretical instability domains and on the theories of interacting between the pulsation and the convection of solar-type stars. Meanwhile, 88 $\delta$ Scuti stars are detected to be the candidates of binary or multiple systems.

About 3196 EA-type binaries (EAs) were observed by LAMOST by June 16, 2017 and their spectral types were derived. Meanwhile stellar atmospheric parameters of 2020 EAs were determined. In the paper, those EAs are catalogued and their physical properties and evolutionary states are investigated. The period distribution of EAs suggests that the period limit of tidal locking for the close binaries is about 6 days. It is found that the metallicity of EAs is higher than that of EWs indicating that EAs are generally younger than EWs and they are the progenitors of EWs. The metallicities of long-period EWs (0.4 < P < 1 days) are the same as those of EAs with the same periods, while their values of Log (g) are usually smaller than those of EAs. These support the evolutionary process that EAs evolve into long-period EWs through the combination of angular momentum loss (AML) via magnetic braking and case A mass transfer. For short-period EWs, their metallicities are lower than those of EAs, while their gravitational accelerations are higher. These reveal that they may be formed from cool short-period EAs through AML via magnetic braking with little mass transfer. For some EWs with high metallicities, they may be contaminated by material from the evolution of unseen neutron stars and black holes or they have third bodies that may help them to form rapidly through a short timescale of pre-contact evolution. The present investigation suggests that the modern EW populations may be formed through the combination of aforementioned mechanisms.

The double-degenerate (DD) model, involving the merging of massive double carbon-oxygen white dwarfs (CO WDs) driven by gravitational wave radiation, is one of the classical pathways for the formation of type Ia supernovae (SNe Ia). Recently, it has been proposed that the WD+He subgiant channel has a significant contribution to the production of massive double WDs, in which the primary WD accumulates mass by accreting He-rich matter from a He subgiant. We evolved about 1800 CO WD+He star systems and obtained a large and dense grid for producing SNe Ia through the DD model. We then performed a series of binary population synthesis simulations for the DD model, in which the WD+He subgiant channel is calculated by interpolations in this grid. According to our standard model, the Galactic birthrate of SNe Ia is about 2.4*10^{-3} yr^{-1} for the WD+He subgiant channel of the DD model; the total birthrate is about 3.7*10^{-3} yr^{-1} for all channels, reproducing that of observations. Previous theoretical models still have deficit with the observed SNe Ia with delay times <1 Gyr and >8 Gyr. After considering the WD+He subgiant channel, we found that the delay time distributions is comparable with the observed results. Additionally, some recent studies proposed that the violent WD mergers are more likely to produce SNe Ia based on the DD model. We estimated that the violent mergers through the DD model may only contribute to about 16% of all SNe Ia.

Two blue-straggler sequences discovered in globular cluster M30 provide a strong constraint on the formation mechanisms of blue stragglers. We study the formation of the blue-straggler binaries through binary evolution, and find that binary evolution can contribute to the blue stragglers in both of the sequences. Whether a blue straggler is located in the blue sequence or red sequence depends on the contribution of the mass donor to the total luminosity of the binary, which is generally observed as a single star in globular clusters. The blue stragglers in the blue sequence have a cool white-dwarf companion, while the majority ($\sim 60\%$) of the objects in the red sequence are the binaries that are still experiencing mass transfer, but there are also some objects that the donors have just finished the mass transfer (the stripped-core stars, $\sim 10\%$) or the blue stragglers (the accretors) have evolved away from the blue sequence ($\sim 30\%$). Meanwhile, W UMa contact binaries found in both sequences may be explained by various mass ratios, that is, W UMa contact binaries in the red sequence have two components with comparable masses (e.g. mass ratio $q\sim$ 0.3-1.0), while those in the blue sequence have low mass ratio (e.g. $q<$ 0.3). However, the fraction of blue sequence in M30 cannot be reproduced by binary population synthesis if we assumed the initial parameters of a binary sample to be the same as that of the field. This possible indicates that dynamical effects on binary systems is very important in globular clusters.

1SWASP\,J162117.36$+$441254.2 was originally classified as an EW-type binary with a period of 0.20785\,days. However, it was detected to have undergone a stellar outburst on June 3, 2016. Although the system was latter classified as a cataclysmic variable (CV) and the event was attributed as a dwarf-nova outburst, the physical reason is still unknown. This binary has been monitored photometrically since April 19, 2016 and many light curves were obtained before, during and after the outburst. Those light and color curves observed before the outburst indicate that the system is a special CV. The white dwarf is not accreting material from the secondary and there are no accretion disks surrounding the white dwarf. By comparing the light curves obtained from April 19 to September 14, 2016, it was found that magnetic activity of the secondary is associated with the outburst. We show strong evidence that the $L_1$ region on the secondary was heavily spotted before and after the outburst and thus quench the mass transfer, while the outburst is produced by a sudden mass accretion of the white dwarf. These results suggest that J162117 is a good astrophysical laboratory to study stellar magnetic activity and its influences on CV mass transfer and mass accretion.

We present the second paper of a series of publications aiming at obtaining a better understanding regarding the nature of type Ia supernovae (SNIa) progenitors by studying a large sample of detached F, G and K main sequence stars in close orbits with white dwarf companions (i.e. WD+FGK binaries). We employ the LAMOST (Large Sky Area Multi-Object Fibre Spectroscopic Telescope) data release 4 spectroscopic data base together with GALEX (Galaxy Evolution Explorer) ultraviolet fluxes to identify 1,549 WD+FGK binary candidates (1,057 of which are new), thus doubling the number of known sources. We measure the radial velocities of 1,453 of these binaries from the available LAMOST spectra and/or from spectra obtained by us at a wide variety of different telescopes around the globe. The analysis of the radial velocity data allows us to identify 24 systems displaying more than 3sigma radial velocity variation that we classify as close binaries. We also discuss the fraction of close binaries among WD+FGK systems, which we find to be ~10 per cent, and demonstrate that high-resolution spectroscopy is required to efficiently identify double-degenerate SNIa progenitor candidates.

Type Ia supernovae (SNe Ia) play a crucial role in studying cosmology and galactic chemical evolution. They are thought to be thermonuclear explosions of carbon-oxygen white dwarfs (CO WDs) when their masses reach the Chandrasekar mass limit in binaries. Previous studies have suggested that He novae may be progenitor candidates of SNe Ia. However, the mass retention efficiencies during He nova outbursts are still uncertain. In this article, we aim to study the mass retention efficiencies of He nova outbursts and to investigate whether SNe Ia can be produced through He nova outbursts. Using the stellar evolution code Modules for Experiments in Stellar Astrophysics, we simulated a series of multicycle He-layer flashes, in which the initial WD masses range from 0.7 to 1.35 Msun with various accretion rates. We obtained the mass retention efficiencies of He nova outbursts for various initial WD masses, which can be used in the binary population synthesis studies. In our simulations, He nova outbursts can increase the mass of the WD to the Chandrasekar mass limit and the explosive carbon burning can be triggered in the center of the WD; this suggests that He nova outbursts can produce SNe Ia. Meanwhile, the mass retention efficiencies in the present work are lower than those of previous studies, which leads to a lower birthrates of SNe Ia through the WD + He star channel. Furthermore, we obtained the elemental abundances distribution at the moment of explosive carbon burning, which can be used as the initial input parameters in studying explosion models of SNe Ia.

He accretion onto carbon-oxygen white dwarfs (CO WDs) plays a fundamental role when studying the formation of type Ia supernovae (SNe Ia). Employing the MESA stellar evolution code, we calculated the long-term evolution of He-accreting CO WDs. Previous studies usually supposed that a WD can grow in mass to the Chandrasekhar limit in the stable He burning region and finally produce a SN Ia. However, in this study we find that off-centre carbon ignition occurs in the stable He burning region if the accretion rate is above a critical value (~2.05*10^{-6}{M$_\odot$}{yr}^-1), resulting in accretion-induced collapse rather than a SN Ia. If the accretion rate is below the critical value, explosive carbon ignition will eventually happen in the centre producing a SN Ia. Taking into account the possibility of off-centre carbon ignition, we have re-determined the initial parameter space that produces SNe Ia in the He star donor channel, one of the promising channels to produce SNe Ia in young populations. Since this parameter space is smaller than was found in the previous study of Wang et al. (2009), the SN Ia rates are also correspondingly smaller. We also determined the chemical abundance profile of the He-accreting WDs at the moment of explosive carbon ignition, which can be used as initial input for SN Ia explosion models.

Context: Although type Ia supernovae (SNe Ia) play a key role in astrophysics, the companions of the exploding carbon-oxygen white dwarfs (CO WDs) are still not completely identified. It has been suggested recently that a He-rich WD (a He WD or a hybrid HeCO WD) merges with a CO WD may produce an SN Ia. This theory was based on the double-detonation model, in which the shock compression in the CO core caused by the surface explosion of the He-rich shell might lead to the explosion of the whole CO WD. However, so far, very few binary population synthesis (BPS) studies have been made on the merger scenario of a CO WD and a He-rich WD in the context of SNe Ia. Aims: We aim to systematically study the Galactic birthrates and delay-time distributions of SNe Ia based on the merger scenario of a CO WD and a He-rich WD. Methods: We performed a series of Monte Carlo BPS simulations to investigate the properties of SNe Ia from the merging of a CO WD and a He-rich WD based on the Hurley rapid binary evolution code. We also considered the influence of different metallicities on the final results. Results: From our simulations, we found that no more than 15% of all SNe Ia stem from the merger scenario of a CO WD and a He-rich WD, and their delay times range from ~110 Myr to the Hubble time. This scenario mainly contributes to SN Ia explosions with intermediate and long delay times. The present work indicates that the merger scenario of a CO WD and a He-rich WD can roughly reproduce the birthrates of SN 1991bg-like events, and cover the range of their delay times. We also found that SN Ia birthrates from this scenario would be higher for the cases with low metallicities.

EX Dra is a long-period eclipsing dwarf nova with $\sim2-3$ mag amplitude outbursts. This star has been monitored photometrically from November, 2009 to March, 2016 and 29 new mid-eclipse times were obtained. By using new data together with the published data, the best fit to the $O-C$ curve indicate that the orbital period of EX Dra have an upward parabolic change while undergoing double-cyclic variations with the periods of 21.4 and 3.99 years, respectively. The upward parabolic change reveals a long-term increase at a rate of $\dot{P}={+7.46}\times10^{-11}{s} {s^{-1}}$. The evolutionary theory of cataclysmic variables (CVs) predicts that, as a CV evolves, the orbital period should be decreasing rather than increasing. Secular increase can be explained as the mass transfer between the secondary and primary or may be just an observed part of a longer cyclic change. Most plausible explanation for the double-cyclic variations is a pair of light travel-time effect via the presence of two companions. Their masses are determined to be $M_{A}sini'_{A}=29.3(\pm0.6) M_{Jup}$ and $M_{B}sini'_{B}=50.8(\pm0.2) M_{Jup}$. When the two companions are coplanar to the orbital plane of the central eclipsing pair, their masses would match to brown dwarfs.

We present the photometric results of the eclipsing cataclysmic variable (CV) WZ Sge near the period minimum ($P_{min}$). Eight new mid-eclipse times were determined and the orbital ephemeris was updated. Our result shows that the orbital period of WZ Sge is decreasing at a rate of $\dot{P}=-2.72(\pm0.23)\times{10^{-13}}\,s s^{-1}$. This secular decrease, coupled with previous detection of its donor, suggest that WZ Sge is a pre-bounce system. Further analysis indicates that the observed period decrease rate is about $1.53$ times higher than pure gravitational radiation (GR) driving. We constructed the evolutionary track of WZ Sge, which predicts that $P_{min}$ of WZ Sge is $\sim77.98 (\pm0.90)$ min. If the orbital period decreases at the current rate, WZ Sge will evolve past its $P_{min}$ after $\sim25.3$ Myr. Based on the period evolution equation we find $\dot{M}_{2}\simeq4.04(\pm0.10)\times10^{-11}M_{\odot}yr^{-1}$, which is compatible with the current concept of CV evolution at ultrashort orbital periods.

DV UMa is an eclipsing dwarf nova with an orbital period of $\sim2.06$ h, which lies just at the bottom edge of the period gap. To detect its orbital period changes we present 12 new mid-eclipse times by using our CCD photometric data and archival data. Combining with the published mid-eclipse times in quiescence, spanning $\sim30$ yr, the latest version of the $O-C$ diagram was obtained and analyzed. The best fit to those available eclipse timings shows that the orbital period of DV UMa is undergoing a cyclic oscillation with a period of $17.58(\pm0.52)$ yr and an amplitude of $71.1(\pm6.7)$ s. The periodic variation most likely arises from the light-travel-time effect via the presence of a circumbinary object because the required energy to drive the Applegate mechanism is too high in this system. The mass of the unseen companion was derived as $M_{3}\sin{i'}=0.025(\pm0.004)M_{\odot}$. If the third body is in the orbital plane (i.e. $i'=i=82.9^{\circ}$) of the eclipsing pair, it would match to a brown dwarf. This hypothetical brown dwarf is orbiting its host star at a separation of $\sim8.6$ AU in an eccentric orbit ($e=0.44$).

New eclipse timings of the Z Cam-type dwarf nova AY Psc were measured and the orbital ephemeris was revised. Based on the long-term AAVSO data, moreover, the outburst behaviors were also explored. Our analysis suggests that the normal outbursts are quasi-periodic, with an amplitude of $\sim2.5(\pm0.1)$ mag and a period of $\sim18.3(\pm0.7)$ days. The amplitude vs. recurrence-time relation of AY Psc is discussed, and we concluded that this relation may represents general properties of dwarf nova (DN) outbursts. The observed standstill ends with an outburst, which is inconsistent with the general picture of Z Cam-type stars. This unusual behavior was considered to be related to the mass-transfer outbursts. Moreover, the average luminosity is brighter during standstills than during outburst cycles. The changes in brightness marks the variations in $\dot{M}_{2}$ due to the disc of AY Psc is nearly steady state. $\dot{M}_{2}$ value was limited to the range from $6.35\times10^{-9}$ to $1.18\times10^{-8}M_{\odot}yr^{-1}$. More detailed examination shows that there are a few small outbursts presence during standstills. These events with amplitudes of $\sim0.5-0.9$ mag are very similar to the stunted outbursts reported in some NLs. We discussed several possible mechanisms and suggested that the most reasonable mechanism for these stunted outbursts is a changing mass-transfer rate.

GSC 4560-02157 is a new eclipsing cataclysmic variable with an orbital period of $0.265359$ days. By using the published $V-$ and $R-$band data together with our observations, we discovered that the $O-C$ curve of GSC 4560-02157 may shows a cyclic variation with the period of $3.51$ years and an amplitude of $1.40$ min. If this variation is caused by a light travel-time effect via the existence of a third body, its mass can be derived as $M_{3}sini'\approx91.08M_{Jup}$, it should be a low-mass star. In addition, several physical parameters were measured. The colour of the secondary star was determined as $V-R=0.77(\pm0.03)$ which corresponds to a spectral type of K2-3. The secondary star's mass was estimated as $M_{2}=0.73(\pm0.02)M_{\odot}$ by combing the derived $V-R$ value around phase 0 with the assumption that it obeys the mass-luminosity relation of the main sequence stars. This mass is consistent with the mass$-$period relation of CV donor stars. For the white dwarf, the eclipse durations and contacts of the white dwarf yield an upper limit of the white dwarf's radius corresponding to a lower limit mass of $M_{1}\approx0.501M_{\odot}$. The overestimated radius and previously published spectral data indicate that the boundary layer may has a very high temperature.

Early polarization observations on Type Ia supernovae (SNe Ia) may reveal the geometry of supernova ejecta, and then put constraints on their explosion mechanism and their progenitor model. We performed a literature search of SNe Ia with polarization measurements and determined the polarization and relative equivalent width (REW) of Si II 635.5-nm absorption feature at -5 days after the maximum light. We found that either the distribution of observed polarization degree is bimodal, i.e. the broad line SNe Ia have systematically higher polarization than all other SNe Ia, or all kind of SNe Ia share the same polarization sequence, i.e. the polarization of Si II 635.5-nm absorption feature increases with the REW. We also discussed the potential meaning of the discovery on the explosion mechanism and progenitor model of SNe Ia.

KIC~8262223 is an eclipsing binary with a short orbital period ($P=1.61$ d). The {\it Kepler} light curves are of Algol-type and display deep and partial eclipses, ellipsoidal variations, and pulsations of Delta Scuti type. We analyzed the {\it Kepler} photometric data, complemented by phase-resolved spectra from the R-C Spectrograph on the 4-meter Mayall telescope at Kitt Peak National Observatory and determined the fundamental parameters of this system. The low mass and oversized secondary ($M_2=0.20M_{\odot}$, $R_2=1.31R_{\odot}$) is the remnant of the donor star that transferred most of its mass to the gainer, and now the primary star. The current primary star is thus not a normal $\delta$ Scuti star but the result of mass accretion from a lower mass progenitor. We discuss the possible evolutionary history and demonstrate with the MESA evolution code that the system can be understood as the result of non-conservative binary evolution similar to that for the formation of EL CVn type binaries. The pulsations of the primary star can be explained as radial and non-radial pressure modes. The equilibrium models from single star evolutionary tracks can match the observed mass and radius ($M_1=1.94M_{\odot}$, $R_1=1.67R_{\odot}$) but the predicted unstable modes associated with these models differ somewhat from those observed. This work presents a preliminary asteroseismic analysis of the abnormal Delta Scuti pulsators, and we discuss the need for better theoretical understanding of such pulsating mass gaining stars.

From the canonical binary scenario, the majority of sdBs are produced from low-mass stars with degenerate cores where helium is ignited in a way of flashes. Due to numerical difficulties, the models of produced sdBs are generally constructed from more massive stars with non-degenerate cores, leaving several uncertainties on the exact characteristics of sdB stars. Employing MESA, we systematically studied the characteristics of sdBs produced from the common envelope (CE) ejection channel, and found that the sdB stars produced from the CE ejection channel appear to form two distinct groups on the effective temperature-gravity diagram. One group (the flash-mixing model) almost has no H-rich envelope and crows at the hottest temperature end of the extremely horizontal branch (EHB), while the other group has significant H-rich envelope and spreads over the whole canonical EHB region. The key factor for the dichotomy of the sdB properties is the development of convection during the first helium flash, which is determined by the interior structure of the star after the CE ejection. For a given initial stellar mass and a given core mass at the onset of the CE, if the CE ejection stops early, the star has a relatively massive H-rich envelope, resulting in a canonical sdB generally. The fact of only a few short-orbital-period sdB binaries being in the flash-mixing sdB region and the lack of He-rich sdBs in short-orbital-period binaries indicate that the flash mixing is not very often in the products of the CE ejection. A falling back process after the CE ejection, similar to that happened in nova, is an appropriate way of increasing the envelope mass, then prevents the flash mixing.

Employing tidally enhanced stellar wind, we studied in binaries the effects of metallicity, mass ratio of primary to secondary, tidal enhancement efficiency and helium abundance on the formation of blue hook (BHk) stars in globular clusters (GCs). A total of 28 sets of binary models combined with different input parameters are studied. For each set of binary model, we presented a range of initial orbital periods that is needed to produce BHk stars in binaries. All the binary models could produce BHk stars within different range of initial orbital periods. We also compared our results with the observation in the Teff-logg diagram of GC NGC 2808 and {\omega} Cen. Most of the BHk stars in these two GCs locate well in the region predicted by our theoretical models, especially when C/N-enhanced model atmospheres are considered. We found that mass ratio of primary to secondary and tidal enhancement efficiency have little effects on the formation of BHk stars in binaries, while metallicity and helium abundance would play important roles, especially for helium abundance. Specifically, with helium abundance increasing in binary models, the space range of initial orbital periods needed to produce BHk stars becomes obviously wider, regardless of other input parameters adopted. Our results were discussed with recent observations and other theoretical models.

Recent observations suggest that some type Ia supernovae (SNe Ia) originate from the merging of two carbon-oxygen white dwarfs (CO WDs). Meanwhile, recent hydrodynamical simulations have indicated that the accretion-induced collapse may be avoided under certain conditions when double WDs merge violently. However, the properties of SNe Ia from this violent merger scenario are highly dependent on a particular mass-accretion stage, the so-called WD + He subgiant channel, during which the primary WD is able to increase its mass by accreting He-rich material from a He subgiant before the systems evolves into a double WD system. In this article, we aim to study this particular evolutionary stage systematically and give the properties of violent WD mergers. By employing the Eggleton stellar evolution code, we followed a large number of binary calculations and obtained the regions in parameter space for producing violent mergers based on the WD + He subgiant channel. According to these simulations, we found that the primary WDs can increase their mass by ~0.10-0.45 Msun during the mass-accretion stage. We then conducted a series of binary population synthesis calculations and found that the Galactic SN Ia birthrate from this channel is about 0.01-0.4*10-3 yr-1. This suggests that the violent WD mergers from this channel may only contribute to ~0.3%-10% of all SNe Ia in our Galaxy. The delay times of violent WD mergers from this channel are >1.7 Gyr, contributing to the SNe Ia in old populations. We also found that the WD + He subgiant channel is the dominent way for producing violent WD mergers that may be able to eventually explode as SNe Ia.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

Theoretical modelling of the evolution of classical and recurrent novae plays an important role in studies of binary evolution, nucleosynthesis and accretion physics. However, from a theoretical perspective the observed statistical properties of novae remain poorly understood. In this paper, we have produced model populations of novae using a hybrid binary population synthesis approach for differing star formation histories (SFHs): a starburst case (elliptical-like galaxies), a constant star formation rate case (spiral-like galaxies) and a composite case (in line with the inferred SFH for M31). We found that the nova rate at 10\;Gyr in an elliptical-like galaxy is $\sim 10-20$ times smaller than a spiral-like galaxy with the same mass. The majority of novae in elliptical-like galaxies at the present epoch are characterized by low mass white dwarfs (WDs), long decay times, relatively faint absolute magnitudes and long recurrence periods. In contrast, the majority of novae in spiral-like galaxies at 10\;Gyr have massive WDs, short decay times, are relatively bright and have short recurrence periods. The mass loss time distribution for novae in our M31-like galaxy is in agreement with observational data for Andromeda. However, it is possible that we underestimate the number of bright novae in our model. This may arise in part due to the present uncertainties in the appropriate bolometric correction for novae.

We carried out a series of binary stellar evolution calculations, in which the effect of tidally enhanced wind on the evolution of WD + RG systems is incorporated. The WDs increase their mass to the Chandrasekhar mass limit, and then explode as SNe Ia. Based on the binary evolution results, we estimated the X-ray/radio (the excess UV) luminosity from the interactions between supernova ejecta and the CSM (the secondary) via some published standard models. We found that the X-ray flux may be high enough to be detected for a nearby SN Ia from a symbiotic system, while the radio flux is more likely to de detected when the companion is an asymptotic giant branch (AGB) star, and for a first giant branch (FGB) companion, the radio flux is generally lower than the detection limit. For two well observed SNe Ia, 2011fe and 2014J, almost all symbiotic systems are excluded by X-ray observations, but WD + FGB systems may not be ruled out by radio observations. The excess UV luminosity resulting from the collision of supernova ejecta with the RG secondary may be high enough to be detected if the secondary fills its Roche lobe at the moment of supernova explosion. The X-ray/radio emissions are more prevalent in SNe Ia from WD + AGB systems, although SNe Ia from such systems are rare. The UV luminosity from the collision of supernova ejecta to RG secondary is high enough, but only 1 in every few hundred SNe Ia manifests the signal from the collision.

We present a catalog of 166 spectroscopically identified hot subdwarf stars from LAMOST DR1, 44 of which show the characteristics of cool companions in their optical spectra. Atmospheric parameters of 122 non-composite spectra subdwarf stars were measured by fitting the profiles of hydrogen (H) and helium (He) lines with synthetic spectra from non-LTE model atmospheres. Most of the sdB stars scatter near the Extreme Horizontal Branch in the $T_{\rm eff}-\log{g}$ diagram and two well defined groups can be outlined. A clustering of He-enriched sdO stars appears near $T_{\rm eff}=45\,000$ K and $\log(g) = 5.8$. The sdB population separates into several nearly parallel sequences in the $T_{\rm eff}-{\rm He}$ abundance diagram with clumps corresponding to those in the $T_{\rm eff}-\log{g}$ diagram. Over $38\,000$ K (sdO) stars show abundance extremes, they are either He-rich or He-deficient and we observe only a few stars in the $ -1 < \log(y) < 0$ abundance range. With increasing temperature these extremes become less prominent and the He abundance approaches to $\log(y)\sim-0.5$. A unique property of our sample is that it covers a large range in apparent magnitudes and galactic latitudes, therefore it contains a mix of stars from different populations and galactic environments. Our results are consistent with the findings of Hirsch (2009) and we conclude that He-rich and He-deficient sdB stars ($\log(y) < 1$) probably origin from different populations. We also find that most sdO and sdB stars lie in a narrow strip in the luminosity and helium abundance plane, which suggests that these atmospheric parameters are correlated.

The correlation between the jet power and accretion disk luminosity is investigated for active galactic nuclei (AGNs) and black hole X-ray binaries (BHXBs) from the literature. The power-law correlation index is steep ($\mu \sim$ 1.0--1.4) for radio loud quasars and the `outliers' track of BHXBs, and it is flatter ($\mu \sim$ 0.3--0.6) for radio loud galaxies and the standard track of BHXBs. The steep-index groups are mostly at higher accretion rates (peaked at Eddington ratio $>$ 0.01) and the flatter-index groups are at relatively low accretion rates (peaked at Eddington ratio $<$ 0.01), implying that the former groups could be dominated by the inner disk accretion of black hole, while the jet in latter groups would be a hybrid production of the accretion and black hole spin. We could still have a fundamental plane of black hole activity for the BHXBs and AGNs with diverse (maybe two kinds of) correlation indices. It is noted that the fundamental plane of black hole activity should be referred to the correlation between the jet power and disk luminosity or equivalently to the correlation between jet power, Eddington ratio and black hole mass, rather than the jet power, disk luminosity and black hole mass.

Supernovae of type Ia (SNe Ia) are believed to be thermonuclear explosions of carbon-oxygen white dwarfs (CO WDs). However, the mass accretion process onto CO WDs is still not completely understood. In this paper, we study the accretion of He-rich matter onto CO WDs and explore a scenario in which a strong wind forms on the surface of the WD if the total luminosity exceeds the Eddington limit. Using a stellar evolution code called modules for experiments in stellar astrophysics (MESA), we simulated the He accretion process onto CO WDs for WDs with masses of 0.6-1.35Msun and various accretion rates of 10^{-8}-10^{-5}Msun/yr. If the contribution of the total luminosity is included when determining the Eddington accretion rate, then a super-Eddington wind could be triggered at relatively lower accretion rates than those of previous studies based on steady-state models. The super-Eddington wind can prevent the WDs with high accretion rates from evolving into red-giant-like He stars. We found that the contributions from thermal energy of the WD are non-negligible, judging by our simulations, even though the nuclear burning energy is the dominating source of luminosity. We also provide the limits of the steady He-burning regime in which the WDs do not lose any accreted matter and increase their mass steadily, and calculated the mass retention efficiency during He layer flashes for various WD masses and accretion rates. These obtained results can be used in future binary population synthesis computations.

Type Ia supernovae luminosities can be corrected to render them useful as standard candles able to probe the expansion history of the universe. This technique was successful applied to discover the present acceleration of the universe. As the number of SNe Ia observed at high redshift increases and analysis techniques are perfected, people aim to use this technique to probe the equation of state of the dark energy. Nevertheless, the nature of SNe Ia progenitors remains controversial and concerns persist about possible evolution effects that may be larger and harder to characterize than the more obvious statistical uncertainties.

Accreting white dwarfs (WDs) with non-degenerate companions are expected to emit in soft X-rays and the UV, if accreted H-rich material burns stably. They are an important component of the unresolved emission of elliptical galaxies, and their combined ionizing luminosity may significantly influence the optical line emission from warm ISM. In an earlier paper we modeled populations of accreting WDs, first generating WD with main-sequence, Hertzsprung gap and red giant companions with the population synthesis code \textsc{BSE}, and then following their evolution with a grid of evolutionary tracks computed with \textsc{MESA}. Now we use these results to estimate the soft X-ray (0.3-0.7keV), H- and He II-ionizing luminosities of nuclear burning WDs and the number of super-soft X-ray sources for galaxies with different star formation histories. For the starburst case, these quantities peak at $\sim 1$ Gyr and decline by $\sim 1-3$ orders of magnitude by the age of 10 Gyr. For stellar ages of $\sim$~10 Gyr, predictions of our model are consistent with soft X-ray luminosities observed by Chandra in nearby elliptical galaxies and He II 4686$\AA/\rm{H}{\beta}$ line ratio measured in stacked SDSS spectra of retired galaxies, the latter characterising the strength and hardness of the UV radiation field. However, the soft X-ray luminosity and He~II~4686$\AA/\rm{H}{\beta}$ ratio are significantly overpredicted for stellar ages of $\lesssim 4-8$ Gyr. We discuss various possibilities to resolve this discrepancy and tentatively conclude that it may be resolved by a modification of the typically used criteria of dynamically unstable mass loss for giant stars.

In the limit of extremely rapid mass transfer, the response of a donor star in an interacting binary becomes asymptotically one of adiabatic expansion. We survey here adiabatic mass loss from Population I stars of mass 0.10 Msun to 100 Msun from the zero age main sequence to the base of the giant branch, or to central hydrogen exhaustion for lower main sequence stars. For intermediate- and high-mass stars, dynamical mass transfer is preceded by an extended phase of thermal time scale mass transfer as the star is stripped of most of its envelope mass. The critical mass ratio qad above which this delayed dynamical instability occurs increases with advancing evolutionary age of the donor star, by ever-increasing factors for more massive donors. Most intermediate- or high-mass binaries with nondegenerate accretors probably evolve into contact before manifesting this instability. As they approach the base of the giant branch, however, and begin developing a convective envelope, qad plummets dramatically among intermediate-mass stars, to values of order unity, and a prompt dynamical instability occurs. Among low-mass stars, the prompt instability prevails throughout main sequence evolution, with q_ad declining with decreasing mass. Our calculated qad agree well with the behavior of timedependent models by Chen & Han (2003) of intermediate-mass stars initiating mass transfer in the Hertzsprung gap. Application of our results to cataclysmic variables, as systems which must be stable against rapid mass transfer, nicely circumscribes the range in qad as a function of orbital period in which they are found. These results are intended to advance the verisimilitude of population synthesis models of close binary evolution.

The mass function of hydrogen-rich atmosphere white dwarfs has been frequently found to reveal a distinctive high-mass excess near 1Msun. However, a significant excess of massive white dwarfs has not been detected in the mass function of the largest white dwarf catalogue to date from the Sloan Digital Sky Survey. Hence, whether a high-mass excess exists or not has remained an open question. In this work we build the mass function of the latest catalogue of data release 10 SDSS hydrogen-rich white dwarfs, including the cool and faint population (i.e. effective temperatures 6,000 <~ Teff <~ 12,000 K, equivalent to 12 mag <~ Mbol <~ 13 mag). We show that the high-mass excess is clearly present in our mass function, and that it disappears only if the hottest (brightest) white dwarfs (those with Teff >~ 12,000 K, Mbol <~ 12 mag) are considered. This naturally explains why previous SDSS mass functions failed at detecting a significant excess of high-mass white dwarfs. Thus, our results provide additional and robust observational evidence for the existence of a distinctive high-mass excess near 1Msun. We investigate possible origins of this feature and argue that the most plausible scenario that may lead to an observed excess of massive white dwarfs is the merger of the degenerate core of a giant star with a main sequence or a white dwarf companion during or shortly after a common envelope event.

Hyper-velocity stars are believed to be ejected out from the Galactic center through dynamical interactions between (binary) stars and the central massive black hole(s). In this paper, we report 19 low mass F/G/K type hyper-velocity star candidates from over one mil- lion stars of the first data release of the LAMOST general survey. We determine the unbound probability for each candidate using a Monte-Carlo simulation by assuming a non-Gaussian proper-motion error distribution, Gaussian heliocentric distance and radial velocity error dis- tributions. The simulation results show that all the candidates have unbound possibilities over 50% as expected, and one of them may even exceed escape velocity with over 90% probabili- ty. In addition, we compare the metallicities of our candidates with the metallicity distribution functions of the Galactic bulge, disk, halo and globular cluster, and conclude that the Galactic bulge or disk is likely the birth place for our candidates.

Symbiotic stars are interacting binary systems with the longest orbital periods. They are typically formed by a white dwarf, a red giant and a nebula. These objects are natural astrophysical laboratories for studying the evolution of binaries. Current estimates of the population of Milky Way symbiotic stars vary from 3000 up to 400000. However, the current census is less than 300. The Large sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST) survey can obtain hundreds of thousands of stellar spectra per year, providing a good opportunity to search for new symbiotic stars. In this work we detect 4 of such binaries among 4,147,802 spectra released by the LAMOST, of which two are new identifications. The first is LAMOST J12280490-014825.7, considered to be an S-type halo symbiotic star. The second is LAMOST J202629.80+423652.0, a D-type symbiotic star.

The Large sky Area Multi-Object Spectroscopic Telescope (LAMOST) General Survey is a spectroscopic survey that will eventually cover approximately half of the celestial sphere and collect 10 million spectra of stars, galaxies and QSOs. Objects both in the pilot survey and the first year general survey are included in the LAMOST First Data Release (DR1). The pilot survey started in October 2011 and ended in June 2012, and the data have been released to the public as the LAMOST Pilot Data Release in August 2012. The general survey started in September 2012, and completed its first year of operation in June 2013. The LAMOST DR1 includes a total of 1202 plates containing 2,955,336 spectra, of which 1,790,879 spectra have observed signal-to-noise S/N >10. All data with S/N>2 are formally released as LAMOST DR1 under the LAMOST data policy. This data release contains a total of 2,204,696 spectra, of which 1,944,329 are stellar spectra, 12,082 are galaxy spectra and 5,017 are quasars. The DR1 includes not only spectra, but also three stellar catalogues with measured parameters: AFGK-type stars with high quality spectra (1,061,918 entries), A-type stars (100,073 entries), and M stars (121,522 entries). This paper introduces the survey design, the observational and instrumental limitations, data reduction and analysis, and some caveats. Description of the FITS structure of spectral files and parameter catalogues is also provided.

The formation mechanism for blue hook (BHk) stars in globular clusters (GCs) is still unclear. Following one of the possible scenario, named late hot flash scenario, we proposed that tidally enhanced stellar wind in binary evolution may provide the huge mass loss on the red giant branch (RGB) and produce BHk stars. Employing the detailed stellar evolution code, Modules for Experiments in Stellar Astrophysics (MESA), we investigated the contributions of tidally enhanced stellar wind as a possible formation channel for BHk stars in GCs. We evolved the primary stars with different initial orbital periods using the binary module in MESA (version 6208) from zero age main-sequence (ZAMS) to post horizontal branch (HB) stage, and obtained their evolution parameters which are compared with the observation. The results are consistent with observation in the color-magnitude diagram (CMD) and the logg-Teff plane for NGC 2808, which is an example GC hosting BHk stars. However, the helium abundance in the surface for our models is higher than the one obtained in BHk stars. This discrepancy between our models and observation is possibly due to the fact that gravitational settling and radiative levitation which are common processes in hot HB stars are not considered in the models as well as the fact that the flash mixing efficiency may be overestimated in the calculations. Our results suggested that tidally enhanced stellar wind in binary evolution is able to naturally provide the huge mass loss on the RGB needed for late hot flash scenario and it is a possible and reasonable formation channel for BHk stars in GCs.

Hypervelocity stars (HVS) travel with velocities so high, that they exceed the escape velocity of the Galaxy. Several acceleration mechanisms have been discussed. Only one HVS (US 708, HVS 2) is a compact helium star. Here we present a spectroscopic and kinematic analysis of US\,708. Travelling with a velocity of $\sim1200\,{\rm km\,s^{-1}}$, it is the fastest unbound star in our Galaxy. In reconstructing its trajectory, the Galactic center becomes very unlikely as an origin, which is hardly consistent with the most favored ejection mechanism for the other HVS. Furthermore, we discovered US\,708 to be a fast rotator. According to our binary evolution model it was spun-up by tidal interaction in a close binary and is likely to be the ejected donor remnant of a thermonuclear supernova.

The super-Eddington wind scenario has been proposed as an alternative way for producing type Ia supernovae (SNe Ia). The super-Eddington wind can naturally prevent the carbon--oxygen white dwarfs (CO WDs) with high mass-accretion rates from becoming red-giant-like stars. Furthermore, it works in low-metallicity environments, which may explain SNe Ia observed at high redshifts. In this article, we systematically investigated the most prominent single-degenerate WD+MS channel based on the super-Eddington wind scenario. We combined the Eggleton stellar evolution code with a rapid binary population synthesis (BPS) approach to predict SN Ia birthrates for the WD+MS channel by adopting the super-Eddington wind scenario and detailed mass-accumulation efficiencies of H-shell flashes on the WDs. Our BPS calculations found that the estimated SN Ia birthrates for the WD+MS channel are ~0.009-0.315*10^{-3}{yr}^{-1} if we adopt the Eddington accretion rate as the critical accretion rate, which are much lower than that of the observations (<10% of the observed SN Ia birthrates). This indicates that the WD+MS channel only contributes a small proportion of all SNe Ia. The birthrates in this simulation are lower than previous studies, the main reason of which is that new mass-accumulation efficiencies of H-shell flashes are adopted. We also found that the critical mass-accretion rate has a significant influence on the birthrates of SNe Ia. Meanwhile, the results of our BPS calculations are sensitive to the values of the common-envelope ejection efficiency.

The origin of the 6.67 hr period X-ray source, 1E161348-5055, in the young supernova remnant RCW 103 is puzzling. We propose that it may be the descendant of a Thorne-Zytkow Object (TZO). A TZO may at its formation have a rapidly spinning neutron star as a core, and a slowly rotating envelope. We found that the core could be braked quickly to an extremely long spin period by the coupling between its magnetic field and the envelope, and that the envelope could be disrupted by some powerful bursts or exhausted via stellar wind. If the envelope is disrupted after the core has spun down, the core will become an extremely long-period compact object, with a slow proper motion speed, surrounded by a supernova-remnant-like shell. These features all agree with the observations of 1E161348-5055. TZOs are expected to have produced extraordinary high abundances of lithium and rapid proton process elements that would remain in the remnants and could be used to test this scenario.

Ca-rich transients (CRTs, 2005E-like supernovae) exhibit unusually strong Ca features during their nebular phase, and their distribution in their host galaxies indicates that they belong to a metal-poor old population. We carried out a series of binary population synthesis and presented four different channels of forming CO WD + He WD pairs. We selected the systems fulfilling the constraints of the old population and the birth rate from all the CO WD + He WD pairs by constraining the component mass of the WD pairs. For the four channels, the stable Roche lobe overflow (RLOF) could have a significant influence on the formation of the WD + WD pairs. Based on their position on the $M_{\rm CO}$-$M_{\rm He}$ plane, the mass-transfer between the components for most of the CO WD + He WD pairs is neither always unstable nor always stable. We found that it is necessary that the CO WDs are less massive than 0.6 $M_{\odot}$ and the He WDs are less massive than 0.25 $M_{\odot}$ if CO WD + He WD pairs fulfill the constraints from both the old population and the birth rate of CRTs. However, the He WD mass is lower than, while the total mass of the low-mass WD pairs is larger than{\bf ,} the ejecta mass of the CRTs derived from observations. Our results imply that the CO WDs participate in CRT explosions and at the same time, a bound remnant could be left after the CRT explosion if the low-mass WD pairs are the progenitors of CRTs. Therefore, whether or not the helium detonation on a low-mass CO WD may lead to the second detonation in the center of the CO WDs, and whether or not the binding remnant is left after the thermonuclear explosion should be examined carefully in the future.

Radiation fields emitted by O, B-type stars or young stellar populations (SPs) are generally considered as significant central ionizing sources (CISs) of classic H II regions. In our previous studies, we show that the inclusion of binary interactions in stellar population synthesis models can significantly increase the ultraviolet spectrum hardness and the number of ionizing photons of intermediate-age (IA, 7<~log(t/yr)<~8) SPs. In this work, we present photoionization models of H II regions ionized by radiation fields emitted by IA SPs, and show that radiation fields of IA SPs, including binary systems, are in theory possible candidates of significant CISs of classic H II regions. When radiation fields of IA SPs comprising binary systems are used as the CISs of classic H II regions, the theoretical strengths of a number of lines (such as [O III]4959', [S II]6716', etc.), weaker than observations, are raised; the border /selection-criterion lines, between star-forming galaxies and AGNs in the diagnostic diagrams (for example, [N II]6583/Ha~[O III]5007/Hb), move into the region occupied originally by AGNs; and He II 1640 line, observed in Lyman break and high-redshift gravitationally lensed galaxies, also can be produced.

From Liu and Han (2014), the accretion-dominated jet power has a linear proportionality with the accretion rate, whereas the power law index is <=0.5 at lower accretion rate. Attributing the jet power in low accretion rate AGN to the black hole spin, it implies that the jet power has a flatter spectrum than the accretion-dominated jet versus the accretion rate. The black hole must be spinning rapidly for producing such jet power efficiently, and this may allow us to find high-spin black holes in the radio-loud low-luminosity AGN.

The hybrid CONe white dwarfs (WDs) have been suggested to be possible progenitors of type Ia supernovae (SNe Ia). In this article, we systematically studied the hybrid CONe WD + He star scenario for the progenitors of SNe Ia, in which a hybrid CONe WD increases its mass to the Chandrasekhar mass limit by accreting He-rich material from a non-degenerate He star. According to a series of detailed binary population synthesis simulations, we obtained the SN Ia birthrates and delay times for this scenario. The SN Ia birthrates for this scenario are ~0.033-0.539*10^(-3)yr^(-1), which roughly accounts for 1-18% of all SNe Ia. The estimated delay times are ~28Myr-178Myr, which are the youngest SNe Ia predicted by any progenitor model so far. We suggest that SNe Ia from this scenario may provide an alternative explanation of type Iax SNe. We also presented some properties of the donors at the point when the WDs reach the Chandrasekhar mass. These properties may be a good starting point for investigating the surviving companions of SNe Ia, and for constraining the progenitor scenario studied in this work.

Accreting, nuclear-burning white dwarfs have been deemed to be candidate progenitors of type Ia supernovae, and to account for supersoft X-ray sources, novae, etc. depending on their accretion rates. We have carried out a binary population synthesis study of their populations using two algorithms. In the first, we use the binary population synthesis code \textsf{BSE} as a baseline for the "rapid" approach commonly used in such studies. In the second, we employ a "hybrid" approach, in which we use \textsf{BSE} to generate a population of white dwarfs (WD) with non-degenerate companions on the verge of filling their Roche lobes. We then follow their mass transfer phase using the detailed stellar evolution code \textsf{MESA}. We investigate the evolution of the number of rapidly accreting white dwarfs (RAWDs) and stably nuclear-burning white dwarfs (SNBWDs), and estimate the type Ia supernovae (SNe Ia) rate produced by "single-degenerate" systems (SD). We find significant differences between the two algorithms in the predicted numbers of SNBWDs at early times, and also in the delay time distribution (DTD) of SD SNe Ia. Such differences in the treatment of mass transfer may partially account for differences in the SNe Ia rate and DTD found by different groups. Adopting 100\% efficiency for helium burning, the rate of SNe Ia produced by the SD-channel in a Milky-way-like galaxy in our calculations is $2.0\times10^{-4}\rm{yr}^{-1}$, more than an order of magnitude below the observationally inferred value. In agreement with previous studies, our calculated SD DTD is inconsistent with observations.

The potential importance of the angular momentum which is gained by accreting white dwarfs (WDs) has been increasingly recognized in the context of type Ia supernova (SN Ia) single-degenerate model. The expectation that the spin of the WD can delay the explosion should help the single-degenerate model to be consistent with the observed properties of most SNe Ia, in particular by avoiding hydrogen contamination. In this article, we attempt to study the most prominent single-degenerate supersoft (WD + MS) channel when the rotation of accreting WDs is considered. We present a detailed binary population synthesis study to examine the predicted population of SNe Ia for this channel. For our standard model, we find that 77% of these SNe Ia explode with WD masses which are low enough to be supported by solid-body rotation (1.378-1.5Msun); this is a substantially higher proportion than found by previous work. Only 2% have WD explosion masses >2.0Msun; these require the initial WD mass to be larger than 1.0 Msun. We further discuss the possible origin of the diversity of SNe Ia from the pre- and post- accretion properties of the WDs in this population. We also suggest that some SN Ia progenitors with substantial circumstellar hydrogen, including some apparent type IIn SNe, might be related to WDs which required support from differential rotation to avoid explosion, since these can still be accreting from hydrogen-rich donors with a relatively high mass-transfer rate at the time of the SN explosion.

A model for the relation between radio jet power and the product of central black hole (BH) mass and Eddington ratio of AGN is proposed, and the model is examined with data from the literature. We find that radio jet power positively correlates but not linearly with the product of BH mass ($m$ in solar mass) and Eddington ratio ($\lambda$), and the power law indices ($\mu$) are significantly less than unity for relatively low accretion ($\lambda<0.1$) AGN, $P_{j}\propto (\lambda m)^{\mu}$, in the radio galaxies and the Seyfert galaxies. This leads to a negative correlation between radio loudness and $\lambda m$ for the low luminosity AGN, i.e. $R\propto (\lambda m)^{\rho}$ with $\rho=(7/6)\mu-1<0$, which may be attributed to a contribution of BH spin to total jet power assuming that the spin induced jet is gradually suppressed as the accretion rate increases. Whereas, for the high-z quasars which often show the slope $\mu\geq1$, a positive correlation between the radio loudness and disc luminosity is predicted. We discuss that the jet powers of the high-z FRII quasars are likely dominated by the accretion disc rather than by the BH spin.

We present a newly developed version of BayeSED, a general Bayesian approach to the spectral energy distribution (SED) fitting of galaxies. The new BayeSED code has been systematically tested on a mock sample of galaxies. The comparison between estimated and inputted value of the parameters show that BayeSED can recover the physical parameters of galaxies reasonably well. We then applied BayeSED to interpret the SEDs of a large Ks-selected sample of galaxies in the COSMOS/UltraVISTA field with stellar population synthesis models. With the new BayeSED code, a Bayesian model comparison of stellar population synthesis models has been done for the first time. We found that the model by Bruzual & Charlot (2003), statistically speaking, has larger Bayesian evidence than the model by Maraston (2005) for the Ks-selected sample. Besides, while setting the stellar metallicity as a free parameter obviously increases the Bayesian evidence of both models, varying the IMF has a notable effect only on the Maraston (2005) model. Meanwhile, the physical parameters estimated with BayeSED are found to be generally consistent with those obtained with the popular grid-based FAST code, while the former exhibits more natural distributions. Based on the estimated physical parameters of galaxies in the sample, we qualitatively classified the galaxies in the sample into five populations that may represent galaxies at different evolution stages or in different environments. We conclude that BayeSED could be a reliable and powerful tool for investigating the formation and evolution of galaxies from the rich multi-wavelength observations currently available. A binary version of MPI parallelized BayeSED code is publicly available at https://bitbucket.org/hanyk/bayesed.

Orbital decay by gravitational-wave radiation will cause some close-binary white dwarfs (WDs) to merge within a Hubble time. The results from previous hydrodynamical WD-merger simulations have been used to guide calculations of the post-merger evolution of carbon-oxygen + helium (CO+He) WD binaries. Our models include the formation of a hot corona in addition to a Keplerian disk. We introduce a 'destroyeddisk' model to simulate the effect of direct disk ingestion into the expanding envelope. These calculations indicate significant lifetimes in the domain of the rare R Coronae Borealis (RCB) stars, before a fast evolution through the domain of the hotter extreme helium (EHe) stars. Surface chemistries of the resulting giants are in partial agreement with the observed abundances of RCB and EHe stars. The production of 3He, 18O and 19F are discussed. Evolutionary timescales combined with binary white-dwarf merger rates from binary-star population synthesis are consistent with present-day numbers of RCBs and EHes, provided that the majority come from relatively recent (< 2Gyr) star formation. However, most RCBs should be produced by CO-WD + low-mass HeWD mergers, with the He-WD having a mass in the range 0.20-0.35 Msolar. Whilst, previously, a high He-WD mass (> 0.40 Msolar) was required to match the carbon-rich abundances of RCB stars, the 'destroyed-disk' model yields a high-carbon product with He-WD mass > 0.30 Msolar, in better agreement with population synthesis results.

Tidally enhanced stellar wind may affect horizontal branch (HB) morphology in globular clusters (GCs) by enhancing the mass loss of primary star during binary evolution. Lei et al. (2013a, 2013b) studied the effect of this kind of wind on HB morphology in details, and their results indicated that binary is a possible secondparameter (2P) candidate in GCs. Binary fraction is an very important fact in the tidally-enhanced-stellar-wind model. In this paper, we studied the effect of binary fraction on HB morphology by removing the effects of metallicity and age. Five different binary fractions (i.e., 10%, 15%, 20%, 30% and 50%) are adopted in our model calculations. The synthetic HB morphologies with different binary fractions are obtained at different metallicities and ages. We found that, due to the great influence of metallicity and age, the effect of binary fraction on HB morphology may be masked by these two parameters. However, when the effects of metallicity and age are removed, the tendency that HB morphologies become bluer with increasing of binary fractions is clearly presented. Furthermore, we compared our results with the observation by Milone et al. (2012). Our results are consistent well with the observation at metalrich and metal-poor GCs. For the GCs with intermediate metallicity, when the effect of age on HB morphology is removed, a weak tendency that HB morphologies become bluer with increasing of binary fractions is presented in all regions of GCs, which is consistent with our results obtained in this metallicity range.

Observations revealed the presence of multiple stellar populations in globular clusters (GCs) that exhibit wide abundance variations and multiple sequences in Hertzsprung-Russell (HR) diagram. We present a scenario for the formation of multiple stellar populations in GCs. In this scenario, initial GCs are single-generation clusters, and our model predicts that the abundance anomalous stars observed in GCs are the merged stars and the accretor stars produced by binary interactions, which are rapidly rotating stars at the moment of their formation and are more massive than normal single stars in the same evolutionary stage. We find that due to their own evolution, these rapidly rotating stars have different surface abundances, effective temperatures and luminosities from normal single stars in the same evolutionary stage. The stellar population with binaries can reproduce two important observational evidences of multiple stellar populations, the Na-O anticorrelation and the multiple sequences in HR diagram. This suggests that binary interactions may be a possible scenario for the formation of multiple stellar populations in GCs.

The detached-binary channel is an important channel for the formation of contact binaries, according to which a detached binary might evolve into contact by evolutionary expansion of the components, or angular momentum loss through the effect of magnetic braking (MB). We have carried out a detailed binary population synthesis (BPS) study of this channel, and obtained the parameter regions for detached binaries to evolve into contact. Combining the observations from the $\emph{Kepler}$ satellite with our results, we found that the ratio of the birth rate of the progenitors of contact binaries to that of contact binaries is greater than about 1.2. This suggests that for the detached-binary channel, the progenitors can be sufficient to produce observed contact binaries. In this channel, we find that the distribution of orbital period for contact binaries has a peak at about 0.25 days and a tail extending to longer periods, and the formation timescale of contact binaries has a large range from $\sim1$\,Myr to 15\,Gyr. These results show that the detached-binary channel could explain satisfactorily the main observational characteristics of contact binaries, in particular the distribution of orbital period shown by the $\emph{Kepler}$ observations and the existence of very young contact binaries.

The most favored progenitor scenarios for Type Ia supernovae (SNe Ia) involve the single-degenerate (SD) scenario and the double-degenerate scenario. The absence of stripped hydrogen (H) in the nebular spectra of SNe Ia challenges the SD progenitor models. Recently, it was shown that pure deflagration explosion models of Chandrasekhar-mass white dwarfs ignited off-center reproduce the characteristic observational features of 2002cx-like SNe Ia very well. In this work we predict, for the first time, the amount of stripped H for the off-center pure deflagration explosions. We find that their low kinetic energies lead to inefficient H mass stripping (less than 0.01 M_sun), indicating that the stripped H may be hidden in (observed) late-time spectra of SN 2002cx-like SNe Ia.

The accretion of hydrogen-rich material onto carbon-oxygen white dwarfs (CO WDs) is crucial for understanding type Ia supernova (SN Ia) from the single-degenerate model, but this process has not been well understood due to the numerical difficulties in treating H and He flashes during the accretion. For the CO WD masses from 0.5 to $1.378\,{M}_\odot$ and accretion rates in the range from $10^{-8}$ to $10^{-5}\,{M}_\odot\,\mbox{yr}^{-1}$, we simulated the accretion of solar-composition material onto CO WDs using the state-of-the-art stellar evolution code of {\sc MESA}. For comparison with the steady-state models (e.g \citet{nskh07}), we firstly ignored the contribution from nuclear burning to the luminosity when determining the Eddington accretion rate and found that the properties of H burning in our accreting CO WD models are similar to those from the steady-state models, except that the critical accretion rates at which the WDs turn into red giants or H-shell flashes occur on their surfaces are slightly higher than those from the steady-state models. However, the super-Eddington wind is triggered at much lower accretion rates, than previously thought, when the contribution of nuclear burning to the total luminosity is included. This super-Eddington wind naturally prevents the CO WDs with high accretion rates from becoming red giants, thus presenting an alternative to the optically thick wind proposed by \cite{hkn96}. Furthermore, the super-Eddington wind works in low-metallicity environments, which may explain SNe Ia observed at high redshifts.

It has recently been proposed that one sub-class of type Ia supernovae (SNe Ia) is sufficiently both distinct and common to be classified separately from the bulk of SNe Ia, with a suggested class name of "type Iax supernovae" (SNe Iax), after SN 2002cx. However, their progenitors are still uncertain. We study whether the population properties of this class might be understood if the events originate from a subset of sub-Chandrasekhar mass explosions. In this potential progenitor population, a carbon--oxygen white dwarf (CO WD) accumulates a helium layer from a non-degenerate helium star; ignition of that helium layer then leads to ignition of the CO WD. We incorporated detailed binary evolution calculations for the progenitor systems into a binary population synthesis model to obtain rates and delay times for such events. The predicted Galactic event rate of these explosions is ~1.5\times10^{-3}{yr}^{-1} according to our standard model, in good agreement with the measured rates of SNe Iax. In addition, predicted delay times are ~70Myr-800Myr, consistent with the fact that most of SNe Iax have been discovered in late-type galaxies. If the explosions are assumed to be double-detonations -- following current model expectations -- then based on the CO WD masses at explosion we also estimate the distribution of resulting SN brightness (-13 \gtrsim M_{bol} \gtrsim -19mag), which can reproduce the empirical diversity of SNe Iax. We speculate on why binaries with non-degenerate donor stars might lead to SNe Iax if similar systems with degenerate donors do not. We suggest that the high mass of the helium layer necessary for ignition at the lower accretion rates typically delivered from non-degenerate donors might be necessary to produce SN 2002cx-like characteristics, perhaps even by changing the nature of the CO ignition.

Eclipsing binary millisecond pulsars (the so-called black widows and redbacks) can provide important information about accretion history, pulsar irradiation of their companion stars and the evolutionary link between accreting X-ray pulsars and isolated millisecond pulsars. However, the formation of such systems is not well understood, nor the difference in progenitor evolution between the two populations of black widows and redbacks. Whereas both populations have orbital periods between $0.1-1.0\;{\rm days}$ their companion masses differ by an order of magnitude. In this paper, we investigate the formation of these systems via evolution of converging low-mass X-ray binaries by employing the MESA stellar evolution code. Our results confirm that one can explain the formation of most of these eclipsing binary millisecond pulsars using this scenario. More notably, we find that the determining factor for producing either black widows or redbacks is the efficiency of the irradiation process, such that the redbacks absorb a larger fraction of the emitted spin-down energy of the radio pulsar (resulting in more efficient mass loss via evaporation) compared to that of the black widow systems. We argue that geometric effects (beaming) is responsible for the strong bimodality of these two populations. Finally, we conclude that redback systems do not evolve into black widow systems with time.

Evolutionary population synthesis (EPS) models play an important role in many studies on the formation and evolution of galaxies. Most EPS models are still poorly calibrated for certain stellar evolution stages, especially for the treatment of binary stars, which are very different from single stars. We aim to present color-magnitude (C-M) and color-color (C-C) relations for passive model galaxies in the redshift z$\sim0.0-3.0$ and to study the effect of binary interactions on these relations for high-redshift passive galaxies. Assuming exponentially declining star formation rate, we used a set of theoretical galaxy templates obtained from Yunnan EPS models (with and without binary interactions) to present the C-M and C-C relations for passive galaxies via Monte Carlo simulation. In Yunnan EPS models with binary interactions, various processes are included, such as mass transfer, mass accretion, common-envelope evolution, collisions, supernova kicks, tidal evolution, and all angular momentum loss mechanisms. In these models, approximately 50 per cent of the stellar systems are binary systems with orbital periods less than 100\,yr. This fraction is a typical value for the Milky Way. We find that the inclusion of binary interactions in the model galaxies' spectra can dramatically alter the predicted C-M and C-C relations and their evolution with redshift. For $z\sim 0.0$ and $1.0$, the binary interactions have a minor effect on the C-M and C-C relations, but at $z\sim 2.0$ and $3.0$ the binary interactions have a major effect on the C-M and C-C relations. Especially for the redshift $z\sim 2.0$, the $g-$band magnitude becomes smaller by $1.5$\,mag, the $g-r$ color becomes bluer by $1.0$\,mag, and the $u-g$ color becomes redder by $1.0$\,mag when binary interactions are included.

In the version of the SD scenario of SNe Ia studied here, a CO WD explodes close to the Chandrasekhar limit after accreting material from a non-degenerate He companion. In the present study, we employ the Stellar GADGET code to perform 3D hydrodynamical simulations of the interaction of the SN Ia ejecta with the He companion taking into account its orbital motion and spin. It is found that only 2%--5% of the initial companion mass are stripped off from the outer layers of He companions due to the SN impact. The dependence of the unbound mass (or the kick velocity) on the orbital separation can be fitted in good approximation by a power law for a given companion model. After the SN impact, the outer layers of a He donor star are significantly enriched with heavy elements from the low-expansion-velocity tail of SN Ia ejecta. The total mass of accumulated SN-ejecta material on the companion surface reaches about > 10e-3 M_sun for different companion models. This enrichment with heavy elements provides a potential way to observationally identify the surviving companion star in SN remnants. Finally, by artificially adjusting the explosion energy of the W7 explosion model, we find that the total accumulation of SN ejecta on the companion surface is also dependent on the explosion energy with a power law relation in good approximation.

It has recently been proposed that one sub-class of type Ia supernovae (SNe Ia) is sufficiently both distinct and common to be classified separately from the bulk of SNe Ia, with a suggested class name of "type Iax supernovae" (SNe Iax), after SN 2002cx. We show that the population properties of this class can be understood if the events originate from helium double-detonation sub-Chandrasekhar mass explosions, in which a carbon--oxygen white dwarf (CO WD) accumulates a helium layer from a non-degenerate helium star. We have incorporated detailed binary evolution calculations for the progenitor systems into a binary population synthesis model to obtain birthrates and delay times for such events. The predicted Galactic event rate is $\sim$$0.6-1.8\times 10^{-3}\,{\rm yr}^{-1}$, in good agreement with the measured rates of SNe Iax. In addition, predicted delay times are $\sim$70\,Myr$-$710\,Myr, consistent with the fact that SNe Iax have so far only been discovered in late-type galaxies. Based on the CO WD mass at explosion and previous detonation models, we also estimate the distribution of resulting SN brightness ($-13 \gtrsim M_{\rm bol} \gtrsim -19$ mag), which can reproduce the empirical diversity of SNe Iax. We speculate on why binaries with non-degenerate donor stars might lead to SNe Iax if similar systems with degenerate donors do not and suggest that the higher mass of the helium layer necessary for ignition at the lower accretion rates typically delivered from non-degenerate donors might provide the explanation.

More than about 50% stars are in binaries, but most stellar population studies take single star stellar population (ssSSP) models, which do not take binary interactions into account. In fact, the integrated peculiarities of ssSSPs are various from those of stellar populations with binary interactions (bsSSPs). Therefore, it is necessary to investigate the effects of binary interactions on the Lick indices and colours of populations detailedly. We show some formulae for calculating the difference between the Lick indices and colours of bsSSPs, and those of ssSSPs. Twenty-five Lick indices and 12 colours are studied in the work. The results can be conveniently used for calculating the effects of binary interactions on stellar population studies and for adding the effects of binary interactions into present ssSSP models. The electrical data and fortran procedures of the paper can be obtained on request to the authors.

Long-orbital-period subdwarf B (sdB) stars with main-sequence companions are believed to be the product of stable Roche Lobe overflow (RLOF), a scenario challenged by recent observations. Here we represent the results of a systematic study of the orbital-period distribution of sdB binaries in this channel using detailed binary evolution calculations. We show that the observed orbital-period distribution of long-period sdB binaries can be well explained by this scenario. Furthermore, we find that, if the progenitors of the sdB stars have initial masses below the helium flash mass, the sdB binaries produced from stable RLOF follow a unique mass -- orbital period relation for a given metallicity $Z$; increasing the orbital period from $\sim 400$ to $\sim 1100$\,d corresponds to increasing the mass of the sdB star from $\sim 0.40$ to $\sim 0.49\,M_\odot$ for $Z=0.02$. We suggest that the longest sdB binaries (with orbital period $> 1100$\,d) could be the result of atmospheric RLOF. The mass -- orbital period relation can be tested observationally if the mass of the sdB star can be determined precisely, e.g.\ from asteroseismology. Using this relation, we revise the orbital period distribution of sdB binaries produced by the first stable RLOF channel for the best fitting model of Han et al (2003), and show that the orbital period has a peak around 830\,d.

The calibrations of star formation rate (SFR) are prone to be affected by many factors, such as metallicity, initial mass function (IMF), evolutionary population synthesis (EPS) models and so on. In this paper we will discuss the effects of binary interactions, metallicity, EPS models and IMF on several widely used SFR calibrations based on the EPS models of Yunnan with and without binary interactions, BC03, SB99, PEGASE and POPSTAR. The inclusion of binary interactions makes these SFR conversion coefficients smaller (less than 0.2dex), and these differences increase with metallicity. The differences in the calibration coefficient between SFR and the luminosity of $\rm H\alpha$ recombination line (C$_{\rm H\alpha}$) and that between SFR and the ultraviolet (UV) fluxes at 1500 and 2800\,$\rm \AA$ (C$_{i, {\rm UV}}$), caused by IMF, are independent of metallicity (0.03-0.33\,dex) except $\Delta$C$_{\rm H\alpha, IMF}$ when using the POPSTAR and $\Delta$C$_{i, {\rm UV, IMF}}$ when using the PEGASE models. Moreover, we find that $L_{\rm 2800}$ is not suitable to the linear calibration of SFR at low metallicities. At last, we compare the effects of these several factors on the SFR calibrations considered in this paper. The effects of metallicity/IMF and EPS models on the C$_{\rm H\alpha}$ and C$_{\rm FIR}$ (the conversion coefficient between SFR and the far-infrared flux) are the largest among these factors, respectively. For the calibration between SFR and C$_{i, {\rm UV}}$, the effects of these several factors are comparable.

It is believed that, except for metallicity, some other parameters are needed to explain the horizontal branch (HB) morphology of globular clusters (GCs). Furthermore, these parameters are considered to be correlated with the mass loss of the red giant branch (RGB) stars. In our previous work, we proposed that tidally enhanced stellar wind during binary evolution may affect the HB morphology by enhancing the mass loss of the red giant primary. As a further study, we now investigate the effects of metallicity and age on HB morphology by considering tidally enhanced stellar winds during binary evolution. We incorporated the tidally enhanced-stellar-wind model into Eggleton's stellar evolution code to study the binary evolution. To study the effects of metallicity and age on our final results, we conducted two sets of model calculations: (i) for a fixed age, we used three metallicities, namely Z=0.0001, 0.001, and 0.02. (ii) For a fixed metallicity, Z=0.001, we used five ages in our model calculations: 14, 13, 12, 10, and 7 Gyr. We found that HB morphology of GCs becomes bluer with decreasing metallicity, and old GCs present bluer HB morphology than young ones. These results are consistent with previous work. Although the envelope-mass distributions of zero-age HB stars produced by tidally enhanced stellar wind are similar for different metallicities, the synthetic HB under tidally enhanced stellar wind for Z=0.02 presented a distinct gap between red and blue HB. However, this feature was not seen clearly in the synthetic HB for Z=0.001 and 0.0001. We also found that higher binary fractions may make HB morphology become bluer, and we discussed the results with recent observations.

Type Ia supernovae (SN Ia) are the most important standard candles for measuring the expansion history of the universe. The thermonuclear explosion of a white dwarf can explain their observed properties, but neither the progenitor systems nor any stellar remnants have been conclusively identified. Underluminous SN Ia have been proposed to originate from a so-called double-detonation of a white dwarf. After a critical amount of helium is deposited on the surface through accretion from a close companion, the helium is ignited causing a detonation wave that triggers the explosion of the white dwarf itself. We have discovered both shallow transits and eclipses in the tight binary system CD-30 11223 composed of a carbon/oxygen white dwarf and a hot helium star, allowing us to determine its component masses and fundamental parameters. In the future the system will transfer mass from the helium star to the white dwarf. Modelling this process we find that the detonation in the accreted helium layer is sufficiently strong to trigger the explosion of the core. The helium star will then be ejected at so large a velocity that it will escape the Galaxy. The predicted properties of this remnant are an excellent match to the so-called hypervelocity star US 708, a hot, helium-rich star moving at more than 750 km/s, sufficient to leave the Galaxy. The identification of both progenitor and remnant provides a consistent picture of the formation and evolution of underluminous type Ia supernovae.

In the SD scenario of SNe Ia the companion survives the SN explosion and thus should be visible near the center of the SN remnant and may show some unusual features. A promising approach to test progenitor models of SNe Ia is to search for the companion in SNRs. Here we present the results of 3D hydrodynamics simulations of the interaction between the SN Ia blast wave and a MS companion taking into consideration its orbital motion and spin. The primary goal of this work is to investigate the rotation of surviving companions after SN Ia explosions in the WD+MS scenario. We use Eggleton's code including the optically thick accretion wind model to obtain realistic models of companions. The impact of the SN blast wave on these companions is followed in 3D hydrodynamic simulations employing the SPH code GADGET3. We find that the rotation of the companion does not significantly affect the amount of stripped mass and the kick velocity caused by the SN impact. However, in our simulations, the rotational velocity of the companion is significantly reduced to about 14% to 32% of its pre-explosion value due to the expansion of the companion and the fact that 55%-89% of the initial angular momentum is carried away by the stripped matter. Compared with the observed rotational velocity of the presumed companion star of Tycho's SN, Tycho G, of 6 km/s the final rotational velocity we obtain is still higher by at least a factor of two. Whether this difference is significant, and may cast doubts on the suggestion that Tycho G is the companion of SN 1572, has to be investigated in future studies. Based on binary population synthesis results we present, for the first time, the expected distribution of rotational velocities of companions after the explosion which may provide useful information for the identification of the surviving companion in observational searches in other historical SNRs.

Using Eggleton's stellar evolution code, we study the minimum mass ratio ($q_{\rm min}$) of W Ursae Majoris (W UMa) binaries that have different primary masses. It is found that the minimum mass ratio of W UMa binaries decreases with increasing mass of the primary if the primary's mass is less than about 1.3$M_{\rm \odot}$, and above this mass the ratio is roughly constant. By comparing the theoretical minimum mass ratio with the observational data, it is found that the existence of low-$q$ systems can be explained by the different structure of the primaries with different masses. This suggests that the dimensionless gyration radius ($k_1^2$) and thus the structure of the primary is very important in determining the minimum mass ratio. In addition, we investigate the mass loss during the merging process of W UMa systems and calculate the rotation velocities of the single stars formed by the merger of W UMa binaries due to tidal instability. It is found that in the case of the conservation of mass and angular momentum, the merged single stars rotate with a equatorial velocity of about $588\sim819$ km s$^{-1}$, which is much larger than their break-up velocities ($v_{\rm b}$). This suggests that the merged stars should extend to a very large radius (3.7$\sim$5.3 times the radii of the primaries) or W UMa systems would lose a large amount of mass (21$\sim$33 per cent of the total mass) during the merging process. If the effect of magnetic braking is considered, the mass loss decreases to be 12$\sim$18 per cent of their total masses. This implies that the significant angular momentum and mass might be lost from W UMa systems in the course of the merging process, and this kind of mass and angular momentum loss might be driven by the release of orbital energy of the secondaries, which is similar to common-envelope evolution.

VFTS102 has a projected rotational velocity ($>$500\,km s$^{-1}$) and would appear to be the fastest rotating O-type star. We show that its high rotational velocity could be understood within the framework of the binary merger. In the binary merger channel, the progenitor binary of VFTS102 would evolve into contact while two components are still on the main sequence, and then merge into a rapidly rotating single star. Employing Eggleton's stellar evolution code, we performed binary stellar evolution calculations and mapped out the initial parameters of the progenitor of VFTS102 in the orbital period-mass ratio ($P-q$) plane. We found that the progenitor binary of VFTS102 with initial mass ratio $q_0 \la 0.7$ should have an initial orbital period shorter than $3.76-4.25$\,days, while above this mass ratio it should have an initial orbital period shorter than $1.44-1.55$\,days. The progenitor of VFTS102 would evolve into contact during the rapid mass transfer phase or during the subsequent slow mass transfer phase, and might ultimately merge into a rapidly rotating massive star. In addition, we performed Monte Carlo simulations to investigate the binary merger channel. We estimated the fraction of binaries that would merge into single stars and the fraction of single stars that might be produced from the binary merger channel. It is found that about 8.7% of binaries would evolve into contact and merge into rapidly rotating single stars, and about 17.1% of single stars might be produced from the binary merger channel and should have similar properties to VFTS102. This suggests that the binary merger channel might be one of the main channels for the formation of rapidly rotating massive stars like VFTS102.

Metallicity is the first parameter to influence the horizontal branch (HB) morphology of globular clusters (GCs). It has been found, however, that some other parameters may also play an important role in affecting the morphology. While the nature of these important parameters remains unclear, they are believed to be likely correlated with wind mass-loss of red giants, since this mass loss determines their subsequent locations on the HB. Unfortunately, the mass loss during the red giant stages of the stellar evolution is poorly understood at present. The stellar winds of red giants may be tidally enhanced by companion stars if they are in binary systems. We investigate evolutionary consequences of red giants in binaries by including tidally enhanced stellar winds, and examine the effects on the HB morphology of GCs. We find that red, blue, and extreme horizontal branch stars are all produced under the effects of tidally enhanced stellar wind without any additional assumptions on the mass-loss dispersion. Furthermore, the horizontal branch morphology is found to be insensitive to the tidal enhancement parameter, Bw. We compare our theoretical results with the observed horizontal branch morphology of globular cluster NGC 2808, and find that the basic morphology of the horizontal branch can be well reproduced. The number of blue horizontal branch stars in our calculations, however, is lower than that of NGC 2808.

This work aims to present our current best physical understanding of common-envelope evolution (CEE). We highlight areas of consensus and disagreement, and stress ideas which should point the way forward for progress in this important but long-standing and largely unconquered problem. Unusually for CEE-related work, we mostly try to avoid relying on results from population synthesis or observations, in order to avoid potentially being misled by previous misunderstandings. As far as possible we debate all the relevant issues starting from physics alone, all the way from the evolution of the binary system immediately before CEE begins to the processes which might occur just after the ejection of the envelope. In particular, we include extensive discussion about the energy sources and sinks operating in CEE, and hence examine the foundations of the standard energy formalism. Special attention is also given to comparing the results of hydrodynamic simulations from different groups and to discussing the potential effect of initial conditions on the differences in the outcomes. We compare current numerical techniques for the problem of CEE and also whether more appropriate tools could and should be produced (including new formulations of computational hydrodynamics, and attempts to include 3D processes within 1D codes). Finally we explore new ways to link CEE with observations. We compare previous simulations of CEE to the recent outburst from V1309 Sco, and discuss to what extent post-common-envelope binaries and nebulae can provide information, e.g. from binary eccentricities, which is not currently being fully exploited.

We build the Yunnan-III evolutionary population synthesis (EPS) models by using the MESA stellar evolution code, BaSeL stellar spectra library and the initial mass functions (IMFs) of Kroupa and Salpeter, and present colours and integrated spectral energy distributions (ISEDs) of solar-metallicity stellar populations (SPs) in the range of 1Myr-15 Gyr. The main characteristic of the Yunnan-III EPS models is the usage of a set of self-consistent solar-metallicity stellar evolutionary tracks (the masses of stars are from 0.1 to 100Msun). This set of tracks is obtained by using the state-of-the-art MESA code. MESA code can evolve stellar models through thermally pulsing asymptotic giant branch (TP-AGB) phase for low- and intermediate-mass stars. By comparisons, we confirm that the inclusion of TP-AGB stars make the V-K, V-J and V-R colours of SPs redder and the infrared flux larger at ages log(t/yr)>7.6 (the differences reach the maximum at log(t/yr)~8.6, ~0.5-0.2mag for colours, ~2 times for K-band flux). The stellar evolutionary tracks, isochrones, colours and ISEDs can be obtained on request from the first author or from our website (http://www1.ynao.ac.cn/~zhangfh/). Using the isochrones, you can build your EPS models. Now the format of stellar evolutionary tracks is the same as that in the STARBURST99 code, you can put them into the STARBURST99 code and get the SP's results. Moreover, the colours involving other passbands or on other systems (for example, HST $F439W-F555W$ colour on AB system) can also be obtained on request.

The identity of the progenitor systems of SNe Ia is still uncertain. In the single-degenerate (SD) scenario, the interaction between the SN blast wave and the outer layers of a main sequence (MS) companion star strips off H-rich material which is then mixed into the ejecta. Strong contamination of the SN ejecta with stripped material could lead to a conflict with observations of SNe Ia. This constrains the SD progenitor model. In this work, our previous simulations based on simplified progenitor donor stars have been updated by adopting more realistic progenitor-system models that result from fully detailed, state-of-the-art binary evolution calculations. We use Eggleton's stellar evolution code including the optically thick accretion wind model and the possibility of the effects of accretion disk instabilities to obtain realistic models of companions for different progenitor systems. The impact of the SN blast wave on these companion stars is followed in three-dimensional hydrodynamic simulations employing the SPH code GADGET3. We find that the stripped masses range from 0.11 to 0.18 M_sun. The kick velocity is between 51 and 105 km/s. We find that the stripped mass and kick velocity depend on the ratio of the orbital separation to the radius of a companion. They can be fitted by a power law for a given companion model. However, the structure of the companion star is also important for the amount of stripped material. With more realistic companion star models than in previous studies, our simulations show that the H masses stripped from companions are inconsistent with the best observational limits (< 0.01 M_sun) derived from nebular spectra. However, a rigorous forward modeling based on impact simulations with radiation transfer is required to reliably predict observable signatures of the stripped H and to conclusively assess the viability of the considered SN Ia progenitor scenario.

We construct a parameterized model to explore the main properties of the star formation history of M33. We assume that the disk originates and grows by the primordial gas infall and adopt the simple form of gas accretion rate with one free parameter, the infall time-scale. We also include the contribution of gas outflow process. A major update of the model is that we adopt a molecular hydrogen correlated star formation law and calculate the evolution of the atomic and molecular gas separately. Comparisons between the model predictions and the observational data show that the model predictions are very sensitive to the adopted infall time-scale, while the gas outflow process mainly influences the metallicity profile. The model adopting a moderate outflow rate and an inside-out formation scenario can be in good agreement with most of observed constraints of M33 disk. We also compare the model predictions based on the molecular hydrogen correlated star formation law and that based on the Kennicutt star formation law. Our results imply that the molecular hydrogen correlated star formation law should be preferred to describe the evolution of the M33 disk, especially the radial distributions of both the cold gas and the stellar population.

Recent theoretical and observational studies both argue that the merging of double carbon-oxygen white dwarfs (WDs) is responsible for at least some Type Ia supernovae (SNe Ia). Previous (standard) studies of the anticipated SN birthrate from this channel have assumed that the merger process is conservative and that the primary criterion for explosion is that the merged mass exceeds the Chandrasekhar mass. Han & Webbink (1999) demonstrated that mass transfer and merger in close double WDs will in many cases be non-conservative. Pakmor et al. (2011) further suggested that the merger process should be violent in order to initiate an explosion. We have therefore investigated how the SN Ia birthrate from the double-degenerate (DD) channel is affected by these constraints. Using the binary-star population-synthesis method, we have calculated the DD SN Ia birthrate under conservative and non-conservative approximations, and including lower mass and mass-ratio limits indicated by recent smoothed-particle-hydrodynamic calculations. The predicted DD SN Ia rate is significantly reduced by all of these constraints. With dynamical mass loss alone (violent merger) the birthrate is reduced to 56% of the conservative rate. Requiring that the mass ratio $q>2/3$ further reduces the birthrate to 18% that of the standard assumption. An upper limit of 0.0061 SNuM, or a Galactic rate of $4.6 \times 10^{-4}{\rm yr}^{-1}$, might be realistic.

The weather at Xinglong Observing Station, where the Guo Shou Jing Telescope (GSJT) is located, is strongly affected by the monsoon climate in north-east China. The LAMOST survey strategy is constrained by these weather patterns. In this paper, we present a statistics on observing hours from 2004 to 2007, and the sky brightness, seeing, and sky transparency from 1995 to 2011 at the site. We investigate effects of the site conditions on the survey plan. Operable hours each month shows strong correlation with season: on average there are 8 operable hours per night available in December, but only 1-2 hours in July and August. The seeing and the sky transparency also vary with seasons. Although the seeing is worse in windy winters, and the atmospheric extinction is worse in the spring and summer, the site is adequate for the proposed scientific program of LAMOST survey. With a Monte Carlo simulation using historical data on the site condition, we find that the available observation hours constrain the survey footprint from 22h to 16h in right ascension; the sky brightness allows LAMOST to obtain the limit magnitude of V = 19.5mag with S/N = 10.

We describe the footprint and input catalog for bright nights in the LAMOST Pilot Survey, which began in October 2011. Targets are selected from two stripes in the north and south Galactic Cap regions, centered at $\alpha$= 29$^\circ$, with 10$^\circ$ width in declination, covering right ascension of 135$^\circ-290^\circ$ and -30$^\circ$ to 30$^\circ$ respectively. We selected spectroscopic targets from a combination of the SDSS and 2MASS point source catalogs. The catalog of stars defining the field centers (as required by the Shack-Hartmann wavefront sensor at the center of the LAMOST field) consists of all V < 8m stars from the Hipparcos catalog. We employ a statistical selection algorithm that assigns priorities to targets based on their positions in multidimensional color/magnitude space. This scheme overemphasizes rare objects and de-emphasizes more populated regions of magnitude and color phase space, while ensuring a smooth, well-understood selection function. A demonstration of plate design is presented based on the Shack-Hartmann star catalog and an input catalog that was generated by our target selection routines.

We describe the current plans for a spectroscopic survey of millions of stars in the Milky Way galaxy using the Guo Shou Jing Telescope (GSJT, formerly the Large Area Multi-Object Spectroscopic Telescope - LAMOST). The survey will obtain spectra for 2.5 million stars brighter than $r<19$ during dark/grey time, and 5 million stars brighter than $r<17$ or $J<16$ on nights that are moonlit or have low transparency. The survey will begin in fall of 2012, and will run for at least four years. The telescope design constrains the optimal declination range for observations to $10^\circ<\delta<50^\circ$, and site conditions lead to an emphasis on stars in the direction of the Galactic anticenter. The survey is divided into three parts with different target selection strategies: disk, anticenter, and spheroid. The resulting dataset will be used to study the merger history of the Milky Way, the substructure and evolution of the disks, the nature of the first generation of stars through identification of the lowest metallicity stars, and star formation through study of open clusters and the OB associations. Detailed design of the LEGUE survey will be completed after a review of the results of the pilot survey in summer 2012.

We outline the design of the dark nights portion of the LAMOST Pilot Survey, which began observations in October 2011. In particular, we focus on Milky Way stellar candidates that are targeted for the LEGUE (LAMOST Experiment for Galactic Understanding and Exploration) survey. We discuss the regions of sky in which spectroscopic candidates were selected, and the motivations for selecting each of these sky areas. Some limitations due to the unique design of the telescope are discussed, including the requirement that a bright (V < 8) star be placed at the center of each plate for wavefront sensing and active optics corrections. The target selection categories and scientific goals motivating them are briefly discussed, followed by a detailed overview of how these selection functions were realized. We illustrate the difference between the overall input catalog - Sloan Digital Sky Survey (SDSS) photometry - and the final targets selected for LAMOST observation.

We describe the target selection algorithm for the low latitude disk portion of the LAMOST Pilot Survey, which aims to test systems in preparation for the LAMOST spectroscopic survey. We use the PPMXL (Roeser et al. 2010) astrometric catalog, which provides positions, proper motions, B/R/I magnitudes (mostly) from USNO-B (Monet et al. 2003) and J/H/Ks from The Two Micron All Sky Survey (2MASS, see Skrutskie et al. 2006) as well. We chose 8 plates along the Galactic plane, in the region $0^\circ<\alpha<67^\circ$ and $42^\circ<\delta<59^\circ$, that cover 22 known open clusters with a range of ages. Adjacent plates may have small overlapping. Each plate covers an area $2.5^\circ$ in radius,with central star (for Shack-Hartmann guider) brighter than $\sim8^{\rm th}$ magnitude. For each plate, we create an input catalog in the magnitude range $11.3<Imag<16.3$ and $Bmag$ available from PPMXL. The stars are selected to satisfy the requirements of the fiber positioning system and have a uniform distribution in the $I$ vs. $B-I$ color-magnitude diagram. Our final input catalog consists of 12,000 objects on each of 8 plates that are observable during the winter observing season in Xinglong Station of the National Astronomical Observatory of China.

We describe a general target selection algorithm that is applicable to any survey in which the number of available candidates is much larger than the number of objects to be observed. This routine aims to achieve a balance between a smoothly-varying, well-understood selection function and the desire to preferentially select certain types of targets. Some target-selection examples are shown that illustrate different possibilities of emphasis functions. Although it is generally applicable, the algorithm was developed specifically for the LAMOST Experiment for Galactic Understanding and Exploration (LEGUE) survey that will be carried out using the Chinese Guo Shou Jing Telescope. In particular, this algorithm was designed for the portion of LEGUE targeting the Galactic halo, in which we attempt to balance a variety of science goals that require stars at fainter magnitudes than can be completely sampled by LAMOST. This algorithm has been implemented for the halo portion of the LAMOST pilot survey, which began in October 2011.

Type Ia supernovae (SNe Ia) play an important role in astrophysics and are crucial for the studies of stellar evolution, galaxy evolution and cosmology. They are generally thought to be thermonuclear explosions of accreting carbon-oxygen white dwarfs (CO WDs) in close binaries, however, the nature of the mass donor star is still unclear. In this article, we review various progenitor models proposed in the past years and summarize many observational results that can be used to put constraints on the nature of their progenitors. We also discuss the origin of SN Ia diversity and the impacts of SN Ia progenitors on some fields. The currently favourable progenitor model is the single-degenerate (SD) model, in which the WD accretes material from a non-degenerate companion star. This model may explain the similarities of most SNe Ia. It has long been argued that the double-degenerate (DD) model, which involves the merger of two CO WDs, may lead to an accretion-induced collapse rather than a thermonuclear explosion. However, recent observations of a few SNe Ia seem to support the DD model, and this model can produce normal SN Ia explosion under certain conditions. Additionally, the sub-luminous SNe Ia may be explained by the sub-Chandrasekhar mass model. At present, it seems likely that more than one progenitor model, including some variants of the SD and DD models, may be required to explain the observed diversity of SNe Ia.

We present a detailed comparison between the 2-10 keV hard X-ray and infrared (IR) luminosity function (LF) of active galactic nuclei (AGN). The composite X-ray to IR spectral energy distributions (SEDs) of AGN used for connecting the hard X-ray LF (HXLF) and IR LF (IRLF) are modeled with a simple but well tested torus model based on the radiative transfer and photoionization code CLOUDY. Four observational determinations of the evolution of 2-10 keV HXLF and six evolution models of the obscured type-2 AGN fraction ($f_2$) have been considered. The 8.0 and 15 \micron LFs for the total, unobscured type-1 and obscured type-2 AGN are predicted from the HXLFs, and then compared with the measurements currently available. We find that the IRLFs predicted from HXLFs tend to underestimate the number of the most IR-luminous AGN. This is independent of the choices of HXLF and $f_2$, and even more obvious for the HXLFs recently measured. We show that the discrepancy between the HXLFs and IRLFs can be largely resolved when the anticorrelation between the UV to X-ray slope $\alpha_{\mathrm{ox}}$ and UV luminosity $L_{\rm UV}$ is appropriately considered. We also discuss other possible explanations for the discrepancy, such as the missing population of Compton-thick AGN and possible contribution of star-formation in the host to the mid-IR. Meanwhile, we find that the HXLFs and IRLFs of AGN can be more consistent with each other if the obscuration mechanisms of quasars and Seyferts are assumed to be different, corresponding to their different triggering and fueling mechanisms. More accurate measurements of the IRLFs of AGN, especially that determined at smaller redshift bins and more accurately separated to that for type-1 and type-2, are very helpful for clarifying these interesting issues.

We present BayeSED, a general purpose tool for doing Bayesian analysis of SEDs by using whatever pre-existing model SED libraries or their linear combinations. The artificial neural networks (ANNs), principal component analysis (PCA) and multimodal nested sampling (MultiNest) techniques are employed to allow a highly efficient sampling of posterior distribution and the calculation of Bayesian evidence. As a demonstration, we apply this tool to a sample of hyperluminous infrared galaxies (HLIRGs). The Bayesian evidences obtained for a pure Starburst, a pure AGN, and a linear combination of Starburst+AGN models show that the Starburst+AGN model have the highest evidence for all galaxies in this sample. The Bayesian evidences for the three models and the estimated contributions of starburst and AGN to infrared luminosity show that HLIRGs can be classified into two groups: one dominated by starburst and the other dominated by AGN. Other parameters and corresponding uncertainties about starburst and AGN are also estimated by using the model with the highest Bayesian evidence. We found that the starburst region of the HLIRGs dominated by starburst tends to be more compact and has a higher fraction of OB star than that of HLIRGs dominated by AGN. Meanwhile, the AGN torus of the HLIRGs dominated by AGN tend to be more dusty than that of HLIRGs dominated by starburst. These results are consistent with previous researches, but need to be tested further with larger samples. Overall, we believe that BayeSED could be a reliable and efficient tool for exploring the nature of complex systems such as dust-obscured starburst-AGN composite systems from decoding their SEDs.

The number of methods used to study the properties of galaxies is increased, and testing these methods is very important. Galactic globular clusters (GCs) provide an excellent medium for such test, because they can be considered as simple stellar populations. We present ages and metallicities for 40 Galactic GCs as determined from three publicly available techniques, including colour, Lick-index and spectrum-fitting methods, based on Bruzual & Charlot evolutionary population synthesis (EPS) models. By comparing with the ages obtained from colour-magnitude diagrams (CMDs) and metallicities obtained from spectra of stars, we are able to estimate the ability of 'these methods on determination of GCs parameters, which is absolutely necessary. As a result, we find that: (i) for the metallicity, our derived metallicities agree with those derived from the spectra of stars, Lick-index method is suitable to study metallicity for the stellar population systems in the range of -1.5=<[Fe/H]=<-0.7 and spectrum- fitting method is suitable to study metallicity for the stellar population systems in the range of -2.3=<[Fe/H]=<-1.5; (ii) for the age, these three methods have difficulties in age determination, our derived ages are smaller (about 2.0 Gyr, on average) than the results of CMDs for all these three methods. We use Vazdekis and Maraston models to analyze whether our results are dependent on EPS models, and find that the tendency of these two models is the same as that of Bruzual & Charlot models. Our results are independent of the EPS models. In addition, our test is based on the old GCs and our conclusions may hold for old stellar population systems. The whole abstract can be found in my PDF version.

Using the evolutionary population synthesis (EPS) models with and without binary interactions (BIs), we present L_Ha, L_[OII], L_{i,UV} and L_FIR for Burst, S0, Sa-Sd and Irr galaxies, and present the calibrations of star formation rate (SFR) in terms of these diagnostics. By comparison, we find that BIs lower the SFR.vs.L_Ha and SFR.vs.L_[OII] conversion factors by ~0.2 dex. Moreover, BIs do not significantly vary the calibrations of SFR in terms of L_{i,UV}. In addition, BIs have little effect on the flux at 2800A. At last, the calibration of SFR from L_FIR is almost unaffected by BIs. This is caused by the fact that BIs almost do not affect the bolometric magnitudes of SPs. We also discuss the effects of initial mass function (IMF), gas-recycle assumption (GRA) and EPS models on these SFR calibrations. Comparing the results by using Salpeter (S55) IMF with those by using Miller & Scalo (MS79) IMF, we find that the SFR.vs.L_Ha and SFR.vs.L_[OII] conversion factors by using S55 IMF are greater by 0.4 and 0.2 dex than those by using MS79 IMF for the models with and without BIs, respectively. The SFR.vs.L_{i,UV} and SFR.vs.L_FIR conversion factors by using S55 IMF are larger by an amount of 0.2 dex than the corresponding ones by using MS79 IMF. The inclusion of GRA only lowers these SFR calibrations at faint SFR. Moreover, comparing the results when using different EPS models, we find that the differences in the SFR.vs.L_Ha and SFR.vs.L_[OII] conversion factors reach ~ 0.7 and 0.9 dex, the difference in the SFR.vs.L_FIR conversion factor reaches 0.4 and 0.8 dex, and the differences in the SFR.vs.L_{i,UV} conversion factors reach 0.3 and 0.2 dex when using S55 and NON-S55 IMF (including Cha03, K01, K93' and MS79 IMFs, partly caused by the difference in the IMF), respectively. At last, we give the conversion coefficients between SFR and these diagnostics for all models.

The stability of mass transfer is important in the formation of contact binaries from detached binaries when the primaries of the initially detached binaries fill their Roche lobes. Using Eggleton's stellar evolution code, we investigate the formation and the short-period limit of contact binaries by considering the effect of the instability of mass transfer. It is found that with decreasing initial primary mass from 0.89M$_{\rm \odot}$ to 0.63M$_{\rm \odot}$, the range of the initial mass ratio decreases for detached binaries that experience stable mass transfer and evolve into contact. If the initial primary mass is less than 0.63M$_{\rm \odot}$, detached binaries would experience dynamically unstable mass transfer when the primaries of detached binaries fill their Roche lobes. These systems would evolve into a common envelope situation and probably then to a complete merger of two components on a quite short timescale. This results in a low mass limit at about 0.63M$_{\rm \odot}$ for the primary mass of contact binaries, which might be a main reason why the period distribution of contact binaries has a short limit of about 0.22 days. By comparing the theoretical period distribution of contact binaries with the observational data, it is found that the observed contact binaries are above the low mass limit for the primary mass of contact binaries and no observed contact binaries are below this limit. This suggests that the short-period limit of contact binaries can be explained by the instability of the mass transfer that occurs when the primaries of the initially detached binaries fill their Roche lobes.

Employing Eggleton's stellar evolution code with the optically thick wind assumption, we systematically studied the He star donor channel of type Ia supernovae (SNe Ia), in which a carbon-oxygen white dwarf (WD) accretes material from a He main-sequence star or a He subgiant to increase its mass to the Chandrasekhar mass. We mapped out the initial parameters for producing SNe Ia in the orbital period--secondary mass plane for various WD masses from this channel. Based on a detailed binary population synthesis approach, we find that this channel can produce SNe Ia with short delay times (\sim100Myr) implied by recent observations. We obtained many properties of the surviving companions of this channel after SN explosion, which can be verified by future observations. We also find that the surviving companions from the SN explosion scenario have a high spatial velocity (>400km/s), which could be an alternative origin for hypervelocity stars (HVSs), especially for HVSs such as US 708.

In the symbiotic (or WD+RG) channel of the single degenerate scenario for type Ia supernovae (SNe Ia) the explosions occur a relatively long time after star formation. The birthrate from this channel would be too low to account for all observed SNe Ia were it not for some mechanism to enhance the rate of accretion on to the white dwarf. A tidally enhanced stellar wind, of the type which has been postulated to explain many phenomena related to giant star evolution in binary systems, can do this. Compared to mass-stripping this model extends the space of SNe Ia progenitors to longer orbital periods and hence increases the birthrate to about $0.0069\,{\rm yr^{-1}}$ for the symbiotic channel. Two symbiotic stars, T CrB and RS Oph, considered to be the most likely progenitors of SNe Ia through the symbiotic channel, are well inside the period--companion mass space predicted by our models.

We present high quality long slit spectra along the major and minor axes out to 1.5 effective radius ($R_e$) of the massive galaxy NGC4636 taken by Hobby-Eberly Telescope (HET). Using Fourier Correlation Quotient (FCQ) method, we measured the stellar line-of-sight velocity distribution along the axes. Furthermore, six Lick/IDS indices ($H\beta,Mgb,Fe_{5015},Fe_{5270},Fe_{5335},Fe_{5406}$) are derived from the clean spectrum. By comparing the measured absorption line strengths with the predictions of Simple Stellar Populations (SSP) models, we derived ages, total metallicity and $\alpha$ abundance profiles of the galaxy. This galaxy presents old and $[\alpha/Fe]$ over abundant stellar populations. Indeed, using the SSP model, we obtained the broadband color profiles. The theoretical colors match well with the measured colors and present red sharp peaks at the galaxy center. The sharp peaks of the colors are mainly shaped by the high metallicity in the galaxy center. Interestingly, the galaxy has steep negative metallicity gradients, but trend flattens outwards. This result likly suggests that the center and outer regions of the galaxy formed through different formation process.

Recent discovery of several overluminous type Ia supernovae (SNe Ia) indicates that the explosive masses of white dwarfs may significantly exceed the canonical Chandrasekhar mass limit. Rapid differential rotation may support these massive white dwarfs. Based on the single-degenerate scenario, and assuming that the white dwarfs would differentially rotate when the accretion rate $\dot{M}>3\times 10^{-7}M_{\odot}\rm yr^{-1}$, employing Eggleton's stellar evolution code we have performed the numerical calculations for $\sim$ 1000 binary systems consisting of a He star and a CO white dwarf (WD). We present the initial parameters in the orbital period - helium star mass plane (for WD masses of $1.0 M_{\odot}$ and $1.2 M_{\odot}$, respectively), which lead to super-Chandrasekhar mass SNe Ia. Our results indicate that, for an initial massive WD of $1.2 M_{\odot}$, a large number of SNe Ia may result from super-Chandrasekhar mass WDs, and the highest mass of the WD at the moment of SNe Ia explosion is 1.81 $M_\odot$, but very massive ($>1.85M_{\odot}$) WDs cannot be formed. However, when the initial mass of WDs is $1.0 M_{\odot}$, the explosive masses of SNe Ia are nearly uniform, which is consistent with the rareness of super-Chandrasekhar mass SNe Ia in observations.

The nature of type Ia supernovae (SNe Ia) is still unclear. Employing Eggleton's stellar evolution code with the optically thick wind assumption, we systematically studied the He star donor channel of SNe Ia, in which a carbon-oxygen white dwarf accretes material from a He main-sequence star or a He subgiant to increase its mass to the Chandrasekhar mass. We mapped out the initial parameters for producing SNe Ia in the orbital period--secondary mass plane for various WD masses from this channel. According to a detailed binary population synthesis approach, we find that this channel can produce SNe Ia with short delay times (~100Myr) implied by recent observations. We obtained many properties of the surviving companions of this channel after SN explosion, which can be verified by future observations. We also find that the surviving companions from the SN explosion scenario have a high spatial velocity (>400km/s), which could be an alternative origin for hypervelocity stars (HVSs), especially for HVSs such as US 708.

We propose a formation channel for the previously unexplained helium-rich subdwarf O (He-rich sdO) stars in which post-subdwarf B (sdB) stars (i.e. hybrid COHe white dwarfs) reignite helium burning in a shell after gaining matter from their helium white-dwarf (WD) companions. Such short-period binaries containing post-sdB WDs and helium WDs are predicted by one of the major binary formation channels for sdB stars. In the majority of cases, mass transfer is expected to lead to a dynamically unstable merger event, leaving a single-star remnant. Calculations of the evolution of these stars show that their properties are consistent with the observed He-rich sdO stars. The luminosity of these stars is about an order of magnitude higher than that of canonical sdB stars. We also suggest that binary systems such as PG 1544+488 (Ahmad et al. 2004) and HE 0301-3039 (Lisker et al. 2004), which each contain two hot subdwarfs, could be the outcome of a double-core common-envelope phase. Since this favours intermediate-mass progenitors, this may also explain why the subdwarfs in these systems are He-rich.

Using Yunnan evolutionary population synthesis (EPS) models, we present integrated colours, integrated spectral energy distributions (ISEDs) and absorption-line indices defined by the Lick Observatory image dissector scanner (Lick/IDS) system, for an extensive set of instantaneous-burst binary stellar populations (BSPs) with interactions. By comparing the results for populations with and without interactions we show that the inclusion of binary interactions makes the appearance of the population substantially bluer. This effect raises the derived age and metallicity of the population. To be used in the studies of modern spectroscopic galaxy surveys at intermediate/high spectral resolution, we also present intermediate- (3A) and high-resolution (~0.3A) ISEDs and Lick/IDS absorption-line indices for BSPs. To directly compare with observations the Lick/IDS absorption indices are also presented by measuring them directly from the ISEDs.

Subdwarf B stars (sdBs) can significantly change the ultraviolet spectra of populations at age t~1Gyr, and have been even included in the volutionary population synthesis (EPS) models by Han et al. (2007). In this study we present the spectral energy distributions (SEDs) of binary stellar populations (BSPs) by combining the EPS models of Han et al. (2007) and those of the Yunnan group (Zhang et al. 2004, 005), which have included various binary interactions (except sdBs) in EPS models. This set of SEDs is vailable upon request from the authors. Using this set of SEDs of BSPs we build the spectra of Burst, E, S0-Sd and Irr types of galaxies by using the package of Bruzual & Charlot (2003, BC03). Combined with the photometric data (filters and magnitudes), we obtain the photometric redshifts and morphologies of 1502 galaxies by using the Hyperz code of Bolzonella et al. (2000). This sample of galaxies is obtained by removing those objects, mismatched with the SDSS/DR7 and GALEX/DR4, from the catalogue of Fukugita et al. (2007). By comparison the results with the SDSS spectroscopic redshifts and the morphological index of Fukugita et al. (2007), we find that the photo-zs fluctuate with the SDSS spectroscopic redshifts, while the Sa-Sc galaxies in the catalogue of Fukugita et al. (2007) are classified earlier as Burst-E galaxies.

Using the Hyperz code and a template spectral library which consists of 4 observed galaxy spectra from Coleman, Wu & Weedman (CWW, 1980) and 8 spectral families built with evolutionary population synthesis models, we present photometric redshift estimates (photo-z) for a spectroscopic sample of 6,531 galaxies, and morphologies for a morphological sample of 1,502 bright galaxies. All galaxies are matched with the SDSS DR7 and GALEX DR4. The inclusion of Fuv or Nuv or both photometry decreases the number of catastrophic identifications (CIs, |z_phot -z_spec| > 1.0). If CIs are removed, the inclusion of both Fuv and Nuv photometry mainly increases the number of non-CIs in the low redshift, g-r < 0.8 and fainter r-magnitude regions. The inclusion of binary interactions (BIs) mainly increases the number of non-CIs and decreases the deviations in the 0.3 < g-r < 0.8 region in the case of only using optical photometry. The inclusion of UV photometry would decrease and increase the probability that early types are classified as Burst and E types, respectively, and increase that late types are classified as CWW-Sbc and CWW-Scd types. If CIs are excluded, the inclusion of UV data mainly raises the identifications of late types in all redshift, bluer g-r and r > 14 regions. Moreover, BIs mainly affect the determinations of E and S0 types. Nuv -u = 1.94 and 5.77-1.47(u-r) = Fuv discriminators can be used as morphology selection indicators. These two criteria have comparable reliability and completeness for selecting early- and late-type galaxies to C=2.6 criterion and higher completeness for early-type selection than u-r=2.22 criterion.

In our evolutionary population synthesis models, the samples of binaries are reproduced by the 'patched' Monte Carlo simulation and the stellar masses, integrated J, H, K, L, L2 and M magnitudes, mass-to-light ratios and broad colours involving infrared bands are presented, for an extensive set of instantaneous-burst binary stellar populations. In addition, the fluctuations in the integrated colours, which have been given by Zhang et al. (2005), are reduced. By comparing the results for binary stellar populations with (Model A) and without (Model B) binary interactions we show that the inclusion of binary interactions makes the stellar mass of a binary stellar population smaller (~3.6-4.5% during the past 15Gyr); magnitudes greater (except U, ~0.18mag at the most); colours smaller (~0.15mag for V-K at the most); the mass-to-light ratios greater (~0.06 for K-band) except those in the U and B passbands at higher metallicities. And, Binary interactions make the V magnitude less sensitive to age, R and I magnitudes more sensitive to metallicity. Given an age, the absolute values of the differences in the stellar mass, magnitudes, mass-to-light ratios (except those in the U and B bands) between Models A and B reach the maximum at Z=0.0001, i.e., the effects of binary interactions on these parameters reach the maximum, while the differences in some colours reach the maximum at Z ~0.01-0.0004. On the contrary, the absolute value of the difference in the stellar mass is minimal at Z=0.03, those in the U,B,V magnitudes and the mass-to-light ratios in the U and B bands reach the minimum at Z ~0.01-0.004.

We present the integrated J, H, K, L, M and N magnitudes and the colours involving infrared bands, for an extensive set of instantaneous-burst binary stellar populations (BSPs) by using evolutionary population synthesis (EPS). By comparing the results for BSPs WITH and WITHOUT binary interactions we show that the inclusion of binary interactions makes the magnitudes of populations larger (fainter) and the integrated colours smaller (bluer) for t > 1Gyr. Also, we compare our model magnitudes and colours with those of Bruzual & Charlot (2003, hereafter BC03) and Maraston (2005, hereafter M05). At last, we compare these model broad colours with Magellanic Clouds globular clusters (GCs) and Milky Way GCs. In (V-R)-[Fe/H] and (V-I)-[Fe/H] diagrams it seems that our models match the observations better than those of BC03 and M05.

Using the Hyperz code (Bolzonella et al. 2000) we present photometric redshift estimates for a random sample of galaxies selected from the SDSS/DR7 and GALEX/DR4, for which spectroscopic redshifts are also available. We confirm that the inclusion of ultraviolet photometry improves the accuracy of photo-zs for those galaxies with g*-r* < 0.7 and z_spec < 0.2. We also address the problem of how binary interactions can affect photo-z estimates, and find that their effect is negligible.

Roche-lobe overflow and common envelope evolution are very important in binary evolution, which is believed to be the main evolutionary channel to hot subdwarf stars. The details of these processes are difficult to model, but adiabatic expansion provides an excellent approximation to the structure of a donor star undergoing dynamical time scale mass transfer. We can use this model to study the responses of stars of various masses and evolutionary stages as potential donor stars, with the urgent goal of obtaining more accurate stability criteria for dynamical mass transfer in binary population synthesis studies. As examples, we describe here several models with the initial masses equal to 1 Msun and 10 Msun, and identify potential limitations to the use of our results for giant-branch stars.

The asymptotic response of donor stars in interacting binary systems to very rapid mass loss is characterized by adiabatic expansion throughout their interiors. In this limit, energy generation and heat flow through the stellar interior can be neglected. We model this response by constructing model sequences, beginning with a donor star filling its Roche lobe at an arbitrary point in its evolution, holding its specific entropy and composition profiles fixed as mass is removed from the surface. The stellar interior remains in hydrostatic equilibrium. Luminosity profiles in these adiabatic models of mass-losing stars can be reconstructed from the specific entropy profiles and their gradients. These approximations are validated by comparison with time-dependent binary mass transfer calculations. We describe how adiabatic mass loss sequences can be used to quantify threshold conditions for dynamical time scale mass transfer, and to establish the range of post-common envelope binaries that are allowed energetically.

Among the initial results from Kepler were two striking lightcurves, for KOI 74 and KOI 81, in which the relative depths of the primary and secondary eclipses showed that the more compact, less luminous object was hotter than its stellar host. That result became particularly intriguing because a substellar mass had been derived for the secondary in KOI 74, which would make the high temperature challenging to explain; in KOI 81, the mass range for the companion was also reported to be consistent with a substellar object. We re-analyze the Kepler data and demonstrate that both companions are likely to be white dwarfs. We also find that the photometric data for KOI 74 show a modulation in brightness as the more luminous star orbits, due to Doppler boosting. The magnitude of the effect is sufficiently large that we can use it to infer a radial velocity amplitude accurate to 1 km/s. As far as we are aware, this is the first time a radial-velocity curve has been measured photometrically. Combining our velocity amplitude with the inclination and primary mass derived from the eclipses and primary spectral type, we infer a secondary mass of 0.22+/-0.03 Msun. We use our estimates to consider the likely evolutionary paths and mass-transfer episodes of these binary systems.

HD 49798 is a hydrogen depleted subdwarf O6 star and has an X-ray pulsating companion (RX J0648.0-4418). The X-ray pulsating companion is a massive white dwarf. Employing Eggleton's stellar evolution code with the optically thick wind assumption, we find that the hot subdwarf HD 49798 and its X-ray pulsating companion could produce a type Ia supernova (SN Ia) in future evolution. This implies that the binary system is a likely candidate of SN Ia progenitors. We also discussed the possibilities of some other WD + He star systems (e.g. V445 Pup and KPD 1930+2752) for producing SNe Ia.

We present high quality long slit spectra along the major and minor axes out to 1.5-2 Re (14-22 kpc) of three bright elliptical galaxies (NGC1600, NGC4125, NGC7619) obtained at the Hobby-Eberly Telescope (HET). We derive stellar kinematic profiles and Lick/IDS indices (Hbeta, Mgb, Fe5015, Fe5270, Fe5335, Fe5406). Moreover, for NGC4125 we derive gas kinematics and emission line strengths. We model the absorption line strengths using Simple Stellar Populations models that take into account the variation of [\alpha/Fe] and derive ages, total metallicity and element abundances. Overall, we find that the three galaxies have old and [\alpha/Fe] overabundant stellar populations with no significant gradients. The metallicity is supersolar at the center with a strong negative radial gradient. For NGC4125, several pieces of evidence point to a recent dissipational merger event. We calculate the broad band color profiles with the help of SSP models. All of the colors show sharp peaks at the center of the galaxies, mainly caused by the metallicity gradients, and agree well with the measured colors. Using the Schwarzschild's axisymmetric orbit superposition technique, we model the stellar kinematics to constrain the dark halos of the galaxies. We use the tight correlation between the Mgb strength and local escape velocity to set limits on the extent of the halos by testing different halo sizes. Logarithmic halos - cut at 60 kpc -minimize the overall scatter of the Mgb-Vesc relation. Larger cutoff radii are found if the dark matter density profile is decreasing more steeply at large radii.

Type Ia supernovae (SNe Ia) play an important role in the study of cosmic evolution, especially in cosmology. There are several progenitor models for SNe Ia proposed in the past years. In this paper, by considering the effect of accretion disk instability on the evolution of white dwarf (WD) binaries, we performed detailed binary evolution calculations for the WD + red-giant (RG) channel of SNe Ia, in which a carbon-oxygen WD accretes material from a RG star to increase its mass to the Chandrasekhar mass limit. According to these calculations, we mapped out the initial and final parameters for SNe Ia in the orbital period--secondary mass plane for various WD masses for this channel. We discussed the influence of the variation of the duty cycle value on the regions for producing SNe Ia. Similar to previous studies, this work also indicates that the long-period dwarf novae offer a possible ways for producing SNe Ia. Meanwhile, we find that the surviving companion stars from this channel have a low mass after SN explosion, which may provide a means for the formation of the population of single low-mass WDs (<0.45Msun).

Context: The nature of type Ia supernovae (SNe Ia) is still unclear. Metallicities may have an important effect on their properties. Aims: In this paper, we study the He star donor channel towards SNe Ia comprehensively and systematically at various metallicities. Methods: Employing Eggleton's stellar evolution code with the optically thick wind assumption, we calculated about 10 000 WD + He star systems and obtained SN Ia production regions of the He star donor channel with metallicities Z=0.03, 0.02, 0.004 and 0.0001. According to a detailed binary population synthesis approach, we also obtained SN Ia birthrates at various metallicities. Results: Our study shows that both the initial mass of the He donor star and the initial orbital period for SNe Ia increase with metallicity, while the minimum initial mass of the carbon--oxygen white dwarfs producing SNe Ia decreases with metallicity. For a constant star-formation galaxy, SN Ia birthrates increase with metallicity. If a single starburst is assumed, SNe Ia occur systemically earlier and the peak value of the birthrate is larger for a high Z, and the He star donor channel with different metallicities can produce the young SNe Ia with delay times $\sim$45-220Myr.

Recent investigations of the WD + MS channel of Type Ia supernovae (SNe Ia) imply that this channel may be the main contribution to the old population (>1Gyr) of SNe Ia. In the WD + MS channel, the WD could accrete material from a main-sequence or a slightly evolved star until it reaches the Chandrasekhar mass limit. The companions in this channel would survive after SN explosion and show distinguishing properties. In this Letter, based on SN Ia production regions of the WD + MS channel and three formation channels of WD + MS systems, we performed a detailed binary population synthesis study to obtain the properties of the surviving companions. The properties can be verified by future observations. We find that the surviving companions of the old SNe Ia have a low mass, which provides a possible way to explain the formation of the population of single low-mass WDs (<0.45Msun).

With more and more exoplanets being detected, it is paid closer attention to whether there are lives outside solar system. We try to obtain habitable zones and the probability distribution of terrestrial planets in habitable zones around host stars. Using Eggleton's code, we calculate the evolution of stars with masses less than 4.00 \mo. We also use the fitting formulae of stellar luminosity and radius, the boundary flux of habitable zones, the distribution of semimajor axis and mass of planets and the initial mass function of stars. We obtain the luminosity and radius of stars with masses from 0.08 to 4.00 \mo, and calculate the habitable zones of host stars, affected by stellar effective temperature. We achieve the probability distribution of terrestrial planets in habitable zones around host stars. We also calculate that the number of terrestrial planets in habitable zones of host stars is 45.5 billion, and the number of terrestrial planets in habitable zones around K type stars is the most, in the Milky Way.

Ultraviolet radiation is a double-edged sword to life. If it is too strong, the terrestrial biological systems will be damaged. And if it is too weak, the synthesis of many biochemical compounds can not go along. We try to obtain the continuous ultraviolet habitable zones, and compare the ultraviolet habitable zones with the habitable zones of host stars. Using the boundary ultraviolet radiation of ultraviolet habitable zone, we calculate the ultraviolet habitable zones of host stars with masses from 0.08 to 4.00 \mo. For the host stars with effective temperatures lower than 4,600 K, the ultraviolet habitable zones are closer than the habitable zones. For the host stars with effective temperatures higher than 7,137 K, the ultraviolet habitable zones are farther than the habitable zones. For hot subdwarf as a host star, the distance of the ultraviolet habitable zone is about ten times more than that of the habitable zone, which is not suitable for life existence.

A star will become brighter and brighter with stellar evolution, and the distance of its habitable zone will become farther and farther. Some planets outside the habitable zone of a host star during the main sequence phase may enter the habitable zone of the host star during other evolutionary phases. A terrestrial planet within the habitable zone of its host star is generally thought to be suited to life existence. Furthermore, a rocky moon around a giant planet may be also suited to life survive, provided that the planet-moon system is within the habitable zone of its host star. Using Eggleton's code and the boundary flux of habitable zone, we calculate the habitable zone of our Solar after the main sequence phase. It is found that Mars' orbit and Jupiter's orbit will enter the habitable zone of Solar during the subgiant branch phase and the red giant branch phase, respectively. And the orbit of Saturn will enter the habitable zone of Solar during the He-burning phase for about 137 million years. Life is unlikely at any time on Saturn, as it is a giant gaseous planet. However, Titan, the rocky moon of Saturn, may be suitable for biological evolution and become another Earth during that time. For low-mass stars, there are similar habitable zones during the He-burning phase as our Solar, because there are similar core masses and luminosities for these stars during that phase.

Type Ia supernovae (SNe Ia) play an important role in diverse areas of astrophysics, from the chemical evolution of galaxies to observational cosmology. However, the nature of the progenitors of SNe Ia is still unclear. In this paper, according to a detailed binary population synthesis study, we obtained SN Ia birthrates and delay times from different progenitor models, and compared them with observations. We find that the Galactic SN Ia birthrate from the double-degenerate (DD) model is close to those inferred from observations, while the birthrate from the single-degenerate (SD) model accounts for only about 1/2-2/3 of the observations. If a single starburst is assumed, the distribution of the delay times of SNe Ia from the SD model is a weak bimodality, where the WD + He channel contributes to the SNe Ia with delay times shorter than 100Myr, and the WD + MS and WD + RG channels to those with age longer than 1Gyr.

We selected brown dwarf candidates from the seventh Data Release of the Sloan Digital Sky Survey (SDSS DR7) with new photometric selectioncriteria based on a parameteriaztion of well-known L and T dwarfs. Then we confirmed their status with SDSS spectra. The candidates without SDSS spectra are cross matched in the Two Micron All Sky Survey (2MASS) and the Fourth Data Release of the UKIRT Infrared Deep Sky Survey (UKIDSS DR4). With the help of colors based on SDSS, 2MASS and UKIDSS, we are able to estimate spectral types of our candidates. We obtain reliable proper motions using positional and epoch information downloaded direct from the survey databases.

We identify 806 ultra-cool dwarfs from their SDSS riz photometry (of which 34 are newly discovered L dwarfs) and obtain proper motions through cross matching with UKIDSS and 2MASS. Proper motion and distance constraints show that nine of our ultra-cool dwarfs are members of widely separated binary systems; SDSS 0101 (K5V+M9.5V), SDSS 0207 (M1.5V+L3V), SDSS 0832 (K3III+L3.5V), SDSS 0858 (M4V+L0V), SDSS 0953 (M4V+M9.5V), SDSS 0956 (M2V+M9V), SDSS 1304 (M4.5V+L0V), SDSS 1631 (M5.5V+M8V), SDSS 1638 (M4V+L0V). One of these (SDSS 0832) is shown to be a companion to the bright K3 giant Eta Cancri. Such primaries can provide age and metallicity constraints for any companion objects, yielding excellent benchmark objects. Eta Cancri AB is the first wide ultra-cool dwarf + giant binary system identified. We present new observations and analysis that constrain the metallicity of Eta Cancri A to be near solar, and use recent evolutionary models to constrain the age of the giant to be 2.2-6.1 Gyr. If Eta Cancri B is a single object, we estimate its physical attributes to be; mass = 63-82 M_Jup, T_eff = 1800+/-150 K, log g = 5.3-5.5, [M/H] = 0.0+/-0.1. Its colours are non typical when compared to other ultra-cool dwarfs, and we also assess the possibility that Eta Cancri B is itself an unresolved binary, showing that the combined light of an L4 + T4 system could provide a reasonable explanation for its colours.

Employing Eggleton's stellar evolution code with an optically thick wind assumption, we have systematically studied the WD + He star channel of Type Ia supernovae (SNe Ia), in which a carbon-oxygen WD accretes material from a He main-sequence star or a He subgiant to increase its mass to the Chandrasekhar mass. We mapped out the parameter spaces for producing SNe Ia. According to a detailed binary population synthesis approach, we find that the Galactic SN Ia birthrate from this channel is $\sim$$0.3\times 10^{-3} {\rm yr}^{-1}$, and that this channel can produce SNe Ia with short delay times ($\sim45-$140 Myr). We also find that the surviving companion stars in this channel have a high spatial velocity ($>$400 km/s) after SN explosion, which could be an alternative origin for hypervelocity stars (HVSs), especially for HVSs such as US 708.

{Context} Recent investigations of the white dwarf (WD) + He star channel of Type Ia supernovae (SNe Ia) imply that this channel can produce SNe Ia with short delay times. The companion stars in this channel would survive and be potentially identifiable. {Aims} In this Letter, we study the properties of the companion stars of this channel at the moment of SN explosion, which can be verified by future observations. {Methods} According to SN Ia production regions of the WD + He star channel and three formation channels of WD + He star systems, we performed a detailed binary population synthesis study to obtain the properties of the surviving companions. {Results} We obtained the distributions of many properties of the companion stars of this channel at the moment of SN explosion. We find that the surviving companion stars have a high spatial velocity (>400 km/s) after SN explosion, which could be an alternative origin for hypervelocity stars (HVSs), especially for HVSs such as US 708.

Most stars are members of binaries, and the evolution of a star in a close binary system differs from that of an ioslated star due to the proximity of its companion star. The components in a binary system interact in many ways and binary evolution leads to the formation of many peculiar stars, including blue stragglers and hot subdwarfs. We will discuss binary evolution and the formation of blue stragglers and hot subdwarfs, and show that those hot objects are important in the study of evolutionary population synthesis (EPS), and conclude that binary interactions should be included in the study of EPS. Indeed, binary interactions make a stellar population younger (hotter), and the far-ultraviolet (UV) excess in elliptical galaxies is shown to be most likely resulted from binary interactions. This has major implications for understanding the evolution of the far-UV excess and elliptical galaxies in general. In particular, it implies that the far-UV excess is not a sign of age, as had been postulated prviously and predicts that it should not be strongly dependent on the metallicity of the population, but exists universally from dwarf ellipticals to giant ellipticals.

Blue stragglers have been found in all populations. These objects are important in both stellar evolution and stellar population synthesis. Much evidence shows that blue stragglers are relevant to primordial binaries. Here, we summarize the links of binary evolution and blue stragglers, describe the characteristics of blue stragglers from different binary evolutionary channels, and show their consequences for binary population synthesis, such as for the integrated spectral energy distribution, the colour-magnitude diagram, the specific frequency, and the influences on colours etc..

The nature of the progenitors of Type Ia supernovae (SNe Ia) is still unclear. In this paper, by considering the effect of the instability of accretion disk on the evolution of white dwarf (WD) binaries, we performed binary evolution calculations for about 2400 close WD binaries, in which a carbon--oxygen WD accretes material from a main-sequence star or a slightly evolved subgiant star (WD + MS channel), or a red-giant star (WD + RG channel) to increase its mass to the Chandrasekhar (Ch) mass limit. According to these calculations, we mapped out the initial parameters for SNe Ia in the orbital period--secondary mass ($\log P^{\rm i}-M^{\rm i}_2$) plane for various WD masses for these two channels, respectively. We confirm that WDs in the WD + MS channel with a mass as low as $0.61 M_\odot$ can accrete efficiently and reach the Ch limit, while the lowest WD mass for the WD + RG channel is $1.0 \rm M_\odot$. We have implemented these results in a binary population synthesis study to obtain the SN Ia birthrates and the evolution of SN Ia birthrates with time for both a constant star formation rate and a single starburst. We find that the Galactic SN Ia birthrate from the WD + MS channel is $\sim$$1.8\times 10^{-3} {\rm yr}^{-1}$ according to our standard model, which is higher than previous results. However, similar to previous studies, the birthrate from the WD + RG channel is still low ($\sim$$3\times 10^{-5} {\rm yr}^{-1}$). We also find that about one third of SNe Ia from the WD + MS channel and all SNe Ia from the WD + RG channel can contribute to the old populations ($\ga$1 Gyr) of SN Ia progenitors.

Masses are a fundamental parameter, but they are not well known for most hot subdwarfs. In general, the mass of a hot subdwarf is derived with asteroseismology or dynamical methods, for which it is often difficult to obtain the necessary data from observations. We intend to find an approach to deriving the masses of hot subdwarfs from observational data in the literature. We presented full evolutionary calculations for hot subdwarfs in a wide mass range (0.33 $M_\odot$ to 1.4 $M_\odot$) for a Population I metallicity of $Z$=0.02, and obtained a relation between $M_{\rm p}$ and $\log (\frac{T_{\rm eff}^4}{g})$, where $M_{\rm p}$, $T_{\rm eff}$, and $g$ are the most probable mass, effective temperature, and gravity. This relation is used to study the masses of some observed hot subdwarfs. We proposed a method of determining the masses of hot subdwarfs. Using this method, we studied the masses of hot subdwarfs from the ESO supernova Ia progenitor survey and Hamburg quasar survey. The study shows that most of subdwarf B stars have masses between 0.42 and 0.54 $M_\odot$, whilst most sdO stars are in the range 0.40 $\sim$ 0.55 $M_\odot$. Comparing our study to the theoretical mass distributions of Han et al. (2003), we found that sdO stars with mass less than $\sim$ 0.5 $M_\odot$ may evolve from sdB stars, whilst most high-mass($>$ 0.5 $M_\odot$) sdO stars result from mergers directly.

Single Degenerate model is the most widely accepted progenitor model of type Ia supernovae (SNe Ia), in which a carbon-oxygen white dwarf (CO WD) accretes hydrogen-rich material from a main sequence or a slightly evolved star (WD +MS) to increase its mass, and explodes when its mass approaches the Chandrasekhar mass limit. During the mass transfer phase between the two components, an optically thick wind may occur and the material lost as the wind may exist as circumstellar material (CSM). Searching the CSM around progenitor star is helpful to discriminate different progenitor models of SNe Ia. Meanwhile, the CSM is a source of color excess.The purpose of this paper is to study the color excess produced from the single-degenerate progenitor model with optically thick wind, and reproduce the distribution of color excesses of SNe Ia. Meng et al. (2009) systemically carried out binary evolution calculation of the WD +MS systems for various metallicities and showed the parameters of the systems before Roche lobe overflow and at the moment of supernova explosion in Meng & Yang (2009). With the results of Meng et al. (2009), we calculate the color excesses of SNe Ia at maximum light via a simple analytic method.We reproduces the distribution of color excesses of SNe Ia by our binary population synthesis approach if the velocity of the optically thick wind is taken to be of order of magnitude of 10 km s$^{\rm -1}$. However, if the wind velocity is larger than 100 km s$^{\rm -1}$, the reproduction is bad.

The single-degenerate model is currently a favourable progenitor model for Type Ia supernovae (SNe Ia). Recent investigations on the WD + He star channel of the single-degenerate model imply that this channel is noteworthy for producing SNe Ia. In this paper we studied SN Ia birthrates and delay times of this channel via a detailed binary population synthesis approach. We found that the Galactic SN Ia birthrate from the WD + He star channel is $\sim 0.3\times 10^{-3} {\rm yr}^{-1}$ according to our standard model, and that this channel can explain SNe Ia with short delay times ($\sim4.5\times10^7 - 1.4\times10^8$ yr). Meanwhile, these WD + He star systems may be related to the young supersoft X-ray sources prior to SN Ia explosions.

The energy transfer of W UMa contact binaries and its effects on the secondary in W UMa contact binaries are investigated. Relations are given between the mass ratio (q) and the relative energy transfer rates, i.e. U1, the ratio of the transferred luminosity to the surface luminosity of the primary and U2, the ratio of the transferred luminosity to the nuclear luminosity of the secondary. The theoretical curves(U1 vs q and U2 vs q) are derived. Although these curves can reflect the distribution of U1 vs q and U2 vs q, some observational systems are significantly deviated from these curves. It is mainly resulted from the difference in the effective temperatures of the components in W UMa systems.The radius and the density of the secondary are related to the relative energy tran-sfer rate U2: the higher is U2, the greater is the expansion and the lower is the density of the secondaries in W UMa systems. The temperature difference of W UMa binary components is correlated with the relative energy transfer rate U1 and decreases with increasing U1. This might suggest that there is a thermal coupling between two components in W UMa contact binaries, and that the classif-ication of W UMa contact binaries into A- or W-types depends on the energy tran-sfer from the primary to the secondary. The temperature difference of W UMa bin-ary components is poorly correlated with the mass of the primary. This suggests that the properties of the common envelope of W UMa contact binaries might not have a significant effect on the energy transfer between two components.

In this paper, we systematically studied blue stragglers produced from primordial binary evolution via a binary population synthesis approach, and examined their contribution to the integrated spectral energy distributions of the host clusters. The mass transfer efficiency, $\beta$, is an important parameter for the final products (then blue stragglers) after mass transfer, and it is set to be 0.5 except for case A binary evolution. The study shows that primordial binary evolution may produce blue stragglers at any given times and that different evolutionary channels are corresponding for blue stragglers in different visual magnitude regions (in V band) on the colour-magnitude diagram (CMD) of clusters. The specific frequency of blue stragglers obtained from primordial binary evolution decreases with time first, and then increases again when the age is larger than 10Gyr, while that from angular momentum loss induced by magnetic braking in low-mass binaries increases with time and exceeds that of primordial binary evolution in a population older than 3 Gyr. Meanwhile, blue stragglers resulting from primordial binary evolution are dominant contributors to the ISEDs in ultraviolet and blue bands in a population between 0.3 and 2.0 Gyr. The value of $\beta$ significantly affects on the final results. For old open clusters, the assumption of $\beta =1$ when the primary is in HG at the onset of mass transfer matches the observations better than that of $\beta =0.5$ from the locations of BSs on the CMDs. Our study also shows that, for most Galactic open clusters, the specific frequency of blue stragglers obtained from our simulations is lower than that of observations, which is puzzling.

Aims. We try to identify ultra-cool dwarfs from the seventh Data Release of the Sloan Digital Sky Survey (SDSS DR7) with SDSS i-z and r-z colors. We also obtain proper motion data from SDSS, 2MASS, and UKIDSS and improve spectral typing from SDSS and 2MASS photometric colors. Methods. We selected ultra-cool dwarf candidates from the SDSS DR7 with new photometric selection criteria, which are based on a parameterization study of known L and T dwarfs. The objects are then cross-identified with the Two Micron All Sky Survey and the Fourth Data Release of the UKIRT Infrared Deep Sky Survey (UKIDSS DR4). We derive proper motion constraints by combining SDSS, 2MASS, and UKIDSS positional information. In this way we are able to assess, to some extent, the credence of our sample using a multi epoch approach, which complements spectroscopic confirmation. Some of the proper motions are affected by short baselines, but, as a general tool, this method offers great potential to confirm faint L dwarfs as UKIDSS coverage increases. In addition we derive updated color-spectral type relations for L and T dwarfs with SDSS and 2MASS magnitudes. Results. We present 59 new nearby M and L dwarfs selected from the imaging catalog of the SDSS DR7, including proper motions and spectral types calculated from the updated color-spectral type relations. and obtain proper motions from SDSS, 2MASS, and UKIDSS for all of our objects.

The single-degenerate channel for the progenitors of type Ia supernovae (SNe Ia) are currently accepted, in which a carbon-oxygen white dwarf (CO WD) accretes hydrogen-rich material from its companion, increases its mass to the Chandrasekhar mass limit, and then explodes as a SN Ia. Incorporating the prescription of \citet{HAC99a} for the accretion efficiency into Eggleton's stellar evolution code and assuming that the prescription is valid for \emph{all} metallicities, we performed binary stellar evolution calculations for more than 25,000 close WD binaries with metallicities $Z=0.06, 0.05, 0.04, 0.03, 0.02, 0.01, 0.004, 0.001, 0.0003$ and 0.0001. For our calculations, the companions are assumed to be unevolved or slightly evolved stars (WD + MS). As a result, the initial parameter spaces for SNe Ia at various $Z$ are presented in orbital period-secondary mass ($\log P_{\rm i}, M_{\rm 2}^{\rm i}$) planes. Our study shows that both the initial mass of the secondary and the initial orbital period increase with metallicity. Thus, the minimum mass of the CO WD for SNe Ia decreases with metallicity $Z$. The difference of the minimum mass may be as large as 0.24 $M_{\odot}$ for different Z. Adopting the results above, we studied the birth rate of SNe Ia for various $Z$ via a binary population synthesis approach. If a single starburst is assumed, SNe Ia occur systemically earlier and the peak value of the birth rate is larger for a high $Z$. The Galactic birth rate from the WD + MS channel is lower than (but comparable to) that inferred from observations. Our study indicates that supernovae like SN2002ic would not occur in extremely low-metallicity environments, if the delayed dynamical-instability model in \citet{HAN06} is appropriate.

Type Ia supernovae (SNe Ia) play an important role in astrophysics, especially in the study of cosmic evolution. There are several progenitor models for SNe Ia proposed in the past years. In this paper, we have carried out a detailed study of the He star donor channel, in which a carbon-oxygen white dwarf (CO WD) accretes material from a He main sequence star or a He subgiant to increase its mass to the Chandrasekhar mass. Employing Eggleton's stellar evolution code with an optically thick wind assumption, and adopting the prescription of Kato & Hachisu (2004) for the mass accumulation efficiency of the He-shell flashes onto the WDs, we performed binary evolution calculations for about 2600 close WD binary systems. According to these calculations, we mapped out the initial parameters for SNe Ia in the orbital period--secondary mass ($\log P^{\rm i}-M^{\rm i}_2$) plane for various WD masses from this channel. The study shows that the He star donor channel is noteworthy for producing SNe Ia (i.e. $\sim 1.2\times10^{-3} {\rm yr}^{-1}$ in the Galaxy), and that the progenitors from this channel may appear as supersoft X-ray sources. Importantly, this channel can explain SNe Ia with short delay times ($\la 10^{8}$ yr), which is consistent with recent observational implications of young populations of SN Ia progenitors.

There is still considerable debate over the progenitors of type Ia supernovae (SNe Ia). Likewise, it is not agreed how single white dwarfs with masses less than ~0.5 Msun can be formed in the field, even though they are known to exist. We consider whether single low-mass white dwarfs (LMWDs) could have been formed in binary systems where their companions have exploded as a SN Ia. In this model, the observed single LMWDs are the remnants of giant-branch donor stars whose envelopes have been stripped off by the supernova explosion. We investigate the likely remnants of SNe Ia, including the effects of the explosion on the envelope of the donor star. We also use evolutionary arguments to examine alternative formation channels for single LMWDs. In addition, we calculate the expected kinematics of the potential remnants of SNe Ia. SN Ia in systems with giant-branch donor stars can naturally explain the production of single LMWDs. It seems difficult for any other formation mechanism to account for the observations, especially for those single LMWDs with masses less than ~0.4 Msun. Independent of those results, we find that the kinematics of one potentially useful population containing single LMWDs is consistent with our model. Studying remnant white-dwarf kinematics seems to be a promising way to investigate SN Ia progenitors. The existence of single LMWDs appears to constitute evidence for the production of SNe Ia in binary systems with a red-giant donor star. Other single white dwarfs with higher space velocities support a second, probably dominant, population of SN Ia progenitors which contained main-sequence or subgiant donor stars at the time of explosion. The runaway stars LP400-22 and US 708 suggest the possibility of a third formation channnel for some SNe Ia in systems where the donor stars are hot subdwarfs.

The excess of far-ultraviolet (far-UV) radiation in elliptical galaxies has remained one of their most enduring puzzles. In contrast, the origin of old blue stars in the Milky Way, hot subdwarfs, is now reasonably well understood: they are hot stars that have lost their hydrogen envelopes by various binary interactions. Here, we review the main evolutionary channels that produce hot subdwarfs in the Galaxy and present the results of binary population synthesis simulations that reproduce the main properties of the Galactic hot-subdwarf population. Applying the same model to elliptical galaxies, we show how this model can explain the main observational properties of the far-UV excess, including the far-UV spectrum, without the need to invoke ad hoc physical processes. The model implies that the UV excess is not a sign of age, as has been postulated previously, and predicts that it should not be strongly dependent on the metallicity of the population.

Well-determined physical parameters of 130 W UMa systems have been collected from the literature. Based on these data, the evolutionary status and dynamical evolution of W UMa systems are investigated. It is found that there is no evolutionary difference between W- and A-type systems in $M-J$ diagram which is consistent with the results derived from the analysis of observed spectral type, $M-R$ and $M-L$ diagrams of W UMa systems. $M-R$ and $M-L$ diagrams of W- and A-type systems indicate that a large amount of energy should be transferred from the more massive to the less massive component so that they are not in thermal equilibrium and undergo thermal relaxation oscillation (TRO). Meanwhile, the distribution of angular momentum, together with the distribution of mass ratio, suggests that the mass ratio of the observed W UMa systems is decreased with the decrease of their total mass. This could be the result of the dynamical evolution of W UMa systems which suffer angular momentum loss (AML) and mass loss due to magnetic stellar wind (MSW). Consequently, the tidal instability forces these systems towards the lower q values and finally to fast rotating single stars.

We have carried out a detailed study of the single-degenerate channel for the progenitors of type Ia supernovae (SNe Ia). In the model, a carbon-oxygen white dwarf (CO WD) accretes hydrogen-rich material from an unevolved or a slightly evolved non-degenerate companion to increase its mass to Chandrasekhar mass limit. Incorporating the prescription of \cite{HAC99a} for the accretion efficiency into Eggleton's stellar evolution code and assuming that the prescription is valid for all metallicities, we performed binary stellar evolution calculations for more than 25,000 close WD binary systems with various metallicities. The initial parameter spaces for SNe Ia are presented in an orbital period-secondary mass ($\log P_{\rm i}, M_{\rm 2}^{\rm i}$) plane for each $Z$. Adopting the results above, we studied the birth rate of SNe Ia for various $Z$ via binary population synthesis. From the study, we see that for a high $Z$, SNe Ia occur systemically earlier and the peak value of the birth rate is larger if a single starburst is assumed. The Galactic birth rate from the channel is lower than (but comparable to) that inferred from observations. We also showed the distributions of the parameters of the binary systems at the moment of supernova explosion and the distributions of the properties of companions after supernova explosion. The former provides physics input to simulate the interaction between supernova ejecta and its companion, and the latter is helpful to search for the companions in supernova remnants.

In this talk, we present the general principles of binary evolution and give two examples. The first example is the formation of subdwarf B stars (sdBs) and their application to the long-standing problem of ultraviolet excess (also known as UV-upturn) in elliptical galaxies. The second is for the progenitors of type Ia supernovae (SNe Ia). We discuss the main binary interactions, i.e., stable Roche lobe overflow (RLOF) and common envelope (CE) evolution, and show evolutionary channels leading to the formation of various binary-related objects. In the first example, we show that the binary model of sdB stars of Han et al. (2002, 2003) can reproduce field sdB stars and their counterparts, extreme horizontal branch (EHB) stars, in globular clusters. By applying the binary model to the study of evolutionary population synthesis, we have obtained an ``a priori'' model for the UV-upturn of elliptical galaxies and showed that the UV-upturn is most likely resulted from binary interactions. This has major implications for understanding the evolution of the UV excess and elliptical galaxies in general. In the second example, we introduce the single degenerate channel and the double degenerate channel for the progenitors of SNe Ia. We give the birth rates and delay time distributions for each channel and the distributions of companion stars at the moment of SN explosion for the single degenerate channel, which would help to search for the remnant companion stars observationally.

Single-star stellar population (ssSSP) models are usually used for stellar population studies. However, more than 50% of stars are in binaries and evolve differently from single stars. This suggests that the effects of binary interactions should be considered when modeling the stellar populations of galaxies and star clusters. Via a rapid spectral stellar population synthesis (RPS) model, we give detailed studies of the effects of binary interactions on the Lick indices and colours of stellar populations, and on the determination of the stellar ages and metallicities of populations. The results show that binary interactions make stellar populations less luminous, bluer, with larger age-sensitive Lick index (Hbeta) and less metallicity-sensitive indices (e.g., Mgb, Fe5270 and Fe5335) compared to ssSSPs. It also shows that when ssSSP models are used to determine the ages and metallicities of stellar populations, lower ages or metallicities will be obtained, taking two line indices (Hbeta and [MgFe]) and two colours (e.g., u-R and R-K), respectively. Some relations for linking the stellar-population parameters obtained by ssSSPs to those obtained by binary-star stellar populations (bsSSPs) are presented in the work. This can help us to get some absolute values for stellar-population parameters and is useful for absolute studies. However, it is found that the relative luminosity-weighted stellar ages and metallicities obtained via ssSSPs and bsSSPs are similar. This suggests that ssSSPs can be used for most stellar population studies, except in some special cases.

Following Paczy\'{n}ski & Zi\'{o}lkowski (1968) and Han et al. (1994), we assume that the envelope of an asymptotic giant branch (AGB) or a first giant branch (FGB) star is lost when the binding energy of the envelope is equal to zero ($\Delta W=0$) and the core mass of the AGB star or the FGB star at the point ($\Delta W=0$) is taken as the final mass. Using this assumption, we calculate the IFMRs for stars of different metallicities.We find that the IFMRs depends strongly on the metallicity, i.e. $Z=0.0001, 0.0003, 0.001, 0.004, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.08$ and 0.1. From $Z=0.04$, the final mass of the stars with a given initial mass increases with increasing or decreasing metallicity. The difference of the final mass due to the metallicity may be up to 0.4 $M_{\odot}$. A linear fit of the initial-final mass relationship in NGC 2099 (M37) shows a potential evidence of the effect of metallicity on the IFMR. The IFMR for stars of $Z=0.02$ obtained in the paper matches well with those inferred observationally in the Galaxy. For $Z\geq 0.02$, helium WDs are obtained from the stars of $M_{\rm i}\leq 1.0 M_{\odot}$ and this result is upheld by the discovery of numerous low-mass WDs in NGC 6791 which is a metal-rich old open cluster. Using the IFMR for stars of $Z=0.02$ obtained in the paper, we have reproduced the mass distribution of DA WDs in Sloan DR4 except for some ultra-massive white dwarfs. The trend that the mean mass of WDs decreases with effective temperature may originate from the increase of the initial metallicities of stars. We predict that metal-rich low-mass stars may become under-massive white dwarfs.

We have carried out a study of the chemical abundances of $^1$H, $^4$He, $^{12}$C, $^{13}$C, $^{14}$N, $^{15}$N, $^{16}$O, $^{17}$O, $^{20}$Ne and $^{22}$Ne in symbiotic stars (SSs) by means of a population synthesis code. We find that the ratios of the number of O-rich SSs to that of C-rich SSs in our simulations are between 3.4 and 24.1, depending on the third dredge-up efficiency $\lambda$ and the terminal velocity of the stellar wind $v(\infty)$. The fraction of SSs with $extrinsic$ C-rich cool giants in C-rich cool giants ranges from 2.1% to 22.7%, depending on $\lambda$, the common envelope algorithm and the mass-loss rate. Compared with the observations, the distributions of the relative abundances of $^{12}$C/$^{13}$C vs. [C/H] of the cool giants in SSs suggest that the thermohaline mixing in low-mass stars may exist. The distributions of the relative abundances of C/N vs. O/N, Ne/O vs. N/O and He/H vs. N/O in the symbiotic nebulae indicate that it is quite common that the nebular chemical abundances in SSs are modified by the ejected materials from the hot components. Helium overabundance in some symbiotic nebulae may be relevant to a helium layer on the surfaces of white dwarf accretors.

The binary fraction among extreme horizontal branch (EHB) stars in Galactic globular clusters (GCs) is an order of magnitude lower than the binary fraction among their counterparts, field hot subdwarfs. This casts serious doubt on their formation channels. In this {\em Letter}, I explain the difference between the field and the cluster EHB stars with the binary model of Han et al. (2002, 2003) for the formation of EHB stars. With the binary population synthesis code of Han et al. (2002, 2003), I follow the evolution of simple stellar populations resulting from single star bursts (note that Han et al. 2002, 2003, adopted a constant star formation rate over the last 15 Gyr for the production of field EHB stars), and obtain EHB stars at different stellar population ages. I found that the binary fraction among EHB stars decreases with the stellar population age. The fraction of EHB binaries with orbital periods $P<5 {\rm d}$ is $\sim 2.5%$ for a stellar population of 10 Gyr from the standard simulation set. The binary model of Han et al. (2002, 2003) is able to explain the lack of EHB binaries in globular clusters. I also propose that the precise determination of the physical parameters of close EHB binaries in GCs can lead to the strictest constraint on common-envelope ejection efficiency.

Binary population synthesis shows that mass transfer from a giant star to a main-sequence (MS) companion may account for some observed long-orbital period blue stragglers. However, little attention {\bf is paid to this blue straggler formation scenario} as dynamical instability often happens when the mass donor is a giant star. In this paper, we have studied the critical mass ratio, $q_{\rm c}$, for dynamically stable mass transfer from a giant star to a MS companion using detailed evolution calculations. The results show that a more evolved star is generally less stable for Roche lobe overflow. Meanwhile, $q_{\rm c}$ almost linearly increases with the amount of the mass and angular momentum {\bf lost} during mass transfer, but has little dependance on stellar wind. To conveniently use the result, we give a fit of $q_{\rm c}$ as a function of the stellar radius at the onset of Roche lobe overflow and of the mass transfer efficiency during the Roche lobe overflow. To examine the formation of blue stragglers from mass transfer between giants and MS stars, we have performed Monte Carlo simulations with various $q_{\rm c}$. {\bf The simulations show that some binaries with the mass donor on the first giant branch may contribute to blue stragglers with $q_{\rm c}$ obtained in this paper but will not from previous $q_{\rm c}$. Meanwhile, from our $q_{\rm c}$, blue stragglers from the mass transfer between an AGB star and a MS companion may be more numerous and have a wider range of orbital periods than those from the other $q_{\rm c}$.

The flux excess of elliptical galaxies in the far-ultraviolet can be reproduced by population synthesis models when accounting for the population of old hot helium-burning subdwarf stars. This has been achieved by Han and coworkers through a quantitative model of binary stellar evolution. Here, we compare the resulting evolutionary population synthesis model to the GALEX far-near ultraviolet colors (FUV-NUV) of Virgo cluster early-type galaxies that were published by Boselli and coworkers. FUV-NUV is reddest at about the dividing luminosity of dwarf and giant galaxies, and becomes increasingly blue for both brighter and fainter luminosities. This behavior can be easily explained by the binary model with a continuous sequence of longer duration and later truncation of star formation at lower galaxy masses. Thus, in contrast to previous conclusions, the GALEX data do not require a dichotomy between the stellar population properties of dwarfs and giants. Their apparently opposite behavior in FUV-NUV occurs naturally when the formation of hot subdwarfs through binary evolution is taken into account.

The WD+MS channel of the single-degenerate scenario is currently favourable for progenitors of type Ia supernovae (SNe Ia). Incorporating the results of detailed binary evolution calculations for this channel into the latest version of a binary population synthesis code, I obtained the distributions of many properties of the companion stars at the moment of SN explosion. The properties can be verified by future observations.

We first presented an isochrone database that can be widely used for stellar population synthesis studies and colour-magnitude diagram (CMD) fitting. The database consists of the isochrones of both single star and binary star simple stellar populations (ss-SSPs and bs-SSPs). The ranges for the age and metallicity of populations are 0--15 Gyr and 0.0001--0.03, respectively. All data are available for populations with Salpeter IMF and Chabrier IMF. Then based on the isochrone database, we built a rapid stellar population synthesis (RPS) model and calculated the high-resolution (0.3 AA) integrated spectral energy distributions (SEDs), Lick indices and colour indices for bs-SSPs and ss-SSPs. In particular, we calculated the UBVRIJHKLM colours, ugriz colours and some composite colours that consist of magnitudes on different systems. As an example for applying the isochrone database for CMD fitting, we fitted the CMDs of star clusters M67 and NGC1868 and obtained their distance moduli, colour excesses, stellar metallicities and ages. The results showed that the isochrones of bs-SSPs are closer to those of real star clusters. The isochrone database, colours, Lick indices of populations can be downloaded via http://www.ynao.ac.cn/~bps/zhongmu/download.htm

The amount of $^{56}$Ni produced in type Ia supernova (SN Ia) explosion is probably the most important physical parameter underlying the observed correlation of SN Ia luminosities with their light curves. Based on an empirical relation between the $^{56}$Ni mass and the light curve parameter $\triangle m_{15}$, we obtained rough estimates of the $^{56}$Ni mass for a large sample of nearby SNe Ia with the aim of exploring the diversity in SN Ia. We found that the derived $^{56}$Ni masses for different SNe Ia could vary by a factor of ten (e.g., $M_{\rm Ni}=0.1 - 1.3$ $M_{\odot}$), which cannot be explained in terms of the standard Chandrasekhar-mass model (with a $^{56}$Ni mass production of 0.4 -- 0.8 $M_{\odot}$). Different explosion and/or progenitor models are clearly required for various SNe Ia, in particular, for those extremely nickel-poor and nickel-rich producers. The nickel-rich (with $M_{\rm Ni}$ $>$ 0.8 $M_{\odot}$) SNe Ia are very luminous and may have massive progenitors exceeding the Chandrasekhar-mass limit since extra progenitor fuel is required to produce more $^{56}$Ni to power the light curve. This is also consistent with the finding that the intrinsically bright SNe Ia prefer to occur in stellar environments of young and massive stars. For example, 75% SNe Ia in spirals have $\Delta m_{15} < 1.2$ while this ratio is only 18% in E/S0 galaxies. On the other hand, the nickel-poor SNe Ia (with $M_{\rm Ni}$ $<$ 0.2 $M_{\odot}$) may invoke the sub-Chandrasekhar model, as most of them were found in early-type E/S0 galaxies dominated by the older and low-mass stellar populations. This indicates that SNe Ia in spiral and E/S0 galaxies have progenitors of different properties.

More than about 50% stars of galaxies are in binaries, but most stellar population studies take single star-stellar population (ss-SSP) models, which do not take binary interactions into account. In fact, the integrated peculiarities of ss-SSPs are various from those of stellar populations with binary interactions (bs-SSPs). Therefore, it is necessary to investigate the effects of binary interactions on the Lick indices and colours of populations detailedly. We show some formulae for calculating the difference between the Lick indices and colours of bs-SSPs, and those of ss-SSPs. Twenty-five Lick indices and 12 colours are studied in the work. The results can be conveniently used for estimating the effects of binary interactions on stellar population studies and for adding the effects of binary interactions into existing ss-SSP models. The results and a few procedures can be obtained on request to the authors or via http://www.ynao.ac.cn/~bps/zhongmu/download.htm

We describe a kinematic method which is capable of determining the overall mass scale in SUSY-like events at a hadron collider with two missing (dark matter) particles. We focus on the kinematic topology in which a pair of identical particles is produced with each decaying to two leptons and an invisible particle (schematically, $pp\to YY+jets$ followed by each $Y$ decaying via $Y\to \ell X\to \ell\ell'N$ where $N$ is invisible). This topology arises in many SUSY processes such as squark and gluino production and decay, not to mention $t\anti t$ di-lepton decays. In the example where the final state leptons are all muons, our errors on the masses of the particles $Y$, $X$ and $N$ in the decay chain range from 4 GeV for 2000 events after cuts to 13 GeV for 400 events after cuts. Errors for mass differences are much smaller. Our ability to determine masses comes from considering all the kinematic information in the event, including the missing momentum, in conjunction with the quadratic constraints that arise from the $Y$, $X$ and $N$ mass-shell conditions. Realistic missing momentum and lepton momenta uncertainties are included in the analysis.

The presence of blue straggler stars (BSs) in star clusters has proven a challenge to conventional simple stellar population (SSP) models. Conventional SSP models are based on the evolution theory of single stars. Meanwhile, the typical locations of BSs in the colour-magnitude diagram of a cluster are brighter and bluer than the main sequence turn-off point. Such loci cannot be predicted by single-star evolution theory. However, stars with such properties contribute significantly to the integrated light of the cluster. In this paper, we reconstruct the integrated properties of the Large Magellanic Cloud cluster ESO 121-SC03, based on a detailed exploration of the individual cluster stars, and with particular emphasis on the cluster's BSs. We find that the integrated light properties of ESO 121-SC03 are dramatically modified by its BS component. The integrated spectral energy distribution (ISED) flux level is significantly enhanced toward shorter wavelengths, and all broad-band colours become bluer. When fitting the fully integrated ISED of this cluster based on conventional SSP models, the best-fitting values of age and metallicity are significantly underestimated compared to the true cluster parameters. The age underestimate is $\sim40$ per cent if we only include the BSs within the cluster's half-light radius and $\sim60$ per cent if all BSs are included. The corresponding underestimates of the cluster's metallicity are $\sim30$ and $\sim60$ per cent, respectively. The populous star clusters in the Magellanic Clouds are ideal objects to explore the potential importance of BSs for the integrated light properties of more distant unresolved star clusters in a statistically robust manner, since they cover a large range in age and metallicity.

Using a widely used stellar population synthesis model, we study the ability of using pairs of AB system colours to break the well-known stellar age--metallicity degeneracy and give constraints on the two stellar-population parameters (age and metallicity). The relative age and metallicity sensitivities of AB system colours that relate to u, B, g, V, r, R, i, I, z, J, H, and K bands are presented, and the abilities of various colour pairs for breaking the age--metallicity degeneracy are quantified by an uncertainty parameter (UP) method. Our results suggest that a few pairs of colours can be used to constrain the two above stellar-population parameters. This will be very useful for exploring the stellar populations of distant galaxies. In detail, colour pairs [(r-K), (u-R)] and [(r-K), (u-r)] are shown to be the best pairs for estimating stellar ages and metallicities. They can constrain two stellar-population parameters on average with age uncertainties less than 3.89 Gyr and metallicity uncertainties less than 0.34 dex for typical uncertainties in colours. Some other colour pairs, such as [(R-K), (u-R)], [(I-K), (u-R)], [(R-K), (u-r)] and [(i-J), (u-R)], can possibly be used for constraining the two parameters, too. As a whole, our results suggest that colours relating to both UBVRIJHK and ugriz magnitudes are much better than either UBVRIJHK colours or ugriz colours for breaking the well-known degeneracy. The results also show that the stellar ages and metallicities of galaxies observed by the Sloan Digital Sky Survey (SDSS) and the Two-Micron All-Sky Survey (2MASS) can be estimated via photometry data. It is also shown that the colours can be used in conjunction with line indices to measure stellar-population parameters.

We constructed some main-sequence mergers from case A binary evolution and studied their characteristics via Eggleton's stellar evolution code. Both total mass and orbital angular momentum are conservative in our binary evolutions. Some mergers might be on the left of the ZAMS as defined by normal surface composition on a CMD because of enhanced surface helium content. The study also shows that central hydrogen content of the mergers is independent of mass. As a consequence, we fit the formula of magnitude and B-V of the mergers when they return back to thermal equilibrium with maximum error 0.29 and 0.037, respectively. Employing the consequences above, we performed Monte Carlo simulations to examine our models in NGC 2682 and NGC 2660. In NGC 2682, binary mergers from our models cover the region with high luminosity, but its importance is much less than that of AML. Our results are well-matched to the observations of NGC2660 if there is about 0.5Mo of mass loss in the merger process.

Some studies show that there are recent star formations in galaxies including early type galaxies. This suggests that we possibly did not obtain the accurate ages and metallicities of the dominating stellar populations of galaxies from colors. We study how young populations (YSPs) in composite stellar populations affect the colors of star systems and analyze how the stellar ages and metallicities derived from colors possibly deviate from those of the dominating populations of composite star systems. Some possible distributions of the differences between stellar ages and metallicities determined by colors and those of the dominating stellar populations of composite stellar systems are presented in this work. Furthermore, some possible distributions of the differences between colors of composite stellar populations and those of their dominating stellar populations are also shown in this paper. In addition, we compare the stellar ages and metallicities determined respectively by a color pair and two line strength indices. The results obtained by a color pair and those obtained by a few independent colors are also compared in the paper. We also discuss about how the YSPs may affect the Fundamental Plane and Kormendy relation of early type galaxies. The main results can be downloaded from the CDS.

The discovery of an excess of light in the far-ultraviolet (UV) spectrum in elliptical galaxies was a major surprise in 1969. While it is now clear that this UV excess is caused by an old population of hot helium-burning stars without large hydrogen-rich envelopes rather than young stars, their origin has remained a mystery. Here we show that these stars most likely lost their envelopes because of binary interactions, similar to the hot subdwarf population in our own Galaxy. This has major implications for understanding the evolution of the UV excess and of elliptical galaxies in general. In particular, it implies that the UV excess is not a sign of age, as had been postulated previously, and predicts that it should not be strongly dependent on the metallicity of the population.

The discovery of a flux excess in the far-ultraviolet (UV) spectrum of elliptical galaxies was a major surprise in 1969. While it is now clear that this UV excess is caused by an old population of hot helium-burning stars without large hydrogen-rich envelopes, rather than young stars, their origin has remained a mystery. Here we show that these stars most likely lost their envelopes because of binary interactions, similar to the hot subdwarf population in our own Galaxy. We have developed an evolutionary population synthesis model for the far-UV excess of elliptical galaxies based on the binary model developed by Han et al (2002, 2003) for the formation of hot subdwarfs in our Galaxy. Despite its simplicity, it successfully reproduces most of the properties of elliptical galaxies with a UV excess: the range of observed UV excesses, both in $(1550-V)$ and $(2000-V)$, and their evolution with redshift. We also present colour-colour diagrams for use as diagnostic tools in the study of elliptical galaxies. The model has major implications for understanding the evolution of the UV excess and of elliptical galaxies in general. In particular, it implies that the UV excess is not a sign of age, as had been postulated previously, and predicts that it should not be strongly dependent on the metallicity of the population, but exists universally from dwarf ellipticals to giant ellipticals.

Using a simple analytic method, we calculate the impact effect between the ejecta of a SN Ia and its companion to survey the influence of initial parameters of the progenitor's system, which is useful for searching the companion in a explosion remnant. The companion models are obtained from Eggleton's evolution code.The results are divided into two groups based on mass transfer stage. For a given condition, more hydrogen-rich material is stripped from the envelope of a Hertzsprung-gap companion than that of a main-sequence companion, while a larger kick velocity and a larger luminosity are gained for a main-sequence companion. The kick velocity is too low to significantly affect the final spatial velocity of the companion, which is mainly affected by the initial parameters of the progenitor systems. The spatial velocity of the stripped material has an upper limit within the range of 8000 - 9500 km/s, which only depends on the total kinetic energy of the explosion.The stripped mass, the ratio of the stripped mass to the companion mass and the kick velocity of the companion all significantly depend on the initial companion mass and orbital period. Our model may naturally explain the spatial velocity of the star G in the remnant of Tycho's supernova, while an energy-loss mechanism is needed to interpret its luminosity.

The discovery of an excess of light in the far-ultraviolet (UV) in 1969 in elliptical galaxies was a major surprise. While it is now clear that this UV excess (UV-upturn) is probably caused by an old population of helium-burning stars. Han et al (2002, 2003) proposed a binary model for the formation of hot subdwarfs (helium burning stars) and the model can reproduce the observations in our Galaxy. By applying the binary model to the study of evolutionary population synthesis, we have obtained an {\it a priori} model for the UV-upturn of elliptical galaxies. The model shows that the UV-upturn is most likely resulted from binary interactions and it is universal (not very much metallicity-dependant) in ellipticals. This has major implications for understanding the evolution of the UV-upturn and elliptical galaxies in general; contrary to previous postulates, it implies that the UV-upturn is not a sign of age, but could be a potentially powerful indicator for a recent minor burst of star-forming activity.

Colors are usually not used for constraining stellar populations because they are thought to have the well-known agemetallicity degeneracy, but some recent works show that colors can also be used. A simple stellar population synthesis model is widely used, but there is no analysis for its colors. We try to find colors that can potentially be used to determine the age and metallicity of stellar populations by the standard model. Principal component analysis and relative sensitive parameter techniques are used in this work. As a result, U-K, U-H, U-J, B-K, B-H, U-I, B-J, and V-K are found to be more important for studying populations than others. Pairs of colors such as B-K and B-V are found to be able to disentangle the stellar age-metallicity degeneracy via the high-resolution model, while pairs such as U-K and R-I may be used instead when the low-resolution model is used. Furthermore, the u-g and r-i colors of the low-resolution model seem to have the same potential, but there are no such colors for the high-resolution one. In conclusion, some colors have been shown to have the potential to determine the age and metallicity of stellar populations, but relative metallicity and age sensitivities of colors in different stellar population synthesis models are usually different. In addition, minor star formations will make star systems look younger and more metal rich than their dominating populations.

Blue straggler stars present as secure members in the Galactic open clusters form a major challenge to the conventional picture of evolutionary population synthesis based on the stellar evolution theory of single stars, as illustrated in our previous work. Expansion of our sample in the current work to include younger age clusters provides a larger data base to expose the question raised for the simple stellar population model. The working sample now includes 97 Galactic open clusters of ages ranging from 0.1 to several Gyrs. The contributions of blue straggler stars to the integrated light of the host clusters are calculated on an individual cluster base. A data base of observational constrained simple stellar population model is made which has a larger age coverage than our previous work. It is shown in this work that the general existence of blue stragglers in star clusters of our sample dramatically altered the predictions of convectional stellar population model in terms of spectral energy distribution. The integrated spectral energy distributions of the synthetic spectra of the clusters are enhanced towards shorter wavelengths, therefore the results of the present work will cast new lights in understanding the properties of stellar populations.

We have carried out a detailed study of symbiotic stars with white dwarf accretors by means of a population synthesis code. We estimate the total number of symbiotic stars with white dwarf accretors in the Galaxy as 1,200 - 15,000. This range is compatible with observational estimates. Two crucial physical parameters that define the birthrate and number of symbiotic stars are the efficiency of accretion by white dwarfs (which greatly depends on the separation of components after common envelope stage and stellar wind velocity) and the mass of the hydrogen layer which the white dwarf can accumulate prior to the hydrogen ignition. The theoretical estimate of the Galactic occurrence rate of symbiotic novae ranges from about 1.3 to about 13.5 per year out of which weak symbiotic novae comprise about 0.5 to 6.0 per year, depending on the model assumptions. We simulate the distributions of symbiotic stars over orbital periods, masses of components, mass-loss rates of cool components, mass-accretion rates of hot components and luminosity of components. Agreement with observations is reasonable.

We determine relative stellar ages and metallicities mainly for about 80 elliptical galaxies in low and high density environments via the latest binary stellar population (BSP) synthesis model and test a latest hierarchical formation model of elliptical galaxies which adopted the new $\Lambda$CDM cosmology for the first time. The stellar ages and metallicities of galaxies are estimated from two high-quality published spectra line indices, i.e. H$\beta$ and [MgFe]. The results show that elliptical galaxies have stellar populations older than 3.9 Gyr and more metal rich than 0.02. Most of our results are in agreement with predictions of the model: First, elliptical galaxies in denser environment are redder and have older populations than field galaxies. Second, elliptical galaxies with more massive stellar components are redder while have older and more metal rich populations than less massive ones. Third, the most massive galaxies are shown to have the oldest and most metal rich stars. However, some of our results are found to be different with predictions of the galaxy formation model, i.e. the metallicity distributions of low- and high-density elliptical galaxies and the relations relating to cluster-centric distance.

In this paper, we present the simulations for the primordial blue stragglers in the old open cluster M67 based on detailed modelling of the evolutionary processes. The principal aim is to discuss the contribution of mass transfer between the components of close binaries to the blue straggler population in M67. First, we followed the evolution of a binary of 1.4M$_\odot$+0.9M$_\odot$. The synthetic evolutionary track of the binary system revealed that a primordial blue straggler had a long lifetime in the observed blue straggler region of color-magnitude diagram. Second, a grid of models for close binary systems experiencing mass exchange were computed from 1Gyr to 6Gyr in order to account for primordial blue-straggler formation in a time sequence. Based on such a grid, Monte-Carlo simulations were applied for the old open cluster M67. Adopting appropriate orbital parameters, 4 primordial blue stragglers were predicted by our simulations. This was consistent with the observational fact that only a few blue stragglers in M67 were binaries with short orbital periods. An upper boundary of the primordial blue stragglers in the color-magnitude diagram (CMD) was defined and could be used to distinguish blue stragglers that were not formed via mass exchange. Using the grid of binary models, the orbital periods of the primordial BSs could be predicted. Compared with the observations, it is clear that the mechanism discussed in this work alone cannot fully predict the blue straggler population in M67. There must be several other processes also involved in the formation of the observed blue stragglers in M67.

We re-assess the question of a systematic time delay between the formation of the progenitor and its explosion in a type Ia supernova (SN Ia) using the Hubble Higher-z Supernova Search sample (Strolger et al. 2004). While the previous analysis indicated a significant time delay, with a most likely value of 3.4 Gyr, effectively ruling out all previously proposed progenitor models, our analysis shows that the time-delay estimate is dominated by systematic errors, in particular due to uncertainties in the star-formation history. We find that none of the popular progenitor models under consideration can be ruled out with any significant degree of confidence. The inferred time delay is mainly determined by the peak in the assumed star-formation history. We show that, even with a much larger Supernova sample, the time delay distribution cannot be reliably reconstructed without better constraints on the star-formation history.

Supernova 2002ic was an atypical Type Ia supernova (SN Ia) with evidence for substantial amounts of hydrogen associated with the system. Contrary to previous claims, we show that its unusual properties can be understood within the framework of one of the most favoured progenitor models, the so-called supersoft channel. This requires that the donor star was initially relatively massive ($\sim 3\Ms$) and that the system experienced a delayed dynamical instability, leading to a large amount of mass loss from the system in the last few $10^4 $yr before the explosion. This can produce the inferred hydrogen-rich circumstellar environment, most likely with a disc-like geometry. However, to apply these models requires a larger accretion efficiency onto the white dwarf than is assumed in present parameterizations. If this is confirmed, it would most likely increase estimates for the frequency of the single-degenerate channel. Based on population synthesis simulations we estimate that not more than 1 in 100 SNe Ia should belong to this subgroup of SNe Ia.

Thanks to a stellar evolution code able to compute through the C-flash we link the binary population synthesis of single degenerate progenitors of type Ia supernovae (SNe Ia) to their physical condition at the time of ignition. We show that there is a large range of possible ignition densities and we detail how their probability distribution depends on the accretion properties. The low density peak of this distribution qualitatively reminds of the clustering of the luminosities of Branch-normal SNe Ia. We tighten the possible range of initial physical conditions for explosion models: they form a one-parameter family, independent of the metallicity. We discuss how these results may be modified if we were to relax our hypothesis of a permanent Hachisu wind or if we were to include electron captures.

Using EPS we present high resolution (0.3 \AA, HRes) ISEDs from 3000 to 7000 \AA and Lick/IDS absorption-line indices, for an extensive set of instantaneous burst binary stellar populations (BSPs) with binary interactions. The ages of the BSPs are in the range 1 - 15 Gyr and the metallicities are in the range 0.004 - 0.03. This HRes synthesis results can satisfy the needs of modern spectroscopic galaxy surveys, and are available on request. By comparing the synthetic continuum of BSPs at high and low resolution (LRes) we show that there is a good agreement for Z=0.02 and a tolerable disagreement for non-solar metallicity. The strength of the Balmer lines at HRes is greater than that at LRes for all metallicities. The comparison of Lick/IDS absorption-line indices at LRes and HRes, both of which are obtained by the fitting functions (FFs), shows that the discrepancies in all indices except for TiO_1 and TiO_2 are insignificant for BSPs with Z=0.004 and Z=0.02. The HRes Ca4227, Fe5015 and Mg_b indices are redder than the corresponding LRes one for BSPs with Z=0.01 and Z=0.03, this effect lowers the derived age and metallicity of the BSP. The high resolution Mg_1, Fe5709 and Fe5782 indices are bluer than those at LRes, it raises the age and metallicity. The discrepancy in these six indices is greater for BSPs with Z=0.03 in comparison to Z=0.01. At HRes we compare the Lick/IDS spectral absorption indices obtained by using the FFs with those measured directly from the synthetic spectra, and see that Ca4455, Fe4668, Mg_b and Na D indices obtained by the use of the FFs are redder for all metallicities, Fe5709 is redder at Z=0.03 and becomes to be bluer at Z=0.01 and 0.004, and other indices are bluer for all metallicities than the corresponding values measured directly from the synthetic spectra.

In a preceding paper, using Eggleton's evolution code we have discussed the structure and evolution of low-mass W UMa type contact binaries without angular momentum loss (AML). The models exhibit cyclic behavior about a state of marginal contact on a thermal time-scale. Part of the time of each cycle is spent in contact and part in a semi-detached state. According to observations, W UMa systems suffer AML. We present the models of low-mass contact binaries with AML due to gravitational wave radiation (GR) or magnetic stellar wind (MSW) are presented. We find that gravitational radiation cannot prevent the cyclic evolution of W UMa systems, and the effect of gravitational radiation on the cyclic behavior of contact binary evolution is almost negligible. We also find that the most likely AML mechanism for W UMa systems is magnetic braking, and that magnetic braking effects can increase the period of the cyclic evolution, and shorten the fraction of the time spent in the poor thermal contact state. If W UMa stars do not undergo cyclic evolution, and their angular momentum loss is caused simultaneously by MSW of both components, we find that the value of the parameter, $\lambda$, should be taken a larger value in comparison with those derived from observations of single stars. This indicates that the AML efficiency in W UMa systems may be lowered in comparison with non-contact stars because of the less mass contained in the convective envelopes of the components in W UMa systems. If W UMa systems lose their angular momentum at a constant rate. An angular momentum rate of $\frac{{\rm dln}J}{{\rm d}t}\approx 1.6\times 10^{-9} {\rm yr^{-1}}$ can prevent the cyclic behaviour of the model, and the model can keep in good contact with an essentially constant depth of contact.

Using evolutionary population synthesis we present integrated colours, integrated spectral energy distributions and absorption-line indices defined by the Lick Observatory image dissector scanner (referred to as Lick/IDS) system, for an extensive set of instantaneous burst binary stellar populations with and without binary interactions. The ages of the populations are in the range $1-15 $Gyr and the metallicities are in the range $0.0001-0.03$. By comparing the results for populations with and without binary interactions we show that the inclusion of binary interactions makes the integrated $\rm U-B$, $\rm B-V$, $\rm V-R$ and $\rm R-I$ colours and all Lick/IDS spectral absorption indices (except for $\rm H_\beta$) substantially smaller. In other words binary evolution makes a population appear bluer. This effect raises the derived age and metallicity of the population. We calculate several sets of additional solar-metallicity binary stellar populations to explore the influence of input parameters to the binary evolution algorithm (the common-envelope ejection efficiency and the stellar wind mass-loss rate) on the resulting integrated colours. We also look at the dependence on the choice of distribution functions used to generate the initial binary population. The results show that variations in the choice of input model parameters and distributions can significantly affect the results. However, comparing the discrepancies that exist between the colours of various models, we find that the differences are less than those produced between the models with and those without binary interactions. Therefore it is very necessary to consider binary interactions in order to draw accurate conclusions from evolutionary population synthesis work.

We present the analysis of a high-resolution, high-quality sample of optical spectra for 76 subdwarf B (sdB) stars from the ESO Supernova Ia Progenitor Survey (SPY). Effective temperature, surface gravity, and photospheric helium abundance are determined simultaneously by fitting the profiles of hydrogen and helium lines using synthetic spectra calculated from LTE and NLTE model atmospheres. We perform a detailed comparison of our measurements with theoretical calculations, both for single star evolution and for binary population synthesis models of close binary evolution. The luminosity evolution given by the standard EHB evolutionary tracks from Dorman et al. (1993) shows an overall agreement in shape with our observations, although a constant offset in luminosity exists. The various simulation sets for binary formation channels of sdB stars calculated by Han et al. (2003) are compared individually to our data for testing our current understanding of sdB formation processes and the physical effects involved. The best-matching sets manage to reproduce the observed sdB distribution in the temperature-gravity-plane well. However, they do not match the observed cumulative luminosity function, indicating that theoretical improvement is necessary. We also investigate composite-spectrum objects showing clear signatures of a cool companion with optical and infrared photometry. These stars have cool main sequence companions of spectral types F to K. Typical helium abundances of composite and non-composite sdB stars do not differ.

We studied the effects of surface composition and thermohaline mixing caused by secular instability on the accreting components for low--mass binaries and applied the results on a short-orbital-period BS F190 in the old cluster M67. The results indicate no distinction in surface composition between the models with and without thermohaline mixing during Roche lobe overflow, but we still see the divergences of evolutionary tracks on HRD and CMD. The change of surface composition makes the gainer bluer and smaller than the ones with original surface composition while thermohaline mixing lessens the effect slightly. If thermohaline mixing were to act instantaneously, the effect would be lessened more. Our calculation shows that case A and case B mass transfer may produce BSs in short- or relatively short-orbital-period binaries (including Algol systems), and that CNO abundance abnormalities could be observed in these products. Our simulation of F190 shows that the primary's mass $M_{\rm 1i}$ of the appropriate models is located in the range of 1.40 to 1.45$M_{\odot}$ with initial mass ratio $q_{\rm i}=1.5$ and initial orbital period $P_{\rm i}=0.8$ days, indicating that case A is a more likely evolutionary channel than case B to form this object. The simulation also shows that it is very likely that F190 is still in a slow stage of mass transfer. As a consequence, obvious CNO abundance abnormalities should be observed for the object.

The structure and evolution of low mass W UMa type contact binaries are discussed by employing Eggleton's stellar evolution code. Assuming that these systems completely satisfy Roche geometry, we calculate the relative radii of both components of contact binaries in different contact depth between inner and outer Roche lobes. We obtain a radius grid of contact binaries, and can ensure the surfaces of two components lying on an equipotential surface by interpolation using this radius grid when we follow the evolution of the contact binaries. We assume that the energy transfer takes place in the different regions of the common envelope to investigate the effects of the region of energy transfer on the structure and evolution of contact binaries. We find that the region of energy transfer has significant influence on the structure and evolution of contact binaries, and conclude that the energy transfer may occur in the outermost layers of the common convective envelope for W-type systems, and this transfer takes place in the deeper layers of the common envelope for A-type systems. Meanwhile, if we assume that the energy transfer takes place in the outermost layers for our model with low total mass, and find that our model steadily evolves towards a system with a smaller mass ratio and a deeper envelope, suggesting that some A-type W UMa systems with low total mass could be considered as the later evolutionary stages of W-subtype systems, and that the surface temperature of the secondary excesses that of the primary during the time when the primary expands rapidly, or the secondary contracts rapidly, suggesting that W-subtype systems may be caused by expansion of the primary, or by the contraction of the secondary.

Using evolutionary population synthesis, we present integrated spectral energy distributions and absorption-line indices defined by the Lick Observatory image dissector scanner (referred to as Lick/IDS) system, for an extensive set of instantaneous burst single stellar populations (SSPs). The ages of the SSPs are in the range 1-19 Gyr and the metallicities [Fe/H] are in the range -2.3 - 0.2. Our models use the rapid single stellar evolution algorithm of Hurley, Pols and Tout for the stellar evolutionary tracks, the empirical and semi-empirical calibrated BaSeL-2.0 model of Lejeune, Cuisinier and Buser for the library of stellar spectra and the empirical fitting functions of Worthey, Faber, Gonzalez and Burstein for the Lick/IDS spectral absorption-feature indices. Applying our synthetic Lick/IDS absorption-line indices to the merit function, we obtain the age and the metallicity of the central region of M32, which can be well explained by an instantaneous SSP with an age of 6.5 Gyr and a metallicity similar to solar. Applying the derived age and the metallicity from the merit function to a number of index-index diagrams, we find that the plots of Hbeta-Fe5015 and Hbeta-Fe5782 are the best index-index diagrams from which we can directly obtain reasonable age and metallicity.

We present integrated colours, integrated spectral energy distributions, and absorption-line indices, for instantaneous burst solar-metallicity binary stellar populations with ages in the range 1 - 15Gyr. By comparing the results for populations with and without binary interactions we show that the inclusion of binary interactions makes the appearance of the population substantially bluer -- this is the case for each of the quantities we have considered. This effect raises the derived age and metallicity of the population. Therefore it is necessary to consider binary interactions in order to draw accurate conclusions from evolutionary population synthesis work.

We have carried out a detailed study of one of the most popular evolutionary channels for the production of Type Ia supernova (SN Ia) progenitors, the semi-degenerate channel (CO+MS), where a carbon/oxygen (CO) white dwarf (WD) accretes matter from an unevolved or slightly evolved non-degenerate star until it reaches the Chandrasekhar mass limit. Employing Eggleton's stellar evolution code and adopting the prescription of Hachisu et al. (1999) for the accretion efficiency, we have carried out full binary evolution calculations for about 2300 close WD binary systems and mapped out the initial parameters in the orbital period -- secondary mass ($P$--$M_2$) plane (for a range of WD masses) which lead to a successful Type Ia supernova. We obtained accurate, analytical fitting formulae to describe this parameter range which can be used for binary population synthesis (BPS) studies. The contours in the $P$--$M_2$ plane differ from those obtained by Hachisu et al. (1999) for low-mass CO WDs, which are more common than massive CO WDs. We show that white dwarfs with a mass as low as $0.67M_\odot$ can accrete efficiently and reach the Chandrasekhar limit. We have implemented these results in a BPS study to obtain the birthrates for SNe Ia and the evolution of birthrates with time of SNe Ia for both a constant star formation rate and a single star burst. The birthrates are somewhat lower than (but comparable to) those inferred observationally.

The present work is designed to explore the evolutionary and pulsational properties of low-mass white dwarfs with carbon/oxygen cores. In particular, we follow the evolution of a 0.33 Msun white dwarf remnant in a self-consistent way with the predictions of nuclear burning, element diffusion and the history of the white dwarf progenitor. Attention is focused on the occurrence of hydrogen shell flashes induced by diffusion processes during cooling phases. The evolutionary stages prior to the white dwarf formation are also fully accounted for by computing the conservative binary evolution of an initially 2.5-Msun Pop. I star with a 1.25 Msun companion, and period P_i= 3 days. Evolution is followed down to the domain of the ZZ Ceti stars on the white dwarf cooling branch. We find that chemical diffusion induces the occurrence of an additional hydrogen thermonuclear flash which leads to stellar models with thin hydrogen envelopes. As a result, a fast cooling is encountered at advanced stages of evolution. In addition, we explore the adiabatic pulsational properties of the resulting white dwarf models. As compared with their helium-core counterparts, low-mass oxygen-core white dwarfs are characterized by a pulsational spectrum much more featured, an aspect which could eventually be used for distinguishing both types of stars if low-mass white dwarfs were in fact found to pulsate as ZZ Ceti-type variables. Finally, we perform a non-adiabatic pulsational analysis on the resulting carbon/oxygen low-mass white dwarf models.

We present the results of a systematic study of the formation and evolution of binaries containing black holes and normal-star companions with a wide range of masses. We first reexamine the standard formation scenario for close black-hole binaries, where the spiral-in of the companion in the envelope of a massive star causes the ejection of the envelope. We estimate the formation rates for different companion masses and different assumptions about the common-envelope structure and other model parameters. We find that black-hole binaries with intermediate- and high-mass secondaries can form for a wide range of assumptions, while black-hole binaries with low-mass secondaries can only form with apparently unrealistic assumptions (in agreement with previous studies). We then present detailed binary evolution sequences for black-hole binaries with secondaries of 2 to 17 Msun and demonstrate that in these systems the black hole can accrete appreciably even if accretion is Eddington limited (up to 7 Msun for an initial black-hole mass of 10 Msun) and that the black holes can be spun up significantly in the process. We discuss the implications of these calculations for well-studied black-hole binaries (in particular GRS 1915+105), ultra-luminous X-ray sources and Cygnus X-1. Finally, we discuss how some of the assumptions in the standard model could be relaxed to allow the formation of low-mass, short-period black-hole binaries which appear to be very abundant in Nature. (Abstract abridged)

We have carried out a detailed binary populations synthesis (BPS) study of the formation of subdwarf B (sdB) stars and related objects (sdO, sdOB stars) using the latest version of the BPS code developed by Han et al.(1994, 1995a, 1995b, 1998, 2001). We systematically investigate the importance of the five main evolutionary channels in which the sdB stars form after one or two common-envelope (CE) phases, one or two phases of stable Roche-lobe overflow (RLOF) or as the result of the merger of two helium white dwarfs (WD) (see Han et al. 2002, Paper I). Our best BPS model can satisfactorily explain the main observational characteristics of sdB stars, in particular their distributions in the orbital period - minimum companion mass diagram and in the effective temperature - surface gravity diagram, their distributions of orbital period, log (g theta^4), and mass function, their binary fraction and the fraction of sdB binaries with WD companions, their birthrates and their space density. We obtain a Galactic formation rate, a total number in the Galaxy, the intrinsic binary fraction for sdB stars. We also predict a distribution of masses for sdB stars that is wider than is commonly assumed and that some sdB stars have companions of spectral type as early as B. The percentage of A type stars with sdB companions can in principle be used to constrain some of the important parameters in the binary evolution model. We conclude that (a) the first RLOF phase needs to be more stable than is commonly assumed; (b) mass transfer in the first stable RLOF phase is non-conservative, and the mass lost from the system takes away a specific angular momentum similar to that of the system; (c) common-envelope ejection is very efficient.

Subdwarf B (sdB) stars (and related sdO/sdOB stars) are believed to be helium core-burning objects with very thin hydrogen-rich envelopes. In recent years it has become increasingly clear from observational surveys that a large fraction of these objects are members of binary systems. To better understand their formation, we here present the results of a detailed investigation of the three main binary evolution channels that can lead to the formation of sdB stars: the common envelope (CE) ejection channel, the stable Roche lobe overflow (RLOF) channel and the double helium white dwarfs (WDs) merger channel. We obtained the conditions for the formation of sdB stars from each of these channels using detailed stellar and binary evolution calculations where we modelled the detailed evolution of sdB stars and carried out simplified binary population synthesis simulations. The observed period distribution of sdB stars in compact binaries strongly constrains the CE ejection parameters. We also present the distribution of sdB stars in the $T_{\rm eff}$ - $\log g$ diagram, the Hertzsprung-Russell diagram and the distribution of mass functions.

We present the results of a systematic study of cataclysmic variables (CVs) and related systems, combining detailed binary-population synthesis (BPS) models with a grid of 120 binary evolution sequences calculated with a Henyey-type stellar evolution code. In these sequences, we used 3 masses for the white dwarf (0.6, 0.8, 1.0 Msun) and seven masses for the donor star in the range of 0.6-1.4 Msun. The shortest orbital periods were chosen to have initially unevolved secondaries, and the longest orbital period for each secondary mass was taken to be just longer than the bifurcation period (16 - 22 hr), beyond which systems evolve towards long orbital periods. These calculations show that systems which start with evolved secondaries near the end or just after their main-sequence phase become ultra-compact systems with periods as short as 7 min. These systems are excellent candidates for AM CVn stars. Using a standard BPS code, we show how the properties of CVs at the beginning of mass transfer depend on the efficiency for common-envelope (CE) ejection and the efficiency of magnetic braking. In our standard model, where CE ejection is efficient, some 10 per cent of all CVs have initially evolved secondaries (with a central hydrogen abundance X_c < 0.4) and ultimately become ultra-compact systems (implying a Galactic birthrate for AM CVn-like stars of 10^{-3} yr^{-1}). Almost all CVs with orbital periods longer than 5 hr are found to have initially evolved or relatively massive secondaries. We show that their distribution of effective temperatures is in good agreement with the distribution of spectral types obtained by Beuermann et al. (1998).

Using Eggleton's stellar evolution code, we carry out 150 runs of Pop I binary evolution calculations, with the initial primary mass between 1 and 8 solar masses the initial mass ratio between 1.1 and 4, and the onset of Roche lobe overflow (RLOF) at an early, middle, or late Hertzsprung-gap stage. We assume that RLOF is conservative in the calculations, and find that the remnant mass of the primary may change by more than 40 per cent over the range of initial mass ratio or orbital period, for a given primary mass. This is contrary to the often-held belief that the remnant mass depends only on the progenitor mass if mass transfer begins in the Hertzsprung gap. We fit a formula, with an error less than 3.6 per cent, for the remnant (white dwarf) mass as a function of the initial mass of the primary, the initial mass ratio, and the radius of the primary at the onset of RLOF. We also find that a carbon-oxygen white dwarf with mass as low as 0.33 solar masses may be formed if the primary's initial mass is around 2.5 solar masses.

WD1704+481 is a visual binary in which both components are white dwarfs. We present spectra of the H-alpha line of both stars which show that one component (WD1704+481.2 = Sanduleak B = GR 577) is a close binary with two white dwarf components. Thus, WD1704+481 is the first known triple degenerate star. From radial velocity measurements of the close binary we find an orbital period of 0.1448d, a mass ratio, q=Mbright/Mfaint of q=0.70+-0.03 and a difference in the gravitational redshifts of 11.5+-2.3km/s. The masses of the close pair of white dwarfs predicted by the mass ratio and gravitational redshift difference combined with theoretical cooling curves are 0.39+-0.05 solar mass and 0.56+-0.07 solar masses. WD1704+481 is therefore also likely to be the first example of a double degenerate in which the less massive white dwarf is composed of helium and the other white dwarf is composed of carbon and oxygen.

We present a simple and efficient, yet reasonably accurate, equation of state, which at the moderately low temperatures and high densities found in the interiors of stars less massive than the Sun is substantially more accurate than its predecessor by Eggleton, Faulkner & Flannery. Along with the most recently available values in tabular form of opacities, neutrino loss rates, and nuclear reaction rates for a selection of the most important reactions, this provides a convenient package of input physics for stellar modelling. We briefly discuss a few results obtained with the updated stellar evolution code.