Spin Evolution of Neutron Stars in OB/X-ray Binaries

We have investigated the relation between the orbital period Porb and the spin period Ps of neutron stars in OB/X-ray binaries. By simulating the time-development of the mass loss rate and radius expansion of a 20M☉ donor star, we have calculated the detailed spin evolution of the neutron star before steady wind accretion occurs (that is, when the break spin period is reached), or when the OB star begins evolving off the main sequence or has filled its Roche lobe. Our results are compatible with the observations of OB/X-ray binaries. We find that in relatively narrow systems with orbital periods less than tens of days, neutron stars with initial magnetic field B0 stronger than about 3×1012 G can reach the break spin period to allow steady wind accretion in the main sequence time, whereas neutron stars with B0 < 3×1012 G and/or in wide systems would still be in one of the pulsar, rapid rotator or propeller phases when the companion evolves off the main sequence or fills its Roche lobe. Our results may help understand the various characteristics of the observed OB/neutron star binaries along with their distributions in the Ps-Porb diagram.

Light curve solutions of six eclipsing binaries at the lower limit of periods for W UMa stars

Photometric observations are presented in V and I bands of six eclipsing binaries at the lower limit of the orbital periods for W UMa stars. Three of them are newly discovered eclipsing systems. The light curve solutions reveal that all shortperiod targets are contact or overcontact binaries and six new binaries are added to the family of short-period systems with estimated parameters. Four binaries have com- ponents that are equal in size and a mass ratio near 1. The phase variability shown by the V-I colors of all targets may be explained by lower temperatures on their back surfaces than those on their side surfaces. Five systems exhibit the O＇Connell effect that can be modeled by cool spots on the side surfaces of their primary components. The light curves of V1067 Her in 2011 and 2012 are fitted by diametrically opposite spots. Applying the criteria for subdivision of W UMa stars to our targets leads to ambiguous results.

Periodic Variation Studies of the Two Short Period W UMa-type Eclipsing Binaries: LX Lyn and V0853 Aur

In this paper,new light curves(LCs) of contact eclipsing binary(CEB) systems LX Lyn and V0853 Aur are presented and analyzed by using the 2015 version of the Wilson-Devinney(W-D) code.In order to explain their asymmetric LCs,cool starspots on the components were employed.It is suggested that their fill-out degrees are f=12.0%(LX Lyn) and f=26.3%(V0853 Aur).At the same time,we found that LX Lyn is a W-type eclipsing binary(EB) with an orbital inclination of i=84°.88 and a mass ratio of q=2.31.V0853 Aur is also a W-type CEB with a mass ratio of q=2.77 and an orbital inclination of i= 79°.26.Based on all available times of light minimum,their orbital period changes are studied by using the O-C method.The O-C diagram of LX Lyn reveals a cyclic oscillation with a period of about 14.84 yr and an amplitude of 0.0019 days,which can be explained by the light-travel time effect(LTTE) due to the presence of a third body with a minimum mass of0.06M_⊙.For V0853 Aur,it is discovered that the O-C diagram of the system also shows a cyclic oscillation with a period of 9.64 yr and an amplitude of 0.03365 days.The cyclic oscillation of V0853 Aur can be attributed to the LTTE by means of a third body with a mass no less than 3.77M_⊙.The third body may play an important role in the formation and evolution of these systems.

Variable Stars in the 50BiN Open Cluster Survey.Ⅲ.NGC 884

Open clusters are the basic building blocks that serve as a laboratory for the study of young stellar populations in the Milky Way.Variable stars in open clusters provide a unique way to accurately probe the internal structure,temporal and dynamical evolutionary stages of individual stars and the host cluster.The most powerful tool for such studies is time-domain photometric observations.This paper follows the route of our previous work,concentrating on a photometric search for variable stars in NGC 884.The target cluster is the companion of NGC869,forming the well-known double cluster system that is gravitationally bound.From the observation run in 2016 November,a total of 9247 B-band CCD images and 8218Ⅴ-band CCD images were obtained.We detected a total of 15 stars with variability in visual brightness,including five Be stars,three eclipsing binaries,and seven of unknown types.Two new variable stars were discovered in this work.We also compared the variable star content of NGC 884 with its companion NGC 869.

