Optical-uv spectrum and proper motion of the middle-aged pulsar b1055-52

Draft version July 20, 2010

Preprint typeset using LAT^X style cmulatcapj v. 08/22/09

OPTICAL-UV SPECTRUM AND PROPER MOTION OF THE MIDDLE-AGED PULSAR B1055-52

R. P. Mignani 1 , G. G. Pavlov 2 , O. Kargaltsev 3

Draft version July 20, 2010

o

Oh-
6

CO

>

O

&\
(N

o
o

*

ABSTRACT

PSRB1055— 52 is a middle-aged (r = 535 kyr) radio, X-ray, and 7-ray pulsar showing X-ray thermal
emission from the neutron star (NS) surface. A candidate optical counterpart to PSRB1055— 52 was
proposed by Mignani and coworkers based on Hubble Space Telescope (HST) observations performed
in 1996, in one spectral band only. We report on HST observations of this field carried out in 2008, in
four spectral bands. The astrometric and photometric analyses of these data confirm the identification
of the proposed candidate as the pulsar's optical counterpart. Similarly to other middle-aged pulsars,
its optical-UV spectrum can be described by the sum of a power-law (PLo) component (F u oc v^ a °) 1
presumably emitted from the pulsar magnetosphere, and a Rayleigh- Jeans (RJ) component emitted
from the NS surface. The spectral index of the PLq component, oq = 1.05 ± 0.34, is larger than
for other pulsars with optical counterparts. The RJ component, with the brightness temperature
To = (0.66 ± 0.10)d3 50 i?o 13 MK~ (where ^350 and Ro,i3 are the distance to the pulsar in units of
350 pc and the radius of the emitting area in units of 13 km), shows a factor of 4 excess with respect
to the extrapolation of the X-ray thermal component into the UV-optical. This hints that the RJ
component is emitted from a larger, colder area, and suggests that the distance to the pulsar is smaller
than previously thought. From the absolute astrometry of the HST images we measured the pulsar
coordinates with a position accuracy of 0'.'15. From the comparison with previous observations we
measured the pulsar proper motion, /i = 42 ± 5 mas yr _1 , which corresponds to a transverse velocity
U t = (70±8)d 35 okms- 1 .
Subject headings: pulsars: individual (PSRB1055— 52) — stars: neutron

1. INTRODUCTION

The radio pulsar B1055— 52 (also known as
PSRJ1057-5226) was discovered by Vaughan &
Large (1972). The period, P = 197 ms, and its
derivative, P = 5.83 x 10~ 15 s s _1 , correspond to the
spin-down age r = P/2P = 535 kyr, rotational energy
loss rate (spin-down power) E — 3.0 x 10 34 erg s^ 1 , and
surface magnetic field B = 1.1 x 10 12 G. Neither proper
motion nor parallax measurements have been reported
for this pulsar so far. Based on the pulsar's dispersion
measure, DM=30.1 pc cm -3 , the distance was estimated
as d = 1.53 kpc for the Galactic free electron density
model by Taylor k Cordes (1993), while the NE2001
model by Cordes & Lazio (2002) yields d = 0.73 ± 0.15
kpc.

PSRB1055— 52 was one of the few radio pulsars de-
tected by the Einstein X-ray observatory (Cheng &
Helfand 1983). Observations with EXOSAT showed
a thermal nature of its soft X-ray emission, presum-
ably originating from the neutron star (NS) surface
(Brinkmann & Ogelman 1987), while ROSAT observa-
tions found that at least two spectral components were
needed to fit the spectrum (e.g., a blackbody [BB] and a
power-law [PL]) and discovered X-ray pulsations at the

radio period (Ogelman & Finley 1993). The similarities

Mullard Space Science Laboratory, University College Lon-
don, Holmbury St. Mary, Dorking, Surrey, RH5 6NT, UK;
rm2@mssl.ucl.ac.uk

2 Department of Astronomy and Astrophysics, Pennsylvania
State University, PA 16802, USA; pavlov@astro.psu.edu

3 Department of Astronomy, University of Florida, FL 32611,
USA; oykl00Oastro.ufi.edu

in the ages and X-ray properties with two nearby middle-
aged pulsars, PSRB0656+14 (r = 110 kyr) and Geminga
(r = 340 kyr), for which ROSAT also detected thermal
X-ray emission, earned these three pulsars the nickname
of "The Three Musketeers" (Becker & Triimper 1997).

The X-ray spectrum and pulsations of PSRB1055-
52 were also studied with ASCA (Greiveldinger et al.
1996; Wang et al. 1998) and BeppoSAX (Mineo et al.
2002). Recent observation with the more sensitive Chan-
dra (Pavlov et al. 2002) and XMM-Newton (De Luca
et al. 2005) have shown that the phase-integrated X-
ray spectrum, similar to those of the other two "Mus-
keteers", can be fitted with three spectral components:
a cold BB (T c ~ 0.8 MK), a hot BB (T H « 1.8 MK),
and a PL component (photon index T w 1.7, luminosity
L nonth _ 8 x io30rf2 50 er g s" 1 in the 0.5-10 keV band,
where d^ 50 = d/(750pc)), presumably emitted from the
pulsar magnetosphere. Phase-resolved spectroscopy of
the Chandra and XMM-Newton data suggests that the
changing projected emitting area of the hot component
is responsible for the bulk of the phase variations (Pavlov
et al. 2002; de Luca et al. 2005).

PSRB1055— 52 was one of the seven 7-ray pulsars de-
tected by the EGRET instrument aboard the Compton
Gamma-ray Observatory (CGRO) (Fierro et al. 1993;
Thompson et al. 1996). It has also been studied by
the Fermi LAT detector (Abdo et al. 2010). Its 0.1-100
GeV luminosity, I 7 ~2x lO 34 ^ 2 -^ erg s _1 (assuming an
approximately isotropic emission), implies an extremely
high 7-ray efficiency, n^ = L 1 jE ~ 0.6d^ 50 .

To study and interpret the magnetospheric and ther-
mal emission from PSRB1055— 52, the X-ray and 7-ray

Mignani, Pavlov, & Kargaltsev

data should be supplemented by optical data. How-
ever, optical observations with ground-based telescopes
(Manchester et al. 1978; Cheng & Hclfand 1983; Big-
nami et al. 1988) failed to detect PSRB1055-52 because
of the presence of a bright (V w 14) F-type field star
w 4" from the radio pulsar position (dubbed Star A by
Manchester et al. 1978). Thanks to the sharp resolution
and high sensitivity of the Faint Object Camera (FOC)
aboard the Hubble Space Telescope (HST), Mignani et
al. (1997; hereafter M+97) were able to identify a faint
(U ~ 25) probable counterpart at the radio pulsar po-
sition. Since the observation was taken in only one fil-
ter, the shape of the optical spectrum of the candidate
counterpart remained unknown, and even its identifica-
tion with the pulsar was not certain. To confirm the
identification, measure the spectral shape (in particular,
separate the magnetospheric and thermal components,
similar to the two other Musketeers - see, e.g. Kargal-
tev & Pavlov 2007; KP07 hereafter), and measure the
optical pulsations, new HST observations were obviously
required. Our proposal to reobserve PSRB1055— 52 with
HST was accepted for Cycle 16 (program 11154). These
observations are described in Section 2, while the results
are presented in Section 3, discussed in Section 4, and
summarized in Section 5.

2. OBSERVATIONS AND DATA ANALYSIS

