Free Access
Issue
A&A
Volume 587, March 2016
Article Number A37
Number of page(s) 9
Section Astrophysical processes
DOI https://doi.org/10.1051/0004-6361/201527286
Published online 12 February 2016

© ESO, 2016

1. Introduction

The compact radio source Sagittarius A* (Sgr A*) is commonly assumed to be associated with the supermassive black hole (SMBH) of ~4.0 × 106M at a solar distance of ~8.0 kpc in the center of the Milky Way (Eckart et al. 2002; Schödel et al. 2002, 2003; Ghez et al. 2003, 2008; Gillessen et al. 2009). Despite the fact that Sgr A* is extremely dim in terms of Eddington luminosity, based on the correlation between black hole (BH) mass and the velocity dispersion, most if not all Galactic nuclei contain an SMBH at their centers (e.g., Richstone et al. 1998; Gebhardt et al. 2000; Ferrarese & Merritt 2000). Because it is about one hundred times closer than the second nearest Galactic nucleus (M31) and because it has the largest projected Schwarzschild radius on the sky, Sgr A* is the most interesting target in which to study the physics of these objects.

thumbnail Fig. 1

Uniformly weighted clean LCP maps of Sgr A* observed on a) May 16 2012 (7:55–12:19 UT) with a flux of 1.11 ± 0.02. b) May 17 2012 (6:04–12:19 UT) with a flux of 1.34 ± 0.03. c) May 18 2012 (7:32–12:19 UT) with a flux of 1.17 ± 0.02. The maps are restored with a beam of 2.74 × 1.12 at 1.76°. The contour levels are 1.73%, 3.46%, 6.93%, 13.9%, 27.7%, and 55.4%.

The intensity of Sgr A* shows spontaneous flux density outbursts at radio to X-ray frequencies, commonly referred to as flares (Mauerhan et al. 2005; Marrone et al. 2006; Yusef-Zadeh et al. 2008; Eckart et al. 2008a,b,c, 2009, 2012; Lu et al. 2011; Miyazaki et al. 2012). These intensity irregularities appear on timescales ranging from 1–2 h (main flares) down to 7–10 min (sub-flares) (Eckart et al. 2006a) with stronger activity at shorter wavelengths (Baganoff et al. 2001; Genzel et al. 2003; Ghez 2004; Eckart et al. 2006b,c). The short timescales suggest that this feature originates in a very compact region, possibly located close to the event horizon of the BH. Even though the detection of structures on these scales has not been achieved yet, Doeleman et al. (2008, 2009) and Fish et al. (2011) have detected an intrinsic source structure of ~40 μas on event-horizon scales using Very Long Baseline Interferometry (VLBI) at 0.7 mm and 1.3 mm.

Additionally, its apparent size is frequency dependent because of the scatter broadening of the interstellar medium (Bower et al. 2006; Krichbaum et al. 2006). It is expected that the λ2-dependency breaks at around 43 GHz and the intrinsic source size is no longer dominated by scatter broadening and thus reveales its true structure (Krichbaum et al. 1998, 2006; Lo et al. 1998; Doeleman et al. 2001; Bower et al. 2004, 2006; Shen et al. 2005), while Bower et al. (2014) reported that strong interstellar scattering is still present at 7 mm. Therefore, observations of Sgr A* with 7 mm VLBI are of particular interest to study its source size and morphology.

The current literature provides several models that try to explain the observed variability by a relativistically aberrated accretion disk emission, hot spots, inhomogeneities in the accretion disk, or a jet (Falcke et al. 2000; Falcke & Markoff 2000; Broderick & Loeb 2006; Huang et al. 2007; Broderick et al. 2009). The hot spot model proposed by Yuan et al. (2009) describes in analogy to the coronal mass ejection of the Sun that magnetic flux ropes can be formed in the accretion disk by Parker instabilities. If this system looses its equilibrium stability, the material is rapidly expelled and the flux rope is propelled away from the accretion disk. Such hot spots orbiting at detectable radii would alter the morphology of Sgr A* at timescales corresponding to their orbit periodicity. In this case, closure phases would periodically deviate from zero. Another currently discussed model is a temporary jet anchored at Sgr A* (e.g., Markoff et al. 2007). A jet can be induced by a higher accretion rate than during quiescent states or by other instabilities in the accretion flow, and temporary jets would also alter the morphology of Sgr A*, resulting in non-zero closure phases.

Relativistic magnetohydrodynamic simulations predict a constant size and shape of the Sgr A* emission region (e.g., Chan et al. 2009). However, if adiabatic expansion were the cause of these flares, a change of the morphology and/or a full width at half-maximum (FWHM) would be observed. An orbiting or asymmetrically located expanding feature would, on the other hand, be detectable by a position deviation of the emission center caused by the change of its structure. The associated closure phases would also deviate from from zero as a result of the increased asymmetry of the source. Therefore, the best properties to investigate the nature of the Sgr A* flare activity are the position, the morphology, and the FWHM Gaussian size of the compact radio source Sgr A*.

While the flux variability of Sgr A* has been studied by many authors (e.g., Zhao et al. 2003; Mauerhan et al. 2005; Marrone et al. 2006, 2008; Doeleman et al. 2008; Eckart et al. 2008a,b,c, 2009, 2012; Yusef-Zadeh et al. 2008, 2009; Li et al. 2009; Kunneriath et al. 2010; Sabha et al. 2010; Zamaninasab et al. 2010; Fish et al. 2011; Miyazaki et al. 2012), size measurements during flares have not been performed as frequently (e.g., Rogers et al. 1994; Krichbaum et al. 1998; Lo et al. 1998; Doeleman et al. 2001, 2008; Shen et al. 2005, 2006; Huang et al. 2007; Lu et al. 2011; Bower et al. 2014). Recently, Bower et al. (2014) performed triggered multi-frequency VLBI observations of Sgr A* at NIR, X-ray, and 43 GHz that revealed an elliptical intrinsic source size of 35.4 × 12.6Rs with an rms variation of ~5% and maximum peak-to-peak change of \hbox{$\lessapprox$}15%. Akiyama et al. (2013) found even stronger fluctuations of its intrinsic size of 19% at 43 GHz based on VERA observations of the emitting region. Lu et al. (2011) also reported a tendency for the minor axis to increase during higher flux periods. The positional change of Sgr A* has been investigated by Reid et al. (2008), who found an average centroid change of Sgr A* at 7 mm of 71 ± 45μas for timescales between 50 and 100 min and 113 ± 50μas for 100–200 min, which rules out hot spots that contribute more than 30% of the total flux with orbital radii above 15 GMSgrAc280\hbox{$\frac{G M_\mathrm{Sgr\,A^{*}}}{c^{2}}\approx80$}μas.

