EDP Sciences
Free Access
Issue
A&A
Volume 576, April 2015
Article Number L16
Number of page(s) 4
Section Letters
DOI https://doi.org/10.1051/0004-6361/201525969
Published online 16 April 2015

© ESO, 2015

1. Introduction

With a mass of M ≈ 4.3 × 106M and located at the center of the Milky Way at a distance of R0 ≈ 8 kpc, Sagittarius A* (Sgr A*) is the nearest supermassive black hole (see, e.g., Genzel et al. 2010 for a review). Thanks to its proximity, the Galactic center serves as an excellent laboratory for the astrophysics of galactic nuclei; for the given mass and distance, 1 m ≡ 8 a.u. ≡ 94 RS, with RS denoting the Schwarzschild radius. Sgr A* is characterized by its low luminosity compared to active galactic nuclei (AGN), with a bolometric luminosity 2 × 10-8 of its Eddington luminosity (see Narayan et al. 1998 and references therein). Sgr A* shows a slightly inverted radio spectrum peaking at mm-to-submm wavelengths (Zylka et al. 1992; Falcke et al. 1998; Melia & Falcke 2001; Bower et al. 2015).

The emission mechanism of Sgr A* is a matter of ongoing debate. On the one hand, radiatively inefficient accretion flows (RIAFs; e.g., Yuan et al. 2003) successfully reproduce the observed fluxes at various wavelengths. On the other hand, multiple studies (Falcke & Markoff 2000; Markoff et al. 2001, 2007; Yuan et al. 2002; Falcke et al. 2009) have pointed out that parts of the emission need to originate from jets, which is supported by observational evidence for the presence of jets (Yusef-Zadeh et al. 2012; Li et al. 2013). The current lack of structural information on Sgr A* prevents an unambiguous decision between competing theories. At the short-wavelength side of the submm-bump, where the emission becomes optically thin, instruments with the necessary spatial resolution are not available. At wavelengths of more than a few millimeters, the source structure is washed out by interstellar scattering (e.g., Bower et al. 2006). Lu et al. (2011) found evidence that the shape and orientation of the elliptical Gaussian changes with frequency; this could be interpreted as an intrinsic structure that is slightly misaligned with the scattering disk, which shines through toward higher frequencies. In addition, recent Very Long Baseline Interferometry (VLBI) observations at 43 GHz have been able to resolve an intrinsic elliptical structure with a preferred geometrical axis (Bower et al. 2014). This structure might indicate an existence of jets, but could also be the result of an elongated accretion flow such as a black hole crescent (e.g., Kamruddin, & Dexter 2013).

Recently, a gas cloud labeled G2 was observed to move toward Sgr A* on a nearly radial orbit (Gillessen et al. 2012). So far, mainly two possible structures of G2 have been discussed, with the first scenario being that G2 is a localized overdense region within an extended gas streamer. This agrees with observations reporting that G2 is composed of a compact head and a more widespread tail (Gillessen et al. 2013b; Pfuhl et al. 2015). The two components are on approximately the same orbit and are connected by a faint bridge in position-velocity diagrams, in Gillessen et al. (2013a) and lasted over one year,while G2 has been stretched substantially along its orbit by tidal shearing caused by the gravitational potential of Sgr A* (Pfuhl et al. 2015). Test particle simulations have provided a good explanation for the dynamics of G2; the results showed that hydrodynamic effects have not been significant (e.g., Pfuhl et al. 2015; see also Schartmann et al. 2012). The second possibility under discussion is that G2 is a circumstellar cloud around a star that provides stabilizing gravity and continuously replenishes the gas. Studies supporting this scenario highlighted that the tail might only be a fore- or background feature and might not be physically connected to G2 (Phifer et al. 2013; Valencia et al. 2015). Recent observations have also shown that G2 survived its pericenter passage as a compact source, which might be a clue that G2 is a star enshrouded by gas and/or dust (Witzel et al. 2014; Valencia et al. 2015). This scenario was modeled analytically (Miralda-Escudé 2012; Scoville & Burkert 2013) and explored with the help of numerical simulations (Ballone et al. 2013; De Colle et al. 2014; Zajaček et al. 2014).

