Source author record

Xiang-Dong Li

Xiang-Dong Li appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

80works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

80 published item(s)

preprint2026arXiv

Formation of Recycled Pulsars in Common Envelope Binaries

We present a systematic study of the evolution of low- and intermediate-mass X-ray binaries (L/IMXBs) consisting of a $1.4\,M_{\odot}$ neutron star (NS) and a donor star of mass $1-8\,M_{\odot}$. Using grids of detailed MESA simulations, we show that for donor masses of $2-8\,M_{\odot}$, mass transfer may be dynamically unstable, leading to a common envelope (CE) phase. By adopting CE ejection efficiencies in the range $α_{\rm CE} = 0.3-3.0$, we find that post-CE binaries frequently experience a CE decoupling phase (CEDP), which plays a critical role in determining their final orbital and compositional properties. Systems with initial donor masses $\gtrsim 3.5\,M_{\odot}$ predominantly evolve into NS binaries with carbon-oxygen or oxygen-neon white dwarfs (WDs) with masses between $0.5\,M_{\odot}$ and $1.4\,M_{\odot}$. Comparison with the observed population of binary pulsars with a WD companion shows better agreement with higher CE ejection efficiencies ($α_{\rm CE} = 3.0$). Furthermore, we demonstrate that NSs can accrete a sufficient amount of matter ($\gtrsim 0.01\,M_{\odot}$) during the CEDP and subsequent Case BA/BB/BC mass transfer phases to be effectively recycled into millisecond pulsars. We identify two distinct evolutionary channels capable of reproducing the observed characteristics of the millisecond pulsar PSR J1928+1815 with a helium-star companion. Our results highlight the importance of the CEDP in the formation of recycled pulsars and provide constraints on the CE ejection efficiency during binary evolution.

preprint2026arXiv

How Beaming Shapes the Demographics of Ultraluminous X-ray Sources?

Ultraluminous X-ray sources (ULXs) are off-nuclear compact objects with apparent luminosities above 10^39 erg/s, often exceeding the Eddington limit for stellar-mass black holes. Beaming is a commonly invoked mechanism to explain their extreme brightness, and the dependence of the beaming factor on accretion rate is a critical parameter. In this work, we investigate how different beaming prescriptions affect the predicted properties of ULX populations. Using binary population synthesis, we construct synthetic X-ray luminosity functions (XLFs) for both classical and log-modified beaming models at solar and sub-solar metallicities. The classical model predicts a larger intrinsic number of bright ULXs, but strong beaming reduces their observable fraction, resulting in fewer visible ULXs compared to the log-modified model. The log-modified prescription yields a shallower slope at high-luminosity, aligning better with observed XLFs, and increases the fraction of observable neutron star ULXs above 10^39 erg/s. These results underscore the significant role of the beaming law in shaping ULX statistical distributions and assessing neutron star contributions to the population.

preprint2026arXiv

Pulsar Gleaners: Discovery of 19 Pulsars in FAST Archival Data at $|b|<5°$ and Decl.$<-5°$

We report the discovery of 19 new pulsars identified from archival observations of the Five-hundred-meter Aperture Spherical radio Telescope (FAST) within Galactic latitudes $|b|<5°$ and declinations ${\rm Decl.}<-5°$. The dataset was recorded using FAST's $L$-band 19-beam receiver and covered $\sim 3.6~{\rm deg}^2$ with a cumulative integration time of $\sim$ 500 hr and a total raw data volume of $\sim$ 700 TB. Our search employed fast Fourier transform (FFT)-based and fast folding algorithm (FFA)-based periodic searches, and the single-pulse search. These new pulsars have spin periods range from 0.03 to 5.54 s. Two have periods under 0.1 s, suggesting they are likely young pulsars or mildly recycled pulsars. Four pulsars exhibit dispersion measures (DMs) exceeding $1000~{\rm pc~cm^{-3}}$ with PSR J1839$-$0558t having the highest value in our sample at $\sim 1271~{\rm pc~cm^{-3}}$, providing valuable samples for pulsar studies in the high-DM regime. Two rotating radio transients, PSRs J1836$-$0552t and J1847$-$0624t, were detected by FFA and single-pulse searches but failed with the FFT-based searches. In addition, three faint pulsars that were also missed by FFT-based searches were successfully detected using FFA. These discoveries demonstrate the critical role of FFA in uncovering faint, long-period, and sporadic pulsars, and highlight the significant potential of FAST archival data, especially when combined with longer integration times and complementary search techniques, to reveal rare and weak pulsar populations.

preprint2022arXiv

Formation of Mass-gap Black Holes from Neutron Star X-ray Binaries with Super-Eddington Accretion

Electromagnetic and gravitational wave observations indicate that there is dearth of compact objects with mass $\sim 2.5-5~{\rm M_\odot}$. This so-called "mass gap" may be linked to the supernova explosion mechanisms that produce neutron stars (NSs) and black holes (BHs). However, the existence of a few mass-gap compact objects, some of which have been confirmed to be BHs, poses a challenge to the traditional theory of black hole formation. In this work we investigate the possible formation channel of BHs from accretion-induced collapse (AIC) of NSs in X-ray binaries. In particular, we consider the influence of super-Eddington accretion of NSs. Recent observations of ultraluminous X-ray pulsars suggest that their apparent luminosities may reflect the true accretion luminosities of the accreting NSs, even exceeding the Eddington limit by a factor of $\gtrsim 100$. Thus, NSs accreting at a super-Eddington accretion rate may rapidly grow into BHs in intermediate/low-mass X-ray binaries. Based on the super-Eddington accretion disk models, we have investigated the evolution of NSs in intermediate/low-mass X-ray binaries by combining binary population synthesis and detailed stellar evolutionary calculations. We show that super-Eddington accretion plays a critical role in mass growth of NSs, and the final masses of the descendant BHs are heavily dependent on the NS magnetic fields, the metallicity of the donor star, and the bifurcation period of the binaries. AIC of NSs may account for some of the observed mass-gap BHs like GRO J0422+32. We also present the parameter distributions of the potential mass-gap BHs in a Milky Way-like galaxy, and point out that future space-based gravitational wave observations may provide important test of or constraints on the formation of mass-gap BHs from the AIC channel.

preprint2022arXiv

Impact of Asymmetrical Mass Ejection from Proto-White Dwarfs on the Properties of Binary Millisecond Pulsars

The standard formation theory of binary millisecond pulsars (BMSPs) predicts efficient orbital circularization due to tidal interaction during the previous mass transfer phase. Therefore, BMSPs are expected to have a circular orbit. However, the discovery of several eccentric BMSPs (eBMSPs) with a white dwarf (WD) companion has challenged this picture. In particular, recent observation reveals that the spin angular momentum of the eBMSP J0955$-$6150 is tilted at an angle $>4.8^{\rm \degree}$ from the orbital angular momentum. This is the first time that a tilt angle is deduced for eBMSPs, which provides an important clue to their formation mechanism. Both the orbital eccentricity and tilt angle could be qualitatively accounted for by asymmetrical mass ejection during thermonuclear flashes from proto-WDs (so-called the thermonuclear rocket model), but detailed studies are still lacking. In this paper, we simulate the impact of the kick caused by asymmetrical mass ejection on the properties of BMSPs. We find that the thermonuclear rocket model can potentially explain the observational characteristics of both eBMSPs and normal BMSPs under reasonable input parameters. In addition, our results predict a wide range of the orbital period (from less than one day to more than several hundred days) for eBMSPs, which can be tested by future observations.

preprint2022arXiv

Stable mass transfer can explain massive binary black hole mergers with a high spin component

Recent gravitational wave observations showed that binary black hole (BBH) mergers with massive components are more likely to have high effective spins. In the model of isolated binary evolution, BH spins mainly originate from the angular momenta of the collapsing cores before BH formation. Both observations and theories indicate that BHs tend to possess relatively low spins, the origin of fast-spinning BHs remains a puzzle. We investigate an alternative process that stable Case A mass transfer may significantly increase BH spins during the evolution of massive BH binaries. We present detailed binary evolution calculations and find that this process can explain observed high spins of some massive BBH mergers under the assumption of mildly super-Eddington accretion.

preprint2021arXiv

A Study of Magnetized White Dwarf + Helium Star Binary Evolution to Type Ia Supernovae

The white dwarf (WD) + helium (He) star binary channel plays an important role in the single degenerate scenario for the progenitors of type Ia supernovae (SNe Ia). Previous studies on the WD + main sequence star evolution have shown that the magnetic fields of WDs may significantly influence their accretion and nuclear burning processes. In this work we focus on the evolution of magnetized WD + He star binaries with detailed stellar evolution and binary population synthesis (BPS) calculations. In the case of magnetized WDs, the magnetic fields may disrupt the inner regions of the accretion disk, funnel the accretion flow onto the polar caps, and even confine helium burning within the caps.We find that, for WDs with sufficiently strong magnetic fields, the parameter space of the potential SN Ia progenitor systems shrinks toward shorter orbital periods and lower donor masses compared with that in the non-magnetized WD case. The reason is that the magnetic confinement usually works with relatively high mass transfer rates, which can trigger strong wind mass loss from the WD, thus limiting the He-rich mass accumulation efficiency. The surviving companion stars are likely of low-mass at the moment of the SN explosions, which can be regarded as a possible explanation for the non-detection of surviving companions after the SNe or inside the SN remnants. However, the corresponding birthrate of Galactic SNe Ia in our high-magnetic models is estimated to be ~(0.08-0.13) * 10^{-3} yr^{-1}( ~0.17-0.28 * 10^{-3}yr^{-1} for the non-magnetic models), significantly lower than the observed Galactic SN Ia birthrate.

preprint2021arXiv

Are there magnetars in high-mass X-ray binaries?

Magnetars form a special population of neutron stars with strong magnetic fields and long spin periods. About 30 magnetars and magnetar candidates known currently are probably isolated. But the possibility that magnetars are in binaries hasn't been excluded. In this work, we perform spin evolution of neutron stars with different magnetic fields in wind-fed high-mass X-ray binaries and compare the spin period distribution with observations, aiming to find magnetars in binaries. Our simulation shows that some of the neutron stars, which have long spin periods or in wide-separation systems, need strong magnetic fields to explain their spin evolution. This implies that there are probably magnetars in high-mass X-ray binaries. Moreover, this can further provide a theoretical basis for some unclear astronomical phenomena, such as the possible origin of periodic fast radio bursts from magnetars in binary systems.

