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

Most dynamically confirmed stellar-mass black holes and the candidates were originally selected from X-ray outbursts. In the present work, we search for black hole candidates in the LAMOST survey by using the spectra along with photometry from the ASAS-SN survey, where the orbital period of the binary may be revealed by the periodic light curve, such as the ellipsoidal modulation type. Our sample consists of 9 binaries, where each source contains a giant star with large radial velocity variation ($\Delta V_{\rm R} > 70~{\rm km~s^{-1}}$) and periods known from light curves. We focus on the 9 sources with long periods ($T_{\rm ph} > 5$ days) and evaluate the mass $M_2$ of the optically invisible companion. Since the observed $\Delta V_{\rm R}$ from only a few repeating spectroscopic observations is a lower limit of the real amplitude, the real mass $M_2$ can be significantly higher than the current evaluation. It is likely an efficient method to place constraints on $M_2$ by combining $\Delta V_{\rm R}$ from LAMOST and $T_{\rm ph}$ from ASAS-SN, particularly by the ongoing LAMOST Medium Resolution Survey.

We propose a method to search for stellar-mass black hole (BH) candidates with giant companions from spectroscopic observations. Based on the stellar spectra of LAMOST Data Release 6, we obtain a sample of seven giants in binaries with large radial velocity variation $\Delta V_R > 80~{\rm km~s^{-1}}$. With the effective temperature, surface gravity, and metallicity provided by LAMOST, and the parallax given by {\it Gaia}, we can estimate the mass and radius of the giant, and therefore evaluate the possible mass of the optically invisible star in the binary. We show that the sources in our sample are potential BH candidates, and are worthy of dynamical measurement by further spectroscopic observations. Our method may be particularly valid for the selection of BH candidates in binaries with unknown orbital periods.

Stars of spectral type Oe are very rare. To date, only 13 Oe stars have been identified within our Galaxy. In this paper, we present six new Oe stars and four new B0e stars found in LAMOST DR5. Repeated spectral observations of the same Oe stars show some emission line variability. The H$\beta$ emission of TYC 4801-17-1 shows rapid V/R variation. Phase lags in the V/R ratio of TYC 4801-17-1 spectra are also seen. We found the unusual O4.5 star RL 128 is an Oe star with variable H$\alpha$ intensity and its Ca II triplet emission appears when H$\alpha$ emission reaches maximum intensity. These newly identified early type Oe and B0e stars significantly increase the known sample.

With Rossi X-ray Timing Explorer (RXTE) data, we systematically study the cross-correlation and time lag in all spectral states of black hole X-ray binary (BHXB) XTE J1859+226 in detail during its entire 1999-2000 outburst that lasted for 166 days. Anti-correlations and positive correlations and their respective soft and hard X-ray lags are only detected in the first 100 days of the outburst when the luminosity is high. This suggests that the cross-correlations may be related to high luminosity. Positive correlations are detected in every state of XTE J1859+226, viz., hard state, hard-intermediate state (HIMS), soft-intermediate state (SIMS) and soft state. However, anti-correlations are only detected in HIMS and SIMS, anti-correlated hard lags are only detected in SIMS, while anti-correlated soft lags are detected in both HIMS and SIMS. Moreover, the ratio of the observations with anti-correlated soft lags to hard lags detected in XTE J1859+226 is significantly different from that in neutron star low-mass X-ray binaries (NS LMXBs). So far, anti-correlations are never detected in the soft state of BHXBs but detected in every branch or state of NS LMXBs. This may be due to the origin of soft seed photons in BHXBs is confined to the accretion disk and, for NS LMXBs, from both accretion disk and the surface of the NS. We notice that the timescale of anti-correlated time lags detected in XTE J1859+226 is similar with that of other BHXBs and NS LMXBs. We suggest that anti-correlated soft lag detected in BHXB may result from fluctuation in the accretion disk as well as NS LMXB.