Interactions with the accretion flows toward Sgr A* might cause G2 to lose angular momentum, and as a result, parts of it may be accreted by the black hole (Anninos et al. 2012; Burkert et al. 2012; Schartmann et al. 2012). Accordingly, an increased radio luminosity (e.g., Mahadevan 1997; Mościbrodzka et al. 2012) as well as an increase in source size might be expected (e.g., Mościbrodzka et al. 2012). In addition, radio-bright outflows such as jets or wind-like outflows related to RIAFs (Yuan et al. 2003; Mościbrodzka et al. 2012; Liu & Wu 2013) might become observable on spatial scales of 1 mas. To search for such transient structures, we performed VLBI observations with four European stations of the Global Millimeter VLBI Array (GMVA) at 86 GHz, providing an angular resolution down to 0.3 mas. Our observing frequency of 86 GHz is in a region where scattering vanishes and is weaker in the images of Sgr A* at frequencies lower than about 43 GHz (Bower et al. 2004, 2006; Shen et al. 2005).

2. Observations and data analysis

We observed Sgr A* at 86 GHz on 1 October 2013 using four GMVA stations of the European VLBI Network (EVN): Effelsberg (EF), Pico Veleta (PV), Plateau de Bure (PB), and Yebes (YS). A combination of the low declination of Sgr A* (−29°), the high latitude of the stations (the minimum latitude being 37° for Pico Veleta), and a technical problem at Yebes station limited our observing times to around 1.5 h for Yebes and 2.5 h for the other stations. We observed both circular polarizations except for Yebes station, where only left circular polarization (LCP) data were obtained. The data were recorded with the Mark 5 VLBI system (two-bit sampling) using the digital baseband converter (DBBC) in polyphase-filter bank mode with a bandwidth of 32 channels, each 16 MHz wide (total bandwidth 512 MHz). All data were correlated with the DifX VLBI correlator of the MPIfR (Bonn, Germany). Our calibrators (fringe-finding sources) were NRAO 530 and 1633+38.

We followed the standard procedures for initial phase and amplitude calibration using the AIPS software package (Greisen 2003) and applied phase self-calibration using the Caltech Difmap package (Shepherd 1997) with a point-source model. The geometry of our observation resulted in a very elongated beam (point spread function) with full widths at half maximum (FWHM) of 3.02 mas and 0.33 mas, respectively, at a position angle of −22°. For the flux density of Sgr A* we found a value of 1.4 ± 0.3 Jy).

We used the evolution of closure phases to search for asymmetric, extended emission around Sgr A*. A closure phase is the sum of the interferometric phases of the three baselines in a closed triangle of stations; it is free from antenna-based phase errors (e.g., Thompson et al. 2001). The closure phase for a centrally symmetric brightness distribution is zero. The amount of deviation from zero and the timescales of fluctuations of the closure phases can be used to probe the structure of an asymmetric source even if it cannot be imaged (due to lack of flux or insufficient uv plane coverage). This technique was previously used by Krichbaum et al. (2006) and Lu et al. (2011), who found the closure phases of Sgr A* to be consistent with zero throughout their observations in October 2005 and May 2007.

thumbnail Fig. 1

Closure phases vs. time for three independent VLBI triangles. Each data point denotes a closure phase measurement averaged over one scan of six minutes; error bars correspond to 1σ errors. Dotted lines indicate weighted averages of all data points, dashed lines represent the corresponding standard errors of mean. Average phase values are noted in each diagram. Phases for the triangle EF-PB-PV are extracted from Stokes I data, phases for the other triangles from LL data.

Open with DEXTER

We extracted the closure phases from the visibility data for three independent triangles of the VLBI stations, initially binned into 10-second time bins. We flagged data with large errors (larger than the standard deviation of the data for a given triangle) and obvious outliers. The fraction of flagged data is 8.9%, 7.9%, and 6.0% for the triangles EF-PB-PV, EF-PV-YS, and PB-PV-YS, respectively; the difference between the results obtained with and without flagging is insignificant, however. We took the weighted average of the remaining values for each scan of 6 min; the resulting data set is shown in Fig. 1. We obtained a combined reduced χ2 value of 1.34 for all the closure phases for all the independent triangles based on the null hypothesis that the closure phase is zero all the time. The corresponding false-alarm probability (p-value) is 0.136; this value is clearly too high to reject the null hypothesis (this outcome did not change when choosing bin sizes other than 6 min). Consequently, we conclude that the closure phases agree with zero during the observing time.

3. Discussion

thumbnail Fig. 2

Outflow orientations assumed in our simulation. The gray ellipse (with FWHM and position angle as noted) indicates the beam, dashed lines indicate four outflow orientations: along the major axis, along the minor axis, along the jet direction claimed by Li et al. (2013), and along the direction claimed by Yusef-Zadeh et al. (2012).

Open with DEXTER