This work presents a multi-epoch measurement campaign using 7 mm Very Long Baseline Array (VLBA) data triggered by NIR observations to investigate the accretion physics and multi-wavelength variability properties of Sgr A*.

2. Observations and data analysis

2.1. Radio and NIR data

The presented 7 mm VLBA data of Sgr A* were observed on three consecutive days on May 16–18, 2012 for ~6 h each day in dual polarization mode (see Fig. 1, stations: Fort Davis, Hancock, Kitt Peak, Los Alamos, Mauna Kea, Owens Valley, Pie Town, Brewster and North Liberty). These observations were triggered by a preceding NIR flare observed at the Very Large Telescope (VLT). Dual circular polarization was recorded at each station at an aggregate bit rate of 2 Gbps (8 intermediate frequency (IF) channels at 32 MSamples/s of 2 bit) using the VLBA/Mk5c correlator.

The science target (Sgr A*) and the compact extragalactic sources (J1745283, J1748291) were observed in a duty cycle of 60 s rapidly switching to Sgr A* between each pointing (J1745283Sgr A*J1748291Sgr A*J1745283). The individual scan length on each source was ~6 s.

The calibrators 3C 345 and NRAO 530 were observed every hour for 2 min each; they served as fringe tracers and amplitude calibrators. The third calibrator 3C 279 was targeted only at the beginning of each day for 2 min between one (first and third day) and three times (second day). The data were correlated at the VLBA correlator in Socorro, NM, USA with an integration time of 1 s at 512 spectral bands per baseband channel.

The trigger was performed on a single-day 2.2 μm NIR dataset. These observations were performed on May 17 4:55–7:51 UT at the VLT with the NACO instrument. Further information on NIR observations can be gained from Shahzamanian et al. (2015), who provide a detailed discussion of this dataset.

2.2. Calibration and analysis

Calibration and editing of the data was performed using the Astronomical Image Processing System (AIPS). This software package offers standard algorithms for phase and delay calibration and fringe fitting. We also used the vlbautil package, which includes additional calibration scripts.

The amplitude calibration was performed based on measured antenna system temperatures and gain elevation curves for each station. The atmospheric opacity corrections were determined by plotting the system temperatures against the air mass and fitting the variations (AIPS task: APCAL). Furthermore, corrections for the total electron content (AIPS script: vlbatecr), Earth orientation parameters (AIPS script: vlbaeop), and parallactic angle (AIPS script: vlbapang) were applied. Fringe fitting was made using each target itself as calibrator source and performing a two-point interpolation of delays and delay rates with a solution interval of 30 s.

Images of Sgr A* and the calibrators were produced using the standard hybrid mapping and CLEAN methods of AIPS and DIFMAP (Shepherd et al. 1994) by averaging over 15 s at a map resolution of 0.3 mas. The beam was enlarged during the iterative amplitude self-calibration process in three steps from one-third over one-half up to its total initial size. Sgr A* has only been detected on short baselines between 110–340 Mega-λ (May 16: 220 Mega-λ (FD, KP, LA, OV, PT); May 17: 110 Mega-λ (FD, KP, LA, PT); May 18: 340 Mega-λ (BR, KP, NL, OV, PT)). Marginal detections at longer baselines (MK) were also acquired with fluxes at expected values, but these visibilities could not be modeled accurately and their fluxes increased to implausible values during the self-calibration procedure and were therefore disbanded. The corresponding left-handed circular polarization maps covering a complete observing day of 6 h are shown in Fig. 1.

2.3. Constraints on the data

Observing with the VLBA at 43 GHz at very low elevations is limited by several effects. The array was designed to operate at longer centimeter wavelengths, steeper gain curves, higher residual pointing, and focus errors. These and other effects are the cause for a deterioration in the data quality (see, e.g., Lu et al. 2011). Additionally, variable weather conditions and higher atmospheric opacities have a stronger effect at these frequencies than at other radio bands. Even with most cautious observing efforts, Sgr A* is susceptible to residual calibration inaccuracies because of its low elevation above the horizon for the VLBA sites located on the northern hemisphere. The low declination of the Galactic center also results in an elliptical UV-coverage and beam size (see Fig. 2), which causes lower positional accuracies along the major beam axis.

thumbnail Fig. 2

UV-coverage of Sgr A* on May 17 2012 (6:04–12:19 UT) for the stations FD, KP, LA, and PT.

The presented dataset BE061 suffered from many problems of deteriorating data quality. A variety of technical difficulties during the observations (NL: not working AC, BR: Wideband upgrade, OV: Saturated system temperatures) left only up to six stations to be used on each date for data acquisition. While the only stations present on all three days are KP, PT, and MK, baselines to the latter unfortunately did not provide Sgr A* detections. Because of changing system temperature values and other uncorrectable station- and scan-based flux errors, the calibrators 3C 345 and 3C 279 had to be excluded from most of the analysis. The remaining calibrator NRAO 530 provided acceptable data quality and shows a constant flux of (2.42 ± 0.04) Jy over the whole experiment.

thumbnail Fig. 3