Sky subtraction is the key technique in data reduction of multi-fiber spectra. Knowledge of the related instrument character is necessary to determine the method adopted in sky subtraction. In this study, we described the sky subtraction method designed for LAMOST(Large sky Area Multi-Object fiber Spectroscopic Telescope) survey. The method has been intergrated into LAMOST 2D Pipeline v2.6 and applied to data of LAMOST DR3 and later. For LAMOST, sky emission line calibration is used to alleviate the position-dependent (thus time-dependent) ~4% fiber throughput uncertainty and the small wavelength instability (0.1\AA ) during observation. PCA (Principal Component Analysis) sky subtraction further reduces 25% of the sky line residual of the OH lines in the red part of the LAMOST spectra after the mater sky spectrum, which is derived from a B-spline fit of 20 sky fibers in each spectrograph, is adjusted by sky emission line and subtracted from each fiber. Further analysis shows that our wavelength calibration accuracy is about 4.5km/s, and the average sky subtraction residuals are about 3% for sky emission lines and 3% for continuum region. The relative sky subtraction residuals vary with the moon light background brightness, could reach as low as 1.5% for the sky emission line regions in the dark night. Tests on the F stars of both similar sky emission line strength and similar object continuum intensity show that the sky emission line residual of LAMOST is smaller than those of SDSS survey.

We present candidate members of the Pal 5, GD-1, Cetus Polar, and Orphan tidal stellar streams found in LAMOST DR3, SDSS DR9 and APOGEE catalogs. In LAMOST DR3, we find 20, 4, 24 high confidence candidates of tidal streams GD-1, Cetus Polar and Orphan respectively. We also list from the SDSS DR9 spectroscopic catalog 59, 118, 10 high confidence candidates of tidal streams Cetus Polar, Orphan and Pal 5, respectively. Furthermore, we find 7 high confidence candidates of the Pal 5 tidal stream in APOGEE data. Compared with SDSS, the new candidates from LAMOST DR3 are brighter, so that together, more of the color-magnitude diagram, including the giant branch can be explored. Analysis of SDSS data shows that there are 3 metallicity peaks of the Orphan stream and also shows some spatial separation. LAMOST data confirms multiple metallicities in this stream. The metallicity, given by the higher resolution APOGEE instrument, of the Pal 5 tidal stream is [Fe/H] $\sim -1.2$, higher than that given earlier by SDSS spectra. Many previously unidentified stream members are tabulated here for the first time, along with existing members, allowing future researchers to further constrain the orbits of these objects as they move within the Galaxy's dark matter potential.

We present the second release of value-added catalogues of the LAMOST Spectroscopic Survey of the Galactic Anticentre (LSS-GAC DR2). The catalogues present values of radial velocity $V_{\rm r}$, atmospheric parameters --- effective temperature $T_{\rm eff}$, surface gravity log$g$, metallicity [Fe/H], $\alpha$-element to iron (metal) abundance ratio [$\alpha$/Fe] ([$\alpha$/M]), elemental abundances [C/H] and [N/H], and absolute magnitudes ${\rm M}_V$ and ${\rm M}_{K_{\rm s}}$ deduced from 1.8 million spectra of 1.4 million unique stars targeted by the LSS-GAC since September 2011 until June 2014. The catalogues also give values of interstellar reddening, distance and orbital parameters determined with a variety of techniques, as well as proper motions and multi-band photometry from the far-UV to the mid-IR collected from the literature and various surveys. Accuracies of radial velocities reach 5kms$^{-1}$ for late-type stars, and those of distance estimates range between 10 -- 30 per cent, depending on the spectral signal-to-noise ratios. Precisions of [Fe/H], [C/H] and [N/H] estimates reach 0.1dex, and those of [$\alpha$/Fe] and [$\alpha$/M] reach 0.05dex. The large number of stars, the contiguous sky coverage, the simple yet non-trivial target selection function and the robust estimates of stellar radial velocities and atmospheric parameters, distances and elemental abundances, make the catalogues a valuable data set to study the structure and evolution of the Galaxy, especially the solar-neighbourhood and the outer disk.