Searching forγ-ray counterparts to very faint X-ray transient neutron star binaries

Very faint X-ray transients （VFXTs） are a group of X-ray binaries with low luminosities, displaying peak X-ray luminosities during their outbursts of only 1034-1036 erg s^-1. Using γ-ray data obtained with the Large Area Telescope （LAT） onboard the Fermi Gamma-Ray Space Telescope （Fermi）, we investigate their possible nature of containing rotation-powered pulsars, or more specifically being transitional millisecond pulsars （MSPs）. Among more than 40 known VFXTs, we select 12 neutron star systems. We analyze the LAT data for the fields of 12 VFXTs in the energy range 0.2-300GeV, but do not find any counterparts likely detected by Fermi. We obtain luminosity upper limits for the 12 sources. While the distances to the sources are largely uncertain, the upper limits are comparable to the luminosities of two transitional systems, PSR J1023-0038 and XSS J12270-4859. From our study, we conclude that no evidence is found at γ-rays for the suggestion that some VFXTs could contain rotation-powered MSPs （or be transitional MSP systems）.

The Binarity of Early-type Stars from LAMOST medium-resolution Spectroscopic Survey

Massive binaries play significant roles in many fields.Identifying massive stars,particularly massive binaries,is of great importance.In this paper,by adopting the technique of measuring the equivalent widths of several spectral lines,we identified 9382 early-type stars from the LAMOST medium-resolution survey and divided the sample into four groups,T1(~O-B4),T2(~B5),T3(~B7)and T4(~B8-A).The relative radial velocities RVrelwere calculated using Maximum Likelihood Estimation.The stars with significant changes of RVreland at least larger than 15.57 km s-1 were identified as spectroscopic binaries.We found that the observed spectroscopic binary fractions for the four groups are 24.6%±0.5%,20.8%±0.6%,13.7%±0.3%and 7.4%±0.3%,respectively.Assuming that orbital period(P)and mass ratio(q)have intrinsic distributions as f(P)∝Pπ(1<P<1000 days)and f(q)∝qκ(0.1<q<1),respectively,we conducted a series of Monte-Carlo simulations to correct observational biases for estimating the intrinsic multiplicity properties.The results show that the intrinsic binary fractions for the four groups are 68%±8%,52%±3%,44%±6%and 44%±6%,respectively.The best estimated values forπare-1±0.1,-1.1±0.05,-1.1±0.1 and-0.6±0.05,respectively.Theκcannot be constrained for groups T1 and T2 and is-2.4±0.3 for group T3 and-1.6±0.3 for group T4.We confirmed the relationship of a decreasing trend in binary fractions toward late-type stars.No correlation between the spectral type and orbital period distribution has been found yet,possibly due to the limitation of observational cadence.

The formation of neutron star systems through accretion-induced collapse in white-dwarf binaries

The accretion-induced collapse(AIC)scenario was proposed 40 years ago as an evolutionary end state of oxygen-neon white dwarfs(ONe WDs),linking them to the formation of neutron star(NS)systems.However,there has been no direct detection of any AIC event so far,even though there exists a lot of indirect observational evidence.Meanwhile,the evolutionary pathways resulting in NS formation through AIC are still not thoroughly investigated.In this article,we review recent studies on the two classic progenitor models of AIC events,i.e.,the single-degenerate model(including the ONe WD+MS/RG/He star channels and the CO WD+He star channel)and the double-degenerate model(including the double CO WD channel,the double ONe WD channel and the ONe WD+CO WD channel).Recent progress on these progenitor models is reviewed,including the evolutionary scenarios leading to AIC events,the initial parameter space for producing AIC events and the related objects(e.g.,the pre-AIC systems and the post-AIC systems).For the single-degenerate model,the pre-AIC systems(i.e.,the progenitor systems of AIC events)could potentially be identified as supersoft X-ray sources,symbiotics and cataclysmic variables(such as classical novae,recurrent novae,Ne novae and He novae)in the observations,whereas the post-AIC systems(i.e.,NS systems)could potentially be identified as low-/intermediate-mass X-ray binaries,and the resulting low-/intermediate-mass binary pulsars,most notably millisecond pulsars.For the double-degenerate model,the pre-AIC systems are close double WDs with short orbital periods,whereas the post-AIC systems are single isolated NSs that may correspond to a specific kind of NS with peculiar properties.We also review the predicted rates of AIC events,the mass distribution of NSs produced via AIC and the gravitational wave(GW)signals from double WDs that are potential GW sources in the Galaxy in the context of future spacebased GW detectors,such as LISA,TianQin,Taiji,etc.Recent theoretical and observational constraints on the detection of AIC events are summarized.In order to confirm the existence of the AIC process,and resolve this long-term issue presented by current stellar evolution theories,more numerical simulations and observational identifications are required.

