We present Giant Metrewave Radio Telescope (GMRT) 323 MHz radio continuum observations toward 13 radio-loud quasars at $z>5$, sampling the low-frequency synchrotron emission from these objects. Among the 12 targets successfully observed, we detected 10 above $4\sigma$ significance, while 2 remain undetected. All of the detected sources appear as point sources. Combined with previous radio continuum detections from the literature, 9 quasars have power-law spectral energy distributions throughout the radio range; for some the flux density drops with increasing frequency while it increases for others. Two of these sources appear to have spectral turnover. For the power-law-like sources, the power-law indices have a positive range between 0.18 and 0.67 and a negative values between $-0.90$ and $-0.27$. For the turnover sources, the radio peaks around $\sim1$ and $\sim10$ GHz in the rest frame, the optically thin indices are $-0.58$ and $-0.90$, and the optically thick indices are 0.50 and 1.20. A magnetic field and spectral age analysis of SDSS J114657.59+403708.6 at $z=5.01$ may indicate that the turnover is not caused by synchrotron self-absorption, but rather by free-free absorption by the high-density medium in the nuclear region. Alternatively, the apparent turnover may be an artifact of source variability. Finally, we calculated the radio loudness $R_{2500\rm\, \AA}$ for our sample, which spans a very wide range from 12$^{+13}_{-13}$ to 4982$^{+279}_{-254}$.

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.

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}$).

The Beijing-Arizona Sky Survey (BASS) is a wide and deep imaging survey to cover a 5400 deg$^2$ area in the Northern Galactic Cap with the 2.3m Bok telescope using two filters ($g$ and $r$ bands). The Mosaic $z$-band Legacy Survey (MzLS) covers the same area in $z$ band with the 4m Mayall telescope. These two surveys will be used for spectroscopic targeting of the Dark Energy Spectroscopic Instrument (DESI). The BASS survey observations were completed in 2019 March. This paper describes the third data release (DR3) of BASS, which contains the photometric data from all BASS and MzLS observations between 2015 January and 2019 March. The median astrometric precision relative to {\it Gaia} positions is about 17 mas and the median photometric offset relative to the PanSTARRS1 photometry is within 5 mmag. The median $5\sigma$ AB magnitude depths for point sources are 24.2, 23.6, and 23.0 mag for $g$, $r$, and $z$ bands, respectively. The photometric depth within the survey area is highly homogeneous, with the difference between the 20\% and 80\% depth less than 0.3 mag. The DR3 data, including raw data, calibrated single-epoch images, single-epoch photometric catalogs, stacked images, and co-added photometric catalogs, are publicly accessible at \url{http://batc.bao.ac.cn/BASS/doku.php?id=datarelease:home}.

We report Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO $(8-7)$, $(9-8)$, $\rm H_{2}O (2_{0,2}-1_{1,1})$ and $\rm OH^{+} (1_{1}-0_{1})$ and NOrthern Extended Millimeter Array (NOEMA) observations of CO $(5-4)$, $(6-5)$, $(12-11)$ and $(13-12)$ towards the $z = 6.003$ quasar SDSS J231038.88+185519.7, aiming to probe the physical conditions of the molecular gas content of this source. We present the best sampled CO spectral line energy distribution (SLED) at $z = 6.003$, and analyzed it with the radiative transfer code MOLPOP-CEP. Fitting the CO SLED to a one-component model indicates a kinetic temperature $T_{\rm kin} = 228 \ \rm K$, molecular gas density $log (n(\rm H_{2})/\rm cm^{-3}$ )=4.75, and CO column density $log(N(\rm CO)/\rm cm^{-2}) =17.5$, although a two-component model better fits the data. In either case, the CO SLED is dominated by a "warm" and "dense" component. Compared to samples of local (Ultra) Luminous Infrared Galaxies ((U)LIRGs), starburst galaxies and high redshift Submillimeter Galaxies (SMGs), J2310+1855 exhibits higher CO excitation at ($J \geq 8$), like other high redshift quasars. The high CO excitation, together with the enhanced $L_{\rm H_{2}O}/ L_{IR} $, $L_{\rm H_{2}O}/ L_{CO} $ and $L_{OH^{+}}/L_{\rm H_{2}O} $ ratios, suggests that besides the UV radiation from young massive stars, other mechanisms such as shocks, cosmic rays and X-rays might also be responsible for the heating and ionization of the molecular gas. In the nuclear region probed by the molecular emissions lines, any of these mechanisms might be present due to the powerful quasar and the starburst activity.

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

We present ALMA observations of the [CII] 158 $\mu$m fine structure line and dust continuum emission from two quasars, SDSS J104433.04-012502.2 and SDSS J012958.51-003539.7, at z=5.78. The ALMA observations at 0.2'' resolution map the dust and gas on kpc scales. The spatially resolved emission shows a similar trend of decreasing [CII]-FIR ratios with increasing FIR surface brightnesses as was found in the infrared luminous galaxies with intense star formation. We confirm the velocity gradients of [CII] emission found previously in SDSS J0129-0035. No clear evidence of order motion is detected in SDSS J1044-0125. The velocity maps and PV-diagrams also suggest turbulent gas clumps in both objects. We tentatively detect a [CII] peak offset 4.9 kpc to the East of SDSS J1044-0125. This may be associated with an infalling companion, or node of gas outflow. All these results suggest significant dynamical evolution of the ISM in the nuclear region of these young quasar-starburst systems. We fit the velocity map of the [CII] emission from SDSS J0129-0035 with a rotating disk model. The result suggests a face-on system with an inclination angle of 16 +/- 20 degree and constrains the lower limit of the host galaxy dynamical mass to be 2.6x10^10 Msun within the [CII] emitting region. It is likely that SDSS J0129-0035, as well as other young quasars with super massive black hole masses on order of 10^7 to 10^8 Msun, falls close to the black hole and host galaxy mass relation defined by local galaxies.

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} $.