preprint2020arXiv

A peculiar cyclotron line near 16 keV detected in the 2015 outburst of 4U 0115+63?

In 2015 October, the Be/X-ray binary 4U 0115+63 underwent a type II outburst, reaching an X-ray luminosity of $\sim$ {\DP{38}\es}. During the outburst, Nuclear Spectroscopic Telescope Array ({\it NuSTAR}) performed two Target of Opportunity observations. Using the broad-band spectra from {\it NuSTAR} (3$-$79 keV), we have detected multiple cyclotron lines of the source, i.e., $\sim$ 12, 16, 22 and 33/35 keV. Obviously, the 16 keV line is not a harmonic component of the 12 keV line. As described by the phase-dependent equivalent widths of these cyclotron lines, the 16 keV and 12 keV lines are two different fundamental lines. In our work, we apply the two-poles cyclotron line model to the observation, i.e., the two line sets are formed at the same altitude ($\sim$ 0.2 km over the NS surface) of different magnetic poles, with $\sim$ {1.1\TDP{12}} and {1.4\TDP{12}} G in two poles, respectively.

preprint2020arXiv

On the formation of PSR J1640+2224: a neutron star born massive?

PSR J1640+2224 is a binary millisecond pulsar (BMSP) with a white dwarf (WD) companion. Recent observations indicate that the WD is very likely to be a $\sim 0.7\,M_{\odot}$ CO WD. Thus the BMSP should have evolved from an intermediate-mass X-ray binary (IMXB). However, previous investigations on IMXB evolution predict that the orbital periods of the resultant BMSPs are generally $<40$ days, in contrast with the 175 day orbital period of PSR J1640+2224. In this paper, we explore the influence of the mass of the neutron star (NS) and the chemical compositions of the companion star on the formation of BMSPs. Our results show that, the final orbital period becomes longer with increasing NS mass, and the WD mass becomes larger with decreasing metallicity. In particular, to reproduce the properties of PSR J1640+2224, the NS was likely born massive ($>2.0\,M_{\odot}$).

preprint2020arXiv

On the Formation of SS433

SS433 is an extraordinary X-ray binary which is ejecting bipolar jets with $26\%$ of the speed of light. Associated with the supernova-like shell W50, SS433 is also probably one of the youngest X-ray binaries with an age of $\lesssim 10^5$ yr. However, the masses of the two components in SS433 and even the nature of the compact object are still under debate. In this work, assuming that the compact object is a black hole, we employ a binary population synthesis method to study the formation of SS433. We use previous estimates of the age of W50 and the duration of the jet activity to constrain the evolutionary history. Our calculations suggest that SS433 likely harboured a Hertzsprung gap star at the beginning of the current Roche-lobe overflow phase. The masses of the black hole and the optical/donor star in SS433 suggested by the simulations are around $8\ M_{\odot}$ and $24\ M_{\odot}$, respectively. Future measurement of the donor mass and radius can help infer not only the orgin of the binary but also the nature of the nebula W50.

preprint2020arXiv

On the Renyi entropy power and the Gagliardo-Nirenberg-Sobolev inequality on Riemannian manifolds

In this paper, we prove the concavity of the Renyi entropy power for nonlinear diffusion equation (NLDE) associated with the Laplacian and the Witten Laplacian on compact Riemannian manifolds with non-negative Ricci curvature or $CD(0,m)$-condition and on compact manifolds equipped with time dependent metrics and potentials. Our results can be regarded as natural extensions of a result due to Savaré and Toscani \cite{ST} on the concavity of the Renyi entropy for NLDE on Euclidean spaces. Moreover, we prove that the rigidity models for the Renyi entropy power are the Einstein or quasi-Einstein manifolds and a special $(K,m)$-Ricci flow with Hessian solitons. Inspired by Lu-Ni-Vazquez-Villani \cite{LNVV}, we prove the Aronson-Benilan estimates for NLDE on compact Riemannian manifolds with $CD(0,m)$-condition. We also prove the NIW formula which indicates an intrinsic relationship between the second order derivative of the Renyi entropy power $N_p$, the $p$-th Fisher information $I_p$ and the time derivative of the $W$-entropy associated with NLDE. Finally, we prove the entropy isoperimetric inequality for the Renyi entropy power and the Gagliardo-Nirenberg-Sobolev inequality on complete Riemannian manifolds with non-negative Ricci curvature or $CD(0, m)$-condition and maximal volume growth condition.

preprint2020arXiv

Population Synthesis of Black Hole X-Ray Binaries

We present a systematic study of the X-ray binaries (XRBs) containing a black hole (BH) and a non-degenerate companion, in which mass transfer takes place via either capturing the companion's wind or Roche lobe overflow (RLO). As shown in our previous work that focusing on the formation and evolution of detached BH binaries, our assumed models relevant to BH's progenitors predicted significantly different binary properties \citep{sl19}. In this paper, we further follow the evolutionary paths of the BH systems that appearing as XRBs. By use of both binary population synthesis and detailed binary evolution calculations, we can obtain the potential population of BH XRBs. Distributions at the current epoch of various binary parameters have been computed. The observed sample of wind-fed XRBs can be well reproduced under assumption of all our models. The wind-fed XRBs are expected to be so rare ($ \lesssim 100 $) that only a couple of such systems have been detected. Comparison of known RLO XRBs with the calculated distributions of various binary parameters indicates that only the models assuming relatively small masses for BH progenitors can roughly match the observations. Accordingly we estimate that there are hundreds of RLO XRBs in the Milky Way, of which the majority are low-mass XRBs. The RLO systems may become ultraluminous X-ray sources (ULXs) if the BH accretes at a very high rate, we expect that about a dozen ULXs with a BH accretor may exist in a Milky Way$-$like galaxy.

preprint2020arXiv

Population Synthesis of Galactic Be-star Binaries with A Helium-star Companion

LB-1 was originally suggested to harbour a very massive ($ \sim 70 M_\odot $) black hole, but was recently suggested to be a post-mass transfer binary containing a Be star and a helium (He) star. In this paper, we use the binary population synthesis method to simulate the potential population of the Be$ - $He binaries in the Milky Way. Mass transfer process during the progenitor binary evolution plays a vital role in determining the possible properties of the Be$ - $He binary population. By constructing a range of physical models with significantly different mass-transfer efficiencies, we obtain the predicted distributions at the current epoch of the component masses and the orbital periods for the Be$ - $He binaries. In particular, we show that, LB-1 very likely has evolved through non-conservative mass transfer if it is indeed a Be$ - $He system. We estimate that there are more than $ 10^{3} $ Be$ - $He binaries with V-band apparent magnitudes brighter than LB-1.

preprint2019arXiv

A Population of Neutron Star Ultraluminous X-ray Sources with A Helium Star Companion

It was recently proposed that a significant fraction of ultraluminous X-ray sources (ULXs) actually host a neutron star (NS) accretor. We have performed a systematic study on the NS ULX population in Milky Way-like galaxies, by combining binary population synthesis and detailed stellar evolution calculations. Besides a normal star, the ULX donor can be a helium star (the hydrogen envelope of its progenitor star was stripped during previous common envelope evolution) if the NS is accreting at a super-Eddington rate via Roche lobe overflow. We find that the NS$-$helium star binaries can significantly contribute the ULX population, with the overall number of about several in a Milky Way-like galaxy. Our calculations show that such ULXs are generally close systems with orbital period distribution peaked at $ \sim 0.1 $ day (with a tail up to $ \sim100 $ days), and the helium stars have relatively low masses distributing with a maximum probability at $ \sim 1M_{\odot} $.

preprint2019arXiv

A strange star scenario for the formation of isolated millisecond pulsars

According to the recycling model, neutron stars in low-mass X-ray binaries were spun up to millisecond pulsars (MSPs), which indicates that all MSPs in the Galactic plane ought to be harbored in binaries. However, about $20\%$ Galactic field MSPs are found to be solitary. To interpret this problem, we assume that the accreting neutron star in binaries may collapse and become a strange star when it reaches some critical mass limit. Mass loss and a weak kick induced by asymmetric collapse during the phase transition (PT) from neutron star to strange star can result in isolated MSPs. In this work, we use a population-synthesis code to examine the PT model. The simulated results show that a kick velocity of $\sim60~{\rm km~s}^{-1}$ can produce $\sim6\times10^3$ isolated MSPs and birth rate of $\sim6.6\times10^{-7} {\rm ~yr}^{-1}$ in the Galaxy, which is approximately in agreement with predictions from observations. For the purpose of comparisons with future observation, we also give the mass distributions of radio and X-ray binary MSPs, along with the delay time distribution.

preprint2019arXiv

Gravitational Waves from Newborn Accreting Millisecond Magnetars

Two accretion columns have been argued to form over the surface of a newborn millisecond magnetar for an extremely high accretion rate $\gtrsim1.8\times10^{-2}M_\odot\ {\rm s^{-1}}$ that may occur in the core-collapse of a massive star. In this paper, we investigate the characteristics of these accretion columns and their gravitational wave (GW) radiation. For a typical millisecond magnetar (surface magnetic field strength $B\sim10^{15}$ G and initial spin period $P\sim1$ ms), we find (1) its accretion columns are cooled via neutrinos and can reach a height $\sim1$ km over the stellar surface; (2) its column-induced characteristic GW strain is comparable to the sensitivities of the next generation ground-based GW detectors within a horizon $\sim1$ Mpc; (3) the magnetar can survive only a few tens of seconds; (4) during the survival timescale, the height of the accretion columns increases rapidly to the peak and subsequently decreases slowly; (5) the column mass, characteristic GW strain, and maximum GW luminosity have simultaneous peaks in a similar rise-fall evolution. In addition, we find that the magnetar's spin evolution is dominated by the column accretion torque. A possible association with failed supernova is also discussed.

preprint2019arXiv

On the rapid orbital expansion in the compact low-mass X-ray binary 2A 1822$-$371