We analyze the relation between the emission radii of twin kilohertz quasi-periodic oscillations (kHz QPOs) and the co-rotation radii of the 12 neutron star low mass X-ray binaries (NS-LMXBs) which are simultaneously detected with the twin kHz QPOs and NS spins. We find that the average co-rotation radius of these sources is r_co about 32 km, and all the emission positions of twin kHz QPOs lie inside the corotation radii, indicating that the twin kHz QPOs are formed in the spin-up process. It is noticed that the upper frequency of twin kHz QPOs is higher than NS spin frequency by > 10%, which may account for a critical velocity difference between the Keplerian motion of accretion matter and NS spin that is corresponding to the production of twin kHz QPOs. In addition, we also find that about 83% of twin kHz QPOs cluster around the radius range of 15-20 km, which may be affected by the hard surface or the local strong magnetic field of NS. As a special case, SAX J1808.4-3658 shows the larger emission radii of twin kHz QPOs of r about 21-24 km, which may be due to its low accretion rate or small measured NS mass (< 1.4 solar mass).

We present preliminary results of the quasar survey in Large Sky Area Multi- Object Fiber Spectroscopic Telescope (LAMOST) first data release (DR1), which includes pilot survey and the first year regular survey. There are 3921 quasars identified with reliability, among which 1180 are new quasars discovered in the survey. These quasars are at low to median redshifts, with highest z of 4.83. We compile emission line measurements around the H{\alpha}, H{\beta}, Mg II, and C IV regions for the new quasars. The continuum luminosities are inferred from SDSS photo- metric data with model fitting as the spectra in DR1 are non-flux-calibrated. We also compile the virial black hole mass estimates, and flags indicating the selec- tion methods, broad absorption line quasars. The catalog and spectra for these quasars are available online. 28% of the 3921 quasars are selected with optical- infrared colours independently, indicating that the method is quite promising in completeness of quasar survey. LAMOST DR1 and the on-going quasar survey will provide valuable data in the studies of quasars.

We infer the emission positions of twin kilohertz quasi-periodic oscillations (kHz QPOs) in neutron star low mass X-ray binaries (NS-LMXBs) based on the Alfven wave oscillation model (AWOM). For most sources, the emission radii of kHz QPOs cluster around a region of 16-19 km with the assumed NS radii of 15 km. Cir X-1 has the larger emission radii of 23-38 km than those of the other sources, which may be ascribed to its large magnetosphere-disk radius or strong NS surface magnetic field. SAX J1808.4-3658 is also a particular source with the relative large emission radii of kHz QPOs of 20 - 23 km, which may be due to its large inferred NS radius of 18 - 19 km. The emission radii of kHz QPOs for all the sources are larger than the NS radii, and the possible explanations of which are presented. The similarity of the emission radii of kHz QPOs (16-19 km) for both the low/high luminosity Atoll/Z sources is found, which indicates that both sources share the similar magnetosphere- disk radii.

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

LAMOST has released more than two million spectra, which provide the opportunity to search for double-peaked narrow emission line (NEL) galaxies and AGNs. The double-peaked narrow-line profiles can be well modeled by two velocity components, respectively blueshifted and redshifted with respect to the systemic recession velocity. This paper presents 20 double-peaked NEL galaxies and AGNs found from LAMOST DR1 using a search method based on multi-gaussian fit of the narrow emission lines. Among them, 10 have already been published by other authors, either listed as genuine double-peaked NEL objects or as asymmetric NEL objects, the remaining 10 being first discoveries. We discuss some possible origins for double-peaked narrow-line features, as interaction between jet and narrow line regions, interaction with companion galaxies and black hole binaries. Spatially resolved optical imaging and/or follow-up observations in other spectral bands are needed to further discuss the physical mechanisms at work.

Spectral type recognition is an important and fundamental step of large sky survey projects in the data reduction for further scientific research, like parameter measurement and statistic work. It tends out to be a huge job to manually inspect the low quality spectra produced from massive spectroscopic survey, where the automatic pipeline may not provide confident type classification results. In order to improve the efficiency and effectiveness of spectral classification, we develop a semi-automated toolkit named ASERA, A Spectrum Eye Recognition Assistant. The main purpose of ASERA is to help the user in quasar spectral recognition and redshift measurement. Furthermore it can also be used to recognize various types of spectra of stars, galaxies and AGNs (Active Galactic Nucleus). It is an interactive software allowing the user to visualize observed spectra, superimpose template spectra from the Sloan Digital Sky Survey (SDSS), and interactively access related spectral line information. It is an efficient and user-friendly toolkit for accurate classification of spectra observed by LAMOST (the Large Sky Area Multi-object Fiber Spectroscopic Telescope). The toolkit is available in two modes: a Java standalone application and a Java applet. ASERA has a few functions, such as wavelength and flux scale setting, zoom in and out, redshift estimation, spectral line identification, which helps user to improve the spectral classification accuracy especially for low quality spectra and reduce the labor of eyeball check. The function and performance of this tool is displayed through the recognition of several quasar spectra and a late type stellar spectrum from the LAMOST Pilot survey. Its future expansion capabilities are discussed.

