A&A 469, 993-1002 (2007)
DOI: 10.1051/0004-6361:20066265
K. Muzic1,2 - A. Eckart1,2 - R. Schödel1 - L. Meyer1 - 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
Received 18 August 2006 / Accepted 23 April 2007
Abstract
Context. L'-band (3.8 m) images of the Galactic center show a large number of thin filaments in the mini-spiral, located west of the mini-cavity and along the inner edge of the Northern Arm. One possible mechanism that could produce such structures is the interaction of a central wind with the mini-spiral. Additionally, we identify similar features that appear to be associated with stars.
Aims. We present the first proper motion measurements of the thin dust filaments observed in the central parsec around SgrA* and investigate possible mechanisms that could be responsible for the observed motions.
Methods. The observations have been carried out using the NACO adaptive optics system at the ESO VLT. The images have been transformed to a common coordinate system and features of interest were extracted. Then a cross-correlation technique could be performed in order to determine the offsets between the features with respect to their position in the reference epoch.
Results. We derive the proper motions of a number of filaments and 2 cometary shaped dusty sources close (in projection) to SgrA*. We show that the shape and the motion of the filaments does not agree with a purely Keplerian motion of the gas in the potential of the supermassive black hole at the position of SgrA*. Therefore, additional mechanisms must be responsible for their formation and motion. We argue that the properties of the filaments are probably related to an outflow from the disk of young mass-losing stars around SgrA*. In part, the outflow may originate from the black hole itself. We also present some evidence and theoretical considerations that the outflow may be collimated.
Key words: Galaxy: center - infrared: ISM - ISM: jets and outflows - accretion, accretion disks
At a distance of 8 kpc (Reid 1993; Eisenhauer et al. 2005,2003)
the center of the Milky Way
shows a broad variety of structures.
Near-infrared diffraction limited imaging over the past 10 years
(Schödel et al. 2003; Ghez et al. 1998,2005; Eckart & Genzel 1996; Eckart et al. 2002; Genzel et al. 1997,2000; Ghez et al. 2000; Genzel et al. 2003)
has resulted in convincing evidence for a
black hole being located at the center of a dense stellar cluster.
This finding is supported by the discovery of a variable X-ray and NIR source
at the position of SgrA* (Eckart et al. 2004; Baganoff et al. 2001,2003).
![]() |
Figure 1:
Left: L'-band image of the central part of the Milky Way (![]() ![]() ![]() |
Open with DEXTER |
High angular resolution spectroscopic observations have confirmed
the presence of a dense clumpy molecular ring (also called circum nuclear
disk, CND) of warm dust (Zylka et al. 1995) and neutral gas
(Güsten et al. 1987; Herrnstein & Ho 2002; Wright et al. 2001; Jackson et al. 1993)
surrounding a short-lived central cavity of much lower mean gas density.
The central cavity is about 3 pc in diameter (see Scoville et al. 2003; Rieke et al. 1989; Chan et al. 1997)
and contains
the so-called mini-spiral, mainly composed of ionized gas and dust. Vollmer & Duschl (2000)
have re-examined the mini-spiral gas dynamics traced by the
[Ne II] line emission
(
m; Lacy et al. 1991)
and the H92
radio emission data (Roberts & Goss 1993)
and discuss the kinematic structure of the inner
pc
of the Galaxy. Vollmer & Duschl (2000)
derive a three-dimensional kinematic model of gas streams that describes
the bulk gas motion in three different planes.
Most of the material is located in a main plane which
is connected to the inner edge of the CND.
Paumard et al. (2004) observed the mini-spiral with integral field spectroscopy at
2.06 m (He I) and 2.16
m (Br
)
covering roughly
around SgrA*.
They analyze the kinematics of the mini-spiral and describe nine distinguishable structures.
From their analysis, the Northern Arm consists of a weak, continuous,
triangular-shaped surface, drawn out into a narrow stream in the vicinity
of SgrA* where it shows a strong velocity gradient, and a bright western rim.
In the previously published observations of the mini-spiral, a number of narrow,
filamentary structures
have been reported. The observations include high-resolution
radio data (Zhao & Goss 1998; Yusef-Zadeh et al. 1998) and NIR
Pa
(Scoville et al. 2003),
Br
(Morris 2000) and He I (Paumard et al. 2001) emission line maps.
Here we present L'-band (3.8
m; see Fig. 1) observations of the Galactic Center
performed with the NACO adaptive optics system at the ESO VLT UT4
.
These data allow us not only to identify features as narrow as the
100 mas angular
resolution of the observations and in much
larger number than in any previous dataset, but also to measure them clearly with high precision.
Some of the filaments, particularly those
associated with the Northern Arm and one detected stellar bow shock, have earlier been reported
by Clénet et al. (2004) from observations at the same wavelength.
The L'-band flux from the mini-spiral is primarily coming from the
thermal dust emission, which is confirmed by the observations at mid-infrared wavelengths (Viehmann et al. 2006).
The detection of some of the filaments in Pa
and Br
tells us that the emission is as well associated with the ionized gas component of the ISM.
In the L'-band, there is probably no significant contribution of the emission from the light scattered
by dust grains. In this case one would expect to observe a significant
amount of the polarized emission also at the shorter wavelength K-band - which is not
observed (Eckart et al. 1995; Ott et al. 1999).
Furthermore, there is no bright NIR source at the Galactic Center that
could be identified as the source of large amounts of scattered light.
Polarized emission from the mini-spiral observed at longer infrared wavelengths is due
to thermal emission from elongated dust grains (Aitken et al. 1998,1991).
The central parsec of the Milky Way harbours a cluster of massive stars
(Clénet et al. 2001; Blum et al. 1995; Paumard et al. 2006; Genzel et al. 1996,2000; Eckart et al. 1999; Krabbe et al. 1995) that supply
-3
yr-1 of gas in the form
of stellar winds to the center (Najarro et al. 1997).
The amount of gas available for accretion onto SgrA* at the Bondi radius
is estimated to be
10-5
yr-1 (Quataert 2003; Baganoff et al. 2003).
The actual mass accreted onto SgrA* is up to 2 orders of magnitude lower, as determined
by the linear polarization measurements at 230 GHz (Bower et al. 2003; Marrone et al. 2006).
We can see that about 99% of the material from stellar winds
does not even get close to the Bondi radius and must therefore
escape the central arcseconds in form of a wind.
Blandford & Begelman (1999) point out that for accretion rates well below the Eddington rate most of
the gas accretes with positive energy and may escape. Therefore they propose a solution
(ADIOS; adiabatic inflow-outflow solution) in
which most of the gas is lost through winds and only a small fraction of the gas
is actually accreted onto the black hole.
This is also
in agreement with the RIAF theory (radiatively inefficient accretion flow; see e.g. Yuan 2006)
that postulates strong outflows from SgrA*.
The outflow may be partially collimated in form of a jet.
Jet and shock models
(see Markoff & Falcke 2003; Melia & Falcke 2001)
as well as accretion from an in-falling wind due to the surrounding cluster
of mass-losing hot stars (e.g. Coker et al. 1999; Yuan et al. 2002) are discussed
to explain the radio to X-ray properties of SgrA*.
There is some observational evidence for outflows in the central parsec as well.
The "mini-cavity'' region 3.5
southwest of the compact radio source SgrA*
was first pointed out on cm-radio maps by Yusef-Zadeh et al. (1990).
The strong Fe[III] line emission seen toward that region
(Lutz et al. 1993; Eckart et al. 1992) is consistent with gas excited by
a collision with a fast (
1000 km s-1) wind from a
source within the central few arcseconds
(Yusef-Zadeh & Wardle 1993; Yusef-Zadeh & Melia 1992).
Most recently Schödel et al. (2007) present an
extinction map of the central parsec.
In this map an extended region of low extinction,
centered on SgrA* runs in NE-SW direction across SgrA* and
is continued in the mini-cavity. It
supports the assumption of an outflow from the central arcseconds.
For one of the most prominent dust embedded sources, IRS 3, Viehmann et al. (2005)
find that the extended L- and M-band continuum emission has a characteristic bow-shock shape.
A possible explanation for its asymmetric appearance
is that IRS 3 consists of a mass-losing star
surrounded by a very thick, extended dust shell,
which is pushed northwest by a wind
from the direction of the IRS 16 cluster and SgrA*
(Pott et al. 2007).
Similarly to the north, the mass-losing envelope of IRS 7 appears to be
influenced by a strong central wind (Yusef-Zadeh & Melia 1992).
The narrow filaments observed in the L'-band are the main concern of this paper: we identify them, present proper motion measurements, discuss their nature and present an interpretation in the light of the postulated outflow from the center. In Sect. 2 we present the observations. In Sect. 3 we identify thin filaments and explain proper motion measurements in detail. In Sects. 4 and 5 we describe the proper motions and discuss the nature of the detected thin filaments. In order to validate our proper motion results we compare the K-band and L'-band stellar proper motions in Appendix A. In Sect. 6 we present the probable wind sources in the central few arcseconds. In Appendix C we present the outflow model and discuss the overall 3-dimensional orientation. The summary and conclusions are given in Sect. 7.
![]() |
Figure 2:
Identification of the thin features observed in the the L'-band. SgrA* is located at the center of the image.
We have marked the line of nodes for disk (i) at
a position angle of 28![]() |
Open with DEXTER |
The L' (3.8 m) and
-band (2.1
m)
images were taken with the NAOS/CONICA adaptive optics assisted imager/spectrometer
(Rousset et al. 1998; Lenzen et al. 1998; Brandner et al. 2002)
at the UT4 (Yepun) at the ESO VLT.
The L'-band data set includes images from 5 epochs
(2002.660, 2003.356, 2004.320, 2005.364 and 2006.408) with a resolution of
100 mas.
The
-band data set includes the images from epochs 2002.339, 2003.356, 2004.512, 2004.521 and
2006.413 with a resolution of
56 mas.
Data reduction (bad pixel correction, sky subtraction, flat field correction)
and formation of final mosaics was performed using the DPUSER software for
astronomical image analysis (T. Ott; see also Eckart & Duhoux 1990).
The absolute positions of sources in our AO images were derived by comparison to the VLA positions of IRS 10EE, 28, 9, 12N, 17, 7 and 15NE as given by Reid (1993). The radio positions and the positions in the L'-band image agree to within less than a single 27 mas pixel i.e. a quarter of the beam. However, our proper motion measurements were based on relative astrometry between frames in order to avoid errors that may arise from uncertainties of the absolute positions.
The L'-band images (Fig. 1, left) clearly show stars, bright dust embedded
sources as well as the emission of the mini-spiral.
In addition, one can distinguish a
number of previously unknown thin filaments.
In order to highlight these structures we produced high-pass filtered
maps (Fig. 1, right) by subtracting a smoothed version of the images from themselves,
and then smoothing the final images.
As a smoothing function in both cases we used a 4 pixel Gaussian (1 pixel corresponds
to 0.027
).
This procedure enhances all structures that have significant
power at spatial frequencies that correspond to the diffraction
limited beam size.
The thin filaments that are apparent in the L'-band maps are located
close to the mini-cavity, along the inner edge of the Northern Arm and,
in some cases, in the vicinity of stars.
The width of the filaments is 100 mas, the diffraction limit of NACO L'-band images.
Figure 2
shows a schematic representation of the thin features in the central parsec of the Milky Way as
derived from the data shown in the right panel of Fig. 1.
![]() |
Figure 3: High-pass filtered image with the narrow features labeled. |
Open with DEXTER |
Table 1: Proper motions of thin filaments.
According to their appearance we can divide the filaments in two classes: the first class
represents features most probably associated with stars (denominated with X); the second class
features are more elongated and associated with the mini-spiral. Those located west and south from
SgrA* are denominated SW and those associated with the Northern Arm and located east
from SgrA* we label NE.
They cannot be associated with stars in a straightforward way.
We find that the second class of filaments to the east and northeast
of SgrA* are located at the inner side of the Northern Arm and curved with their
convex sides eastward, whereas filaments to the west and southwest
of SgrA* are curved with their convex sides westward and elongated almost perpendicularly to the Bar.
Several of the more luminous filaments to the west and southwest are
located on the western side of the mini-cavity close (in projection) to
the IRS13/IRS2 complex or just south of it (most prominent are SW5, SW6, SW7, SW8).
One can also notice that some of these filaments are located
up to about 5
SSW of SgrA*, almost opposite to and approximately at the same
distance as the most distant filament in the Northern Arm.
Concerning the X features, we notice no preferred direction towards which they are curved.
We discuss their possible nature in Sect. 4.1.
In order to measure the positional offsets of a feature at different epochs, entire
image frames had to be transformed to a common coordinate system. This was done using IDL
image transformation routines. Transformation to second order was performed providing
correction of all translations, rotations and possible distortions between two images.
All the frames were transformed to the coordinate system of the 2003.356 epoch using
a total of 34 stars to calculate the transformation matrix. The stars used for
the transformation were chosen to be uniformly distributed across the field. The positions
were corrected for the stellar proper motions as derived from the -band images.
In order to validate our method and to demonstrate that reliable proper motions can be
obtained from the lower resolution L'-band data we compare
-band and L'-band stellar
proper motions in
Appendix A.
In order to measure proper motions of the observed thin filaments we calculated offsets at all epochs with respect to the reference epoch (2003.356). First, we extracted an image section containing the entire feature, or the part of the feature that is not directly contaminated by stars. The position of the center and dimensions of the image sections are given in Table 1. Then we re-binned those image sections by a factor of 10 (since all the offsets are on sub-pixel scale), masked residual stars if present in the frame, and finally subtracted the background. The resulting small frames were cross-correlated with a reference frame in order to calculate the offset of the feature with respect to the reference position. We computed the inverse of the sum of the squares of the difference between overlapping pixel values. This resulted in a 2D cross-correlation function which was then fitted with a 2D Gaussian to determine the exact position of the peak that represents the best overlap of two features. There are two main sources of error that occur during the determination of the offset. The first is caused by the uncertainties due to the cross-correlation method. In order to determine this uncertainty for each feature and epoch, we shift it by different known sub-pixel offsets in all directions, run the cross-correlation routine and then calculate the standard deviation obtained from the imposed and derived shifts. The second is resulting from the frame transformation procedure. From the initial sample of 34 stars we have randomly chosen 50 different sub-samples of 20 stars each. For all sub-samples we then repeated the transformation of all frames and calculated the positional uncertainty. The error bars in diagrams in Fig. 6 contain contributions of both sources of uncertainty.
![]() |
Figure 4:
Left: L'-band ( up) and ![]() ![]() |
Open with DEXTER |
It is important to note that the extraction of stars from our images via PSF fitting leaves significant residuals at the positions of the brightest stars. Thus, our approach was to extract sub-images containing the filaments from the original L'-band images if possible. Due to the presence of a large number of stars in L'-band images, some of the identified features could not be extracted at all to perform the cross-correlation (e.g. SW6, SW8, NE2). In some cases we were only able to extract one part of the filament, rather then the entire structure. Besides the proximity of stellar sources, there is another reason for this: many of the filaments are parts of more extended structures and could not be completely isolated. This means that the positional data along the feature do not result in a measurable proper motion. Therefore, the velocities of all NE and SW features, except that of the compact features SW7 and NE4, are measured approximately perpendicular to their extent. Some of the features are very faint in the L'-band and therefore unsuitable for the extraction. They can be more clearly identified in the smooth-subtracted images (features SW7, SW1, X4 and X6). In this case, only the SW7 appearance is sharp enough to give reliable results. On the other hand, some of the features are very bright (X5 and SW2) and may be associated with embedded stellar sources.
Using the Starfinder code for stellar field analysis (Diolaiti et al. 2000) we also produced low-pass images containing the background emission. The result of the point source extraction largely depends on the size of the box (in units of a FWHM of the PSF) used for the background estimate. For small values of the box size (<10), the background estimate is still contaminated by stars, and for values larger than 15 the compact features become diffuse. The brightest of the Northern Arm filaments, NE4, has a high enough S/N to be clearly detected in the background images. The background image is free from stars and therefore allows us to extract the NE4 feature without confusion. However, other NE and SW filaments cannot be clearly defined against the smooth background. Thus, only for the NE4 filament we additionally give the proper motion component along the extent of the feature (see also Appendix B). We performed the NE4 proper motion analysis for three different sizes of the box used for the background estimate (10, 13 and 15). The values given in Table 1 represent the average of the obtained results and the corresponding diagram in Fig. 6 shows only the result for the box with a width of 13 PSF FWHM values.
Along the line that connects SgrA* and the feature X3
(3.3
southwest from SgrA*) at the position (-0.56
, -0.55
)
from SgrA*,
one can notice feature X7, with a shape similar to that of X3.
The faint star at this position in the
-band image is represented
by a much more extended structure in the L'-band image (Fig. 4).
Both X3 and X7 are approximately aligned (in projection) with SgrA*.
The line of symmetry of the X7 "V''-shape is oriented
40
east of north and to
within less than 5
it is passing through the position of SgrA*.
The line of symmetry of X3 is slightly displaced from the line that connects it with SgrA* (
5
)
and is oriented
37
east of north. The opening angle toward SgrA* is
20
for X7 and
10
for X3.
The orientation of both sources towards SgrA* suggests an interaction with a possibly existing
outflow from the central region (see Sect. 6.1.2).
They may represent bow-shock-like features pointing back to the position of SgrA* as a common origin of such an outflow.
![]() |
Figure 5: L'-band image of the Galactic Center. Boxes mark thin filaments with measurable proper motions. Note that boxes in this image are different from those used for measurements (stated in Table 1). The arrows show proper motions of the thin filaments obtained in our study: light blue arrows stand for the features with measurable proper motion in both directions, while black arrows show only the proper motion component perpendicular to the feature (see the text for the explanation). The insignificant motion of the Northern Arm filament NE1 is marked with a circle rather than an arrow. The cross marks the position of SgrA*. |
Open with DEXTER |
![]() |
Figure 6:
Proper motions of thin filaments. The error bars show the 1![]() |
Open with DEXTER |
Table 1 summarizes the proper motions of thin filaments obtained in our study. The proper motions have been successfully determined in 11 regions. In order to transform the angular motions into velocities we assume a distance to SgrA* of 8 kpc.
Boxes in Fig. 5 mark the filaments that were analyzed.
The directions of motion in the plane of the sky are represented by arrows.
To emphasize that some of the motions shown in Fig. 5 do
not contain the complete information about the motion in the plane of the sky, we
represent them by different colors.
The light blue arrows in Fig. 5 show
the full proper motion of the corresponding feature, while the black ones represent
only the component perpendicular to the feature of interest.
A filled circle rather than an arrow at the position of the NE1 feature marks the
insignificant motion of this feature.
Note that for display purpose the dimensions of the boxes shown in Fig. 5 are
different from those used for the extraction of the feature (as they are stated in Table 1).
Diagrams in Fig. 6 contain detailed information about the results.
The velocity difference v in Table 1
represents the 1
uncertainty of the linear fit to the various epochs. Each offset is weighted by
the corresponding uncertainty. The escape velocities
were determined using the projected distance from SgrA* with a mass of
6
.
Feature X1 is associated with a star located (2.95
, 3.46
)
from SgrA*.
The proper motion of the associated star is
= (
) km s-1,
= (-46
5) km s-1
as derived from our L'-band data.
We find that the feature is moving approximately in the same direction as the star apparently
associated with it - which confirmes the bow-shock interpretation of Clénet et al. (2004).
Yusef-Zadeh et al. (1998) presented radio continuum observations at
cm and calculated
proper motions of ionized gas at the Galactic Center. Although no thin filaments similar to those
detected in NACO L'-band images are detected in the radio continuum images, some comparison can be done.
The filamentary ionized structure connecting IRS 13 and IRS 2 (Box 11 in Yusef-Zadeh et al. 1998)
coincides with the SW5 feature in our images (in this work divided in two parts, SW5a and SW5b,
because of the presence of a star). A westward direction of motion is obtained from
both measurements, but the value of this component obtained from the radio continuum is larger
(
= (-329
56) km s-1). The reason for this discrepancy is not clear. It is possible
that at different wavelengths we are actually probing different material.
Zhao & Goss (1998) also performed radio continuum observations at 7 and 13 mm and calculated
proper motions of 57 compact H II components. As in the previous case, the only comparable
structure is SW5 ("loop'' in their notation), but the results are not in good agreement with ours.
This feature, in the measurements of Zhao & Goss (1998),
is divided in several non-consistently moving components,
which is probably consequence of the granular substructure
seen in the radio image
(probably due to a low signal to noise ratio).
In our IR continuum images the structure is well defined.
Yusef-Zadeh et al. (1998) report an anomalous high-velocity feature, so-called "bullet'',
at the offset (-3.39
,
2.065
)
from SgrA* in the 1990 epoch. Approximately at this position
in the L'-band images one can notice a dusty feature (see Fig. 5). However,
given the high velocity quoted for the bullet it should have moved by
(
,
) = (0.30
, 0.34
)
to the northwest
from 1990 to 2006, which is the epoch of the image
in Fig. 5. No similar feature can be found at the corresponding
position in the 2006 epoch, or close to it.
Although the feature we see in the L'-band cannot be disentangled from the rest of the material
and therefore we are not able to measure its proper motion, a simple inspection
of images at different epochs shows that it does not seem to have a significant
velocity compared to the bullet. Therefore, we cannot identify this feature as the bullet.
The overall kinematic properties of the Northern Arm and Bar section
can be modeled as a single gas streamer orbiting SgrA*.
Such a model can be constrained via the radial velocity information from
the H92
(Roberts & Goss 1993) and [Ne II] line
(Lacy et al. 1991; Vollmer & Duschl 2000) as well as the overall distribution of the gas
and dust along these features.
Proper motions of ionized gas were in detail studied by Yusef-Zadeh et al. (1998).
They additionally give the radial velocities as measured by Roberts & Goss (1993) and
Roberts et al. (1996). According to Yusef-Zadeh et al. (1998) the predominant component
of the motion in the plane of the sky is from east to west. The flow of gas in the Northern Arm
follows its elongated shape: from the north-east with slightly redshifted
velocities it follows the orbital trajectory to
the southwest as it crosses the plane of the sky.
The ionized gas passes south of SgrA* before it moves to the northwest.
The complex velocity structure of the
mini-cavity, however, fails to be explained by such a simple model.
Similar results are reported by Paumard et al. (2004) based on He I and Br
observations. At the position of the NE1 filament the line emission is
slightly red shifted (40-60 km s-1) while at the positions of the other
Northern Arm filaments the radial velocity becomes negative. At the position
of NE3 the line emission is blue shifted by 40-60 km s-1 and close to NE4
by more than 130 km s-1.
The Bar material at the position of SgrA* and within few arcseconds south of it
has radial velocities close
to 0 km s-1. The emission becomes blue shifted as we move further to the west
(where SW3 and SW4 are located).
Paumard et al. (2004) analyzed the Northern Arm as a Keplerian system in the gravitational field of SgrA*. They fit a bundle of Keplerian orbits in order to model the full velocity field of the Northern Arm. The authors interpret the rim of the Northern Arm in terms of line-of-sight orbit crowding. The western edge facing SgrA* is the region that shows the highest density of crowded orbits (see Fig. 8 in Paumard et al. 2004). In order to construct Fig. 7 we used the part of the Fig. 8 in Paumard et al. (2004) having the same field of view as the L'-band images presented in this work.
![]() |
Figure 7:
The Pa ![]() |
Open with DEXTER |
The nature of the observed thin features, especially the SW and NE filaments, is unclear. One possible interpretation is that they are the consequence of the interaction of a central wind originating from the hot stars and/or possibly from SgrA*. While the origin and the mechanisms of such a wind will be discussed later in the text, here we concentrate on the effects it would have on the mini-spiral material in combination with the dominating gravitational potential of the central black hole.
We favor a scenario in which
the proposed wind compresses gas and dust in the mini-spiral and forms the observed filaments.
For the SW filaments this interaction could have taken place at the western edge of the
mini-cavity. The thin Northern Arm filaments may have been generated by a wind in the
opposite direction of the mini-cavity.
For the mini-spiral Vollmer & Duschl (2001) quote a sound speed of
km s-1.
This velocity is very small in comparison with the overall motion of
the Northern Arm (Paumard et al. 2004; Yusef-Zadeh et al. 1998).
Also, once formed, the filaments may still be
influenced by wind from the central few arcseconds
that may keep them confined.
The magnetic field of the Northern Arm, as reported by Aitken et al. (1991) and
Aitken et al. (1998),
could also play a role in confining ionized gas and charged dust grains along the
filaments.
If we assume a temperature of T=7000 K for the mini-spiral material
(Roberts & Goss 1993) then the average velocity of a hydrogen atom associated
with thermal energy is v = (3kT/m)
km s-1.
At some position in the Northern Arm it takes approximately
250 yr for it to cross the distance equivalent to the width of a filament
that we can resolve (
100 mas).
During this time the filament can drift
2.5
soutwards due to the motion of
the mini-spiral material in the central gravitational potential.
For dust particles it takes even longer to cross the
equivalent width of a filament due to their thermal motion.
This leads to the conclusion that a filament formed in the Northern Arm
of the mini-spiral will have a lifetime that is long enough to drift to the
southwest together with the gas stream.
Once formed, the dominant source for their
destruction will be the interaction with their surrounding turbulent
medium. Judged from the line widths of the H92
emission
(Roberts & Goss 1993) the velocity dispersion of the gas within the
mini-spiral is of the order of 50 km s-1.
Vollmer & Duschl (2000) assume a maximum value for the velocity
dispersion of 40% of the Keplerian velocity.
![]() |
Figure 8:
Results of our preliminary 3D dynamical model calculation.
Three particles (marked +, ![]() ![]() |
Open with DEXTER |
In order to investigate the kinematic behaviour of
linear features formed in a central potential
we carried out orbit simulations.
We assumed a dominating central potential
from a
6
black hole
at the position of SgrA*.
Schödel et al. (2007) estimate
that the gravitational field at the Galactic Center is dominated by the point mass of
the black hole inside a projected radius of about 6
.
Since all of the measured filaments
are either inside or just at the border of this circle, we do not take into account
the possible contribution of the
stellar cluster to the gravitational potential.
In agreement to similar calculations by
Vollmer & Duschl (2000)
and
Paumard et al. (2004)
we chose 3 test particles with velocities of
= 20 km s-1,
= -150 km s-1 and
= -100 km s-1launched from locations around
mpc = 5.6
,
mpc = 6.4
and
= 150 mpc
(see Fig. 8).
Here z is the distance behind the plane of the sky including SgrA*.
The corresponding orbits show that
elongated thin features "imprinted'' into the Northern Arm
are subject to a weak shear but to first
order keep their original structure while drifting to the south.
In our current picture Northern Arm
filaments drift south (while being stretched by less than a
factor of 2 within the Northern Arm)
until they reach the mini-cavity region where
they may be destroyed. Filaments created at the western edge of
the mini-spiral drift westward along the Bar.
From Fig. 7 one can clearly see that the Northern Arm filaments
cannot be simply explained as the consequence of orbit-crowding as suggested by Paumard et al. (2004).
A comparison with the velocity field obtained from their model shows a good agreement for the
perpendicular velocity component of the filaments NE1 and NE3 (within a 1 error for NE1 and
2
for NE3). This supports the assumption that the observed filaments are really
drifting southwards
together with the Northern Arm material.
However, none of the orbits actually follows the shape of the filaments and therefore
the proper motions of the gas as stated in Paumard et al. (2004) are not aligned with the extent
of the filaments. Thus one needs to include an additional
force that could form the filaments, keep them confined
and prevent them from being destroyed by the bulk gas flow.
This is especially valid for the filament NE3.
The exact position of the filament NE4 does not coincide with any of the Keplerian orbits
fitted by Paumard et al. (2004). They find that the Northern Arm is split into
two layers as it reaches the mini-cavity. The mini-cavity itself is a part of the first layer,
while the second layer is deflected northward from it and seems to contain the filament NE4.
Concerning the SW filaments, our study shows that the proper motion components perpendicular to their extent
are consistent with the motion in the dominating central potential,
but we are still facing the question of their formation.
As proposed earlier in this section, they could have been formed by an interaction at the western
edge of the mini-cavity, and then drift
westwards together with the rest of the mini-spiral material.
Concerning the possibility that the filaments are the consequence of an
interaction between a fast wind and the mini-spiral material, we can try to
predict the width of a shock using formula (8) in Hartigan et al. (1987).
Assuming shock velocities 100-200 km s-1 and a density of the pre-shock medium 103 cm-3,
we get the cooling distance of the shock to be
30 AU
4 mas, much smaller then
the resolution of the current observations.
Since the motions of the Northern Arm filaments fail to be completely explained simply by gravitational influence, one has to include additional forces that could drive the observed motions. In the next section we consider the winds at the Galactic Center as the main candidates.
As already discussed in the introduction, less than 1% of the material supplied by massive stars in the central parsec of the Galaxy is actually available for the accretion onto SgrA*. Most of the material never gets close to the BH, i.e. it must be blown away from the center and therefore can interact with the material of the mini-spiral.
There is a body of observational evidence that supports the interaction of a wind with the local ISM as well as with stars in the central parsec. Some of those were already presented in the introduction (mini-cavity, IRS 3, IRS 7, low extinction region). Then, the features presented in this work: the narrow filaments could be the contact surfaces between the wind and the mini-spiral material. The "V''-shaped sources X3 and X7 are supporting the SgrA* wind scenario (see Sect. 6.1.2).
The main candidates for the source of the central outflow are the high-mass-losing stars at the center. Alternatively, the outflow could be due to a jet-like structure directly associated with the black hole at the position of SgrA*. Some of the theoretical models that successfully fit the observational data include such a jet (e.g. Yuan et al. 2002; Falcke & Markoff 2000). In Appendix C we consider the possibility that the outflow could, in fact, be partially collimated.
Most of the above mentioned bright and massive He-stars are located
within a single plane and most of them within a radius of 0.1 pc (2.5
)
from SgrA*
(Genzel et al. 2000; Levin & Beloborodov 2003; Genzel et al. 2003),
forming a disk of clockwise rotating stars (CWS).
Paumard et al. (2006) largely expand the sample of stars
belonging to the CWS disk, and also propose
the existence of a second, less populated stellar disk containing counter-clockwise rotating
stars (CCWS; see also Genzel et al. 2003).
According to Paumard et al. (2006) most of the bright stars in the IRS 16 complex are part of
the clockwise system, including IRS 16C and IRS 16SW.
Since the bright He-stars in the CWS are also those that show clear indications for P Cygni profiles
we assume that the stars in this disk are also the predominant source of mass loss in the central parsec.
Since the geometry of the CWS plane is such that the rotational axis
is inclined towards the observer, the bright He-stars to the north and
northeast of SgrA* are located on the far
side of the disk (i.e. IRS 16C, IRS 16CC, E29, E38; Beloborodov et al. 2006).
Therefore they must be as close to the
Northern Arm - which is passing by the IRS 16 cluster from behind -
as the IRS 16 cluster is with respect to the mini-cavity.
It is therefore likely that their stellar wind may have a similarly strong effect on
the Northern Arm as it may have
towards the mini-cavity.
A particularly interesting feature in this context is feature X7. Despite
the high velocity component to the northwest, the cometary or bow-shock-like shape of the
feature is pointing approximately toward SgrA*.
The wind from the direction of SgrA* could be responsible for this shape.
If the star is moving within an area with no significant amounts of ISM then any
supersonic wind from the direction of SgrA* could easily produce such a shape.
If the amount of the ISM is not negligible then the wind
should have a much higher velocity than the observed velocity of the feature.
Taking into account the slight displacement
of the line of symmetry of X7 from the line that connects it with SgrA*, a simple
calculation tells us that the velocity of the wind from SgrA* should be
at least one order of magnitude higher than the velocity of the feature.
In addition, approximately along the same line that connects SgrA* and X7,
lies the feature X3 with a same type of a cometary appearance (see Fig. 4).
The alignment of X3 and X7 with SgrA*, and the shape of both features, give a
strong support to the central outflow hypothesis.
Recent VLT NACO -band polarization measurements (Meyer et al. 2006b) have
shown that source X7 is also polarized by
% along a position angle
of -34
10
.
Assuming that the polarization is due to Mie scattering of a non-spherical
dust distribution then an overall direction of the bow-shock's symmetry axis
is 56
10
.
This includes the direction towards SgrA*
(see the description of data acquisition and calibration in Eckart et al. 2006;
and Meyer et al. 2006b).
We have presented L'-band (3.8 m) adaptive optics observations of an
approximately
arcsec2 (
pc2) region around SgrA*
obtained with the NAOS/CONICA system at ESO VLT.
The images reveal a large number of thin filaments associated with the mini-spiral,
that could be due to the interaction of a fast wind with an ambient ISM.
The previous detection of some of the filaments in the radio regime, as well as in
Br
and Pa
lines, shows that
the filaments are not only associated with the continuum emission by dust, but also
with the ionized gas component of the ISM.
High-resolution spectroscopy is needed that could show if there are shock-diagnostic lines
coming primarily from the filaments in order to confirm our interpretation. However, there
are a number of observational challenges involved (e.g. extinction, stellar contamination, resolution etc.)
that may preclude a clear answer at this time.
We performed proper motion measurements of these features on the time basis of approximately 4 years. The filaments located west from SgrA* are curved with their convex sides westwards and show a significant proper motion to the west. The comparison with the Keplerian orbit fitting of the Northern Arm presented in Paumard et al. (2004), tells us that the process that drives the motion of the filaments cannot be purely due to the gravitational potential of the supermassive black hole at the center.
We propose the following scenario: An at least partially collimated outflow is emanating from a combination from the disk of high mass-losing He-stars and an outflow associated with the black hole at the position of SgrA*. This scenario represents a unified model that would explain the mini-cavity, the two cometary shaped sources X3 and X7, as well as the northern and southern lobes in the NIR H2 emission (see Appendix C). It would also be in agreement with the recent NIR polarization data that is indicative for a temporal accretion disk around SgrA* and would be supportive for a collimated outflow that is often invoked to explain the properties of SgrA*.
Acknowledgements
We thank P. Pejovic for providing his K-band stellar proper motion measurements. Part of this work was supported by the Deutsche Forschungsgemeinschaft (DFG) via SFB 494. K. Muzic and L. Meyer were supported for this research through a stipend from the International Max Planck Research School (IMPRS) for Radio and Infrared Astronomy at the Universities of Bonn and Cologne.
Table A.1:
Comparison of stellar proper motions in - and L'-band. All velocities are in km s-1.
In Sect. 3.2 we explain how L'-band image frames from different epochs have been transformed to
the coordinate system of a reference frame using IDL image transformation routines.
In order to validate our method and to demonstrate that (i) reliable proper motions can be
obtained from the lower resolution L'-band data and (ii)
that the resulting transformed images therefore are indeed astrometric and suitable
for the measurements conducted in this work,
we compare -band and L'-band stellar proper motions.
The PSF fitting was done using the StarFinder code (Diolaiti et al. 2000).
While the L'-band positions were obtained from transformed images
we derived the
-band positions from the original image frames and applied
the corresponding transformation to the resulting lists.
So the basic difference between two sets of proper motions is in the transformation of
data to the reference coordinate system.
As can be seen from Fig. A.1, 97
of our
- and L'-band stellar proper motion
results are in agreement to within 3
.
Also, there is no systematic offset between
proper motion difference when projected in RA and Dec (see histograms in Fig. A.2).
This demonstrates that the L'-band image transformation we used produces astrometric frames.
As an example, we list in Table A.1. proper motions for some of the sources in the central parsec obtained in our study
(
- and L'). We also add the values published in Genzel et al. (2000). The uncertainties of the proper motions
represent the 1
deviation of the linear fit.
Burton & Allen (1992); Gatley et al. (1986,1984) and Yusef-Zadeh et al. (2001)
report H2(1-0) S(1) molecular hydrogen line emission from the CND that is
concentrated on the northern and southern lobes.
Based on HST NICMOS and ground based data and adopting a lower limit
to the extinction toward the lobes of AK=3 (Scoville et al. 2003),
Yusef-Zadeh et al. (2001) derive the extinction-corrected peak line
intensity from the lobes as approximately
to be about a factor
higher than the mean surface brightness of the CND.
The source of excitation of the H2 emission associated with the
CND is unclear.
Most commonly it is assumed that the
70
inclination
of the 0.5 pc thick CND torus toward the line of sight will
produce the north and south lobes by
limb brightening of the
0.5 pc torus
(e.g. Yusef-Zadeh et al. 2001; Gatley et al. 1986; Burton & Allen 1992).
As summarized by Yusef-Zadeh et al. (2001), neither the
intense UV filed combined with the high densities, nor the shock processes
due to isotropic wind originating from the hot stars at the central parsec, can
explain such a high intensity of molecular hydrogen emission on the CND lobes.
The latter consideration shows that the ram
pressure of an isotropic wind from IRS 16 cluster would be at least
a factor
too low compared to what is required.
In this section we consider the possibility that the outflow originating from the central part of the Galaxy, presumably from the cluster of high mass losing stars, with a possible contribution from the massive black hole itself, might actually be partially collimated. In the case of the existence of a wind from SgrA*, the collimation would have to be seen as a consequence of the yet unexplained accretion process onto the black hole but could be linked to the plane in which the mass-losing He-stars are moving in. This motion within a single plane implies that the viscous friction of gas being blown into that plane is larger than that for gas being blown out perpendicular to it. On the one hand that results in a pressure gradient perpendicular to the plane, on the other hand it will result in an accretion flow within that plane onto the central black hole and may give rise to a corresponding outflow perpendicular to it. So both the ordered motion of the mass-losing stars as well as the resulting structured accretion flow may result in a combined, possibly collimated outflow.
This scenario would also explain the bow shock "tail'' of the IRS 7
and the northern and southern lobes
within the CND to be associated with the outflow.
If we assume that the H2 line emission from the lobes is dominated
by the interaction with a collimated wind rather than an approximately
isotropic stellar wind from the central few arcseconds
then - due to the limited size of the two lobes in the CND -
a geometrical concentration of the outflow over an isotropic flow by
a factor
has to be taken into account.
We note that such an outflow could also have a grazing incident onto the CND
(see Fig. C.2).
In this case the opening angle could be larger and the compression factor smaller.
We also infer the lower limit on the collimated outflow velocity from the fact that
NE1, the northern most filament, has an insignificant proper motion
perpendicular to its extent and that filaments
NE3 and NE4 which are likely to be more affected by the IRS 16 cluster wind have
proper motion velocities of
km s-1.
When generated in the Northern Arm at angular distances of 3 to 7 arcseconds
from the IRS 16 cluster or SgrA* and taking
as an upper limit of the
proper motion inferred by a collimated outflow, then a lower limit to the velocity of that flow
with a total opening angle of
-30
may be of the order of
/tan(
)
500-900 km s-1.
The total required mass-loss carried by the collimated outflow is then
![]() |
(C.1) |
Whether in the presence of the deep central potential
provided by the
6
black hole
the wind contribution is dominated by the stars or by SgrA*
has to be investigated via detailed hydro-dynamical model calculations.
In the context of the mini-cavity such calculations have been
carried out by Melia et al. (1996) and most recently by Cuadra et al. (2006).
The authors have shown that the structure, kinematics, and luminosity
of the mini-cavity are consistent with
a model in which a few 106
black hole
interacts with an
ambient 500-700 km s-1 Galactic Center wind originating within
the IRS 16 cluster.
It accretes a small portion of it and then may expel the rest
in form of a collimated flow toward the mini-spiral streamer
(Melia et al. 1996).
The filaments along the Northern Arm presented
here, the shape of the Northern Arm, as well as the approximate alignments
between the rotation axis of the He-star disk and the connecting
line between the northern and southern lobes in the CND
(see references given above) present evidence for a similar collimated
flow to the northeast.
![]() |
Figure C.2: Positions of main features at the Galactic Center as seen by an observer at the center of the Galaxy. Two crosses indicate the normal to the clockwise stellar disk plane: disk A as found by Levin & Beloborodov (2003) and disk B by Paumard et al. (2006). The large light-gray-shaded area south of the galactic plane represents the results from NIR polarization measurements as presented by Eckart et al. (2006) and Meyer et al. (2006b,a). Here we only show the direction towards the Northern Arm. The dark-gray-shaded area close to the center of the plot outlines the approximate direction of a possible collimated outflow that could explain out new observations plus the H2 results by Burton & Allen (1992); Gatley et al. (1986,1984) and Yusef-Zadeh et al. (2001). Such an outflow would also be in agreement with mass-loss from the the stellar disks and the recent polarization data (see references above). |
In the following we discuss the space orientation of the main structures at the Galactic Center, as well as the orientation of the proposed outflow.
The observed geometry in the plane of the sky is such that the Northern Arm filaments and the
mini-cavity are on almost opposite sides and on a well defined
connection line through SgrA*. This line also coincides with the strip of low extinction
presented by Schödel et al. (2007). The wind that we consider responsible for
the filaments along the Northern Arm and
the formation of the mini-cavity as well,
is strongly interacting with the gas in the gas disk (i)presented by Vollmer & Duschl (2000). This disk contains the
Northern Arm and the part of the Bar, and connects well with the CND at an inclination
of 30
to the line of sight. The disk (i) is inclined
25
to
the line of sight and has
the position angle of 28
on the sky (N to E).
Both CND and the Northern Arm are at a
relatively low inclination
to the plane of the Galaxy (Paumard et al. 2004; Vollmer & Duschl 2000).
Thus we conclude that the wind could in fact be
collimated with an opening angle of the order of 30
or less.
Another indication for an outflow to the northeast
comes from the wavy structure of the Northern Arm.
Here Yusef-Zadeh & Wardle (1993) already noted that
the distortions of the arm are produced by ram
pressure of the wind incident on the surface of
the arm and inducing there Rayleigh-Taylor or
Kelvin-Helmholtz instabilities.
The position angle of the outflow projected on the sky is close to that of the galactic plane. This would also be in good agreement with the extent of the northern and southern lobes in the H2 emission of the CND. These lobes are located at the outer tips of the projected ellipse of the rotating CND. Therefore the lobes must be located approximately within the plane of the sky and within the Galactic plane.
This comparison also shows that the rotation axis of the hot He-stars
(Paumard et al. 2006; Genzel et al. 2000,2003; Levin & Beloborodov 2003)
and the axis of the outflow may in fact be aligned to within less then 30.
Both results for the disk of clockwise-moving stars from Levin & Beloborodov 2003 (disk A)
and Paumard et al. 2006 (disk B) are shown in the Fig. C.2. The rotational axis of
the disk A is closer to the direction of the postulated outflow than the axis of the
disk B. The reason for this behavior could be the choice of the stellar sample in both
measurements.
The sample in Levin & Beloborodov (2003) includes 10 of the 13 brightest
brightest He stars, most of them designated as Ofpe/WN9 stars by Paumard et al. (2006).
Disk A therefore contains the largest number of heavily mass-loosing stars.
Eckart et al. (2006) obtained polarization measurements of SgrA* at two observing epochs
(July 2005 and June 2004) and report the mean position angle of the E-vector to be about 60
with a swing in polarization angle of up to 40-50
.
Meyer et al. (2006a) find that the inclination of the possible accretion disk around SgrA* has
a lower limit of about 30
.
It is interesting to point out
that the range over which polarization angle varies on the plane of the sky, centered on SgrA*
(see Fig. 8 in Eckart et al. 2006), approximately coincides with the opening angle of the central
outflow proposed in this work. These two phenomena could be connected under the assumption that the accretion disk,
or the possible jet, may be polarized and the
E-vector is positioned perpendicularly to the accretion disk, or along the jet
(see the discussion in Eckart et al. 2006; and Meyer et al. 2006a).
Fig C.2 shows the orientations of the discussed structures
for an observer at the center of the Galaxy.