The neutron star low-mass X-ray binary 2A 1822$-$371 has an orbital period of 5.57 hr. Mass transfer in such short-period binaries is thought to be driven by magnetic braking with orbital shrinking. However, 2A 1822$-$371 shows a very rapid orbital expansion, implying that mass transfer occurs rapidly in this system. The accretion rate of the neutron star is observationally estimated to be higher than the Eddington limit, which is also hard to be explained by the standard magnetic braking mechanism. In this work, we construct a model to account for the peculiar properties of 2A 1822$-$371. We assume that the donor star possesses a relatively strong magnetic field, which is coupled with the stellar winds excited by the X-ray radiation from the neutron star. This would generate efficient angular momentum loss, leading to a high mass transfer rate and hence orbital expansion. We provide possible evolutionary tracks of 2A 1822$-$371 and study how the input parameters affect the results. The observational implications of the irradiation-driven mass loss are also briefly discussed in the context of evolution of low-mass X-ray binaries and millisecond pulsars.

preprint2016arXiv

A timing view of the heartbeat state of GRS 1915+105

We present a timing analysis of two Rossi X-ray Timing Explorer observations of the microquasar GRS 1915+105 during the heartbeat state. The phase-frequency-power maps show that the intermediate-frequency aperiodic X-ray variability weakens as the source softens in the slow rise phase, and when the quasi-periodic oscillation disappears in the rise phase of the pulse of the double-peaked class its sub-harmonic is still present with a hard phase lag. In the slow rise phase, the energy-frequency-power maps show that most of the aperiodic variability is produced in the corona, and may also induce the aperiodic variability observed at low energies from an accretion disk, which is further supported by the soft phase lag especially in the intermediate-frequency range (with a time delay up to 20 ms). In the rise phase of the pulse, the low-frequency aperiodic variability is enhanced significantly and there is a prominent hard lag (with a time delay up to 50 ms), indicating that the variability is induced by extension of the disk toward small radii as implied by the increase in flux and propagates into the corona. However, during the hard pulse of the double-peaked class, the variability shows no significant lag, which may be attributed to an optically thick corona. These timing results are generally consistent with the spectral results presented by Neilsen et al. (2011, 2012) which indicated that the slow rise phase corresponds to a local Eddington limit and the rise phase of the pulse corresponds to a radiation pressure instability in the disk.

preprint2016arXiv

Can the Subsonic Accretion Model Explain the Spin Period Distribution of Wind-fed X-ray Pulsars?

Neutron stars in high-mass X-ray binaries (HMXBs) generally accrete from the wind matter of their massive companion stars. Recently Shakura et al. (2012) suggested a subsonic accretion model for low-luminosity ($<4\times 10^{36}$ ergs$^{-1}$), wind-fed X-ray pulsars. To test the feasibility of this model, we investigate the spin period distribution of wind-fed X-ray pulsars with a supergiant companion star, using a population synthesis method. We find that the modeled distribution of supergiant HMXBs in the spin period - orbital period diagram is consistent with observations provided that the winds from the donor stars have relatively low terminal velocities ($\lesssim 1000$ kms$^{-1}$). The measured wind velocities in several supergiant HMXBs seem to favor this viewpoint. The predicted number ratio of wind-fed X-ray pulsars with persistent X-ray luminosities higher and lower than $4\times 10^{36}$ ergs$^{-1}$ is about $1:10$.

preprint2016arXiv

Evolution of Low-Mass X-ray Binaries: the Effect of Donor Evaporation

Millisecond pulsars (MSPs) are thought to originate from low-mass X-ray binaries (LMXBs). The discovery of eclipsing radio MSPs including redbacks and black widows indicates that evaporation of the donor star by the MSP's irradiation takes place during the LMXB evolution. In this work, we investigate the effect of donor evaporation on the secular evolution of LMXBs, considering different evaporation efficiencies and related angular momentum loss. We find that for widening LMXBs, the donor star leaves a less massive white dwarf than without evaporation; for contracting systems, evaporation can speed up the evolution, resulting in dynamically unstable mass transfer and possibly the formation of isolated MSPs.

preprint2016arXiv

Expanding molecular bubble surrounding Tycho's supernova remnant (SN 1572) observed with IRAM 30 m telescope: evidence for a single-degenerate progenitor

Whether the progenitors of SNe Ia are single-degenerate or double-degenerate white dwarf (WD) systems is a highly debated topic. To address the origin of the Type Ia Tycho's supernova remnant (SNR), SN 1572, we have carried out a $^{12}$CO J=2-1 mapping and a 3-mm line survey toward the remnant using the IRAM 30 m telescope. We show that Tycho is surrounded by a clumpy molecular bubble at a local standard of rest velocity of $\sim 61$ km s$^{-1}$ which expands at a speed of $\sim 4.5$ km s$^{-1}$ and has a mass of $\sim 220$ M$_\odot$ (at the distance of 2.5 kpc). Enhanced $^{12}$CO J=2-1 line emission relative to $^{12}$CO J=1-0 emission and possible line broadenings (in velocity range $-64$ to $-60$ km s$^{-1}$) are found at the northeastern boundary of the SNR, where the shell is deformed and decelerated. These features, combined with the morphological correspondence between the expanding molecular bubble and Tycho, suggest that the SNR is associated with the bubble at the velocity range $-66$ to $-57$ km s$^{-1}$. The most plausible origin for the expanding bubble is the fast outflow (with velocity of hundreds km s$^{-1}$) driven from the vicinity of a WD as it accreted matter from a nondegenerate companion star. The SNR has been expanding in the low-density wind-blown bubble, and the shock wave has just reached the molecular cavity wall. The SNR has been expanding in the low-density wind-blown bubble, and the shock wave has just reached the molecular cavity wall. This is the first unambiguous detection of an expanding bubble driven by the progenitor of a Type-Ia SNR, which constitutes evidence for a single-degenerate progenitor for this SN Ia.

preprint2016arXiv

Exploration of Spin-down Rate of Neutron Star in High Mass X-ray Binaries

We use the evolutionary population synthesis method to investigate the statistical properties of the wind-fed neutron star (NS) compact ($P_{\rm orb}<10$ days) high-mass X-ray binaries (HMXBs) in our Galaxy, based on different spin-down models. We find that the spin-down rate in the supersonic propeller phase given \textbf{by assuming that the surrounding material is treated as forming a quasi-static atmosphere} or \textbf{by assuming that the characteristic velocity of matter and the typical Alfv$\acute{e}$n velocity of material in the magnetospheric boundary layer are comparable to the sound speed in the external medium} is too low to produce the observed number of compact HMXBs. We also find that the models suggested \textbf{by assuming that the infalling material is ejected with the corotation velocity at the magnetospheric radius when the magnetospheric radius is larger than the corotation radius} and \textbf{by simple integration of the magnetic torque over the magnetosphere} with a larger spin-down rate than that given by \citet{dav81} or \citet{is75} can predict a reasonable number of observed wind-fed NS compact HMXBs. Our calculated results indicate that subsonic propeller phase may not exist at all by comparing with the observed particular distributions of wind-fed NS compact HMXBs in the $P_s-P_{orb}$ diagram. However, the spin-down rate suggested by \citet{wan85,dai06,jia05} and that given by \citet{dav73} both seem reasonable to produce the observed distribution of wind-fed NS compact HMXBs in the $P_s-P_{orb}$ diagram. We cannot find which spin-down rate seems more reasonable from our calculations.

preprint2016arXiv

Fe Line Diagnostics of Cataclysmic Variables and Galactic Ridge X-ray Emission

The properties of the Galactic Ridge X-ray Emission (GRXE) observed in the 2-10 keV band place fundamental constraints on various types of X-ray sources in the Milky Way. Although the primarily discrete origin of the emission is now well established, the responsible populations of these sources remain uncertain, especially at relatively low fluxes. To provide insights into this issue, we systematically characterize the Fe emission line properties of the candidate types of the sources in the solar neighborhood and compare them with those measured for the GRXE. Our source sample includes 6 symbiotic stars (SSs), 16 intermediate polars (IPs), 3 polars, 16 quiescent dwarf novae (DNe) and 4 active binaries (ABs). We find that the mean equivalent width ($EW_{6.7}$) of the 6.7-keV line and the mean 7.0/6.7-keV line ratio are $107\pm16.0$ eV and $0.71\pm 0.04$ for intermediate polars and $221\pm 135$ eV and $0.44\pm 0.14$ for polars, respectively, which are all substantially different from those ($490\pm15 $~eV and $0.2\pm 0.08$) for the GRXE. Instead, the GRXE values are better agreed by the $EW_{6.7}$ ($438\pm 84.6$~eV) and the ratio ($0.27\pm 0.06$) observed for the DNe. We further find that the $EW_{6.7}$ is strongly correlated with the 2-10-keV luminosity of the DNe, which can be characterized by the relation $EW_{6.7}=(438\pm95 {\rm~eV}) (L/10^{31}~{\rm~ergs~s^{-1} })^{(-0.31\pm0.15)}$. Accounting for this correlation, the agreement can be improved further, especially when the contributions from other classes sources to the GRXE are considered, which all have low $EW_{6.7}$ values. We conclude that the GRXE mostly consists of typically faint, but numerous DNe, plus ABs, while magnetic CVs are probably mainly the high-flux representatives of the responsible populations and dominate the GRXE only in harder energy bands.

preprint2016arXiv

Harnack inequalities and $W$-entropy formula for Witten Laplacian on Riemannian manifolds with $K$-super Perelman Ricci flow

In this paper, we prove logarithmic Sobolev inequalities and derive the Hamilton Harnack inequality for the heat semigroup of the Witten Laplacian on complete Riemannian manifolds equipped with $K$-super Perelman Ricci flow. We establish the $W$-entropy formula for the heat equation of the Witten Laplacian and prove a rigidity theorem on complete Riemannian manifolds satisfying the $CD(K, m)$ condition, and extend the $W$-entropy formula to time dependent Witten Laplacian on compact Riemannian manifolds with $(K, m)$-super Perelman Ricci flow, where $K\in \mathbb{R}$ and $m\in [n, \infty]$ are two constants. Finally, we prove the Li-Yau and the Li-Yau-Hamilton Harnack inequalities for positive solutions to the heat equation $\partial_t u=Lu$ associated to the time dependent Witten Laplacian on compact or complete manifolds equipped with variants of the $(K, m)$-super Ricci flow.