Using all the observations from Rossi X-ray Timing Explorer for accreting neutron star XTE J1701-462, we carry out a systematic study on the cross-correlation between its soft and hard light curves. The anti-correlations appear at the hard vertex and on the upper normal branch (NB) in the Cyg-like interval, occur on the horizontal branch (HB) and/or upper NB in the first and second Sco-like intervals, and display at the upper flaring branch (FB) in the third Sco-like interval. We suggest that the anti-correlation might evolve with luminosity, because with decreasing of luminosity, the source evolved from a Cyg-like Z source, via a Sco-like Z source, to an atoll source. In the Cyg-like interval, the positive correlations are mostly distributed on the HB, which is not consistent with that of the Cyg-like Z source GX 5-1 and Cyg X-2 whose HBs host ambiguous correlations and anti-correlations. From our spectral analyses, it is found that the ratio of the hard emission to the soft emission basically keeps unvaried from the HB to the hard vertex, which might result in the positive correlation. Additionally, we compare the spectra in the third Sco-like interval with those of atoll source 4U 1735-44 and find that their spectral evolution along the tracks in the color-color diagrams are similar, indicating that in this interval, XTE J1701-462 approaches an atoll source. Both hard time lags and soft time lags are detected and the possible origins of them are discussed.

We collect the data of twin kilohertz quasi-periodic oscillations (kHz QPOs) published before 2012 from 26 neutron star (NS) low-mass X-ray binary (LMXB) sources, then we analyze the centroid frequency (\nu) distribution of twin kHz QPOs (lower frequency \nu_1 and upper frequency \nu_2) both for Atoll and Z sources. For the data without shift-and-add, we find that Atoll and Z sources show the different distributions of \nu_1, \nu_2 and \nu_2/\nu_1, but the same distribution of \Delta\nu (difference of twin kHz QPOs), which indicates that twin kHz QPOs may share the common properties of LXMBs and have the same physical origins. The distribution of \Delta\nu is quite different from constant value, so is \nu_2/\nu_1 from constant ratio. The weighted mean values and maxima of \nu_1 and \nu_2 in Atoll sources are slightly higher than those in Z sources. We also find that shift-and-add technique can reconstruct the distribution of \nu_1 and \Delta\nu. The K-S test results of \nu_1 and \Delta\nu between Atoll and Z sources from data with shift-and-add are quite different from those without it, and we think that this may be caused by the selection biases of the sample. We also study the properties of the quality factor (Q) and the root-mean-squared (rms) amplitude of 4U 0614+09 with the data from the two observational methods, but the errors are too big to make a robust conclusion. The NS spin frequency (\nu_s) distribution of 28 NS-LMXBs show a bigger mean value (about 408Hz) than that (about 281 Hz) of the radio binary millisecond pulsars (MSPs), which may be due to the lack of the spin detections from Z sources (systematically lower than 281 Hz). Furthermore, on the relations between the kHz QPOs and NS spin frequency \nu_s, we find the approximate correlations of the mean values of \Delta\nu with NS spin and its half, respectively.

We test the relativistic precession model (RPM) and the MHD Alfven wave oscillation model (AWOM) for the kHz QPOs by the sources with measured NS masses and twin kHz QPO frequencies. For RPM, the derived NS mass of Cyg X-2 (SAX J1808.4-3658 and 4U 1820-30) is 1.96 +/- 0.10 solar masses (2.83 +/- 0.04 solar masses and 1.85 +/- 0.02 solar masses), which is 30% (100% and 40%) higher than the measured result 1.5 +/- 0.3 solar masses (< 1.4 solar masses and 1.29 + 0.19 / - 0.07 solar masses). For AWOM, where the free parameter of model is the density of star, we infer the NS radii to be around 10 - 20 km for the above three sources, based on which we can infer the matter compositions inside NSs with the help of the equations of state (EOSs). In particular, for SAX J1808.4-3658, AWOM shows a lower mass density of its NS than those of the other known kHz QPO sources, with the radius range of 17 - 20 km, which excludes the strange quark matter inside its star.

