Issue |
A&A
Volume 517, July 2010
|
|
---|---|---|
Article Number | A46 | |
Number of page(s) | 14 | |
Section | Astrophysical processes | |
DOI | https://doi.org/10.1051/0004-6361/200913613 | |
Published online | 03 August 2010 |
Coordinated NIR/mm observations of flare emission from Sagittarius A*
D. Kunneriath1,2 - G. Witzel1 - A. Eckart1,2 - M. Zamaninasab2,1 - R. Gießübel2,1 - R. Schödel3 - F. K. Baganoff4 - M. R. Morris5 - M. Dovciak6 - W. J. Duschl7,8 - M. García-Marín1 - V. Karas6 - S. König1 - T. P. Krichbaum2 - M. Krips12 - R.-S. Lu2,1 - J. Mauerhan16 - J. Moultaka9 - K. Muzic1 - N. Sabha1 - F. Najarro10 - J.-U. Pott11 - K. F. Schuster12 - L. O. Sjouwerman13 - C. Straubmeier1 - C. Thum12 - S. N. Vogel14 - P. Teuben14 - A. Weiss2 - H. Wiesemeyer15 - J. A. Zensus2,1
1 - I.Physikalisches Institut, Universität zu Köln, Zülpicher Str. 77,
50937 Köln, Germany
2 - Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121
Bonn, Germany
3 - Instituto de Astrofísica de Andalucía (CSIC), Camino Bajo de Huétor
50, 18008 Granada, Spain
4 - MKI, Massachusetts Institute of Technology, Cambridge,
MA 02139-4307, USA
5 - Department of Physics and Astronomy, University of California, Los
Angeles, CA 90095-1547, USA
6 - Astronomical Institute, Academy of Sciences, Bocní II, 14131
Prague, Czech Republic
7 - Institut für Theoretische Physik und Astrophysik,
Christian-Albrechts-Universität zu Kiel, Leibnizstr. 15, 24118 Kiel,
Germany
8 - Steward Observatory, The University of Arizona, 933 N. Cherry Ave.
Tucson, AZ 85721, USA
9 - LATT, Université de Toulouse, CNRS, 14 Avenue Edouard Belin, 31400
Toulouse, France
10 - Departamento de Astrofísica, Centro de Astrobiología, CSIC-INTA,
Ctra. Torrejón a Ajalvir km 4, 28850 Torrejn de Ardoz, Spain
11 - Max-Planck-Institut für Astronomie, Königstuhl 17, 69117
Heidelberg, Germany
12 - Institut de Radio Astronomie Millimetrique, Domaine Universitaire,
38406 Saint Martin d'Heres, France
13 - National Radio Astronomy Observatory, PO Box 0, Socorro, NM 87801,
USA
14 - Department of Astronomy, University of Maryland, College Park, MD
20742-2421, USA
15 - IRAM, Avenida Divina Pastora 7, Núcleo Central, 18012 Granada,
Spain
16 - IPAC, California Institute of Technology, 770 South Wilson Avenue,
Pasadena, CA 91125, USA
Received 5 November 2009 / Accepted 29 March 2010
Abstract
Context. We report on a successful, simultaneous
observation and modelling of the millimeter (mm) to near-infrared (NIR)
flare emission of the Sgr A* counterpart associated with the
supermassive (
) black hole at the Galactic centre (GC). We present a mm/sub-mm light
curve of Sgr A* with one of the highest quality continuous
time coverages.
Aims. We study and model the physical processes
giving rise to the variable emission of Sgr A*.
Methods. Our non-relativistic modelling is based on
simultaneous observations carried out in May 2007 and 2008, using the
NACO adaptive optics (AO) instrument at the ESO's VLT and the mm
telescope arrays CARMA in California, ATCA in Australia, and the
30 m IRAM telescope in Spain. We emphasize the importance of
multi-wavelength simultaneous fitting as a tool for imposing adequate
constraints on the flare modelling. We present a new method for
obtaining concatenated light curves of the compact mm-source
Sgr A* from single dish telescopes and interferometers in the
presence of significant flux density contributions from an extended and
only partially resolved source.
Results. The observations detect flaring activity in
both the mm domain and the NIR. Inspection and modelling of the light
curves show that in the case of the flare event on 17 May 2007, the mm
emission follows the NIR flare emission with a delay of h.
On 15 May 2007, the NIR flare emission is also followed by elevated
mm-emission. We explain the flare emission delay by an adiabatic
expansion of source components. For two other NIR flares, we can only
provide an upper limit to any accompanying mm-emission of about
0.2 Jy. The derived physical quantities that describe the
flare emission give a source component expansion speed of
,
source sizes of about one Schwarzschild radius, flux densities of a few
Janskys, and spectral indices of
to 1.3. These source components peak in the THz regime.
Conclusions. These parameters suggest that either
the adiabatically expanding source components have a bulk motion
greater than
or the expanding material contributes to a corona or disk, confined to
the immediate surroundings of Sgr A*. Applying the flux
density values or limits in the mm- and X-ray domain to the observed
flare events constrains the turnover frequency of the synchrotron
components that are on average not lower than about 1 THz,
such that the optically thick peak flux densities at or below these
turnover frequencies do not exceed, on average, about
1 Jy.
Key words: black hole physics - infrared: general - accretion, accretion disks - Galaxy: center - galaxies: nuclei - radio continuum: general
1 Introduction
Sgr A*, the compact non-thermal radio and infrared source at
the centre of the Milky Way galaxy (8 kpc away) is known to be associated
with a supermassive black hole (SMBH) of mass
(Ghez
et al. 2009; Eckart & Genzel 1996; Schödel
et al. 2002; Gillessen et al. 2009;
Eisenhauer
et al. 2005; Eckart et al. 2002; Schödel
et al. 2003; Ghez et al. 2000,1998; Genzel
et al. 1997; Ghez et al. 2004b; Genzel
et al. 2000; Eisenhauer et al. 2003;
Ghez
et al. 2005,2004a). The close proximity of
Sgr A* makes it ideal for studying the evolution and physics
of SMBHs located in the nuclei of galaxies. The SMBH radiates far below
its Eddington luminosity at all wavelengths, partly because of its low
observed accretion rate. For Sgr A*, we assume that
as,
being one
Schwarzschild radius and
the gravitational radius of the SMBH.
Evidence of flaring activity occurring from a few hours to
days has
been found from variability studies ranging from the radio to sub-mm
wavelengths (Mauerhan
et al. 2005; Zhao et al. 2004; Bower
et al. 2002; Herrnstein et al. 2004;
Zhao
et al. 2003). There is also evidence that variations
in radio/sub-mm emission are linked to NIR/X-ray flares, with the
radio/sub-mm flares occurring after a delay of 100 min after the NIR/X-ray flares (Eckart
et al. 2004,2008b; Yusef-Zadeh
et al. 2008; Eckart et al. 2006a;
Marrone
et al. 2008).
These flares have been explained with a synchrotron self Compton (SSC) model that involves up-scattered sub-mm photons from a compact (inferred from short flare timescales) source component (e.g., Eckart et al. 2004,2006a). The X-ray emission is caused by the inverse Compton scattering of the THz-peaked flare spectrum by relativistic electrons. Both synchrotron and SSC mechanisms contribute to NIR flux density. Adiabatic expansion of the source components then gives rise to flares in the mm/sub-mm regimes (Yusef-Zadeh et al. 2006a,b; Marrone et al. 2008; Eckart et al. 2006a; Yusef-Zadeh et al. 2008).
Based on the assigned VLT time, we organized an extensive
multi-frequency
campaign in May 2007 and May 2008 which included millimeter to NIR
observations at
single telescopes and interferometers around the world. We present data
from observations of Sgr A* using CARMA(Bock et al. 2006), ATCA
, and the MAMBO bolometer
at
the IRAM
30 m
telescope in the mm regime and NIR data from the ESO VLT. We detected
simultaneous emission in the NIR and mm-regimes using the CARMA and ESO
VLT telescopes. The main results are: 2 bright NIR flares
(16 mJy and 10 mJy), a
0.4 Jy mm flare, and a possible
weaker third flare in the NIR and mm in the combined ATCA, CARMA,
IRAM 30 m mm/sub-mm light curve from May 2007, and a
bright NIR flare in May 2008 also covered by the CARMA 3 mm
observations. We present updated versions of light curves first
presented in Kunneriath
et al. (2008), a more detailed description of the
methods used to obtain them, and physical models to explain the flaring
activity. The observations and data reduction are described in
Sect. 2, an outline of the flare modelling in Sect. 3
and a discussion and summary follow in Sects. 4
and 5.
2 Observations and data reduction
Interferometric observations in the mm/sub-mm wavelength domain are especially well suited to differentiating the flux density contribution of Sgr A* from the thermal emission of the circumnuclear disk (CND, a ring-like structure of gas and dust surrounding the Galactic centre at a distance of about 1.5-4 pc (see, e.g., Guesten et al. 1987; or Christopher et al. 2005). In May 2007 and 2008, global coordinated multiwavelength observations were carried out in the NIR and mm regimes to study the variability of Sgr A*.
2.1 The mm data
We observed the GC at 100 and 86 GHz (3 and 3.5 mm
wavelength) with the two mm-arrays CARMA and ATCA, respectively. In
addition, we observed with the MAMBO 2 bolometer at the IRAM
30 m-telescope at a wavelength of 1.2 mm. CARMA is
located in Cedar Flat, Eastern California, and consists
of 15 antennas (
m and
m).
The quasar source 3C273 was used for bandpass calibration while
1733-130 and Uranus were used for phase and amplitude calibration of
Sgr A* data, respectively. The Australia Telescope Compact
Array (ATCA), at the Paul Wild Observatory, is an array of six 22-m
telescopes located in Australia. Calibrator sources 1253-055, 1921-293,
and Uranus were used for bandpass, phase, and amplitude calibration.
The interferometer data were mapped using the Miriad
interferometric data reduction package. Details of the observation are
given in Table 1.
The Max-Planck Millimeter Bolometer (MAMBO 2) array is installed at the IRAM 30 m telescope on Pico Veleta, Spain. The 37 channel array of the precursor instrument MAMBO has been successfully used by many observers since the end of 1998. The bolometer data was reduced using the bolometer array data reduction, analysis, and handling software package, the BoA (Bolometer Data Analysis).
Table 1: Log of the mm- and sub-mm observations.
To correct for extended flux contributions in the
interferometer data, we extracted visibilities from two orthogonal
pairs of the longest baselines for the two arrays, and subtracted for
each baseline the median baseline and time dependent visibility trend
![]() |
(1) |
from each visibility data set



![]() |
(2) |
and
![]() |
(3) |
Here the operator


![]() |
Figure 1: Combined differential light curve of Sgr A* in the mm/sub-mm domain for the May 2007 observing run. The MAMBO 2 bolometer at the IRAM 30 m-telescope was operated at a central wavelength of 1.2 mm (250 GHz). The CARMA data were centred on 100 GHz and the ATCA data on 86 GHz. We also show the daily 7 mm and 3 mm flux density averages from simultaneous VLBA observations (Lu et al. 2008,2009). In the lower part of the plot, we show the combined NIR light curve consisting of K-band and L'-band light curves (see Appendix for individual light curves). The time axis is labelled with UT hours starting at 00 h on May 15. |
Open with DEXTER |
![]() |
Figure 2: Difference maps at 3 mm of a 40'' diameter region (shown by the dashed white line) centred on Sgr A*, obtained from the difference between full synthesis maps of the individual days of CARMA observations and the full CARMA data set as described in the text. The maps have been plotted with the same color-coding table. The figure shows that the flux density variations that are evident from the differential light curves (see Fig. 1) can also be seen in the maps constructed from the corresponding data. |
Open with DEXTER |
The influence of flux spectral index variations can be calculated in
the following way.
If S1 and S2
are the flux densities at the frequencies
and
,
then a change from spectral index
to
results
in a new flux density value S'1
and with
![]() |
(4) |
![]() |
(5) |
we find
![]() |
(6) |
or for S'1 =f S1
![]() |
(7) |
We estimate a spectral index variation of




The light curve in Fig. 1 contains two peaks, on May 15 and 17 (there is a weaker, third possible peak on May 19). In Fig. 1, we also show the daily flux density averages of the 7 mm VLBA observations that were conducted in parallel (Lu et al. 2008; Kunneriath et al. 2008; Lu et al. 2009). The VLBA data follow the overall trend of the combined CARMA/ATCA/30 m light-curve very well. The NIR coverage in the K and L'-band is given in the lower panel of Fig. 1 (see Appendix for a description of the scaling and the individual light curves for each day).
To verify the mm/sub-mm flux density variations and in particular the flare detected on 17 May by CARMA, we produced residual maps from the four individual tracks obtained with the CARMA array, since CARMA has the best combination of high resolution and signal-to-noise ratio of all the mm-telescopes we use. These maps shown in Fig. 2 were computed by subtracting the mean of all 4 maps from the maps of individual epochs. The rms noise in these differential maps is of the order of 0.1 Jy per beam. This procedure clearly represents the trend shown in the combined light curve and detects the excess flux density of 0.4 Jy detected on May 17 and a slightly positive flux density on May 19. On the 16th and 18th May 2007, a mixture of negative and positive flux densities in the difference maps corresponds to the non detection of flaring activity in the differential light curve shown in Fig. 1.
Table 2: Log of the near-infrared observations.
![]() |
Figure 3:
The 3 mm CARMA and ATCA (top red and blue, respectively) and K-
and L'-band NIR (middle red and blue, respectively)
light curves of the 17 May 2007 flare, with the cross-correlation
(bottom) between the two showing that the mm flare follows the NIR
flare with a time lag of |
Open with DEXTER |
![]() |
Figure 4: 3 mm CARMA Sgr A* ( top) and calibrator (1733-130) ( middle) light curves from 17 May 2007, with the cross-correlation ( bottom) between the two showing zero correlation between the calibrator and the source fluxes. |
Open with DEXTER |
2.2 The NIR data
The NIR data were taken in the -band (2.18
m, FWHM
0.35
m)
and in the L'-band (3.8
m, FWHM
0.62
m)
at ESO's Very Large Telescope (VLT) with the NACO infrared camera and
adaptive optics (AO) system (Lenzen et al. 2003; Rousset
et al. 2003) in Chile. For the 3.8
m
observations, the integration times were NDIT
DIT =
s
(DIT is Detector Integration Time, and NDIT is the number of DITs). At
2.2
m,
the integration times were NDIT
DIT =
s.
Additional details of the observation are given in Table 2.
The diffraction limit in K-band is about
60 mas. The infrared wavefront sensor of NAOS was used to lock
the AO loop on the NIR bright (K-band magnitude
6.5)
supergiant
IRS 7, located about 5.6'' north of Sgr A*. Therefore
the AO was able to provide a stable correction with a high Strehl ratio
(of the order of 50%). After correcting the images for bad pixels, and
subtracting the sky and flat fielding, the point spread function (PSF)
was extracted from each individual image using
(Diolaiti et al. 2000)
and a Lucy-Richardson (LR) deconvolution was applied. A Gaussian beam
with FWHM corresponding to the respective wavelength (
60 mas
at 2.2
m
and
104 mas
at 3.8
m)
was used to restore the beam.
To obtain a flux density in Jansky, we perform aperture photometry of
the sources with circular apertures of radius 52 mas, apply extinction
corrections of AK=2.8
and AL'=1.8,
and calibrate the flux with flux densities of known sources, which
resulted in K- and L'-band flux
densities of S2 of
mJy
and
mJy,
respectively. The measurement uncertainties were obtained from the
reference star S2.
2.3 Flare events
2.3.1 The 15 May 2007 flare event
The first NIR flare on May 15 detected during the multiwavelength campaign (see also Eckart et al. 2008a) preceded the combined mm/sub-mm monitoring and the first maximum detected therein by nearly 3 h. The May 15, 2007 flare shows 2 polarized bright sub-flares centred on about 07:29 UT and 07:51 UT, respectively, after the start of observations (see Fig. A.1 in Appendix). We define sub-flares to be the shorter flux excursions superimposed on the longer main underlying flare. The time difference between the two sub-flares is 22 min which is fully consistent with previously reported sub-flare separations. Zamaninasab et al. (2010) demonstrated that these highly polarized sub-flares are part of the flare structure that is significant compared to the randomly polarized red-noise.
The second sub-flare shows substructure that is interpreted by Eckart et al. (2008a) in the framework of spot evolution. The flux density between the 2 sub-flares does not reach the emission level well before and after (i.e., <70 and >170 min into the observations) the flare, probably because of fore- or background stellar flux density contributions (star S17). Additional details of the May 15, 2007 NIR flare were given by Eckart et al. (2008a).
For the flux density limit reached in low luminosity states of Sgr A*, the reader is referred to Sabha et al. (2010) and Do et al. (2009). The brief L'-band observations at 10 UT show a NIR signal that is a factor of 2 times stronger than at other times (see Sect. 3.3.1).
2.3.2 The 17 May 2007 flare event
The NIR flare on 17 May had overlap with the CARMA observations, which indicate that a 0.4 Jy 3 mm flare followed the NIR flare.
Figure 3
gives the updated differential light curve
obtained from the four longest baselines of CARMA and ATCA, along with
the NIR light curve and a cross-correlation between the two showing a
time lag of h,
with a significant, broad and positive peak corresponding to the flares
found at both wavelengths.
The presence of negative power in the cross-correlation is most likely
due to a combination
of two effects: (1) the signal-to-noise ratio of the differential ATCA
data is lower than that of the CARMA data; and (2) the negative power
may also be indicative of
residual power on timescales of the length of the (sub-)mm-datasets.
Some amount of negative power is expected since we are dealing with
differential light curves that measure the variable flux with respect
to a median flux level over the
length of the corresponding dataset. A significant, broad, and positive
cross-correlation peak centred on a time lag of
1.5 h is also present if the correlation
is only performed with the K-band and CARMA data.
Figure 4
shows that the calibrator and source fluxes are not correlated for the
CARMA data.
Here the cross-correlation was carried out only for the CARMA data
since the ATCA array did not observe the same calibrator interleaved
with the Sgr A* data.
The NIR K-band data show variable emission, with four sub-flares starting from about 7:00 UT and lasting until the end of the observations at 9:12 UT. The mm-data started at about 8:00 UT and had ended by 10:00 UT. In the mm domain, the variable emission is dominated by a single flare that is slightly wider than the individual NIR sub-flares.
2.3.3 The 19 May 2007 flare event
The flux density excursions on 19 May are only partially constrained by the CARMA millimeter measurements. Two flare events of about 3.5 mJy and 1.1 mJy peak flux density are separated by about 3 h. The gap in the observations at a time during which the first event reached its peak flux density (see Appendix for individual light curve) implies that the true peak was missed and that the flare was probably brighter than 3.5 mJy. There is an overlap between the NIR and the mm-measurements of about 1.5 h with no detection of a flare event in the mm-wavelength domain brighter than 0.1-0.2 Jy.
2.3.4 The 26 May 2008 flare event
The 26 May 2008 NIR flare in the L'-band was covered by the CARMA observations, as shown in Fig. 5. The L'-band flare shows a flux density increase of about 70 mJy, followed by a plateau with almost constant flux density at a 40 mJy level. After a total flare duration of 200 min, the flux density reaches the original level of 20 mJy. The CARMA data were taken under moderate to poor weather conditions. From 07:10 UT to 11:10 UT, the source visibility was comparable to the levels reached in the 2007 observing session. No flare was detected with a flux density limit of about 0.2 Jy. Starting at 11:10 UT, the weather conditions deteriorated resulting in strong coherence losses.
![]() |
Figure 5:
3.8 |
Open with DEXTER |
3 Flare analysis
The seven coordinated Sgr A* measurements reported so far that
include sub-mm/mm data (Yusef-Zadeh
et al. 2009; Eckart et al. 2008b;
Yusef-Zadeh
et al. 2006b; Eckart et al. 2006a;
Marrone
et al. 2008,
including this work) have shown that the observed sub-millimeter/mm
flares follow the largest event observed at shorter wavelengths
(NIR/X-ray; see detailed discussion in Marrone
et al. 2008) with a delay of h.
We therefore assume that the millimeter flares presented here are
related to the observed IR flare events. This flaring activity can be
explained using a multi-component relativistic spot/disk model and an
adiabatic expansion model in combination with a synchrotron self
Compton (SSC) formalism. While this model has been presented before
(see references above),
it is essential to test new datasets against it to improve
our understanding of the variable emission from Sgr A*.
3.1 Relativistic spot/disk modelling of NIR flares
Multiple hot spots revolving around a black hole in an accretion disk in AGNs or galactic black holes can give rise to light curves in the X-ray and NIR regimes whose power spectral density (PSD) is described by a broken power law with a slope similar to that of red noise processes (Armitage & Reynolds 2003; Pechácek et al. 2008; Do et al. 2008). Simulations show that polarized light curves exhibit behaviour associated with lensing of hot spots orbiting around a SMBH (Zamaninasab et al. 2010). Previously observed NIR light curves of Sgr A* have been successfully modelled by a multi-component hot spot model involving source components orbiting around the SMBH in a temporary accretion disk (Zamaninasab et al. 2008b; Meyer et al. 2006b,a; Eckart et al. 2006b; Zamaninasab et al. 2008a; Meyer et al. 2007).
This relativistic model consists of source components that are
a mix of synchrotron and SSC components with an optically thin spectral
index
and relativistic electrons with
,
using the formalism described by Gould
(1979) and Marscher
(1983). The KY-code by Dovciak
et al. (2004) produces light amplification curves
for each individual component orbiting the SMBH taking into account
relativistic effects, and by combining these curves with the SSC model,
we can estimate the SSC X-ray and NIR flux densities and the magnetic
field (Eckart et al.
2008a). Figure 6 presents
simulations of the accretion disk with multiple spots revolving around
the SMBH, and Fig. 7
gives the light curves produced by these simulations. Zamaninasab et al. (2010)
describe the modelling and simulations shown in Figs. 6 and 7 in greater
detail.
These light curves may be the result of a single hot spot that dominates the disk for several orbital periods. If the spot were to sink (as may occur) towards the centre while orbiting, this would cause a quasi-periodic modulation of the light curve for the time the emission of the spot dominates the red-noise emission of the remaining disk. This modulation may also be produced by several different hot spots that have individual lifetimes of shorter than one orbital period close to the last stable orbit. In this case, the light-curve will be modulated quasi-periodically at a rate close to the frequency of the last stable orbit (see discussion in Eckart et al. 2008a; and Eckart et al. 2006b). If the hot spots expand within the accretion disk or are the source of a short outflow above the disk, this quasi-periodic signal will be smeared out. In particular, it may be undetectable at lower (sub-)mm observing frequencies (see below). In the case of our experiment, only the NIR data on 15 May shows a strong sub-flare modulation (as explained below) that may be linked to quasi-periodicity. All other NIR flares (with the possible exception of the steeply rising and falling flanks of the L'-band flare on 26 May 2008) cannot be associated with this phenomenon. However, in the case of a source expansion this model can be linked to the adiabatic expansion model for synchrotron sources that is described in the following section.
![]() |
Figure 6:
Simulated images of multiple spots revolving around a black hole,
projected on the observer image plane ( |
Open with DEXTER |
![]() |
Figure 7: Light curve obtained from the multiple spot simulations shown in Fig. 6, the solid black line representing the overall light curve and the different coloured dashed lines showing the contributions of individual spots. |
Open with DEXTER |
3.2 Adiabatically expanding source components
Adiabatic expansion of the synchrotron components can explain the
apparent time difference of h between the NIR
and sub-mm/mm flares (Yusef-Zadeh
et al. 2006a; Eckart et al. 2008b,2006a).
We assume a uniformly expanding blob of relativistic electrons with a
power-law energy spectrum,
,
threaded by a magnetic field that declines as R-2,
and of energy and density that decline as R-1
and R-3, respectively, as a
result of expansion of the blob (van
der Laan 1966). If R0,
S0 and
0
are the size, flux density, and optical depth of the source component
at the peak frequency (
0)
of the synchrotron spectrum, the optical depth and flux density at a
given frequency
scale as
and
respectively. Here


For instance, if p ranges from 1 to 3,














Table 3: Source component parameters for the combined SSC and adiabatic expansion model of the 15, 17, and 19 May 2007 and the 26 May 2008 flares.
3.3 Adiabatic expansion modelling
Table 3
summarises the properties of the model, such as the adiabatic expansion
velocity ,
the optically thin spectral index
,
cutoff frequency
,
and flux density
of the source components, using the same nomenclature as Eckart
et al. (2009,2006a). The optically thin
NIR flux density is represented by
,
while
and
give the upper limits to the flux densities of upscattered SSC
components.
A global variation in a single parameter by the value listed in the
corresponding column in Table 3 results in an
increase of
(from the reduced
value). Here global variation means adding to a single model parameter
for all source components the 1
uncertainty, such that a maximum positive or negative flux density
deviation is reached. Alternatively, a variation by the listed
uncertainty for only a single source component results in a variation
of the model
predicted NIR and X-ray flux density by more than 30%.
Judging from the
based on the mm-data alone,
the global uncertainties for
,
,
and R0 may then be doubled.
When performing the reduced
fit for N source components, we used N times
4 (
,
,
R0,
)
plus one common expansion velocity
and
time offset (leaving the time differences between the components
fixed),
i.e., 4N+2 degrees of freedom.
Unfortunately, the model parameters may not all be considered
independent, e.g., the width and peak of a light curve signature
depends to a varying extent on all 4 parameters
,
,
R0, and
.
Therefore, for all models we stayed with the minimum number of source
components to estimate the degrees of freedom. Since in the case of the
May 2007 flare the VLT data consist of
10 times the number of CARMA data
points, we weighted the squared CARMA flux deviations and number of
data points by an additional factor of 19. The
test was then carried out using the sum of the squared flux deviations
and data points of the VLT and CARMA datasets.
In general, we aimed to develop a model at a lower and a higher expansion speed labelled A and B (Table 3). In comparison to earlier modelling results (Eckart et al. 2008a,2009; Marrone et al. 2008; Yusef-Zadeh et al. 2008), these speeds were chosen to be 0.007c and 0.017cfor the 15 May and 19 May 2007 flares. For the 17 May 2007 and the 26 May 2008 data, a violation of the observed flare width and the expected range for the magnetic field strength restricted the choice of expansion speeds. For comparison, we listed the model components in Table 3. The models were developed by minimizing the number of free parameters (and maximizing the description of significant flare features in the observed light curves).
SSC modelling as an additional constraint:
The SSC model described in Sect. 3.1
allows us to estimate the SSC contribution of the X-ray and NIR flux
densities, and the magnetic field. Higher SSC NIR emission would
violate the assumption in our adiabatic expansion model that the
dominant source of NIR flux density is synchrotron emission of the THz
peaked expanding component, as well as require us to model the unknown
X-ray flux density. Magnetic fields of up to 70 Gauss (Eckart et al. 2008b;
Marrone
et al. 2008; Eckart et al. 2006a;
Yusef-Zadeh
et al. 2008) have been obtained for previous models
via
,
where
and
are the synchrotron turnover frequency and flux density, respectively.
These results should also not be violated by the model.
The expansion speed:
The model yields expansion velocities from 0.005c-0.017c
(see Table 3),
which agree well with previously published results (Yusef-Zadeh et al. 2008).
These velocities are lower than expected for relativistic sound speed
in orbital velocity near the SMBH. This may be caused by the bulk
velocity of the source components being larger than the expansion
velocity
or because the expanding material is confined to the immediate region
surrounding Sgr A* in the form of a disk or corona. In this
case, shearing caused by differential rotation within the accretion
disk is responsible for the ``expansion'' and the observed low
expansion velocities (Pechácek et al. 2008;
Zamaninasab
et al. 2008b; Eckart et al. 2008b).
Owing to the expansion or the presence of several strong spots in the
disk or corona, the short time modulation often seen in the NIR will
also
be smeared out in the (sub-)mm-wavelength domain - implying that it
may not be observable at all at these wavelengths.
By comparing other flare events that were detected
simultaneously
at X-ray and sub-mm/mm-wavelengths, we imposed the additional
constraint that the 100 GHz flux density should not be much
stronger than about 0.2 Jy, h after the NIR
event, and the X-ray flux density should be limited to a few
10 nJy. The following sections describe the modelling of the
individual flares.
3.3.1 Expansion modelling of the 15 May flare emission
The 15 May 2007 flare can only be constrained by means of its NIR
lightcurve (Fig. 8)
since the mm data begins nearly 3 h after the first peak in
the NIR. According to the adiabatic expansion model, 3-4 h
after the NIR flare, the flux density in the mm band was insufficiently
strong to be significant compared to the median flux density variations
in the data. We modelled a background component
and two components for the two strongest sub-flares
and
.
The expansion speed could not be constrained at all. In model A, we
used a value of 0.007c, close to the speed we found
for most of the modelled flares. In model B, we used a higher speed of
0.017c, which required lower
spectral indices and a weaker background component to obtain a good fit
to the NIR data
and to fullfil the upper limits at 100 GHz and in the X-ray
domain.
![]() |
Figure 8: Models A and B of the 15 May 2007 NIR flare: a) fitted NIR light curve, with the black curve corresponding to the model and the red curve corresponding to the data, b) contribution of individual source components to NIR light curve, and c) contribution of individual source components to mm light curve. |
Open with DEXTER |
![]() |
Figure 9: Models A and B of the 17 May 2007 NIR flare: a) fitted NIR light curve, with the black curve corresponding to the model and the red curve corresponding to the data; b) contribution of individual source components to NIR light curve; c) fitted mm light curve, with the black curve corresponding to the model and the red curve corresponding to the data; and d) contribution of individual source components to mm light curve. |
Open with DEXTER |
![]() |
Figure 10: Models A and B of the 19 May 2007 NIR flare: a) fitted NIR light curve, with the black curve corresponding to the model and the red curve corresponding to the data; b) contribution of individual source components to NIR light curve; c) fitted mm light curve, with the black curve corresponding to the model and the red curve corresponding to the data; and d) contribution of individual source components to mm light curve. |
Open with DEXTER |
![]() |
Figure 11: Models A and B of the 26 May 2008 NIR flare: a) fitted NIR light curve, with the black curve corresponding to the model and the red curve corresponding to the data; b) contribution of different source components to NIR light curve; and c) Contribution of different source components to mm light curve. |
Open with DEXTER |
A correlation between the L'-band data at 10 UT and the radio data cannot be excluded, but the L'-band data light curve is not long enough to justify including this data in a full modelling of the NIR/mm-light curves in the framework of the adiabatic expansion model. However, high NIR flux density levels precede a positive flux density excursion in the mm-domain, as it is expected to in the framework of the adiabatic expansion model. Although the L'-band flux density is higher than that at other times, the event is not linked to a particularly strong flux density excursion in the mm-domain. However, the high L'-band flux may also be caused by a transient steeper infrared spectral index. This may indicate either a significantly higher turnover frequency (extrapolating into the sub-mm domain) or the presence of stronger synchrotron losses. Both effects may explain a lack of strong mm-emission following the NIR event at 10 UT.
3.3.2 Expansion modelling of the 17 May flare emission
The 17 May 2007 flare is constrained by the NIR and CARMA mm-data
(Fig. 9).
This time lag of h
between the flare events mainly determines the expansion velocity. A
constant flux density of 0.003 Jy was subtracted from the NIR
data to model the individual flare events. We present two modelling
approaches that lead to a successful description of the flare. In model
A, we used only 4 components describing 4 sub-flares
through
.
Motivated by the relativistic disk modelling referred to in
Sect. 3.1,
we used in model B a background component
located close to the centre of the overall NIR flare. We then modelled
the 4 sub-flares with additional components. In comparison to model A,
this
results in lower values for
and the spectral index
,
smaller source component sizes R0,
and higher cutoff frequencies
for components
through
.
The quality of the fit is comparable for models A
and B.
3.3.3 Expansion modelling of the 19 May flare emission
For the 19 May 2007 flare, the slower adiabatic expansion velocity of
0.007c allows us to fit the first event (centred on
the NIR at 5:30 UT)
with a single component (Fig. 10).
For higher velocities, a larger number of source components has to be
used to model the NIR data.
For both velocities, the second, later flare component (centred on the
NIR at
8:20 UT)
is constrained by the 3 mm CARMA data. For the model
parameters used to match the NIR data, the
predicted 3 mm flux densities of both flare components (at
5:30 UT
and
8:20 UT)
are at or below the flux density limit of 0.1-0.2 Jy, above
which no millimeter wavelength
flare event has been detected.
3.3.4 Expansion modelling of the 26 May flare emission
Modelling the 26 May 2008 flare (Fig. 11) requires a satisfactory description of an initially bright flare with sharply rising and falling flanks followed by a 1.5 h plateau. This can only be achieved with several source components each covering a maximum portion of time of about a 30 min duration. A smaller number of source components would require larger source component sizes resulting in higher magnetic fields and a violation of the 100 GHz flux limit. Given these difficulties, the number of source components is of course not unique and the models proposed in Table 3 can only be taken as an example.
The overall flare structure is quite reminiscent of the flare that was observed on 3 June 2008 simultaneously using APEX and the VLT, as described by Eckart et al. (2008b), with the exception that the 26 May flare is more continuous than the 3 June 2008 flare. An overall flare structure that starts out with a bright event followed by a single or a few fainter events was also reported by Eckart et al. (2009,2006a). These events may relate to a common physical structure. An explanation could be a disk structure that is expanded by differential rotation within the disk. Alternative explanations may involve a common history in terms of either accretion or magnetic field instability - possibly also within a jet structure (see Eckart et al. 2008a,b and Sect. 4).
4 Alternative short jet model
Emission from a jet and an underlying accretion process may also account for the spectrum of Sgr A*, with the plasma in the jet starting to become optically thin at longer wavelengths with increasing distance from the centre. Thus, radiation at different radio wavelengths probe different sections of the jet, leading to a correlation with emission in other wavelength regimes. Very close to the SMBH, it is difficult to distinguish between emission from the jet and the accretion flow. For this jet model, flux variations may then be caused by accretion processes or jet instabilities rather than be the result of modulation from an orbiting spot, followed by an adiabatic expansion of jet components. Details of a compact, weak jet structure are discussed in Markoff et al. (2007; also Markoff et al. 2005).
5 Summary
We have presented results from global coordinated multiwavelength observational campaigns carried out in 2007 and 2008 using NACO at VLT in the NIR K- and L'-bands, and CARMA, ATCA, and IRAM in the mm regime.We have presented a new method to obtain concatenated light curves of the compact mm-source Sgr A* from single dish telescopes and interferometers in the presence of significant flux density contributions from an extended and only partially resolved source. The method requires several consecutive datasets to be available for each participating observatory. From these data, we have evidence of four different flaring events in the NIR, with three of these events being covered later in the mm regime.
We have modelled the flare emission with a model involving
adiabatic expansion of synchrotron source components, and obtained
spectral index values, expansion velocities, and other parameters that
are consistent with previously published variability studies. The
mm-flares were found to occur on average about h after the NIR
flare emission, and the expansion velocities for the various flares
ranged from 0.005c-0.017c.
Modelling the NIR-flares with synchrotron components that
become
optically thick in the sub-mm wavelength regime and obey the flux
density values or limits in the mm- and X-ray domain
constrains their turnover frequency and flux density.
For the described global experiment, the low mm-regime flux density
limits imply that the turnover frequencies are not lower than about
THz
and that the optically thick peak flux densities
at or below this turnover frequencies do not exceed about
1 Jy
on average. However, in other experiments the source showed higher flux
density
flare emission in the mm-domain (e.g., 0.5-1.5 Jy in the
sub-mm;
Yusef-Zadeh et al. 2009; Marrone et al. 2008; Eckart et al. 2008,
2006b)
Further monitoring of the variability of Sgr A* in different wavelength regimes (including X-ray, NIR, sub-mm, and mm radio) and in polarised NIR/radio emission is required to improve our understanding of the adiabatic expansion model and improve statistics describing the flaring activity. Current mm-interferometers such as CARMA, PdBI, and ATCA and future telescopes such as ALMA which are able to distinguish between emission from Sgr A* and the thermal emission of the CND and mini-spiral surrounding the SMBH, and NIR telescopes with large apertures such as VLT, Keck, and LBT, which can distinguish Sgr A* from surrounding stars, may be combined to provide us with high quality data to study the evolution of synchrotron components in greater detail. With future mm-VLBI measurements at frequencies of 230 GHz and above, imaging of the central region of Sgr A* may become possible, enabling a deeper study of the accretion physics and testing of current theories of emission in the region (Doeleman et al. 2009,2008).
AcknowledgementsPart of this work was supported by the German Deutsche Forschungsgemeinschaft, DFG via grant SFB 494. M. Zamaninasab, D. Kunneriath, and R.-S. Lu, are members of the International Max Planck Research School (IMPRS) for Astronomy and Astrophysics at the MPIfR and the Universities of Bonn and Cologne. R.S. acknowledges support by the Ramón y Cajal programme by the Ministerio de Ciencia y Innovación of the government of Spain. Macarena Garcia-Marin is supported by the German federal department for education and research (BMBF) under the project numbers: 50OS0502 & 50OS0801. N. Sabha is a member of the Bonn Cologne Graduate School (BCGS) of Physics and Astronomy.
Appendix A: NIR light curves
The individual NIR K and L'-band
light curves from Fig. 1
are shown in Figs. A.1
to A.5.
A spectral index of -0.9 (as we preferentially used in our modelling,
see Table 3
and Eckart et al. 2008b, 2009; Yusef-Zadeh et al. 2008) gives us a
scaling factor of 0.6, which we use to scale the L'-band
data, and then combine the K and L'-band
data as shown in Fig. 1
and Figs. A.1
to A.5.
In the remaining figures, no scaling has been applied. The
uncertainties in flux measurements were obtained from the reference
star S2, which has known flux densities of mJy and
mJy
in the K and L'-bands,
respectively.
![]() |
Figure A.1: Observing time-K-band: 05:29:55 to 08:31:48, L'-band: 10:05:48 to 10:26:45. |
Open with DEXTER |
![]() |
Figure A.2: Observing time-K-band: 04:47:22 to 07:54:41, L'-band: 08:34:27 to 10:41:46. |
Open with DEXTER |
![]() |
Figure A.3: Observing time-K-band: 04:42:14 to 09:34:40, L'-band: 10:16:24 to 10:27:06. |
Open with DEXTER |
![]() |
Figure A.4: Observing time-L'-band: 06:03:26 to 10:26:00. |
Open with DEXTER |
![]() |
Figure A.5: Observing time-K-band: 04:55:00 to 09:28:22. |
Open with DEXTER |
References
- Armitage, P. J., & Reynolds, C. S. 2003, MNRAS, 341, 1041 [NASA ADS] [CrossRef] [Google Scholar]
- Bock, D., Bolatto, A. D., Hawkins, D. W., et al. 2006, in SPIE Conf. Ser., 6267 [Google Scholar]
- Bower, G. C., Falcke, H., Sault, R. J., & Backer, D. C. 2002, ApJ, 571, 843 [NASA ADS] [CrossRef] [Google Scholar]
- Christopher, M. H., Scoville, N. Z., Stolovy, S. R., & Yun, M. S. 2005, ApJ, 622, 346 [NASA ADS] [CrossRef] [Google Scholar]
- Diolaiti, E., Bendinelli, O., Bonaccini, D., et al. 2000, in SPIE Conf. Ser. 4007, ed. P. L. Wizinowich, 879 [Google Scholar]
- Do, T., Ghez, A. M., Morris, M. R., et al. 2008, J. Phys. Conf. Ser., 131, 012003 [NASA ADS] [CrossRef] [Google Scholar]
- Do, T., Ghez, A. M., Morris, M. R., et al. 2009, ApJ, 691, 1021 [NASA ADS] [CrossRef] [Google Scholar]
- Doeleman, S. S., Weintroub, J., Rogers, A. E. E., et al. 2008, Nature, 455, 78 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Doeleman, S. S., Fish, V. L., Broderick, A. E., Loeb, A., & Rogers, A. E. E. 2009, ApJ, 695, 59 [NASA ADS] [CrossRef] [Google Scholar]
- Dovciak, M., Karas, V., & Yaqoob, T. 2004, ApJS, 153, 205 [NASA ADS] [CrossRef] [Google Scholar]
- Eckart, A., & Genzel, R. 1996, Nature, 383, 415 [NASA ADS] [CrossRef] [Google Scholar]
- Eckart, A., Genzel, R., Ott, T., & Schödel, R. 2002, MNRAS, 331, 917 [NASA ADS] [CrossRef] [Google Scholar]
- Eckart, A., Baganoff, F. K., Morris, M., et al. 2004, A&A, 427, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Eckart, A., Baganoff, F. K., Schödel, R., et al. 2006a, A&A, 450, 535 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Eckart, A., Schödel, R., Meyer, L., et al. 2006b, A&A, 455, 1 [Google Scholar]
- Eckart, A., Baganoff, F. K., Zamaninasab, M., et al. 2008a, A&A, 479, 625 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Eckart, A., Schödel, R., García-Marín, M., et al. 2008b, A&A, 492, 337 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Eckart, A., Baganoff, F. K., Morris, M. R., et al. 2009, A&A, 500, 935 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Eisenhauer, F., Schödel, R., Genzel, R., et al. 2003, ApJ, 597, L121 [NASA ADS] [CrossRef] [Google Scholar]
- Eisenhauer, F., Genzel, R., Alexander, T., et al. 2005, ApJ, 628, 246 [NASA ADS] [CrossRef] [Google Scholar]
- Genzel, R., Eckart, A., Ott, T., & Eisenhauer, F. 1997, MNRAS, 291, 219 [NASA ADS] [CrossRef] [Google Scholar]
- Genzel, R., Pichon, C., Eckart, A., Gerhard, O. E., & Ott, T. 2000, MNRAS, 317, 348 [NASA ADS] [CrossRef] [Google Scholar]
- Ghez, A. M., Klein, B. L., Morris, M., & Becklin, E. E. 1998, ApJ, 509, 678 [NASA ADS] [CrossRef] [Google Scholar]
- Ghez, A. M., Morris, M., Becklin, E. E., Tanner, A., & Kremenek, T. 2000, Nature, 407, 349 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Ghez, A. M., Hornstein, S. D., Bouchez, A., et al. 2004a, in BAAS, 36, 1384 [Google Scholar]
- Ghez, A. M., Wright, S. A., Matthews, K., et al. 2004b, ApJ, 601, L159 [CrossRef] [Google Scholar]
- Ghez, A. M., Salim, S., Hornstein, S. D., et al. 2005, ApJ, 620, 744 [NASA ADS] [CrossRef] [Google Scholar]
- Ghez, A., Morris, M., Lu, J., et al. 2009, in Astronomy, AGB Stars and Related Phenomenastro2010: The Astronomy and Astrophysics Decadal Survey, 2010, 89 [Google Scholar]
- Gillessen, S., Eisenhauer, F., Trippe, S., et al. 2009, ApJ, 692, 1075 [NASA ADS] [CrossRef] [Google Scholar]
- Gould, R. J. 1979, A&A, 76, 306 [NASA ADS] [Google Scholar]
- Guesten, R., Genzel, R., Wright, M. C. H., et al. 1987, ApJ, 318, 124 [NASA ADS] [CrossRef] [Google Scholar]
- Herrnstein, R. M., Zhao, J., Bower, G. C., & Goss, W. M. 2004, AJ, 127, 3399 [NASA ADS] [CrossRef] [Google Scholar]
- Kunneriath, D., Eckart, A., Vogel, S., et al. 2008, J. Phys. Conf. Ser., 131, 012006 [NASA ADS] [CrossRef] [Google Scholar]
- Lenzen, R., Hartung, M., Brandner, W., et al. 2003, in SPIE Conf. Ser. 4841, ed. M. Iye, & A. F. M. Moorwood, 944 [Google Scholar]
- Lu, R., Krichbaum, T. P., Eckart, A., et al. 2008, J. Phys. Conf. Ser., 131, 012059 [NASA ADS] [CrossRef] [Google Scholar]
- Lu, R., Krichbaum, T. P., Eckart, A., et al. 2009, A&A, submitted [Google Scholar]
- Markoff, S., Falcke, H., Yuan, F., & Biermann, P. L. 2001, A&A, 379, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Markoff, S., Nowak, M. A., & Wilms, J. 2005, ApJ, 635, 1203 [NASA ADS] [CrossRef] [Google Scholar]
- Markoff, S., Bower, G. C., & Falcke, H. 2007, MNRAS, 379, 1519 [NASA ADS] [CrossRef] [Google Scholar]
- Marrone, D. P., Baganoff, F. K., Morris, M. R., et al. 2008, ApJ, 682, 373 [NASA ADS] [CrossRef] [Google Scholar]
- Marscher, A. P. 1983, ApJ, 264, 296 [NASA ADS] [CrossRef] [Google Scholar]
- Mauerhan, J. C., Morris, M., Walter, F., & Baganoff, F. K. 2005, ApJ, 623, L25 [NASA ADS] [CrossRef] [Google Scholar]
- Meyer, L., Eckart, A., Schödel, R., et al. 2006a, A&A, 460, 15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Meyer, L., Schödel, R., Eckart, A., et al. 2006b, A&A, 458, L25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Meyer, L., Schödel, R., Eckart, A., et al. 2007, A&A, 473, 707 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pechácek, T., Karas, V., & Czerny, B. 2008, A&A, 487, 815 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rousset, G., Lacombe, F., Puget, P., et al. 2003, in SPIE Conf. Ser. 4839, ed. P. L. Wizinowich, & D. Bonaccini, 140 [Google Scholar]
- Sabha, N., Witzel, G., Eckart, A., et al. 2010, A&A, 512, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Schödel, R., Ott, T., Genzel, R., et al. 2002, Nature, 419, 694 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Schödel, R., Ott, T., Genzel, R., et al. 2003, ApJ, 596, 1015 [NASA ADS] [CrossRef] [Google Scholar]
- van der Laan, H. 1966, Nature, 211, 1131 [NASA ADS] [CrossRef] [Google Scholar]
- Yusef-Zadeh, F., Bushouse, H., Dowell, C. D., et al. 2006a, ApJ, 644, 198 [NASA ADS] [CrossRef] [Google Scholar]
- Yusef-Zadeh, F., Roberts, D., Wardle, M., Heinke, C. O., & Bower, G. C. 2006b, ApJ, 650, 189 [NASA ADS] [CrossRef] [Google Scholar]
- Yusef-Zadeh, F., Wardle, M., Heinke, C., et al. 2008, ApJ, 682, 361 [NASA ADS] [CrossRef] [Google Scholar]
- Yusef-Zadeh, F., Bushouse, H., Wardle, M., et al. 2009, ApJ, 706, 348 [NASA ADS] [CrossRef] [Google Scholar]
- Zamaninasab, M., Eckart, A., Kunneriath, D., et al. 2008a, Mem. Soc. Astron. It., 79, 1054 [NASA ADS] [Google Scholar]
- Zamaninasab, M., Eckart, A., Meyer, L., et al. 2008b, J. Phys. Conf. Ser., 131, 012008 [NASA ADS] [CrossRef] [Google Scholar]
- Zamaninasab, M., Eckart, A., Witzel, G., et al. 2010, A&A, 510, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Zhao, J., Young, K. H., Herrnstein, R. M., et al. 2003, ApJ, 586, L29 [NASA ADS] [CrossRef] [Google Scholar]
- Zhao, J., Herrnstein, R. M., Bower, G. C., Goss, W. M., & Liu, S. M. 2004, ApJ, 603, L85 [NASA ADS] [CrossRef] [Google Scholar]
Footnotes
- ... CARMA
- Support for CARMA construction was derived from the states of California, Illinois, and Maryland, the Gordon and Betty Moore Foundation, the Kenneth T. and Eileen L. Norris Foundation, the Associates of the California Institute of Technology, and the National Science Foundation. Ongoing CARMA development and operations are supported by the National Science Foundation under a cooperative agreement, and by the CARMA partner universities.
- ... ATCA
- ATCA is operated by the Australia Telescope National Facility, a division of CSIRO, which also includes the ATNF Headquarters at Marsfield in Sydney, the Parkes Observatory and the Mopra Observatory near Coonabarabran.
- ... IRAM
- The IRAM 30 m millimeter telescope is operated by the Institute for Radioastronomy at millimeter wavelengths - Granada, Spain, and Grenoble, France.
All Tables
Table 1: Log of the mm- and sub-mm observations.
Table 2: Log of the near-infrared observations.
Table 3: Source component parameters for the combined SSC and adiabatic expansion model of the 15, 17, and 19 May 2007 and the 26 May 2008 flares.
All Figures
![]() |
Figure 1: Combined differential light curve of Sgr A* in the mm/sub-mm domain for the May 2007 observing run. The MAMBO 2 bolometer at the IRAM 30 m-telescope was operated at a central wavelength of 1.2 mm (250 GHz). The CARMA data were centred on 100 GHz and the ATCA data on 86 GHz. We also show the daily 7 mm and 3 mm flux density averages from simultaneous VLBA observations (Lu et al. 2008,2009). In the lower part of the plot, we show the combined NIR light curve consisting of K-band and L'-band light curves (see Appendix for individual light curves). The time axis is labelled with UT hours starting at 00 h on May 15. |
Open with DEXTER | |
In the text |
![]() |
Figure 2: Difference maps at 3 mm of a 40'' diameter region (shown by the dashed white line) centred on Sgr A*, obtained from the difference between full synthesis maps of the individual days of CARMA observations and the full CARMA data set as described in the text. The maps have been plotted with the same color-coding table. The figure shows that the flux density variations that are evident from the differential light curves (see Fig. 1) can also be seen in the maps constructed from the corresponding data. |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
The 3 mm CARMA and ATCA (top red and blue, respectively) and K-
and L'-band NIR (middle red and blue, respectively)
light curves of the 17 May 2007 flare, with the cross-correlation
(bottom) between the two showing that the mm flare follows the NIR
flare with a time lag of |
Open with DEXTER | |
In the text |
![]() |
Figure 4: 3 mm CARMA Sgr A* ( top) and calibrator (1733-130) ( middle) light curves from 17 May 2007, with the cross-correlation ( bottom) between the two showing zero correlation between the calibrator and the source fluxes. |
Open with DEXTER | |
In the text |
![]() |
Figure 5:
3.8 |
Open with DEXTER | |
In the text |
![]() |
Figure 6:
Simulated images of multiple spots revolving around a black hole,
projected on the observer image plane ( |
Open with DEXTER | |
In the text |
![]() |
Figure 7: Light curve obtained from the multiple spot simulations shown in Fig. 6, the solid black line representing the overall light curve and the different coloured dashed lines showing the contributions of individual spots. |
Open with DEXTER | |
In the text |
![]() |
Figure 8: Models A and B of the 15 May 2007 NIR flare: a) fitted NIR light curve, with the black curve corresponding to the model and the red curve corresponding to the data, b) contribution of individual source components to NIR light curve, and c) contribution of individual source components to mm light curve. |
Open with DEXTER | |
In the text |
![]() |
Figure 9: Models A and B of the 17 May 2007 NIR flare: a) fitted NIR light curve, with the black curve corresponding to the model and the red curve corresponding to the data; b) contribution of individual source components to NIR light curve; c) fitted mm light curve, with the black curve corresponding to the model and the red curve corresponding to the data; and d) contribution of individual source components to mm light curve. |
Open with DEXTER | |
In the text |
![]() |
Figure 10: Models A and B of the 19 May 2007 NIR flare: a) fitted NIR light curve, with the black curve corresponding to the model and the red curve corresponding to the data; b) contribution of individual source components to NIR light curve; c) fitted mm light curve, with the black curve corresponding to the model and the red curve corresponding to the data; and d) contribution of individual source components to mm light curve. |
Open with DEXTER | |
In the text |
![]() |
Figure 11: Models A and B of the 26 May 2008 NIR flare: a) fitted NIR light curve, with the black curve corresponding to the model and the red curve corresponding to the data; b) contribution of different source components to NIR light curve; and c) Contribution of different source components to mm light curve. |
Open with DEXTER | |
In the text |
![]() |
Figure A.1: Observing time-K-band: 05:29:55 to 08:31:48, L'-band: 10:05:48 to 10:26:45. |
Open with DEXTER | |
In the text |
![]() |
Figure A.2: Observing time-K-band: 04:47:22 to 07:54:41, L'-band: 08:34:27 to 10:41:46. |
Open with DEXTER | |
In the text |
![]() |
Figure A.3: Observing time-K-band: 04:42:14 to 09:34:40, L'-band: 10:16:24 to 10:27:06. |
Open with DEXTER | |
In the text |
![]() |
Figure A.4: Observing time-L'-band: 06:03:26 to 10:26:00. |
Open with DEXTER | |
In the text |
![]() |
Figure A.5: Observing time-K-band: 04:55:00 to 09:28:22. |
Open with DEXTER | |
In the text |
Copyright ESO 2010
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.