NIR and 7 mm light curves on May 17 6:04–12:19 h UT 2012. a)Ks-band (2.2 μm) light curve of Sgr A* observed in polarimetry mode on 17 May 2012. The light curve shown is produced by combining pairs of orthogonal polarization channels: 0° and 90° (for further details see Shahzamanian et al. 2015). b) 7 mm DFT of NRAO 530 and Sgr A* on May 17. Sgr A* was averaged over 15 min. NRAO 530 was averaged over 2 min. c) Detrended 7 mm DFT of Sgr A* on May 17. d) 7 mm light curve derived from the two-hour maps observed on May 17 (see Fig. 4). Fluxes were obtained from delta-component maps (green) and model fitting (blue). Errors represent the formal errors of 1.7% derived from NRAO 530. Not shown are the systematic errors of ~18%, which should be corrected for by calibration.

On May 18, the flux values of Sgr A* dropped by 50% on all baselines after 10:00 UT. These visibilities are therefore considered to be faulty and were excluded from this discussion.

3. Results

3.1. Flux density

On May 17, 2012 the VLT detected a NIR flare peaking at 5:30 UT (Fig. 3a). We have plotted the discrete Fourier transform (DFT) of the complex visibilities for Sgr A* and the calibrator NRAO 530 provided by the AIPS task DFTPL (Fig. 3b). On this date an increase of its flux density from 0.97 Jy to 1.17 Jy at ~9:30 UT has been observed (see Fig. 3c). This intensity was derived from uncalibrated data and is free of any errors that might arise during calibration. This therefore is the most unbiased way to determine relative flux trends that are introduced by observational effects. The light curve of Sgr A* was corrected for the intensity fluctuations of NRAO 530 by determining the difference of every data point from the mean value in percent and subtracting these as percentages of its mean value from Sgr A*. This produces the detrended light curve of Sgr A* presented in Fig. 3c. Based on the flux measurements of NRAO 530, we adopted an error of 1.7% for all amplitudes within this dataset. It would be better to use the flux densities of closer calibrators such as the phase reference sources, but these are too weak to be used for an accurate flux estimation and are therefore not regarded in this work.

For further analysis the presented 7 mm data on this date were split into five overlapping 2 h parts and then mapped individually (see Fig. 4). Because of technical problems during the observations, only LCP data provided acceptable data quality and therefore RCP had to be disbanded. We observe increasing flux densities at 7 mm peaking between 9:00–11:00 UT (see Fig. 3d). The flux density on May 16 and May 18 remained constant at values of (1.14 ± 0.02) Jy and (1.0 ± 0.2) Jy, respectively. May 18 was only considered until 10:00 h UT because of observational problems. While Sgr A* shows constant flux densities during its quiescent states on May 16 and 18, its intensity increases from 1.23 Jy (6:00–8:00 UT) to 1.41 Jy (9:00–11:00 UT) on May 17, which agrees well with the fluxes observed at the shortest baselines before self-calibration and is well above the quiescent flux values reported on May 16 and 18. This intensity variability is seen on independent baselines, which excludes the possibility that station-based errors affect the measured flux densities. Each map was produced using the same beam and map size to achieve consistent and comparable images. The calibrated fluxes of NRAO 530 show systematic station-based correction factors of ~18% compared to its uncalibrated values, which agrees well with the values found for Sgr A*. We therefore considered them to be corrected by the presented calibration. For this analysis, we excluded some poor visibilities from PT at the start and end of the experiment. The FD station inherits an uncorrectable offset with respect to the other stations and was therefore completely excluded.

Additionally, Fig. 4 shows Sgr A* as a point-like source during the initial quiescent state and then developing an extended feature of ~(1.7 ± 0.3) mas toward the southeast during its flaring state. The apex of this change of morphology appears shortly (8:00–10:00 h UT) before the peak of the flux density at 9:00–11:00 h UT. To test the validity of this feature and rule out clean-window-biased and faulty detected components, we performed different cleaning methods by changing the beam sizes and areas covered by clean windows on all maps during the hybrid mapping process, which all developed a somewhat pronounced extended feature at the same position. Model fitting by placing circular Gaussian components to the UV-data at the regions of expected flux, as suggested by the clean maps, and solving for size, radius, and positional degree yields a reliable test of the observed source morphology. Figure 5 shows the resulting best-fit model components. The presented maps were modeled using components well above the scattering size at 7 mm and favor a two-component model during the periods of higher flux densities, as intended by the clean maps. The χ2-value of this model fit represents the goodness of the model fitting of the UV-data and is therefore a measure for the quality of the input model. Trying to model a single component for the time of the most asymmetric maps (8:00–10:00 h UT) produces higher reduced χ2 values than a two-component model (two components: χ2 = 0.34, single component: χ2 = 0.54).

Analyzing the corresponding right-handed circular polarization map reveals a similar feature. The positions of the detected delta components in the left- and right-handed circular polarization maps correspond to a position of the secondary component at a radial separation of (1.7 ± 0.3) mas (LCP) and (1.8 ± 0.2) mas (RCP) at position angles of (126 ± 16)° (LCP) and (138 ± 9)° (RCP). These values are consistent within their error limits and are another indication that the secondary component that is present at this position is a real feature and not a data artifact.

thumbnail Fig. 4

Two-hour LCP maps of Sgr A* observed on May 17 2012. a) May 17 6–8 h (6:04–08:00 UT) with a flux of 1.23 ± 0.03. b) May 17 (7:10–09:00 UT) with a flux of 1.26 ± 0.03. c) May 17 8–10 h (8:14–09:58 UT) with a flux of 1.31 ± 0.03. d) May 17 9–11 h (9:11–10:56 UT) with a flux of 1.41 ± 0.03. e) May 17 10–12 h (10:09–12:19 UT) with a flux of 1.40 ± 0.03. All maps are restored with a beam of 2.74 × 1.12 at 1.76°. The plotted contour levels are 1.73%, 3.46%, 6.93%, 13.9%, 27.7%, and 55.4%.

thumbnail Fig. 5

