arXiv:astro-ph/0609676v1 25 Sep 2006 G 21+03i oae ntefil fteSRG78.2+2.1 SNR (SNRs). as starform- the of remnants in known field the supernova rich (commonly in located very and is is binaries, J2018+4043 region IGR X-ray Cygnus sites, The ing band. keV 40 bevtoswith Observations ul oe edo iw(CO)osrainyeddthe ks yielded 320 (1 observation The (FCFOV) of view flux 2004). of source al. field et coded (Bykov fully J2018+4043 IGR source, u6Aû 7 â 5,40 ig,Belgium Liège, 4000 B5c, Bât 17, Août 6 du s mn paetysed ndnie GE sources. EGRET unidentified bright- the steady is apparently 1999), al. among et Hartman est 1995; Dermer & (Sturner vteCdx France Cedex, Yvette ore C 7+,i the in 078+2, 2CG source, hs htnflxis J2020+4023), GEV flux photon = whose J2020+4017 3EG = J2020+4026 (2EG rila 0,N-08S mtra,TeNetherlands The Amsterdam, SJ NL-1098 403, Kruislaan o h u bv 0 e,idctn ra ntespec- the the the that be in suggested could (2004) sources break ISGRI al. and a et EGRET Bykov GeV. indicating few GeV, a 300 above above trum flux the for AUrct h Netherlands The Utrecht, CA ihthe with 98 aeetbihda pe ii f2.2 of limit upper an established have 1998) tory 48Beo ie,Argentina Aires, Buenos 1428 Whipple oneprso h aesuc in source same the of counterparts A16802 PA 901 [email protected] 194021; einwt h BSIGIiae Lbu ta.2003) al. et (Lebrun imager IBIS-ISGRI the aboard the with region .B H. 8 7 6 5 4 3 2 rpittpstuigL using typeset Preprint 1 h rthr -a mgn bevto fteCygnus the of observation imaging X-ray hard first The T ..B A.M. ntttdAtohsqee eGohsqe nvriéde Université Géophysique, de France et Paris, d’Astrophysique 75231 Institut M.Berthelot, Place 11 7164, APC-UMR E-aly S/ANASriedAtohsqe 19 G 91191 d’Astrophysique, DSM/DAPNIA/Service CEA-Saclay, srnmclIsiue“no anke” nvriyo A of University Pannekoek”, “Anton Institute Astronomical RNNtelnsIsiuefrSaeRsac,Sorbonnela Research, Space for Institute Netherlands SRON nttt eAtooí íiadlEpco(AE,C 7 S 67, CC (IAFE), Espacio del Física y Astronomía de Instituto enyvnaSaeUiest,55DvyLbrtr,Univ Laboratory, Davey 525 University, State Pennsylvania ..IfeIsiueo hsc n ehooy t Petersb St. Technology, and Physics of Institute Ioffe A.F. HE LOEMEN ( INTEGRAL A STROPHYSICAL aaw on lgtyetne n ihyasre object absorbed highly and extended slightly a found we data wt 4 a (with ks 5 A nosue crtn bet ..a G raGlci X-ray Galactic bright a the or to AGN an J2018+4043 IGR i.e. of f connection object, distance accreting J20 beam-width obscured IGR a of an at spectra multiwavelength peaked and source morphology observed GHz 1.4 VLA extended index photon with model law 3.4 of flux keV 0.5–10 snwcndnl eetdas nte4-0kVbn,wt fl a with band, keV 40-80 the in also sour detected The confidently field. now source is the of observations multi-wavelength ujc headings: Subject EGRET. γ O B COS YKOV rytlsoeosrain eg,Bclye al. et Buckley (e.g., observations telescope -ray efudavr ieycutratt h eetydiscovere recently the to counterpart likely very a found We nentoa am-a srpyisLabora- Astrophysics Gamma-Ray International 4 .H W. , 1 ..K A.M. , Swift . ike ta.20)hsrvae new a revealed has 2003) al. et Winkler ; 7 A . ′′ aelt eg,Saebr ta.1981). al. et Swanenburg (e.g., satellite T CGRO E netit) h combined The uncertainty). 2 ± ERMSEN tl mltajv 08/13/06 v. emulateapj style X J 1. ∼ OURNAL 0 bevto fteIRJ0844 edrvae adpoint hard a revealed field J2018+4043 IGR the of observation γ . 4) INTRODUCTION 1.2 CgiSR.Ahigh-energy A SNR). -Cygni RASSILCHTCHIKOV × S:idvda ( individual ISM: NTENTR FTEHR -A OREIRJ2018+4043 IGR SOURCE X-RAY HARD THE OF NATURE THE ON GE hwdta hssource this that showed EGRET 4,5 × 10 γ , .K J. , 10 V Cgifil a discovered was field -Cygni 4,L1L4 2006 L21-L24, 649, . − − + − 11 0 0 . . 6 8 7 γ × AASTRA cm rscm ergs Cgi,eg,aGeminga- a e.g., , -Cygni 10 − 2 − Γ γ 12 s ryadhr X-ray hard and -ray − 1 = 1 r cm erg 4 1 − h srpyia ora,v 4,L1L4 2006 L21-L24, 649, v. Journal, Astrophysical The .L F. , × Y , γ 2 bv 0 MeV 100 above . Cgi)Xry:idvda IRJ0844;3GJ2020+40 3EG J2018+4043; (IGR individual )—X-rays: -Cygni 3 10 s U − ± 1 A U .A. − EBRUN − 11 0 nte25- the in 2 ig,Allée Liège, r,Russia, urg, riyPark, ersity . Swift n2 3584 2, an cm and 2 s msterdam, − VAROV 1 c 28, uc. γ − γ 6 if-sur- 9%cndnelvl at level) confidence (90% 2 rysuc E 22+03dtce by detected J2020+4023 GEV source -ray .R M. , -ray - ABSTRACT s Integral − N 1 H 1 ..K J.A. , 6 = ENAUD . ido h al yeOVsa D132,lctdjs 7 just located 193322, HD with star J2018+4043, O9V IGR type from early the interactio of the wind particular, In accelerate ejecta/shell. particles supernova from the emission extended or pulsar like INTEGRAL the of sitioning i betNC63B(e asn ta.20;Bkve al. et Bykov 2005; al. et extragala Bassani obscured (see 2006). 6334B highly NGC a object is likely tic source most The that of 6334. the NGC identification onto region projected starforming J17204–3554 a Galactic IGR such rich source of neces- example the not An is does confusion connection. the field physical hand, SNR their imply other the sarily within the position On source ISGRI acceleration. particle of mechanism G,addsustepsiiiyta h GE oreis source EGRET the a that likely possibility its source, indeed the accreting discuss obscured and strongly AGN, a is source th suggest the we observations, § multiwavelength these and on (§2.3 Based data archival in counterparts radio/IR/optical nytesinewnosweetecutrt a within was taking rate criteria count selection the the where applied windows ± We science the 16). only Nov Dec (2002 2005 19–255 – revolutions 9 at taken data FCFOV archival a wthdofframr cuaeflxetmto.Broad estimation. flux spreading accurate Pixel more a the 2003). for with off al. switched et performed (Goldwurm was was software reduction 5.1 standard OSA The exposure. n h eeto facutrato G 21+03i the in (§2.1) J2018+4043 with IGR ISGRI band of 1.5 with keV counterpart recent a source 0.5–10 of of the detection results the of and the observations here FCFOV report (4–5 We Ms weak relatively instruments. the other of identification σ reliable a 2003), bu 1 about pcrmcnb ecie ya bobdpower- absorbed an by described be can spectrum 1 SR orerqie ihrslto bevtoswith observations high-resolution requires source ISGRI ) − + 1 ie h oia 12 nominal the Given h edof field The 2 3 ENNEA 6 e rgnlydsoee nte2-0kVband, keV 20-40 the in discovered originally ce, . . σ .T R. , 2 2 844 r ossetwt hs xetdfor expected those with consistent are 18+4043 × rudtema au eemndb h hl . Ms 1.5 whole the by determined value mean the around tthe at adXrysuc G 21+03i the in J2018+4043 IGR source X-ray hard d iay h dnicto ugsspossible suggests identification The binary. 10 ′ 2 ERRIER on orelclzto fIGI(erne al. et (Lebrun ISGRI of localization source point o h pia n -a oiin.The positions. X-ray and optical the rom 22 γ ..P G.G. , xo (1 of ux 2. rycutrat(§3). counterpart -ray Swift cm SR.W sd15M fIGIgo time good ISGRI of Ms 1.5 used We ISGRI. BEVTOSADDT ANALYSIS DATA OBSERVATIONS AND γ − α 7 Swift Cgihsbe rqetyosre with observed frequently been has -Cygni 2 .D M. , nacia pia n infrared and optical archival In . AVLOV 20 = orepsto.Teei loan also is There position. source 2.1. . 4 lk orewt h observed the with source -like ± h orealwdu ofidispossible its find to us allowed source SR observation ISGRI 18 E Swift 0 2 B ′ ..D G.M. , . m 4) ECKER nua eouin(WM and (FWHM) resolution angular 38 × γ R §.) h cuaepo- accurate The (§2.2). XRT . s 55, 10 Cgiwudb plausible a be would -Cygni 8 O B COS − .R G. , UBNER 11 δ = rscm ergs + AUW 40 3 and ..G E.B. , ◦ 17) 41 8 .P S J.-P. , CGRO − ′ 2 00 s . − ′′ IACANI 1 4 . WINGS fthe of n 2.4). 3 in d , 8 c- at n ′ 2 A.M.Bykov et al.

′ ′ FIG. 1.— 40-80 keV ISGRI significance map (≥ 3.5σ) of the γ-Cygni field, with VLA 1.4 GHz contours overlaid. The light and dark 5 × 5 ISGRI pixels correspond to 3.6 and 3.8σ significance, respectively. The VLA contours run from 9 to 73 mJy/beam on a linear scale. The Swift source position is marked by the white cross. The smoothed Swift XRT 2-8 keV image is shown in the inset. The bright O9V star HD193322 was removed from the XRT image. energy bands were used to increase S/N. The source was de- event files were processed with the xselect tool (grades 0– tected with a flux of (1.2±0.2) ×10−11 erg cm−2 s−1 in the 12 were used). The XSPEC package (ver. 12.1.1) was used 20-40 keV band, and (1.4±0.4)×10−11 erg cm−2 s−1 in the for the spectral analysis. A point-like source was detected in 40-80 keV band, obtained by cross calibration with the Crab the XRT image at a position consistent with that of the IS- (see, e.g., Lubinski et al. 2004). The flux in 20- 40 keV band GRI source (see Fig. 1): α=20h18m38.s55, δ=+40◦41′00.′′4 is below that quoted in Bykov et al. (2004) due to the new (J2000), with a position uncertainty of 4.′′2 (errors here and background model used in OSA 5 and more conservative se- below in this subsection are at the 90% confidence level). lection criteria that we applied now to the much longer ex- We extracted 105 source counts from an aperture of 47.′′2 posure dataset available. The 40-80 keV ISGRI significance radius and estimated the observed 0.5–10 keV flux FX = +0.7 × −12 −2 −1 map is presented in Figure 1. 3.4−0.8 10 erg cm s . We fit the spectrum with the ab- To search for ISGRI flux variations, we grouped all events sorbed power law model using the C-statistics and obtained Γ +0.9 in three time bins (the maximum numberallowed by the count the photon index =1.1−0.8 and hydrogen column density statistics), but found no significant flux variations. With the +3.3 × 22 −2 NH =5.2−2.4 10 cm . The unabsorbed 0.5-10 keV flux scarce statistics available, we can only exclude flux variations +2.0 × −12 −2 −1 is FX =5.4−1.3 10 erg cm s . greater than a factor of 4 on a timescale of about a year. The combined Swift-ISGRI spectrum can be fitted with a Γ +0.2 +3.2 × 22 single power-law model with =1.3−0.2, NH =6.1−2.2 10 2.2. Swift observation cm−2. Based on the positional coincidence and spectral com- The field of IGR J2018+4043 was observed with Swift patibility, we conclude that we detected a Swift counterpart of XRT (Burrowset al. 2000) on 2006 March 26, 29, and 30, IGR J2018+4043; we will call this source J2018 hereafter. with a total exposure of 4.8 ks (Kennea et al. 2006). The XRT Because of the presence of the bright (mV =5.82) O9V star was operated in photon counting mode. The data were re- HD193322 in the FOV, only W2 and M2 ultraviolet filters duced with the standard HEAsoft package (ver. 6.0.5). The were used with the Swift/UVOT camera. The source was not Hard X-ray source IGR J2018+4043 3 detected in these filters, the upper limits being not constrain- object in the γ-Cygni SNR (d ≈ 1.5 kpc, expected NH ∼ 1– ing given the strong ISM extinction in this direction (see §3). 3 × 1021 cm−2; Landecker et al. 1980). J2018 is likely a more distant object, not associated with γ-Cygni . 2.3. VLA data analysis A number of ISGRI sourceshave been identified as Galactic The radio image at 1.4 GHz was produced from VLA X-ray binaries, mostly HMXBs (e.g., Filliarte & Chaty 2004; archival data obtained in D configurationon 1996 Sep 19. The Tomsick et al. 2006, and references therein). The J2018’s X- 35 2 data were processed with the MIRIAD software following ray spectrum (§2.2) and luminosity, LX ∼ 3×10 (d/10kpc) standard procedures. The source 1411+522 was used as a pri- ergs s−1 in 1–80 keV range, as well as the magnitudes of mary flux density calibrator (S1.4GHz =22.8 Jy) and 1924+334 the IR-optical counterpart (§2.4) do not contradict such an ′′ ′′ as a phase calibrator. The synthesized beam is 54 × 53 interpretation. The apparent ∼10′′ extension of the optical ◦ (beamwidth at half power), P.A. = 33 , and the rms noise is and 2MASS counterparts could be associated with a powerful about 5 mJy/beam. wind of an HMXB’s O or B companion. Alternatively, the ex- Figure 1 shows a contour image of the region around J2018, tended IR-optical source could be interpreted as a galaxy that at 1.4 GHz. The J2018 position is within an extended radio hosts an AGN detected in X-rays. ′ ′′ feature, of about 5 in size, ≈ 70 south-east of the apparent An extrapolation of the X-ray spectrum of J2018 into the position of radio intensity peak. The total flux density of the EGRET range reaches the flux level of GEV J2020+4023, radio feature is about 480 mJy. The resolution in D config- which suggests that these sources could be X-ray and γ- uration is too low to resolve a possible point source within ray counterparts of the same object. An extensive dedicated the structured extended feature. Overall, we cannot rule out search for a counterpart to 3EG J2020+4017 with Chandra the possibilty that the radio source is a counterpart of J2018, and the Green Bank radio telescope in a close vicinity of its although the connection is not very certain yet. nominal position failed to show any clear candidate object (Becker et al. 2004; Weisskopf et al. 2006). However, J2018 2.4. Optical and infrared data analysis is located outside the fields of these searches. Having the position of J2018 well constrained by the Swift The position of J2018 is 30.′6 (35.′5) offset from the nomi- observations, we looked for counterparts of the X-ray source nal position of GEV J2020+4023 (3EG J2020+4017),exceed- in archival optical and IR data. We found a likely counterpart, ing the 8′ (9.′6) 95% localization accuracy in those catalogs with the flux of 88±9 mJy at 8.28 µm, in the 6th MSX cat- (Lamb & Macomb 1997; Hartman et al. 1999). However, ◦ ′ ′′ alog (Egan et al. 2003) at α = 20h18m39.s1, δ = +40 40 56 such offsets are not unusual among EGRET sources. The ′′ (about 2 positional accuracy). We also found an extended position accuracy of relatively faint EGRET sources in the ◦ ′ ′′ NIR source centered at α = 20h18m38.s73, δ = +40 41 00. 1 Galactic plane was estimated by Gehrels et al. (2001) as ∼1◦. ′′ (2. 6 from the Swift position) in the 2MASS XSC catalog Moreover, even the very bright Vela and Geminga pulsars are (Skrutskie et al. 2006), with magnitudes J = 13.0, H = 11.5, located outside their EGRET 99% confidence contours (see Ks = 10.7 (Fig. 2). These magnitudes are brighter than re- Hartman et al. 1999). 3EG J2020+4017 is a bright source, ported by Kennea et al. (2006) because those authors did not but it is located in a very crowded region of the plane, with account for source extension. enhanced, highly inhomogeneous diffuse emission from the The optical images, taken from the POSS-II survey γ-Cygni SNR and nearby molecular clouds, and at least two (Reid et al. 1991), are also shown in Figure 2. A counter- other EGRET sources (3EG J2021+3716 and J2016+3657) part to J2018 can be seen in the POSS-Red (5878-7121 Å) within the 3◦ vicinity, which may have affected the source localization. In such an environment, systematic uncertain- and POSS-IR (6939-9030 Å) filters at the same position as ◦ the 2MASS source, and with a similar spatial extent. ties can increase the localization error up to ∼0. 5. Thus, de- We estimated the fluxes of the extended optical source, us- spite the positional discrepancy, we believe that one can con- ing the magnitudes of two nearby stars from the GSC 2.3 sider the possibility that Swift-ISGRI and EGRET detected catalog (McLean et al. 2000), marked with blue crosses in the same source. Future GLAST (e.g., Gehrels & Michelson panels b and c of Figure 2. We estimated m = 18.1 and 1999) observation could resolve the issue. Red ∼ mIR = 17.3 via aperture count-rate calculations. The absence It is also worth mentioning the 3σ detection of a pos- of the source in the POSS-II Blue image yields an upper limit sible Fe Kα feature at 6.2 keV in the RXTE spectrum of γ- mBlue > 22.5. The aperture boundaries were set at the count Cygni (Bykov et al. 2004). The feature, which could be at- rate corresponding to half a difference between the bright- tributed to a putative hidden accreting source (e.g., an AGN), est source pixel value and the local background level. Two does not appear in the ASCA spectra of the γ-Cygni regions reference stars in each band were used to estimate the zero- observed by this satellite, which did not include the J2018 point and the count-rate vs. magnitude dependence. The position, while this position was well within the RXTE colli- same local background was subtracted from all the count mator field. rates, normalized by the aperture size. To convert the mag- To investigate the nature of J2018, we constructed its spec- nitudes into fluxes, we used the mAB and POSS-II magni- tral energy distribution (SED; Fig. 3) using the ISGRI and tudes of the standard secondary spectrophotometric objects Swift observations and the archival data described above. It HZ4 and G191B2B (Oke 1990) and obtained FBlue < 2.5µJy, should be noted that the data are noncontemporaneous,which FRed =0.27 mJy, FIR = 0.76 mJy for the optical source. We is a point of some concern because the source could be vari- estimate the error of the cross-calibration, about 20%, to be able. The Galactic absorption (AV ≈ 9.8) was accounted for mainly due to the extended nature of the source. in the unabsorbed NIR-optical spectrum using the reddening maps by Schlegel et al. (1998)9. Being rather accurate for 3. DISCUSSION high latitude sources, the method often overestimates the red- The Swift-ISGRI spectrum of J2018 shows that it is a non- dening by a factor of up to 1.5 for strongly reddenned objects thermal, highly obscured source. The hydrogen column den- sity estimated from the Swift data (§2.2) looks too high for an 9 See also http://irsa.ipac.caltech.edu/applications/DUST/ 4 A.M.Bykov et al.

′′ FIG. 2.— Archival POSS-II (top) and 2MASS (bottom) images of the J2018 field. The Swift XRT position is indicated with the 4. 2 radius error circle. The green and blue crosses denote point-like sources from the GSC 2.3 catalog.

(Arce & Goodman 1999). The optical emission is apparently cesses responsible for such a spectrum. On the other hand, extended, but with the data available it is not possible to disen- the leptonic jet model for a microquasar, such as discussed by tangle the core component, corresponding to the point source Dermer & Böttcher (2006) for LS 5039, also provides a SED detected by Swift, from the extended halo. Therefore, the generally consistent with that given in Figure 3. optical points can be lower than shown in Figure 3, and the NIR and optical spectroscopy could distinguish between the extinction corrected (with AV =9.8) optical fluxes in Figure two interpretations by measuring the redshift of the galaxy 3 should be regarded as upper limits, because of a possible that hosts the putative AGN. Also, microquasars are usually overestimation (up to 3 mag) in the AV used. variable on minutes scale, while AGNs are variable on longer The hypothesis that the 1.4GHz VLA excess is due to a time scales making a variability study conclusive. point-like source at the Swift source position (e.g., an AGN) The analysis given above allows us to conclude that the projected onto an extended radio source is not in conflict with properties of the source discovered by Swift are consistent the data available. Therefore, we conditionally include the with those of IGR J2018+4043 . However, that does not ex- VLA point, as well as the EGRET data, in Figure 3. Within a clude a sizable extended hard X-ray flux from the possible in- more conservative approach, they should be regarded as upper teraction of the γ-Cygni SNR with stellar winds as discussed limits. by Bykov et al. (2004). The extended VLA feature seen The radio through hard X-ray SED in Figure 3 is consistent in Figure 1 can also be attributed to the SNR-wind interac- with that expected for a nearby Seyfert AGN or a radio-galaxy tion. High-resolution observations with Chandra and XMM- (e.g., Dermer & Gehrels 1995). For instance, at z =0.02 the Newton, as well as simultaneous multiwavelength studies, are source luminosity below 100 keV would be about 1043 ergs needed to separate the componentsand to make a firm conclu- s−1. sion on the nature of IGR J2018+4043. If GEV J2020+4023 is indeed associated with J2018, then the radio through γ-ray source could be a blazar. Given the AV uncertainties discussed above, the multiwavelength spectrum We thank the anonymous referee whose useful remarks resembles the two-peak synchrotron self-Compton SED of S5 helped us to improve the presentation of the results. Sup- 0716+714, a blazar studied by Foschini et al. (2006). Si- port for this work was provided by the NASA through grant multaneous multiwavelength observations are required, how- NAG5-10865, by RBRF grants 06-02-16844, 04-02-16595, ever, to verify the blazar hypothesis and understand the pro- and by France-Argentina ECOS A04U03.

REFERENCES Arce, H. G., & Goodman, A. A. 1999, ApJ, 512, L135 Buckley, J.H., et al. 1998, A&A, 329, 639 Bassani, L., et al. 2005, ApJ, 634, L21 Burrows, D.N., et al. 2000, Proc. SPIE, 4140, 64 Becker, W., et al. 2004, ApJ, 615, 875 Bykov, A.M., et al. 2004, A&A, 427, L21. Hard X-ray source IGR J2018+4043 5

FIG. 3.— Observed (top) and extinction-corrected (bottom) multiwavelength spectra of J2018. The NIR-optical and soft X-ray points in the bottom panel are 22 −2 corrected for absorption assuming AV = 9.8 (see §3) and NH = 6.1 × 10 cm , respectively. Note that the observational data are noncontemporaneous. The EGRET data points for 3EG J2020+4017 are included, although the association with J2018 is not firmly established (see §3).

Bykov, A.M., et al. 2006, A&A, 449, 917 McLean, B.J., Greene, G.R., Lattanzi, M.G., Pirenne, B. 2000, ADASS IX, Dermer, C.D. & Gehrels, N. 1995, ApJ, 447, 103 ASP Conf. Proc., 216, eds. N. Manset, C. Veillet & D. Crabtree, 145 Dermer, C.D. & Böttcher, M. 2006, ApJ, 643, L13 Oke, J.B. 1990, AJ, 99, 1621 Egan, M.P., Price, S.D., Kraemer, K.E. 2003, AAS Meeting 203, #58.07 Reid, I.N. et al. 1991, PASP, 103, 661 Filliatre, P. & Chaty, S. 2004, ApJ, 616, 469 Schlegel, D.J., Finkbeiner, D.P., Davis, M. 1998, ApJ, 500, 525 Foschini, L et al. 2006, A&A (in press), astro/ph 0604600. Skrutskie, M.F. et al. 2006, AJ, 131, 1163 Gehrels, N., Macomb, D.J., Bertsch, D.L., Thompson, D.J., Hartman, R.C., Sturner, S.J., Dermer, C.D. 1995, A&A, 293, L17 Mattson, B.J., 2001, in: The Nature of Unidentified Galactic High Energy Swanenburg, B.N. et al. 1981, ApJ, 243, L69 Gamma-ray Sources, eds. A.Carrami˜nana et al. (Kluwer: Dordrecht), 81 Tomsick, J., Chaty, S., Rodriguez, J., Foschini, L., Walter, R., & Kaaret, P. Gehrels, N., Michelson, P. 1999, Astropart. Phys., 11, 277 2006, astro-ph/0603810 Goldwurm, A., et al. 2003, A&A, 411, L223 Weisskopf, M.C. et al. 2006, ApJ (submitted), astro-ph/0606596 Hartman, R.C., et al. 1999, ApJS, 123, 79 Winkler, C., et al. 2003, A&A, 411, L1 Kennea, J., et al. 2006, The Astronomer’s Telegram, 788. Landecker, T., Roger, R., Higgs, L. 1980, A&AS, 39, 133 Lamb, R. C. & Macomb, D.J. 1997, ApJ, 488, 880 Lebrun, F., et al. 2003, A&A, 411, L141 Lubinski, P., Dubath, P., Kretschmar, P., Pottschmidt, K., Walter, R. 2004, Proc. 5th INTEGRAL Workshop, ESA SP-552, eds. V. Schönfelder, G. Lichti & C. Winkler (ESA: Munich), 871