Formation and Destiny of White Dwarf and Be Star Binaries

The binary systems consisting of a Be star and a white dwarf(Be WDs) are very interesting.They can originate from the binaries composed of a Be star and a subdwarf O or B star(Besd OBs),and they can merge into red giants via luminous red nova or can evolve into double WD potentially detected by the LISA mission.Using the method of population synthesis,we investigate the formation and the destiny of Be WDs,and discuss the effects of the metallicity(Z) and the common envelope evolution parameters.We find that Besd OBs are significant progenitors of Be WDs.About 30%(Z = 0.0001)-50%(Z = 0.02) of Be WDs come from Besd OBs.About 60%(Z = 0.0001)-70%(Z = 0.02) of Be WDs turn into red giants via a merger between a WD and a non-degenerated star.About 30%(Z = 0.0001)-40%(Z = 0.02) of Be WDs evolve into double WDs which are potential gravitational waves of the LISA mission at a frequency band between about 3 × 10^(-3)and 3 × 10^(-2)Hz.The common envelope evolution parameter introduces an uncertainty with a factor of about 1.3 on Be WD populations in our simulations.

Fractions of Compact Object Binaries in Star Clusters:Theoretical Predictions

The binary population in field stars and star clusters contributes to the formation of gravitational wave(GW)sources.However,the fraction of compact-object binaries(CBs),which is an important feature parameter of binary populations,is still difficult to measure and very uncertain.This paper predicts the fractions of important CBs and semi-compact object binaries(SCBs) making use of an advanced stellar population synthesis technique.A comparison with the result of N-body simulation is also presented.It is found that most CBs are formed within about 500 Myr after the starburst.The fractions of CBs and SCBs are demonstrated to correlate with stellar metallicity.The higher the metallicity becomes,the smaller the fraction of black hole binaries(BHBs),neutron star binaries(NSBs) and SCBs.This suggests that the GW sources of BHBs and NSBs are more likely to form in metal-poor environments.However,the fraction of black hole-neutron star binaries is shown to be larger for metalrich populations on average.

Quark-novae in neutron star - white dwarf binaries: a model for luminous (spin-down powered) sub-Chandrasekhar-mass Type Ia supernovae?

We show that, by appealing to a Quark-Nova （QN） in a tight binary system containing a massive neutron star and a CO white dwarf （WD）, a Type Ia explosion could occur. The QN ejecta collides with the WD, driving a shock that triggers carbon burning under degenerate conditions （the QN-Ia）. The conditions in the compressed low-mass WD （MwD 〈 0.9 M） in our model mimic those of a Chandrasekhar mass WD. The spin-down luminosity from the QN compact remnant （the quark star） pro- vides additional power that makes the QN-Ia light-curve brighter and broader than a standard SN-Ia with similar 56Ni yield. In QNe-Ia, photometry and spectroscopy are not necessarily linked since the kinetic energy of the ejecta has a contribution from spin-down power and nuclear decay. Although QNe-Ia may not obey the Phillips relationship, their brightness and their relatively ＂normal looking＂ light-curves mean they could be included in the cosmological sample. Light-curve fitters would be con- fused by the discrepancy between spectroscopy at peak and photometry and would correct for it by effectively brightening or dimming the QNe-Ia apparent magnitudes, thus over- or under-estimating the true magnitude of these spin-down powered SNe-Ia. Contamination of QNe-Ia in samples of SNe-Ia used for cosmological analyses could systematically bias measurements of cosmological parameters if QNe-Ia are numerous enough at high-redshift. The strong mixing induced by spin-down wind combined with the low 56Ni yields in QNe-Ia means that these would lack a secondary maximum in the/-band despite their luminous nature. We discuss possible QNe-Ia progenitors.