Optical observations of normal-stars in binary systems with massive unseen objects have been proposed to search for candidate black holes (BHs) and provide a direct measurement of their dynamical masses. In this paper, we have performed binary population synthesis calculations to simulate the potential population of detached binaries containing BHs and normal-star companions in the Galaxy. We focus on the influence of the BH progenitors. In the traditional model, BHs in binaries evolve from stars more massive than $ \sim25M_{\odot} $. However, it is difficult for this model to produce BH low-mass X-ray binaries. Recent investigations on massive star evolution suggest that the BH progenitors may have masses as low as $ \sim15M_{\odot} $. Based on this model, we provide the expected distributions of various parameters for detached BH binaries with normal-star companions, including the component masses, the orbital parameters of the binary systems, the radial velocity semi-amplitudes, and the astrometric signatures of the optical companions. %the detached BH systems with low-mass companions can be easily produced in our simulations. Our calculations show that there are more than thousands of such detached binaries in the Galaxy, and hundreds of them are potentially observable systems with luminous companions brighter than 20 mag. In addition, detached BH binaries are dominated by those with main-sequence companions and only a few percent of them are expected to have giant companions.

V1082 Sgr is a cataclysmic variable with accretion luminosity above $10^{34}$erg s$^{-1}$, indicating a mass transfer rate above $10^{-9}M_{\odot}$yr$^{-1}$. However, its K type companion was suggested to be under-filling its Roche lobe (RL), making the high mass transfer rate a mystery. In this work we propose a possible model to explain this discrepancy. The system is proposed to be an intermediate polar, with its K type companion filling its RL. The mass of the white dwarf star is evaluated to be $0.77\pm0.11M_{\odot}$ from both X-ray continuum fitting and Fe line flux ratio measurements. We make numerical simulations to search for the possible progenitors of the system. The results show that a binary with an initial 1.5--2.5$M_{\odot}$ companion in a 1--2 day orbit (or an initial 1.0--1.4$M_{\odot}$ companion in a 3.2--4.1 day orbit) may naturally evolve to a cataclysmic variable with a $\sim 0.55 \pm 0.11M_{\odot}$, Roche-lobe filling companion in a 0.86 day orbit. The effective temperature of the donor star, the mass transfer rate, and the derived V band magnitude are all consistent with previous observations.