As no timing capabilities were available in the HST Cy-
cle 16, we could only image the target in several filters.
To observe PSRB1055-52 in the far-UV (FUV), most
useful for studying the thermal component, we employed
the Solar Blind Channel (SBC) detector of the Advanced
Camera for Surveys (ACS). The SBC is a Csl Multi-
Anode MicroChannel Array (MAMA) photon-counting
detector, which is operated in ACCUM mode producing
time-integrated images. The detector is sensitive to ra-
diation in the 1115-1700 A wavelength range, and it pro-
vides a nominal field of view (FOV) of 34'.'6 x 30'.T, with a
pixel scale of 0'.'034 x 0'.'030. To filter out the geocoronal
lines, we used the longpass filter F140LP (pivot wave-
length A = 1528 A; equivalent gaussian FWHM AA =
297 A), which cuts off radiation at A < 1350 A. The
pulsar was observed on 2008 February 13 over two con-
secutive spacecraft orbits for a total integration time of
5569 s. The data were reduced and calibrated (including
the linearity correction, dark subtraction, flat-fielding)
through the ACS data reduction pipeline (CALACS), using
the closest-in-time calibration frames. Geometric distor-
tion correction, cosmic-ray filtering, and exposure stack-
ing were applied during the pipeline processing using the
Multidrizzle task. This task also applied a re-sampling
of the drizzled image to an even pixel size of 0'.'025.

We also observed PSRB1055— 52 in three optical bands
with the Wide Field and Planetary Camera 2 ( WFPC2) 4
on 2008 March 8. To exploit the maximum spatial reso-
lution, the pulsar was centred in the Planetary Camera
(PC) chip, which has a pixel scale of 0'.'045 and a FOV
of 35" x 35". To maximize the spectral coverage, we
observed through the F450W (A = 4557 A; AA = 951

4 We originally proposed the observations in the optical bands
for the higher throughput ACS Wide Field Channel (WFC), but
we had to switch to the WFPC2 since the WFC was turned off in
January 2007 because of problems with the ACS electronics.

A), F555W (A = 5442 A; AA = 1229 A), and F702W
(A = 6917 A; AA = 1381 A) filters. Observations were
carried out during four consecutive orbits for the total
integration times of 3600 s in the F450W filter and 1800
s in each of the F555W and F702W filters. For each
filter, the observations were split in shorter, dithered ex-
posures of 600 s to filter out cosmic ray hits and bad
pixels. The data were processed through the WFPC2
CALWP2 reduction pipeline for bias, dark, and flat-field
correction and flux calibration. For each filter, we then
combined single exposures with the IRAF task drizzle
to produce co-added and cosmic-ray-free images.

We also used the archival FOC observations by M+97.
They were carried out on 1996 May 11 in the F342W
filter (A = 3402 A; AA = 519 A) during three spacecraft
orbits with the exposure times of 2771 s for the first orbit
and 3062 s for the second and third orbits, corresponding
to the total integration time of 8895 s. The data were
obtained with the high-resolution F/96 relay, with the
FOV of 7" x 7" and pixel scale of 07014. The data were
retrieved from the HST science data archive 5 and on-the-
fly recalibrated through the CALF0C reduction pipeline
with the updated reference files.

3. RESULTS

3.1. Detection of the candidate counterpart

The most recent published radio coordinates of
PSRB1055-52 (Newton et al. 1981),

a = 10 h 57 m 58?84 ± 0?03, 5 = -52°26'5673 ± 073, (1)

are for the epoch of 1978.13 (equinox J2000). These co-
ordinates, without proper motion values, are reported
in the ATNF radio pulsar catalogue 6 (Manchester et al.
2005), and they were also used by Weltevrede et al.
(2010) to compute the PSRB1055— 52 radio ephemeris
for the timing of the Fermi 7-ray observations. Since
pulsars are known to have high velocities, PSRB1055— 52
could move a few arcseconds in the 30 year span between
the epochs of the radio position and our HST observa-
tions (2008.12 and 2008.18). For instance, for an aver-
age radio pulsar velocity of 400 km s _1 (e.g., Hobbs et
al. 2005), we would expect a proper motion of 110 dj 50
mas yr _1 and a position shift of 374 d? 50 in 30 years, in
an unknown direction. Such an uncertainty hampers a
direct identification based on the pulsar's position only.
Fortunately, the presence of Star A in the immediate
vicinity of the pulsar can be used for the pulsar identifi-
cation.

The 10" x 10" cutout of the SBC image of the
PSRB1055-52 field (left panel of Figure 1) shows the
only two objects detected in the entire SBC FOV, sepa-
rated by a distance of 4741. From the comparison of this
image with that obtained by M+97 with FOC, it is nat-
ural to assume that the northwestern and southeastern
objects are Star A and the candidate pulsar counterpart,
respectively.

Indeed, the nominal SBC coordinates of the northwest-
ern source are a = 10 h 57 m 58?790 and 5 = -52°26'53"04.
The most accurate coordinates of Star A, cxa —

See http://archive.stsci.edu.
6 See http://www.atnf.csiro.au/research/pulsar/psrcat/.

HST observations of PSRB1055-52

10 h 57 m 58 i :728±0!004 and 5 A = -52°26'52'.'45±0'/06, at
the epoch of our HST observation, can be derived from
the Third US Naval Observatory CCD Astrograph Cata-
log (UCAC3; Zacharias et al. 2010), with account for the
Star A's proper motion, \x a = —8.9 ±1.3 mas yr _1 and
fig = 7.8 ± 2.2 mas yr -1 , given in the same catalog and
independently verified by us (see Appendix). As the dif-
ference between these and the nominal SBC coordinates
(0'.'57 and —0759 in right ascension and declination, re-
spectively) is within the error budget of the HST abso-
lute astrometry, we conclude that the SBC northwestern
source is indeed Star A.

Applying the boresight correction to the SBC image,
such that the coordinates of the northwestern source co-
incide with the UCAC3 coordinates of Star A, we obtain
the following coordinates of the southeastern source

10 h 57 m 58!954,

-52°26'56'/37,

(2)

with a nominal position uncertainty of about 0'.'07.
This uncertainty is mostly determined by the uncer-
tainty of the UCAC3 absolute astrometry, thanks to
the small angular distance (4'.'41) between Star A and
PSRB1055-52, and the - 0.1° accuracy on the HST
roll angle with respect to the equatorial reference frame.