Variable binaries and variables in binaries in the Binary star DataBase

The Binary star DataBase(BDB, http://bdb.inasan.ru) combines data from catalogs of binary and multiple stars of all observational types. There is a number of ways for variable stars to form or to be a part of binary or multiple systems. We describe how such stars are represented in the database.

Flaring activity from quiescent states in neutron-star X-ray binaries

We examine systematically the observed X-ray luminosity jumps(or flares) from quiescent states in millisecond binary pulsars(MSBPs) and high-mass X-ray binary pulsars(HMXBPs). We rely on the published X-ray light curves of seven pulsars: four HMXBPs, two MSBPs and the ultraluminous X-ray pulsar M82 X-2. We discuss the physics of their flaring activities or lack thereof, paying special attention to their emission properties when they are found on the propeller line, inside the Corbet gap or near the light-cylinder barrier. We provide guiding principles for future interpretations of faint X-ray observations, as well as a method of constraining the propeller lines and the dipolar surface magnetic fields of pulsars using a variety of quiescent states. In the process, we clarify some disturbing inaccuracies that have made their way into the published literature.

A Brief Review for Quasi-Periodic Oscillations in Neutron-Star Low-Mass X-Ray Binaries

In accreting neutron star (NS) low-mass X-ray binaries (LMXBs), the turbulent flow in accretion disk may show magnetic structures. Its emission will vary in time due to inhomogeneous motions through and with the accretion flow. These emissions contribute to considerable X-ray variability on a wide range of timescales in all wavelengths, and down to milliseconds. In this article, we give a brief review for quasi-periodic oscillations (QPOs), one of a periodic X-ray variability, in NS/ LMXBs. Firstly, we give a brief introduction to NS/LMXBs and the fruitful QPO components. As an example, the energy dependence of normal branch oscillations in Scorpius X-1 is discussed. We mostly focus on the properties and mechanism of kilohertz QPOs—the fastest variability components that have the same order as the dynamical timescales of the innermost regions of accretion flow. Finally, we discuss the success and questions for theoretical interpretations and present the possible entry for investigation of nature of QPOs.

Physical properties and catalog of EW-type eclipsing binaries observed byLAMOST

EW-type eclipsing binaries(hereafter called EWs)are strong interacting systems in which both component stars usually fill their critical Roche lobes and share a common envelope.Numerous EWs were discovered by several deep photometric surveys and there were about 40 785 EW-type binary systems listed in the international variable star index(VSX)by 2017 March 13.7938 of them were observed with LAMOST by 2016 November 30 and their spectral types were identified.Stellar atmospheric parameters of 5363 EW-type binary stars were determined based on good spectroscopic observations.In the paper,those EWs are cataloged and their properties are analyzed.The distributions of orbital period(P),effective temperature(T),gravitational acceleration(log(g)),metallicity([Fe/H])and radial velocity(RV)are presented for these observed EW-type systems.It is shown that about 80.6% of sample stars have metallicity below zero,indicating that EW-type systems are old stellar populations.This is in agreement with the conclusion that EW binaries are formed from moderately close binaries through angular momentum loss via magnetic braking that takes a few hundred million to a few billion years.The unusually high metallicities of a few percent of EWs may be caused by contamination of material from the evolution of unseen neutron stars or black holes in the systems.The correlations between orbital period and effective temperature,gravitational acceleration and metallicity are presented and their scatters are mainly caused by(i)the presence of third bodies and(ii)sometimes wrongly determined periods.It is shown that some EWs contain evolved component stars and the physical properties of EWs mainly depend on their orbital periods.It is found that extremely short-period EWs may be older than their long-period cousins because they have lower metallicities.This reveals that they have a longer timescale of pre-contact evolution and their formation and evolution are mainly driven by angular momentum loss via magnetic braking.

A catalog of early-type emission-line stars and Hα line profiles from LAMOST DR2

We present a catalog including 11 204 spectra of 10436 early-type emission-line stars from LAMOST DR2, among which 9752 early-type emission-line spectra are newly discovered. For these early- type emissionqine stars, we discuss the morphological and physical properties of their low-resolution spectra. In this spectral sample, the Hα emission profiles display a wide variety of shapes. Based on the Hα line profiles, these spectra are categorized into five distinct classes： single-peak emission, single-peak emission in absorption, double-peak emission, double-peak emission in absorption, and P-Cygni profiles. To better understand what causes the Hc~ line profiles, we divide these objects into four types from the perspective of physical classification, which include classical Be stars, Herbig Ae/Be stars, close binaries and spectra contaminated by H n regions. The majority of Herbig Ae/Be stars and classical Be stars are identified and separated using a （H-K, K-W1） color-color diagram. We also discuss 31 binary systems that are listed in the SIMBAD on-line catalog and identify 3600 spectra contaminated by HⅡ regions after cross-matching with positions in the Dubout-Crillon catalog. A statistical analysis of line profiles versus classifications is then conducted in order to understand the distribution of Hc~ profiles for each type in our sample. Finally, we also provide a table of 172 spectra with Fen emission lines and roughly calculate stellar wind velocities for seven spectra with P-Cygni profiles.

NSVS 2569022:a peculiar binary among W UMa stars with extremely small mass ratios

Photometric observations of the W UMa binary NSVS 2569022 are presented. The light curve solution reveals that both components are of F spectral type （temperatures T1 = T2 = 6100 K）. NSVS 2569022 undergoes a total eclipse of W subtype and the mass ratio is well-determined. Its extremely small value of only 0.077 implies that the target will probably experience instability and a possible merger. This value ranks NSVS 2569022 in sixth place among binaries with the smallest mass ratio. Based on an empirical relation of ＂period - total mass＂ for low mass-ratio binaries, we estimate the global parameters of NSVS 2569022： masses M1 = 1.17 M⊙ and M2 = 0.09 M⊙; radii R1 = 1.19 R⊙ and R2 = 0.38 R⊙; luminosities L1 = 1.73 L⊙ and L2 = 0.17L⊙. An analysis of the characteristics of binaries with extremely low-mass ratios is made. NSVS 2569022 turns out to be a peculiar binary among W UMa stars with extremely small mass ratios due to its unexpectedly small fill-out factor of only 0.014 （slightly overcontact configuration）.

Photometric investigation and orbital period analyses of the W UMa binaries FP Lyn, FV CVn and V354 UMa

New light curves and photometric solutions of FP Lyn,FV CVn and V354 UMa are presented.We found that these three systems are W-subtype shallow contact binaries.In addition,it is obvious that the light curves of FP Lyn and V354 UMa are asymmetric.Therefore,a hot spot was added on the primary star of FP Lyn and a dark spot was added on the secondary star of V354 UMa.At the same time,we added a third light to the photometric solution of FP Lyn for the final result.The obtained mass ratios and fill-out factors are q = 1.153 and f = 13.4% for FP Lyn,q = 1.075 and f = 4.6% for FV CVn,and q = 3.623 and f = 10.7% for V354 UMa respectively.The investigations of orbital period for these three systems indicate that the periods are variable.FP Lyn and V354 UMa were discovered to have secularly increasing components with rates of dp/dt = 4.19 ×10^-7 dyr^-1 and dp/dt = 7.70 ×10^-7 dyr^-1 respectively,which are feasibly caused by conservative mass transfer from the less massive component to the more massive component.In addition,some variable components were discovered for FV CVn,including a rate of dp/dt =-1.13 ×10^-6 dyr^-1 accompanied by a cyclic oscillation with amplitude and period of 0.0069 d and 10.65 yr respectively.The most likely explanation for the long-term decrease is angular momentum loss.The existence of an additional star is the most plausible explanation for the periodic variation.

The W UMa binaries USNO-A2.0 1350-17365531,V471 Cas,V479 Lac and V560 Lac:light curve solutions and global parameters based on Gaia distances

We present photometric observations in Sloan filters g′, i′of the eclipsing W UMa stars USNOA2.0 1350-17365531, V471 Cas, V479 Lac and V560 Lac. The sinusoidal-like O-C diagram of V471 Cas indicates the presence of a third body with mass 0.12 M_⊙(a red dwarf) at distance 897 R_⊙. The O-C diagram of V479 Lac reveals a period decrease of d P/dt =-1.69 × 10-6d yr-1. The results of the light curve solutions are:(i) the targets are overcontact binaries with small fill-out factors;(ii) their components are F–K stars, comparable in size, whose temperature differences are below 80 K;(iii) all targets undergo partial eclipses and to limit the possible mass ratios we carried out two-step q-search analysis. The target global parameters(luminosities, radii, masses) were obtained on the basis of their Gaia distances and the results of our light curve solutions. The obtained total mass of V560 Lac turns out to be smaller than the lower mass limit for presently known W UMa binaries of 1.0-1.2 M_⊙, i.e. this target is a peculiar overcontact system.

Contact binaries at different evolutionary stages

Contact binaries consist of two strongly interacting component stars where they are filling their critical Roche lobes and sharing a common envelope.Most of them are main-sequence stars,but some of them are post main-sequence systems.They are good astrophysical laboratories for studying several problems such as the merging of binary stars,evolution of the common envelope,the origin of luminous red nova outbursts and the formation of rapidly rotating single stars with possible planetary systems.A large number of contact binary candidates were detected by several photometric surveys around the world and many of them were observed by the LAMOST spectroscopic survey.Based on follow-up observations,the evolutionary states and geometrical structures of some systems were understood well.In this review,we will introduce and catalog new stellar atmospheric parameters(i.e.,the effective temperature(Teff),the gravitational acceleration(log(g)),metallicity([Fe/H])and radial velocity(Vr))for 9149 EW-type contact binaries that were obtained based on low-and medium-resolution spectroscopic surveys of LAMOST.Then we will focus on several groups of contact binary stars,i.e.,marginal contact binary systems,deep and low-mass ratio contact binary stars,binary systems below the short-period limit of contact binaries and evolved contact binaries.Marginal contact binaries are at the beginning of the contact stage,while deep and low-mass ratio contact binary stars are at the final evolutionary stage of tidally locked binaries.Several statistical relations including the period-temperature relation are determined well by applying LAMOST data and their formation and evolutionary states are reviewed.The period-color relation of M-type binaries reveals that there are contact binaries below the short-period limit.Searching for and investigating contact binaries near and below this limit will help us to understand the formation of contact binary systems and a new prediction for the short-period limit is about 0.15 d.Some evolved contact binaries were detected by the LAMOST survey where both components are sub-giants or giants.They provide a good opportunity to investigate evolution of the common envelope and are the progenitors of luminous red novae like V1309 Sco.

Physical properties of γ Doradus pulsating stars and their relationship with long-period δ Scuti variables

We present LAMOST data on 168 γ Doradus(γ Dor) pulsating stars including stellar atmospheric parameters of 137 variables and spectral types for all of the samples. The distributions of period(P), temperature(T), gravitational acceleration(log(g)) and metallicity [Fe/H] are shown. It is found that most γ Dor variables are main-sequence stars with early F spectral types and temperatures from 6880 K to7280 K. They are slightly more metal poor than the Sun with a metallicity range from-0.4 to 0. On the H-R and log g-T diagrams, both the γ Dor and δ Scuti(δ Sct) stars occupy in the same region and some are beyond the borders predicted by current stellar pulsation theories. It is discovered that the physical properties of γ Dor stars are similar to those of long-period δ Sct(P > 0.3 d) stars. The stellar atmospheric parameters are all correlated with the pulsation period for short-period δ Sct variables(P < 0.3 d), but there are no such relations for γ Dor or long-period δ Sct stars. These results reveal that γ Dor and long-period δ Sct are the same group of pulsating stars and they are different from short-period δ Sct variables. Meanwhile, 33γ Dor stars are identified as candidates of binary or multiple systems.