preprint2016arXiv

Monte Carlo population synthesis of post-common-envelope white dwarf binaries and type Ia supernova rate

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

preprint2016arXiv

Nonconservative Mass Transfer in Massive Binaries and the Formation of Wolf-Rayet+O Binaries

The mass transfer efficiency during the evolution of massive binaries is still uncertain. We model the mass transfer processes in a grid of binaries to investigate the formation of Wolf-Rayet+O (WR+O) binaries, taking into account two kinds of non-conservative mass transfer models: Model I with rotation-dependent mass accretion and Model II of half mass accretion. Generally the mass transfer in Model I is more inefficient, with the average efficiency in a range of $\sim0.2-0.7$ and $ \lesssim0.2 $ for Case A and Case B mass transfer, respectively. We present the parameter distributions for the descendant WR+O binaries. By comparing the modeled stellar mass distribution with the observed Galactic WR+O binaries, we find that highly non-conservative mass transfer is required.

preprint2016arXiv

On the Formation of Galactic Black Hole Low-Mass X-ray Binaries

Currently, there are 24 black hole (BH) X-ray binary systems that have been dynamically confirmed in the Galaxy. Most of them are low-mass X-ray binaries (LMXBs) comprised of a stellar-mass BH and a low-mass donor star. Although the formation of these systems has been extensively investigated, some crucial issues remain unresolved. The most noticeable one is that, the low-mass companion has difficulties in ejecting the tightly bound envelope of the massive primary during the spiral-in process. While initially intermediate-mass binaries are more likely to survive the common envelope (CE) evolution, the resultant BH LMXBs mismatch the observations. In this paper, we use both stellar evolution and binary population synthesis to study the evolutionary history of BH LMXBs. We test various assumptions and prescriptions for the supernova mechanisms that produce BHs, the binding energy parameter, the CE efficiency, and the initial mass distributions of the companion stars. We obtain the birthrate and the distributions of the donor mass, effective temperature and orbital period for the BH LMXBs in each case. By comparing the calculated results with the observations, we put useful constraints on the aforementioned parameters. In particular, we show that it is possible to form BH LMXBs with the standard CE scenario if most BHs are born through failed supernovae.

preprint2016arXiv

The Binding Energy Parameter for Common Envelope Evolution

The binding energy parameter $λ$ plays a vital role in common envelope evolution. Though it is well known that $λ$ takes different values for stars with different masses and varies during stellar evolution, it has been erroneously adopted as a constant in most of the population synthesis calculations. We have systematically calculated the values of $λ$ for stars of masses $1-60\,M_{\odot}$ by use of an updated stellar evolution code, taking into account contribution from both gravitational energy and internal energy to the binding energy of the envelope. We adopt the criterion for the core-envelope boundary advocated by \citet{Ivanova2011}. A new kind of $λ$ with the enthalpy prescription is also investigated. We present fitting formulae for the calculated values of various kinds of $λ$, which can be used in future population synthesis studies.

preprint2015arXiv

A Population of Ultraluminous X-ray Sources with An Accreting Neutron Star

Most ultraluminous X-ray sources (ULXs) are believed to be X-ray binary systems, but previous observational and theoretical studies tend to prefer a black hole rather than a neutron star accretor. The recent discovery of 1.37 s pulsations from the ULX M82 X-2 has established its nature as a magnetized neutron star. In this work we model the formation history of neutron star ULXs in an M82- or Milky Way-like galaxy, by use of both binary population synthesis and detailed binary evolution calculations. We find that the birthrate is around $10^{-4}\, \rm yr^{-1}$ for the incipient X-ray binaries in both cases. We demonstrate the distribution of the ULX population in the donor mass - orbital period plane. Our results suggest that, compared with black hole X-ray binaries, neutron star X-ray binaries may significantly contribute to the ULX population, and high-mass and intermediate-mass X-ray binaries dominate the neutron star ULX population in M82- and Milky Way-like galaxies, respectively.

preprint2015arXiv

A Strange Star Scenario for the Formation of Eccentric Millisecond Pulsar/Helium White Dwarf Binaries

According to the recycling scenario, millisecond pulsars (MSPs) have evolved from low-mass X-ray binaries (LMXBs). Their orbits are expected to be circular due to tidal interactions during the binary evolution, as observed in most of the binary MSPs. There are some peculiar systems that do not fit this picture. Three recent examples are PSRs J2234$+$06, J1946$+$3417 and J1950$+$2414, all of which are MSPs in eccentric orbits but with mass functions compatible with expected He white dwarf companions. It has been suggested these MSPs may have formed from delayed accretion-induced collapse of massive white dwarfs, or the eccentricity may be induced by dynamical interaction between the binary and a circumbinary disk. Assuming that the core density of accreting neutron stars in LMXBs may reach the density of quark deconfinement, which can lead to phase transition from neutron stars to strange quark stars, we show that the resultant MSPs are likely to have an eccentric orbit, due to the sudden loss of the gravitational mass of the neutron star during the transition. The eccentricities can be reproduced with a reasonable estimate of the mass loss. This scenario might also account for the formation of the youngest known X-ray binary Cir X$-$1, which also possesses a low-field compact star in an eccentric orbit.

preprint2015arXiv

Evolution of Transient Low-Mass X-ray Binaries to Redback Millisecond Pulsars

Redback millisecond pulsars (hereafter redbacks) are a sub-population of eclipsing millisecond pulsars in close binaries. The formation processes of these systems are not clear. The three pulsars showing transitions between rotation- and accretion-powered states belong to both redbacks and transient low-mass X-ray binaries (LMXBs), suggesting a possible evolutionary link between the them. Through binary evolution calculations, we show that the accretion disks in almost all LMXBs are subject to the thermal-viscous instability during certain evolutionary stages, and the parameter space for the disk instability covers the distribution of known redbacks in the orbital period - companion mass plane. We accordingly suggest that the abrupt reduction of the mass accretion rate during quiescence of transient LMXBs provides a plausible way to switch on the pulsar activity, leading to the formation of redbacks, if the neutron star has been spun up to be an energetic millisecond pulsar. We investigate the evolution of redbacks, taking into account the evaporation feedback, and discuss its possible influence on the formation of black widow millisecond pulsars.

preprint2015arXiv

Formation and Evolution of Galactic Intermediate/Low-Mass X-ray Binaries

We investigate the formation and evolutionary sequences of Galactic intermediate- and low-mass X-ray binaries (I/LMXBs) by combining binary population synthesis (BPS) and detailed stellar evolutionary calculations. Using an updated BPS code we compute the evolution of massive binaries that leads to the formation of incipient I/LMXBs, and present their distribution in the initial donor mass vs. initial orbital period diagram. We then follow the evolution of I/LMXBs until the formation of binary millisecond pulsars (BMSPs). We find that the birthrate of the I/LMXB population is in the range of $ 9\times10^{-6} - 3.4\times10^{-5} \, {\rm yr^{-1}}$, compatible with that of BMSPs which are thought to descend from I/LMXBs. We show that during the evolution of I/LMXBs they are likely to be observed as relatively compact binaries with orbital periods $ \lesssim $ 1 day and donor masses $\lesssim 0.3 M_{\odot}$. The resultant BMSPs have orbital periods ranging from less than 1 day to a few hundred days. These features are consistent with observations of LMXBs and BMSPs. We also confirm the discrepancies between theoretical predications and observations mentioned in the literature, that is, the theoretical average mass transfer rates ($ \sim 10^{-10} M_{\odot} $\,yr$^{-1}$) of LMXBs are considerably lower than observed, and the number of BMSPs with orbital periods $\sim 0.1-10$ day is severely underestimated. These discrepancies imply that something is missing in the modeling of LMXBs, which is likely to be related to the mechanisms of the orbital angular momentum loss.

preprint2015arXiv

Formation of Black Hole Low-Mass X-ray Binaries

The majority of known Galactic black holes reside in low-mass X-ray binaries. They are rare and fascinating objects, providing unique information on strong gravity, accretion disc physics, and stellar and binary evolution. There is no doubt that our understanding of the formation of black hole low-mass X-ray binaries has significantly advanced in the past decade. However, some key issues are still unresolved. In this paper we briefly summarize the observational clues and theoretical progress on the formation of black hole low-mass X-ray binaries.

preprint2015arXiv

Formation of the Double Neutron Star System PSR J1930$-$1852

The spin period (185 ms) and period derivative ($1.8\times10^{-17}\,\rm s\,s^{-1}$) of the double neutron star (DNS) system PSR J1930$-$1852 recently discovered indicate that the pulsar was mildly recycled through the process of Roche-lobe overflow. This system has the longest orbital period (45 days) of the known DNS systems, and can be formed from a helium star-NS binary if the initial mass of the helium star was $ \lesssim 4.0M_{\odot} $; otherwise the helium star would never fill its Roche-lobe \citep{t15}. At the moment of the supernova explosion, the mass of the helium star was $ \lesssim3.0M_{\odot} $. We find that the probability distribution of the velocity kick imparted to the new-born neutron star has a maximum at about $30 \,\rm km\,s^{-1}$ (and a tail up to $ 260 \,\rm km\,s^{-1}$), indicating that this NS most probably received a low kick velocity at birth.

preprint2015arXiv

Low braking index of PSR J1734-3333: an interaction between fall-back disk and magnetic field?

Recent timing observation reported that the radio pulsar PSR J1734 - 3333 with a rotating period $P=1.17~\rm s$ is slowing down with a period derivative $\dot{P}=2.28\times 10^{-12}\rm s\,s^{-1}$. Its derived braking index $n=0.9 \pm 0.2$ is the lowest value among young radio pulsars with the measured braking indices. In this Letter, we attempt to investigate the influence of the braking torque caused by the interaction between the fall-back disk and the strong magnetic field of the pulsar on the spin evolution of PSR J1734 - 3333. Analytical result show that this braking torque is obviously far more than that by magnetic dipole radiation for pulsars with spin period of $> 0.1$ s, and play an important role during the spin-down of the pulsars. Our simulated results indicate that, for some typical neutron star parameters, the braking index and the period derivative approximately in agreement with the measured value of PSR J1734 - 3333 if the material inflow rate in the fallback disk is $2 \times 10^{17} \rm g\,s^{-1}$. In addition, our scenario can account for the measured braking indices of four young pulsars. However, our predicted X-ray luminosity are 1 -2 order of magnitude higher than the observation. We proposed that this discrepancy may originate from the instability of fall-back disk.