We present Karl G. Jansky Very Large Array (VLA) observations of the CO (2$-$1) line emission towards three far-infrared luminous quasars at $z\sim6$: SDSS J231038.88$+$185519.7 and SDSS J012958.51$-$003539.7 with $\sim0\farcs6$ resolution and SDSS J205406.42$-$000514.8 with $\sim2\farcs1$ resolution. All three sources are detected in the CO (2$-$1) line emission -- one source is marginally resolved, and the other two appear as point sources. Measurements of the CO (2$-$1) line emission allow us to calculate the molecular gas mass even without a CO excitation model. The inferred molecular gas masses are (0.8$-$4.3) $\times$ 10$^{10}$ $M_{\odot}$. The widths and redshifts derived from the CO (2$-$1) line are consistent with previous CO (6$-$5) and [\ion{C}{2}] measurements. We also report continuum measurements using the Herschel for SDSS J231038.88$+$185519.7 and SDSS J012958.51$-$003539.7, and for SDSS J231038.88+185519.7, data obtained at $\sim140$ and $\sim300$ GHz using the Atacama Large Millimeter/submillimeter Array (ALMA). In the case of SDSS J231038.88+185519.7, we present a detailed analysis of the spectral energy distribution and derive the dust temperature ($\sim40$ K), the dust mass ($\sim10^{9}$ $M_{\odot}$), the far-infrared luminosity (8$-$1000 $\mu$m; $\sim10^{13}$ $ L_{\odot}$) and the star formation rate (2400$-$2700 $M_{\odot}$ yr$^{-1}$). Finally, an analysis of the photo-dissociation regions associated with the three high redshift quasars indicates that the interstellar medium in these sources has similar properties to local starburst galaxies.

We perform simulations of the capabilities of the next generation Very Large Array in the context of imaging low order CO emission from typical high redshift star forming galaxies at ~ 1 kpc resolution. We adopt as a spatial and dynamical template the CO 1-0 emission from M 51, scaled accordingly for redshift, transition, and total gas mass. The molecular gas masses investigated are factors of 1.4, 3.5, and 12.5 larger that of M 51, at z = 0.5, 2, and 4.2, respectively. The z = 2 galaxy gas mass is comparable to the lowest mass galaxies currently being discovered in the deepest ALMA and NOEMA cosmological CO line surveys, corresponding to galaxies with star formation rates ~ 10 to 100 $M_\odot$ yr$^{-1}$. The ngVLA will perform quality imaging at 1kpc resolution of the gas distribution and dynamics over this disk. We recover the overall rotation curve, galaxy orientation properties, and molecular ISM internal velocity dispersion. The model at z = 4.2 corresponds to a massive star forming main sequence disk (SFR ~ 130 $M_\odot$ yr$^{-1}$). The ngVLA can obtain 1kpc resolution images of such a system in a reasonable integration time, and recover the basic galaxy orientation parameters, and, asymptotically, the maximum rotation velocity. We compare the ngVLA results with capabilities of ALMA and the Jansky VLA. ALMA and the VLA can detect the integrated low order CO emission from these galaxies, but lack the sensitivity to perform the high resolution imaging to recover the dynamics at 1kpc scales. To do so would require of order 1000 hrs per galaxy with these current facilities. We investigate a 'minimal' ngVLA configuration, removing the longest baselines and much of the very compact core, and find good imaging can still be performed at 1 kpc resolution.

In this work we focus on a group of Galactic double neutron star (DNS) systems with long orbital periods of $ \gtrsim 1$ day and low eccentricities of $\lesssim 0.4$. The feature of these orbital parameters is used to constrain the evolutionary processes of progenitor binaries and the supernova (SN) kicks of the second born NSs. Adopting that the mass transfer during primordial binary evolution is highly non-conservative (rotation-dependent), the formation of DNS systems involves a double helium star binary phase, the common envelope (CE) evolution initiates before the first NS formation. During the CE evolution the binary orbital energy is obviously larger when using a helium star rather than a NS to expel the donor envelope, this can help explain the formation of DNS systems with long periods. SN kicks at NS birth can lead to eccentric orbits and even the disruption of binary systems, the low eccentricities require that the DNSs receive a small natal kick at the second collapse. Compared with the overall distribution of orbital parameters for observed DNS binaries, we propose that the second born NSs in most DNS systems are subject to small natal kicks with the Maxwellian dispersion velocity of less than $ 80 \,\rm km\,s^{-1} $, which can provide some constraints on the SN explosion processes. The mass distribution of DNS binaries is also briefly discussed. We suggest that the rotation-dependent mass transfer mode and our results about SN kicks should be applied to massive binary evolution and population synthesis studies.

