A&A 439, L9-L13 (2005)
Y. Clénet1 - D. Rouan1 - D. Gratadour1 - O. Marco2 - P. Léna1 - N. Ageorges2 - E. Gendron1
1 - Observatoire de Paris, LESIA, 5 place Jules Janssen, 92195 Meudon Cedex, France
2 - European Southern Observatory (ESO), Alonso de Cordova 3107, Vitacura, Casilla 19001, Santiago 19, Chile
Received 28 April 2005 / Accepted 26 June 2005
We have collected in 2004 adaptive optics corrected images of the Galactic Center region with NAOS-CONICA at VLT. A strong variability was observed as well as a correlation between the photocenter positions and fluxes of the -band counterpart of Sgr A*. It is interpreted as the combined emission of a point-like flaring source at the position of Sgr A*/IR itself and an extended dust structure, 75 mas south west of Sgr A*/IR, which we name Sgr A*-f. We examine the different possible mechanisms to explain this dual Sgr A* emission and conclude it is likely a flaring emission associated to energetic events in the close environment of the black hole plus a quiescent emission resulting from the collision of Sgr A*-f by a jet from Sgr A*.
Key words: Galaxy: center - infrared: stars - instrumentation: adaptive optics
We report the detection of the -band counterpart of Sgr A* (hereafter Sgr A*/) as a source whose intensity and position are time dependent. It results from the combined emission of Sgr A*/IR itself and a nearby extended structure we name Sgr A*-f. We first detail the observations and data reduction, next explain the photometry and astrometry methods, then examine the position-flux correlation of Sgr A*/ and finally explore different mechanisms for this Sgr A*/ emission.
On-source (A) and on-sky (B) 10241024 frames, each one resulting from the mean of 60 subintegrations of 0.175 s, have been alternatively acquired following an ABBA cycle, with a random jitter within a box 6 wide every other cycle. The sky position was 3 away from the on-source position. The total on-source integration times were 27.3 min, 16.8 min, 49 min, 10.5 min and 9.8 min for Nights 1, 2, 3, 4, 5 respectively.
After sky subtraction, flat fielding and bad pixel correction, the recentering has been done in two steps: (1) the relative offsets between individual images have been estimated by computing their cross-correlation function, leading to an accuracy of 0.5 pixel; (2) a subpixel recentering has been performed with a maximum likelihood based algorithm (Gratadour et al. 2005), leading to a final accuracy of 0.01 pixel i.e. 0.3 mas.
For each night, to allow for a study of the time evolution of the Sgr A*/ flux with a sufficient signal to noise, we have averaged the recentered individual images into groups of four, leading to a cube of averaged images with a typical sampling of 265 s and a 42 s integration time for each averaged image.
|Figure 1: Astrometric/photometric errors and their standard deviations. Crosses are for Night1, stars for Night2, diamonds for Night3, triangles for Night4, and squares for Night5. Errors are estimated for each magnitude bin of the luminosity function computed by StarFinder. The translation to get the absolute photometry can be done with the zero-point of each night (4.5). To avoid confusion, the different symbols are horizontally shifted with respect to the bin centers.|
|Open with DEXTER|
The astrometry and the photometry of all averaged images have been obtained with StarFinder (Diolaiti et al. 2000), a PSF-fitting software designed for AO data.
To estimate the StarFinder accuracy for a given averaged image, we have added, in each magnitude bin of the luminosity function computed by StarFinder, 20% of synthetic stars at random positions and with random magnitudes and re-analyzed this new image with StarFinder. For each night, this procedure has been applied to three different averaged images.
Up to the relative magnitude bin [18, 19] ( [13.5, 14.5]), the error dispersion between the different nights is small, as well as the standard deviation associated to each night error value (Fig. 1). This is no longer the case for the last bin [19, 20] ( [14.5, 15.5]), most probably because of the lower accuracy of StarFinder for these fainter stars.
The absolute photometry has been computed as in Clénet et al. (2004b), using the photometry of the non variable stars IRS 16C, IRS 29N and IRS 33SE published in Blum et al. (1996). Given the accuracy of this zero-point computation (0.08 mag) and the StarFinder accuracy (Fig. 1), we finally estimate the photometric error for the 14 stars to be 0.18 mag for Night1, 0.17 mag for Night2, and 0.16 mag for Night3, Night4 and Night5.
Dereddened fluxes have been computed as in Clénet et al. (2004b) with an -band extinction and a zero magnitude value F0() = 248 Jy.
Following Clénet et al. (2004b), the astrometry has been calculated with respect to Sgr A* taken as the dynamical center (DC) of the S2 star orbit. The offsets between S2 and DC are given in Table 1 and the final astrometric errors in Table 2.
Table 1: Offsets between S2 and the dynamical center Sgr A*, computed with the S2 orbit parameters of Eisenhauer et al. (2003). Errors are computed by simulating several S2 orbits with gaussian distributions of the orbit parameters (cf. Clénet et al. 2004b).
Table 2: Astrometric errors for the 14 stars, given the DC positioning error (Table 1) and the StarFinder accuracy (Fig. 1).
|Figure 2: Upper line: Sgr A*/ (continuous line/crosses) and S2 (dash line/filled circles) light curves. Second line: distance in mas between Sgr A*/L and DC. Third line: opposite of the Sgr A*/ polar angle with respect to DC. The curves of the different nights are separated by a vertical line. All curves have been smoothed with a boxcar average of 2. Bottom line: 1.4 1.4 Night 3 images of the central stellar cluster during Sgr A*/ flare (couple of images on the left) and quiescent (couple of images on the right) states. Each couple of images is made of a 105 s integration time reduced image and the deconvolved corresponding one, obtained with 40 iterations of the maximum likelihood deconvolution IDL procedure. The deconvolved images are presented here only to highlight the sources position/extension and have not been used elsewhere for any quantitative analysis.|
|Open with DEXTER|
|Figure 3: Left: distance of Sgr A*/ from DC in mas as a function of the Sgr A*/ magnitude. Crosses, stars, diamonds, triangles, squares and filled circles are Night1, 2, 3/part1, 3/part2, 4, 5 data points respectively. Right: Sgr A*/ positions with respect to DC. Crosses represent positions for (Sgr A*/) < 12.6, diamonds for (Sgr A*/) < 13.0, stars for (Sgr A*/) < 13.4, triangles for (Sgr A*/) < 13.8.|
|Open with DEXTER|
The photocenter positions and fluxes of Sgr A*/ were measured on each averaged image using StarFinder. Firstly, as described in Clénet et al. (2004b), Sgr A*/ has exhibited variability on time scales of about one hour: flares with an amplification factor 2 during Night 1 and Night 3/part2; only decays or rises, with an amplification factor 3, for the other nights.
Secondly, we confirm the offset of Sgr A*/ in quiescence with respect to DC reported in Clénet et al. (2004b) at and ; we also observe a clear correlation between the Sgr A*/ flux and a) the DC-Sgr A*/ distance, b) the Sgr A*/ position angle (PA) with respect to DC. We show that these correlations are reproducible from one night to another: Sgr A*/ and DC are all the more distant and Sgr A*/ is all the more shifting toward DC than Sgr A*/ is bright (Figs. 2, 3).
Hence, our observations point at a dual mechanism responsible for the 3-5 m emission in Sgr A*: one traced by the flaring state, unresolved and associated to Sgr A*/IR itself; another one associated to the quiet state, clearly resolved with a typical extension of about 130 mas, occurring at the south-west of the dynamical centre (PA -130) and at a distance of 75 mas typically (600 AU) from it. A simple argument on the light velocity limit comforts the likeliness of this dual mechanism from two distinct sources where the emission originates from two distinct sources that NACO has marginally separated spatially, thanks to the extension of its quiescent part.
The flaring source compactness and the typical variation times measured at do not conflict with the interpretation of Eckart et al. (2004): both near-IR and X- emissions would be related to synchrotron self-Compton emission of electrons. The flare could be mainly due to Doppler boosted emission from material moving toward the observer, either matter at the base of a jet, or matter bound to the black hole on one of the innermost stable orbits. No simultaneous observations at several IR wavelengths have yet been done to allow one to derive the IR spectral index and we can only guess that flares are rather reproducible phenomena, as indicated by the similar behaviour of rises and decays in the observed -band events.
Though, another process should be invoked to explain a quiescent emission as far as 600 AU, corresponding to about 104 Schwarzschild radius (). A stellar emission could hardly account for this quiescent emission since: a) no K-band source can be observed at this location in the images of Eckart et al. (2004, see their Fig. 2); b) the inferred color (larger than 3) would be too large; c) the source appears extended.
An estimate of the IR color temperature
corresponding to the Sgr A*/ quiescent emission can be computed by assuming the quiescent K-band dereddened flux density found by Eckart et al. (2004) ( = 1.9 mJy), a quiescent -band dereddened flux density of 3.0 mJy and by resolving the optically thin case equation:
To estimate a maximum equilibrium temperature
of dust grains heated by Sgr A* itself, we equal the powers received and radiated by the grain:
We have also explored the possibility that the quiescent extended emission could be due to synchrotron emission from a jet. The absorption coefficient of the synchrotron emission is , where p and N0 describe the power-law dependence of the electrons energy distribution ( ) and B is the magnetic field. Using the jet model described in Falcke & Markoff (2000), completed by the analytical expression of the quantity established in Falcke (1996), we have set the parameters at the jet nozzle according to Yuan et al. (2002), i.e. so that to fit the millimetric observations. As expected, the resulting optical depth at a typical distance of 600 AU from the black hole is extremely low, so that the quiescent extended emission cannot be explained by the high frequency wing of synchrotron radiation from a jet.
If we now consider the mass loss rate in the jet found by Yuan et al. (2002), yr-1, and translate it in kinetic power delivered by the jet (assuming an equipartition between ions and electrons), we get , where , , are the mass of electron, the mass of proton and the electron to light speed ratio, respectively. With (Yuan et al. 2002), we obtain erg/s = 32 .
On the other hand, the total luminosity, if due to grains, that would correspond to the observed -band luminosity
In this description, the observed quiescent state extended source would be a small dust condensation, close to Sgr A* and heated through a colliding jet. The observed position-magnitude correlation would then result from Starfinder measuring the position/magnitude of the combined emission of Sgr A*/IR itself and this dust structure.
We note that the Sgr A*/IR naming should be reserved to the flaring location since the role of the black hole is obvious for this state. We propose to name Sgr A*-f (standing for Sgr A* "flake'') the extended emission located south west of Sgr A*/IR. Sgr A*/ is then the observed combined emission of Sgr A*/IR and Sgr A*-f.
We cannot totally exclude that Sgr A*-f belongs to the mini-spiral and is then located at a fraction of parsec from Sgr A* along the line of sight. Though, the physical association of Sgr A*/IR and Sgr A*-f is supported by the presence of a chain of gaseous clouds, probably associated to Sgr A* through a ridge of emission, observed in the radio continuum between Sgr A and the mini-cavity (Yusef-Zadeh et al. 1990). These clouds would be driven by fast winds towards the gas streamers orbiting Sgr A* (the "Minispiral") and create by shocks the mini-cavity.
|Figure 4: Encircled H II components positions on a 2.4 2.4 NACO image. From left to right: H II 5D, H II 5C, source , source (Yusef-Zadeh et al. 1998; Zhao & Goss 1999). The circle diameter is the 0.1 radio resolution. ID18 is the bright extended source just south to the H II 5D circle.|
|Open with DEXTER|
The origin of these radio features is still discussed: Wardle & Yusef-Zadeh (1992) and Melia et al. (1996) proposed they are formed by the gravitational collimation of the IRS 16cluster stellar winds by Sgr A*. Alternatively, Lutz et al. (1993) invoked a jet from Sgr A*, similarly to our scenario, as the wind source.
The position of these gaseous H II components, taking into account their estimated proper motions and their fairly extended sizes (Yusef-Zadeh et al. 1998; Zhao & Goss 1999), is very close to very red IR extended sources (Fig. 4), e.g. H II 5D in Zhao & Goss (1999) and ID18 in Clénet et al. (2004a). This suggests a common origin, if not an identification, between the sources at the different wavelengths. Sgr A*-f could then related to this chain of structures, possibly to the unnamed feature in Fig. 3 upper panel of Zhao & Goss (1999), located, in 1998, 0.35 south west to Sgr A*: a proper motion similar to the other H II structures (25 mas yr-1) could have shifted this feature very close to the 2004 position of Sgr A*-f.
In the future, to constrain the aforementioned physical model of a jet interacting with Sgr A*-f, simultaneous observations at different IR wavelengths of Sgr A*/IR in both quiescent and flare states will be essential. An instrument such as the NACO prism disperser, which covers simultaneously the IR spectrum from J (1.3 m) to M (4.8 m), could be efficiently dedicated to this task.