preprint2015arXiv

On the evolution of a fossil disk around neutron stars originating from merging white dwarfs

Numerical simulations suggest that merging double white dwarfs (WDs) may produce a newborn neutron star surrounded by a fossil disk. We investigate the evolution of the fossil disk following the coalescence of double WDs. We demonstrate that the evolution can be mainly divided into four phases: the slim disk phase (with time $\lesssim$ 1 yr), the inner slim plus outer thin disk phase ($\sim 10-\DP{6}$ yr), the thin disk phase ($\sim \DP{2}-\DP{7}$ yr), and the inner advection-dominated accretion flow plus outer thin disk phase, given the initial disk mass $\sim 0.05-0.5\,M_{\sun}$ and the disk formation time $10^{-3}-1$ s. Considering possible wind mass loss from the disk, we present both analytic formulae and numerically calculated results for the disk evolution, which is sensitive to the condition that determines the location of the outer disk radius. The systems are shown to be very bright in X-rays in the early phase, but quickly become transient within $\lesssim$ 100 yr, with peak luminosities decreasing with time. We suggest that they might account for part of the very faint X-ray transients around the Galactic center region, which generally require a very low mass transfer rate.

preprint2015arXiv

On the Law of Large Numbers for the empirical measure process of Generalized Dyson Brownian motion

We study the generalized Dyson Brownian motion (GDBM) of an interacting $N$-particle system with logarithmic Coulomb interaction and general potential $V$. Under reasonable condition on $V$, we prove the existence and uniqueness of strong solution to SDE for GDBM. We then prove that the family of the empirical measures of GDBM is tight on $\mathcal {C}([0,T],\mathscr{P}(\mathbb{R}))$ and all the large $N$ limits satisfy a nonlinear McKean-Vlasov equation. Inspired by previous works due to Biane and Speicher, Carrillo, McCann and Villani, we prove that the McKean-Vlasov equation is indeed the gradient flow of the Voiculescu free entropy on the Wasserstein space of probability measures over $\mathbb{R}$. Using the optimal transportation theory, we prove that if $V"\geq K$ for some constant $K\in \mathbb{R}$, the McKean-Vlasov equation has a unique weak solution. This proves the Law of Large Numbers and the propagation of chaos for the empirical measures of GDBM. Finally, we prove the longtime convergence of the McKean-Vlasov equation for $C^2$-convex potentials $V$.

preprint2015arXiv

Orbital period decay of compact black hole x-ray binaries: the influence of circumbinary disks?

Recently, compact black hole X-ray binaries XTE J 1118+480 and A0620-00 have been reported to be experiencing a fast orbital period decay, which is two orders of magnitude higher than expected with gravitational wave radiation. Magnetic braking of an Ap/Bp star has been suggested to account for the period change when the surface magnetic field of the companion star $B_{\rm s}\ga 10^{4}$ G. However, our calculation indicates that anomalous magnetic braking cannot significantly contribute to the large orbital period decay rates observed in these two sources even if $B_{\rm s}\ga 10^{4}$ G. Observations have provided evidence that circumbinary disks around two compact black hole X-ray binaries may exist. Our analysis shows that, for some reasonable parameters, tidal torque between the circumbinary disk and the binary can efficiently extract the orbital angular momentum from the binary, and result in a large orbital period change rate. Based on the circumbinary disk model, we simulate the evolution of XTE J 1118+480 via a stellar evolution code. Our computations are approximatively in agreement with the observed data (the masses of two components, donor star radius, orbital period, and orbital period derivative). The mass transfer rate and circumbinary disk mass are obviously far greater than the inferred values from observations. Therefore, it seems that the circumbinary disk is unlikely to be the main cause of the rapid orbital decay observed in some compact black hole X-ray binaries.

preprint2015arXiv

Super-Strong Magnetic Fields of Neutron Stars in Be X-Ray Binaries Estimated with New Torque and Magnetosphere Models

We re-estimate the surface magnetic fields of neutron stars (NSs) in Be X-ray binaries (BeXBs) with different models of torque, improved beyond Klus et al. (2014). In particular, a new torque model (Dai \& Li 2006) is applied to three models of magnetosphere radius. Unlike the previous models, the new torque model does not lead to divergent results for any fastness parameter. The inferred surface magnetic fields of these NSs for the two compressed magnetosphere models are much higher than that for the uncompressed magnetosphere model. The new torque model using the compressed-magnetosphere radius (Shi, Zhang \& Li 2014) leads to unique solutions near spin equilibrium in all cases, unlike other models that usually give two branches of solutions. Although our conclusions are still affected by the simplistic assumptions about the magnetosphere radius calculations, we show several groups of possible surface magnetic field values with our new models when the interaction between the magnetosphere and the infalling accretion plasma is considered. The estimated surface magnetic fields for NSs BeXBs in the Large Magellanic Cloud, the Small Magellanic Cloud and the Milk Way are between the quantum critical field and the maximum "virial" value by the spin equilibrium condition.

preprint2015arXiv

The orbital period evolution of the supersoft X-ray source CAL 87

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

preprint2014arXiv

$W$-entropy formula for the Witten Laplacian on manifolds with time dependent metrics and potentials

In this paper, we develop a new approach to prove the $W$-entropy formula for the Witten Laplacian via warped product on Riemannian manifolds and give a natural geometric interpretation of a quantity appeared in the $W$-entropy formula. Then we prove the $W$-entropy formula for the Witten Laplacian on compact Riemannian manifolds with time dependent metrics and potentials, and derive the $W$-entropy formula for the backward heat equation associated with the Witten Laplacian on compact Riemannian manifolds equipped with Lott's modified Ricci flow. We also extend our results to complete Riemannian manifolds with negative $m$-dimensional Bakry-Emery Ricci curvature, and to compact Riemannian manifolds with $K$-super $m$-dimensional Bakry-Emery Ricci flow. As application, we prove that the optimal logarithmic Sobolev constant on compact manifolds equipped with the $K$-super $m$-dimensional Bakry-Emery Ricci flow is decreasing in time.

preprint2014arXiv

Common Envelope Mechanisms: Constraints from the X-ray Luminosity Function of High Mass X-ray Binaries

We use the measured X-ray luminosity function (XLF) of high-mass X-ray binaries (HMXBs) in nearby star-forming galaxies to constrain the common envelope (CE) mechanisms, which play a key role in governing the binary evolution. We find that the XLF can be reproduced quite closely under both CE mechanisms usually adopted, i.e., the $α_{\rm CE}$ formalism and the $γ$ algorithm, with a reasonable range of parameters considered. Provided that the parameter combination is the same, the $γ$ algorithm is likely to produce more HMXBs than the $α_{\rm CE}$ formalism, by a factor of up to $\sim$ 10. In the framework of the $α_{\rm CE}$ formalism, a high value of $α_{\rm CE}$ is required to fit the observed XLF, though it does not significantly affect the global number of the HMXB populations. We present the detailed components of the HMXB populations under the $γ$ algorithm and compare them with those in Zuo et al. and observations. We suggest the distinct observational properties, as well as period distributions of HMXBs, may provide further clues to discriminate between these two types of CE mechanisms.

preprint2014arXiv

Discovery of the transient magnetar 3XMM J185246.6+003317 near supernova remnant Kesteven 79 with XMM-Newton

We report the serendipitous discovery with XMM-Newton that 3XMM J185246.6+003317 is an 11.56 s X-ray pulsar located 1' away from the southern boundary of supernova remnant Kes 79. The spin-down rate of 3XMM J185246.6+003317 is $<1.1\times 10^{-13}$ s s$^{-1}$, which, together with the long period P=11.558714(2) s, indicates a dipolar surface magnetic field of $<3.6\times 10^{-13}$ G, a characteristic age of $>1.7$ Myr, and a spin-down luminosity of $<2.8\times 10^{30}$ erg s$^{-1}$. The X-ray spectrum of the source is best-fitted with a resonant Compton scattering model, and can be also adequately described by a blackbody model. The observations covering a seven month span from 2008 to 2009 show variations in the spectral properties of the source, with the luminosity decreasing from $2.7\times 10^{34}$ erg s$^{-1}$ to $4.6 \times 10^{33}$ erg s$^{-1}$, along with a decrease of the blackbody temperature from $kT\approx 0.8$ keV to $\approx0.6$ keV. The X-ray luminosity of the source is higher than its spin-down luminosity, ruling out rotation as a power source. The combined timing and spectral properties, the non-detection of any optical or infrared counterpart, together with the lack of detection of the source in archival X-ray data prior to the 2008 XMM-Newton observation, point to this source being a newly discovered transient low-B magnetar undergoing an outburst decay during the XMM-Newton observations. The non-detection by Chandra in 2001 sets an upper limit $4\times 10^{32}$ erg s$^{-1}$ to the quiescent luminosity of 3XMM J185246.6+003317. Its period is the longest among currently known transient magnetars. The foreground absorption toward 3XMM J185246.6+003317 is similar to that of Kes 79, suggesting a similar distance of $\sim$7.1 kpc.

preprint2014arXiv

Formation of Binary Millisecond Pulsars by Accretion-Induced Collapse of White Dwarfs under Wind-Driven Evolution

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

preprint2014arXiv

Formation of Millisecond Pulsars with Low-Mass Helium White Dwarf Companions in Very Compact Binaries

Binary millisecond pulsars (BMSPs) are thought to have evolved from low-mass X-ray binaries (LMXBs). If the mass transfer in LMXBs is driven by nuclear evolution of the donor star, the final orbital period is predicted to be well correlated with the mass of the white dwarf (WD), which is the degenerate He core of the donor. Here we show that this relation can be extended to very small WD mass ($\sim 0.14-0.17\,M_\odot$) and narrow orbital period (about a few hours), mainly depending on the metallicities of the donor stars. There is also discontinuity in the relation, which is due to the temporary contraction of the donor when the H-burning shell crosses the hydrogen discontinuity. BMSPs with low-mass He WD companions in very compact binaries can be accounted for if the progenitor binary experienced very late Case A mass transfer. The WD companion of PSR J1738+0333 is likely to evolve from a Pop~II star. For PSR J0348+0432, to explain its extreme compact orbit in the Roche lobe-decoupling phase, even lower metallicity ($Z=0.0001$) is required.