JINGLE is a new JCMT legacy survey designed to systematically study the cold interstellar medium of galaxies in the local Universe. As part of the survey we perform 850um continuum measurements with SCUBA-2 for a representative sample of 193 Herschel-selected galaxies with M*>10^9Msun, as well as integrated CO(2-1) line fluxes with RxA3m for a subset of 90 of these galaxies. The sample is selected from fields covered by the Herschel-ATLAS survey that are also targeted by the MaNGA optical integral-field spectroscopic survey. The new JCMT observations combined with the multi-wavelength ancillary data will allow for the robust characterization of the properties of dust in the nearby Universe, and the benchmarking of scaling relations between dust, gas, and global galaxy properties. In this paper we give an overview of the survey objectives and details about the sample selection and JCMT observations, present a consistent 30 band UV-to-FIR photometric catalog with derived properties, and introduce the JINGLE Main Data Release (MDR). Science highlights include the non-linearity of the relation between 850um luminosity and CO line luminosity, and the serendipitous discovery of candidate z>6 galaxies.

The Beijing-Arizona Sky Survey (BASS) is a new wide-field legacy imaging survey in the northern Galactic cap using the 2.3m Bok telescope. The survey will cover about 5400 deg$^2$ in the $g$ and $r$ bands, and the expected 5$\sigma$ depths (corrected for the Galactic extinction) in the two bands are 24.0 and 23.4 mag, respectively. BASS started observations in January 2015, and has completed about 41% of the whole area as of July 2016. The first data release contains both calibrated images and photometric catalogs obtained in 2015 and 2016. The depths of single-epoch images in the two bands are 23.4 and 22.9 mag, and the full depths of three epochs are about 24.1 and 23.5 mag, respectively.

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

We have performed population synthesis calculation on the formation of binaries containing a black hole (BH) and a neutron star (NS) in the Galactic disk. Some of important input parameters, especially for the treatment of common envelope evolution, are updated in the calculation. We have discussed the uncertainties from the star formation rate of the Galaxy and the velocity distribution of NS kicks on the birthrate ($ \sim 0.6-13 \rm\, Myr^{-1}$) of BH/NS binaries. From incident BH/NS binaries, by modelling the orbital evolution duo to gravitational wave radiation and the NS evolution as radio pulsars, we obtain the distributions of the observable parameters such as the orbital period, eccentricity and pulse period of the BH/pulsar binaries. We estimate that there may be $\sim 3-80 $ BH/pulsar binaries in the Galactic disk and around 10\% of them could be detected by the Five-hundred-meter Aperture Spherical radio Telescope.

Until recently, determining the rotational properties of galaxies in the early universe (z>4, Universe age <1.5Gyr) was impractical, with the exception of a few strongly lensed systems. Combining the high resolution and sensitivity of ALMA at (sub-) millimeter wavelengths with the typically high strength of the [CII] 158um emission line from galaxies and long-developed dynamical modeling tools raises the possibility of characterizing the gas dynamics in both extreme starburst galaxies and normal star forming disk galaxies at z~4-7. Using a procedure centered around GIPSY's ROTCUR task, we have fit tilted ring models to some of the best available ALMA [CII] data of a small set of galaxies: the MS galaxies HZ9 & HZ10, the Damped Lyman-alpha Absorber (DLA) host galaxy ALMA J0817+1351, the submm galaxies AzTEC/C159 and COSMOS J1000+0234, and the quasar host galaxy ULAS J1319+0950. This procedure directly derives rotation curves and dynamical masses as functions of radius for each object. In one case, we present evidence for a dark matter halo of O(10^11) solar masses. We present an analysis of the possible velocity dispersions of AzTEC/C159 and ULAS J1319+0950 based on matching simulated observations to the integrated [CII] line profiles. Finally, we test the effects of observation resolution and sensitivity on our results. While the conclusions remain limited at the resolution and signal-to-noise ratios of these observations, the results demonstrate the viability of the modeling tools at high redshift, and the exciting potential for detailed dynamical analysis of the earliest galaxies, as ALMA achieves full observational capabilities.