With {\it RXTE} data, we analyzed the cross-correlation function between the soft and hard X-rays of the transient atoll source 4U 1608-52. We found anti-correlations in three outbursts occurred in 1998, 2002 and 2010, and significant time lags of several hundreds of seconds in the latter two outbursts. Our results show no correlation between the soft and hard X-rays in the island state, and a dominated positive correlation in the lower banana state. Anti-correlations are presented at the upper banana state for the outburst of 2010 and at the lower left banana states for the other two outbursts. So far for atoll sources the cross-correlation has been studied statistically only for 4U 1735-44, where anti-correlations showed up in the upper banana state. Here our investigation upon 4U 1608-52 provides a similar result in its 2010 outburst. In addition, we notice that the luminosities in the upper banana of 1998 and 2002 outbursts are about 1.5 times that of 2010 outburst whose luminosity in the upper banana is close to that of 4U 1735-44. The results suggest that the states in color-color diagram of a source could be correlated with the luminosity of the source. A further spectral analysis shows that, during the 2010 outburst, although an anti-correlation presents at the highest fluxes, the contemporary spectrum is not the softest one along the outburst evolution. This suggests that the observed anti-correlation may be relevant to the transition between the hard and soft states, which is consistent with the previous results on 4U 1735-44 and several black hole X-ray binaries that anti-correlations are observed during the transition states.

With a sky brightness monitor in Xinglong station of National Astronomical Observatories of China (NAOC), we collected data from 22 dark clear nights and 90 lunar nights. We first measured the sky brightness variation with time in dark nights, found a clear correlation between the sky brightness and human activity. Then with a modified sky brightness model of moon night and data from moon night, we derived the typical value for several important parameters in the model. With these results, we calculated the sky brightness distribution under a given moon condition for Xinglong station. Furthermore, we simulated the moon night sky brightness distribution in a 5 degree field of view telescope (such as LAMOST). These simulations will be helpful to determine the magnitude limit, exposure time as well as the survey design for LAMOST at lunar night.

We present new quasars discovered in the vicinity of the Andromeda and Triangulum galaxies with the LAMOST during the 2010 and 2011 observational seasons. Quasar candidates are selected based on the available SDSS, KPNO 4 m telescope, XSTPS optical, and WISE near infrared photometric data. We present 509 new quasars discovered in a stripe of ~135 sq. deg from M31 to M33 along the Giant Stellar Stream in the 2011 pilot survey datasets, and also 17 new quasars discovered in an area of ~100 sq. deg that covers the central region and the southeastern halo of M31 in the 2010 commissioning datasets. These 526 new quasars have i magnitudes ranging from 15.5 to 20.0, redshifts from 0.1 to 3.2. They represent a significant increase of the number of identified quasars in the vicinity of M31 and M33. There are now 26, 62 and 139 known quasars in this region of the sky with i magnitudes brighter than 17.0, 17.5 and 18.0 respectively, of which 5, 20 and 75 are newly-discovered. These bright quasars provide an invaluable collection with which to probe the kinematics and chemistry of the ISM/IGM in the Local Group of galaxies. A total of 93 quasars are now known with locations within 2.5 deg of M31, of which 73 are newly discovered. Tens of quasars are now known to be located behind the Giant Stellar Stream, and hundreds behind the extended halo and its associated substructures of M31. The much enlarged sample of known quasars in the vicinity of M31 and M33 can potentially be utilized to construct a perfect astrometric reference frame to measure the minute PMs of M31 and M33, along with the PMs of substructures associated with the Local Group of galaxies. Those PMs are some of the most fundamental properties of the Local Group.