The derived position is consistent with the candi-
date pulsar counterpart position in the FOC observa-
tion (a = 10 h 57 m 58?832, 5 = -52°26'56"28), within the
uncertainty of the FOC position (w l'.'O) 7 and the un-
certainty caused by an unknown proper motion since the
epoch of the FOC observation (1996.36). This suggests
that we indeed detected the PSRB1055— 52 candidate
counterpart of M+97 in our SBC image, while the ex-
tremely blue color (obvious from the comparison of the
FUV (SBC) image with the optical (FOC and WFPC2)
images; see Figure 1 and also Section 3.2) virtually proves
that this object is indeed PSRB1055-52.

Using the offsets between the pulsar counterpart and
Star A measured in the SBC image (2702 and -3792 in
right ascension and declination, respectively), we looked
for the pulsar counterpart at the corresponding posi-
tion in the WFPC2 images. We detected it in both the
F702W image (seen in the circle labeled "SBC" in the
right panel of Figure 1) and in the F555W image but not
in the F450W one.

3.2. Photometry

We measured the magnitudes and fluxes of the
PSRB1055-52 counterpart in the WFPC2, ACS/SBC,
and FOC images through customized aperture photome-
try using the IRAF package digiphot. For the WFPC2
F555W and F702W bands, we used a small source aper-
ture with the radius r = 3 pixels (07135) to maximize the
signal-to- noise ratio (S/N) . We sampled the background
in a surrounding annulus with the inner radius rt, = 5
pixels and radial width 5r = 10 pixels; the distance of
2 pixels from the source aperture was chosen to exclude
the brighter portion of the point spread function (PSF)

7 Before 2002, the pair of guide stars used to compute the as-
tromctric solution in the HST focal plane were selected from the
GSC1.0 (Lasker et al. 1990), which had a mean positional error
of about l'.'O. After 2002, the guide stars arc selected from the
GSC2.3 (Lasker et al. 2008), which has an improved mean posi-
tional error of about 0'.'3.

wings. We then applied the aperture correction using the
fractional encircled energy (FEE) coefficients (Holtzman
et al. 1995), and the correction to compensate for the
time- and position-dependent charge transfer efficiency
(CTE) losses of the WFPC2 detector (Dolphin 2009).
For estimating the 3 a flux upper limit in the F450W
band, we used the same r = 3 pixels source aperture
but sampled the background in a wider annulus, Sr = 30
pixels, around the expected pulsar position.

For the SBC data, we chose a source aperture of r = 6
pixels (072) and a background annulus with r\> = 20 pix-
els and Sr = 10 pixels. The large radial distance of 14
pixels between the source aperture and the background
annulus was chosen to exclude the count excess over the
PSF wings, which is seen from the comparison of the
count distributions around the pulsar and around Star
A (see the insets in the left panel of Figure 1). The
3 a excess (e.g., 114 ± 38 counts in the 073-077 annulus)
cannot be ascribed to different PSFs of the pulsar and
Star A because the latter is only one magnitude brighter
than the pulsar in the SBC F140LP band, and it is sep-
arated from the pulsar by only 4741. From these data
alone, we cannot exclude the possibility that the excess
is a compact pulsar wind nebula (PWN), but, as the
excess is not seen in other bands, it can also be just a
small-scale enhancement of the detector background at
the pulsar position. Since the brightness profile of Star
A is closer to that of the model PSF, we used it as a ref-
erence to compute the aperture correction by fitting its
growth curve (i.e., the number of source counts within
an aperture as a function of its radius).

We also measured anew the magnitude and flux of the
counterpart in the FOC image. Following the approach
described by Pavlov et al. (1997), we chose an optimal
source aperture from the measured growth curves for
each of the three orbits. We found that the growth curves
for the first and third orbits were consistent with each
other and compatible with the nominal values of FEE
at w 3400 A, while the growth curve for the second or-
bit showed a peculiar behavior (e.g., no saturation up to
r ~ 20 pixels, which means an unusually broad PSF).
Therefore, we excluded the second orbit data from the
analysis and chose r — 8 pixels (0711), r& = 10 pixels,
and Sr = 10 pixels.

In all cases, we applied the countrate-to-flux conver-
sion by computing the photometric zero points in the
STmag system using the image keywords PH0TFLAM and
PH0TZPT, derived by the HST data reduction pipelines.
To calculate the extinction-corrected fluxes, we esti-
mated the interstellar reddening toward the pulsar,
E(B — V) = 0.07, using the hydrogen column density
N H = 2.7 x 10 20 cm" 2 , derived by De Luca et al. (2005)
from the X-ray fits, and using the correlation between
N H and E(B - V) found by Predehl & Schmitt (1995).
We then calculated the extinction coefficients using the
optical extinction curves by Fitzpatrick (1999) for the
WFPC2 and FOC passbands and the UV extinction
curves by Seaton (1979) for the SBC passband. The mag-
nitudes and the measured, F° hs , and dereddened, F^ er ,
spectral fluxes are presented in the last two columns of
Table 1.

We plotted the dereddened spectral fluxes F^ er ver-
sus frequency v in Figure 2. The shape of the broad-

Mignani, Pavlov, & Kargaltsev

band spectrum and, in particular, the relatively high
FUV brightness of the object are inconsistent with it
being a non-collapsed star in the Galaxy or an extra-
galactic object. Therefore, the photometry proves un-
equivocally that the source is the optical counterpart of
PSRB1055-52.

It is clear from Figure 2 that a simple one-component
spectral model cannot fit the optical-UV spectral energy
distribution (SED). However, similar to the optical-UV
spectra of other middle-aged pulsars (e.g., KP07), the
SED can be fitted by a model consisting of an (optical)
power-law (PLo) component and a Rayleigh- Jeans (RJ)
component:

iC r = k

VQ

27r R l hT 2

(3)

where i?o is the radius of equivalent emitting sphere,
and To is the brightness temperature (both as seen by
a distant observer). At E(B - V) = 0.07, the fitting
parameters are a = 1.05 ±0.34, f = 0.112 ±0.003 /zJy
at i/ = 1 x 10 15 Hz, and R 2 T = (510 ± 81)d? 50 km 2
MK (or T = (3.02 ± 0.48) df 50 Rq 2 13 MK > where #0,13
is the radius in units of 13 km). The PLo component,
which can be interpreted as magnetospheric radiation,
dominates in the optical (A > 3000 A). The energy flux
of this component is J-o° nth w 1.3 x 10 -15 erg cm -2 s _1
in the 3000-9000 A band. The RJ component dominates
in the FUV and is likely emitted from the NS surface.
We will discuss implications of these results in Section
4.1.

3.3. Astrometry
3.3.1. Absolute position of the pulsar

The pulsar counterpart position given by Equation (2)
was obtained from the SBC astromery using the Star A
position in the UCAC3 catalog. Of course, using just
one reference star may introduce a significant bias in the
absolute target position. Fortunately, the WFPC2 has
a much larger FOV than the SBC, allowing astrometry
with larger numbers of reference stars.

For the WFPC2 astrometry, we used two catalogs,
GSC2.3 (Lasker et al. 2008) and 2MASS (Skrutskie et
al. 2006). Although the astrometric accuracy of these
catalogs is somewhat lower than that of UCAC3, they
contain more stars per given area, which improves the
statistics of the astrometric fits. Moreover, most UCAC3
stars identified in the WFPC2 FOV are saturated, which
prevents their use for the astrometric calibration. We
used the F702W image, in which PSRB1055-52 was de-
tected with a higher S/N than in the F555W image. We
produced the WFPC2 mosaic image with the STSDAS
task wmosaic, which also corrects for the geometric dis-
tortion. Then, we measured the centroids of the reference
stars (13 stars in GSC2.3 and 16 stars in 2MASS) in the
WFPC2 pixel coordinates through gaussian fitting with
the Graphical Astronomy and Image Analysis (GAIA)
tool 8 and used the catalog sky coordinates of these stars
to compute the pixel-to-sky coordinate transformation
with the code ASTROM 9 . This yielded astrometric fits with

8 See star-www.dur.ac.uk/~pdraper/gaia/gaia.html.

9 See http://star-www.rl.ac.uk/Software/software.htm.

the rms values a r — 0'.'17 for both the GSC2.3 and
2MASS stars. To these we added in quadrature the un-
certainties of the registration of the WFPC2 image on
the catalog reference frames, Ctr = 0'.'14 and 0'/08, for
the GSC2.3 and 2MASS, respectively. Following Lat-
tanzi et al. (1997), these uncertainties were estimated
as er t r = (3/ J V cat ) 1 / 2 CT cat , where 3 1 / 2 accounts for 3 free
parameters in the astrometric fit, <7 ca t is the mean posi-
tional error of the catalog coordinates (0'.'3 and 0'.'2 for
the GSC2.3 and 2MASS, respectively), and A cat is the
number of catalog stars used to compute the astrometric
solution. Accounting for the 0'.'15 and 0'.'015 uncertain-
ties on the link of the GSC2.3 and 2MASS to the Inter-
national Celestial Reference Frame (ICRF), the overall
uncertainties of the astrometry (1 a position errors) are
S r = 0726 and 0719, respectively.

We applied these astrometric solutions to compute the
sky coordinates of PSRB1055— 52 from its pixel coordi-
nates. As it has been done for the reference stars, we
computed the pixel coordinates through a gaussian fit-
ting, with an uncertainty of < 0701 for the pulsar cen-
troid, negligible in comparison with the overall uncer-
tainty of our absolute astrometry calibration. We then
obtained the GSC2.3-based pulsar coordinates

a = 10 h 57 m 58?961, 8 = -52°26'56717, (4)

with the position error <5 r = 0726, and the 2MASS-based
coordinates

a = 10 h 57 m 58!967, 8 = -52°26'56731,

(5)

with 8 r = 0719. These two positions are consistent with
each other and with the UCAC3-based position derived
by using Star A as the only reference star in the SBC
astrometry (Equation 2).

Since the coordinates in Equations (4) and (5) have
been computed using different catalogs and different ref-
erence stars, they can be considered independent and
hence can be averaged. The calculation of the weighted
means yields

10 h 57 m 58!965, <5 = -52°26'56"26,

(6)

with S r = 0715. Thus, the WFPC2 astrometry, re-
calibrated with the GSC2.3 and 2MASS, has provided
the absolute pulsar coordinates at the epoch of 2008.18,
which is separated by about 30 years from the epoch of
the only published radio position 10 .

3.3.2. Proper motion of the pulsar

We used the optical coordinates of PSRB1055— 52 to
measure its proper motion from the comparison with
the radio coordinates at the 1978.13 epoch (Equation
1). We note that, since the GSC2.3 and 2MASS, as
well as UCAC3, are linked to the ICRF, there is no a
systematic offset between the optical and radio coordi-
nates. Using Equation (6), we infer the displacement
Aa cos 5 = +1714 ±0729 and AS = +0704±0732, which
corresponds to the proper motion

/j, a = +38 ± 10 masyr -1 , /j,g = +1 ± 11 masyr -1 (7)

10 After completing our astrometric analysis, we became aware
that an unpublished radio timing position was recently obtained
(R. N. Manchester, private communication), which we found con-
sistent with our own results.

HST observations of PSRB1055-52

in right ascenscion and declination. Using the UCAC3-
based positions (Equation 2) results in virtually the same
proper motion (/x Q = 35 ±9 mas yr" 1 , [j,$ = — 2± 10 mas
yr -1 ), whose errors are dominated by the errors of the
radio position (Equation 1).

The pulsar proper motion can be independently mea-
sured through relative HST astrometry, by comparing its
position with respect to Star A in the SBC image (epoch
2008.12) with that in the FOC image (epoch 1996.36).
Of course, measuring a relative proper motion with only
one reference star comes with caveats. Firstly, Star A
is slightly saturated in the FOC image and part of its
PSF is outside the FOC FOV (see Figure 1 of M+97),
which lowers the precision of the star centroid determina-
tion. To be very conservative, we assumed an uncertainty
of 5 pixels (0'.'07). The pulsar centroids were measured
with a much higher precision of < 0.1 pixels in both
the SBC and FOC images (< 0"0025 and < O'.'OOW, re-
spectively) , and the residual geometric distortion (about
/ /004 and 0'.'007 rms for the SBC and FOC, respectively
- see Apellaniz & Cox 2008 and Nota et al. 1996, respec-
tively) is also much smaller than the FOC centroiding
uncertainty for Star A.

Secondly, with only one reference star, we have to rely
upon the nominal values of the angles between the de-
tector axes and the directions of right ascension and dec-
lination, which may introduce an error in the direction
of the proper motion. However, this error is very small
in our case, thanks to the accuracy of the detector posi-
tion angle and the small separation of Star A from the
pulsar (see §3.1). Last but not least, the uncertainty of
the proper motion of the reference star should be small
enough not to hamper the measurement of the relative
pulsar proper motion. We carefully verified it for Star A
(see Appendix).

Using the value of the pixel scale for the two detec-
tors and the nominal orientations of the two images with
respect to the equatorial reference frame, we found the
displacement of the pulsar with respect to Star A during
the time span of 11.76 years between the two HST obser-
vations: Aa cos<5 = 0'/61±0'.'14 and AS = -0'.'19±0'/14.
The uncertainty of the displacement is dominated by the
large uncertainty of the Star A centroid in the FOC im-
age. This displacement corresponds to the proper motion
fi a = 52 ± 6 mas yr -1 and jj,s = —16 ± 6 mas yr . Cor-
recting these values for the proper motion of star A as
given in the UCAC3 catalog (see § 3.1), we obtained

Ha = 43 ± 6 mas yr

Us = —5 ± 6 mas yr

(8)

This proper motion agrees with that measured from the
comparison of the radio and optical absolute coordinates.
The weighted mean of the proper motion values deter-
mined by the two methods is

/j, a = 42 ± 5 mas yr 1 , /is = — 3 ± 5 mas yr

-i

or

^ = 42 ±5 masyr -1 , P.A. = 94° ± 7 C

(9)

(10)

for the total proper motion and position angle (counted
east of north). The large proper motion provides addi-
tional evidence that the M+07 candidate is indeed the
pulsar counterpart.

4. DISCUSSION
4.1. Pulsar spectrum

We have shown in Section 3.2 that the optical-UV SED
for PSRB1055— 52 can be described by a sum of PLo
(magnctosphcric) and RJ (thermal) spectra (see Equa-
tion 3). It is interesting to compare these optical com-
ponents with the corresponding spectral components at
higher energies, particularly in soft X-rays where both
the thermal and magnetospheric components are seen.
As we have mentioned in Section 1, the pulsar's X-ray
spectrum can be described by a model that consists of
cold BB (BBc), hot BB (BB H ), and X-ray power-law
(PLx) components. The BBc an d BBh components
are presumably emitted from the bulk of the NS sur-
face and polar caps, respectively, while the PLx com-
ponent, F v = fx(E/Eo)~ a *, represents the magncto-
sphcric X-ray emission. The parameters of the three
components, as inferred from the XMM-Newton data by
De Luca et al. (2005), are the following: T c = 0.79 ±0.03
MK, R c = 12.3to 5 7 d 750 km; T H = 1.79 ± 0.06 MK,
Ru = (0.46 ± 0.06) d 7 50 km; and a x = 0.7 ± 0.1,

fx = L3i"o.i x 10 ~ 31 cr g cm ~ 2 s_1 Hz_1 at S = 1
keV, respectively. These components and their extrapo-
lations into the optical-UV domain are shown in the left
panel of Figure 3.

We see from Figure 3 that the extension of the best-fit
PLx component into the optical overshoots the optical
fluxes by a factor of ~ 3, similar to most of other pulsars
detected in the optical (e.g., KP07; Mignani et al. 2010).
However, in contrast to the majority of those pulsars,
the PLo component is apparently steeper than the PLx
component (ao = 1.05±0.34vs. ax = 0.7±0.1), and per-
haps steeper than the PLo components in other pulsars
(e.g., a = 0.41 ±0.08 and 0.46±0.12 in PSR B0656+14
and Geminga, respectively; see KP07). Moreover, the ex-
trapolation of the best-fit PLo component into the X-ray
range docs not intersect the PLx spectrum, being at least
a factor of ~ 10 lower at E > 1 keV. Such behavior is
unusual (a notable exception is PSR B0540-69; Mignani
et al. 2010), and it might suggest that different popula-
tions of ultrarelativistic electrons are responsible for the
optical and X-ray emission of PSRB1055— 52. However,
taking into account the large uncertainties of the spectral
slopes, it seems more plausible that ao is substantially
smaller than its best-fit value (e.g., ao J$ 0.5), in which
case the extrapolation of the PLo component can, in
fact, be smoothly connected with the PLx spectrum.

Interestingly, the slope a 7 = 0.06 ± 0.1 of the Fermi
LAT spectrum, fit by a PL with an exponential cut-
off (F v ex E- a - ■ cxp(- E/E cut ); Abdo ct al. 2010), is
even flatter than the X-ray slope. On the other hand,
as we see from the right panel of Figure 3, the optical
and 7-ray points can be connected by a PL spectrum
with the slope ao 7 ~ 0.46, which, however, goes slightly
above the PL tail of the X-ray spectrum. We cannot rule
out the possibility that further observations (especially
in the IR-optical), supplemented by a joint multiwave-
lcngth analysis, will show the magnetospheric spectrum
of PSRB1055— 52 to be similar to those of other pulsars
with optical, X-ray, and 7-ray counterparts, including
the other two Musketeers.

Another property of PSRB1055— 52 to compare with

6

Mignani, Pavlov, & Kargaltsev

those of other pulsars detected in the optical is the re-
lationship between the optical and X-ray luminosities of
the magnetospheric emission. For instance, for the wave-
length range 4000-9000 A, the flux in the PLq compo-

nent, Jq

•nonth

0.94 x 10 15 erg cm 2 s 1 , corresponds

to the luminosity Lq

nonth

47rd 2 J"g onth = 6.3 x 10 28 d? 50
erg s" 1 and optical efficiency r/o = L^^/E = 2.1 x
10~ 6 <i 2 50 . At the adopted distance of 750 pc, both the
luminosity and efficiency of PSRB1055— 52 are consider-
ably higher than those of PSR B0656+14 and Geminga,

4.2 x 10" 7 (i 2

for which L^ onth = 1.6 x 10 28 d 2

290 J ^O

<290-

and L^ onth = 4.8 x lO 27 ^, tj = 1.5 x lO" 7 ^, re-
spectively (KP07). This hints that the actual distance
to PSRB1055— 52 is smaller, which is supported by the
analysis of the thermal emission below. On the other
hand, the ratio of the optical (4000-9000 A) and X-
ray (1-10 kcV) nonthermal luminosities, L™ nth /L™ nth =
0.9 x 10~ 2 (which does not depend on the distance), is
within the same range of ~ 10~ 3 -10 -2 as for all other
pulsars with optical counterparts (Zavlin & Pavlov 2004;
Zharikov et al. 2006).

Let us now compare the optical-UV and X-ray ther-
mal components of the PSRB1055— 52 spectrum. The
BBh component is not seen in the optical because of
the small emitting area, while the extrapolation of the
BBc into the optical-UV is a factor of « 4 below the
observed RJ component. Such an "optical-UV excess"
of the RJ component has been seen in the nearby radio-
quiet isolated NSs (RQINSs), such as RXJ1856. 5-3754
and RX J0720.4-3125 (e.g., Kaplan 2008, and references
therein), and in the millisecond pulsar J0437— 4715 (Kar-
galtsev et al. 2004). It, however, has not been observed
in middle-aged pulsars. In particular, the extrapolated
BBc component lies slightly above the optical-UV RJ
component in Geminga (Kargaltsev et al. 2005), while
these components virtually coincide with each other in
PSR B0656+14 (KP07).

The nature of the optical-UV excess is not fully under-
stood; a popular hypothesis is that it is due to a nonuni-
form distribution of the temperature over the NS surface,
such that the BBc component originates from a smaller
and hotter region than the optical-UV RJ component
(e.g., Pavlov et al. 2002). Such an assumption is crudely
equivalent to adding a "very cold" thermal component,
BBycj which dominates in the UV but makes very lit-
tle contribution in X-rays. In other words, the product
RqTo in the second term of Equation (3) is approxi-
mately equal to RqTq + Ry C Tyc- Taking into account

SlOc^o

that RlT c « 119c&= n km 2 MK and R^Ti

l c J c

';-,{)

km 2 MK, we obtain R? r ~T*

vc-'ve

391d? 50 km 2 MK (or
Tyc ~ 2.31 d 2 50 i?y^ 13 MK). Unfortunately, we cannot
measure Tyc independent of d /Ry C in the RJ regime.
We can, however, estimate an upper limit on Tyc from
the requirement that the BByc spectral flux extrapo-
lated to soft X-ray energies becomes so small, in com-
parison with the BBc, that it does not affect the X-ray
fit. We found that it occurs at T vc < 0.45 MK (at this
temperature the contribution of the BByc flux is < 10%
above 0.3 keV - see Figure 3, left).

Such a limit on Tyc implies a considerably smaller dis-
tance to the pulsar than that estimated from the pulsar's

dispersion measure. Indeed, from the above estimate on
-Ryc^vc we obtain d < 25.4(i? vc /lkm)pc < 330i? N s,i3
pc, where we took into account that i?vc is smaller than
the NS radius, .Rns = 13-Rns,13 km. Since the NS ra-
dius (as seen from infinity) can hardly exceed 20 km, the
distance is < 500 pc, and it can even be substantially
smaller than this upper limit.

We should note that such estimates are, of course,
model dependent. For instance, the temperature distri-
bution over the NS surface can be smooth (i.e., there
are no three distinct regions as we implicitly assumed
above). Moreover, the spectrum of the thermal emission
from a NS surface region can differ from the BB, which
may lead to an excess of the actually emitted optical-UV
thermal component over the extrapolation of the X-ray
thermal component even if the surface temperature is
uniform (e.g., in the case of a light-element, H or He,
NS atmosphere; Pavlov et al. 1996). However, fitting the
model light-clement atmosphere spectra to the observed
X-ray spectrum always leads to a lower temperature and
a larger value of R/d (i.e., to a smaller distance at a
given NS radius) than those obtained from the BB fits
(Pavlov et al. 1995). Therefore, it seems hard to avoid
the conclusion on a smaller distance 11 , although its value
cannot be determined accurately from the data available.
We believe 200 pc and 500 pc are reasonable lower and
upper limits, and we will scale the distance to 350 pc
below.

The conclusion that the distance is smaller than previ-
ously thought has important implications. For instance,
the radius of the BBc emitting region, Rq = 5.7^0.3^350
km, becomes substantially smaller than a plausible NS
radius. It means that T c = 0.79 ± 0.03 MK is not the
temperature of the entire NS surface, and it should notbe
used for comparisons with the NS cooling models. This
removes the problem of PSR B1055— 52 being too hot and
luminous (£c,boi = 4.2 x 10 32 d 2 50 erg s^ 1 ) for its age, in
comparison with other middle-aged NSs and the predic-
tions of standard NS cooling models (e.g., Yakovlev &
Pethick 2004). It also means that PSRB1055-52 is not
a "very slowly cooling NS" (as suggested by Kaminker et
al. 2002), and there is no need to invoke an unusually low
NS mass to explain its thermal emission (see Yakovlev &
Pethick 2004 for references). To infer the true thermal
luminosity and average surface temperature, and com-
pare them with the predictions of the NS cooling mod-
els, the parallax-based distance to the pulsar should be
measured, and the distribution of the temperature over
the NS surface should be determined from the phase-
resolved spectral analysis of the pulsar's soft X-ray and
FUV emission.

Another implication of the smaller distance is that the
magnetospheric luminosities and efficiencies in various
spectral bands are lower than previously thought, being
close to the luminosities and efficiencies of other pulsars.
For instance, if d ~ 350 pc, its 7-ray efficiency, ?7 7 ~
0.13 d 2 35Q in the 0.1-100 GeV band (Abdo et al. 2010),

11 Note a similar discrepancy between the DM-based distance
estimate (750 pc, in the model by Taylor & Cordes 1993) and
the distance determined from parallax measurements (288_ 27 pc;
Brisken et al. 2003) for PSR B0656+14. Interestingly, Anderson
et al. (1993) estimated the distance to be 280_ 50 pc using the NS
atmosphere models by Shibanov et al. (1993) and ROSAT data.

HST observations of PSRB1055-52

7

although still rather high, becomes not so different from
the efficiencies of other pulsars detected with Fermi LAT,
and the optical efficiency, 770 = 4.6 x 10~ 7 <i3 50 , becomes
comparable with those of PSRB0656+14 and Geminga.

4.2. Pulsar astrometry and kinematics

PSRB1055— 52 is the third rotation-powered pulsar,
after Geminga (Caraveo et al. 1998) and PSRB0540-69
(Mignani et al. 2010), for which an absolute position has
been determined with a sub-arcsecond accuracy using
optical astrometry techniques. Knowing the astrometric
position will be useful for the precise timing analysis of
the pulsar.

The proper motion we measured for the pulsar's optical
counterpart (Equation 11) corresponds to a tangential
velocity V t — (70 ± 8)^350 km s . This value is well
below the average 400 km s _1 for radio pulsars (Hobbs
et al. 2005), but it is close to that of the Vela pulsar
(V t ~ 65 km s" 1 ; Caraveo et al. 2001).

The proper motion in Galactic coordinates,

[ii = 39 ± 5 masyr , /if, = 15 ± 5 masyr

(11)

shows that the pulsar, whose current coordinates are
I = 285.984° and b = +6.649°, is moving away from the
Galactic plane. We attempted to use the proper motion
to locate the pulsar's birth place, identify a parent star
cluster (OB association), and estimate the actual age of
the pulsar that can differ from its spin-down age. How-
ever, to calculate the pulsar's trajectory back in time in
the Galactic potential, one should know the present dis-
tance (which is poorly known), the tangential velocity
(proportional to the distance), and the radial velocity
V r (which is unknown). Therefore, one has to calculate
a large set of pulsar trajectories for various values of d
and V r . Using the code developed by Vande Putte et
al. (2010), we calculated the trajectories on a grid of 21
distances (in the range of 200 to 1200 pc) and 21 radial
velocities (—500 to +500 km s _1 ). We also selected a
sample of young (< 100 Myr) clusters/associations (de
Zeeuw et al. 1999; Dias et al. 2002) in a 25° radius cir-
cle around the present position of the pulsar, extrapo-
lated back in time their motion using the same code,
and looked for closest approaches of the pulsar and clus-
ter trajectories. The results, however, turned out to be
rather ambiguous because different values of d and (es-
pecially) V r resulted in quite different candidate parent
associations and pulsar ages (e.g., we found 10 candi-
dates for the age range of 400-700 kyr) , and the uncer-
tainty of V t propagated back in time resulted in very
large uncertainties on the minimum separations. These
uncertainties becomes even larger if one also accounts for
the errors on the cluster proper motions, distances, and
radial velocities. To obtain a more certain solution, we
have to wait for the measurement of the pulsar's parallax
and more accurate measurements of the proper motion,
together with observational constraints on the pulsar's
radial velocity, which could come, e.g. from the mod-
elling of a possible bow-shock PWN created in the ISM
by the supersonically moving pulsar (see, e.g. Pellizza et
al. 2005).

5. SUMMARY

Using the HST observations, we have confirmed that
the candidate proposed by M+97 is indeed the optical-
UV counterpart of PSRB1055— 52, the tenth rotation-
powered pulsar with a secured optical identification (see
Mignani 2010 for a recent review). From multi-band pho-
tometry we found that, similar to the middle-aged pul-
sars PSR B0656+14 and Geminga, its spectrum can be
described by the combination of a PLo component, with
the spectral index ao = 1.05+0.34, and a RJ component,
with brightness temperature To = (0.66+0.10) d§ 50 i?Q 2 13
MK. The PLo component is steeper than those of
PSRB0656+14 and Geminga and, possibly, of all the
other rotation-powered pulsars with known optical coun-
terparts. Moreover, it is possibly steeper than the X-ray
PL component, which might suggest that different pop-
ulations of relativistic electrons are responsible for the
X-ray and optical-UV magnetospheric emission. To ver-
ify this assumption, IR observations of PSRB1055— 52
would be particularly useful, supplemented by a joint
spectral fit of the IR-optical, X-ray, and 7-ray compo-
nents, using multiwavelength spectral models for magne-
tospheric pulsar emission.

The observed RJ component exceeds by a factor of «4
the extrapolation of the X-ray thermal component into
the UV-optical, unlike PSRB0656+14 and Geminga but
similar to RQINSs. It indicates that the temperature
distribution over the NS surface is nonuniform, and the
X-ray thermal component comes from a relatively small
hotter region, while the main contribution to the optical-
UV RJ component comes from a larger and colder area.
It implies that the distance to the pulsar is considerably
smaller than estimated from the pulsar's dispersion mea-
sure and the models for Galactic electron distribution. It
also means that the conclusions inferred from the com-
parison of the X-ray temperature and thermal luminosity
of PSRB1055— 52 with the NS cooling models require a
revision. New HST observations in the UV and near-
IR, including phase-resolved spectroscopy, are required
to separate the thermal and magnetospheric components
and infer the temperature distribution, while the mea-
surement of the radio parallax is crucial to determine
the actual distance and the NS radius.

By recalibrating the astrometry of the HST images
with the GSC2.3, 2MASS and UCAC3 catalogs, we
have measured the PSRB1055— 52 absolute coordinates
with a radial position accuracy of 0'/15 at the epoch of
2008.18, about 30 years after the measurement of the
only published pulsar's coordinates. From relative as-
trometry between the HST images taken 12 years apart,
we obtained the first measurement of the pulsar's proper
motion, /1 = 42 ± 5 mas yr _1 with a position angle of
94° ± 7°, which corresponds to the transverse velocity
V t = (70 ± 8)d 350 km s" 1 . Further HST and radio ob-
servations will allow one to measure the proper motion
with a higher accuracy, to determine the distance from
the parallax measurement, and find the true pulsar age
and the parent cluster/association.

Support for Program number 11154 was provided by
NASA through a grant from the Space Telescope Science
Institute, which is operated by the Association of Uni-
versities for Research in Astronomy, Incorporated, under
NASA contract NAS5-26555. RPM thanks A. C. Jackson

Mignani, Pavlov, & Kargaltsev

for her contribution in the reduction of the VLT data.
We thank R. N. Manchester for communicating the radio
timing position of PSRB1055— 52 before publication, A.
Dolphin for useful discussions of the corrections for CTE
losses in WFPC2, D. Vande Putte for making available
to us his code for Galactic orbit simulations, and M. van
Kerkwijk for his valuable advice on the search for the

pulsar's birthplace. We also thank the anonymous ref-
eree for the very careful reading of the manuscript and
the useful suggestions that helped us to improve the pre-
sentation of the results. The work by GGP and OK was
partially supported by NASA grant NNX09AC84G.
Facilities: HST (ACS, WFPC2),VLT (FORS1)

APPENDIX
STAR A PROPER MOTION

Since the magnitude of the Star A proper motion is crucial for determining the pulsar's absolute position and proper
motion (see Sections 3.1 and 3.3.1), we verified the UCAC3 value (/i Q = —8.9 ± 1.3 mas yr _1 ; [is = 7.8 ± 2.2 mas
yr _1 ) against other astrometric catalogs. Unfortunately, star A is too faint (V « 14.6) to be in the Hipparcos catalog
(V < 12.4; Perryman et al. 1997), while it is strangely not in the deeper Tycho-2 (V < 15.2; H0g et al. 2000) and
UCAC2 (R < 16; Zacharias et al. 2004). It is, however, in the USNO-B1.0 catalog (B < 22; Monet et al. 2003), which
reports a proper motion of fi a = +30 ± 4 mas yr _1 and fig = ± 6 mas yr _1 (much larger than in UCAC3), and
in the PPMXL catalog (Roeser et al. 2010), which gives fi a = —13.8 ± 8.2 mas yr _1 and us = +4.5 ± 8.2 mas yr _1 ,
consistent with the UCAC3 value (albeit less accurate).

As an independent test, we measured the relative proper motion of Star A. We used archival Very Large Telescope
(VLT) images of the PSRB1055-52 field taken with FOcal Reducer/low dispersion Spectrograph ( FORS1) on 2002
March 21, which yields a time baseline of about 6 years with respect to the WFPC2 observations. The images were
taken in high resolution mode (O'.'l/pixel) but with a windowing of the CCD (1.'6 x lf6 FOV). A sequence of eighteen
140 s dithered exposures was obtained through the Bcsscl B filter (A = 4290A; A A = 880 A). We applied standard CCD
reduction with the MIDAS package and performed the astrometric calibration as described in Section 3.3 (position
error S r = 0'.'3). Unfortunately, the image quality of the FORS1 data (seeing ~ l'/2) does not allow one to detect the
pulsar in the halo of Star A.

We then compared the relative position of Star A measured in the FORS1 (2002.21) and WFPC2 (2008.18) images.
We used the F450W WFPC2 image, where Star A is not saturated and the pivot wavelength is similar to that of
the Bessel B filter. We registered the WFPC2 image onto the FORS1 image, aligned along right ascension and
declination to better than ~ 0.05°, by computing the transformation between the pixel coordinates of 15 stars seen
in both images. This yielded an rms of 0'.'015 along each of the two axes, larger than the object centroid uncertainty
(0'.'003). We thus determined the offset of Star A to be Aa cos S = —2.2 ± 15 mas and AS = 24.6 ± 15 mas, which
corresponds to the 3cr upper limit on the proper motion of ~ 15 mas yr _1 , consistent with both the UCAC3 and
PPMXL values.

REFERENCES

Abdo, A. A., ct al. 2010, ApJS, 187, 460

Anderson, S. B., Cordova, F. A., Pavlov, G. C, Robinson, C. R.,

& Thompson, R. J., Jr. 1993, ApJ, 414, 867
Apcllancz, J. M., & Cox, C. 2008, A New Geometric Distortion

Solution for the ACS/SBC, Instr. Science Rep. ACS 2008-2,

Space Telescope Science Institute
Becker, W., & Triimper, J. 1997, A&A, 326, 682
Bignami, G. F., Caraveo, P. A., & Vacanti, G. 1988, A&A, 196,

191
Brisken, W. F., Thorsett, S. E., Golden, A., & Goss, W. M. 2003,

ApJ, 593, L98
Brinkmann, W. & Ogelman, H. 1987, A&A, 182, 71
Caraveo, P. A., Lattanzi, M. G., Massone, G., Mignani, R. P.,

Makarov, V. V., Perryman, M. A. C, Bignami, G. F. 1998,

A&A, 329, LI
Caraveo, P. A., De Luca, A., Mignani, R. P. & Bignami, G. F.

2001, ApJ, 561, 930
Cheng, A. F., & Hclfand, D. J. 1983, ApJ, 271, 271
Cordes, J. M., & Lazio, T. J. W. 2002, preprint

(astro-ph/0207156)
De Luca, A., Caraveo, P. A., Mereghetti, S., Negroni, M.,

Bignami, G. F. 2005, ApJ, 623, 1051
de Zeeuw, P. T., Hoogerwerf, R., de Bruijne, J. H. J., Brown, A.

C. A., & Blaauw, A. 1999, AJ, 117, 354
Dias, W. S., Alcssi, B. S., Moitinho, A., & Lepine, J. R. D. 2002,

A&A, 389, 871
Dolphin, A. E., 2009 PASP, 121, 655
Fierro, J. M., et al. 1993, ApJ, 413, L27
Fitzpatrick, E. L., 1999, PASP, 111, 63
Greiveldinger, C, et al. 1996, ApJ, 465, L35
Giullcmot, L. 2009, in Fermi Symposium 2009, on-line talks

Lyne, A. G., & Kramer, M. 2005,

Hobbs, G., Lorimer, D. R.

MNRAS, 360, 974
Hog, E., et al. 2000, A&A, 355, L27
Kaminkcr, A. D., Yakovlev, D. G., & Gnedin, O. Y. 2002, A&A,

383, 1076
Kaplan, D. L. 2008, Proc. of 40 Years of Pulsars: Millisecond

Pulsars, Magnctars and More, AIP Conf. Proc, 983, 331
Kargaltsev, O., Pavlov, G. G., & Romani, R. W. 2004, ApJ, 602.

327
Kargaltsev, O., Pavlov, G. C, Zavlin, V. E., & Romani, R. W.

2005, ApJ, 625, 307
Kargaltsev, O. & Pavlov, G. G. 2007, Ap&SS, 308, 287 (KP07)
Lasker, B. M., Sturch, C. R., McLean, B. J., Russel, J. L.,

Jcnkner, H., & Shara, M. M. 1990, AJ, 99, 2019
Lasker, B. M., ct al. 2008, AJ, 136, 735
Lattanzi, M. C, Capetti, A., & Macchctto, F. D. 1997, A&A,

318, 997
Manchester, R. N., et al. 1978, MNRAS, 184, 159
Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M. 2005,

AJ, 129, 1993
Mignani, R., Caraveo, P. A., & Bignami, G. F. 1997, ApJ, 474,

L51 (M+97)
Mignani, R. P., Zharikov, S.

891
Mignani, R. P., Sartori, A., De Luca, A., Rudack, B.,

Slowikowska, A., Kanbach, G., Caraveo, P. A. 2010, A&A, 515,

110
Mignani, R. P. 2010, ASpR, in press (arXiv:0912.2931)
Mineo, T., Massaro, E., Cusumano, G., & Becker, W. 2002, A&A,

392, 181
Monet, D. G., et al. 2003, AJ, 125, 984

Caraveo, P. A. 2007, A&A, 473,

HST observations of PSRB1055-52

Newton, L. M., Manchester, R. N., & Cooke, D. J. 1981,

MNRAS, 194, 841
Nota, A., et al., 1996, FOC Instrument Handbook Ver. 7.0

(Baltimore: STScI)
Ogelman, H., & Finley, J. P. 1993, ApJ, 413, L31
Pavlov, G. G., Shibanov, Yu. A., Zavlin, V. E., & Meyer, R. D.

1995, in The Lives of the Neutron Stars, eds. M. A. Alpar, U.

Kiziloglu, & J. van Paradijs (Kluwcr: Dordrecht), p. 71
Pavlov, G. G., Zavlin, V. E., Triimper, J., & Neuhauser, R. 1996,

ApJ, 472, L33
Pavlov, G. G., Welty, A. D., & Cordova, F. A. 1997, ApJ, 489,

L75
Pavlov, G. G., Zavlin, V. E., & Sanwal, D. 2002, in Neutron Stars,

Pulsars and Supernova Remnants, eds. W. Becker, H. Lesch, &

J. Triimper, MPE-Report 278, 273 (astro-ph/0206024)
Pellizza, L. J., Mignani, R. P., Grenier, I. A., & Mirabel, I. F.

2005, A&A, 435, 625
Perryman M. A. C, ct al. 1997, A&A, 323, L49
Predehl, P., & Schmitt, J. H. M. M. 1995, A&A, 293, 889

Roeser, S., Demleitner, M., & Schilbach, E. 2010, AJ, 139, 2440

Seaton, M.J. 1979, MNRAS, 187, 73

Shibanov, Yu. A., Zavlin, V. E., Pavlov, G. G., & Ventura, J.

1993, A&A, 266, 313
Skrutskie, M. F., et al. 2006, AJ, 131, 1163
Taylor, J. H., & Cordes, J. M. 1993, ApJ, 411, 674
Thompson, D. J., et al. 1999, ApJ, 516, 297
Vande Putte, D. W., Gamier, T. P., Ferreras, I., Mignani, R. P.,

&, Cropper, M. S. 2010, MNRAS, in press
Vaughan, A. E. & Large, M. L. 1972, MNRAS, 156, 27P
Weltevrede, P., et al., 2010, PASA, 27, 64
Yakovlev, D. G., & Pethick, C. J. 2004, ARA&A, 42,169
Zacharias, N., ct al. 2004, AJ,127, 3043
Zacharias, N., et al. 2010, AJ, 139, 2184
Zavlin, V. E., & Pavlov, G. G. 2004, ApJ, 616, 452
Zharikov, S., Shibanov, Yu., Komarova, V. 2006, AdSpR, 37, 1979

10

Mignani, Pavlov, & Kargaltscv

TABLE 1

HST PHOTOMETRY OF THE PSRB1055 — 52 COUNTERPART.

Filter

log V

t

r

C r "

m r

Am r

AniCTE

wiob B

pobs

pdor

[Hz]

s

pix

counts

A*Jy

A*Jy

WFPC2 F450W

14.818 (0.045)

600

3

<13.7

>25.61

-0.27

-0.37

>24.97

< 0.26

< 0.34

WFPC2 F555W

14.741 (0.049)

600

3

27.3 ±5.0

25.90 ±0.20

-0.24

-0.23

25.43

0.241 ± 0.044

0.295

WFPC2 F702W

14.637 (0.043)

600

3

30.3 ±3.5

26.46 ±0.13

-0.26

-0.13

26.08

0.214 ± 0.024

0.249

SBC F140LP

15.293 (0.042)

5569

6

525 ± 24

22.83 ± 0.05

-0.19

22.64

0.248 ±0.011

0.417

FOC F342W

14.945 (0.033)

2916

8

242 ± 13

25.18 ±0.06

-0.18

25.00

0.140 ± 0.008

0.190

Note. — Columns report the camera and the filter, the pivot frequency and the band width (in parentheses) in logarithmic units (as derived
from the SYNPHOT tables), the integration time of the average-combined frames, the photometry aperture radius r in detector pixels, the total
number of counts (C r ) measured within the aperture r, the corresponding magnitude in the STmag system (m r ), the aperture correction (Am r ),
the CTE correction (AmcTE), the observed magnitude m b s , and the observed (F° fa ) and extinction-corrected (F^ CI ) spectral fluxes.

a For the WFPC2 and FOC observations, counts have been measured on the average-combined frames,
average of the first and second exposure (see § 3.2).

For the FOC, we used only the

Fig. 1.— SBC F140LP (FUV band; left panel) and WFPC2 F702W (red band; right panel) images of the PSRB1055-52 field (10" X 10").
North is to the top, East to the left. The two objects in the SBC image are the pulsar (in the smaller thin circle of 0"3 radius, corresponding
to a nominal uncertainty of SBC position) and Star A (4'.'4 northwest of the pulsar). Zoomed images of the pulsar and Star A, demonstrating
the count distributions in and around these sources, are shown in the insets; the radii of the circles in the insets are 0'.'3 and 0"7. The
small thin circle in the WFPC2 panel, labeled "SBC" is at the same position with respect to Star A as in the left (SBC) panel. In each of
the panels the smaller thick circle of 0"45 radius is the error circle around the radio pulsar position at the epoch of 1978.13 (Newton et al.
1981), while the larger thick circle of 3" 4 radius encloses possible positions of the pulsar at the epoch of our HST observations if the pulsar
transverse velocity (in an unknown direction) were 400 km s — 1 , at d = 750 pc. The thin circle labelled "FOC" corresponds to the position
of the PSRB1055— 52 candidate counterpart of M+97; its radius of l'.'O corresponds to the uncertainty of a nominal FOC position. See
§3.1 for more detail. Notice the much higher relative brightness of the pulsar in the SBC image.

HST observations of PSRB1055-52

11

0.5

0.0 -

-0.5 -

o

-1.0 -

-1.5 -

-j— i— i— i— j— i—
E(B-V) = 0.07

14.2 14.4 14.6 14.8 15.0

Log v (Hz)

15.2

Fig. 2. — Optical-UV spectrum of PSRB1055— 52. The spectral fluxes have been corrected for interstellar extinction assuming E(B — V) =
0.07. The dashed lines correspond to the Rayleigh- Jeans (RJ; red) and power-law (PLo; blue) components of the best fit, the solid line
shows the sum of these components. The shaded areas show lcr uncertainties of the fit.

12

Mignani, Pavlov, & Kargaltsev

' « ' ' l—l—l I— I I— I I— I I— I I— I I ■ ' ■ ' I— I I— I I— I I— I I— I I I I I !■■ I I

15

16
Log v (Hz)

17

18

^>

~i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r~

BB, ..•••*'

Chandra
XMM-Newton

RJ<

PL

-I I I I I I l_ I I I I I I I I—

—I I I—

16

18 20

Log v (Hz)

22

24

Fig. 3. — Left: Comparison of the optical-UV and X-ray spectral models for PSRB1055— 52. The X-ray fit for the BBc + BB H +
PLx model is taken from De Luca et al. (2005). The excess of the RJ component over the extrapolation of the X-ray BBc component
into the optical-UV range could be explained by the presence of a "very cold" thermal component, BBvo The dotted blue curve shows
this component for the maximum allowed temperature Tvc = 0.45 MK. The shaded areas show la incertainties of the fits. Right:
Multiwavelength SED for PSRB1055-52, from optical to 7-rays, with spectral models for separate ranges. The 7-ray points are from
Guillemot (2009), the spectral model from Abdo et al. (2010). The line v~ °- 46 shows a PL spectrum approximately connecting the optical
and 7-ray bands. See Section 4.1 for more details.