A&A 480, 623-628 (2008)
DOI: 10.1051/0004-6361:20066806
O. Löhmer1 - A. Jessner1 - M. Kramer2 - R. Wielebinski1 - O. Maron3,1
1 - Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany
2 -
University of Manchester, Jodrell Bank Centre for Astrophysics,
Manchester M13 9PL, UK
3 - J. Kepler Institute of Astronomy, University of Zielona Gora, Poland
Received 24 November 2006 / Accepted 11 December 2007
Abstract
Aims. The behaviour of the pulsar spectrum at high radio frequencies can provide decisive information about the nature of the radio emission mechanism.
Methods. We report recent observations of a selected sample of pulsars at mm (32 GHz) with the 100-m Effelsberg radio telescope.
Results. Three pulsars, PSR B0144+59, PSR B0823+26, and PSR B2022+50, were detected for the first time at this frequency. We confirm the earlier flux density measurements for a sample of six pulsars, and we are able to place upper flux density limits for another 12 pulsars. We find that all pulsar spectra have a simple form that can be described using only three parameters, one of which is the lifetime of short nano-pulses in the emission region. The study of the transition region from coherent to incoherent emission needs further and more sensitive observations at even higher radio frequencies.
Key words: stars: pulsars: individual: PSR B0144+59 - stars: pulsars: individual: PSR B0823+26 - stars: pulsars: individual: - PSR B2022+50 - radiation mechanisms: non-thermal - radio continuum: stars
Almost 40 years after the discovery of pulsars, it is still not possible to give a detailed explanation of the radio emission process, and a large number of fundamental questions remain (e.g. Lorimer & Kramer 2005) to be answered. Advances in our theoretical understanding of the radio emission and the associated conditions in the pulsar magnetosphere (e.g. Petrova 2001; Melrose et al. 2006) are still largely driven by observational progress (e.g. Kramer et al. 2006).
In recent years, most observational studies pertaining to emission processes have concentrated on polarisation properties (e.g. Karastergiou & Johnston 2004) or on the drift behaviour of sub-pulses (e.g. Weltevrede et al. 2006). However, as past studies of other discrete radio sources have demonstrated, the investigation of radio frequency spectra can be crucial for identifying of the underlying emission process (e.g. Conway et al. 1963). Although we may expect that similar detailed studies of pulsar spectra can be helpful in solving the ``pulsar problem'', such studies (e.g. Malofeev et al. 1994) have so far not been successful in establishing a relationship between pulsar spectra and other observed properties. It turned out that significant correlations between spectral properties and pulsar spin parameters could not be found or did not survive more detailed scrutiny (Maron et al. 2000). It is therefore more likely that spectral studies of separate profile components (Kramer 1994) or even of individual pulses (Kramer et al. 2003; Karastergiou & Johnston 2006) are better suited to addressing questions about the emission process, as the average pulse profile seen at different frequencies is dominated by different components with different spectral behaviours. Differences in the flux density spectra of separate components can indeed explain certain features in observed spectra of the average profiles, such as breaks in power laws, but may also contribute to the observed polarisation evolution with frequency (Karastergiou et al. 2005). Other authors have, however, claimed that the general decrease in polarised emission at high frequencies is due to a frequency-dependent decrease in the coherence of pulsar radio emission (e.g. Xilouris et al. 1996).
The steepness of pulsar spectra (mean spectral index of
with indices between
and
for power
law spectra,
,
Maron et al. 2000)
may be caused
by the mechanisms that enforce the coherence of pulsar radiation as has
been suggested by Michel (1991).
This concept is based on a model in which an initially incoherent radiation process
produces a coherent output by means of antenna mechanisms or resonances.
The radiation produced by other considered processes, especially
maser mechanisms and plasma radiation, is however inherently coherent.
Should the first class of models be applicable, then pulsar
radiation would be a combination of coherent and incoherent emission, where
the coherent emission is seen to dominate at radio frequencies,
establishing a typical coherence length
beyond which
incoherent emission should dominate. Interestingly, the non-thermal
emission of pulsars observed at infrared and optical frequencies is
indeed believed to be due to incoherent synchrotron emission
(Crusius-Waetzel & Lesch 2002).
The study of pulsars at very high radio frequencies, where the coherent emission process might become less dominant, is an obvious procedure for investigating a possible intrinsic coherence length and the contribution of an additional incoherent emission component. So far, the flux density spectrum of pulsar emission in the transition region between the coherent radio emission and the incoherent infrared emission is unknown. We therefore aim to investigate the exact position and nature of this transition from coherent to incoherent emission through access to the spectral region where the spectral index is flattening or perhaps showing an upturn.
The information gap between high radio and infrared frequencies was first pointed out after early observations of the Crab pulsar (e.g. Smith 1977) at radio and optical wavelengths. The Crab pulsar, among a few others, also provides direct evidence about the fine details of the radio-emission process in the form of giant radio pulses (GRPs). These have been observed up to 15 GHz for the Crab pulsar (Hankins private communication, Kondratiev et al., in preparation) and occur during all rotation phases where the average pulse profile shows radio emission (Jessner et al. 2005). It is worthwhile investigating the possibility of the spectrum, or parts of it, being the natural result of a superposition of a very large number of very short pulses, similar to the ``nano-shots'' (Eilek & Hankins 2006) observed for the Crab pulsar. One can show that a simple form of the spectrum will be the result, depending - apart from a scale factor - on at most two more parameters, one of which is the typical decay time of the nano-shots.
The filling of the spectral gap, hence the attempt to learn more about the radio emission process by studying this crucial frequency range, was one of the motivations for the technically challenging observations of pulsars at mm-wavelengths (Wielebinski et al. 1993; Kramer et al. 1996). Only at high frequencies will propagation effects become negligible for the majority of pulsars, their intrinsic emission properties may reveal themselves only here. This is commonly appreciated for the interstellar medium, but it may even be important for propagation effects in the pulsar magnetosphere (e.g. Hoensbroech & Lesch 1999).
As we summarise in the following, the first observations of pulsars at mm-wavelengths provided some unexpected results, which triggered further follow-up studies. The results of such a study are described in this paper, confirming the results for the known detected pulsars but with an additional enlargement of the numbers of observed pulsars to provide better statistics.
In this paper we briefly review the main properties of pulsar spectra, before we describe new observations of pulsars at mm-wavelengths and present their results. We then compare these new results with previous measurements before we discuss the implications of our results and our conclusions are drawn.
Prior to the observations by Wielebinski et al. (1993),
it was commonly believed that the generally steep pulsar spectrum
continues to the highest radio frequencies (e.g. Manchester & Taylor
1977; Malofeev et al. 1994). On the other
hand, one can even
find claims in the literature that pulsar spectra show a ``cut-off''
at about 10 GHz or so (e.g. Gurevich et al. 1993). This
erroneous impression was caused by the lack of large sensitive telescopes
and high-frequency receivers at that time.
The absence of any cut-off was demonstrated by the first
detection of four pulsars at 33.8 GHz (
mm Wielebinski et al. 1993).
In fact, two
pulsars were stronger at these very high radio frequencies than was
previously thought. In contrast to predictions based on the
extrapolation of the spectra measured up to
24 GHz, the
measured flux densities of these two sources (PSRs B1929+10 and
B2021+51) were factors of five to ten higher than expected. This
apparent spectral turn-up was further investigated by Kramer et al. (1996). With their observations between
27 GHz and
35 GHz, they confirmed the results for B1929+10 and B2021+51,
but also established a normal continuation of the known spectrum for six
additional sources. The four pulsars strongest at 32 GHz were
later also successfully observed at 43 GHz (
mm),
underlining the impression that pulsar spectra show
unusual behaviour at high radio frequencies; i.e. they may become flatter
or even possibly show a turn-up at mm-wavelengths (Kramer et al. 1997). The first and so far only
detection
of a radio pulsar at 87 GHz (
mm) with the 30-m Pico
Veleta telescope (Morris et al. 1997) shows a
continuation of the relatively flat spectrum of the detected pulsar
B0355+54.
However, recent reports suggest that turn-ups in pulsar
spectra may even occur at much lower radio frequencies (Kijak & Maron 2004).
Clearly, more studies of more pulsars
at mm-wavelengths are required.
We selected a sample of 15 pulsars from the catalogue of Maron et al. (2000), from the list of brightest pulsars, and from the list of pulsars with unusual frequency evolution given by von Hoensbroech & Lesch (1999). All observations were made with the 100-m Effelsberg radio telescope of the Max-Planck-Institut für Radioastronomie in Bonn, Germany. We observed the sources at three separate epochs with good to excellent weather conditions: 2002 March, 2002 September, and 2003 January. Integration times of the pulsars were between 30 and 180 min, depending on the expected flux density and visibility of the sources.
Table 1: Results of observations at 32 GHz.
We used an HEMT receiver installed in the secondary focus with a fixed
central frequency of 32 GHz (
mm). The total bandwidth of
the system was 2 GHz, with a receiver noise temperature of about
60 K. The telescope gain at this frequency is
1 K/Jy.
The receiver provided LHC and RHC signals that were digitised
and independently sampled every P/1024 s and synchronously folded
with the topocentric pulse period P, using the Effelsberg Pulsar
Observation System (EPOS, Jessner 1996). Records of 15 s
sub-integrations were transferred to disk for off-line data analysis.
For flux-calibration purposes, we made use of an internal noise diode fed directly into the waveguide following the feed horn. EPOS allows us to switch this noise diode synchronously with the pulse period, occupying the first 50 phase bins of the pulse period in the measured profile. To check that the actual pulsar signal, expected to be relatively weak at 32 GHz, does not overlap with the strong, artificial calibration signal, we also observed each pulsar briefly at 4.85 GHz and 8.35 GHz, making sure that it is positioned correctly, and then switched to the high frequency receiver in less than 30 s that way preserving the pulse phase for all frequencies. The strong and clearly visible pulsar signals observed at the two lower frequencies could then be directly related to the possible signal of the high frequency emission, greatly facilitating confirmation of a possible detection.
Flux densities were obtained for all three frequencies by comparing the
pulse profiles to the corresponding signal of the noise diode, which
itself was calibrated by observing known reference sources during
regular pointing observations. The continuum references
(NGC 7027, 3C 286, and NRAO 150) were observed approximately
every hour, checking both pointing and focus stability. Typical
pointing errors were rms
5'', compared to a beamwidth of 23''at 32 GHz. Flux densities for the chosen reference sources were
obtained from the catalogue of Peng et al. (2000).
Uncertainties in the resulting pulsar flux densities are estimated
to be about 20% for a single observation.
![]() |
Figure 1: Time-aligned profiles of the three PSRs B0144+59, B0823+26 and B2022+50 at centre frequencies of 4.85, 8.35 and 32 GHz. |
Open with DEXTER |
![]() |
Figure 2: Resulting flux density spectra for PSRs B0144+59, B0823+26, and B2022+50. Low-frequency data have been taken from available literature (see text). The dashed line shows the model fit with parameters listed in Table 2. |
Open with DEXTER |
In total, we observed a sample of 21 pulsars. The results of the measurements are summarised in Table 1. We were able to detect PSRs B0144+59, B0823+26, and B2022+50 at 32 GHz for the first time, determining their flux densities at their so far highest observation frequency. The time-aligned profiles of these sources measured at 4.85, 8.35, and 32 GHz are shown in Fig. 1, while their resulting spectra are displayed in Fig. 2. Flux densities shown for lower frequencies are taken from the literature (Maron et al. 2000) or from unpublished observations.
We also successfully observed all those pulsars that were already detected by Wielebinski et al. (1993), Kramer (1995), and Kramer et al. (1997). The measured flux densities are in very good agreement with previous observations.
For those pulsars that were not detected, we estimate upper flux limits
following the procedure described by Sieber, & Wielebinski (1987).
The upper limit is based on an estimate of five times the RMS of the noise
signal multiplied by the equivalent pulse width. For our sample we used the
pulse width W50 at the highest known frequency as given in the EPN pulsar
data base or by Kramer et al. (1994) or Kijak et al. (1998).
Our results confirm all previously published measurements of flux densities of pulsars at mm-wavelengths. We also present spectral information for three newly detected pulsars at mm-wavelengths. No new turn-up in the spectrum was found. The spectra of the newly detected pulsars follow the trends determined from observations at lower radio frequencies. While on one hand this gives us great confidence in the reliability of the adopted calibration procedure, it does suggest, on the other, that a spectral turn-up at a frequency as low as 30 GHz is the exception rather than the rule. We note that the spectrum of PSR B0144+59 shows a peculiar ``kink'' at 3 GHz to 10 GHz, which may be misinterpreted as a ``turn-up'' or definite ``flattening'' unless data above 10 GHz are considered (cf. Kijak & Maron 2004). We believe that this is caused by the peculiar high-frequency profile evolution of this source, which will be studied elsewhere. In any case, it adds to the notion that some pulsar spectra may not be adequately described by simple power laws. Moreover, in view of the evidence that individual pulsar spectra are quite different (differing not only in shape but also e.g. in the location of the low-frequency turn-over, spectral index, and the possible existence of a spectral steepening (Maron et al. 2000), it would be rather naive to assume that all pulsars show identical behaviour in the narrow frequency range probed by our mm-observations. Hence we should not expect to detect a spectral turn-up for all pulsars at 32 GHz. Instead, the known diversity of pulsar properties may indeed eventually reveal a turn-up or flattening of pulsar spectra at higher or even at lower frequencies (Kijak & Maron 2004).
Although pulsar spectra are traditionally described by a power law
with
about -1.6 (Sieber 1973),
many of the observed spectra turned out to fit badly to that concept.
A number of additional characteristic features are traditionally employed to describe
more complex pulsar spectra : some of them
exhibit a turn-over at low
,
some have a ``broken power law'', some a
turn-up or flattening, etc. The original rationale was the analogy of
synchrotron radiation, which has a power-law spectrum, so that radio astronomers
initially tried to model most observations in terms of power laws. But one of
the implicit assumptions was that one wanted to model a temporarily continuous
and spatially large-scale radiation mechanism. The observations have, however, come up
with strong variability on timescales from days down to nanoseconds (GRPs).
A few attempts have been made to fit another function to a pulsar
spectrum, most notable Ochelkov & Usov (1984),
who proposed a six-parameter model equation:
![]() |
(1) |
Do we know of other natural processes, involving flowing charges, that produce spectra of the observed form? Indeed we do, solid-state currents in everyday materials exhibit a non-thermal low frequency noise, often called ``flicker noise'', ``1/f-noise'' or simply lf-noise. The spectrum is very similar to a pulsar spectrum, but at frequencies lower by a factor of 107-108. That in itself is not surprising, as the conditions on the pulsar surface and in the magnetosphere are also many orders of magnitude different from our experience with everyday materials. People had already been studying that kind of noise in solids and on surface boundaries in the 1950's (Pfeifer 1959; Bess 1953). It is indeed surprising how far one can get using this very simple one-parameter theory.
We already know from observations that pulsar radio emission varies on very short
timescales. Very high-resolution observations of giant radio
pulses have so far not found a shortest timescale for the emission, although
they have reached nanosecond to sub nanosecond resolution (Hankins et al. 2003;
Jessner et al. 2005; Eilek & Hankins 2006)
for the Crab pulsar. Eilek & Hankins (2006)
report an apparent increase in micro-burst duration
,
nano-shots with
ns, and they infer that each
giant pulse is a coherent superposition of a number of powerful nano-pulses.
Based on these theoretical studies and on the observational evidence, we assume
that all pulsar radio emission is a superposition of short-lived
(
)
elementary emission processes. The observed shape of the
spectrum may then be used to provide constraints on the timescales of the elementary
processes.
A shot noise model for the emission due to short-lived non-Gaussian elementary processes in ordinary pulsars has been proposed by Jenet et al. (2001). However, in this study we use a different elementary process, one that has been commonly used to model flicker noise in solids and vacuum discharges, and evaluate the expected observable broad band radio spectrum.
![]() |
Figure 3:
Spectral fits for PSR B2021+51 using additional data from O. Maron et al.
Solid: fit with a spectrum according to Eq. (3) with S0=1.61 Jy,
![]() ![]() ![]() |
Open with DEXTER |
Let us assume that we
have a stochastic formation of radio emission centres that are characterised
by a localised potential
with a lifetime (observer system) of i.e.
s, very similar to the nano pulses proposed i.e. by
Weatherall (1998) and also evident in the
PIC-simulations of Jaroschek & Lesch (2006). The decay of the electric
potential of a nano-pulse may be described by
with (t>0). A random superposition of short pulses has the simple spectrum
![]() |
(2) |
![]() |
(3) |
Table 2: Parameters of spectral fits for a number of pulsars.
We have presented the first detection of three more pulsars at mm-wavelengths, while presenting upper limits of flux densities for 12 other pulsars. None of these upper limits is yet low enough to establish the nature of the flux density spectrum of those sources, the continuation of the previous spectral trends, or a possible flattening or turn-up at high frequencies. Indeed, while the three newly detected pulsars show a continuation of the lower frequency spectra, we confirmed the flux densities of six previously detected pulsars. This sample of re-detected pulsars includes PSRs B1929+10 and B2021+51, thus giving further credibility to their published flux density measurements, which led to the suggestion of a spectral turn-up at mm-wavelengths for these pulsars. Establishing such a trend could provide a possible link between the spectrum of the coherent radio emission and the incoherent higher frequency radiation, but more sensitive observations at even higher frequencies are needed to substantiate any conclusion in this direction.
We have also found evidence that a heuristic model of pulsar radio emission,
consisting of a superposition of a large number of short pulses of only
nano-second duration describes observed pulsar radio spectra
quite naturally. Apart from intensity scaling, a model fit needs only two parameters,
one being the typical nano pulse life time
,
the other (n)
related to the geometry of the individual emission process. These processes
already play a role, with much greater magnitude for the individual nano-pulse,
in known giant pulse emission of several pulsars. Different emission
heights for different frequencies would have to be due to propagation effects (refraction)
in the pulsar magnetosphere. A detailed appraisal of these issues is in preparation.
To achieve final reliable proof of these findings, it is important to
obtain even higher radio frequency data of sufficient sensitivity. Future observations
of single pulses at very high frequencies combined with very high sensitivity should
reveal the elementary emission processes we propose here in ordinary
pulsars. Our present observations extend the spectra of observed pulsars into the long mm
wavelength range.
So far only the detection of one object, PSR 0355+54, has been made at 87 GHz with
the Pico Veleta 30-m telescope (Morris et al. 1997) at the
level. Observations of nine pulsars at
mm (345 GHz) have not been
successful (Löhmer et al. 2004). The advent of ALMA should change this.
ALMA (www.eso.org/projects/alma/) will have observational capabilities in
ten frequency bands. Initially the frequency bands 84-116 GHz, 275-373 GHz, and
602-720 GHz will be available. An IF bandwidth (in 2 polarisations) of 8 GHz
will be provided. The flux sensitivity of ALMA is expected to be sub-mJy (about 0.1 mJy) range
after 10 min of integration. Considering a 4-h integration, a detection level
of roughly 0.05 mJy should be reached. This would allow the study of the flux densities of numerous pulsars even if they do not show a turn-up in their spectra. By these observations the gap between radio and optical (IR) flux detections could be closed, giving us new information about
the emission mechanism inherent to pulsars.
For a proper assessment of the range of coherent emission from pulsars one needs
more information about the spectral
behaviour of pulsars at frequencies below 100 MHz where observational data is
comparatively rare. Here pulsars are
usually fairly strong radio sources with flux densities in the range of several Jy.
The new generation of large low-frequency interferometers (LOFAR and possibly LWA)
will be able to make detailed spectral measurements over several octaves on metre and dekametric wavelengths (Eilek et al. 2006; Stappers et al. 2006). The knowledge about the upper
and lower limits of coherent pulsar radio emission is expected
to become a decisive constraint on pulsar radio emission theories.
Acknowledgements
The authors would like to thank the staff of the radio-observatory Effelsberg for their help and support during the observations. We would also like to thank Harald Lesch, Jan Kuijpers and Francis Graham-Smith for encouraging discussions about the shape of pulsar spectra. We are grateful to the anonymous referee for helpful advice.