preprint2014arXiv

Hamilton's Harnack inequality and the $W$-entropy formula on complete Riemannian manifolds

In this paper, we prove Hamilton's Harnack inequality and the gradient estimates of the logarithmic heat kernel for the Witten Laplacian on complete Riemainnian manifolds. As applications, we prove the $W$-entropy formula for the Witten Laplacian on complete Riemannian manifolds, and prove a family of logarithmic Sobolev inequalities on complete Riemannian manifolds with natural geometric condition.

preprint2014arXiv

On the Common Envelope Efficiency

In this work, we try to use the apparent luminosity versus displacement (i.e., $L_{\rm X}$ vs. $R$) correlation of high mass X-ray binaries (HMXBs) to constrain the common envelope (CE) efficiency $α_{\rm CE}$, which is a key parameter affecting the evolution of the binary orbit during the CE phase. The major updates that crucial for the CE evolution include a variable $λ$ parameter and a new CE criterion for Hertzsprung gap donor stars, both of which are recently developed. We find that, within the framework of the standard energy formula for CE and core definition at mass $X=10$\%, a high value of $α_{\rm CE}$, i.e., around 0.8-1.0, is more preferable, while $α_{\rm CE}< \sim 0.4$ likely can not reconstruct the observed $L_{\rm X}$ vs. $R$ distribution. However due to an ambiguous definition for the core boundary in the literature, the used $λ$ here still carries almost two order of magnitude uncertainty, which may translate directly to the expected value of $α_{\rm CE}$. We present the detailed components of current HMXBs and their spatial offsets from star clusters, which may be further testified by future observations of HMXB populations in nearby star-forming galaxies.

preprint2014arXiv

On the formation of Be stars through binary interaction

Be stars are rapidly rotating B type stars. The origin of their rapid rotation is not certain, but binary interaction remains to be a possibility. In this work we investigate the formation of Be stars resulting from mass transfer in binaries in the Galaxy. We calculate the binary evolution with both stars evolving simultaneously and consider different possible mass accretion histories for the accretor. From the calculated results we obtain the critical mass ratios $q_{\rm cr}$ that determine the stability of mass transfer. We also numerically calculate the parameter $λ$ in common envelope evolution, and then incorporate both $q_{\rm cr}$ and $λ$ into the population synthesis calculations. We present the predicted numbers and characteristics of Be stars in binary systems with different types of companions, including helium stars, white dwarfs, neutron stars, and black holes. We find that in Be/neutron star binaries the Be stars can have a lower limit of mass $ \sim 8 M_{\odot}$ if they are formed by stable (i.e., without the occurrence of common envelope evolution) and nonconservative mass transfer. We demonstrate that isolated Be stars may originate from both mergers of two main-sequence stars and disrupted Be binaries during the supernova explosions of the primary stars, but mergers seem to play a much more important role. Finally the fraction of Be stars which have involved binary interactions in all B type stars can be as high as $ \sim 13%-30% $, implying that most of Be stars may result from binary interaction.

preprint2014arXiv

The Magneto Hydro Dynamical Model of KHz Quasi Periodic Oscillations in Neutron Star Low Mass X-ray Binaries (II)

We study the kilohertz quasi-periodic oscillations (kHz QPOs) in neutron star low mass X-ray binaries (LMXBs) with a new magnetohydrodynamics (MHD) model, in which the compressed magnetosphere is considered. The previous MHD model (Shi \& Li 2009) is re-examined and the relation between the frequencies of the kHz QPOs and the accretion rate in LMXBs is obtained. Our result agrees with the observations of six sources (4U 0614+09, 4U 1636--53, 4U 1608--52, 4U 1915--15, 4U 1728--34, XTE 1807--294) with measured spins. In this model the kHz QPOs originate from the MHD waves in the compressed magnetosphere. The single kHz QPOs and twin kHz QPOs are produced in two different parts of the accretion disk and the boundary is close to the corotation radius. The lower QPO frequency in a frequency-accretion rate diagram is cut off at low accretion rate and the twin kHz QPOs encounter a top ceiling at high accretion rate due to the restriction of innermost stable circular orbit.

preprint2013arXiv

From the Boltzmann $H$-theorem to Perelman's $W$-entropy formula for the Ricci flow

In 1870s, L. Boltzmann proved the famous $H$-theorem for the Boltzmann equation in the kinetic theory of gas and gave the statistical interpretation of the thermodynamic entropy. In 2002, G. Perelman introduced the notion of $W$-entropy and proved the $W$-entropy formula for the Ricci flow. This plays a crucial role in the proof of the no local collapsing theorem and in the final resolution of the Poincaré conjecture and Thurston's geometrization conjecture. In our previous paper \cite{Li11a}, the author gave a probabilistic interpretation of the $W$-entropy using the Boltzmann-Shannon-Nash entropy. In this paper, we make some further efforts for a better understanding of the mysterious $W$-entropy by comparing the $H$-theorem for the Boltzmann equation and the Perelman $W$-entropy formula for the Ricci flow. We also suggest a way to construct the "density of states" measure for which the Boltzmann $H$-entropy is exactly the $W$-entropy for the Ricci flow.

preprint2013arXiv

Generalized Dyson Brownian motion, McKean-Vlasov equation and eigenvalues of random matrices

Using Itô's calculus and the mass optimal transportation theory, we study the generalized Dyson Brownian motion (GDBM) and the associated McKean-Vlasov evolution equation with an external potential $V$. Under suitable condition on $V$, we prove the existence and uniqueness of strong solution to SDE for GDBM. Standard argument shows that the family of the process of empirical measures $L_N(t)$ of GDBM is tight and every accumulative point of $L_N(t)$ in the weak convergence topology is a weak solution of the associated McKean-Vlasov evolution equation, which can be realized as the gradient flow of the Voiculescu free entropy on the Wasserstein space over $\mathbb{R}$. Under the condition $V''\geq -K$ for some constant $K\geq 0$, we prove that the McKean-Vlasov equation has a unique solution $μ(t)$ and $L_N(t)$ converges weakly to $μ(t)$ as $N\rightarrow \infty$. For $C^2$ convex potentials, we prove that $μ(t)$ converges to the equilibrium measure $μ_V$ with respect to the $W_2$-Wasserstein distance on $\mathscr{P}_2(\mathbb{R})$ as $t\rightarrow \infty$. Under the uniform convexity or a modified uniform convexity condition on $V$, we prove that $μ(t)$ converges to $μ_V$ with respect to the $W_2$-Wasserstein distance on $\mathscr{P}_2(\mathbb{R})$ with an exponential rate as $t\rightarrow \infty$. Finally, we discuss the double-well potentials and raise some conjectures.

preprint2013arXiv

On the $L^p$-estimates for Beurling-Ahlfors and Riesz transforms on Riemannian manifolds

In our previous papers \cite{Li2008, Li2011}, we proved some martingale transform representation formulas for the Riesz transforms and the Beurling-Ahlfors transforms on complete Riemannian manifolds, and proved some explicit $L^p$-norm estimates for these operators on complete Riemannian manifolds with suitable curvature conditions. In this paper we correct a gap contained in \cite{Li2008, Li2011} and prove that the $L^p$-norm of the Riesz transforms $R_a(L)=\nabla(a-L)^{-1/2}$ can be explicitly bounded by $C(p^*-1)^{3/2}$ if $Ric+\nabla^2ϕ\geq -a$ for $a\geq 0$, and the $L^p$-norm of the Riesz transform $R_0(L)=\nabla(-L)^{-1/2}$ is bounded by $2(p^*-1)$ if $Ric+\nabla^2ϕ=0$. We also prove that the $L^p$-norm estimates for the Beurling-Ahlfors transforms obtained in \cite{Li2011} remain valid. Moreover, we prove the time reversal martingale transform representation formulas for the Riesz transforms and the Beurling-Ahlfors transforms on complete Riemannian manifolds.

preprint2013arXiv

On the $L^p$-estimates of Riesz transforms on forms over complete Riemanian manifolds

In our previous paper \cite{Li2010}, we proved a martingale transform representation formula for the Riesz transforms on forms over complete Riemannian manifolds, and proved some explicit $L^p$-norm estimates for the Riesz transforms on complete Riemannian manifolds with suitable curvature conditions. In this paper we correct a gap contained in \cite{Li2010} and prove that the main result obtained in \cite{Li2010} on the $L^p$-norm estimates for the Riesz transforms on forms remain valid. Moreover, we prove a time reversal martingale transform representation formula for the Riesz transforms on forms. Finally, we extend our approach and result to the Riesz transforms acting on Euclidean vector bundles over complete Riemannian manifolds with suitable curvature conditions.

preprint2013arXiv

On the formation of the peculiar low-mass X-ray binary IGR J17480$-$2446 in Terzan 5

IGR J17480$-$2446 is an accreting X-ray pulsar in a low-mass X-ray binary harbored in the Galactic globular cluster Terzan 5. Compared with other accreting millisecond pulsars, IGR J17480$-$2446 is peculiar in its low spin frequency (11 Hz), which suggests that it might be a mildly recycled neutron star at the very early phase of mass transfer. However, this model seems to be in contrast with the low field strength deduced from the kiloHertz quasi-periodic oscillations observed in IGR J17480$-$2446. Here we suggest an alternative interpretation, assuming that the current binary system was formed during an exchange encounter either between a binary (which contains a recycled neutron star) and the current donor, or between a binary and an isolated, recycled neutron star. In the resulting binary, the spin axis of the neutron star could be parallel or anti-parallel with the orbital axis. In the later case, the abnormally low frequency of IGR J17480$-$2446 may result from the spin-down to spin-up evolution of the neutron star. We also briefly discuss the possible observational implications of the pulsar in this scenario.

preprint2013arXiv