The absence of a deviation of the closure phases from zero – within typical errors of a few degrees – implies the absence of asymmetric structure around Sgr A* at the time of our observation. A priori, bipolar, symmetric jets can account for zero closure phases. However, if the jet axis is oriented partly along the line of sight, the flux observed from the jet approaching the observer (f+) is amplified while that from the receding jet (f) is dimmed by the Doppler effect. If the jets on both sides of the black hole have the same luminosity intrinsically, the observed ratio of the two fluxes is (1)(e.g., Beckmann & Shrader 2012), where β is the jet speed in units of speed of light, θ is the angle between the jet axis and the line of sight, and α is the spectral index defined via the flux density Sννα. Various observational constraints (α ≈ 0.5, β ≲ 0.1; Lu et al. 2011) and model predictions for transient jets (θ ≈ 60°, β ≈ 0.7; Falcke et al. 1998, 2009) place the flux ratios to be expected at between 1 and 6. Accordingly, even intrinsically symmetric jets from Sgr A* might appear asymmetric in the observer frame.

To constrain the size of outflows that we could have missed, we performed simple simulations. We placed artificial, unipolar secondary sources next to a primary point source model representing Sgr A* and compared the closure phases obtained from the resulting artificial visibility data with the observations. We considered two geometries: a single-point source and a jet composed of ten equally spaced point sources (knots) with equal fluxes. We probed four orientations for the simulated outflows (see Fig. 2): along the major axis, along the minor axis of the beam, the jet direction claimed by Li et al. (2013), and the jet direction claimed by Yusef-Zadeh et al. (2012). We used total fluxes of 0.2 Jy and 0.55 Jy for the artificial sources; these values ensure that our simulated outflows are sufficiently faint to not violate the constraints given by the known recent brightness evolution of Sgr A* (0.2 Jy from the mean variability of 15% from June 2013 to February 2014 at 41 GHz, with 0.55 Jy corresponding to the strongest variation in the same period, Chandler & Sjouwerman 2014). For each simulation setup, we measured the average of absolute values of the closure phases for each triangle. We varied the distances of the model sources (for the jet model: the largest distance) from Sgr A* until we found a critical distance at which the absolute values of the simulated closure phases exceeded those of the observations by more than the 1σ error at all triangles. We summarize our results in Table 1. As expected, the critical distances are smaller for brighter outflows. Jet-like structures lead to larger critical distances than equally luminous single, compact sources. As a consequence of the very elongated beam, the critical distances for sources located along the major axis of the beam are larger by a factor of 7 than for those located along the minor axis. In a few cases (denoted “N/A” in Table 1), the absolute values of the simulated closure phases were similar to those of the observations for all distances of the model sources, meaning that we were unable to identify a critical distance. Overall, our observations limit the extension of asymmetric (in the observer frame) jet-like outflows from Sgr A* to projected distances of 2.5 mas along the major axis and 0.4 mas along the minor axis.

Table 1

Size constraints for outflows from Sgr A*.

Our analysis limits the (projected) extension of linear outflows to about 232 and 35 Schwarzschild radii, respectively, for outflows with fluxes of about 0.2 Jy; obviously, outflows substantially fainter could be larger and still remain undetected. Unfortunately, the resolution of our observations is not sufficient to probe structures in accretion flows that are expected to occur on scales 10 RS (cf., e.g., Broderick et al. 2011); those observations will probably have to await the Event Horizon Telescope (cf. Fish et al. 2014). When referring to radio galaxies for comparison, especially to M 87, which has a central black hole with small angular diameter (10μas) second only to Sgr A*, a distance of tens of Schwarzschild radii appears to be critical for the formation of AGN jets. Recent VLBI observations of the jet of M 87 find a transition in the collimation geometry at a distance of about 100 RS, with the jet opening angle being smaller outside this boundary; this has been interpreted as 100 RS being the characteristic distance for the conversion of magnetic to kinetic energy in a magnetically launched jet (Hada et al. 2013). Accordingly, Sgr A* is potentially a very important test case for AGN jet physics if a jet is ever detected.

Our nondetection of outflows is consistent with earlier null results from VLBI observations of Sgr A* at 86 GHz (Krichbaum et al. 1998; Lo et al. 1998; Shen et al. 2005; Lu et al. 2011) even though we observed Sgr A* at an epoch of potentially increased accretion. Our results are in line with other recent observations finding that Sgr A* has been quiescent from radio to X-rays in 2013 and 2014 (Akiyama et al. 2013; Brunthaler & Falcke 2013; Chandler & Sjouwerman 2014; Degenaar et al. 2014). In addition, the zero closure phases within errors can set some constraints on the substructure in scattering disks of Sgr A* (e.g., Gwinn et al. 2014); this substructure must remain symmetric on spatial scales of from submas to a few mas depending on an axis in the geometry of the substructure.