We analyze the cross-correlation function between the soft and hard X-rays of atoll source 4U 1735-44 with RXTE data, and find the anti-correlated soft and hard time lags of about hecto-second. On the island state, the observations do not show any obvious correlations, and most observations of banana branch show positive correlation. However, the anti-correlations are detected on the upper banana branch. These results are different from those of Z sources (Cyg X-2, GX 5-1), where the anti-correlation is detected in the low luminosity states, then the lag timescales of both this atoll and Z sources are found to be similar, at the magnitude of several tens to hundreds of seconds. As a comparison, it is noted that the anti-correlated lags of thousand-second have been reported from the several black hole candidates in their intermediate states. Finally, we compare the correspondent results of atoll source 4U 1735-44 with those observed in black hole candidates and Z sources, and the possible origins of the anti-correlated time lags are discussed.

In this paper, the rms-flux (root mean square-flux) relation along the Z-track of the bright Z-Source Cyg X-2 is analyzed using the observational data of Rossi X-ray Timing Explorer (RXTE). Three types of rms-flux relations, i.e. positive, negative, and 'arch'-like correlations are found in different branches. The rms is positively correlated with flux in normal branch (NB), but anti-correlated in the vertical horizontal branch (VHB). The rms-flux relation shows an 'arch'-like shape in the horizontal branch (HB). We also try to explain this phenomenon using existing models.

A direct WIMP (Weakly Interacting Massive Particle) detector with a neutron veto system is designed to better reject neutrons. Two experimental configurations are studied in the present paper: one is for an Xe detector enclosed in a Gd-loaded scintillator and the other one is for an Xe detector placed inside a reactor neutrino detector. The Gd-doped liquid scintillator (or the neutrino detector) is used as a neutron veto device. The neutron backgrounds for the two experimental designs have been estimated using Geant4 simulations. The results show that the neutron backgrounds can decrease to O(0.1) events per year per tonne of liquid Xenon. We calculate the sensitivities to spin-independent WIMP-nucleon elastic scattering. An exposure of one tonne $\times$ year could reach a cross-section of about 6$\times$$10^{-11}$ pb.

The high mass X-ray binary 4U 1901+03 was reported to have the pulse profile evolving with the X-ray luminosity and energy during its outburst in February-July 2003: the pulse peak changed from double to single along with the decreasing luminosity. We have carried out a detailed analysis on the contemporary phase-resolved energy spectrum of 4U 1901+03 as observed by Rossi X-ray Timing Explorer (RXTE). We find that, both the continuum and the pulse spectra are phase dependent. The optical depth derived from the pulse spectrum is in general larger than that from the continuum. Fe Ka emission line is only detected in the spectrum of the continuum and is missing in the pulse spectrum. This suggests an origin of Fe emission from the accretion disk but not the surface of the neutron star.

Based the Alfven wave oscillation model (AWOM) and relativistic precession model (RPM) for twin kHz QPOs, we estimate the emission positions of most detected kHz QPOs to be at r=18+-3 km (R/15km) except Cir X-1 at r = 30\+-5 km (R/15km). For the proposed Keplerian frequency as an upper limit to kHz QPO, the spin effects in Kerr Spacetime are discussed, which have about a 5% (2%) modification for that of the Schwarzchild case for the spin frequency of 1000 (400) Hz.The application to the four typical QPO sources, Cir X-1, Sco X-1, SAX J1808.4-3658 and XTE 1807-294, is mentioned.

We report the detections of the anti-correlated soft and hard X-rays, and the time lags of $\sim$ hecto-second from the neutron star low-mass X-ray binary Cyg X-2, a well-known Z-type luminous source. Both the anti-correlation and the positive correlation were detected during the low-intensity states, while only the latter showed up during high-intensity states. Comparing with the lower part of normal branch and flaring branch, more observations located on the horizontal and the upper normal branches are accompanied with the anti-correlation, implying the occurrence of the anti-correlation under circumstance of a low mass accretion rate. So far the anti-correlated hard lag of thousand-second timescale are only reported from the Galactic black hole candidates in their hard states. Here we provide the first evidence that a similar feature can also establish in a neutron-star system like Cyg X-2. Finally, the possible origins of the observed time lags are discussed under the current LMXB models.

