Issue |
A&A
Volume 499, Number 2, May IV 2009
|
|
---|---|---|
Page(s) | L17 - L20 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/200912327 | |
Published online | 29 April 2009 |
LETTER TO THE EDITOR
Discovery of a bright radio transient in M 82: a new radio supernova?
A. Brunthaler1 - K. M. Menten1 - M. J. Reid2 - C. Henkel1 - G. C. Bower3 - H. Falcke4,5
1 - Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69,
53121 Bonn, Germany
2 -
Harvard-Smithsonian Center for Astrophysics, 60 Garden Street,
Cambridge, MA 02138, USA
3 -
UC Berkeley, 601 Campbell Hall, Astronomy Department & Radio
Astronomy Lab, Berkeley, CA 94720, USA
4 -
Department of Astrophysics, Radboud Universiteit
Nijmegen, Postbus 9010, 6500 GL Nijmegen, The Netherlands
5 -
ASTRON, Postbus 2, 7990 AA Dwingeloo, The Netherlands
Received 14 April 2009 / Accepted 21 April 2009
Abstract
In this Letter, we report the discovery of a new bright radio transient in
82. Using the Very Large Array, we observed the nuclear region of M 82 at
several epochs at 22 GHz and detected a new bright radio source in this
galaxy's central region. We find a flux density for this flaring source
that is
times larger than the upper limits determined in
previous observations. The flare must have
started between 2007 October 29 and 2008 March 24. Over the past year,
the flux density of this new source has decreased from
100 mJy to
11 mJy. The lightcurve (based on only three data points) can be
fitted better with an exponential decay than with a power law. Based on
the current data we cannot identify the nature of this transient
source. However, a new radio supernova seems to be the most natural
explanation. With its flux density of more than 100 mJy, it is
at least 1.5 times brighter than SN1993J in M 81 at the peak of its
lightcurve at 22 GHz.
Key words: stars: supernovae: general - radio continuum: general - galaxies: individual: M 82
1 Introduction
M 82 is a nearby (3.6 Mpc based on a Cepheid distance to M 81 by Freedman et al. 1994) irregular (I0) galaxy with a very active starburst in its nuclear region. It harbors many bright supernova remnants in its central region, which have been studied extensively for decades (Fenech et al. 2008; Beswick et al. 2006; Muxlow et al. 1994). van Buren & Greenhouse (1994) estimate a supernova rate of 0.1 year-1. However, a new radio supernova has not been discovered. Singer et al. (2004) reported a supernova in M 82 (SN2004am) that was classified as type-II (Mattila et al. 2004), but it has not been detected at radio wavelengths (Beswick et al. 2004).
Kronberg & Sramek (1985) and Kronberg et al. (2000) monitored
the flux densities of 24 radio sources in M 82 from 1980 until 1992. Most
sources (75%) remained surprisingly constant. There is some controversy about
how the fluxes of these compact radio sources can be stable.
Models of supernova remnants expanding into a dense medium may explain this
(Chevalier & Fransson 2001). Seaquist & Stankovic (2007) argue that the
radio emission could arise from wind-driven bubbles. Studying the evolution
of a young source could be very important for understanding these models.
The strongest source, 41.9+58.0, shows an exponential decay with a decay
rate of
years (Kronberg et al. 2000). Another source,
41.5+597, which was detected in 1981 with a flux density of
10 mJy,
faded within a few months to a flux density below 1 mJy
(Kronberg et al. 2000; Kronberg & Sramek 1985). The nature of this
strongly variable source was never clarified.
Radio supernova are rare events. So far only about two dozen have been detected (Weiler et al. 2002) and most of them were quite distant and rather weak. This makes it difficult to study them in great detail. One notable exception is SN1993J (Schmidt et al. 1993) in M 81, which has been studied extensively (Bartel et al. 2007; Bietenholz et al. 2001; Marcaide et al. 1997; Bartel et al. 2002; Bietenholz et al. 2003; Pérez-Torres et al. 2002; Marcaide et al. 2009; Pérez-Torres et al. 2001). Thus, the detection of a new nearby supernova would be highly desirable.
M 82 is part of the M 81 group of galaxies and it shows clear signs of tidal
interaction with M 81 and NGC 3077 (Yun et al. 1994). This makes the M 81
group an ideal system for studying galaxy interaction in great detail. With
this
motivation, we have initiated a project to measure the proper motions of M 81
and M 82 with VLBI astrometry. We are observing M 81*, the nuclear radio source
in M 81, bright water masers in M 82, and three compact background
quasars. Based on our experience with measurements of proper motions in the
Local Group (Brunthaler et al. 2007,2005b), we
expect a detection of the tangential motions of M 81 and M 82 relative to the
Milky Way within a few years. So far, we have observed at three epochs at 22
GHz with the High Sensitivity Array (including the Very Long Baseline Array,
the Very Large Array, and the Greenbank and Effelsberg telescopes).
Here, we report the detection of a new transient source in M 82 based on the
data from the NRAO Very
Large Array (VLA).
2 Observations and data reduction
M 82 was observed with the Very Large Array (VLA) as part
of the High Sensitivity Array observation under projects BB229
and BB255 on 2007 January 28, 2008 May 03, and 2009 April 08. The total observing time in each epoch was 12 h. We
used M 81* as phase calibrator and switched between M 81*, M 82, and 3
extragalactic background quasars every 50 s
in the cycle M 81* - 0945+6924 - M 81* - 0948+6848 - M 81*
- M 82 - 1004+6936 - M 81*, yielding an integration time of
100 min on M 82. We observed with two frequency
bands of 50 MHz, each in dual circular polarization.
The data reduction was performed with the Astronomical
Image Processing System (AIPS) and involved standard steps.
On 2007 January 28 and 2009 April 08, we used 3 C48 as
flux density calibrator. M 81* was used as a gain and phase calibrator and
one round of phase self-calibration was performed on M 82.
Unfortunately, no flux density calibrator was observed on 2008
May 03. Here, we assumed a flux density of 150 mJy for the highly variable
source M 81*. This value was chosen since it was in the range of typical values
at cm wavelenghts (e.g. Brunthaler et al. 2006) and it yielded flux
densities for 0945+6924, 1004+6936, and 44.0+59.6 in M 82 that were consistent
with their flux densities at the other two epochs.
M 82 was also observed on 2008 March 24 with the VLA at 22 GHz for 10 min
(6 min integration time on M 82) in spectral line mode. 3C 48 was
used as flux density calibrator, and 1048+717 was used as gain and phase
calibrator. A total bandwidth of 9.18 MHz was observed. We also analyzed the
archival VLA data of 2007 October 29 at 4.8 GHz, which is the latest available
data before our observations. Here, 3C 286 was used as flux density
calibrator and 1048+717 was used as phase calibrator. Eighteen minutes of
integration time on M 82 was spread over 1.5 h.
3 Results
A new bright source was detected on 2008 May 03 (see Fig. 1,
middle and Fig. 2, bottom). With its flux density of
90 mJy, it was clearly the brightest radio source in the field (at least five
times brighter than the second brightest source). Almost one year later, on
2009 April 08, the source had faded to a flux density of
11 mJy,
i.e. losing almost 90% of its flux density (see Fig. 1, bottom). We
estimated its position to be
55
51.551
40
45.792
.
Following the
traditional convention of previous papers naming a source after the
position offset to
55
and
,
our new detection would be 42.82+59.54. The
positions in the two observations are consistent and the uncertainty was
estimated by comparing the position of another clearly identified source
(44.0+59.6) with the positions from observations with the Multi-Element Radio
Linked Interferometer Network (MERLIN) by Muxlow et al. (1994).
![]() |
Figure 1:
Top: VLA B-configuration image at 4.8 GHz of M 82 on 2007
October 29. Contours start at 3 mJy and increase with factors of 2. The beamsize is
|
Open with DEXTER |
![]() |
Figure 2: Contours and grey scale represent 22.2 GHz images of the central region of M 82 taken on 2007 January 28 ( top) and 2008 May 03 ( bottom). The new bright source is very conspicuous in the latter image. Contour values are -4, 4, 8, 16, 32, 64, 128, 256, and 512 times 0.2 mJy beam-1, roughly the (comparable) rms noise level in both images. The star symbols mark the positions of the X-ray sources discussed by Matsumoto et al. (2001) and in the top panel bear these authors' nomenclature. The asterisk gives the position of the radio kinematic center (R.K.C.) of M 82 as determined by Weliachew et al. (1984). Both images were restored with a circular beam of 1.5 arcsec. Since the observation on 2007 January 28 was made in C configuration, there was much more data from short baselines, which results in better sensitivity to extended emission. |
Open with DEXTER |
Table 1: Details of the detected compact sources in our VLA observations with statistical errors on the 22.2 GHz flux densities.
To constrain the start time of the flare, we reduced two earlier VLA
observations of M 82. The spectral line observation on 2008 March 24 had lower
quality, since we observed with a narrow bandwidth and had only 6 min
integration time on M 82. Nevertheless, we could easily confirm the detection
of the new source at the same position with a similar but slightly higher flux
density as in the 2008 May 03 observation. Next, we analyzed 4.8 GHz data from
2007 October 29. Here, no bright point source was discovered (see
Fig. 1, top), so we conclude that the start of the flare was
between 2007 October 29 and 2008 March 24. However, it is possible that the
source could be highly self-absorbed early in its development and thus might
have been detectable at high frequencies, but not at 4.8 GHz. In this case,
the flare might have started earlier, but not earlier than 2007 January 28
(the day of our first 22 GHz observation, see Fig. 2, top).
The comparision of flux densities from different epochs in M 82 is difficult due to the diffuse emission, in particular when comparing observations with different spatial resolutions. In Table 1, we list the flux densities of the new source along with the flux densities of one additional source in M 82 (44.0+59.6), M 81*, and three background quasars.
4 Pre-flare observations at other wavelengths
Matsumoto et al. (2001) present high-resolution (
)
X-ray
imaging of the central
(
kpc) region of M 82
with the High-Resolution Camera (HRC) aboard the Chandra X-ray Observatory.
These images, taken on 1999 October 28 and 2000 January 20, show a total of
9 sources within this area, some of which are highly variable. They do
not find a counterpart to our new radio source. Neither do Kong et al. (2007) in
a total of 12 datasets of a similar sized region taken with the Chandra HRC and
Advanced CCD Imaging Spectrometer Array (ACIS-1 and -2) taken between 1999
September 20 and 2005 August 18. These authors also present Hubble Space
Telescope (HST) H-band (
m) imaging of the region with
the Near-Infrared Camera and Multi-Object Spectrometer (NICMOS) that, again,
does not show a counterpart to our variable radio source.
Körding et al. (2005) observed M 82 with the VLA at
8.4 GHz eight times between 2003 June and October in order to detect possible
radio flares
from several ultraluminous X-ray sources (ULX). No significant emission was
found at the position of our new detection above a noise level of 70 Jy
at each epoch. A very sensitive eight-day integration of M 82 at 5 GHz
wavelength was performed with MERLIN between 2002 April 1 and 28 (Fenech et al. 2008). The resulting
40 millarcsec resolution images show no
significant emission at the position of our new source above the rms noise
level of
Jy beam-1. Tsai et al. (2009) present a 7 mm map from
the VLA from 2005 April 22, with no detection of more than 0.3 mJy at
the position of our source.
5 Discussion
The most straightforward explanation for this new source is a new radio
supernova. A quantitative analysis of the lightcurve is difficult due
to the small number of data points (3) we obtained with different angular
resolutions, the complication of diffuse emission in M 82, and the uncertain
absolute flux scale in our observation of 2008 May 03. We extracted the
flux densities of our source by restricting the interferometer (u, v)-data to
>30 k
to remove most of the extended emission
(Table 1).
We fitted the lightcurve of the source with an exponential decay:
.
This yields a decay timescale of
yr
(Fig. 3). Since we do not know the exact time of the onset of the
flare, t0, we get values of S0 between 110 and 270 mJy for t0 between
2008 March 24 and 2007 October 29. For the fit, we added in quadrature the
difference between peak and integrated flux density as one error estimate and
systematic (flux density scale) errors of 5% for the first and third epochs,
and 10% for the second epoch, where no
proper flux density calibrator was observed. The resulting
was 0.5. Based on this fit, we predict that the source flux density will drop to
5 mJy by mid 2009 and 1.5 mJy in early 2010. We also fitted a simple exponential
decay to the 22 GHz lightcurve of SN199J in M 81. Here we used the first year
of data after the peak in the lightcurve published in
Weiler et al. (2007). The fit yields
yr,
indicating that the new source in M 82 decays faster than SN1993J.
A power-law fit (
)
is not consistent with the
lightcurve of the new source in M 82 (
,
for t0 between
2007.8 and 2008). If one allows an earlier t0 (e.g., if the source is highly
self-absorbed at lower frequencies), the power-law fit improves
(
,
for
t0=2007.08,
one day after our 2007 January 28
observation).
![]() |
Figure 3: Lightcurve of the new radio transient from our observations (squares). Also shown are an exponential-decay fit (solid line) and power-law fits (dashed lines), assuming different peak times of the flare: 2008.0 (lower line), 2007.8 (middle line), and 2007.08 (upper line). |
Open with DEXTER |
The position of the new source is marginally consistent with the position of
the kinematic center of M 82 (Weliachew et al. 1984). This raises
the possibility that, rather than a supernova, we could have detected a
flare from a supermassive black hole in the center of M 82. Flares in AGN
sources can often be fitted with exponential decays
(Valtaoja et al. 1999). For
example, a strong radio flare in 1999 in the Seyfert galaxy III Zw 2 has a
decay rate of
year (Brunthaler et al. 2005a).
Other explanations, such as luminous flares from quiescent supermassive black
holes induced by a close passage of a star that is torn apart by tidal forces,
are also possible. However, such flares show a power-law decay with
(Ayal et al. 2000; Evans & Kochanek 1989; Gezari et al. 2009), which is not
consistent with our lightcurve.
However, M 82 has never shown evidence of a nuclear supermassive black hole (which would be surprising for a smallirregular galaxy). Since the progenitor for our flare showed no X-ray emission, a stellar or intermediate mass black hole is also not probable. Thus, based on the current data, a new radio supernova seems to be the most likely explanation.
Acknowledgements
We thank the referee Dr. K. Weiler for critically reading the manuscript.
References
- Ayal, S., Livio, M., & Piran, T. 2000, ApJ, 545, 772 [NASA ADS] [CrossRef]
- Bartel, N., Bietenholz, M. F., Rupen, M. P., et al. 2002, ApJ, 581, 404 [NASA ADS] [CrossRef]
- Bartel, N., Bietenholz, M. F., Rupen, M. P., & Dwarkadas, V. V. 2007, ApJ, 668, 924 [NASA ADS] [CrossRef]
- Beswick, R. J., Muxlow, T. W. B., Argo, M. K., & Pedlar, A. 2004, IAU Circ., 8332, 2 (In the text)
- Beswick, R. J., Riley, J. D., Marti-Vidal, I., et al. 2006, MNRAS, 369, 1221 [NASA ADS] [CrossRef]
- Bietenholz, M. F., Bartel, N., & Rupen, M. P. 2001, ApJ, 557, 770 [NASA ADS] [CrossRef]
- Bietenholz, M. F., Bartel, N., & Rupen, M. P. 2003, ApJ, 597, 374 [NASA ADS] [CrossRef]
- Brunthaler, A., Falcke, H., Bower, G. C., et al. 2005a, A&A, 435, 497 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Brunthaler, A., Reid, M. J., Falcke, H., Greenhill, L. J., & Henkel, C. 2005b, Science, 307, 1440 [NASA ADS] [CrossRef]
- Brunthaler, A., Bower, G. C., & Falcke, H. 2006, A&A, 451, 845 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Brunthaler, A., Reid, M. J., Falcke, H., Henkel, C., & Menten, K. M. 2007, A&A, 462, 101 [NASA ADS] [CrossRef] [EDP Sciences]
- Chevalier, R. A., & Fransson, C. 2001, ApJ, 558, L27 [NASA ADS] [CrossRef] (In the text)
- Evans, C. R., & Kochanek, C. S. 1989, ApJ, 346, L13 [NASA ADS] [CrossRef]
- Fenech, D. M., Muxlow, T. W. B., Beswick, R. J., Pedlar, A., & Argo, M. K. 2008, MNRAS, 391, 1384 [NASA ADS] [CrossRef]
- Freedman, W. L., Hughes, S. M., Madore, B. F., et al. 1994, ApJ, 427, 628 [NASA ADS] [CrossRef] (In the text)
- Gezari, S., Heckman, T., Cenko, S. B., et al. 2009, ApJ, in press [arXiv:0904.1596]
- Kong, A. K. H., Yang, Y. J., Hsieh, P.-Y., Mak, D. S. Y., & Pun, C. S. J. 2007, ApJ, 671, 349 [NASA ADS] [CrossRef] (In the text)
- Körding, E., Colbert, E., & Falcke, H. 2005, A&A, 436, 427 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Kronberg, P. P., & Sramek, R. A. 1985, Science, 227, 28 [NASA ADS] [CrossRef] (In the text)
- Kronberg, P. P., Sramek, R. A., Birk, G. T., et al. 2000, ApJ, 535, 706 [NASA ADS] [CrossRef] (In the text)
- Marcaide, J. M., Alberdi, A., Ros, E., et al. 1997, ApJ, 486, L31 [NASA ADS] [CrossRef]
- Marcaide, J. M., Marti-Vidal, I., Alberdi, A., et al. 2009, A&A, submitted [arXiv:0903.3833]
- Matsumoto, H., Tsuru, T. G., Koyama, K., et al. 2001, ApJ, 547, L25 [NASA ADS] [CrossRef] (In the text)
- Mattila, S., Meikle, W. P. S., Groeningsson, P., et al. 2004, IAU Circ., 8299, 2 (In the text)
- Muxlow, T. W. B., Pedlar, A., Wilkinson, P. N., et al. 1994, MNRAS, 266, 455 [NASA ADS]
- Pérez-Torres, M. A., Alberdi, A., & Marcaide, J. M. 2001, A&A, 374, 997 [NASA ADS] [CrossRef] [EDP Sciences]
- Pérez-Torres, M. A., Alberdi, A., & Marcaide, J. M. 2002, A&A, 394, 71 [NASA ADS] [CrossRef] [EDP Sciences]
- Schmidt, B. P., Kirshner, R. P., Eastman, R. G., et al. 1993, Nature, 364, 600 [NASA ADS] [CrossRef] (In the text)
- Seaquist, E. R., & Stankovic, M. 2007, ApJ, 659, 347 [NASA ADS] [CrossRef] (In the text)
- Singer, D., Pugh, H., & Li, W. 2004, IAU Circ., 8297, 2 (In the text)
- Tsai, C.-W., Turner, J. L., Beck, S. C., Meier, D. S., & Ho, P. T. P. 2009, AJ, in press [arXiv:0903.1858] (In the text)
- Valtaoja, E., Lähteenmäki, A., Teräsranta, H., & Lainela, M. 1999, ApJS, 120, 95 [NASA ADS] [CrossRef] (In the text)
- van Buren, D., & Greenhouse, M. A. 1994, ApJ, 431, 640 [NASA ADS] [CrossRef] (In the text)
- Weiler, K. W., Panagia, N., Montes, M. J., & Sramek, R. A. 2002, ARA&A, 40, 387 [NASA ADS] [CrossRef] (In the text)
- Weiler, K. W., Williams, C. L., Panagia, N., et al. 2007, ApJ, 671, 1959 [NASA ADS] [CrossRef] (In the text)
- Weliachew, L., Fomalont, E. B., & Greisen, E. W. 1984, A&A, 137, 335 [NASA ADS] (In the text)
- Yun, M. S., Ho, P. T. P., & Lo, K. Y. 1994, Nature, 372, 530 [NASA ADS] [CrossRef] (In the text)
Footnotes
- ... NRAO
- The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc.
All Tables
Table 1: Details of the detected compact sources in our VLA observations with statistical errors on the 22.2 GHz flux densities.
All Figures
![]() |
Figure 1:
Top: VLA B-configuration image at 4.8 GHz of M 82 on 2007
October 29. Contours start at 3 mJy and increase with factors of 2. The beamsize is
|
Open with DEXTER | |
In the text |
![]() |
Figure 2: Contours and grey scale represent 22.2 GHz images of the central region of M 82 taken on 2007 January 28 ( top) and 2008 May 03 ( bottom). The new bright source is very conspicuous in the latter image. Contour values are -4, 4, 8, 16, 32, 64, 128, 256, and 512 times 0.2 mJy beam-1, roughly the (comparable) rms noise level in both images. The star symbols mark the positions of the X-ray sources discussed by Matsumoto et al. (2001) and in the top panel bear these authors' nomenclature. The asterisk gives the position of the radio kinematic center (R.K.C.) of M 82 as determined by Weliachew et al. (1984). Both images were restored with a circular beam of 1.5 arcsec. Since the observation on 2007 January 28 was made in C configuration, there was much more data from short baselines, which results in better sensitivity to extended emission. |
Open with DEXTER | |
In the text |
![]() |
Figure 3: Lightcurve of the new radio transient from our observations (squares). Also shown are an exponential-decay fit (solid line) and power-law fits (dashed lines), assuming different peak times of the flare: 2008.0 (lower line), 2007.8 (middle line), and 2007.08 (upper line). |
Open with DEXTER | |
In the text |
Copyright ESO 2009
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.