2 h model fits of Sgr A* on May 17 2012. Extended morphology during flare is best fitted by a two component model. a) May 17 6–8 h (6:04–08:00 UT), flux: 1.21 ± 0.03. b) May 17 (7:10–09:00 UT), flux: 1.21 ± 0.03. c) May 17 8–10 h (8:14–09:58 UT), flux: 1.26 ± 0.03. d) May 17 9–11 h (9:11–10:56 UT), flux: 1.29 ± 0.03. e) May 17 10–12 h (10:09–12:19 UT), flux: 1.28 ± 0.03. All maps are restored with a beam of 2.74 × 1.12 at 1.76°. The plotted contour levels are 1.73%, 3.46%, 6.93%, 13.9%, 27.7% and 55.4%.

Table 1

Summary of closure phases for all closure triangles on May 16 and 17.

Table 2

Summary of closure phases for all simulated closure triangles on May 17 8:00–10:00 h UT.

3.2. Closure phases

A good quantity on which to test the symmetry of a compact radio source are the closure phases, which are the sum of three baseline phases in a triangle of antennas. These phases are independent of all station-based phase errors because they are a phase quantity of the complex visibilities. A point-symmetric source would have closure phases of zero, while any deviation from its symmetry would cause non-zero values. In Fig. 7 we plot closure phases for May 17 extracted from 30-min averaged UV-data. A mean value was determined for each set of closure triangles in dual polarizations as well as an average mean phase for all triangles of (0.5 ± 0.2)° (LCP) and (0.0 ± 0.1)° (RCP) (see Table 1). The values for all triangles and polarizations are all equal within their error limits. Because the station arrays changed on May 18, there are no high-quality closure phases that can be compared to the previous days, and they can therefore not be used in this analysis.

Of most interest is the closure triangle FD-KP-PT, which shows non-zero closure phases at 8:00–10:00 h UT on May 17. While no triangle offers values exceeding 3σ, there are several values higher than 2σ. These excessive values coincide with the observed change in morphology at this time (see Fig. 4c). This trend is not reproduced by the other triangles, which also show some closure phases above the two-sigma range, but in a more random fashion.

To further investigate this effect, we tried to simulate the observed closure phase structure using the Caltech VLBI analysis program FAKE (Pearson 1991). FAKE is limited to single polarization and a single IF, but since the presented maps of this work are LCP measurements, and there is no a priori reason for simulated RCP values to be different, and this should not change the result. We simulated several datasets using different input models (see Table 2) of the source and parameters similar to the presented measurement campaign. The summary of all simulated closure phase values can be found in Table 2.

Starting with two components separated by 1.5 mas at 140° east of north (central component: 1.55 Jy, secondary component: 0.02 Jy), as suggested by the presented modelfit, we tried to reproduce the observed closure phase values by changing the position of the secondary component to 0.7 mas, 0.3 mas, and finally a single central component. All of these simulations for the time range of maximum source asymmetry (8–10 h on May 17) produced higher closure phases than observed. Therefore, we changed the method of error application to the data. The applied errors for the simulations were first determined from Tsys, station diameters, station efficiencies, and pointing efficiency, and second by changing the amplitude of an applied additive Gaussian noise (erradd). erradd specifies the amplitude of additive Gaussian noise to a set of stations. The error of stations I and J is computed by (erradd(i)×erradd(j))12\hbox{$(\textsc{erradd(i)} \times \textsc{erradd(j)})^\frac{1}{2}$}. The simulated datasets also show higher errors in the range of 0.8 to 6.7 and are all zero within these limits, while a perfect simulation without any noise produces very low errors and values similar to those that are observed. To understand this behavior, additional datasets were simulated using different seeds, which were used to generate random numbers for error application to the phases and amplitudes in FAKE. These simulations show that depending on the chosen seed, the mean closure phases for the presented triangles can differ by ~10° for each closure triangle, and since the observed effect on closure phases is about ~5°, it is not possible to accurately simulate the phases of the presented data set.

In this context, it can be expected that a secondary component at 140° will produce higher closure phases on closure triangles with an axis of maximum resolution close to this angle. Therefore, FD-KP-LA, FD-KP-PT, and FD-LA-PT should show the highest closure phase values since their maximum resolution angles are at 97°, 101°, and 125° and thus closer to the asymmetric feature than KP-LA-PT (65°). This is well reproduced by the perfect simulation without errors, but it can change significantly depending on the seed value for randomization. Even though for all simulations the highest deviation from zero is present in one of the three favored triangles, a chosen seed was capable of producing a mean value of the closure phase for KP-LA-PT of 7.2°, which is the second highest value among triangles within this simulation. The observed closure phases of our dataset show the highest closure angle for FD-KP-PT of 1.9°, but it also has its second highest value for KP-LA-PT. Even though these values might be similar to the expected trend, as suggested by the angle of maximum resolution and the perfect simulation, with the given error limits, we can not place reliable constraints on this hypothesis.

To test the probability of detecting the observed source structure with the given array, we iterated several simulations with different Gaussian noise values. Starting at similar values to our observations and increasing the noise level in 10% steps revealed that the two-component source structure can be reproduced up to Gaussian noise values of erradd = 0.16 with a χ2 = 1.22. To test the possible effect of noise being the cause of the secondary component in Sgr A*, we performed a similar simulation by generating a single-component dataset and increasing the Gaussian noise to values of up to erradd = 0.80 and checked the resulting visibilities for indications of a possibly appearing secondary component generated by random noise. The resulting visibilities showed no evidence of a secondary component for all simulated values.

By examining the time range of highest flare activity ranging from 8:00–12:00 h UT on May 17, 2012, we obtained a radial plot of the model generated for the clean maps of Sgr A* on May 17 (see Fig. 8). It is noticeable that the most point-like map on 10–12 h (blue) does not inherit a strong secondary component in its models for UV-ranges between 40–80 Mλ, while the two models for the more asymmetric maps during 9–11 h (orange) and especially during 8–10 h (red) clearly show secondary model features. This means that the data were best fit by two components for this UV-range, and it is another indication for the existence of an off-core feature. All presented independent tests fit consistently with the picture that Sgr A* experienced a change in morphology on May 17 8:00–10:00 h UT.

thumbnail Fig. 6

RCP map of Sgr A* on May 17, 2012 (8:00–10:00 h UT). The map was convolved with a beam of 2.74 × 1.12 at 1.76°. Contour levels are 1.73%, 3.46%, 6.93%, 13.9%, 27.7%, and 55.4% of the peak flux density of 1.5 Jy/beam.