As noted above, the radio flux density that can realistically be attributed to accretion of parts of G2 is Sν ≈ 0.2 Jy, translating into a radio luminosity (for ν = 86 GHz). Without making further assumptions, we can state that this value corresponds to an effective accretion rate η = LR/c2 ≈ 1.5 × 109 kg s-1 ≈ 7.7 × 10-9M yr-1, with η ∈ [0,1] being the matter-to-light conversion efficiency and c denoting the speed of light. However, in addition, we have to take into account that accreted matter might not be converted into electromagnetic radiation, but into jets, with the presence of jets being an ad hoc working hypothesis. In this case, a rough quantitative estimate of the accretion rate is provided by the relation of jet power to accretion power of radio galaxies: for highly sub-Eddington accretion (as is the case for Sgr A*), kinetic jet power Pjet and accretion power Paccc2 are related as Pjet ≈ 0.01Pacc (e.g., Allen et al. 2006; Trippe 2014). The powers Pjet and LR in turn are empirically relatedas Pjet ≈ 5.8 × 1036 (LR/ 1033)0.7 W (Cavagnolo et al. 2010). Combining the two relations and using, again, LR ≲ 0.2 Jy implies an accretion rate ≲ 1017 kg s-1 ≈ 0.5 M yr-1. We note that this calculation assumes a highly idealized situation, neglecting interactions between G2 and the accretion flow around Sgr A*. Given the low accretion rates as well as the complexity of the accretion flows, it seems realistic that relatively large amounts of matter could be “peeled off” G2 and driven out of the accretion zone by winds or other noncollimated outflows.

These limits on accretion rates, which are always substantially lower than the total mass of G2 (with details depending on which structure of G2 is assumed), are consistent with the observed kinematics of G2 during its pericenter passage: as noted by several studies, the orbit was purely Keplerian even after the pericenter passage (Witzel et al. 2014; Pfuhl et al. 2015; Valencia et al. 2015). This indicates that G2 did not experience a notable loss of angular momentum and energy, which indicates rather weak interactions with the hot gas around Sgr A*. As suggested by several studies based on hydrodynamical simulations, the viscous timescale can be on the order of years (Burkert et al. 2012; Schartmann et al. 2012; Mościbrodzka et al. 2012), meaning that it could take a few more years to see AGN-like (mostly in radio mode associated with hot accretion flows and jets, see Yuan & Narayan 2014 for a review) activity in Sgr A*. Overall, our analysis suggests that G2 is mostly stable against loss of angular momentum and subsequent (partial) accretion at least on timescales 1 yr.

Acknowledgments

We are grateful to Eric W. Greisen (NRAO, U.S.A.) for technical support. We acknowledge financial support from the Korean National Research Foundation (NRF) via Global Ph.D. Fellowship Grant 2014H1A2A1018695 (J.-H. P.) and Basic Research Grant 2012R1A1A2041387 (S.T.). We made use of data obtained with the Global Millimeter VLBI Array (GMVA), which consists of telescopes operated by the MPIfR, IRAM, Onsala, Metsahovi, Yebes, and the VLBA. The data were correlated at the correlator of the MPIfR in Bonn, Germany. Last but not least, we are gratful to the anonymous referee for a helpful review and valuable suggestions.

References

All Tables

Table 1

Size constraints for outflows from Sgr A*.

All Figures

thumbnail Fig. 1

Closure phases vs. time for three independent VLBI triangles. Each data point denotes a closure phase measurement averaged over one scan of six minutes; error bars correspond to 1σ errors. Dotted lines indicate weighted averages of all data points, dashed lines represent the corresponding standard errors of mean. Average phase values are noted in each diagram. Phases for the triangle EF-PB-PV are extracted from Stokes I data, phases for the other triangles from LL data.

Open with DEXTER
In the text
thumbnail Fig. 2

Outflow orientations assumed in our simulation. The gray ellipse (with FWHM and position angle as noted) indicates the beam, dashed lines indicate four outflow orientations: along the major axis, along the minor axis, along the jet direction claimed by Li et al. (2013), and along the direction claimed by Yusef-Zadeh et al. (2012).

Open with DEXTER
In the text

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.