215 Acta Polytechnica CTU Proceedings 1(1): 215–221, 2014 215 doi: 10.14311/APP.2014.01.0215 Highly Magnetized Accreting Pulsars: Are There Accreting Magnetars? Pere Blay1, Pablo Reig2, Vı́ctor Reglero1 1IPL, University of Valencia, Spain 2University of Crete, Greece Corresponding author: pere.blay@uv.es Abstract 2S 0114+650, GX 301-2, IGR J16358-4726, X Per, 4U 2206+54, SXP 1062, and 3A 1954+319 are thought to possess high magnetic fields. They have recently been named accreting magnetars, or highly magnetized accreting pulsars. In this work their properties are reviewed. Within the context of their observational properties (mainly from INTEGRAL data), and the recent models of accretion onto highly magnetized neutron stars, their similarities and differences are analyzed. The aim is to find a common framework to understand the evolution (in terms of past and present history) of these sources, and to establish the basis of a possible new kind of accreting sources. Two of these sources, namely X Per and 4U 2206+54, contain a massive main-sequence companion, while the rest are supergiant X-ray binaries or symbiotic systems. The variety of astrophysical scenarios represented by this set is wide, therefore the study of these systems is also important in order to establish commonalities between the different types of accreting X-ray pulsars and to study the possible evolutionary links between them. Keywords: stars: binaries: close - stars: magnetars - stars: neutron - X-rays: binaries. 1 Introduction Magnetars are neutron stars with bright persistent X- ray emission (on the order of LX ∼ 1033−36ergs−1) and spin periods on the order of a few seconds. They show a very drastic spin period evolution, with large spin period derivatives and a long-term spin down ten- dency. The fact that they are isolated objects, together with the difficulties in considering their rotational en- ergy as the source of the observed bright X-ray emis- sion, are key issues in order to postulate their highly- magnetized nature (with extreme magnetic fields, on the order or greater than ∼ 1014Gauss). Magnetars are usually grouped into two subclasses, namely Anomalous X-Ray Pulsars (AXPs) and Soft Gamma-Ray Repeaters (SGRs). For a review of magnetars see, for example, Panchapakesan (2003) and Rea & Esposito (2011). High Mass X-ray Binary Systems (HMXRBS) con- tain a massive companion (≥10M�) and a neutron star or a black hole. The evolution of neutron stars in HMXRBS is affected directly by the wind plasma in which they are immersed (originating from the massive companion). The ambient plasma interacts with the magnetic field of the neutron star modifying its spin period evolution. Neutron stars are born as rapid ro- tators. The ambient plasma removes rotational energy from the neutron star and transports it outwards and the neutron star suffers an spin-down process, which goes on until the ambient matter becomes radiative, cools rapidly, and accretion can take place. The time spent in this initial spin-down epoch will determine the final X-Ray pulsar rotational period which will be an equilibrium period which resulting from the balance of two forces: the magnetic field pressure (which stops and drags the in-falling matter), and the inward pressure ex- erted by the ambient plasma. Initial modelings (see, for example, Illarionov and Sunyaev, 1975) predicted that the slowest reachable pulse period of a neutron star in a binary system is on the order of ∼500 s. Revisions of these initial models (see, for example, Davies and Pringle, 1981) relaxed this limit, allowing longer spin periods but only associated with higher neutron star magnetic fields. The possible presence of magnetars in some accret- ing HMXRBS, in which a slowly rotating neutron star has been found, has been proposed recently by a num- ber of authors (see, for example, Li & van den Heuvel 1999, Finger et al. 2010, Reig et al. 2012, Ikhsanov & Beskrovnaya 2010, Doroshenko et al. 2010, Patel et al. 2007, Popov & Turolla 2012, Bozzo et al. 2008). Our goal is to compare the observational properties of these systems and determine which observational parameters are more suitable to derive the strength of the magnetic field. We also want to stress the difficulties in deter- mining some of the properties of these systems when 215 http://dx.doi.org/10.14311/APP.2014.01.0215 Pere Blay, Pablo Reig, Vı́ctor Reglero considering only the long-term analysis (i.e., compar- ing pulse periods from different epochs), or the tempta- tion to identify spectral features as Cyclotron Resonant Scattering Features (CRSF) and derive magnetic field values from their central energy. Despite the small number of systems, the variety of astrophysical scenarios represented by this set is wide, from main sequence accreting binary systems (X Per, SXP 1062 and 4U 2206+54) to supergiants (2S 0114+65 and GX 301-2), and even symbiotic systems (IGR J16358-4726 and 3A 1954+319). In this sense, the study of these systems is important also to establish commonalities between the different types of HMXRBs and to study the possible evolutionary links between them. In this work we will review the observational properties of some of these systems and we will try to establish the scientific rational to develop a deeper analysis of their common properties and evolutionary status. Due to the large inhomogeneity in the energy ranges and sets of data available in the literature, we have started an analysis of these systems by using IN- TEGRAL/IBIS/ISGRI1 data, setting up common en- ergy ranges, spectral resolutions, and performing a ho- mogeneous timing analysis. 2 The Candidates A number of HMXRBs have already been proposed to contain highly magnetized neutron stars, namely 2S 0114+65 (Li & van den Heuvel 1999), 4U 2206+54 (Reig et al. 2012, Finger et al. 2010), X Per (Doroshenko et al. 2012), GX 301-2 (Ikhsanov & Finger 2012), IGR J16358-476 (Patel et al. 2007). There are also sugges- tions that SXP 1062 (Fu & Li 2012) and 3A 1954+319 (Marcu et al. 2011) contain a highly magnetized neu- tron star. Attempts have also been made to try to explain the properties of SFXTs with the presence of a highly-magnetized neutron star with some degree of success (see Bozzo et al. 2008 ) . A complete list of sources has been compiled in Table 1. Despite the fact that there are peculiarities in each of these systems which make difficult to group them as a kind (see, for example, Reig et al. 2012, Blay and Reglero 2013), we can summarize their overall commonalities: • slowly rotating neutron stars (Pspin slower than 500 s), with • high magnetic fields (larger than 1013 Gauss) • most of them show a long-term spin-down ten- dency The neutron star properties in these systems are differ- ent from those of magnetars, for this reason the terms magnetar-like, or magnetar-descendant (making refer- ence to a possible evolution from a neutron star born as a magnetar), have also been used to identify these sources (Reig et al. 2012, Blay and Reglero 2013). Table 1: List of sources. Source Pspin (s) Porb (d) Type 4U 2206+54 5560 19.2 MS donor X Per 864 250 BeXRB 2S 0114+650 9700 11.6 SXRB GX 301-2 680 41.498 SXRB IGR J16358-4726 5880 −− Symbiotic 3A 1954+319 19080 −− Symbiotic SXP 1062 1062 ∼300 BeXRB After 10 years of operations of the INTEGRAL satellite, a very complete and homogeneous long-term database of observations of all these sources has been compiled. We have analysed all public observations in which 4U 2206+54, X Per, 2S 0114+65, IGR J16358- 476, and GX 301-2 were in the field of view of INTE- GRAL/IBIS/ISGRI. The behavior of these source is similar in long time scales, with transient emission in the form of bursts or peaks and more or less continuous detection. They all show the typical behavior of wind accreting systems (al- though X Per is a BeX, it has been proposed that its behavior resembles that of wind-fed systems because of the large orbital distance, which will keep the neutron star away from the Be disk, see Doroshenko et al. 2012). The light curves, in every case, are complete enough to develop both short and long-term analysis of the spec- tral and timing properties of these systems. The magnetic properties of the neutron stars in these systems have been investigated/explained mainly in these terms: • spin period (neutron star evolutionary considera- tions) • spectral features (Cyclotron Resonant Scattering Features, CRSF) 1INTEGRAL (INTERnational Gamma-RAy Laboratory) is an ESA mission with contributions from USA and Russia, with two high-energy imagers (IBIS, working in the energy range 15 kev−10 MeV , and JEM-X sensitive in the 3−35 keV energy range) and one spectrometer (SPI, sensitive in the range 20 keV−8 MeV). IBIS, in turn, has two detector layers, ISGRI, for the lower energy band,15 keV−1 MeV, and PiCSIT working in the 175 keV−20.4 MeV energy range. A detailed description of the mission can be found in Winkler et al. (2003) 216 Highly Magnetized Accreting Pulsars: Are There Accreting Magnetars? With regard to the spin-period history, according to Ikhsanov (2007), the maximum spin period reachable by a neutron star can be related to the intensity of its magnetic field and the quantity of accreted matter (i.e., mass accretion rate). The latter quantity is considered to be directly proportional to the X-Ray luminosity of the source (see, for example, Finger et al. 2010). The relationship between these quantities is: Pmax = 15000µ 16/21 32 M −4/21 NS ( Ṁ 1015gs−1 )−5/7s (1) where µ32 is the dipole magnetic moment (in units of 1032Gcm−3) and MNS the mass of the neutron star (in units 1.5M�), and Ṁ is the mass accretion rate. A more recent calculation of quasi-spherical accretion by Popov and Turolla 2012, yields this relationship, directly in terms of the magnetic field intensity: B12 ∼ 8.1Ṁ 1/3 16 v −11/3 300 ( P1000 Porb,300 )11/2 (2) where B12 denotes the magnetic field in units of 1012 Gauss, Ṁ16 is the mass accretion rate in units of 1016gs−1, v300 is the wind velocity in units of 300kms−1, P1000 is the spin period in units of 1000s, and Porb,300 is the orbital period of the system in units of 300d. On the other hand, a direct measurement of the magnetic field of the neutron star is possible when its high-energy spectrum shows CRSFs. CRSFs are due to the splitting of energy levels of the electron in the presence of a magnetic field. Therefore the scatter- ing of X-ray photons by these electrons is produced at quantized Landau levels producing the absorption-like features seen in the X-Ray spectra of many HMXRBs (see Schönherr et al. 2007). The magnetic field strength can be calculated from the position of the CRSF by the formula: ECRSF = 11.6 B12 (1 + z) (3) where z is the gravitational redshift. It should be noted that this redshift, and consequently the measured mag- netic field, will depend on the height over the neutron star surface where the CRSF is formed. In the next sections, we will review these methods when applied mainly to some of sources listed in Table 1. 3 Magnetic Field Determination via Spin Period Evolution Fig. 1 shows examples of spin period histories of 4U 2206+54, 2S 0114+650, IGR 16358-4726, and X Per. 53600 53700 53800 53900 54000 54100 54200 54300 MJD 5300 5400 5500 5600 5700 5800 P sp in ( s) 4U 2206+54 53000 53500 54000 54500 55000 55500 56000 MJD 2.00 2.25 2.50 2.75 3.00 P sp in (h ) 2S 0114+650 54000 54200 54400 54600 54800 55000 55200 55400 MJD 830 831 832 833 834 835 836 837 838 839 840 841 842 843 P sp in ( s) X Per 53000 53500 54000 54500 55000 55500 56000 MJD 5780 5800 5820 5840 5860 5880 5900 P sp in ( s) IGR J16358-4726 Figure 1: Spin period evolution of 4U 2206+54 (top panel), 2S 0114+650 (second panel from the top), X Per (third panel from the top), and IGR 16358-4726 (lower panel) as obtained from INTEGRAL/IBIS/ISGRI data in the 20-40 keV energy range. 217 Pere Blay, Pablo Reig, Vı́ctor Reglero 0.0001 0.001 0.01 0.1 1 F lu x (c ou nt s -1 k eV -1 ) 20 40 60 80 100 120 140 160 180200 Energy (keV) -0.04 -0.02 0 0.02 0.04 re si du al s 4U 2206+54 0.01 0.1 1 F lu x (c ou nt s s- 1 ke V -1 ) 20 40 60 80 100 120 140 160 180200 Energy (keV) -0.4 -0.2 0 0.2 0.4 re si du al s 2S 0114+650 0.0001 0.001 0.01 0.1 1 F lu x ( co un t s- 1 ke V - 1 ) 20 30 40 50 60 70 80 90 100 200120 140 160 180 Energy (keV) -0.02 -0.01 0 0.01 0.02 re si du al s X Per Figure 2: The 20-200 keV spectra of 4U 2206+54 (top panel), 2S 0114+650 (middle panel) and X Per (bottom panel) extracted from INTEGRAL/IBIS/ISGRI data. The spectrum of 4U 2206+54 has been fitted with a bulk comptonization model with a χ2red ∼ 1.4, no evi- dence of a feature around ∼30 keV can be seen, neither in the spectrum nor in the residuals. For 2S 0114+650 a bulk comptonization model has been used to fit the data with a χ2red ∼ 1.2, again no evidence of a CRSF at ∼30 keV is evident. For the case of X Per, a comp- tonization model and two absorption features at 31.2 keV and 82.3 keV (indicated by arrows) have been used to fit the observed spectrum, with a χ2red ∼ 1.1. In the case of 4U 2206+54, frequency derivative val- ues derived by Finger et al. (2010) and Reig et al. (2012) are 1.5 × 10−14Hz s−1 and 1.7 × 10−14Hz s−1 respectively. Finger et al. 2010 included an older Bep- poSAX measurement. From these spin-down measure- ments, Finger et al. (2010) estimated a magnetic field of B∼ 1014 G. By using the relationship of Ikhsanov (2007) they also derived a magnetic field for 4U 2206+54 on the order of 1014 G. Reig et al. (2012) calculated a magnetic field strength on the same order of magni- tude by considering the relationship of Popov and Tur- olla (2012). The upper panel of Fig. 1 shows the spin period measurements from INTEGRAL/ISGRI data of 4U 2206+54 in the period 53600 − 54300 MJD. During this epoch (see Blay and Reglero 2013) the data is com- patible with no changes in the spin period of the system, as the pulse period measurements show large errors. The same thing happens for X Per and 2S 0114+650 (second and third panels from the top, respectively), in which large errors in the spin period determination make the measurements compatible with no pulse pe- riod change. Therefore we cannot report on magnetic field measurements from spin period changes for these three systems. The luminosity of 2S 0114+65 in the epoch shown in Fig. 1 is 4 × 1037ergs−1 (estimated from the model fit- ted to the spectrum shown in Fig. 2), considering a dis- tance to the source of 7.2 kpc (Hall et al. 2000). There- fore, the Ikhsanov (2007) relationship yields a magnetic moment on the order of 1033 G cm−3, which yields a magnetic field strength of B ∼ 1015 G. The relationship from Popov and Turolla (2012), however, results in an unreasonable number when applied to this source. Doroshenko et al. (2012) argue that, for the case of X Per, according to the torques applied to the neu- tron star magnetosphere in the case of wind accre- tion, and taking into account the theoretical approaches of Doroshenko et al. (2010), Davidson and Ostriker (1973), Davies et al. (1979), and Bisnovatyi-Kogan (1991), the magnetic field in X Per is expected to be on the order of B∼ 1014. We show in the lower panel of Fig. 1 the spin period measurements made with INTE- GRAL/ISGRI data. For this source We can estimate the luminosity to be Lx ∼ 1.09 × 1035ergs−1 in this epoch (estimated from the model fitted to the spectrum shown in Fig. 2 ). As the luminosity of the source is closely related to the accretion rate, we can estimate the field strength of X Per by using the relationship of Popov and Turolla (2012), yielding B∼ 1012 G, which falls within the range of normal neutron star magnetic fields. The relationship from Ikhsanov (2007) results in a magnetic moment of µ ∼ 1030 G cm−3, which implies a magnetic field on the order of B ∼ 1012 G. We find, therefore, a discrepancy with the results of Doroshenko et al. (2012). We cannot reproduce their calculations, however, because of the large errors in the spin period determination, which prevent the determination of a reliable frequency derivative, needed for the determina- tion of the spin torques. For IGR J16358-4726 we see in the lower panel of 218 Highly Magnetized Accreting Pulsars: Are There Accreting Magnetars? Fig. 1 that the source experiences large pulse period changes. It shows a spin-down between MJD∼53250 and MJD∼54000 of ν̇ ∼ 2.4 × 10−14 Hz s−1. According to the Equation (9) from Finger et al. (2010), which relates the magnetic moment of the neutron star with the measured spin-down, we find in IGR J16358-4726 a minimum magnetic moment of µ ∼ 2 × 1032 G cm−3, which implies a magnetic field on the order of or larger than B∼ 1014 G. These calculations rely on a good measurement of the luminosity (or mass accretion rate) or of the changes in the spin period. The major uncertainty comes from the determination of the luminosity of the source, cal- culated from the source flux in a given energy range and the estimated distance to the source, which is usu- ally subject to large uncertainties. With respect to the use of spin period changes, the estimates always rely on long-term or average behavior. On the one hand this is due to the poor knowledge of the short-term or detailed pulse period evolution in many cases, and, on the other hand, to the lack of a detailed and complete theoretical description of the wind accretion scenario. The efforts of Postnov et al. (2013), however, try to fill this lack of good theoretical approximations. We will comment on that in the last section. 4 Magnetic Field Determination via Cyclotron Resonant Scatering Features A more direct measurement of the magnetic field strength of the neutron star in these systems can be obtained by determining the position of the fundamen- tal CRSF seen in their X-Ray spectra. CRSF have been reported in GX 301-02 (La Barbera et al. 2005, Makishima and Mihara 1992), 2S 0114+65 (Wang 2010, Bonning and Falanga 2005), 4U 2206+54 (Blay et al. 2005, , Wang 2009, Torrejón et al. 2004, Masseti et al. 2004), and X Per (this work, see Fig. 2, Doroshenko et al. 2010), see Table 2. Table 2 summarizes the CRSF reported in these sys- tems and the magnetic field strength measurement de- rived. Table 2: CRSFs and B determination. See references in the text. Source ECRSF (keV) B (10 12 Gauss) 4U 2206+54 30 2.6 X Per 31.2 2.7 2S 0114+650 22 1.9 GX 301-2 37 3.2 5 2S 0114+650 22(c) 1.9 We see in all cases that the magnetic field strength derived by this method is on the order of those found in other neutron stars in HMXRBs, i.e., 100 to 1000 times lower than those of the magnetars. We find, therefore, a discrepancy in magnetic field strength determination with respect to the methods shown in the previous sec- tion. One could think that the magnetic fields derived from spin period evolution are not properly calculated, in view of the more direct determination from CRSF measurements. However the location of CRSF is not free of ambiguities. In the case of 4U 2206+54 (see top panel of Fig. 2), for example, this feature has been ob- served only marginally, and only at one epoch (see Blay et al. 2005 and references there in). It has not been ob- served again (see Blay and Reglero 2013). Therefore, although it has been assumed that the absorption ob- served was a CRSF, this may not be the case. A sim- ilar situation applies to 2S 0114+65 (see middle panel of Fig. 2), in which the detection of the possible CRSF has only been marginal (see Wang 2010) and not con- firmed again. In the case of X Per, Doroshenko et al. (2010) show how instead of a CRSF fundamental and its first harmonic (see bottom panel of Fig. 2), the X Per spectra can be very well described by a combina- tion of two comptonization models (compTT+compTT, in XSPEC notation) demonstrating the possibility that the features observed may not be CRSFs. Therefore, at least in some cases, the possibility of finding magnetar-like magnetic fields in these systems is not ruled out by the reported observations of CRSFs. 5 Discussion and Conclusion We have seen that in most cases the determination of the magnetic field strength of the neutron star is am- biguous, or even contradictory, when the result from different approaches are compared. We want to em- phasize how different approaches may seem to result in different magnetic field strength estimates. An homo- geneous and coherent approach is needed in order to determine if there are links between these systems (in terms of their spin period and binary evolution) and to support or discard the different theories trying to explain the puzzling presence of long spin periods in HMXRBs. There have been recent theoretical efforts in order to provide better modeling of the wind-accretion sce- nario. Ikshanov and Beskrovnaya (2013), for example, report on a likely explanation of the large spin-down shown in 4U 2206+54 by taking into account a higher than usual magnetic field of the companion star and, consequently, that the accreting material can be mag- netized. The possibility of an optical companion with a large magnetic field had been already proposed by Blay and Reglero (2011). 219 Pere Blay, Pablo Reig, Vı́ctor Reglero The efforts of Postnov et al. (2013) in model- ing more accurately the wind accretion mechanism in HMXRBs also result in magnetic field estimations for 4U 2206+54 and SXP 1062 on the order of typical mag- netic fields found in accreting X-Ray pulsars (∼ 1012G). They explain the behavior of sources like GX 301-2 with the hypothesis of being older systems which have al- ready reached their equilibrium period. Together with an improvement in the modeling of wind accretion in HMXRBs, a better observational ap- proach is needed in order to understand the behavior of these sources. Multiwavelength campaigns (IR, optical, UV, X-rays) on these objects are needed (as simultane- ous as possible) in order to determine the various pa- rameters involved: obtain more accurate spectral clas- sification and orbital solutions, measure accurate mass loss rates and wind velocity laws, etc. Acknowledgement P. Blay acknowledges funding from the Spanish Ministerio de Economia y Competitividad through project AYA-2011-29936-C05. References [1] Bisnovatyi-Kogan, G. S. 1991, A&A, 245, 528 [2] Blay P., Ribó M., Negueruela I., Torrejn J. M., Reig P., Camero A., Mirabel, I.F., & Reglero V. 2005, A&A 438, 936 [3] Blay P. and Reglero V. 2011, BSRSL, 80, 634 [4] Blay P. and Reglero V. 2013, PoS(INTEGRAL 2012)012 [5] Bonning E. W., Falanga M. 2005, A&A, 436, L31 [6] Bozzo E., Falanga M., Stella L. 2008, ApJ, 683, 1031 doi:10.1086/589990 [7] Davidson, K., and Ostriker, J. P. 1973, ApJ, 179, 585 [8] Davies, R. E., Fabian, A. C., and Pringle, J. E. 1979, MNRAS, 186, 779 [9] Davies R.E., and Pringle J.E. 1981, MNRAS, 196, 209 [10] Doroshenko V., Santangelo A., Suleimanov V., Kreykenbohm I., Staubert R., Ferrigno C., Klochkov D. 2010, A&A 515, A10 [11] Doroshenko V., Santangelo A., Kreykenbohm I., and Doroshenko R. 2012, A&A 540, L1 [12] Finger M.H., Ikhsanov N.R., Wilson-Hodge C.A., & Patel S.K. 2010., ApJ, 709, 1249 doi:10.1088/0004-637X/709/2/1249 [13] Fu L., Li X.D. 2012, ApJ, 757,171 doi:10.1088/0004-637X/757/2/171 [14] Fürst F., Marcu D.M., Pottschmidt K., Grin- berg V., Wilms J., and Cadolle Bel M. 2011, PoS(INTEGRALL 2010)017 [15] Hall T. A., Finley J. P., Corbet R. H. D., Thomas R. C., 2000, ApJ, 536, 450 doi:10.1086/308924 [16] Harbel, F., Sturm, R., & Filipović, M. D. et al. 2012, A&A, 537, L1 [17] Ikhsanov, N. R. 2007, MNRAS, 375, 698 doi:10.1111/j.1365-2966.2006.11331.x [18] Ikhsanov N.R., Finger M.H. 2012, ApJ, 753,1 doi:10.1088/0004-637X/753/1/1 [19] Ikhsanov N.R., Beskrovnaya N.G. 2010, Astro- physics, 53-2, 237 [20] Ikhsanov N. R. and Beskrovnaya N. G.2013, ARep, 57, 287 [21] Illarionov A.F., and Sunyaev R.A. 1975, A&A, 39,185 [22] Koh D.T., Bildsten L., Chakrabarty D., et al. 1997, ApJ, 479, 933 [23] Li X. D., van den Heuvel E. P. J., 1999, ApJ, 513, L45 [24] La Barbera, A., Segreto, A., Santangelo, A., Kreykenbohm, I., & Orlandini, M. 2005, A&A, 438, 617 [25] Makishima K., Mihara T. 1992, in Magnetic Fields of Neutron Stars, ed. T. Tanaka & K. Koyama (Tokyo: Universal Academy Press), 23 [26] Marcu D. M., Fürst F., Pottschmidt K., et al. 2011, ApJ, 742, L11 doi:10.1088/2041-8205/742/1/L11 [27] Masetti N., Dal Fiume D., Amati L., Del Sordo S., Frontera F., Orlandini M., and Palazzi E. 2004, A&A, 423, 311 [28] Panchapakesan N. 2003, Bull. Astr. Soc. India, 31, 19 [29] Patel S.K., Zurita J., Del Santo M. et al. 2007, ApJ, 657, 994 doi:10.1086/510374 [30] Popov S.B., Turolla R. 2012, MNRAS, 421, L127 doi:10.1111/j.1745-3933.2012.01220.x 220 http://dx.doi.org/10.1086/589990 http://dx.doi.org/10.1088/0004-637X/709/2/1249 http://dx.doi.org/10.1088/0004-637X/757/2/171 http://dx.doi.org/10.1086/308924 http://dx.doi.org/10.1111/j.1365-2966.2006.11331.x http://dx.doi.org/10.1088/0004-637X/753/1/1 http://dx.doi.org/10.1088/2041-8205/742/1/L11 http://dx.doi.org/10.1086/510374 http://dx.doi.org/10.1111/j.1745-3933.2012.01220.x Highly Magnetized Accreting Pulsars: Are There Accreting Magnetars? [31] Postnov K.A., Shakura N.I., Kochetkova A.Yu., and Hjalmarsdotter L., 2013, arXiv:1307.3032 [32] Rea N., Esposito P. 2011, in High-Energy Emission from Pulsars and their Systems, ed. D. F. Torres & N. Rea, Astrophysics and Space Science Proceed- ings (Springer Berl in Heidelberg), 247 [33] Reig P., Torrejón J.M., & Blay P. 2012, MNRAS, 425, 595 doi:10.1111/j.1365-2966.2012.21509.x [34] Schönherr G., Wilms J., Kretschmar P., Kreyken- bohm I., Santangelo A., Rothschild R.E., Coburn W., Staubert R. 2007, A&A 472, 353 [35] Torrejn J. M., Kreykenbohm I., Orr A., Titarchuk L., and Negueruela I. 2004, A&A, 423, 301 [36] Wang, W. 2009, MNRAS, 398, 1428 doi:10.1111/j.1365-2966.2009.15200.x [37] Wang W. 2010, A&A 516, A15 [38] Winkler C., Courvoisier T.J.L., Di Cocco G., et al. 2003, A&A, 411, L1 221 http://dx.doi.org/10.1111/j.1365-2966.2012.21509.x http://dx.doi.org/10.1111/j.1365-2966.2009.15200.x Introduction The Candidates Magnetic Field Determination via Spin Period Evolution Magnetic Field Determination via Cyclotron Resonant Scatering Features Discussion and Conclusion