The rms (root mean square) variability is the parameter for understanding the emission temporal properties of X-ray binaries (XRBs) and active galactic nuclei (AGN). The rms-flux relation with Rossi X-ray Timing Explorer (RXTE) data for the dips and nondip of black hole Cyg X-1 has been investigated in this paper. Our results show that there exist the linear rms-flux relations in the frequency range 0.1-10 Hz for the dipping light curve. Moreover, this linear relation still remains during the nondip regime, but with the steeper slope than that of the dipping case in the low energy band. For the high energy band, the slopes of the dipping and nondipping cases are hardly constant within errors. The explanations of the results have been made by means of the ``Propagating Perturbation'' model of Lyubarskii (1997).

We present X-ray spectral analyses of low mass X-ray binary Cir X-1 during X-ray dips, using the Rossi X-ray Timing Explorer (RXTE) data. Each dip was divided into several segments, and the spectrum of each segment was fitted with a three-component blackbody model, in which two components are affected by partial covering and the third one is unaffected. A Gaussian emission line is also included in the spectral model to represent the Fe Ka line at ~ 6.4 keV. The fitted temperatures of the two partially covered components are about 2 keV and 1 keV, while the uncovered component has a temperature of ~ 0.5-0.6 keV. The equivalent blackbody emission radius of the hottest component is the smallest and that of the coolest component is the biggest. During dips, the fluxes of the two hot components are linearly correlated, while that of the third component doesn't show any significant variation. The Fe line flux remains constant within errors during the short dips. However, during the long dips the line flux changes significantly and is positively correlated with the fluxes of the two hot components. These results suggest: (1) the temperature of the X-ray emitting region decreases with radius, (2) the Fe Ka line emitting region is close to the hot continuum emitting region, and (3) the size of the Fe line emitting region is bigger than the size of the obscuring matters causing short dips but smaller than the sizes of those causing long dips.

We studied the correlations between spin frequencies and kilohertz quasi-periodic oscillations (kHz QPOs) in neutron star low mass X-ray binaries. The updated data of kHz QPOs and spin frequencies are statistically analyzed. We found that when two simultaneous kHz QPOs are present in the power spectrum, the minimum frequency of upper kHz QPO is at least 1.3 times larger than the spin frequency, i.e. \nu_{s}<\nu_{2min}/1.3. We also found that the average kHz QPO peak separation in 6 Atoll sources anti-correlates with the spin frequency in the form \lan\dn\ran = -(0.19\pm0.05)\ns+(389.40\pm21.67)Hz. If we shifted this correlation in the direction of the peak separation by a factor of 1.5, this correlation matches the data points of the two accretion powered millisecond X-ray pulsars, SAX J1808.4-3658 and XTE J1807-294.

The recombination history of the Universe provides a useful tool for constraining the annihilation of dark matter particles. Even a small fraction of dark matter particles annihilated during the cosmic dark age can provide sufficient energy to affect the ionization state of the baryonic gas. Although this effect is too small for neutralinos, lighter dark matter particle candidates, e.g. with mass of 1-100 MeV, which was proposed recently to explain the observed excess of positrons in the Galactic Center, may generate observable differences in the cosmic microwave background (CMB) temperature and polarization anisotropies. The annihilations at the era of recombination affects mainly the CMB anisotropy at small angular scales (large $\ell$), and is distinctively different from the effect of early reionization. We perform a multi-parameter analysis of the CMB data, including the both the WMAP first year and three year data, and the ACBAR, Boomerang, CBI, and VSA data. Assuming that the observed excess of $e^+e^-$ pairs in the galactic center region is produced by dark matter annihilation, and that a sizable fraction of the energy produced in the annihilation is deposited in the baryonic gas during recombination, we obtain a %95 dark matter mass limit of $M<8 \MeV$ with the current data set.

This paper presents the simulation of the GLAST high energy gamma-ray telescope. The simulation package, written in C++, is based on the Geant4 toolkit, and it is integrated into a general framework used to process events. A detailed simulation of the electronic signals inside Silicon detectors has been provided and it is used for the particle tracking, which is handled by a dedicated software. A unique repository for the geometrical description of the detector has been realized using the XML language and a C++ library to access this information has been designed and implemented. A new event display based on the HepRep protocol was implemented. The full simulation was used to simulate a full week of GLAST high energy gamma-ray observations. This paper outlines the contribution developed by the Italian GLAST software group.