thumbnail Fig. 7

Closure phases of all closure triangles on May 17 6:00–12:00 h UT with plotted standard deviation levels. a) FD-KP-PT, b) FD-KP-LA, c) FD-LA-PT, d) KP-LA-PT. We show standard errors provided by standard deviation.

thumbnail Fig. 8

Radial plot of the models obtained during the time of highest flaring activity of Sgr A* on May 17, 2012. The colors correspond to the maps of 8:00–10:00 h UT (red), 9:00–011:00 h UT (orange), and 10:00–12:00 h UT (blue).

4. Discussion

Several sources claim quasi-periodic NIR flares appearing on timescales ranging from 1–2 h (main flares) down to 7–10 min (sub-flares) (Genzel et al. 2003; Aschenbach et al. 2004; Eckart et al. 2006c). Since the discovery of these flares, many models have been discussed. These models can be distinguished by a wide range of properties such as their periodicity or the source morphology (see Sect. 1).

Sgr A* shows variable flux densities from the radio to NIR and X-ray wavelengths. First evidence of correlated multi-frequency flare activities has been made in the X-ray and radio bands (Zhao et al. 2004). Related flares at NIR and sub-mm frequencies have been detected by Eckart et al. (2008c) and Yusef-Zadeh et al. (2009). While in the NIR and X-ray regime flares appear almost synchronously, a time delay of 1.5 ± 0.5 h between strong flares in the NIR and at 345 GHz have been observed (Eckart et al. 2008a,c). A typical delay of ~100 min is assumed for high-energy flares to reach the sub-mm/mm regime (Meyer et al. 2008; Marrone et al. 2008; Yusef-Zadeh et al. 2008). Since different wavelengths observe different areas of an opaque object, multi-frequency observations provide a radial view of the source. The traveling time that a NIR event needs to reach the 7 mm area is the observed time delay between these observations and can therefore be used to connect flares in different frequencies to each other.

A cross correlation between NIR and 7 mm data results in a time delay of (4.5 ± 0.5) h with a correlation coefficient of 0.83. The error of 0.5 h represents the range of coefficient values exceeding 68%. This agrees well with the time delay of up to 5.25 h reported by Yusef-Zadeh et al. (2009) between NIR/X-ray and 43 GHz observations. If a jet feature were present, it is commonly accepted to have an appreciable inclination toward the line of sight (Markoff et al. 2007). Therefore we adopt an inclination angle θ of 80° for the following discussion. Lower inclinations are viable as well and will not significantly change the result. The observed radial separation of 1.5 mas for a time delay of (4.5 ± 0.5) h corresponds to a traveling velocity (vt) of

vt=(0.4±0.2)mash-1=(1.3±0.7)×108ms-1=(0.4±0.3)c.$$ v_{\rm t}=(0.4\pm0.2)\mathrm{mas\,h}^{-1}=(1.3\pm0.7)\times 10^{8}\,\mathrm{m\,s}^{-1}=(0.4\pm0.3)\,c. $$Taking superluminal effects into account, the true gas velocity (vg) can be calculated by

vg=vtsinθβcosθ=(0.4±0.3)c.$$ v_{\rm g}=\frac{v_{\rm t}}{\sin{\theta}-\beta \cos{\theta}}=(0.4\pm0.3)\,c. $$The MHD model by Yuan et al. (2009) proposes an expelled gas velocity (vg) of 0.8 c above the accretion disk shortly after the beginning of the flare, and therefore the presented velocity of 0.4 c can easily be achieved. Even though there are other possible cooling processes of jets than adiabatic expansion, the time delay of (4.5 ± 0.5) h agrees well with the current literature and is evidence for a causal connection between the NIR and radio events by this process. A correlation between NIR/X-ray and mm-flares together with a change in source morphology or size provides strong evidence for the adiabatic expansion model.

We have detected a mean closure phase of (0.5 ± 0.2)° (LCP) and (0.0 ± 0.1)° (RCP) for Sgr A* at 7 mm using VLBI. While this mean value is zero within its error limits with no detections above 3σ, there are a few detections with values exceeding the 2σ range. For the closure triangle FD-KP-PT these high values coincide with the change of source structure on May 17 8:00–10:00 h UT (see Fig. 4c). Since the closure phases are a measure for the symmetry, this is another indication for the true existence of the change in the morphology of Sgr A*. This trend is not reproduced by the other triangles, but the observed and simulated values show that the highest deviations from zero-closure phases occur for triangles with an angle of maximum resolution close to the position angle of the detected secondary component.

Fish et al. (2011) detected closure phase values of (0 ± 40)° for Sgr A* using mm-VLBI at 1.3 mm with an Earth-scaled triangle of SMT-JCMT-CARMA. Broderick et al. (2011) showed that accretion flow models can reproduce closure phase values up to ±30° in some cases. We report much lower closure phases at 7 mm of , which are not excluded by theoretically expected values of (45–90)° at 1.3 mm (Broderick et al. 2011).

According to Bower et al. (2014), scattering still plays a significant role at 7 mm. In this context a cause for a secondary component might be refractive noise, which can introduce artificial compact features on long baselines (Gwinn et al. 2014; Johnson & Gwinn 2015). The refractive timescale for Sgr A* at 7 mm is reported to be about three months (Akiyama et al. 2013), which is much longer than the presented two-hour sub-images. We have only detected Sgr A* on May 17 at 7 mm on baselines of up to 110 Mega-λ. The detection of a secondary component on shorter timescales and baselines as would be required for interstellar scintillation exclude this effect as a cause for the observed source structure.