Population synthesis of young isolated neutron stars: the effect of fallback disk accretion and magnetic field evolution

The spin evolution of isolated neutron stars (NSs) is dominatd by their magnetic fields. The measured braking indices of young NSs show that the spin-down mechanism due to magnetic dipole radiation with constant magnetic fields is inadequate. Assuming that the NS magnetic field is buried by supernova fallback matter and re-emerges after accretion stops, we carry out Monte-Carlo simulation of the evolution of young NSs, and show that most of the pulsars have the braking indices ranging from -1 to 3. The results are compatible with the observational data of NSs associated with supernova remnants. They also suggest that the initial spin periods of NSs might occupy a relatively wide range.

preprint2013arXiv

Population synthesis on high-mass X-ray binaries: prospects and constraints from the universal X-ray luminosity function

Using an updated population synthesis code initially developed by Hurley et al. we modeled the synthetic X-ray binary (XRB) populations for direct comparison with the universal, featureless X-ray luminosity function (XLF) of high mass X-ray binaries (HMXBs) in star-forming galaxies. Our main goal is to use the universal XLF to constrain the model parameters, given the current knowledge of binary evolution. We find that the one-dimensional (1D) Maxwellian velocity dispersion of the natal kick can be constrained to be of the order of ~150km/s, supporting earlier findings that neutron stars (NSs) formed in binaries seem to receive significantly smaller natal kicks than the velocities of Galactic single pulsars would indicate. The super-Eddington accretion factor is further confirmed in the framework of stellar mass black holes (BHs), revealing the true origin of the most of the ultraluminous X-ray sources (ULXs) may indeed be the high-luminosity extension of ordinary HMXBs which harbor stellar-mass BHs rather than exotic intermediate-mass BHs or ones. We present the detail properties of the model-predicted present-day HMXBs, which may be investigated by future high-resolution X-ray and optical observations of sources in nearby star-forming galaxies.

preprint2013arXiv

Strange star equation of state fits the refined mass measurement of 12 pulsars and predicts their radii

There are three categories of stars whose masses have been found accurately in recent times: (1) two for which Shapiro delay is used which is possible due to GR light bending as the partner is heavy : PSR J1614-2230 and PSR J1903+0327 (2) six eclipsing stars for which numerical Roche Lobe geometry is used and (3) 3 stars for which spectroscopic methods are used and in fact for these three the mass and radii both are estimated. Motivated by large color (N_c) expansion using a modified Richardson potential, along with density dependent quark masses thereby allowing chiral symmetry restoration, we get compact strange stars fitting all the observed masses.

preprint2013arXiv

Wind-Driven Evolution of White Dwarf Binaries to Type Ia Supernovae

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

preprint2012arXiv

Could SXP 1062 be an Accreting Magnetar?

In this work we explore the possible evolutionary track of the neutron star in the newly discovered Be/X-ray binary SXP 1062, which is believed to be the first X-ray pulsar associated with a supernova remnant. Although no cyclotron feature has been detected to indicate the strength of the neutron star's magnetic field, we show that it may be $\ga 10^{14} $G. If so SXP 1062 may belong to the accreting magnetars in binary systems. We attempt to reconcile the short age and long spin period of the pulsar taking account of different initial parameters and spin-down mechanisms of the neutron star. Our calculated results show that, to spin down to a period $\sim 1000 $s within $10-40 $kyr requires efficient propeller mechanisms. In particular, the model for angular momentum loss under energy conservation seems to be ruled out.

preprint2012arXiv

Formation of Millisecond Pulsars from Intermediate- and Low-Mass X-ray Binaries

We present a systematic study of the evolution of intermediate- and low-mass X-ray binaries consisting of an accreting neutron star of mass $1.0-1.8 M_{\odot}$ and a donor star of mass $1.0-6.0 M_{\odot}$. In our calculations we take into account physical processes such as unstable disk accretion, radio ejection, bump-induced detachment, and outflow from the $L_{2}$ point. Comparing the calculated results with the observations of binary radio pulsars, we report the following results. (1) The allowed parameter space for forming binary pulsars in the initial orbital period - donor mass plane increases with increasing neutron star mass. This may help explain why some MSPs with orbital periods longer than $\sim 60$ days seem to have less massive white dwarfs than expected. Alternatively, some of these wide binary pulsars may be formed through mass transfer driven by planet/brown dwarf-involved common envelope evolution. (2) Some of the pulsars in compact binaries might have evolved from intermediate-mass X-ray binaries with anomalous magnetic braking. (3) The equilibrium spin periods of neutron stars in low-mass X-ray binaries are in general shorter than the observed spin periods of binary pulsars by more than one order of magnitude, suggesting that either the simple equilibrium spin model does not apply, or there are other mechanisms/processes spinning down the neutron stars.

preprint2012arXiv

On the first eigenvalue of the Witten-Laplacian and the diameter of compact shrinking Ricci solitons

We prove a lower bound estimate for the first non-zero eigenvalue of the Witten-Laplacian on compact Riemannian manifolds. As an application, we derive a lower bound estimate for the diameter of compact gradient shrinking Ricci solitons. Our results improve some previous estimates which were obtained by the first author and Y. Sano in [12], and by B. Andrews and L. Ni in [1]. Moreover, we extend the diameter estimate to compact self-similar shrinkers of mean curvature flow.

preprint2011arXiv

Angular Momentum Loss Mechanisms in Cataclysmic Variables below the Period Gap

Mass transfer in cataclysmic variables (CVs) is usually considered to be caused by angular momentum loss (AML) driven by magnetic braking and gravitational radiation (GR) above the period gap, and solely by GR below the period gap. The best-fit revised model of CV evolution recently by \citet{kbp11}, however, indicates that AML rate below the period gap is $2.47(\pm 0.22)$ times the GR rate, suggesting the existence of some other AML mechanisms. We consider several kinds of consequential AML mechanisms often invoked in the literature: isotropic wind from the accreting white dwarfs, outflows from the Langrangian points, and the formation of a circumbinary disk. We found that neither isotropic wind from the white dwarf nor outflow from the $L_1$ point can explain the extra AML rate, while ouflow from the $L_2$ point or a circumbinary disk can effectively extract the angular momentum provided that $\sim (15-45)%$ of the transferred mass is lost from the binary. A more promising mechanism is a circumbinary disk exerting gravitational torque on the binary. In this case the mass loss fraction can be as low as $\lesssim 10^{-3}$.

preprint2011arXiv

Can eccentric binary millisecond pulsars form by accretion induced collapse of white dwarfs?

Binary radio pulsars are generally believed to have been spun up to millisecond periods (i.e. recycling) via mass accretion from their donor stars, and they are the descendants of neutron star low-mass X-ray binaries. However, some studies indicate that the formation of pulsars from the accretion-induced collapse (AIC) of accreting white dwarfs (WDs) cannot be excluded. In this work, we use a population synthesis code to examine if the AIC channel can produce eccentric binary millisecond pulsars (BMSPs) in the Galaxy. Our simulated results indicate that, only when the natal MSPs receive a relatively strong kick ($\ga100\rm km\,s^{-1}$), can the AIC channel produce $\sim 10-180$ eccentric ($e>0.1$) BMSPs in the Galaxy, most of which are accompanied by a Helium star. Such a kick seems to be highly unlikely in the conventional AIC process, hence the probability of forming eccentric BMSPs via the AIC channel can be ruled out. Even if a high kick is allowed, the AIC channel cannot produce eccentric BMSPs with an orbital period of $\ga 20$ days. Therefore, we propose that the peculiar BMSP PSR J1903+0327 cannot be formed by the AIC channel. However, the AIC evolutionary channel may produce some fraction of isolated millisecond pulsars, and even sub-millisecond pulsars if they really exist.

preprint2011arXiv

He star evolutionary channel to intermediate-mass binary pulsar PSR J1802-2124

The intermediate-mass binary pulsars (IMBPs) are characterized by relatively long spin periods (10 - 200 ms) and massive ($\ga 0.4 M_{\odot}$) white dwarf (WD) companions. Recently, precise mass measurements have been performed for the pulsar and the WD in the IMBP PSR J1802-2124. Some observed properties, such as the low mass of the pulsar, the high mass of the WD, the moderately long spin period, and the tight orbit, imply that this system has undergone a peculiar formation mechanism. In this work, we attempt to simulate the detailed evolutionary history of PSR J1802-2124. We propose that a binary system consisting of a neutron star (NS, of mass $1.3 M_{\odot}$) and an He star (of mass $1.0 M_{\odot}$), and with an initial orbital period of 0.5 d, may have been the progenitor of PSR J1802-2124. Once the He star overflows its Roche lobe, He-rich material is transferred onto the NS at a relatively high rate of $\sim 10^{-7}-10^{-6} M_{\odot}\,\rm yr^{-1}$, which is significantly higher than the Eddington accretion rate. A large amount of the transferred material is ejected from the vicinity of the NS by radiation pressure and results in the birth of a mildly recycled pulsar. Our simulated results are consistent with the observed parameters of PSR J1802-2124. Therefore, we argue that the NS + He star evolutionary channel may be responsible for the formation of most IMBPs with orbital periods $\la 3 \rm d$.

preprint2011arXiv

KHz QPOs in LMXBs, relations between different frequencies and compactness of stars

We suggest that the mass of four compact stars SAX J1808.4$-$3658, KS 1731$-$260, SAX J1750.8$-$2900 and IGR J17191$-$2821 can be determined from the difference in the observed kiloHertz quasi periodic oscillations (kHz QPO-s) of these stars. The stellar radius is very close to the marginally stable orbit $R_{ms}$ as predicted by Einstein's general relativity. It may be noted that the first of these stars was suggested to be a strange star more than a decade back by Li \emph{et al.} (1999) from the unique millisecond X-ray pulsations with an accurate determination of its rotation period. It showed kHz QPO-s eight years back and so far it is the only set that has been observed. This is the first time we give an estimate of the mass of the star and of three other compact stars in Low-Mass X-ray Binaries using their observed kHz QPO-s.

preprint2011arXiv

On Pulsar-Driven Mass Ejection in Low-Mass X-ray Binaries