We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations of the dust continuum and [C II] 158 $\mu$m fine structure line emission towards a far-infrared-luminous quasar, ULAS J131911.29$+$095051.4 at $z=6.13$, and combine the new Cycle 1 data with ALMA Cycle 0 data. The combined data have an angular resolution $\sim$ $0.3$, and resolve both the dust continuum and the [C II] line emission on few kpc scales. The [C II] line emission is more irregular than the dust continuum emission which suggests different distributions between the dust and [C II]-emitting gas. The combined data confirm the [C II] velocity gradient that we previously detected in lower resolution ALMA image from Cycle 0 data alone. We apply a tilted ring model to the [C II] velocity map to obtain a rotation curve, and constrain the circular velocity to be 427 $\pm$ 55 km s$^{-1}$ at a radius of 3.2 kpc with an inclination angle of 34$^\circ$. We measure the dynamical mass within the 3.2 kpc region to be 13.4$_{-5.3}^{+7.8}$ $\times 10^{10}\,M_{\odot}$. This yields a black hole and host galaxy mass ratio of 0.020$_{-0.007}^{+0.013}$, which is about 4$_{-2}^{+3}$ times higher than the present-day $M_{\rm BH}$/$M_{\rm bulge}$ ratio. This suggests that the supermassive black hole grows the bulk of its mass before the formation of the most of stellar mass in this quasar host galaxy in the early universe.

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.

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$.

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.

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.

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.

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 $\lambda$ in common envelope evolution, and then incorporate both $q_{\rm cr}$ and $\lambda$ 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.

There is a remarkable correlation between the spin periods of the accreting neutron stars in Be/X-ray binaries (BeXBs) and their orbital periods . Recently Knigge et al. (2011) showed that the distribution of the spin periods contains two distinct subpopulations peaked at $\sim 10$ s and $\sim 200$ s respectively, and suggested that they may be related to two types of supernovae for the formation of the neutron stars, i.e., core-collapse and electron-capture supernovae. Here we propose that the bimodal spin period distribution is likely to be ascribed to different accretion modes of the neutron stars in BeXBs. When the neutron star tends to capture material from the warped, outer part of the Be star disk and experiences giant outbursts, a radiatively-cooling dominated disk is formed around the neutron star, which spins up the neutron star, and is responsible for the short period subpopulation. In BeXBs that are dominated by normal outbursts or persistent, the accretion flow is advection-dominated or quasi-spherical. The spin-up process is accordingly inefficient, leading to longer periods of the neuron stars. The potential relation between the subpopulations and the supernova mechanisms is also discussed.

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.

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}$.

Considering that the universe is filled with the nonrelativistic matter and dark energy and each component is respectively satisfied with its conservation condition in the absence of their interaction, we give the change rate of the fractional density and the density of dark energy from the conservation condition. It is clear that the fractional density of dark energy will monotonously increase and gradually become the dominating contribution to the universe as the redshift becomes low. Combining the evolutional trend of the state equation of dark energy and the change rate of the density of dark energy we find that the density of dark energy will decrease up to a minimum and whereafter it will increase again as the redshift becomes low. This can be regarded as the phase transition of dark energy from the quintessence phase to the phantom phase.

We propose that complex symmetry gravitational theory (CSGT) explain the accelerating expansion of universe. In this paper universe is taken as the double complex symmetric space. Cosmological solution is obtained within CSGT. The conditions of the accelerating expansion of universe are discussed within CSGT. Moreover, the range of equation of state of matter $\omega_\epsilon$ is given in the hyperbolic imaginary space.

We assume that dark matter and dark energy satisfy the unified equation of state: $p=B(z)\rho$, with $p=p_{dE}$, $\rho=\rho_{dm}+\rho_{dE}$, where the pressure of dark matter $p_{dm}=0$ has been taken into account. A special function $B=-\frac{A}{(1+z)^{\alpha}}$ is presented, which can well describe the evolution of the universe. In this model, the universe will end up with a Big Rip. By further simple analysis, we know other choices of the function $B$ can also describe the universe but lead to a different doomsday.