In the case of an orbiting hot spot the deviation of the closure phases from zero depends on all properties affecting its symmetry and periodicity, such as the hot spot orbital size, its inclination, and the flux ratio between disk and hot spot. We have detected a change of morphology at a distance of 1.5 mas from the emission center of Sgr A* (see Fig. 4). Reid et al. (2008) were able to rule out hot spots at radii above ~80 μas that contribute more than 30% of the total 7 mm flux. This limit also applies for variable adiabatically expanding off-center components. Even though fainter hot spots at larger radii may occur, the detected radial separation of 1.5 mas on May 17 is much higher, and therefore a flux component of 0.02 Jy at this scale cannot be related to a hot spot. A jet feature is not excluded by such a change in morphology. Since closure phases and 7 mm maps provide independent tests for the source morphology, the presented consistent picture is a good indication for the correct detection of a secondary component.

5. Summary and conclusions

This work presented observations of Sgr A* acquired in the framework of a three-day global observing campaign. A preceding NIR flare observed at the VLT triggered 7 mm VLBA observations. Analysis of the DFT and mapped fluxes of Sgr A* showed a radio flare following the NIR observations by ~(4.5 ± 0.5) h on May 17, 2012 (9:30 ± 1 h UT). Splitting the data of this date into two-hour bins provided evidence of a possible change in source morphology shortly before the peak of the flare (8:00–10:00 h UT). This source morphology was adequately modeled by two circular components well above the scattering size. The best-fitting model inherits a central component of 1.55 Jy and a secondary 0.02 Jy component separated by 1.5 mas at 140° (E–N). This two-component fit is superior to a single-component model for the specific time frame.

This change in morphology should be detectable by non-zero closure phases because of its asymmetry. The phases of FD-KP-PT show increasing values during 8:00–10:00 h UT with no values above 3σ, but some are in the 2σ range. The other triangles show less prominent trends, but also have a few phases exceeding 2σ. The mean value of all individual closure phases as well as the total mean closure phase (0.5 ± 0.2)° are close to zero for all times and all triangles within their errors. While the effect on closure phases is small and might be hidden within the error limits, we have tried to place constraints on this evidence by simulating several single- and two-component datasets using the Caltech VLBI analysis software FAKE and trying to reproduce the observed values. These simulations are strongly dependent on the seed chosen to generate random noise and can differ by ~10° depending on this parameter. The observed effect on the closure phases is in the order of ~5° and can therefore not be reliably tested by this method. But it can be gained from these simulations that all simulated closure phases are, like the observed values, zero within their error limits for the presented parameters.

Furthermore, it is expected that a secondary component at a position angle of 140° should have the strongest impact on closure triangles with an angle of maximum resolution close to this angle. This is true for all simulations, while depending on the applied random errors, it can change within any of the three triangles FD-KP-LA, FD-KP-PT, and FD-LA-PT. While KP-LA-PT is never the highest, for some cases it can offer the second highest mean closure phase. This is consistent with the presented observations, which show the strongest deviation from zero for FD-KP-PT on May 17. We can therefore not place entirely reliable constraints on this hypothesis, but we showed that it is possible to simulate datasets with a two-component model that can reproduce the observed trends.

A secondary component present at the observed radial separation of ~1.5 mas would rule out hot spots, which can only occur on smaller radii, but it could be produced by an adiabatically expanding feature. The event would need traveling speeds of ~0.4 c to reach 1.5 mas within (4.5 ± 0.5) h, which is easily achievable according to the current literature.

Bower et al. (2014) reported major-axis sizes of Sgr A* as an elliptical Gaussian of 35.4 × 12.6RS at a position angle of 95° east of north. This is much lower than the discussed source morphology due to a secondary component of 0.02 Jy at 150 RS at 140° east of north.

We see indications of the existence of a secondary feature leading to an asymmetric morphology during the flaring state of Sgr A*. We cannot discard that these observed effects may be affected by random observational or weather effects. The observed morphology suggests an adiabatically expanding jet. Further observations are needed to clarify its existence.

Acknowledgments

We would like to thank the anonymous referee for the constructive comments that further improved this paper. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. C.R. was supported for this research through a scholarship from the International Max Planck Research School (IMPRS) for Astronomy and Astrophysics at the Universities of Bonn and Cologne. E.R. was partially supported by the Spanish MINECO project AYA2012-38491-C02-01 and by the Generalitat Valenciana project PROMETEOII/2014/057, as well as by the COST MP0905 action Black Holes in a Violent Universe.