There is accumulating evidence for mass ejection in low-mass X-ray binaries (LMXBs) driven by radio pulsar activity during X-ray quiescence. In this paper we consider the condition for mass ejection by comparing the radiation pressure from a millisecond pulsar, and the gas pressure at the inner Lagrange point or at the surrounding accretion disk. We calculate the critical spin period of the pulsar below which mass ejection is allowed. Combining with the evolution of the mass transfer rate, we present constraints on the orbital periods of the systems. We show that mass ejection could happen in both wide and compact LMXBs. It may be caused by transient accretion due to thermal instability in the accretion disks in the former, and irradiation-driven mass-transfer cycles in the latter.

preprint2011arXiv

Strong field effects on emission line profiles: Kerr black holes and warped accretion disks

If an accretion disk around a black hole is illuminated by hard X-rays from non-thermal coronae, fluorescent iron lines will be emitted from the inner region of the accretion disk. The emission line profiles will show a variety of strong field effects, which may be used as a probe of the spin parameter of the black hole and the structure of the accretion disk. In this paper we generalize the previous relativistic line profile models by including both the black hole spinning effects and the non-axisymmetries of warped accretion disks. Our results show different features from the conventional calculations for either a flat disk around a Kerr black hole or a warped disk around a Schwarzschild black hole by presenting, at the same time, multiple peaks, rather long red tails and time variations of line profiles with the precession of the disk. We show disk images as seen by a distant observer, which are distorted by the strong gravity. Although we are primarily concerned with the iron K-shell lines in this paper, the calculation is general and is valid for any emission lines produced from a warped accretion disk around a black hole.

preprint2011arXiv

The contribution of X-ray binaries to the evolution of late-type galaxies: Evolutionary population synthesis simulations

X-ray studies of normal late-type galaxies have shown that non-nuclear X-ray emission is typically dominated by X-ray binaries, and provides a useful measure of star formation activity. We have modeled the X-ray evolution of late-type galaxies over the $\sim$ 14 Gyr of cosmic history, with an evolutionary population synthesis code developed by Hurley et al. Our calculations reveal a decrease of the X-ray luminosity-to-mass ratio $L_{\rm X}/M$ with time, in agreement with observations (Fig.~7$a$). We show that this decrease is a natural consequence of stellar and binary evolution and mass accumulating process in galaxies. The X-ray-to-optical luminosity ratio $L_{\rm X}/L_{\rm B}$ is found to be fairly constant (around $\sim 10^{30}$ erg\,s$^{-1}$$L_{\rm B,\odot}^{-1}$, Fig.~7$b$), and insensitive to the star formation history in the galaxies. The nearly constant value of $L_{\rm X}/L_{\rm B}$ is in conflict with the observed increase in $L_{\rm X}/L_{\rm B}$ from $z=0$ to 1.4. The discrepancy may be caused by intense obscured star formation activity that leads to nonlinear relationship between X-ray and B-band emission.

preprint2010arXiv

A General Relativistic Magnetohydrodynamic Model of High Frequency Quasi-periodic Oscillations in Black Hole Low-Mass X-ray Binaries

We suggest a possible explanation for the high frequency quasi-periodic oscillations (QPOs) in black hole low mass X-ray binaries. By solving the perturbation general relativistic magnetohydrodynamic equations, we find two stable modes of the Alfven wave in the the accretion disks with toroidal magnetic fields. We suggest that these two modes may lead to the double high frequency QPOs if they are produced in the transition region between the inner advection dominated accretion flow and the outer thin disk. This model naturally accounts for the 3 : 2 relation for the upper and lower frequencies of the QPOs, and the relation between the black hole mass and QPO frequency.

preprint2010arXiv

On the Binding Energy Parameter $λ$ of Common Envelope Evolution

The binding energy parameter $λ$ plays an important role in common envelope (CE) evolution. Previous works have already pointed out that $λ$ varies throughout the stellar evolution, though it has been adopted as a constant in most of the population synthesis calculations. We have systematically calculated the binding energy parameter $λ$ for both Population I and Population II stars of masses $1-20 M_{\sun}$, taking into account the contribution from the internal energy of stellar matter. We present fitting formulae for $λ$ that can be incorporated into future population synthesis investigations. We also briefly discuss the possible applications of the results in binary evolutions.

preprint2010arXiv

On the displacement of X-ray binaries from star clusters in star-burst galaxies

We have modeled the displacement of luminous X-ray binaries from star clusters in star-burst galaxies with an evolutionary population synthesis code developed by Hurley et al. (2000, 2002). In agreement with Kaaret et al. (2004), we find significant spatial offset of X-ray sources from their parent clusters, and the apparent X-ray luminosity vs. displacement correlation can be roughly reconstructed. The correlation is not sensitive to the fundamental properties of the clusters (e.g., initial mass functions of the binary stars) and the kick velocity imparted to the newborn compact stars, except the common envelope parameter. We present the distributions of the main parameters of the current X-ray binaries, which may be used to constrain the models for the formation and evolution of X-ray binaries with future optical observations.

preprint2009arXiv

Formation of ultra-compact X-ray binaries through circum-binary disk-driven mass transfer

The formation of ultra-compact X-ray binaries (UCXBs) has not been well understood. Previous works show that ultra-short orbital periods ($<1$ hr) may be reached through mass transfer driven by magnetic braking in normal low/intermediate-mass X-ray binaries (L/IMXBs) only for an extremely small range of initial binary parameters, which makes it difficult to account for the rather large population of UCXBs. In this paper we report the calculated results on mass transfer processes in L/IMXBs with a circum-binary disk. We show that when the orbital angular momentum loss due to a circum-binary disk is included, ultra-short orbital periods could be reached for a relatively wide range of initial binary parameters. The results of our binary models suggest an alternative formation channel for UCXBs.

preprint2009arXiv

The progenitors of Type Ia supernovae with long delay times

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.

preprint2008arXiv

A Fallback Disk Accretion-Involved Formation Channel to PSR J1903+0327

The discovery of the eccentric binary and millisecond pulsar PSR J1903+03273 has raised interesting questions about the formation mechanisms of this peculiar system. Here we present a born-fast scenario for PSR J1903+03273. We assume that during the supernova (SN) explosion that produced the pulsar, a fallback disk was formed around and accreted onto the newborn neutron star. Mass accretion could accelerate the neutron star's spin to milliseconds, and decrease its magnetic field to $\sim 10^8-10^9$ G, provided that there was sufficient mass ($\sim 0.1 M_{\sun}$) in the fallback disk. The neutron star became a millisecond pulsar after mass accretion terminated. In the meanwhile the binary orbit has kept to be eccentric (due to the SN explosion) for $\sim 10^{9}$ yr. We have performed population synthesis calculations of the evolutions of neutron stars with a fallback disk, and found that there might be tens to hundreds of PSR J1903+03273-like systems in the Galaxy. This scenario also suggests that some fraction of isolated millisecond pulsars in the Galactic disk could be formed through the same channel.

preprint2008arXiv

Spatial distribution of luminous X-ray binaries in spiral galaxies

We have modelled the spatial distribution of luminous X-ray binaries (XRBs) in spiral galaxies that are like the Milky Way using an evolutionary population synthesis code. In agreement with previous theoretical expectations and observations, we find that both high- and low-mass XRBs show clear concentrations towards the galactic plane and bulge.We also compare XRB distributions under the galactic potential with a dark matter halo and the modified Newtonian dynamics potential, and we suggest that the difference may serve as potential evidence to discriminate between these two types of model.

preprint2008arXiv

The bifurcation periods in low-mass X-ray binaries: the effect of magnetic braking and mass loss

The bifurcation period in low-mass X-ray binaries is the initial orbital pe- riod which separates the formation of converging systems (which evolve with decreasing orbital periods until the donor becomes degenerate) from the diverging systems (which evolve with increasing orbital periods until the donor star loses its envelope and a wide detached binary is formed). We calculate systematically the bifurcation periods of binary systems with a 1.4M_\sun neutron star and a 0.5-2M_\sun donor star, taking into account different kinds of magnetic braking and mass loss mechanisms. Our results show that the saturated magnetic braking can considerably decrease the values of bifurcation period compared to the traditional magnetic braking, while the influence of mass loss mechanisms on bifurcation periods is quite weak. We also develop a semi-analytical method to compute the bifurcation period, the result of which agrees well with the numerical method in the leading order.

preprint2006arXiv

Evolution of black-hole intermediate-mass X-ray binaries: the influence of a circumbinary disc

Justham, Rappaport & Podsiadlowski (2006) recently suggested that black-hole low-mass X-ray binaries (BHLMXBs) with short orbital periods may have evolved from black-hole intermediate-mass X-ray binaries (BHIMXBs). In their model the secondaries in BHIMXBs are assumed to possess anomalously high magnetic fields, so that magnetic braking can lead to substantial loss of angular momentum. In this paper we propose an alternative mechanism for orbital angular momentum loss in BHIMXBs. We assume that a small fraction $δ$ of the transferred mass from the donor star form a circumbinary disc surrounding the binary system. The tidal torques exerted by the disc can effectively drain orbital angular momentum from the binary. We have numerically calculated the evolutionary sequences of BHIMXBs, to examine the influence of the circumbinary disc on the binary evolution. Our results indicate when $δ\la 0.01-0.1$ (depending on the initial orbital periods), the circumbinary disc can cause secular orbital shrinking, leading to the formation of compact BHLMXBs, otherwise the orbits always expand during the evolution. This scenario also suggests the possible existence of luminous, persistent BHLMXBs, but it suffers the same problem as in Justham, Rappaport & Podsiadlowski (2006) that, the predicted effective temperatures of the donor stars are significantly higher than those of the observed donor stars in BHLMXBs.

preprint2006arXiv

On the Nature of Part Time Radio Pulsars

The recent discovery of rotating radio transients and the quasi-periodicity of pulsar activity in the radio pulsar PSR B1931$+$24 has challenged the conventional theory of radio pulsar emission. Here we suggest that these phenomena could be due to the interaction between the neutron star magnetosphere and the surrounding debris disk. The pattern of pulsar emission depends on whether the disk can penetrate the light cylinder and efficiently quench the processes of particle production and acceleration inside the magnetospheric gap. A precessing disk may naturally account for the switch-on/off behavior in PSR B1931$+$24.