• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    The Study of Variations of Cyclotron Resonant Scattering Features in GX 301-2 by using the INTEGRAL/IBIS

    2016-07-11 01:35:31YuGuangwenWangWei
    天文研究與技術 2016年1期
    關鍵詞:王偉雙星磁場強度

    Yu Guangwen, Wang Wei

    (1. National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China, Email: wangwei@bao.ac.cn;2. University of Chinese Academy of Sciences, Beijing 100049, China)

    ?

    The Study of Variations of Cyclotron Resonant Scattering Features in GX 301-2 by using the INTEGRAL/IBIS

    Yu Guangwen1,2, Wang Wei1

    (1. National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China, Email: wangwei@bao.ac.cn;2. University of Chinese Academy of Sciences, Beijing 100049, China)

    Abstract:With long-term hard X-ray monitoring, observations on high mass X-ray binary GX 301-2 from 2003-2011 performed by INTEGRAL/IBIS, we first systematically studied its spectral properties in different accreting luminosities and orbital phases. The cyclotron resonant scattering feature (CRSF) at energies from 35-47 keV was detected in the hard X-ray spectra, suggesting the magnetic field of 5×1012Gs. The variations of the CRSF show no relation to the X-ray luminosity while the line centroid energy of the CRSF has a positive correlation to the photon index and spectral cutoff energy, and there also exists a weak correlation between the absorption depth and cutoff energy. These relations support the idea that the spectral cutoff in accreting X-ray pulsars is strongly affected by the cyclotron resonant scattering. The correlation between the ratio of the line width to the centroid energy and absorption depth implies a tall cylindrical accreting column on the surface of the neutron star in GX 301-2. The explanation of the long spin period in GX 301-2 requires a strong surface magnetic field at least higher than 2×1014Gs which is contradict with the measured value based on the cyclotron line energy. The line-forming region with a height at least larger than two neutron star radii is proposed to resolve the contradiction. This scenario is supported by the tall cylindrical accreting column on the neutron star surface in GX 301-2 according to the variation patterns of the CRSF. In addition, the possible evolution scenario of accreting magnetars like GX 301-2 is briefly discussed, and GX 301-2 would be an accreting magnetar in the equivalence phase.

    Key words:GX 301-2; X-rays binaries; Neutron; Magnetic fields

    1Introduction

    The high mass X-ray binaries (HMXBs) are X-ray sources composed of a nearly-type massive star and an accretion compact object (the neutron star or black hole). A major part of HMXBs which harbor X-ray pulsars are believed to be the magnetic neutron stars with a relatively strong magnetic field of 1012Gs[1]. Their X-ray emission can be powered by the accretion disk (disk-fed) or direct accretion from the radiative wind (wind-fed) of the donor star. The high energy luminosity variation is generally due to the changes of the orbital phase and accreting rates.

    GX 301-2 (4U 1223-62) is a high-mass X-ray binary system consisting of a neutron star and an early B-type optical companion Wray 977. Its distance is quite uncertain, with estimates between 1.8 kpc[2]and 5.3 kpc[3]. The latest estimate is about 3 kpc determined by the interstellar absorption[4]. The orbital period of the system is 41.5 d with a high eccentricity of 0.5[5-6]. The neutron star is an accreting X-ray pulsar with a spin period of 680-700 s[7-8]. GX 301-2 is a variable X-ray source whose X-ray luminosity powered by the wind-fed direct accreting process which changes with the orbital phases.

    The cyclotron absorption line feature of GX 301-2 around 37 keV was first discovered by Makishima & Mihara (1992)[9]based on the Ginga observations, and then confirmed by Orlandini (2000)[10]using the BeppoSAX observations. This Cyclotron Resonant Scattering Feature (CRSF) is produced by the resonant scattering of photons by electrons whose energies are quantized into Landau levels by the magnetic fields[11]. Therefore cyclotron resonance scattering features can be used to directly measure the magnetic field of the neutron star, and provide a useful way to understand the accretion geometry and physics near the surface of the neutron star. However, the cyclotron line centroid energy in GX 301-2 shows large variations and deviations by different measurements. Based on the RXTE data, Kreykenbohm et al. (2004)[12]reported the variable cyclotron line with energies from 30-38 keV according to the phase-resolved spectra. La Barbera et al. (2005)[13]detected the CRSFs with the centroid energies from 45-53 keV in different orbital phases with the BeppoSAX observations. Doroshenko et al. (2010)[6]detected a cyclotron scattering line at 46 keV in GX 301-2 with early INTEGRAL observations. Recently Suzaku broadband spectroscopy on GX 301-2 showed the variations of CRSFs from 28-43 keV[14], which is still consistent with early Ginga and RXTE results, but lower than the reported centroid energies by BeppoSAX and INTEGRAL. This discrepancy should result from that previous BeppoSAX and INTEGRAL results used a Gaussian line profile to fit the absorption dips, while other work applied the cyclotron absorption line profile into the data.

    The CRSF centroid energy in GX 301-2 shows a large variation near 25 keV according to previous different observations, which makes it the largest variation scale in all accreting X-ray pulsars with CRSF measurements. Previous different missions reported different centroid energy ranges for CRSFs in GX 301-2 using different line profile models. Then more observations need systematical studies of the centroid energy variation range in GX 301-2 with the same spectral line model in fittings. In addition, what are this variation pattern of CRSFs and the relationship of CRSFs to continuum spectral properties in the supergiant X-ray binary GX 301-2? In the Be X-ray transient and the low-mass X-ray binary Her X-1, the correlations between CRSF centroid energies and X-ray luminosities are discovered[15-17]. GX 301-2 has the different companion star and the different accreting process (the wind-fed direct accretion) from the Be X-ray transient and the low-mass X-ray binary. Thus the relation between the CRSF and the continuum spectral property in GX 301-2 should be studied in detail, which would help to understand the accreting processes and the geometry in supergiant X-ray binaries.

    Using the archived data of GX 301-2 collected by the INTEGRAL observations in 2003-2011, we mainly concentrated on searching for the CRSFs in the hard X-ray spectra of 18-100 keV and studying the correlations between the CRSF and the continuum spectral parameter. In Section.2 of this paper, we show the information of the INTEGRAL observations and briefly introduce the data analysis. In Section.3, we study the spectral properties and search for CRSFs of GX 301-2 in different luminosity ranges. A cyclotron absorption line model cyclabs is applied in our analysis. The study of the possible correlation between the CRSF and the spectral parameter are presented for GX 301-2 in Section.4. We also study the spectral parameter and the CRSF versus the orbital phase in Section.5. We summarize our results in the last section Section.6 and discuss the implication of the results on the nature of GX 301-2.

    2INTEGRAL observations

    Fig.1The IBIS count rate variations of GX 301-2 in the energy range of 20-40 keV from 2003-2011. GX 301-2 is a highly variable source with the mean detected IBIS rate from 9-180 cts/s

    The INTEGRAL is the currently operational space-based hard X-ray/soft gamma-ray telescope by ESA[24]. GX 301-2 is observed with the frequent pointing observational surveys on the Carina region from 2003-2011 by the INTEGRAL satellite. The hard X-ray data are obtained by the INTEGRAL Soft Gamma-Ray Imager[25]which has a 12′ (FWHM)angular resolution and arcmin source location accuracy in the energy band of 15-200 keV. The JEM-X aboard on the INTEGRAL is a small X-ray telescope covering about 4-30 keV, which could help to constrain the continuum of IBIS hard X-ray sources in low energy bands. However, JEM-X has a small field of view (off-source axis angle 5°), cannot detect GX 301-2 in most observational revolutions.

    In this work, we used the available archival data for the IBIS observations where GX 301-2 is within ~10 degrees of the pointing direction with the corrected on-source time longer than 15 ks for each observation. In Table 1, the information of available INTEGRAL observations used in this paper is summarized. The INTEGRAL observations provided a long-term monitoring on GX 301-2 from 2003-2011. The archival data used in our work are available from the INTEGRAL Science Data Center (ISDC). The analysis is made with the standard INTEGRAL off-line scientific analysis (OSA[26]) software, version 10. The individual pointing in each satellite revolution (about 3 days) processed with OSA 10 is mosaicked to create sky images for the source detection. We have used the energy range of 20-40 keV by the IBIS for the source detection, significance levels and the quoted source fluxes for each revolution (see Table 1).

    From Table 1, we can see that the source is variable in hard X-ray bands, which is typical for these class of wind-fed accreting X-ray binaries. The mean IBIS count rate of GX 301-2 (20-40 keV) in each revolution has a large variation range from <10 cts/s to ~180 cts/s. The spectral properties of GX 301-2 would vary with the different accreting luminosities and orbital phases. In following spectral analysis, we will extract the hard X-ray spectrum of GX 301-2 from 18-100 keV for each revolution. The spectral analysis processes then are made on these spectra to probe the variations of spectral properties of GX 301-2 in different luminosity ranges.

    3Hard X-ray spectral properties of GX 301-2

    We have extracted the X-ray spectra from 18-100 keV by the IBIS in all observational revolutions for GX 301-2. These spectra will be fitted with a suitable continuum model, and then we will try to fit the cyclotron absorption line, deriving both the continuum spectral parameters and cyclotron line parameters. With frequent observations by the IBIS, covering a wide luminosity range and different orbital phases, we can study the variation patterns of the cyclotron scattering line feature.

    Table 1INTEGRAL/IBIS observations of the field around GX 301-2 within the off-axis angle 10°. The time intervals of

    observations in the revolution number and the corresponding MJD dates, and the corrected on-source exposure

    time are listed. Mean IBIS count rate in the energy range of 20-40 keV was also shown

    The IBIS covers an energy range from 18-100 keV. A small telescope, the JEM-X covering lower energy ranges from 4-30 keV could also detected the same source but requiring a much smaller off-axis angle. However, we have checked all available JEM-X data, only in a few revolutions we can extract a JEM-X spectrum for GX 301-2. A broader energy band spectrum from 4-100 keV could better constrain the continuum of GX 301-2. However, only very limited number of JEM-X spectra can be used to carried out statistical studies on spectral properties in the followings, for which we mainly used the IBIS data.

    Here we will show one spectrum sample obtained by both the JEM-X and IBIS data to present the spectral properties of GX 301-2 in hard X-rays. This spectrum is derived from the revolution 176 when both JEM-X and IBIS detectors have detected GX 301-2 with a good significance level. In Fig.2, the combined spectrum from 4-100 keV obtained by both JEM-X and IBIS for GX 301-2 is presented. The cross-calibration studies on the JEM-X and IBIS/ISGRI detectors have been made using the Crab observation data, and the calibration between the JEM-X and the IBIS can be good enough within ~ 6% (see samples in [21-22]). In the spectral fittings, the constant factor between the JEM-X and the IBIS is set to be one. The spectral analysis software package used is the XSPEC 12.6.0q.

    Generally, the hard X-ray spectrum of accreting X-ray pulsars like GX 301-2 can be described by a power-law model plus a high energy exponential roll off (hereafter cutoffpl):A(E)=KE-Γexp (-E/Ecutoff). Other simple spectral models like the single power-law model, the thermal bremsstrahlung model, a blackbody plus a power-law model are also applied to fit the spectra of this source. However, these models cannot fit the hard X-ray spectra of GX 301-2 better than the cutoffpl model for the spectra from 4-100 keV, also for the cases from 18-100 keV (i.e., the larger reduced value ofχ2), therefor we just use the cutoff power-law model in the following spectral analysis.

    We also fitted the spectrum from 18-100 keV only obtained by the IBIS. The fitted parameters can be compared with the ones using the combined spectrum from 4-100 keV. In the top panel of Fig.3, the spectrum is fitted with the cutoff power-law model. The absorption dip from 30-40 keV appears in the residuals. In the bottom panel, the cyclotron line feature is added to fit the spectrum. It is seen that no significant features in the residuals can be found.

    From the best fitted parameters of the spectra for two energy bands shown in Table 2, we can compare the spectral fitting results of the combined data and only the IBIS data. For the one revolution, these two data sets give the similar spectral fitting results, for both continuum and cyclotron scattering features. Therefore, the spectral fitting parameters using only the IBIS data can well describe the hard X-ray properties of GX 301-2, specially the parameters of the CRSF. Then we use the cutoff power-law model together with a cyclotron scattering line to fit the spectra derived in all available observational revolutions. The continuum spectral parameters and CRSF parameters are collected together for the next correlation studies.

    Fig.2The hard X-ray spectrum of GX 301-2 from 4-100 keV obtained by both JEM-X (low energy) and IBIS (high energy) for the observation revolution 176. In the Top, the spectrum is fitted with a power-law model anda high energy cutoff (cutoffpl). In the Middle, the spectrum is fitted with the cutoffpl plus a Gaussian line with the fixed Fe line energy at 6.4 keV. A strong absorption line feature from 30-50 keV is detected in residuals. In the Bottom, the spectrum is fitted with the cutoff using the cyclotron line model cyclabs plus a Gaussian line with the fixed Fe line energy at 6.4 keV. See the text for details

    Table 2The hard X-ray spectral properties (according to spectrum samples in Figs.2-3) of the high mass X-ray

    binary GX 301-2. The flux is given in the unit of 10-10erg cm-2s-1within the range of 18-100 keV.

    The iron line fluxKis in the unit of 10-5photons cm-2s-1

    In different accreting luminosities, one cyclotron scattering line feature is detected by the present INTEGRAL/IBIS observations. The detected line centroid energy in GX 301-2 varies from 35-47 keV according to all the available data in different accreting luminosities (also see Fig.4). This energy range is consistent with the early Ginga, RXTE, and Suzaku results. The variation patterns of the CRSF may be related to the continuum spectral properties, which will be discussed in the next section.

    4Correlation studies of cyclotron resonance spectral features

    The variations of CRSFs may be induced by some physical parameters in accreting X-ray pulsars, which can be statistically studied with the relations between the CRSFs and spectral properties. The correlations between the fundamental line energy of the CRSFs and the X-ray luminosity are reported in some X-ray pulsars, like the Be X-ray transients 4U 0115+63[15]and GX 304-1[16]during their outbursts, and a low mass X-ray binary Her X-1[17]. A positive correlation between the energy of the fundamental line and the photon index is discovered in 4U 0115+63 during the giant outburst in 2008[15]. In addition, with collecting observation data for different X-ray pulsars, a positive correlation between the energy of the fundamental CRSF and the cutoff energy in the accreting X-ray pulsar continuum is found[28-30]. It is suggested that the spectral cutoff in accreting X-ray pulsars should be strongly affected by the cyclotron resonant scattering.

    Fig.3The hard X-ray spectrum of GX 301-2 from 18-100 keV obtained by the IBIS for the observation revolution 176. In the Top, the spectrum is fitted with a power-law model and a high energy cutoff (cutoffpl).In the Bottom, the spectrum is fitted with the cutoffpl using the cyclotron line model cyclabs. See the text for details

    Fig.4The derived cyclotron energyEcand the depth of the cyclotron line versus the three spectral parameters: the hard X-ray flux in the range of 18-100 keV, the photon indexΓand the exponential cutoff energyEcutoff

    The production of cyclotron resonant absorption line features near the surface of the neutron star is physically complicated, which will sensitively depend on the accretion states and the accretion geometry (e.g.,[31-32]). Thus the detection of cyclotron scattering line features is not only used to identify a magnetized neutron star in binary systems but also important for studies on the accretion process and the geometry near the surface of a neutron star[32]. The variation pattern of the CRSFs and their physical origins need more studies with more observations in detail. With the derived spectral parameters of GX 301-2 by the IBIS long-term observations, we can first systematically study the variation patterns of the CRSF and probe the accreting geometry in the wind-fed accreting supergiant binary GX 301-2. In the following, we will use the Pearson product-moment correlation coefficient (hereafterr) to describe the correlation between spectral parameters, where -10.5.

    In Fig.4, we have presented the relationships among the centroid energyEc, the depth of the CRSF and three spectral parameters: hard X-ray flux from 18-100 keV, photon indexΓand cutoff energyEcutoff. The centroid energy of the CRSF has no relation to the X-ray flux in GX 301-2, which is different from the Be/X-ray pulsars 4U 0115+63 and GX 304-1. But the line energy has the positive correlations to the photon indexΓ(r=0.609) and the cutoff energyEcutoff(r=0.687). The relation ofEcversusΓis similar to the case in the Be/X-ray pulsar 4U 0115+63[15]. The relation ofEcversusEcutoffis consistent with the common behavior by collecting different accreting X-ray pulsars[28-30]. These correlations support the conclusion that the spectral properties of accreting X-ray pulsars in hard X-ray bands will be affected by the CRSFs.

    The depth of the CRSFs is also changed, however, it seems to have no relation to the X-ray flux. While the depth still show a weak positive correlation to the cutoff energyEcutoff(r=0.509), but no significant relationship to the spectral indexΓ(r=0.114). Hence the spectral cutoff energyEcutoffhas the relationship to both the energy and the depth of the CRSF in GX 301-2. Then the high energy cutoff in the spectrum of the accreting X-ray pulsars GX 301-2 should be mainly induced by the cyclotron resonant scattering.

    In the searching for the CRSFs in different accreting X-ray pulsars with the RXTE, Coburn et al. (2002)[30]also found a correlation between the absorption depth and the ratio of the width to the energy in the phase averaged spectra of the accreting X-ray pulsars. This relation can help to understand the accretion geometry of accreting X-ray pulsars. In an accreting column, one can find[33]: Width/Ec∝(kTecos2θ)1/2, wherekTeas the electron temperature along the magnetic fields is nearly constant,θis the angle between the line of sight and the magnetic field. Generally, we assume that there exist two classes of accreting geometry patterns: a flat coin shape and a tall cylindrical shape. In the flat coin shape, the depth of CRSFs is largest when the line of sight is perpendicular to the direction of the magnetic field[34]. While in a tall cylinder shape, the depth is largest when the line of sight is parallel to the magnetic field[34], which then predicts a positive correlation between the depth and the ratio of the width to the energy. Kreykenbohm et al. (2004)[12]analyzed the RXTE data, derived this relation in the single source GX 301-2. With the present INTEGRAL observations, we derived the variation of the CRSFs in GX 301-2, then presented the relation between the ratio of the width to the energy and the absorption depth in Fig.5. Though large error bars exist, the correlation between the ratio and depth (r=0.837) is confirmed by our analysis with a larger data set from the IBIS. With the properties of the CRSF and its variation pattern, it is suggested that GX 301-2 has a cylindrical column accretion geometry.

    5The CRSFS versus the orbital phase

    Fig.5The derived Width/Ecversus the depth of the cyclotron line for GX 301-2

    GX 301-2 has an orbital period of ~41.5 days, and its X-ray luminosity changes in one binary orbit. The INTEGRAL/IBIS performed a long-term observations on this source, covering the different orbital phases. Thus we can derive the orbital phase resolved spectra for GX 301-2, and study the variation of spectral properties versus the orbital phase.

    Using the orbital ephemeris provided by Koh et al. (1997)[5], we divided one orbital circle into fourteen phase bins, and re-combined the INTEGRAL data (in units of the science windows, each has a duration of ~2000-4000 s) according to the defined phase bins. The data analysis processes were re-made for all data in each orbital phase,so that we can extract the hard X-ray spectrum of GX 301-2 for each orbital phase bin. Similar to the spectral fitting processes in Section 3, we have derived the best fitted continuum spectral parameters, like the X-ray flux,ΓandEcutoff, the centroid energies and the absorption depth of the CRSF in GX 301-2. In Fig.6, we have presented the different spectral parameters versus the orbital phase.

    The Hard X-ray flux variation shows the flare around the orbital phase 0.9, which is still consistent with the orbital profile derived by the RXTE/ASM data in the range of 1.5-12 keV[35]. This flare occurs around the orbital phase just before the periastron passage of the neutron star in GX 301-2. The other spectral parameters also show the variations over the orbital phase.Specially, the mean centroid energy of the CRSF over the orbit varies from 37-44 keV, which is consistent with the variation range derived in different luminosities presented in Section 3. The flux variation has no correlations to the changes of other spectral parameters. But the centroid energy and depth of the CRSF still show the similar variation profiles to those of the photon indexΓand the spectral cutoff energyEcutoffover the whole orbital phase. These correlations are similar to the relations found in Section 3. Hence the relations between the CRSF and continuum spectral parameters in GX 301-2 are confirmed with the orbital resolved spectral analysis.

    Fig.6The derived cyclotron line energyEc, the depth of the cyclotron line and the three spectral parameters: the hard X-ray flux in the range of 18-100 keV, the photon indexΓand the exponential cutoff energyEcutoffversus the orbital phase

    6Summary and discussion

    In this work, we have carried out the data analysis on GX 301-2 with the INTEGRAL/IBIS long-term observations. The cyclotron scattering line at energies from 35-47 keV in the X-ray spectra of GX 301-2 is confirmed by our analysis. The variation range of the line centroid energy is still consistent with the previous reports by Ginga[9], RXTE[12], Suzaku[14].

    The variation patterns of the CRSFs and spectral properties of GX 301-2 are systematically studied over different luminosity ranges and orbital phases. The centroid energy of the CRSF has no relation to the X-ray luminosity, which is different from other accreting X-ray pulsars, like the Be X-ray transients. While the line centroid energy shows the correlation to the photon index and the spectral cutoff energy, and the absorption depth also has a weak correlation to the cutoff energy. These relations support that the spectral profile, specially, the high energy cutoff is strongly affected by the cyclotron resonant scattering in GX 301-2. This characteristic is similar to other accreting X-ray pulsars, suggesting that there exist the similar radiation processes in these X-ray pulsars independent of their companion stars and accreting processes. In addition, a correlation between the ratio of the line width to the centroid energy and absorption depth in GX 301-2 implies a tall cylinder accreting column on the surface of the magnetized neutron star. This special accreting geometry could help us to understand the nature of the neutron star in GX 301-2.

    As shown in the introduction, GX 301-2 could be an accreting magnetar candidate. Since the discovery of the GX 301-2[36], the spin period changes dramatically, and shows both spin-down and spin-up trends in the last thirty years (see [12,14]). Therefore GX 301-2 is located in a torque equivalence phase, hence its equilibrium spin period can be estimated with the surface magnetic field of the neutron star and its average luminosity (the accreting rate). According to a standard model of the accreting disk[18], a critical period is defined by equating the corotational radius of the neutron star to the magnetospheric radius, which gives

    where κ~0.5-1isthegeometricalparameteroftheaccretionflow, μisthedipolemagneticmomentofthe

    法律應當把公共信息生產和服務作為監(jiān)管者的主要職責。這主要包括提供民間融資市場總體供需信息、風險信息和信用信息,以減少融資主體信息盲區(qū),降低信息獲取成本。微觀監(jiān)管作為重要手段雖不可或缺,但須限制在較小范圍內,包括適當?shù)臏嗜肟刂坪秃弦?guī)性運行監(jiān)管。

    However, if assuming that the cyclotron scattering line comes from the surface of the neutron star, the magnetic field of the star can be estimated by

    WhereEcycis the energy of the fundamental line. Here we can take the maximum measured valuesEcyc=47 keV, andz~0.3, the gravitational redshift near the surface of the neutron star. Then we derived the magnetic field of ~5×1012Gsfor GX 301-2 which is much lower than the magnetic field range for common magnetar candidates. Why a large discrepancy in observations and theoretical expectations exists?

    It has been suggested that the line-forming region may reside in an accretion column with a significant height[37,6]above the surface of the neutron star. Then the measured magnetic field based on the cyclotron line energy is just the magnetic field strength near the line-forming region rather than the surface. SinceBcyc=(R/RNS)-3Bsuff, whereRis the radius of the line-forming region. Then assumingBsurf~2×1014Gs, andBcyc~4×1012Gs, one can derive a radius of line emission regionR~3.7 RNS. If the surface magnetic field is near 1015Gs, the requested emission region will be higher. Thus the emission region must be located in a height at least larger than 2RNS. The accreting and emission region will be like a cylinder column with a height >>2RNSin GX 301-2. This geometry is supported by our analysis results on GX 301-2: the correlation between the ratio of the line width to the centroid energy and the absorption depth implies a tall cylinder accreting column on the surface of the highly magnetized neutron star.

    Then the high mass X-ray binary GX 301-2 could be very special, and can belong to a new class of compact objects—accreting magnetars. Accreting magnetars have the quite different properties from the isolated magnetar candidates (see discussions in [21]), anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs), which are the isolated systems without companion stars. Of course, the physical origin of magnetars in both accreting and isolated systems is unknown. An evolution scenario for accreting magnetars has been suggested[21]. Magnetars born in the binaries could spin down faster in the propeller phase than other normal neutron star binaries (fast spin-down trends observed in two systems, 4U 2206+54 and SXP 1062), then its spin period could spin down to slower than 104seconds. After the spin-down phase, the magnetar candidate can transfer into the spin-up process (a long-term spin-up trend is found in 2S 0114+65[22]). This spin-up process will stop till the accreting magnetar is located around the equilibrium spin period range (several hundred seconds). In this scenario, GX 301-2 is the successor of the younger and fast-evolved superslow pulsation X-ray pulsars (like4U 2206+54 and 2S 0114+65), it has become an older accreting magnetar in the torque equivalence phase. Hence GX 301-2 becomes interesting and very unique in the family of accreting magnetars.

    Though the idea of accreting magnetars is proposed in some literatures ([21] and references therein), we need understand the nature of the magnetar definition. In Table 3, we show the present understanding of both isolated and accreting magnetars in observations and theories. In observations both candidates have longer spin period values, larger spin period derivative and younger ages relative to other normal systems. We may need to search for more common features in two classes of magnetars. The well-known magnetar candidates as AXPs and SGRs require aultra-strong magnetic field to produce X-ray bursts, which can reflect magnetic field activity. How to detect the signals of magnetic field activity in accreting magnetar candidates is an intriguing question in studying accreting magnetars in the future. The Chinese hard X-ray telescope-Hard X-ray Modulation Telescope (HXMT) will be launched in the next year, one of main scientific objects will be the monitoring magnetar candidates, specially, their spectral and temporal variations. Further studies on GX 301-2, specially, on the spectral features and flares, even bursts could improve our understanding on the nature of accreting magnetars and the magnetar family.

    Acknowledgements: We are grateful to the referee for the comments to improve the manuscript.This work is based on observations of INTEGRAL, an ESA project with instrument and science data center funded by ESA member states.

    References:

    [1]Bildsten L, Chakrabarty D, Chiu J, et al. Observations of accreting pulsars[J]. The Astrophysical Journal Supplement, 1997, 113(2): 367-408.

    [2]Parkes G E, Culhane J L, Mason K O, et al. A spectral study of Wray 977, the optical counterpart of the binary X-ray pulsar 4U 1223-62[J]. Monthly Notices of the Royal Astronomical Society, 1980, 191(3): 547-558.

    [4]Kaper L, van der Meer A, Najarro F. VLT/UVES spectroscopy of Wray 977, the hypergiant companion to the X-ray pulsar GX301-2[J]. Astronomy and Astrophysics, 2006, 457(2): 595-610.

    [5]Koh D T, Bildsten L, Chakrabarty D, et al. Rapid Spin-Up Episodes in the Wind-fed Accreting Pulsar GX 301-2[J]. The Astrophysical Journal, 1997, 479(2): 933-947.

    [6]Doroshenko V, Santangelo A, Suleimanov V, et al. Is there a highly magnetized neutron star in GX 301-2?[J]. Astronomy and Astrophysics, 2010, 515: A10-A19.

    [7]White N E, Mason K O, Huckle H E, et al. Periodic modulation of three galactic X-ray sources [J]. Astrophysical Journal, 1976, 209(1): L119-L124.

    [8]Evangelista Y, Feroci M, Costa E, et al. Temporal properties of GX 301-2 over a year-long observation with superAGILE[J]. The Astrophysical Journal, 2010, 708(2): 1663-1673.

    [9]Makishima K, Mihara T. Magnetic fields of neutron stars[C]. Tokyo Japan: Frontiers Science Series, 1992: 23.

    [10]Orlandini M, dal Fiume D, Frontera F, et al. BeppoSAX observations of an orbital cycle of the X-ray binary pulsar GX 301-2[J]. Advances in Space Research, 2000, 25(3-4): 417-420.

    [11]Mészáros P. High Energy Radiation from Magnetized Neutron Stars[M]. Chicago: University of Chicago Press, 1992.

    [12]Kreykenbohm I, Wilms J, Coburn W, et al. The variable cyclotron line in GX 301-2[J]. Astronomy and Astrophysics, 2004, 427(3): 975-986.

    [13]La Barbera A, Segreto A, Sanangelo A, et al. A study of an orbital cycle of GX 301-2 observed by BeppoSAX[J]. Astronomy and Astrophysics, 2005, 438(2): 617-632.

    [14]Suchy S, Fürst F, Pottschmidt K, et al. Broadband spectroscopy using two Suzaku observations of the HMXB GX 301-2[J]. The Astrophysical Journal, 2012, 745(2): 124-135.

    [15]Li Jun, Wang Wei, Zhao Yongheng. Cyclotron resonance energies and orbital elements of accretion pulsar 4U 0115+63 during the giant outburst in 2008[J].Monthly Notices of the Royal Astronomical Society, 2012, 423(3): 2854-2867.

    [16]Klochkov D, Doroshenko V, Santangelo A, et al. Outburst of GX 304-1 monitored with INTEGRAL: positive correlation between the cyclotron line energy and flux[J]. Astronomy & Astrophysics, 2012, 542(4): L28-L31.

    [17]Vasco D, Klochkov D, Staubert R. Hercules X-1: the positive cyclotron line energy/luminosity correlation[J]. Astronomy & Astrophysics, 2011, 532: A99-A104.

    [18]Ghosh P, Lamb F K. Accretion by rotating magnetic neutron stars. III-accretion torques and period changes in pulsating X-ray sources[J]. Astrophysical Journal, 1979, 234(15): 296-316.

    [19]Urpin V, Konenkov D, Geppert U. Evolution of neutron stars in high-mass X-ray binaries[J]. Monthly Notices of the Royal Astronomical Society, 1998, 99(1): 73-77.

    [20]Wang Wei. Evidence for a magnetic neutron star in high-mass X-ray binary 4U 2206+54 with INTEGRAL/IBIS observations[J]. Monthly Notices of the Royal Astronomical Society, 2009, 398(3): 1428-1434.

    [21]Wang Wei. Spin and spectral variations of a peculiar high-mass X-ray binary 4U 2206+54[J]. Monthly Notices of the Royal Astronomical Society, 2013, 432(2): 954-966.

    [22]Wang Wei. Long-term hard X-ray monitoring of 2S 0114+65 with INTEGRAL/IBIS[J]. Monthly Notices of the Royal Astronomical Society, 2011, 413(2): 1083-1098.

    [23]Haberl F, Sturm R, Filipovic M D, et al. SXP 1062, a young Be X-ray binary pulsar with long spin period; Implications for the neutron star birth spin[J]. Astronomy & Astrophysics, 2012, 537(1): L1-L5.

    [24]Winkler C, Courvoisier T J L, Di Cocco G, et al. The INTEGRAL mission[J]. Astronomy and Astrophysics, 2003, 411(1): L1-L6.

    [25]Leburn F, Lerav J P, Lavocat P, et al. ISGRI: The INTEGRAL soft gamma-ray imager[J]. Astronomy and Astrophysics, 2003, 411(1): L141-L148.

    [26]Goldwurm A, David P, Foschini L, et al. The INTEGRAL/IBIS scientific data analysis[J]. Astronomy and Astrophysics, 2003, 411(1): L223-L229.

    [27]Mihara T, Makishima K, Ohashi T, et al. New observations of the cyclotron absorption feature in Hercules X-1[J]. Nature, 1990, 346: 250-252.

    [28]Makishima K, Mihara T, Nagase F, et al. Cyclotron resonance effects in two binary X-ray pulsars and the evolution of neutron star magnetic fields[J]. The Astrophysical Journal, 1999, 525(2): 978-994.

    [29]Kreykenbohm I, Coburn W, Wilms J, et al. Confirmation of two cyclotron lines in Vela X-1[J]. Astronomy and Astrophysics, 2002, 395(1): 129-140.

    [30]Coburn W, Heindl W A, Rothschild R E, et al. Magnetic fields of accreting X-ray pulsars with the Rossi X-ray timing explorer[J]. The Astrophysical Journal, 2002, 580(1): 394-412.

    [31]Araya R A, Harding A K. Cyclotron line features from near-critical magnetic fields: the effect of optical depth and plasma geometry[J]. The Astrophysical Journal, 1999, 517(1): 334-354.

    [32]Nishimura O. Superposition of cyclotron lines in accreting X-ray pulsars. I. long spin period[J]. The Astrophysical Journal, 2011, 730(2): 106-116.

    [33]Meszaros P, Nagel W. X-ray pulsar models. I-angle-dependent cyclotron line formation and comptonization[J]. The Astrophysical Journal, 1985, 298: 147-160.

    [34]Isenberg M, Lamb D Q, Wang J C L. Effects of the geometry of the line-forming region on the properties of cyclotron resonant scattering lines[J]. The Astrophysical Journal, 1998, 05(2): 688-714.

    [35]Leahy D A. The RXTE/ASM X-ray light curve of GX301-2[J]. Astronomy and Astrophysics, 2002, 391(1): 219-224.

    [36]Swank J H, Becker R H, Boldt E A, et al. Variations in the spectra and pulse period of GX 301-2[J]. The Astrophysical Journal, 1976, 209: L57-L60.

    [37]Basko M M, Sunyaev R A. The limiting luminosity of accreting neutron stars with magnetic fields[J]. Monthly Notices of the Royal Astronomical Society, 1976, 175(2): 395-417.

    CN 53-1189/PISSN 1672-7673

    用INTEGRAL/IBIS衛(wèi)星研究GX 301-2中回旋共振散射線性質

    于光雯1,2,王偉1

    (1. 中國科學院國家天文臺, 北京100012; 2. 中國科學院大學, 北京100049)

    關鍵詞:GX 301-2; X射線雙星; 中子星; 磁場

    摘要:利用INTEGRAL/IBIS衛(wèi)星2003~2011年對X射線雙星GX 301-2的觀測數(shù)據,系統(tǒng)研究了GX 301-2在不同吸積光度和軌道相位的譜性質,得到GX 301-2的回旋吸收線的性質與X射線光度無關,其能量在35~47 keV之間,且能量與光指數(shù)和截斷能量正相關,回旋吸收線的深度與截斷能量有弱正相關,這表明在吸積X射線脈沖星中,回旋共振散射線對截斷能量有很大影響。而回旋吸收線寬度與能量比值和深度正相關,這表明在GX 301-2表面存在一個吸積柱。由GX 301-2的自旋周期得到表面磁場強度大約為2 × 1014Gs,與通過回旋吸收線得到的磁場強度不符,但可以用吸積柱模型很好地解釋。

    基金項目:國家自然科學基金 (11073030) 資助.

    收稿日期:2015-04-01;修定日期:2015-04-26

    作者簡介:于光雯,女,碩士. 研究方向:高能天體物理. Email: gwyu@bao.ac.cn通訊作者:王偉,男,研究員. 研究方向:高能天體物理. Email: wangwei@bao.ac.cn

    中圖分類號:O572.13

    文獻標識碼:A

    文章編號:1672-7673(2016)01-0011-13

    猜你喜歡
    王偉雙星磁場強度
    CHARACTERIZATION OF RESIDUATED LATTICES VIA MULTIPLIERS*
    雙星啟示錄
    李雙星 一心為民拔“窮根”
    河北畫報(2020年21期)2020-12-14 03:16:08
    小 蝌 蚪 的 尾 巴
    關于醫(yī)用磁共振成像系統(tǒng)(MRI)磁場強度建標
    Convection: a neglected pathway for downward transfer of wind energy in the oceanic mixed layer*
    藝術百家 王偉
    電影文學(2017年10期)2017-12-27 00:47:46
    一種永磁種子磁化機的設計
    超高商業(yè)大廈內部磁場強度的研究與分析
    科技資訊(2017年3期)2017-03-25 08:54:13
    雙星計劃中的雙星
    太空探索(2015年12期)2015-07-18 10:59:19
    成人欧美大片| 亚洲欧美日韩东京热| 伦精品一区二区三区| 伦理电影大哥的女人| 天天一区二区日本电影三级| av.在线天堂| 大型黄色视频在线免费观看| 国产真实伦视频高清在线观看 | 久久久久久久久久成人| 精品不卡国产一区二区三区| 亚洲精品一区av在线观看| 日本色播在线视频| 免费高清视频大片| 久久久久久大精品| 欧美在线一区亚洲| 麻豆国产av国片精品| 亚洲欧美激情综合另类| 美女 人体艺术 gogo| 在线看三级毛片| 精品午夜福利视频在线观看一区| 国产探花在线观看一区二区| 搞女人的毛片| 国产一区二区三区视频了| 此物有八面人人有两片| 热99在线观看视频| 91久久精品国产一区二区三区| 嫩草影院新地址| 一个人观看的视频www高清免费观看| 99久久精品国产国产毛片| 在线免费观看不下载黄p国产 | 最好的美女福利视频网| 嫩草影视91久久| 免费av观看视频| 精品国产三级普通话版| 九九在线视频观看精品| 成人一区二区视频在线观看| 极品教师在线视频| 在现免费观看毛片| 欧美性猛交╳xxx乱大交人| 舔av片在线| 亚洲精品色激情综合| 搡老岳熟女国产| 少妇裸体淫交视频免费看高清| 国产精华一区二区三区| 欧美最黄视频在线播放免费| 亚洲精品亚洲一区二区| 有码 亚洲区| 1024手机看黄色片| 中文字幕高清在线视频| 村上凉子中文字幕在线| 亚洲aⅴ乱码一区二区在线播放| 两人在一起打扑克的视频| 18+在线观看网站| 12—13女人毛片做爰片一| 国产av不卡久久| 春色校园在线视频观看| 毛片女人毛片| 久9热在线精品视频| 观看免费一级毛片| 夜夜看夜夜爽夜夜摸| 精品一区二区免费观看| 日韩亚洲欧美综合| 九色国产91popny在线| 在线观看66精品国产| avwww免费| 免费搜索国产男女视频| 亚洲精品乱码久久久v下载方式| 天天一区二区日本电影三级| 日韩欧美一区二区三区在线观看| 97热精品久久久久久| 日韩欧美在线二视频| 亚洲色图av天堂| 少妇丰满av| 69av精品久久久久久| 91精品国产九色| 一个人看视频在线观看www免费| 久久人人精品亚洲av| 最近视频中文字幕2019在线8| 18禁黄网站禁片午夜丰满| 国产探花在线观看一区二区| 国产av在哪里看| 日本免费a在线| 啦啦啦啦在线视频资源| 看十八女毛片水多多多| 美女高潮喷水抽搐中文字幕| 免费观看的影片在线观看| 亚洲精品456在线播放app | 久久精品国产亚洲av香蕉五月| 99国产极品粉嫩在线观看| 欧美日本亚洲视频在线播放| 国产私拍福利视频在线观看| 我要看日韩黄色一级片| 一夜夜www| 成人无遮挡网站| 最近在线观看免费完整版| 无人区码免费观看不卡| 精品99又大又爽又粗少妇毛片 | 国产蜜桃级精品一区二区三区| 琪琪午夜伦伦电影理论片6080| 琪琪午夜伦伦电影理论片6080| 69av精品久久久久久| 高清日韩中文字幕在线| av在线观看视频网站免费| 欧美日韩黄片免| 丰满的人妻完整版| 日本爱情动作片www.在线观看 | 国产精品三级大全| 校园人妻丝袜中文字幕| а√天堂www在线а√下载| 天美传媒精品一区二区| 日韩欧美精品免费久久| av女优亚洲男人天堂| 蜜桃亚洲精品一区二区三区| 亚洲一区高清亚洲精品| 啦啦啦韩国在线观看视频| 国内精品美女久久久久久| 男人舔女人下体高潮全视频| 精品午夜福利在线看| 有码 亚洲区| 联通29元200g的流量卡| 精品99又大又爽又粗少妇毛片 | 最新中文字幕久久久久| 国产亚洲91精品色在线| 琪琪午夜伦伦电影理论片6080| 69人妻影院| 国内毛片毛片毛片毛片毛片| a在线观看视频网站| 午夜福利视频1000在线观看| 干丝袜人妻中文字幕| 国内久久婷婷六月综合欲色啪| 黄色欧美视频在线观看| 久久久久久久亚洲中文字幕| 日韩在线高清观看一区二区三区 | 国产精品av视频在线免费观看| 久久精品国产亚洲av香蕉五月| 乱人视频在线观看| 直男gayav资源| 国产极品精品免费视频能看的| 少妇的逼好多水| 美女 人体艺术 gogo| 国产伦一二天堂av在线观看| 亚洲一级一片aⅴ在线观看| 麻豆成人午夜福利视频| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品国产三级国产av玫瑰| 麻豆国产av国片精品| 日韩在线高清观看一区二区三区 | 99热6这里只有精品| 午夜日韩欧美国产| avwww免费| 日本欧美国产在线视频| 成人永久免费在线观看视频| 他把我摸到了高潮在线观看| 国产精品美女特级片免费视频播放器| 国产免费男女视频| 国产精品,欧美在线| 国产人妻一区二区三区在| 亚洲无线观看免费| 国产伦精品一区二区三区视频9| 一区二区三区免费毛片| 亚洲国产日韩欧美精品在线观看| 午夜a级毛片| 久久久久国内视频| 国产午夜精品论理片| 男人舔女人下体高潮全视频| 一本一本综合久久| 中国美白少妇内射xxxbb| eeuss影院久久| 国产亚洲精品久久久com| 欧美最黄视频在线播放免费| 精品人妻视频免费看| 色播亚洲综合网| 99热这里只有是精品50| 亚洲av中文av极速乱 | 中文亚洲av片在线观看爽| 国产欧美日韩精品亚洲av| 桃色一区二区三区在线观看| 亚洲人成网站在线播放欧美日韩| 舔av片在线| 别揉我奶头 嗯啊视频| 亚洲一区二区三区色噜噜| 久久久久久大精品| 国产久久久一区二区三区| 动漫黄色视频在线观看| 深爱激情五月婷婷| 国产精品一区二区免费欧美| 日本色播在线视频| 午夜爱爱视频在线播放| 国产精品亚洲一级av第二区| 国产精品三级大全| 亚洲男人的天堂狠狠| 色哟哟哟哟哟哟| 极品教师在线免费播放| 老女人水多毛片| 久久精品人妻少妇| videossex国产| 成人一区二区视频在线观看| 国产不卡一卡二| 亚洲最大成人手机在线| 天天一区二区日本电影三级| 综合色av麻豆| 国产在线精品亚洲第一网站| 久久精品久久久久久噜噜老黄 | 精华霜和精华液先用哪个| 国产成人影院久久av| 天堂av国产一区二区熟女人妻| 国产 一区精品| 久久午夜亚洲精品久久| 色综合色国产| 18禁在线播放成人免费| 啦啦啦韩国在线观看视频| 亚洲欧美日韩无卡精品| av女优亚洲男人天堂| 国产精品一区二区三区四区免费观看 | 午夜a级毛片| 别揉我奶头~嗯~啊~动态视频| 日韩国内少妇激情av| 日本精品一区二区三区蜜桃| 国产69精品久久久久777片| 亚洲色图av天堂| 国产一区二区亚洲精品在线观看| 久久久久久久久久久丰满 | 久久国内精品自在自线图片| 天天一区二区日本电影三级| 亚洲无线在线观看| 十八禁国产超污无遮挡网站| 最新中文字幕久久久久| 97超级碰碰碰精品色视频在线观看| 成年版毛片免费区| 国产精品精品国产色婷婷| 国产高清三级在线| 性欧美人与动物交配| 国产精品久久久久久久电影| 性插视频无遮挡在线免费观看| 欧美潮喷喷水| 久久久久久伊人网av| aaaaa片日本免费| 在线a可以看的网站| 中文字幕熟女人妻在线| 联通29元200g的流量卡| 免费不卡的大黄色大毛片视频在线观看 | 成年女人毛片免费观看观看9| 日韩中文字幕欧美一区二区| 亚洲国产精品成人综合色| 啦啦啦观看免费观看视频高清| 五月伊人婷婷丁香| 热99re8久久精品国产| 国产黄片美女视频| 国产精品女同一区二区软件 | 色综合站精品国产| 国产成人a区在线观看| 国产不卡一卡二| 狂野欧美白嫩少妇大欣赏| 三级男女做爰猛烈吃奶摸视频| 国产精品爽爽va在线观看网站| 直男gayav资源| 国产男靠女视频免费网站| 最后的刺客免费高清国语| 两性午夜刺激爽爽歪歪视频在线观看| 露出奶头的视频| 精品久久久久久久久久久久久| 亚洲欧美日韩无卡精品| 很黄的视频免费| 99视频精品全部免费 在线| 校园人妻丝袜中文字幕| 欧美高清成人免费视频www| 久久中文看片网| 亚洲精品456在线播放app | 91麻豆av在线| 日本免费一区二区三区高清不卡| 久久九九热精品免费| 看十八女毛片水多多多| 最近中文字幕高清免费大全6 | 88av欧美| 国内久久婷婷六月综合欲色啪| 日韩欧美精品v在线| 1024手机看黄色片| 亚洲国产高清在线一区二区三| 69av精品久久久久久| 亚洲国产日韩欧美精品在线观看| 露出奶头的视频| 国产色婷婷99| 国产精品一区二区三区四区久久| 久久久久久久久久久丰满 | 成人av在线播放网站| 黄色配什么色好看| 亚洲性夜色夜夜综合| 亚洲欧美日韩高清在线视频| 国产精品久久久久久精品电影| 99久国产av精品| ponron亚洲| 亚洲精品久久国产高清桃花| 黄色日韩在线| 欧美zozozo另类| 搡老熟女国产l中国老女人| 久久久久久久久久黄片| 黄色丝袜av网址大全| 美女免费视频网站| 久久久国产成人精品二区| 成人特级黄色片久久久久久久| 啪啪无遮挡十八禁网站| 亚洲欧美日韩东京热| 我要看日韩黄色一级片| 久久久久久久久中文| 日韩精品有码人妻一区| 欧美极品一区二区三区四区| 国产三级在线视频| 色视频www国产| 在线观看66精品国产| 亚洲av第一区精品v没综合| 99久久精品一区二区三区| 他把我摸到了高潮在线观看| 欧美国产日韩亚洲一区| 男女做爰动态图高潮gif福利片| xxxwww97欧美| 99国产极品粉嫩在线观看| 欧美色视频一区免费| 国产精品1区2区在线观看.| 一区二区三区四区激情视频 | 啦啦啦韩国在线观看视频| 免费av毛片视频| 91精品国产九色| 午夜福利18| 美女 人体艺术 gogo| 人妻夜夜爽99麻豆av| 日韩欧美精品免费久久| 日本成人三级电影网站| 97热精品久久久久久| 网址你懂的国产日韩在线| 少妇裸体淫交视频免费看高清| 小蜜桃在线观看免费完整版高清| 哪里可以看免费的av片| 在线播放国产精品三级| 精品一区二区免费观看| 精品99又大又爽又粗少妇毛片 | 成人美女网站在线观看视频| 亚洲av免费在线观看| 别揉我奶头 嗯啊视频| 欧美日韩精品成人综合77777| 女人被狂操c到高潮| 国产精品久久久久久亚洲av鲁大| 51国产日韩欧美| 永久网站在线| 1024手机看黄色片| 22中文网久久字幕| av女优亚洲男人天堂| 三级国产精品欧美在线观看| 国产精品嫩草影院av在线观看 | 欧美成人性av电影在线观看| 一个人看视频在线观看www免费| 欧美一区二区亚洲| 日日摸夜夜添夜夜添av毛片 | 少妇高潮的动态图| 一个人观看的视频www高清免费观看| 亚洲国产精品合色在线| 亚洲国产精品成人综合色| 国产精品国产高清国产av| 亚洲av第一区精品v没综合| 尾随美女入室| 亚洲av第一区精品v没综合| 国产精品亚洲美女久久久| 18禁在线播放成人免费| 欧美激情国产日韩精品一区| 嫩草影院新地址| 两性午夜刺激爽爽歪歪视频在线观看| 欧美一区二区亚洲| 国产探花极品一区二区| 亚洲欧美日韩卡通动漫| 亚洲色图av天堂| 国产精品久久久久久av不卡| 欧美最新免费一区二区三区| 午夜老司机福利剧场| 亚洲专区国产一区二区| 久久久久九九精品影院| 我的老师免费观看完整版| 久久午夜亚洲精品久久| 精品人妻偷拍中文字幕| 一进一出抽搐动态| 午夜影院日韩av| 日本黄色片子视频| 无遮挡黄片免费观看| 欧美中文日本在线观看视频| 美女大奶头视频| 干丝袜人妻中文字幕| 色在线成人网| 免费看日本二区| eeuss影院久久| 亚洲av不卡在线观看| 少妇人妻一区二区三区视频| 最新在线观看一区二区三区| 国产不卡一卡二| 欧美成人免费av一区二区三区| 男女下面进入的视频免费午夜| 欧美成人免费av一区二区三区| 97人妻精品一区二区三区麻豆| 亚洲欧美精品综合久久99| 国产 一区精品| 亚洲av电影不卡..在线观看| 最新中文字幕久久久久| 国产综合懂色| 久久人人爽人人爽人人片va| 免费一级毛片在线播放高清视频| 中文字幕熟女人妻在线| 成人高潮视频无遮挡免费网站| 国产精品国产三级国产av玫瑰| 亚洲av电影不卡..在线观看| 久久99热这里只有精品18| 中文字幕人妻熟人妻熟丝袜美| 国内揄拍国产精品人妻在线| 真实男女啪啪啪动态图| 亚洲av.av天堂| 亚洲在线观看片| 亚州av有码| 大又大粗又爽又黄少妇毛片口| 床上黄色一级片| 99精品在免费线老司机午夜| 俄罗斯特黄特色一大片| 日本熟妇午夜| 欧美xxxx黑人xx丫x性爽| 国产毛片a区久久久久| 天美传媒精品一区二区| 国产乱人伦免费视频| 久久久久国内视频| 免费观看在线日韩| 在线观看午夜福利视频| 久久久久国产精品人妻aⅴ院| av天堂在线播放| 亚洲成人久久爱视频| 成人高潮视频无遮挡免费网站| 窝窝影院91人妻| 村上凉子中文字幕在线| 国产淫片久久久久久久久| 欧美不卡视频在线免费观看| 波多野结衣高清作品| 淫妇啪啪啪对白视频| 成人国产麻豆网| 给我免费播放毛片高清在线观看| 国产精品99久久久久久久久| 久久精品影院6| 亚洲国产精品久久男人天堂| 日日干狠狠操夜夜爽| 91麻豆av在线| 嫩草影视91久久| 成年人黄色毛片网站| 精品不卡国产一区二区三区| 日日干狠狠操夜夜爽| 亚洲成av人片在线播放无| 国产亚洲91精品色在线| 日韩欧美国产在线观看| 草草在线视频免费看| 一级av片app| 久久精品夜夜夜夜夜久久蜜豆| 久久精品国产亚洲av天美| 国产精品人妻久久久影院| 国内精品久久久久精免费| 中文字幕熟女人妻在线| 一个人看视频在线观看www免费| 亚洲aⅴ乱码一区二区在线播放| 欧美区成人在线视频| 国产亚洲精品久久久com| 欧美+亚洲+日韩+国产| 欧美丝袜亚洲另类 | 一级黄片播放器| 日本与韩国留学比较| 内地一区二区视频在线| 亚洲18禁久久av| 波多野结衣巨乳人妻| 欧美日韩乱码在线| 在线a可以看的网站| 别揉我奶头~嗯~啊~动态视频| 22中文网久久字幕| 精品日产1卡2卡| 男人舔奶头视频| 久久久国产成人免费| 午夜影院日韩av| 免费黄网站久久成人精品| 我的女老师完整版在线观看| 啦啦啦观看免费观看视频高清| 久久久精品欧美日韩精品| 久久久久久久久久成人| 欧美zozozo另类| 欧美区成人在线视频| 国产美女午夜福利| 亚洲精品日韩av片在线观看| 日本在线视频免费播放| 麻豆国产97在线/欧美| 国产 一区 欧美 日韩| 精品久久久久久久久av| 久久久久免费精品人妻一区二区| 岛国在线免费视频观看| 亚洲色图av天堂| 亚洲av一区综合| 女同久久另类99精品国产91| 一进一出抽搐动态| 日韩中字成人| 九色国产91popny在线| 欧美性猛交黑人性爽| 美女被艹到高潮喷水动态| 极品教师在线免费播放| 欧美最新免费一区二区三区| 亚洲人与动物交配视频| 国产不卡一卡二| 国内精品久久久久精免费| 国产精品精品国产色婷婷| 精品福利观看| 亚洲人与动物交配视频| 一a级毛片在线观看| 亚洲三级黄色毛片| 一个人免费在线观看电影| 婷婷精品国产亚洲av| 能在线免费观看的黄片| 3wmmmm亚洲av在线观看| 欧美最黄视频在线播放免费| 亚洲av第一区精品v没综合| 国产午夜福利久久久久久| 无人区码免费观看不卡| 黄色女人牲交| 亚洲三级黄色毛片| 成人午夜高清在线视频| av在线蜜桃| 中出人妻视频一区二区| 欧美黑人欧美精品刺激| 亚洲性久久影院| 久久久久国内视频| 亚洲真实伦在线观看| 久久草成人影院| 欧美日韩精品成人综合77777| 黄色欧美视频在线观看| 一个人免费在线观看电影| 成年女人看的毛片在线观看| 免费在线观看成人毛片| 国产免费av片在线观看野外av| 亚洲在线自拍视频| 久久久久久久午夜电影| 黄色一级大片看看| 久久午夜亚洲精品久久| 免费看a级黄色片| 亚洲精品乱码久久久v下载方式| 美女高潮喷水抽搐中文字幕| 久久久久久久久久成人| 男插女下体视频免费在线播放| 亚洲成人久久性| 国产v大片淫在线免费观看| 精品午夜福利在线看| 免费看光身美女| 久久久色成人| 色噜噜av男人的天堂激情| 大又大粗又爽又黄少妇毛片口| 亚洲av熟女| 成人av一区二区三区在线看| 内地一区二区视频在线| 国产一区二区激情短视频| 一个人看的www免费观看视频| 精品人妻熟女av久视频| www.色视频.com| 午夜日韩欧美国产| 男人的好看免费观看在线视频| 免费在线观看影片大全网站| 久久久色成人| 日韩中字成人| 别揉我奶头 嗯啊视频| 久久久精品欧美日韩精品| 国产高清不卡午夜福利| 色视频www国产| 97超视频在线观看视频| 一个人看视频在线观看www免费| 国产成人影院久久av| 国产男靠女视频免费网站| 国产aⅴ精品一区二区三区波| 淫秽高清视频在线观看| 夜夜看夜夜爽夜夜摸| 久久久国产成人精品二区| 长腿黑丝高跟| 欧美日韩瑟瑟在线播放| 无人区码免费观看不卡| 久久久久久伊人网av| 伦精品一区二区三区| 亚洲av免费高清在线观看| 色噜噜av男人的天堂激情| 国产激情偷乱视频一区二区| 日韩欧美精品免费久久| 2021天堂中文幕一二区在线观| 国产色爽女视频免费观看| 日日干狠狠操夜夜爽| ponron亚洲| 国产在视频线在精品| 乱人视频在线观看| 国产麻豆成人av免费视频| 久久99热这里只有精品18| 伦精品一区二区三区| 最近在线观看免费完整版| 男女做爰动态图高潮gif福利片| 欧美高清成人免费视频www| 给我免费播放毛片高清在线观看| 亚洲成人中文字幕在线播放| 午夜视频国产福利| 一区二区三区免费毛片| 亚洲国产精品sss在线观看| 联通29元200g的流量卡| 欧美成人a在线观看| 午夜影院日韩av| 亚洲三级黄色毛片| 少妇裸体淫交视频免费看高清| 国产老妇女一区| 很黄的视频免费| 国产不卡一卡二| 国产综合懂色| 男女边吃奶边做爰视频| 亚洲性久久影院| 久9热在线精品视频| www.色视频.com| 欧美激情国产日韩精品一区| 熟女人妻精品中文字幕| 免费av毛片视频| 看免费成人av毛片|