References

  1. Akiyama, K., Takahashi, R., Honma, M., Oyama, T., & Kobayashi, H. 2013, PASJ, 65, 91 [NASA ADS] [Google Scholar]
  2. Aschenbach, B., Grosso, N., Porquet, D., & Predehl, P. 2004, A&A, 417, 71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Baganoff, F. K., Bautz, M. W., Brandt, W. N., et al. 2001, Nature, 413, 45 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  4. Bower, G. C., Falcke, H., Herrnstein, R. M., et al. 2004, Science, 304, 704 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  5. Bower, G. C., Goss, W. M., Falcke, H., Backer, D. C., & Lithwick, Y. 2006, ApJ, 648, L127 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bower, G. C., Markoff, S., Brunthaler, A., et al. 2014, ApJ, 790, 10 [Google Scholar]
  7. Broderick, A. E., & Loeb, A. 2006, MNRAS, 367, 905 [NASA ADS] [CrossRef] [Google Scholar]
  8. Broderick, A. E., Loeb, A., & Narayan, R. 2009, ApJ, 701, 1357 [NASA ADS] [CrossRef] [Google Scholar]
  9. Broderick, A. E., Fish, V. L., Doeleman, S. S., & Loeb, A. 2011, ApJ, 738, 38 [NASA ADS] [CrossRef] [Google Scholar]
  10. Chan, C.-K., Liu, S., Fryer, C. L., et al. 2009, ApJ, 701, 521 [NASA ADS] [CrossRef] [Google Scholar]
  11. Doeleman, S., Shen, Z.-Q., Rogers, A. E. E., et al. 2001, AJ, 121, 2610 [NASA ADS] [CrossRef] [Google Scholar]
  12. Doeleman, S., Weintroub, J., Rogers, A. E. E., et al. 2008, Nature, 455, 78 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  13. Doeleman, S., Fish, V. L., Broderick, A. E., Loeb, A., & Rogers, A. E. E. 2009, ApJ, 695, 59 [NASA ADS] [CrossRef] [Google Scholar]
  14. Eckart, A., Genzel, R., Ott, T., & Schödel, R. 2002, MNRAS, 331, 917 [NASA ADS] [CrossRef] [Google Scholar]
  15. Eckart, A., Schödel, R., Meyer, L., et al. 2006a, The Messenger, 125, 2 [Google Scholar]
  16. Eckart, A., Baganoff, F. K., Schödel, R., et al. 2006b, A&A, 450, 535 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Eckart, A., Schödel, R., Meyer, L., et al. 2006c, A&A, 455, 1 [Google Scholar]
  18. Eckart, A., Baganoff, F. K., Zamaninasab, M. R., et al. 2008a, A&A, 479, 625 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Eckart, A., Baganoff, F. K., Morris, M. R., et al. 2008b, J. Phys. Conf. Ser., 131, 012006 [NASA ADS] [CrossRef] [Google Scholar]
  20. Eckart, A., Schödel, R., García-Marín, M., et al. 2008c, A&A, 492, 337 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Eckart, A., Baganoff, F. K., Morris, M. R., et al. 2009, A&A, 500, 935 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Eckart, A., García-Marín, M., Vogel, S. N., et al. 2012, A&A, 537, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Falcke, H., & Markoff, S. 2000, A&A, 362, 113 [NASA ADS] [Google Scholar]
  24. Falcke, H., Melia, F., & Agol, E. 2000, ApJ, 528, L13 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  25. Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9 [NASA ADS] [CrossRef] [Google Scholar]
  26. Fish, V. L., Doeleman, S. S., Beaudoin, C., et al. 2011, ApJ, 727, L36 [NASA ADS] [CrossRef] [Google Scholar]
  27. Gebhardt, K., Bender, R., Bower, G., et al. 2000, ApJ, 539, L13 [NASA ADS] [CrossRef] [Google Scholar]
  28. Genzel, R., Schödel, R., Ott, T., et al. 2003, Nature, 425, 934 [CrossRef] [PubMed] [Google Scholar]
  29. Ghez, A. M., Duchêne, G., Matthews, K., et al. 2003, ApJ, 586, L127 [NASA ADS] [CrossRef] [Google Scholar]
  30. Ghez, A. M., Wright, S. A., Matthews, K., et al. 2004, ApJ, 601, L159 [CrossRef] [Google Scholar]
  31. Ghez, A. M., Salim, S., Weinberg, N. N., et al. 2008, ApJ, 689, 1044 [NASA ADS] [CrossRef] [Google Scholar]
  32. Gillessen, S., Eisenhauer, F., Trippe, S., et al. 2009, ApJ, 692, 1075 [NASA ADS] [CrossRef] [Google Scholar]
  33. Gwinn, C. R., Kovalev, Y. Y., Johnson, M. D., & Soglasnov, V. A. 2014, ApJ, 749, L14 [NASA ADS] [CrossRef] [Google Scholar]
  34. Huang, L., Cai, M., Shen, Z.-Q., & Yuan, F. 2007, MNRAS, 379, 833 [NASA ADS] [CrossRef] [Google Scholar]
  35. Johnson, M. D., & Gwinn, C. R. 2015, ApJ, 805, 180 [NASA ADS] [CrossRef] [Google Scholar]
  36. Krichbaum, T. P., Graham, D. A., Witzel, A., et al. 1998, A&A, 335, L106 [NASA ADS] [Google Scholar]
  37. Krichbaum, T. P., Graham, D. A., Bremer, M., et al. 2006, J. Phys. Conf. Ser., 54, 328 [NASA ADS] [CrossRef] [Google Scholar]
  38. Kunneriath, D., Witzel, G., Eckart, A., et al. 2010, A&A, 517, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Li, J., Shen, Z.-Q., Miyazaki, A., et al. 2009, ApJ, 700, 417 [NASA ADS] [CrossRef] [Google Scholar]
  40. Lo, K. Y., Shen, Z.-Q., Zhao, J.-H., & Ho, P. T. P. 1998, ApJ, 508, L61 [NASA ADS] [CrossRef] [Google Scholar]
  41. Lu, R. S., Krichbaum, T. P., Eckart, A., et al. 2011, A&A, 525, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Markoff, S., Bower, G. C., & Faclke, H. 2006, MNRAS, 379, 1519 [Google Scholar]
  43. Marrone, D. P., Moran, J. M., Zhao, J.-H., & Rao, R. 2006, J. Phys. Conf. Ser., 54, 354 [Google Scholar]
  44. Marrone, D. P., Baganoff, F. K., Morris, M. R., et al. 2008, ApJ, 682, 373 [NASA ADS] [CrossRef] [Google Scholar]
  45. Mauerhan, J. C., Morris, M., Walter, F., & Baganoff, F. K. 2005, ApJ, 623, L25 [NASA ADS] [CrossRef] [Google Scholar]
  46. Meyer, L., Do, T., Ghez, A., et al. 2008, ApJ, 688, L17 [NASA ADS] [CrossRef] [Google Scholar]
  47. Miyazaki, A., Tsuboi, M., & Tsutsumi, T. 2012, IAU Symp. 290, eds. C. M., Zhang, T., Belloni, M., Mendez, & S. N., Zhang, 1 [Google Scholar]
  48. Pearson, T. J. 1991, Bull. Astron. Soc., 23, 991 [Google Scholar]
  49. Reid, M. J., Broderick, A. E., Loeb, A., & Brunthaler, A. 2008, ApJ, 682, 1041 [NASA ADS] [CrossRef] [Google Scholar]
  50. Richstone, D., Ajhar, E. A., Bender, R., et al. 1998, Nature, 395, A14 [NASA ADS] [Google Scholar]
  51. Rogers, A. E. E., Doeleman, S., Wright, M. C. H., et al. 1994, ApJ, 434, L59 [NASA ADS] [CrossRef] [Google Scholar]
  52. Sabha, N., Eckart, A., Buchholz, R. M., et al. 2010, A&A, 512, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Schödel, R., Ott, T., Genzel, R., et al. 2002, Nature, 419, 694 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  54. Schödel, R., Ott, T., Genzel, R., et al. 2003, ApJ, 596, 1015 [NASA ADS] [CrossRef] [Google Scholar]
  55. Shahzamanian, B., Eckart, A., Valencia-S., M., et al. 2015, A&A, 576, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Shen, Z.-Q. 2006, J. Phys. Conf. Ser., 54, 377 [NASA ADS] [CrossRef] [Google Scholar]
  57. Shen, Z.-Q., Lo, K. Y., Liang, M. C., Ho, P. T. P., & Zhao, J.-H. 2005, Nature, 438, 62 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  58. Shepherd, M. C., Pearson, T. J., & Taylor, G. B. 1994, Bull. Astron. Soc., 26, 987 [Google Scholar]
  59. Yuan, F., Lin, J., Wu, K., & Ho, L. C. 2009, MNRAS, 395, 2183 [NASA ADS] [CrossRef] [Google Scholar]
  60. Yusef-Zadeh, F., Wardle, M., Heinke, C., et al. 2008, ApJ, 682, 361 [NASA ADS] [CrossRef] [Google Scholar]
  61. Yusef-Zadeh, F., Bushouse, H., Wardle, M., et al. 2009, ApJ, 706, 348 [NASA ADS] [CrossRef] [Google Scholar]
  62. Zamaninasab, M., Eckart, A., Witzel, G., et al. 2010, A&A, 510, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Zhao, J.-H., Young, K. H., Herrnstein, R. M., et al. 2003, ApJ, 586, L29 [NASA ADS] [CrossRef] [Google Scholar]
  64. Zhao, J.-H., Bower, G. C., Goss, C., & Liu, S. M. 2004, ApJ, 603, L85 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Summary of closure phases for all closure triangles on May 16 and 17.

Table 2

Summary of closure phases for all simulated closure triangles on May 17 8:00–10:00 h UT.

All Figures

thumbnail Fig. 1

Uniformly weighted clean LCP maps of Sgr A* observed on a) May 16 2012 (7:55–12:19 UT) with a flux of 1.11 ± 0.02. b) May 17 2012 (6:04–12:19 UT) with a flux of 1.34 ± 0.03. c) May 18 2012 (7:32–12:19 UT) with a flux of 1.17 ± 0.02. The maps are restored with a beam of 2.74 × 1.12 at 1.76°. The contour levels are 1.73%, 3.46%, 6.93%, 13.9%, 27.7%, and 55.4%.

In the text
thumbnail Fig. 2

UV-coverage of Sgr A* on May 17 2012 (6:04–12:19 UT) for the stations FD, KP, LA, and PT.

In the text
thumbnail Fig. 3

NIR and 7 mm light curves on May 17 6:04–12:19 h UT 2012. a)Ks-band (2.2 μm) light curve of Sgr A* observed in polarimetry mode on 17 May 2012. The light curve shown is produced by combining pairs of orthogonal polarization channels: 0° and 90° (for further details see Shahzamanian et al. 2015). b) 7 mm DFT of NRAO 530 and Sgr A* on May 17. Sgr A* was averaged over 15 min. NRAO 530 was averaged over 2 min. c) Detrended 7 mm DFT of Sgr A* on May 17. d) 7 mm light curve derived from the two-hour maps observed on May 17 (see Fig. 4). Fluxes were obtained from delta-component maps (green) and model fitting (blue). Errors represent the formal errors of 1.7% derived from NRAO 530. Not shown are the systematic errors of ~18%, which should be corrected for by calibration.

In the text
thumbnail Fig. 4

Two-hour LCP maps of Sgr A* observed on May 17 2012. a) May 17 6–8 h (6:04–08:00 UT) with a flux of 1.23 ± 0.03. b) May 17 (7:10–09:00 UT) with a flux of 1.26 ± 0.03. c) May 17 8–10 h (8:14–09:58 UT) with a flux of 1.31 ± 0.03. d) May 17 9–11 h (9:11–10:56 UT) with a flux of 1.41 ± 0.03. e) May 17 10–12 h (10:09–12:19 UT) with a flux of 1.40 ± 0.03. All maps are restored with a beam of 2.74 × 1.12 at 1.76°. The plotted contour levels are 1.73%, 3.46%, 6.93%, 13.9%, 27.7%, and 55.4%.

In the text
thumbnail Fig. 5

2 h model fits of Sgr A* on May 17 2012. Extended morphology during flare is best fitted by a two component model. a) May 17 6–8 h (6:04–08:00 UT), flux: 1.21 ± 0.03. b) May 17 (7:10–09:00 UT), flux: 1.21 ± 0.03. c) May 17 8–10 h (8:14–09:58 UT), flux: 1.26 ± 0.03. d) May 17 9–11 h (9:11–10:56 UT), flux: 1.29 ± 0.03. e) May 17 10–12 h (10:09–12:19 UT), flux: 1.28 ± 0.03. All maps are restored with a beam of 2.74 × 1.12 at 1.76°. The plotted contour levels are 1.73%, 3.46%, 6.93%, 13.9%, 27.7% and 55.4%.

In the text
thumbnail Fig. 6

RCP map of Sgr A* on May 17, 2012 (8:00–10:00 h UT). The map was convolved with a beam of 2.74 × 1.12 at 1.76°. Contour levels are 1.73%, 3.46%, 6.93%, 13.9%, 27.7%, and 55.4% of the peak flux density of 1.5 Jy/beam.

In the text
thumbnail Fig. 7

Closure phases of all closure triangles on May 17 6:00–12:00 h UT with plotted standard deviation levels. a) FD-KP-PT, b) FD-KP-LA, c) FD-LA-PT, d) KP-LA-PT. We show standard errors provided by standard deviation.

In the text
thumbnail Fig. 8

Radial plot of the models obtained during the time of highest flaring activity of Sgr A* on May 17, 2012. The colors correspond to the maps of 8:00–10:00 h UT (red), 9:00–011:00 h UT (orange), and 10:00–12:00 h UT (blue).

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.