Showing 1–50 of 140 results for author: Wardle, M

Multiwavelength Observations of Sgr A*. II. 2019 July 21 and 26

Authors: Joseph M. Michail, Farhad Yusef-Zadeh, Mark Wardle, Devaky Kunneriath, Joseph L. Hora, Howard Bushouse, Giovanni G. Fazio, Sera Markoff, Howard A. Smith

Abstract: We report on the final two days of a multiwavelength campaign of Sgr A* observing in the radio, submillimeter, infrared, and X-ray bands in July 2019. Sgr A* was remarkably active, showing multiple flaring events across the electromagnetic spectrum. We detect a transient $\sim35$-minute periodicity feature in Spitzer Space Telescope light curves on 21 July 2019. Time-delayed emission was detected… ▽ More We report on the final two days of a multiwavelength campaign of Sgr A* observing in the radio, submillimeter, infrared, and X-ray bands in July 2019. Sgr A* was remarkably active, showing multiple flaring events across the electromagnetic spectrum. We detect a transient $\sim35$-minute periodicity feature in Spitzer Space Telescope light curves on 21 July 2019. Time-delayed emission was detected in ALMA light curves, suggesting a hotspot within the accretion flow on a stable orbit. On the same night, we observe a decreased flux in the submillimeter light curve following an X-ray flare detected by the Chandra X-ray Observatory and model the feature with an adiabatically expanding synchrotron hotspot occulting the accretion flow. The event is produced by a plasma $0.55~R_{\text{S}}$ in radius with an electron spectrum $p=2.84$. It is threaded by a $\sim130$ Gauss magnetic field and expands at $0.6\%$ the speed of light. Finally, we reveal an unambiguous flare in the infrared, submillimeter, and radio, demonstrating that the variable emission is intrinsically linked. We jointly fit the radio and submillimeter light curves using an adiabatically expanding synchrotron hotspot and find it is produced by a plasma with an electron spectrum $p=0.59$, $187$ Gauss magnetic field, and radius $0.47~R_{\text{S}}$ that expands at $0.029c$. In both cases, the uncertainty in the appropriate lower and upper electron energy bounds may inflate the derived equipartition field strengths by a factor of 2 or more. Our results confirm that both synchrotron- and adiabatic-cooling processes are involved in the variable emission's evolution at submillimeter and infrared wavelengths. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 24 pages, 13 figures, accepted to The Astrophysical Journal. Comments welcome! Paper I can be found here: arXiv:2107.09681

The interaction of a large-scale nuclear wind with the high velocity HII region G0.17+0.15

Authors: F. Yusef-Zadeh, Jun-Hui Zhao, R. Arendt, M. Wardle, M. Royster, L. Rudnick, J. Michail

Abstract: We investigate the nature of a Galactic center source, G0.17+0.15, lying along the northern extension of the Radio Arc near l~0.2deg. G0.17+0.15 is an HII region located toward the eastern edge of the radio bubble, embedded within the highly polarized Galactic center eastern Lobe where a number of radio filaments appear to cross through the HII region. We report the detection of hydrogen and heliu… ▽ More We investigate the nature of a Galactic center source, G0.17+0.15, lying along the northern extension of the Radio Arc near l~0.2deg. G0.17+0.15 is an HII region located toward the eastern edge of the radio bubble, embedded within the highly polarized Galactic center eastern Lobe where a number of radio filaments appear to cross through the HII region. We report the detection of hydrogen and helium recombination lines with a radial velocity exceeding 140 km/s based on GBT and VLA observations. The morphology of G0.17+0.15, aided by kinematics, and spectral index characteristics, suggests the presence of an external pressure dragging and shredding the ionized gas. We argue that this ionized cloud is interacting with a bundle of radio filaments and is entrained by the ram pressure of the radio bubble, which itself is thought to be produced by cosmic-ray driven outflows at the Galactic center. In this interpretation, the gas streamers on the western side of G0.17+0.15 are stripped, accelerated from 0 to deltav~35 km/s, over a time scale roughly 8x10^4 years, implying that ablating ram pressure is ~700 eV cm-3, comparable to the ~10^3 eV cm-3 cosmic-ray driven wind pressure in the Galactic center region. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: 23 pages, 12 figures, MNRAS (in press)

G359.13142-0.20005: A steep spectrum radio pulsar candidate with an X-ray counterpart running into the Galactic Center Snake (G359.1-0.2)

Authors: F. Yusef-Zadeh, Jun-Hui Zhao, R. Arendt, M. Wardle, C. O. Heinke, M. Royster, C. Lang, J. Michail

Abstract: The Snake is a remarkable Galactic center radio filament with a morphology characterized by two kinks along its $\sim 20'$ extent. The major and minor kinks are located where the filament is most distorted from a linear magnetized structure running perpendicular to the Galactic plane. We present {\em Chandra}, VLA, and MeerKAT data and report the detection of an X-ray and radio source at the locat… ▽ More The Snake is a remarkable Galactic center radio filament with a morphology characterized by two kinks along its $\sim 20'$ extent. The major and minor kinks are located where the filament is most distorted from a linear magnetized structure running perpendicular to the Galactic plane. We present {\em Chandra}, VLA, and MeerKAT data and report the detection of an X-ray and radio source at the location of the major kink. High-resolution radio images of the major kink reveal a compact source with a steep spectrum with spectral index alpha ~ -2.7 surrounded by extended emission. The radio luminosity and steep spectrum of the compact source are consistent with a pulsar. We also show flattening of the spectrum and enhanced synchrotron emissivity away from the position of the major kink along the Snake, which suggests injection of relativistic particles along the Snake. We argue that the major kink is created by a fast-moving (~500-1000 km/s), object punching into the Snake, distorting its magnetic structure, and producing X-ray emission. X-ray emission pinpoints an active acceleration site where the interaction is taking place. A secondary kink is argued to be induced by the impact of the high-velocity object producing the major kink. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: 12 pages, 8 figures, MNRAS (in press)

Detection of large-scale synchrotron radiation from the molecular envelope of the Sgr B cloud complex at the Galactic center

Authors: F. Yusef-Zadeh, M. Wardle, R. Arendt, J. W. Hewitt, Y. Hu, A. Lazarian, N. Kassim, S. Hyman, I. Heywood

Abstract: We present highly sensitive measurements taken with MeerKAT at 1280 MHz as well as archival GBT, MWA and VLA images at 333, 88 and 74 MHz. We report the detection of synchrotron radio emission from the infrared dark cloud (IRDC) associated with the halo of the Sgr B complex on a scale of ~60 pc. A strong spatial correlation between low-frequency radio continuum emission and dense molecular gas, co… ▽ More We present highly sensitive measurements taken with MeerKAT at 1280 MHz as well as archival GBT, MWA and VLA images at 333, 88 and 74 MHz. We report the detection of synchrotron radio emission from the infrared dark cloud (IRDC) associated with the halo of the Sgr B complex on a scale of ~60 pc. A strong spatial correlation between low-frequency radio continuum emission and dense molecular gas, combined with spectral index measurements, indicates enhanced synchrotron emission by cosmic-ray electrons. Correlation of the FeI 6.4 keV Kalpha line and synchrotron emission provides compelling evidence that the low energy cosmic-ray electrons are responsible for producing the Kalpha line emission. The observed synchrotron emission within the halo of the Sgr B cloud complex has mean spectral index alpha -1+/-1 gives the magnetic field strength ~100 muG for cloud densities nH = 10^4-10^5 cm-3 and estimate cosmic-ray ionization rates between 10^-13 and 10^-14 s^-1. Furthermore, the energy spectrum of primary cosmic-ray electrons is constrained to be E^-3 +/-1 for typical energies of few hundred MeV. The extrapolation of this spectrum to higher energies is consistent with X-ray and gamma-ray emission detected from this cloud. These measurements have important implications on the role that high cosmic-ray electron fluxes at the Galactic center play in production of radio synchrotron emission, the FeI Kalpha line emission at 6.4 keV and ~GeV gamma-ray emission throughout the central molecular zone (CMZ). △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 10 pages, 4 figures, MN (in press)

The Population of the Galactic Center Filaments: Position Angle Distribution Reveal a Degree-scale Collimated Outflow from Sgr A* along the Galactic Plane

Authors: F. Yusef-Zadeh, R. G. Arendt, M. Wardle, I. Heywood

Abstract: We have examined the distribution of the position angle (PA) of the Galactic center filaments with lengths $L > 66''$ and $ < 66''$ as well as their length distribution as a function of PA. We find bimodal PA distributions of the filaments, long and short populations of radio filaments. Our PA study shows the evidence for a distinct population of short filaments with PA close to the Galactic plane… ▽ More We have examined the distribution of the position angle (PA) of the Galactic center filaments with lengths $L > 66''$ and $ < 66''$ as well as their length distribution as a function of PA. We find bimodal PA distributions of the filaments, long and short populations of radio filaments. Our PA study shows the evidence for a distinct population of short filaments with PA close to the Galactic plane. Mainly thermal short radio filaments ($<66''$) have PAs concentrated close to the Galactic plane within $60^\circ < \rm PA <120^\circ$. Remarkably, the short filament PAs are radial with respect to the Galactic center at $l <0^\circ$, and extend in the direction toward Sgr A*. On a smaller scale, the prominent Sgr E HII complex G358.7-0.0 provides a vivid example of the nearly radial distribution of short filaments. The bimodal PA distribution suggests different origin for two distinct filament populations. We argue that alignment of the short filament population results from the ram pressure of a degree-scale outflow from Sgr A* that exceeds the internal filament pressure, and aligns them along the Galactic plane. The ram pressure is estimated to be 2$\times10^6\,$ cm$^{-3}\,$ K at a distance of 300pc, requiring biconical mass outflow rate $10^{-4}$ \msol\, yr$^{-1}$ with an opening angle of $\sim40^\circ$. This outflow aligns not only the magnetized filaments along the Galactic plane but also accelerates thermal material associated with embedded or partially embedded clouds. This places an estimate of $\sim$6 Myr as the age of the outflow. △ Less

Submitted 1 June, 2023; originally announced June 2023.

Comments: 19 pages, 8 figures, ApJL (June 2nd, 2023)

GNOMES II: Analysis of the Galactic diffuse molecular ISM in all four ground state hydroxyl transitions using Amoeba

Authors: Anita Petzler, J. R. Dawson, Hiep Nguyen, Carl Heiles, M. Wardle, M. -Y. Lee, Claire E. Murray, K. L. Thompson, Snezana Stanimirovic

Abstract: We present observations of the four 2 Pi 3/2 J = 3/2 ground-rotational state transitions of the hydroxyl molecule (OH) along 107 lines of sight both in and out of the Galactic plane: 92 sets of observations from the Arecibo telescope and 15 sets of observations from the Australia Telescope Compact Array (ATCA). Our Arecibo observations included off-source pointings, allowing us to measure excitati… ▽ More We present observations of the four 2 Pi 3/2 J = 3/2 ground-rotational state transitions of the hydroxyl molecule (OH) along 107 lines of sight both in and out of the Galactic plane: 92 sets of observations from the Arecibo telescope and 15 sets of observations from the Australia Telescope Compact Array (ATCA). Our Arecibo observations included off-source pointings, allowing us to measure excitation temperature (Tex) and optical depth, while our ATCA observations give optical depth only. We perform Gaussian decomposition using the Automated Molecular Excitation Bayesian line-fitting Algorithm 'AMOEBA' (Petzler, Dawson, and Wardle 2021) fitting all four transitions simultaneously with shared centroid velocity and width. We identify 109 features across 38 sightlines (including 58 detections along 27 sightlines with excitation temperature measurements). While the main lines at 1665 and 1667 MHz tend to have similar excitation temperatures (median Tex(main) difference = 0.6 K, 84% show Tex(main) difference < 2 K), large differences in the 1612 and 1720 MHz satellite line excitation temperatures show that the gas is generally not in LTE. For a selection of sightlines we compare our OH features to associated (on-sky and in velocity) HI cold gas components (CNM) identified by Nguyen et al. (2019) and find no strong correlations. We speculate that this may indicate an effective decoupling of the molecular gas from the CNM once it accumulates. △ Less

Submitted 21 February, 2023; originally announced February 2023.

Comments: Accepted for publication to PASA. 41 pages, 27 figures

Polarized signatures of adiabatically expanding hotspots in Sgr A*'s accretion flow

Authors: Joseph M. Michail, Farhad Yusef-Zadeh, Mark Wardle, Devaky Kunneriath

Abstract: We report 235 GHz linear and circular polarization (LP and CP) detections of Sgr A* at levels of $\sim10\%$ and $\sim-1\%$, respectively, using ALMA. We describe the first full-Stokes modeling of an observed submillimeter flare with an adiabatically-expanding synchrotron hotspot using a polarized radiative transfer prescription. Augmented with a simple full-Stokes model for the quiescent emission,… ▽ More We report 235 GHz linear and circular polarization (LP and CP) detections of Sgr A* at levels of $\sim10\%$ and $\sim-1\%$, respectively, using ALMA. We describe the first full-Stokes modeling of an observed submillimeter flare with an adiabatically-expanding synchrotron hotspot using a polarized radiative transfer prescription. Augmented with a simple full-Stokes model for the quiescent emission, we jointly characterize properties of both the quiescent and variable components by simultaneously fitting all four Stokes parameter light curves. The hotspot has magnetic field strength $71$ G, radius $0.75$ Schwarzschild radii, and expands at speed $0.013$c assuming magnetic equipartition. The magnetic field's position angle projected in the plane-of-sky is $\approx55^\circ$ East of North, which previous analyses reveal as the accretion flow's angular momentum axis and further supports Sgr A* hosting a magnetically-arrested disk. The magnetic field is oriented approximately perpendicular to the line of sight, which suggests repolarization as the cause of the high circular-to-linear polarization ratio observed at radio frequencies. We additionally recover several properties of the quiescent emission, consistent with previous analyses of the accretion flow, such as a rotation measure $\approx-4.22\times10^{5}$ rad m$^{-2}$. Our findings provide critical constraints for interpreting and mitigating the polarized variable emission in future Event Horizon Telescope images of Sgr A*. △ Less

Submitted 24 January, 2023; originally announced January 2023.

Comments: 13 pages, 6 figures, 3 tables. Accepted to Monthly Notices of the Royal Astronomical Society, comments welcome!

Populations of Magnetized Filaments in the Intracluster Medium and the Galactic Center

Authors: F. Yusef-Zadeh, R. G. Arendt, M. Wardle

Abstract: Magnetized radio filaments are found in abundance in the inner few hundred pc of our Galaxy. Progress in understanding this population of filaments has been slow, in part due to a lack of detection elsewhere in the Galaxy or in external galaxies. Recent highly sensitive radio continuum observations of radio galaxies in galaxy clusters have revealed remarkable isolated filamentary structures in the… ▽ More Magnetized radio filaments are found in abundance in the inner few hundred pc of our Galaxy. Progress in understanding this population of filaments has been slow, in part due to a lack of detection elsewhere in the Galaxy or in external galaxies. Recent highly sensitive radio continuum observations of radio galaxies in galaxy clusters have revealed remarkable isolated filamentary structures in the ICM that are linked to radio jets, tails and lobes. The origin of this class of filaments is not understood either. Here, we argue that the underlying physical mechanisms responsible for the creation of the two populations are the same because of their similarity in morphology, spacing between the filaments, aspect ratio, magnetic energy densities to thermal pressure of the medium, and that both populations have undergone synchrotron aging. These similarities provide an opportunity to investigate the physical processes in the ISM and ICM for the first time. We consider that the origin of the filaments in both the GC and ICM is a result of the interaction of a large-scale wind with clouds, or the filaments arise through the stretching and collection of field lines by turbulence in weakly magnetized medium. We examine these ideas toward four radio galaxy filaments associated with four radio galaxies IC 40B, IC 4496, J1333-3141, ESO137-006 and argue that much can be understood in the future by comparing these two populations of filaments. △ Less

Submitted 10 October, 2022; originally announced October 2022.

Comments: 14 pages, 4 figures, ApJL (in press)

The population of Galactic centre filaments III: candidate radio and stellar sources

Authors: F. Yusef-Zadeh, R. G. Arendt, M. Wardle, I. Heywood, W. Cotton

Abstract: Recent MeerKAT radio continuum observations of the Galactic center at 20 cm show a large population of nonthermal radio filaments (NRFs) in the inner few hundred pc of the Galaxy. We have selected a sample of 57 radios ources, mainly compact objects, in the MeerKAT mosaic image that appear to be associated with NRFs. The selected sources are about 4 times the number of radio point sources associat… ▽ More Recent MeerKAT radio continuum observations of the Galactic center at 20 cm show a large population of nonthermal radio filaments (NRFs) in the inner few hundred pc of the Galaxy. We have selected a sample of 57 radios ources, mainly compact objects, in the MeerKAT mosaic image that appear to be associated with NRFs. The selected sources are about 4 times the number of radio point sources associated with filaments than would be expected by random chance. Furthermore, an apparent correlation between bright IR stars and NRFs is inferred from their similar latitude distributions, suggesting that they both co-exist within the same region. To examine if compact radio sources are related to compact IR sources, we have used archival 2MASS, and {\em Spitzer} data to make spectral energy distribution of individual stellar sources coincident or close to radio sources. We provide a catalogue of radio and IR sources for future detailed observations to investigate a potential 3-way physical association between NRFs, compact radio and IR stellar sources. This association is suggested by models in which NRFs are cometary tails produced by the interaction of a large-scale nuclear outflow with stellar wind bubbles in the Galactic center. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Comments: 63 pages, 3 Figures (multiple subfigures), MNRAS (in press)

The structure of jets launched from post-AGB binary systems

Authors: Dylan Bollen, Devika Kamath, Hans Van Winckel, Orsola De Marco, Olivier Verhamme, Jacques Kluska, Mark Wardle

Abstract: We focus on post-asymptotic giant branch (post-AGB) binaries and study the interaction between the different components of these complex systems. These components comprise the post-AGB primary, a main sequence secondary, a circumbinary disk, as well as a fast bipolar outflow (jet) launched by the companion. We obtained well-sampled time series of high resolution optical spectra over the last decad… ▽ More We focus on post-asymptotic giant branch (post-AGB) binaries and study the interaction between the different components of these complex systems. These components comprise the post-AGB primary, a main sequence secondary, a circumbinary disk, as well as a fast bipolar outflow (jet) launched by the companion. We obtained well-sampled time series of high resolution optical spectra over the last decade and these spectra provide the basis of our study. The jet is detected in absorption, at superior conjunction, when the line of sight towards the primary goes through the bipolar cone. Our spectral time series scan the jets during orbital motion. Our spatio-kinematic model is constrained by these dynamical spectra. We complement this with a radiative-transfer model in which the Balmer series are used to derive total mass-loss rates in the jets. The jets are found to be wide and display an angle-dependent density structure with a dense and slower outer region near the jet cone and a fast inner part along the jet symmetry axes. The deprojected outflow velocities confirm that the companions are main sequence companions. The total mass-loss rates are large (10^{-8} and 10^{-5}\,solar mass per year), from which we can infer that the mass-accretion rates onto the companion star must be high as well. The circumbinary disk is likely the main source for the accretion disk around the companion. All systems with full disks that start near the sublimation radius show jets, whereas for systems with evolved transition disks, this lowers to a detection rate of 50%. Objects without an infrared excess do not show jets. We conclude that jet creation in post-AGB binaries is a mainstream process. The interaction between the circumbinary disks and the central binary provide the needed accretion flow, but the presence of a circumbinary disk does not seem to be the only prerequisite to launch a jet. △ Less

Submitted 18 August, 2022; originally announced August 2022.

Comments: 17 pages, 12 figures, accepted by Astronomy and Astrophysics

Journal ref: A&A 666, A40 (2022)

Statistical Properties of the Population of the Galactic Center Filaments II: The Spacing Between Filaments

Authors: F. Yusef-Zadeh, R. G. Arendt, M. Wardle, S. Boldyrev, I. Heywood, W. Cotton, F. Camilo

Abstract: We carry out a population study of magnetized radio filaments in the Galactic center using MeerKAT data by focusing on the spacing between the filaments that are grouped. The morphology of a sample of 43 groupings containing 174 magnetized radio filaments are presented. Many grouped filaments show harp-like, fragmented cometary tail-like, or loop-like structures in contrast to many straight filame… ▽ More We carry out a population study of magnetized radio filaments in the Galactic center using MeerKAT data by focusing on the spacing between the filaments that are grouped. The morphology of a sample of 43 groupings containing 174 magnetized radio filaments are presented. Many grouped filaments show harp-like, fragmented cometary tail-like, or loop-like structures in contrast to many straight filaments running mainly perpendicular to the Galactic plane. There are many striking examples of a single filament splitting into two prongs at a junction, suggestive of a flow of plasma along the filaments. Spatial variations in spectral index, brightness, bending and sharpening along the filaments indicate that they are evolving on a 10^{5-6}-year time scale. The mean spacings between parallel filaments in a given grouping peaks at $\sim16''$. We argue by modeling that the filaments in a grouping all lie on the same plane and that the groupings are isotropically oriented in 3D space. One candidate for the origin of filamentation is interaction with an obstacle, which could be a compact radio source, before a filament splits and bends into multiple filaments. In this picture, the obstacle or sets the length scale of the separation between the filaments. Another possibility is synchrotron cooling instability occurring in cometary tails formed as a result of the interaction of cosmic-ray driven Galactic center outflow with obstacles such as stellar winds. In this picture, the mean spacing and the mean width of the filaments are expected to be a fraction of a parsec, consistent with observed spacing. △ Less

Submitted 21 June, 2022; originally announced June 2022.

Comments: 53 pages, 32 figures, MNRS (in press)

The non-ideal finite Larmor radius effect in the solar atmosphere

Authors: B. P. Pandey, Mark Wardle

Abstract: The dynamics of the partially ionized solar atmosphere is controlled by the frequent collision and charge exchange between the predominant neutral Hydrogen atoms and charged ions. At signal frequencies below or of the order of either of the collision or charge exchange frequencies the magnetic stress is {\it felt} by both the charged and neutral particles simultaneously. The resulting neutral-mass… ▽ More The dynamics of the partially ionized solar atmosphere is controlled by the frequent collision and charge exchange between the predominant neutral Hydrogen atoms and charged ions. At signal frequencies below or of the order of either of the collision or charge exchange frequencies the magnetic stress is {\it felt} by both the charged and neutral particles simultaneously. The resulting neutral-mass loading of the ions leads to the rescaling of the effective ion-cyclotron frequency-it becomes the Hall frequency, and the resultant effective Larmor radius becomes of the order of few kms. Thus the finite Larmor radius (FLR) effect which manifests as the ion and neutral pressure stress tensors operates over macroscopic scales. Whereas parallel and perpendicular (with respect to the magnetic field) viscous momentum transport competes with the Ohm and Hall diffusion of the magnetic field in the photosphere-chromosphre, the gyroviscous effect becomes important only in the transition region between the chromosphere and corona, where it competes with the ambipolar diffusion. The wave propagation in the gyroviscous effect dominated medium depends on the plasma $β$ (a ratio of the thermal and magnetic energies). The abundance of free energy makes gyro waves unstable with the onset condition exactly opposite of the Hall instability. However, the maximum growth rate is identical to the Hall instability. For a flow gradient $\sim 0.1 \,\mbox{s}^{-1}$ the instability growth time is one minute. Thus, the transition region may become subject to this fast growing, gyroviscous instability. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: 17 pages, 12 figures, MNRAS (in press)

Statistical Properties of the Population of the Galactic Center Filaments: The Spectral Index and Equipartition Magnetic Field

Authors: F. Yusef-Zadeh, R. G. Arendt, M. Wardle, I. Heywood, W. D. Cotton, F. Camilo

Abstract: We present high-pass filtered continuum images of the inner $3.5^\circ\times2.5^\circ$ of the Galactic center at 20 cm with $6.4''$ resolution. These mosaic images are taken with MeerKAT and reveal a large number of narrow filaments, roughly an order of magnitude increase in their numbers compared to past measurements. For the first time, we carry out population studies of the spectral index and m… ▽ More We present high-pass filtered continuum images of the inner $3.5^\circ\times2.5^\circ$ of the Galactic center at 20 cm with $6.4''$ resolution. These mosaic images are taken with MeerKAT and reveal a large number of narrow filaments, roughly an order of magnitude increase in their numbers compared to past measurements. For the first time, we carry out population studies of the spectral index and magnetic field of the entire region. The mean spectral indices of the filaments are steeper than supernova remnants (SNRs) (-0.62) with a value of $α\sim-0.83$. The variation in $α$ is much larger than for the SNRs, suggesting that these characteristics have a different origin. A large-scale cosmic-ray driven wind has recently been proposed to explain the origin of filaments and the large-scale 430 pc bipolar radio and X-ray structure. This favors the possibility that the large-scale bipolar radio/X-ray structure is produced by past activity of Sgr A* rather than coordinated burst of supernovae. A trend of steeper indices is also noted with increasing distance from the Galactic plane. This could be explained either by synchrotron cooling or weak shocks accelerating cosmic-ray particles in the context of the cosmic-ray driven wind. The mean magnetic field strengths along the filaments ranges from $\sim100$ to 400 $μ$G depending on the assumed ratio of cosmic-ray protons to electrons. Given that there is a high cosmic ray pressure in the Galactic center, the large equipartition magnetic field implies that the magnetic field is weak in most of the interstellar volume of the Galactic center. △ Less

Submitted 25 January, 2022; originally announced January 2022.

Comments: 15 pages, 9 figures, accepted for publication in Astrophysical Journal Letters

AMOEBA: Automated Molecular Excitation Bayesian Line-Fitting Algorithm

Authors: Anita Petzler, Joanne R Dawson, Mark Wardle

Abstract: The hyperfine transitions of the ground-rotational state of the hydroxyl radical (OH) have emerged as a versatile tracer of the diffuse molecular interstellar medium. We present a novel automated Gaussian decomposition algorithm designed specifically for the analysis of the paired on-source and off-source optical depth and emission spectra of these transitions. In contrast to existing automated Ga… ▽ More The hyperfine transitions of the ground-rotational state of the hydroxyl radical (OH) have emerged as a versatile tracer of the diffuse molecular interstellar medium. We present a novel automated Gaussian decomposition algorithm designed specifically for the analysis of the paired on-source and off-source optical depth and emission spectra of these transitions. In contrast to existing automated Gaussian decomposition algorithms, AMOEBA (Automated MOlecular Excitation Bayesian line-fitting Algorithm) employs a Bayesian approach to model selection, fitting all 4 optical depth and 4 emission spectra simultaneously. AMOEBA assumes that a given spectral feature can be described by a single centroid velocity and full width at half-maximum, with peak values in the individual optical depth and emission spectra then described uniquely by the column density in each of the four levels of the ground-rotational state, thus naturally including the real physical constraints on these parameters. Additionally, the Bayesian approach includes informed priors on individual parameters which the user can modify to suit different data sets. Here we describe AMOEBA and evaluate its validity and reliability in identifying and fitting synthetic spectra with known parameters. △ Less

Submitted 25 July, 2021; originally announced July 2021.

Comments: Submitted to ApJ

Multi-wavelength Observations of Sgr A*. I. 2019 July 18

Authors: Joseph M. Michail, Mark Wardle, Farhad Yusef-Zadeh, Devaky Kunneriath

Abstract: We present and analyze ALMA submillimeter observations from a multi-wavelength campaign of Sgr A* during 18 July 2019. In addition to the submillimeter, we utilize concurrent mid-IR (Spitzer) and X-ray (Chandra) observations. The submillimeter emission lags less than $δt\approx30$ minutes behind the mid-IR data. However, the entire submillimeter flare was not observed, raising the possibility that… ▽ More We present and analyze ALMA submillimeter observations from a multi-wavelength campaign of Sgr A* during 18 July 2019. In addition to the submillimeter, we utilize concurrent mid-IR (Spitzer) and X-ray (Chandra) observations. The submillimeter emission lags less than $δt\approx30$ minutes behind the mid-IR data. However, the entire submillimeter flare was not observed, raising the possibility that the time delay is a consequence of incomplete sampling of the light curve. The decay of the submillimeter emission is not consistent with synchrotron cooling. Therefore, we analyze these data adopting an adiabatically expanding synchrotron source that is initially optically thick or thin in the submillimeter, yielding time-delayed or synchronous flaring with the IR, respectively. The time-delayed model is consistent with a plasma blob of radius $0.8~R_{\text{S}}$ (Schwarzschild radius), electron power-law index $p=3.5$ ($N(E)\propto E^{-p}$), equipartition magnetic field of $B_{\text{eq}}\approx90$ Gauss, and expansion velocity $v_{\text{exp}}\approx0.004c$. The simultaneous emission is fit by a plasma blob of radius $2~R_{\text{S}}$, $p=2.5$, $B_{\text{eq}}\approx27$ Gauss, and $v_{\text{exp}}\approx0.014c$. Since the submillimeter time delay is not completely unambiguous, we cannot definitively conclude which model better represents the data. This observation presents the best evidence for a unified flaring mechanism between submillimeter and X-ray wavelengths and places significant constraints on the source size and magnetic field strength. We show that concurrent observations at lower frequencies would be able to determine if the flaring emission is initially optically thick or thin in the submillimeter. △ Less

Submitted 3 November, 2021; v1 submitted 20 July, 2021; originally announced July 2021.

Comments: 9 pages, 5 figures, Accepted to ApJ. Major updates since original submitted version: 1. Correction of mismatching time coordinate systems, 2. Added text discussing time delay error bars, 3. Made text more model-agnostic with respect to time delayed or simultaneous emission. Second revision edits are grammatic in style with updated references

Detection of a 20 minute time lag observed from Sgr A* between 8 and 10 GHz with the VLA

Authors: Joseph M. Michail, Farhad Yusef-Zadeh, Mark Wardle

Abstract: We report the detection and analysis of a radio flare observed on 17 April 2014 from Sgr A* at $9$ GHz using the VLA in its A-array configuration. This is the first reported simultaneous radio observation of Sgr A* across $16$ frequency windows between $8$ and $10$ GHz. We cross correlate the lowest and highest spectral windows centered at $8.0$ and $9.9$ GHz, respectively, and find the $8.0$ GHz… ▽ More We report the detection and analysis of a radio flare observed on 17 April 2014 from Sgr A* at $9$ GHz using the VLA in its A-array configuration. This is the first reported simultaneous radio observation of Sgr A* across $16$ frequency windows between $8$ and $10$ GHz. We cross correlate the lowest and highest spectral windows centered at $8.0$ and $9.9$ GHz, respectively, and find the $8.0$ GHz light curve lagging $18.37^{+2.17}_{-2.18}$ minutes behind the $9.9$ GHz light curve. This is the first time lag found in Sgr A*'s light curve across a narrow radio frequency bandwidth. We separate the quiescent and flaring components of Sgr A* via flux offsets at each spectral window. The emission is consistent with an adiabatically-expanding synchrotron plasma, which we fit to the light curves to characterize the two components. The flaring emission has an equipartition magnetic field strength of $2.2$ Gauss, size of $14$ Schwarzschild radii, average speed of $12000$ km s$^{-1}$, and electron energy spectrum index ($N(E)\propto E^{-p}$), $p = 0.18$. The peak flare flux at $10$ GHz is approximately $25$% of the quiescent emission. This flare is abnormal as the inferred magnetic field strength and size are typically about $10$ Gauss and few Schwarzschild radii. The properties of this flare are consistent with a transient warm spot in the accretion flow at a distance of $10$-$100$ Schwarzschild radii from Sgr A*. Our analysis allows for independent characterization of the variable and quiescent components, which is significant for studying temporal variations in these components. △ Less

Submitted 24 May, 2021; originally announced May 2021.

Comments: 9 pages, 6 figures, accepted to the Monthly Notices of the Royal Astronomical Society

Revealing the Physical Conditions around Sgr A* using Bayesian Inference -- I. Observations and Radiative Transfer

Authors: Tomas A. James, Serena Viti, Farhad Yusef-Zadeh, Marc Royster, Mark Wardle

Abstract: We report sub-arcsecond ALMA observations between 272 - 375 GHz towards Sgr A*'s Circumnuclear disk (CND). Our data comprises 8 individual pointings, with significant SiO (8(7) - 7(6)) and SO (7 - 6) emission detected towards 98 positions within these pointings. Additionally, we identify H2CS (9(1,9) - 8(1,8)), OCS (25 - 24) and CH3OH (2(1,1) - 2(0,2)) towards a smaller subset of positions. By usi… ▽ More We report sub-arcsecond ALMA observations between 272 - 375 GHz towards Sgr A*'s Circumnuclear disk (CND). Our data comprises 8 individual pointings, with significant SiO (8(7) - 7(6)) and SO (7 - 6) emission detected towards 98 positions within these pointings. Additionally, we identify H2CS (9(1,9) - 8(1,8)), OCS (25 - 24) and CH3OH (2(1,1) - 2(0,2)) towards a smaller subset of positions. By using the observed peak line flux density together with a Bayesian Inference technique informed by radiative transfer models, we systematically recover the physical gas conditions towards each of these positions. We estimate that the bulk of the surveyed gas has temperature T < 500 K and density n $\lessapprox 10^{6}$ cm$^{-3}$, consistent with previous studies of similar positions as traced by HCN clumps. However, we identify an uncharacteristically hot (T $\approx 600$ K) and dense (n $\approx 10^{6}$ cm$^{-3}$) source in the Northeastern Arm. This position is found to be approximately consistent with a gravitationally bound region dominated by turbulence. We also identify a nearby cold (T $\approx 60$ K) and extremely dense (n $\approx 10^{7}$ cm$^{-3}$) position that is again potentially bound and dominated by turbulence. We also determine that the total gas mass contained within the CND is M $\approx 4 \times 10^{4}$ $M_{\odot}$. Furthermore, we qualitatively note that the observed chemical enrichment across large scales within the CND is consistent with bulk grain processing, though multiple desorption mechanisms are plausibly responsible. Further chemical modelling is required to identify the physical origin of the grain-processing, as well as the localised H2CS and OCS emission. △ Less

Submitted 5 May, 2021; v1 submitted 30 April, 2021; originally announced May 2021.

Comments: Accepted by ApJ

Jet parameters for a diverse sample of jet-launching post-AGB binaries

Authors: Dylan Bollen, Devika Kamath, Hans Van Winckel, Orsola De Marco, Mark Wardle

Abstract: Jets are a commonly observed phenomenon in post-asymptotic giant branch (post-AGB) binaries. Due to the orbital motion of the binary, the jet causes variable absorption in the Balmer profiles. In previous work, we have developed spatio-kinematic and radiative transfer models to reproduce the observed Balmer line variability and derive the spatio-kinematic structure of the jet and its mass-loss rat… ▽ More Jets are a commonly observed phenomenon in post-asymptotic giant branch (post-AGB) binaries. Due to the orbital motion of the binary, the jet causes variable absorption in the Balmer profiles. In previous work, we have developed spatio-kinematic and radiative transfer models to reproduce the observed Balmer line variability and derive the spatio-kinematic structure of the jet and its mass-loss rate. Here, we apply our jet model to five post-AGB binaries with distinct Hα line variability and diverse orbital properties. Our models fit the Hα line variations very well. We estimate jet mass-loss rates between 10-8 Mdot yr-1 and 10-4 Mdot yr-1 , from which we deduce accretion rates onto the companion between 10-7 Mdot yr-1 and 10-3 Mdot yr-1 . These accretion rates are somewhat higher than can be comfortably explained with reasonable sources of accretion, but we argue that the circumbinary disc in these systems is most-likely the source feeding the accretion, although accretion from the post-AGB star cannot be ruled out. The diversity of the variability in the five objects is due to their wide ejection cones combined with a range of viewing angles, rather than inherent differences between the objects. The nature of the observations does not let us easily distinguish which jet launching model (stellar jet, disc wind, or X-wind) should be favoured. In conclusion, we show that our jet model includes the physical parameters to successfully reproduce the Hα line variations and retrieve the structure and mass-loss rates of the jet for all five objects that are representative of the diverse sample of Galactic post-AGB binaries. △ Less

Submitted 6 January, 2021; originally announced January 2021.

Comments: 20 pages, 9 figures, accepted for publication in MNRAS

G0.173-0.42: an X-ray and radio magnetized filament near the galactic center

Authors: F. Yusef-Zadeh, M. Wardle, C. Heinke, R. Arendt, M. Royster, I. Heywood, W. Cotton, F. Camilo, J. Michail

Abstract: The recent detection of an X-ray filament associated with the radio filament G0.173-0.42 adds to four other nonthermal radio filaments with X-ray counterparts, amongst the more than 100 elongated radio structures that have been identified as synchrotron-emitting radio filaments in the inner couple of degrees of the Galactic center. The synchrotron mechanism has also been proposed to explain the em… ▽ More The recent detection of an X-ray filament associated with the radio filament G0.173-0.42 adds to four other nonthermal radio filaments with X-ray counterparts, amongst the more than 100 elongated radio structures that have been identified as synchrotron-emitting radio filaments in the inner couple of degrees of the Galactic center. The synchrotron mechanism has also been proposed to explain the emission from X-ray filaments. However, the origin of radio filaments and the acceleration sites of energetic particles to produce synchrotron emission in radio and X-rays remain mysterious. Using MeerKAT, VLA, Chandra, WISE and Spitzer, we present structural details of G0.173-0.42 which consists of multiple radio filaments, one of which has an X-ray counterpart. A faint oblique radio filament crosses the radio and X-ray filaments. Based on the morphology, brightening of radio and X-ray intensities, and radio spectral index variation, we argue that a physical interaction is taking place between two magnetized filaments. We consider that the reconnection of the magnetic field lines at the interaction site leads to the acceleration of particles to GeV energies. We also argue against the synchrotron mechanism for the X-ray emission due to the short $\sim$30 year lifetime of TeV relativistic particles. Instead, we propose that the inverse Compton scattering mechanism is more likely to explain the X-ray emission by upscattering of seed photons emitted from a 10^6 solar luminosity star located at the northern tip of the X-ray filament. △ Less

Submitted 19 October, 2020; originally announced October 2020.

Comments: 13 pages, 7 figures, in press (MNRAS)

Evidence for a jet and outflow from Sgr A*: a continuum and spectral line study

Authors: F. Yusef-Zadeh, M. Royster, M. Wardle, W. Cotton, D. Kunneriath, I. Heywood, J. Michail

Abstract: We study the environment of Sgr A* using spectral and continuum observations with the ALMA and VLA. Our analysis of sub-arcsecond H30alpha, H39alpha, H52alpha and H56alpha line emission towards Sgr A* confirm the recently published broad peak ~500 km/s~spectrum toward Sgr~A*. We also detect emission at more extreme radial velocities peaking near -2500 and 4000 km/s, within 0.2''. We then present b… ▽ More We study the environment of Sgr A* using spectral and continuum observations with the ALMA and VLA. Our analysis of sub-arcsecond H30alpha, H39alpha, H52alpha and H56alpha line emission towards Sgr A* confirm the recently published broad peak ~500 km/s~spectrum toward Sgr~A*. We also detect emission at more extreme radial velocities peaking near -2500 and 4000 km/s, within 0.2''. We then present broad band radio continuum images at multiple frequencies on scales from arcseconds to arcminutes. A number of elongated continuum structures lie parallel to the Galactic plane, extending from ~0.4'' to 10'. We note a nonthermal elongated structure on an arcminute scale emanating from Sgr A* at low frequencies between 1 and 1.4 GHz where thermal emission from the mini-spiral is depressed by optical depth effects. The position angle of this elongated structure and the sense of motion of ionized features with respect to Sgr A* suggest a symmetric, collimated jet emerging from Sgr A* with an opening angle of ~30deg and a position angle of ~60deg punching through the medium before accelerating a significant fraction of the orbiting ionized gas to high velocities. The jet with estimated mass flow rate ~1.4x10^{-5} solar mass/yr emerges perpendicular to the equatorial plane of the accretion flow near the event horizon of Sgr A* and runs along the Galactic plane. To explain a number of east-west features near Sgr A*, we also consider the possibility of an outflow component with a wider-angle launched from the accretion flow at larger radii. △ Less

Submitted 10 August, 2020; originally announced August 2020.

Comments: 26 pages, 14 figures, MNRAS, in press

Determining mass accretion and jet mass-loss rates in post-asymptotic giant branch binary systems

Authors: Dylan Bollen, Devika Kamath, Orsola De Marco, Hans Van Winckel, Mark Wardle

Abstract: Aims. In this study, we determine the morphology and mass-loss rate of jets emanating from the companion in post-asymptotic giant branch (post-AGB) binary stars with a circumbinary disk. In doing so, we also determine the mass accretion rates on to the companion and investigate the source feeding the circum-companion accretion disk. Methods. We perform a spatio-kinematic modelling of the jet of tw… ▽ More Aims. In this study, we determine the morphology and mass-loss rate of jets emanating from the companion in post-asymptotic giant branch (post-AGB) binary stars with a circumbinary disk. In doing so, we also determine the mass accretion rates on to the companion and investigate the source feeding the circum-companion accretion disk. Methods. We perform a spatio-kinematic modelling of the jet of two well-sampled post-AGB binaries, BD+46442 and IRAS19135+3937, by fitting the orbital phased time series of H-alpha spectra. Once the jet geometry, velocity and scaled density structure are computed, we carry out radiative transfer modelling of the jet for the first four Balmer lines to determine the jet densities, thus allowing us to compute the jet mass-loss rates and mass accretion rates. Results. The spatio-kinematic model of the jet reproduces the observed absorption feature in the H-alpha lines. In both objects, the jets have an inner region with extremely low density. Using our radiative transfer model, we find the full three-dimensional density structure of both jets. From these results, we can compute mass-loss rates of the jets, which are of the order of 10^-7 - 10^-5 M_sol/yr. We estimate mass accretion rates onto the companion of 10^-6 - 10^-4 M_sol/yr. Conclusions. Based on the mass accretion rates found for these two objects, we conclude that the circumbinary disk is most likely the source feeding the circum-companion accretion disk. This is in agreement with the observed depletion patterns in post-AGB binaries. The high accretion rates from the circumbinary disk imply that the lifetime of the disk will be short. Mass-transfer from the post-AGB star cannot be excluded in these systems, but it is unlikely to provide a sufficient mass-transfer rate to sustain the observed jet mass-loss rates. △ Less

Submitted 28 July, 2020; originally announced July 2020.

Comments: 16 pages, 13 figures, accepted for publication in A&A

Journal ref: A&A 641, A175 (2020)

Cosmic-ray Particles in the Galactic Center: Blowing in the Wind

Authors: F. Yusef-Zadeh, M. Wardle, I. Heywood, W. Cotton, M. Royster

Abstract: Recent results from multi-wavelength observations of the inner few hundred pc of the Galactic center have added two new characteristics to the ISM in this unique region. One is the cosmic ray ionization rate derived from H$_3^+$ measurements is at least two orders of magnitudes higher than in the disk of the Galaxy. The other is the bipolar thermal X-ray and synchrotron emission from this region,… ▽ More Recent results from multi-wavelength observations of the inner few hundred pc of the Galactic center have added two new characteristics to the ISM in this unique region. One is the cosmic ray ionization rate derived from H$_3^+$ measurements is at least two orders of magnitudes higher than in the disk of the Galaxy. The other is the bipolar thermal X-ray and synchrotron emission from this region, suggesting a relic of past activity. We propose that the high cosmic ray pressure drives a large-scale wind away from the Galactic plane and produces the bipolar emission as well as highly blue-shifted diffuse gas detected in H$_3^+$ absorption studies. We then discuss the interaction of large-scale winds with a number of objects, such as cloudlets and stellar wind bubbles, to explain the unusual characteristics of the ISM in this region including the nonthermal radio filaments. One of the implications of this scenario is the removal of gas driven by outflowing winds may regulate star formation or black hole accretion. △ Less

Submitted 23 December, 2019; originally announced December 2019.

Comments: Four pages, one figure, in the proceedings of GCWS2019 -New Horizons in Galactic Center Astrophysics and Beyond

The hydroxyl satellite-line `flip' as a tracer of expanding HII regions

Authors: Anita Petzler, Joanne R Dawson, Mark Wardle

Abstract: Observations of the four $^{2}Π_{3/2},~J = 3/2$~ground state transitions of the hydroxyl radical (OH) have emerged as an informative tracer of molecular gas in the Galactic ISM. We discuss an OH spectral feature known as the `flip', in which the satellite lines at 1612 and 1720\,MHz flip -- one from emission to absorption and the other the reverse -- across a closely blended double feature. We hig… ▽ More Observations of the four $^{2}Π_{3/2},~J = 3/2$~ground state transitions of the hydroxyl radical (OH) have emerged as an informative tracer of molecular gas in the Galactic ISM. We discuss an OH spectral feature known as the `flip', in which the satellite lines at 1612 and 1720\,MHz flip -- one from emission to absorption and the other the reverse -- across a closely blended double feature. We highlight 30 examples of the flip from the literature, 27 of which exhibit the same orientation with respect to velocity: the 1720\,MHz line is seen in emission at more negative velocities. These same examples are also observed toward bright background continuum, many (perhaps all) show stimulated emission, and 23 of these are coincident in on-sky position and velocity with H\textsc{ii}~radio recombination lines. To explain these remarkable correlations we propose that the 1720\,MHz stimulated emission originates in heated and compressed post-shock gas expanding away from a central H\textsc{ii}~region, which collides with cooler and more diffuse gas hosting the 1612\,MHz stimulated emission. The foreground gas dominates the spectrum due to the bright central continuum, hence the expanding post-shock gas is blue-shifted relative to the stationary pre-shock gas. We employ non-LTE excitation modelling to examine this scenario, and find that indeed FIR emission from warm dust adjacent to the H\textsc{ii}~region radiatively pumps the 1612 MHz line in the diffuse, cool gas ahead of the expanding shock front, while collisional pumping in the warm, dense shocked gas inverts the 1720 MHz line. △ Less

Submitted 29 July, 2020; v1 submitted 10 November, 2019; originally announced November 2019.

Comments: Accepted to MNRAS with significant revisions to the version previously uploaded to arXiv

Cosmic-ray driven outflow from the galactic center and the origin of magnetized radio filaments

Authors: F. Yusef-Zadeh, M. Wardle

Abstract: Radio, X-ray and infrared observations of the inner few hundred pc of the Galactic center have highlighted two characteristics to the ISM. The cosmic ray ionization rate derived from molecular ions such as H$^+_3$, is at least two to three orders of magnitudes higher than in the Galactic disk. The other is bipolar X-ray and radio emission away from the Galactic plane. These features are consistent… ▽ More Radio, X-ray and infrared observations of the inner few hundred pc of the Galactic center have highlighted two characteristics to the ISM. The cosmic ray ionization rate derived from molecular ions such as H$^+_3$, is at least two to three orders of magnitudes higher than in the Galactic disk. The other is bipolar X-ray and radio emission away from the Galactic plane. These features are consistent with a scenario in which high cosmic ray pressure drives large-scale winds away from the Galactic plane. The interaction of such a wind with stellar wind bubbles may explain the energetic nonthermal radio filaments found throughout the Galactic center. Some of the implications of this scenario is the removal of gas driven by outflowing winds, acting as a feedback to reduce the star formation rate in the central molecular zone (CMZ), and the distortion of azimuthal magnetic field lines in the CMZ to vertical direction away from the plane. The combined effects of the wind and vertical magnetic field can explain why most magnetized filaments run perpendicular to the Galactic plane. This proposed picture suggests our Milky Way nucleus has recently experienced starburst or black hole activity, as recent radio and X-ray observations indicate. △ Less

Submitted 16 September, 2019; originally announced September 2019.

Comments: 6 pages, 2 figures, MNRAS, in press

Cosmic snow clouds: self-gravitating gas spheres manifesting hydrogen condensation

Authors: Mark Walker, Mark Wardle

Abstract: We present hydrostatic equilibrium models of spherical, self-gravitating clouds of helium and molecular hydrogen, focusing on the cold, high-density regime where solid- or liquid-hydrogen can form. The resulting structures have masses from 0.1 Msun down to several x 1.e-8 Msun, and span a broad range of radii: 1.e-4 < R(AU) < 1.e7. Our models are fully convective, but all have a two-zone character… ▽ More We present hydrostatic equilibrium models of spherical, self-gravitating clouds of helium and molecular hydrogen, focusing on the cold, high-density regime where solid- or liquid-hydrogen can form. The resulting structures have masses from 0.1 Msun down to several x 1.e-8 Msun, and span a broad range of radii: 1.e-4 < R(AU) < 1.e7. Our models are fully convective, but all have a two-zone character with the majority of the mass in a small, condensate-free core, surrounded by a colder envelope where phase equilibrium obtains. Convection in the envelope is unusual in that it is driven by a mean-molecular-weight inversion, rather than by an entropy gradient. In fact the entropy gradient is itself inverted, leading to the surprising result that envelope convection transports heat inwards. In turn that permits the outer layers to maintain steady state temperatures below the cosmic microwave background. Amongst our hydrostatic equilibria we identify thermal equilibria appropriate to the Galaxy, in which radiative cooling from H2 is balanced by cosmic-ray heating. These equilibria are all thermally unstable, albeit with very long thermal timescales in some cases. The specific luminosities of all our models are very low, and they therefore describe a type of baryonic dark matter. Consequently such clouds are thermally fragile: when placed in a harsh radiation field they will be unable to cool effectively and disruption will ensue as heat input drives a secular expansion. Disrupting clouds should leave trails of gas and H2 dust in their wake, which might make them easier to detect. Our models may be relevant to the cometary globules in the Helix Nebula, and the G2 cloud orbiting Sgr A*. △ Less

Submitted 13 June, 2019; originally announced June 2019.

Comments: 26 pages, 18 figures, ApJ (Accepted)

ALMA Detection of Extreme Blue-Shifted Ionized Gas Within 0.2 pc of Sgr A* from -480 to -300 km/s

Authors: Marc J. Royster, Farhad Yusef-Zadeh, Mark Wardle, Devaky Kunneriath, William Cotton, Doug A. Roberts

Abstract: We have used the capabilities of ALMA to probe the ionized gas towards the Galactic center with the H30alpha mm hydrogen recombination line within 30'' of the Galaxy's dynamical center. The observations are made with spatial and spectral resolutions of 0.46'' X 0.35'' and 3 km/s, respectively. Multiple compact and extended sources are detected in the mini-cavity region with extreme negative radial… ▽ More We have used the capabilities of ALMA to probe the ionized gas towards the Galactic center with the H30alpha mm hydrogen recombination line within 30'' of the Galaxy's dynamical center. The observations are made with spatial and spectral resolutions of 0.46'' X 0.35'' and 3 km/s, respectively. Multiple compact and extended sources are detected in the mini-cavity region with extreme negative radial velocities ranging from -480 to -300 km/s, 2-3'' (0.08-0.12 pc) from Sgr A*. This is the highest radial velocity of ionized gas detected beyond the inner 1'' of Sgr A*. This new component is kinematically isolated from the orbiting ionized gas. The detected gas has a velocity gradient ranging from -50 to -200 km/s per arcsecond located to the southwest of Sgr A* at a position angle of roughly -160 degrees. Previous proper motion measurements in the immediate vicinity of these extreme high-velocity ionized components have been measured and show transverse velocities that range from 313 - 865 km/s. If we assume that they are associated with each other, then the total velocity implies these components are gravitationally unbound. In addition, we report the kinematics of cometary radio and infrared sources. These sources are diagonally distributed with a position angle of roughly 50 degrees within 14'' of Sgr A*. We interpret the highly blue-shifted features to the SW where the mini-cavity is located and red-shifted cometary sources to the NE in terms of the interaction of a collimated outflow with an opening angle of roughly 30 degrees. An expected mass outflow rate of 2E-7 or 4E-5 solar mass per yr is estimated by a relativistic jet-driven outflow or collimated stellar winds, respectively. △ Less

Submitted 11 January, 2019; originally announced January 2019.

Comments: 22 pages, 9 figures, ApJ (in press)

ALMA Detection of Bipolar Outflows: Evidence for Low Mass Star Formation within 1pc of Sgr A*

Authors: F. Yusef-Zadeh, M. Wardle, D. Kunneriath, M. Royster, A. Wootten, D. A. Roberts

Abstract: We report the discovery of 11 bipolar outflows within a projected distance of 1pc from Sgr A* based on deep ALMA observations of $^{13}$CO, H30$α$ and SiO (5-4) lines with sub-arcsecond and $\sim1.3$ km/s, resolutions. These unambiguous signatures of young protostars manifest as approaching and receding lobes of dense gas swept up by the jets created during the formation and early evolution of sta… ▽ More We report the discovery of 11 bipolar outflows within a projected distance of 1pc from Sgr A* based on deep ALMA observations of $^{13}$CO, H30$α$ and SiO (5-4) lines with sub-arcsecond and $\sim1.3$ km/s, resolutions. These unambiguous signatures of young protostars manifest as approaching and receding lobes of dense gas swept up by the jets created during the formation and early evolution of stars. The lobe masses and momentum transfer rates are consistent with young protostellar outflows found throughout the disk of the Galaxy. The mean dynamical age of the outflow population is estimated to be $6.5^{+8.1}_{-3.6}\times10^3$ years. The rate of star formation is $\sim5\times10^{-4}$\msol\,yr$^{-1}$ assuming a mean stellar mass of $\sim0.3$ \msol. This discovery provides evidence that star formation is taking place within clouds surprisingly close to Sgr A*, perhaps due to events that compress the host cloud, creating condensations with sufficient self-gravity to resist tidal disruption by Sgr A*. Low-mass star formation over the past few billion years at this level would contribute significantly to the stellar mass budget in the central few pc of the Galaxy. The presence of many dense clumps of molecular material within 1pc of Sgr A* suggests that star formation could take place in the immediate vicinity of supermassive black holes in the nuclei of external galaxies △ Less

Submitted 28 November, 2017; originally announced November 2017.

Comments: 17 pages, 4 figures

Journal ref: ApJL 850, L30, December 1, 2017

ALMA and VLA observations of emission from the environment of Sgr A*

Authors: F. Yusef-Zadeh, R. Schödel, M. Wardle, H. Bushouse, W. Cotton, M. J. Royster, D. Kunneriath, D. A. Roberts, E. Gallego-Cano

Abstract: We present 44 and 226 GHz observations of the Galactic center within 20$"$ of Sgr A*. Millimeter continuum emission at 226 GHz is detected from eight stars that have previously been identified at near-IR and radio wavelengths. We also detect a 5.8 mJy source at 226 GHz coincident with the magnetar SGR~J1745-29 located 2.39$"$ SE of Sgr A* and identify a new 2.5$"\times1.5"$ halo of mm emission cen… ▽ More We present 44 and 226 GHz observations of the Galactic center within 20$"$ of Sgr A*. Millimeter continuum emission at 226 GHz is detected from eight stars that have previously been identified at near-IR and radio wavelengths. We also detect a 5.8 mJy source at 226 GHz coincident with the magnetar SGR~J1745-29 located 2.39$"$ SE of Sgr A* and identify a new 2.5$"\times1.5"$ halo of mm emission centered on Sgr A*. The X-ray emission from this halo has been detected previously and is interpreted in terms of a radiatively inefficient accretion flow. The mm halo surrounds an EW linear feature which appears to arise from Sgr A* and coincides with the diffuse X-ray emission and a minimum in the near-IR extinction. We argue that the millimeter emission is produced by synchrotron emission from relativistic electrons in equipartition with a $\sim 1.5$mG magnetic field. The origin of these is unclear but its coexistence with hot gas supports scenarios in which the gas is produced by the interaction of winds either from the fast moving S-stars, the photo-evaporation of low-mass YSO disks or by a jet-driven outflow from Sgr A*. The spatial anti-correlation of the X-ray, radio and mm emission from the halo and the low near-IR extinction provides compelling evidence for an outflow sweeping up the interstellar material, creating a dust cavity within 2$"$ of Sgr A*. Finally, the radio and mm counterparts to eight near-IR identified stars within $\sim$10\arcs\ of Sgr A* provide accurate astrometry to determine the positional shift between the peak emission at 44 and 226 GHz. △ Less

Submitted 8 June, 2017; originally announced June 2017.

Comments: 22 pages, 9 figures, MNRAS, in press

Tidal Distortion of the Envelope of an AGB Star IRS 3 near Sgr A*

Authors: F. Yusef-Zadeh, M. Wardle, W. Cotton, R. Schödel, M. J. Royster, D. A. Roberts, D. Kunneriath

Abstract: We present radio and mm continuum observations of the Galactic center taken with the VLA and ALMA at 44 and 226 GHz, respectively. We detect radio and mm emission from IRS 3, lying ~4.5" NW of Sgr A*, with a spectrum that is consistent with the photospheric emission from an AGB star at the Galactic center. Millimeter images reveal that the envelope of IRS 3, the brightest and most extended 3.8$μ$m… ▽ More We present radio and mm continuum observations of the Galactic center taken with the VLA and ALMA at 44 and 226 GHz, respectively. We detect radio and mm emission from IRS 3, lying ~4.5" NW of Sgr A*, with a spectrum that is consistent with the photospheric emission from an AGB star at the Galactic center. Millimeter images reveal that the envelope of IRS 3, the brightest and most extended 3.8$μ$m Galactic center stellar source, consists of two semi-circular dust shells facing the direction of Sgr A*. The outer circumstellar shell at the distance of 1.6$\times10^4$ AU, appears to break up into "fingers" of dust directed toward Sgr A*. These features coincide with molecular CS (5-4) emission and a near-IR extinction cloud distributed between IRS 3 and Sgr A*. The NE-SW asymmetric shape of the IRS 3 shells seen at 3.8 micron and radio are interpreted as structures that are tidally distorted by Sgr A*. Using the kinematics of CS emission and the proper motion of IRS 3, the tidally distorted outflowing material from the envelope after 5000 years constrains the distance of IRS 3 to $\sim$0.7 pc in front of or $\sim$0.5 pc behind Sgr A*. This suggests that the mass loss by stars near Sgr A* can supply a reservoir of molecular material near Sgr A*. We also present dark features in radio continuum images coincident with the envelope of IRS 3. These dusty stars provide examples in which high resolution radio continuum images can identify dust enshrouded stellar sources embedded an ionized medium. △ Less

Submitted 7 February, 2017; originally announced February 2017.

Comments: 16 pages, 5 figures, ApJ (in press)

ALMA and VLA Observations: Evidence for Ongoing Low-mass Star Formation near Sgr A*

Authors: F. Yusef-Zadeh, W. Cotton, M. Wardle, M. J. Royster, D. Kunneriath, D. A. Roberts, A. Wootten, R. Schödel

Abstract: Using the VLA, we recently detected a large number of protoplanetary disk (proplyd) candidates lying within a couple of light years of the massive black hole Sgr A*. The bow-shock appearance of proplyd candidates point toward the young massive stars located near Sgr A*. Similar to Orion proplyds, the strong UV radiation from the cluster of massive stars at the Galactic center is expected to photoe… ▽ More Using the VLA, we recently detected a large number of protoplanetary disk (proplyd) candidates lying within a couple of light years of the massive black hole Sgr A*. The bow-shock appearance of proplyd candidates point toward the young massive stars located near Sgr A*. Similar to Orion proplyds, the strong UV radiation from the cluster of massive stars at the Galactic center is expected to photoevaporate and photoionize the circumstellar disks around young, low mass stars, thus allowing detection of the ionized outflows from the photoionized layer surrounding cool and dense gaseous disks. To confirm this picture, ALMA observations detect millimeter emission at 226 GHz from five proplyd candidates that had been detected at 44 and 34 GHz with the VLA. We present the derived disk masses for four sources as a function of the assumed dust temperature. The mass of protoplanetary disks from cool dust emission ranges between 0.03 -- 0.05 solar mass. These estimates are consistent with the disk masses found in star forming sites in the Galaxy. These measurements show the presence of on-going star formation with the implication that gas clouds can survive near Sgr A* and the relative importance of high vs low-mass star formation in the strong tidal and radiation fields of the Galactic center. △ Less

Submitted 20 January, 2017; originally announced January 2017.

Comments: 13 pages, 3 figures, MNRAS (in press)

Two-fluid dusty shocks: simple benchmarking problems and applications to protoplanetary discs

Authors: Andrew Lehmann, Mark Wardle

Abstract: The key role that dust plays in the interstellar medium has motivated the development of numerical codes designed to study the coupled evolution of dust and gas in systems such as turbulent molecular clouds and protoplanetary discs. Drift between dust and gas has proven to be important as well as numerically challenging. We provide simple benchmarking problems for dusty gas codes by numerically so… ▽ More The key role that dust plays in the interstellar medium has motivated the development of numerical codes designed to study the coupled evolution of dust and gas in systems such as turbulent molecular clouds and protoplanetary discs. Drift between dust and gas has proven to be important as well as numerically challenging. We provide simple benchmarking problems for dusty gas codes by numerically solving the two-fluid dust-gas equations for steady, plane-parallel shock waves. The two distinct shock solutions to these equations allow a numerical code to test different forms of drag between the two fluids, the strength of that drag and the dust to gas ratio. We also provide an astrophysical application of J-type dust-gas shocks to studying the structure of accretion shocks onto protoplanetary discs. We find that two-fluid effects are most important for grains larger than 1 um, and that the peak dust temperature within an accretion shock provides a signature of the dust-to-gas ratio of the infalling material. △ Less

Submitted 29 December, 2016; originally announced December 2016.

Comments: 9 Pages, 8 Figures, submitted to MNRAS, comments welcome

Star Formation Close to Sgr A* and Beyond the Nuclear Cluster

Authors: F. Yusef-Zadeh, M. Wardle

Abstract: Two modes of star formation are involved to explain the origin of young stars near Sgr A*. One is a disk-based mode, which explains the disk of stars orbiting Sgr A*. The other is the standard cloud-based mode observed in the Galactic disk. We discuss each of these modes of star formation and apply these ideas to the inner few parsecs of Sgr A*. In particular, we focus on the latter mode in more d… ▽ More Two modes of star formation are involved to explain the origin of young stars near Sgr A*. One is a disk-based mode, which explains the disk of stars orbiting Sgr A*. The other is the standard cloud-based mode observed in the Galactic disk. We discuss each of these modes of star formation and apply these ideas to the inner few parsecs of Sgr A*. In particular, we focus on the latter mode in more detail. We also discuss how the tidal force exerted by the nuclear cluster makes the Roche density approaching zero and contributes to the collapse of molecular clouds located tens of parsecs away from Sgr A*. △ Less

Submitted 28 November, 2016; originally announced November 2016.

Comments: 4 pages, 2 figures, International Astronomical Union Proceedings Series 324, editors: Andreja Gomboc and Carole Mundell

SHOCKFIND - An algorithm to identify magnetohydrodynamic shock waves in turbulent clouds

Authors: Andrew Lehmann, Christoph Federrath, Mark Wardle

Abstract: The formation of stars occurs in the dense molecular cloud phase of the interstellar medium. Observations and numerical simulations of molecular clouds have shown that supersonic magnetised turbulence plays a key role for the formation of stars. Simulations have also shown that a large fraction of the turbulent energy dissipates in shock waves. The three families of MHD shocks --- fast, intermedia… ▽ More The formation of stars occurs in the dense molecular cloud phase of the interstellar medium. Observations and numerical simulations of molecular clouds have shown that supersonic magnetised turbulence plays a key role for the formation of stars. Simulations have also shown that a large fraction of the turbulent energy dissipates in shock waves. The three families of MHD shocks --- fast, intermediate and slow --- distinctly compress and heat up the molecular gas, and so provide an important probe of the physical conditions within a turbulent cloud. Here we introduce the publicly available algorithm, SHOCKFIND, to extract and characterise the mixture of shock families in MHD turbulence. The algorithm is applied to a 3-dimensional simulation of a magnetised turbulent molecular cloud, and we find that both fast and slow MHD shocks are present in the simulation. We give the first prediction of the mixture of turbulence-driven MHD shock families in this molecular cloud, and present their distinct distributions of sonic and Alfvenic Mach numbers. Using subgrid one-dimensional models of MHD shocks we estimate that ~0.03 % of the volume of a typical molecular cloud in the Milky Way will be shock heated above 50 K, at any time during the lifetime of the cloud. We discuss the impact of this shock heating on the dynamical evolution of molecular clouds. △ Less

Submitted 5 August, 2016; originally announced August 2016.

Comments: 15 pages, 16 figures, accepted by MNRAS

An OH(1720 MHz) Maser and a Nonthermal Radio Source in Sgr B2(M): A SNR--Molecular Cloud Interaction Site?

Authors: F. Yusef-Zadeh, W. Cotton, M. Wardle, H. Intema

Abstract: Sgr B2 is a well-known star forming molecular cloud complex in the Galactic center region showing evidence of high energy activity as traced by the K$α$ neutral FeI line at 6.4 keV, as well as GeV and TeV $γ$-ray emission. Here we present VLA and GMRT observations with respective resolutions of $\approx3.5"\times1.2"$ and 25$"\times25"$ and report the detection of an OH(1720 MHz) maser, with no ac… ▽ More Sgr B2 is a well-known star forming molecular cloud complex in the Galactic center region showing evidence of high energy activity as traced by the K$α$ neutral FeI line at 6.4 keV, as well as GeV and TeV $γ$-ray emission. Here we present VLA and GMRT observations with respective resolutions of $\approx3.5"\times1.2"$ and 25$"\times25"$ and report the detection of an OH(1720 MHz) maser, with no accompanying OH 1665, 1667 and 1612 MHz maser emission. The maser coincides with a 150 MHz nonthermal radio source in Sgr B2(M). This rare class of OH(1720 MHz) masers or the so-called supernova remnant (SNR) masers, with no main line transitions, trace shocked gas and signal the interaction of an expanding SNR with a molecular cloud. We interpret the 150 MHz radio source as either the site of a SNR -- molecular gas interaction or a wind-wind collision in a massive binary system. The interaction of the molecular cloud and the nonthermal source enhances the cosmic-ray ionization rate, allows the diffusion of cosmic rays into the cloud and produces the variable 6.4 keV line, GeV and TeV $γ$-ray emission from Sgr B2(M). The cosmic ray electron interaction with the gas in the Galactic center can not only explain the measured high values of cosmic ray ionization and heating rates but also contribute to nonthermal bremsstrahlung continuum emission, all of which are consistent with observations. △ Less

Submitted 2 March, 2016; originally announced March 2016.

Comments: 15 pages, 4 figures, ApJL (in press)

Sgr A* and its Environment: Low Mass Star Formation, the Origin of X-ray Gas and Collimated Outflow

Authors: F. Yusef-Zadeh, M. Wardle, R. Schödel, D. A. Roberts, W. Cotton, H. Bushouse, R. Arendt, M. Royster

Abstract: We present high-resolution multiwavelength radio continuum images of the region within 150$"$ of Sgr A*, revealing a number of new extended features and stellar sources in this region. First, we detect a continuous 2" east-west ridge of radio emission, linking Sgr A* and a cluster of stars associated with IRS 13N and IRS 13E. The ridge suggests that an outflow of east-west blob-like structures is… ▽ More We present high-resolution multiwavelength radio continuum images of the region within 150$"$ of Sgr A*, revealing a number of new extended features and stellar sources in this region. First, we detect a continuous 2" east-west ridge of radio emission, linking Sgr A* and a cluster of stars associated with IRS 13N and IRS 13E. The ridge suggests that an outflow of east-west blob-like structures is emerging from Sgr A*. We also find arc-like features within the ridge with morphologies suggestive of photoevaporative protoplanetary disks. We use near-IR fluxes to show that the emission has similar characteristics to those of a protoplanetary disk irradiated by the intense radiation field at the Galactic center. This suggests that star formation has taken place within the S cluster 2$"$ from Sgr A*. We suggest that the diffuse X-ray emission associated with Sgr A* is due to an expanding hot wind produced by the mass loss from B-type main sequence stars, and/or the disks of photoevaporation of low mass YSOs at a rate ~10^{-6} solar mass per year. The proposed model naturally reduces the inferred accretion rate and is an alternative to the inflow-outflow style models to explain the underluminous nature of Sgr A*. Second, we detect new cometary radio and nea-IR sources and a striking tower of radio emission suggesting that they are tracing interaction sites of a mildly relativistic jet from Sgr A* with the atmosphere of stars and the nonthermal Sgr A East shell at a PA$\sim50-60^\circ$ with ~10^{-7} solar mass per year, and opening angle 10 degrees. △ Less

Submitted 1 January, 2016; originally announced January 2016.

Comments: 45 pages, 12 figures, ApJ (in press)

Radio Continuum Emission from the Magnetar SGR J1745-2900: Interaction with Gas Orbiting Sgr A*

Authors: F. Yusef-Zadeh, R. Diesing, M. Wardle, L. O. Sjouwerman, M. Royster, W. D. Cotton, D. Roberts, C. Heinke

Abstract: We present radio continuum light curves of the magnetar SGR J1745$-$2900 and Sgr A* obtained with multi-frequency, multi-epoch Very Large Array observations between 2012 and 2014. During this period, a powerful X-ray outburst from SGR J1745$-$2900 occurred on 2013-04-24. Enhanced radio emission is delayed with respect to the X-ray peak by about seven months. In addition, the flux density of the em… ▽ More We present radio continuum light curves of the magnetar SGR J1745$-$2900 and Sgr A* obtained with multi-frequency, multi-epoch Very Large Array observations between 2012 and 2014. During this period, a powerful X-ray outburst from SGR J1745$-$2900 occurred on 2013-04-24. Enhanced radio emission is delayed with respect to the X-ray peak by about seven months. In addition, the flux density of the emission from the magnetar fluctuates by a factor of 2 to 4 at frequencies between 21 and 41 GHz and its spectral index varies erratically. Here we argue that the excess fluctuating emission from the magnetar arises from the interaction of a shock generated from the X-ray outburst with the orbiting ionized gas at the Galactic center. In this picture, variable synchrotron emission is produced by ram pressure variations due to inhomogeneities in the dense ionized medium of the Sgr A West bar. The pulsar with its high transverse velocity is moving through a highly blue-shifted ionized medium. This implies that the magnetar is at a projected distance of $\sim0.1$ pc from Sgr A* and that the orbiting ionized gas is partially or largely responsible for a large rotation measure detected toward the magnetar. Despite the variability of Sgr A* expected to be induced by the passage of the G2 cloud, monitoring data shows a constant flux density and spectral index during this period △ Less

Submitted 10 September, 2015; originally announced September 2015.

Comments: 12 pages, 3 figures, ApJL (in press)

Signatures of fast and slow magnetohydrodynamic shocks in turbulent molecular clouds

Authors: Andrew Lehmann, Mark Wardle

Abstract: The character of star formation is intimately related to the supersonic magnetohydrodynamic (MHD) turbulent dynamics of the molecular clouds in which stars form. A significant amount of the turbulent energy dissipates in low velocity shocks. Fast and slow MHD shocks differ in how they compress and heat the molecular gas, and so their radiative signatures reveal distinct physical conditions. We u… ▽ More The character of star formation is intimately related to the supersonic magnetohydrodynamic (MHD) turbulent dynamics of the molecular clouds in which stars form. A significant amount of the turbulent energy dissipates in low velocity shocks. Fast and slow MHD shocks differ in how they compress and heat the molecular gas, and so their radiative signatures reveal distinct physical conditions. We use a two-fluid model to compare one-dimensional fast and slow MHD shocks propagating at low speeds (a few km/s). Fast shocks are magnetically driven, forcing ion species to stream through the neutral gas ahead of the shock front. This magnetic precursor heats the gas sufficiently to create a large, warm transition zone where all the fluid variables smoothly change in the shock front. In contrast, slow shocks are driven by gas pressure, and neutral species collide with ion species in a thin hot slab that closely resembles an ordinary gas dynamic shock. We consider shocks at velocities $v_s = 2$-$4$ km/s and preshock Hydrogen nuclei densities $n_\mathrm{H} = 10^2$-$10^4$ cm$^{-3}$. We include a simple oxygen chemistry and cooling by CO, H$_2$ and H$_2$O. CO rotational lines above $J = 6\to 5$ are more strongly excited in slow shocks. These slow shock signatures may have already been observed in infrared dark clouds in the Milky Way. △ Less

Submitted 8 October, 2015; v1 submitted 8 July, 2015; originally announced July 2015.

Comments: 12 pages, 8 figures, Accepted to MNRAS

Compact Radio Sources within 30" of Sgr A*: Proper Motions, Stellar Winds and the Accretion Rate onto Sgr A*

Authors: F. Yusef-Zadeh, H. Bushouse, R. Schödel, M. Wardle, W. Cotton, D. A. Roberts, F. Nogueras-Lara, E. Gallego-Cano

Abstract: Recent broad-band 34 and 44 GHz radio continuum observations of the Galactic center have revealed 41 massive stars identified with near-IR counterparts, as well as 44 proplyd candidates within 30" of Sgr A*. Radio observations obtained in 2011 and 2014 have been used to derive proper motions of eight young stars near Sgr A*. The accuracy of proper motion estimates based on near-IR observations by… ▽ More Recent broad-band 34 and 44 GHz radio continuum observations of the Galactic center have revealed 41 massive stars identified with near-IR counterparts, as well as 44 proplyd candidates within 30" of Sgr A*. Radio observations obtained in 2011 and 2014 have been used to derive proper motions of eight young stars near Sgr A*. The accuracy of proper motion estimates based on near-IR observations by Lu et al. and Paumard et al. have been investigated by using their proper motions to predict the 2014 epoch positions of near-IR stars and comparing the predicted positions with those of radio counterparts in the 2014 radio observations. Predicted positions from Lu et al. show an rms scatter of 6 mas relative to the radio positions, while those from Paumard et al. show rms residuals of 20 mas, which is mainly due to uncertainties in the IR-based proper motions. Under the assumption of homogeneous ionized winds, we also determine the mass-loss rates of 11 radio stars, finding rates that are on average $\sim$2 times smaller than those determined from model atmosphere calculations and near-IR data. Clumpiness of ionized winds would reduce the mass loss rate of WR and O stars by additional factors of 3 and 10, respectively. One important implication of this is a reduction in the expected mass accretion rate onto Sgr A* from stellar winds by nearly an order of magnitude to a value of few$\times10^{-7}$ \msol\ yr$^{-1}$. Finally, we present the positions of 318 compact 34.5 GHz radio sources within 30\arcs\ of Sgr A*. At least 45 of these have stellar counterparts in the near-IR $K_s$ (2.18 $μ$m) and $L'$ (3.8$μ$m) bands. △ Less

Submitted 23 June, 2015; originally announced June 2015.

Comments: 30 pages, 4 figures, ApJ (in press)

Signatures of Young Star Formation Activity Within Two Parsecs of Sgr A*

Authors: F. Yusef-Zadeh, M. Wardle, M. Sewilo, D. A. Roberts, I. Smith, R. Arendt, W. Cotton, J. Lacy, S. Martin, M. W. Pound, M. Rickett, M. Royster

Abstract: We present radio and infrared observations indicating on-going star formation activity inside the $\sim2-5$ pc circumnuclear ring at the Galactic center. Collectively these measurements suggest a continued disk-based mode of on-going star formation has taken place near Sgr A* over the last few million years. First, VLA observations with spatial resolution 2.17$"\times0.81"$ reveal 13 water masers,… ▽ More We present radio and infrared observations indicating on-going star formation activity inside the $\sim2-5$ pc circumnuclear ring at the Galactic center. Collectively these measurements suggest a continued disk-based mode of on-going star formation has taken place near Sgr A* over the last few million years. First, VLA observations with spatial resolution 2.17$"\times0.81"$ reveal 13 water masers, several of which have multiple velocity components. The presence of interstellar water masers suggests gas densities that are sufficient for self-gravity to overcome the tidal shear of the 4$\times10^6$ \msol\, black hole. Second, SED modeling of stellar sources indicate massive YSO candidates interior to the molecular ring, supporting in-situ star formation near Sgr A* and appear to show a distribution similar to that of the counter-rotating disks of $\sim$100 OB stars orbiting Sgr A*. Some YSO candidates (e.g., IRS~5) have bow shock structures suggesting that they have have gaseous disks that are phototoevaporated and photoionized by the strong radiation field. Third, we detect clumps of SiO (2-1) and (5-4) line emission in the ring based on CARMA and SMA observations. The FWHM and luminosity of the SiO emission is consistent with shocked protostellar outflows. Fourth, two linear ionized features with an extent of $\sim0.8$ pc show blue and redshifted velocities between $+50$ and $-40$ \kms, suggesting protostellar jet driven outflows with mass loss rates of $\sim5\times10^{-5}$ solar mass yr$^{-1}$. Finally, we present the imprint of radio dark clouds at 44 GHz, representing a reservoir of molecular gas that feeds star formation activity close to Sgr A*. △ Less

Submitted 19 May, 2015; originally announced May 2015.

Comments: 38 pages, 10 figures, ApJ (in press)

Magnetic fields in gaps surrounding giant protoplanets

Authors: Sarah L. Keith, Mark Wardle

Abstract: Giant protoplanets evacuate a gap in their host protoplanetary disc, which gas must cross before it can be accreted. A magnetic field is likely carried into the gap, potentially influencing the flow. Gap crossing has been simulated with varying degrees of attention to field evolution (pure hydrodynamical, ideal, and resistive MHD), but as yet there has been no detailed assessment of the role of th… ▽ More Giant protoplanets evacuate a gap in their host protoplanetary disc, which gas must cross before it can be accreted. A magnetic field is likely carried into the gap, potentially influencing the flow. Gap crossing has been simulated with varying degrees of attention to field evolution (pure hydrodynamical, ideal, and resistive MHD), but as yet there has been no detailed assessment of the role of the field accounting for all three key non-ideal MHD effects: Ohmic resistivity, ambipolar diffusion, and Hall drift. We present a detailed investigation of gap magnetic field structure as determined by non-ideal effects. We assess susceptibility to turbulence induced by the magnetorotational instability, and angular momentum loss from large-scale fields. As full non-ideal simulations are computationally expensive, we take an a posteriori approach, estimating MHD quantities from the pure hydrodynamical gap crossing simulation by Tanigawa et al. (2012). We calculate the ionisation fraction and estimate field strength and geometry to determine the strength of non-ideal effects. We find that the protoplanetary disc field would be easily drawn into the gap and circumplanetary disc. Hall drift dominates, so that much of the gap is conditionally MRI unstable depending on the alignment of the field and disc rotation axes. Field alignment also influences the strong toroidal field component permeating the gap. Large-scale magnetic forces are small in the circumplanetary disc, indicating they cannot drive accretion there. However, turbulence will be key during satellite growth as it affects critical disc features, such as the location of the ice line. △ Less

Submitted 8 May, 2015; originally announced May 2015.

Comments: 14 pages, 6 figures, accepted for publication by MNRAS

Radio Continuum Observations of the Galactic Center: Photoevaporative Proplyd-like Objects near Sgr A*

Authors: F. Yusef-Zadeh, D. A. Roberts, M. Wardle, W. Cotton, R. Schödel, M. J. Royster

Abstract: We present radio images within 30$''$ of Sgr A* based on recent VLA observations at 34 GHz with 7.8 microJy sensitivity and resolution $\sim88\times46$ milliarcseconds (mas). We report 44 partially resolved compact sources clustered in two regions in the E arm of ionized gas that orbits Sgr A*. These sources have size scales ranging between ~50 and 200 mas (400 to 1600 AUs), and a bow-shock appear… ▽ More We present radio images within 30$''$ of Sgr A* based on recent VLA observations at 34 GHz with 7.8 microJy sensitivity and resolution $\sim88\times46$ milliarcseconds (mas). We report 44 partially resolved compact sources clustered in two regions in the E arm of ionized gas that orbits Sgr A*. These sources have size scales ranging between ~50 and 200 mas (400 to 1600 AUs), and a bow-shock appearance facing the direction of Sgr A*. Unlike the bow-shock sources previously identified in the near-IR but associated with massive stars, these 34 GHz sources do not appear to have near-IR counterparts at 3.8 $μ$m. We interpret these sources as a candidate population of photoevaporative protoplanetary disks (proplyds) that are associated with newly formed low mass stars with mass loss rates ~10^{-7} - 10^{-6} solar mass per year and are located at the edge of a molecular cloud outlined by ionized gas. The disks are externally illuminated by strong Lyman continuum radiation from the ~100 OB and WR massive stars distributed within 10'' of Sgr A*. The presence of proplyds implies current in-situ star formation activity near Sgr A* and opens a window for the first time to study low mass star, planetary and brown dwarf formations near a supermassive black hole. △ Less

Submitted 10 February, 2015; originally announced February 2015.

Comments: 13 pages, 4 figures, ApJL (in press)

The Discovery of Radio Stars within 10 arcseconds of Sgr A* at 7mm

Authors: F. Yusef-Zadeh, D. A. Roberts, H. Bushouse, M. Wardle, W. Cotton, M. Royster, G. van Moorsel

Abstract: Very Large Array observations of the Galactic Center at 7 mm have produced an image of the 30 arcseconds surrounding Sgr A* with a resolution of 82x42 milliarcseconds (mas). A comparison with IR images taken simultaneously with the Very Large Telescope (VLT) identifies 41 radio sources with L-band (3.8 microns) stellar counterparts. The well-known young, massive stars in the central Sgr A* cluster… ▽ More Very Large Array observations of the Galactic Center at 7 mm have produced an image of the 30 arcseconds surrounding Sgr A* with a resolution of 82x42 milliarcseconds (mas). A comparison with IR images taken simultaneously with the Very Large Telescope (VLT) identifies 41 radio sources with L-band (3.8 microns) stellar counterparts. The well-known young, massive stars in the central Sgr A* cluster (e.g., IRS 16C, IRS 16NE, IRS 16SE2, IRS 16NW, IRS 16SW, AF, AFNW, IRS 34W and IRS 33E) are detected with peak flux densities between 0.2 and 1.3 mJy. The origin of the stellar radio emission in the central cluster is discussed in terms of ionized stellar winds with mass-loss rates in the range 0.8-5x10^{-5} solar mass per year. Radio emission from eight massive stars is used as a tool for registration between the radio and infrared frames with mas precision within a few arcseconds of Sgr A*. This is similar to the established technique of aligning SiO masers and evolved stars except that radio stars lie within a few arcseconds of Sgr A*. Our data show a scatter of ~6.5 mas in the positions of the eight radio sources that appear in both the L-band and 7 mm images. Lastly, we use the radio and IR data to argue that members of IRS 13N are Young Stellar Objects rather than dust clumps, supporting the hypothesis that recent star formation has occurred near Sgr A*. △ Less

Submitted 15 July, 2014; originally announced July 2014.

Comments: 15 pages, 3 figures, ApJL (in press)

Accretion in giant planet circumplanetary disks

Authors: Sarah L. Keith, Mark Wardle

Abstract: During the final growth phase of giant planets, accretion is thought to be controlled by a surrounding circumplanetary disk. Current astrophysical accretion disk models rely on hydromagnetic turbulence or gravitoturbulence as the source of effective viscosity within the disk. However, the magnetically-coupled accreting region in these models is so limited that the disk may not support inflow at al… ▽ More During the final growth phase of giant planets, accretion is thought to be controlled by a surrounding circumplanetary disk. Current astrophysical accretion disk models rely on hydromagnetic turbulence or gravitoturbulence as the source of effective viscosity within the disk. However, the magnetically-coupled accreting region in these models is so limited that the disk may not support inflow at all radii, or at the required rate. Here, we examine the conditions needed for self-consistent accretion, in which the disk is susceptible to accretion driven by magnetic fields or gravitational instability. We model the disk as a Shakura-Sunyaev $α$ disk and calculate the level of ionisation, the strength of coupling between the field and disk using Ohmic, Hall and Ambipolar diffusevities for both an MRI and vertical field, and the strength of gravitational instability. We find that the standard constant-$α$ disk is only coupled to the field by thermal ionisation within $30\,R_J$ with strong magnetic diffusivity prohibiting accretion through the bulk of the midplane. In light of the failure of the constant-$α$ disk to produce accretion consistent with its viscosity we drop the assumption of constant-$α$ and present an alternate model in which $α$ varies radially according to the level magnetic turbulence or gravitoturbulence. We find that a vertical field may drive accretion across the entire disk, whereas MRI can drive accretion out to $\sim200\,R_J$, beyond which Toomre's $Q=1$ and gravitoturbulence dominates. The disks are relatively hot ($T\gtrsim800\,$K), and consequently massive ($M_{\text{disk}}\sim0.5\,M_J$). △ Less

Submitted 5 February, 2014; originally announced February 2014.

Comments: 19 pages, 5 figures, 2 tables; accepted for publication by MNRAS

On the origin of the central 1" hole in the stellar disk of Sgr A* and the Fermi gamma-ray bubbles

Authors: Mark Wardle, Farhad Yusef-Zadeh

Abstract: The supermassive black hole Sgr A* at the center of the Galaxy is surrounded by two misaligned disks of young, massive stars extending from ~0.04 to 0.4 pc. The stellar surface density increases as ~ r^-2 towards Sgr A* but is truncated within 1" (0.04pc). We explore the origin of this annulus using a model in which star formation occurs in a disk of gas created through the partial capture of a ga… ▽ More The supermassive black hole Sgr A* at the center of the Galaxy is surrounded by two misaligned disks of young, massive stars extending from ~0.04 to 0.4 pc. The stellar surface density increases as ~ r^-2 towards Sgr A* but is truncated within 1" (0.04pc). We explore the origin of this annulus using a model in which star formation occurs in a disk of gas created through the partial capture of a gas cloud as it sweeps through the inner few parsecs of the galaxy and temporarily engulfs Sgr A*. We identify the locations within which star formation and/or accretion onto Sgr A* take place. Within 0.04 pc the disk is magnetically active and the associated heating and enhanced pressure prevents the disk from becoming self gravitating. Instead, it forms a magneto-turbulent disk that drains onto Sgr A* within 3 Myr. Meanwhile, fragmentation of the gas beyond the central 0.04 pc hole creates the observed young stellar disk. The two large scale bubbles of gamma-ray emission extending perpendicular to the Galactic plane may be created by a burst of accretion of ~10^5 Msun of gas lying between 0.01 and 0.03 pc. The observed stellar ages imply that this capture event occurred ~10^6.5 yr ago, thus such events occurring over the life time of the Galaxy could have significantly contributed to the current mass of Sgr A* and to the inner few parsec of the nuclear star cluster. We suggest that these events also occur in extragalactic systems. △ Less

Submitted 17 April, 2014; v1 submitted 31 January, 2014; originally announced January 2014.

Comments: ApJL accepted; eq 1 for surface density profile generalised from previous version plus minor changes

Transport and Accretion in Planet-Forming Disks

Authors: N. J. Turner, S. Fromang, C. Gammie, H. Klahr, G. Lesur, M. Wardle, X. -N. Bai

Abstract: Planets appear to form in environments shaped by the gas flowing through protostellar disks to the central young stars. The flows in turn are governed by orbital angular momentum transfer. In this chapter we summarize current understanding of the transfer processes best able to account for the flows, including magneto-rotational turbulence, magnetically-launched winds, self-gravitational instabili… ▽ More Planets appear to form in environments shaped by the gas flowing through protostellar disks to the central young stars. The flows in turn are governed by orbital angular momentum transfer. In this chapter we summarize current understanding of the transfer processes best able to account for the flows, including magneto-rotational turbulence, magnetically-launched winds, self-gravitational instability and vortices driven by hydrodynamical instabilities. For each in turn we outline the major achievements of the past few years and the outstanding questions. We underscore the requirements for operation, especially ionization for the magnetic processes and heating and cooling for the others. We describe the distribution and strength of the resulting flows and compare with the long-used phenomenological $α$-picture, highlighting issues where the fuller physical picture yields substantially different answers. We also discuss the links between magnetized turbulence and magnetically-launched outflows, and between magnetized turbulence and hydrodynamical vortices. We end with a summary of the status of efforts to detect specific signatures of the flows. △ Less

Submitted 28 January, 2014; originally announced January 2014.

Comments: Refereed review chapter accepted for publication in Protostars and Planets VI, University of Arizona Press (2014), eds. H. Beuther, R. Klessen, C. Dullemond & Th. Henning. 24 pages, 7 figures. The talk presented at the Protostars and Planets VI conference is available at http://youtube.com/watch?v=tEgw0PXwkGE

An analysis of HCN observations of the Circumnuclear Disk at the galactic centre

Authors: Ian L Smith, Mark Wardle

Abstract: The Circumnuclear Disk (CND) is a torus of dust and molecular gas rotating about the galactic centre and extending from approximately 1.6pc to 7pc from the central massive black hole, SgrA*. Large Velocity Gradient modelling of the intensities of the HCN 1-0, 3-2 and 4-3 transitions is used to infer hydrogen density and HCN optical depth. From HCN observations we find the molecular hydrogen densit… ▽ More The Circumnuclear Disk (CND) is a torus of dust and molecular gas rotating about the galactic centre and extending from approximately 1.6pc to 7pc from the central massive black hole, SgrA*. Large Velocity Gradient modelling of the intensities of the HCN 1-0, 3-2 and 4-3 transitions is used to infer hydrogen density and HCN optical depth. From HCN observations we find the molecular hydrogen density ranges from 0.1 to 2 $\times$ 10$^{6}$ cm$^{-3}$, about an order of magnitude less than inferred previously. The 1-0 line is weakly inverted with line-centre optical depth approx $-$0.1, in stark contrast to earlier estimates of 4. The estimated mass of the ring is approximately 3 $-$ 4 $\times$ 10$^{5}$M$_{\odot}$ consistent with estimates based on thermal dust emission. The tidal shear in the disk implies that star formation is not expected to occur without some significant triggering event. △ Less

Submitted 31 October, 2013; originally announced October 2013.

Comments: 17 pages, 8 figures, 7 tables (including 3 landscape tables generated with 2 separate Latex files)

Constraints on Common Envelope Magnetic Fields from Observations of Jets in Planetary Nebulae

Authors: James Tocknell, Orsola De Marco, Mark Wardle

Abstract: The common envelope (CE) interaction describes the swallowing of a nearby companion by a growing, evolving star. CEs that take place during the asymptotic giant branch phase of the primary and may lead to the formation of a planetary nebula (PN) with a post-CE close binary in the middle. We have used published observations of masses and kinematics of jets in four post-CE PN to infer physical chara… ▽ More The common envelope (CE) interaction describes the swallowing of a nearby companion by a growing, evolving star. CEs that take place during the asymptotic giant branch phase of the primary and may lead to the formation of a planetary nebula (PN) with a post-CE close binary in the middle. We have used published observations of masses and kinematics of jets in four post-CE PN to infer physical characteristics of the CE interaction. In three of the four systems studied, Abell 63, ETHOS 1 and the Necklace PN, the kinematics indicate that the jets were launched a few thousand years before the CE and we favour a scenario where this happened before Roche lobe overflow, although better models of wind accretion and wind Roche lobe overflow are needed. The magnetic fields inferred to launch pre-CE jets are of the order of a few Gauss. In the fourth case, NGC 6778, the kinematics indicate that the jets were launched about 3000 years after the CE interaction. Magnetic fields of the order of a few hundreds to a few thousands Gauss are inferred in this case, approximately in line with predictions of post-CE magnetic fields. However, we remark that in the case of this system, it is impossible to find a reasonable scenario for the formation of the two jet pairs observed: the small orbital separation would preclude the formation of even one accretion disk able to supply the necessary accretion rate to cause the observed jets. △ Less

Submitted 13 January, 2014; v1 submitted 22 August, 2013; originally announced August 2013.

Comments: 9 pages, accepted by MNRAS, Key words: Magnetic Fields, ISM: Jets and Outflows, Planetary Nebulae: Individual: Necklace, Planetary Nebulae: Individual: Abell 63, Planetary Nebulae: Individual: ETHOS 1, Planetary Nebulae: Individual: NGC 6778

74 MHz Nonthermal Emission from Molecular Clouds: Evidence for a Cosmic Ray Dominated Region at the Galactic Center

Authors: F. Yusef-Zadeh, M. Wardle, D. Lis, S. Viti, C. Brogan, E. Chambers, M. Pound, M. Rickert

Abstract: We present 74 MHz radio continuum observations of the Galactic center region. These measurements show nonthermal radio emission arising from molecular clouds that is unaffected by free-free absorption along the line of sight. We focus on one cloud, G0.13--0.13, representative of the population of molecular clouds that are spatially correlated with steep spectrum (alpha^{74MHz}_{327MHz}=1.3\pm0.3)… ▽ More We present 74 MHz radio continuum observations of the Galactic center region. These measurements show nonthermal radio emission arising from molecular clouds that is unaffected by free-free absorption along the line of sight. We focus on one cloud, G0.13--0.13, representative of the population of molecular clouds that are spatially correlated with steep spectrum (alpha^{74MHz}_{327MHz}=1.3\pm0.3) nonthermal emission from the Galactic center region. This cloud lies adjacent to the nonthermal radio filaments of the Arc near l~0.2^0 and is a strong source of 74 MHz continuum, SiO (2-1) and FeI Kalpha 6.4 keV line emission. This three-way correlation provides the most compelling evidence yet that relativistic electrons, here traced by 74 MHz emission, are physically associated with the G0.13--0.13 molecular cloud and that low energy cosmic ray electrons are responsible for the FeI Kalpha line emission. The high cosmic ray ionization rate ~10-13 s-1 H-1 is responsible for heating the molecular gas to high temperatures and allows the disturbed gas to maintain a high velocity dispersion. LVG modeling of multi-transition SiO observations of this cloud implies H2 densities ~104-5 cm-3 and high temperatures. The lower limit to the temperature of G0.13-0.13 is ~100K, whereas the upper limit is as high as 1000K. Lastly, we used a time-dependent chemical model in which cosmic rays drive the chemistry of the gas to investigate for molecular line diagnostics of cosmic ray heating. When the cloud reaches chemical equilibrium, the abundance ratios of HCN/HNC and N2H+/HCO+ are consistent with measured values. In addition, significant abundance of SiO is predicted in the cosmic ray dominated region of the Galactic center. We discuss different possibilities to account for the origin of widespread SiO emission detected from Galactic center molecular clouds. △ Less

Submitted 5 May, 2013; originally announced May 2013.

Comments: 48 pages, 10 figures, The Journal of Physical Chemistry, part A, Oka Festschrift: Celebrating 45 Years of Astrochemistry (in press)

Signatures of an Encounter Between the G2 Cloud and a Jet from Sgr A*

Authors: F. Yusef-Zadeh, M. Wardle

Abstract: The recent discovery of the G2 cloud of dense, ionized gas on a trajectory toward Sgr A*, the black hole at the dynamical center of the Galaxy, offers a unique opportunity to observe an accretion event onto a massive black hole as well as to probe its immediate environment. Simulations and models predict increased X-ray and radio variability resulting from increased accretion driven by drag on an… ▽ More The recent discovery of the G2 cloud of dense, ionized gas on a trajectory toward Sgr A*, the black hole at the dynamical center of the Galaxy, offers a unique opportunity to observe an accretion event onto a massive black hole as well as to probe its immediate environment. Simulations and models predict increased X-ray and radio variability resulting from increased accretion driven by drag on an atmosphere of hot, X-ray emitting gas surrounding Sgr A*. Here, we present X-ray and radio light curves of the emission resulting from the potential encounter of the G2 cloud with a relativistic jet from Sgr A*. This interaction would violently shock a portion of the G2 cloud to temperatures $\sim 10^8$\,K resulting in bright X-ray emission from the dense, shocked gas as it adiabatically expands. The 2-10\, keV luminosity may reach $\sim$10 times the quiescent X-ray flux of Sgr A*. approximately $3\,\lsol$ is emitted above 10\,keV at the peak of the light curve, with significant softening of the spectrum occurring as the gas subsequently cools. Observations with NuSTAR would therefore be able to confirm such an event as well as determine the cloud speed. At radio wavelengths, the associated synchrotron radio emission may reach levels of a few Jy. △ Less

Submitted 24 May, 2013; v1 submitted 16 March, 2013; originally announced March 2013.

Comments: 14 pages, 4 figures, ApJL (in press), a new added figure shows the probability of the encounter between G2 and a jet from five different measurements

ALMA Observations of the Galactic Center: SiO Outflows and High Mass Star Formation near Sgr A*

Authors: F. Yusef-Zadeh, M. Royster, M. Wardle, R. Arendt, H. Bushouse, D. C. Lis, M. W. Pound, D. A. Roberts, B. Whitney, A. Wootten

Abstract: ALMA observations of the Galactic center with spatial resolution $2.61"\times0.97"$ resulted in the detection of 11 SiO (5-4) clumps of molecular gas within 0.6pc (15$"$) of Sgr A*, interior to the 2-pc circumnuclear molecular ring. The three SiO (5-4) clumps closest to Sgr A* show the largest central velocities, $\sim150$ \kms, and broadest asymmetric linewidths with full width zero intensity (FW… ▽ More ALMA observations of the Galactic center with spatial resolution $2.61"\times0.97"$ resulted in the detection of 11 SiO (5-4) clumps of molecular gas within 0.6pc (15$"$) of Sgr A*, interior to the 2-pc circumnuclear molecular ring. The three SiO (5-4) clumps closest to Sgr A* show the largest central velocities, $\sim150$ \kms, and broadest asymmetric linewidths with full width zero intensity (FWZI) $\sim110-147$ \kms. The remaining clumps, distributed mainly to the NE of the ionized mini-spiral, have narrow FWZI ($\sim18-56$ \kms). Using CARMA SiO (2-1) data, LVG modeling of the the SiO line ratios for the broad velocity clumps, constrains the column density N(SiO) $\sim10^{14}$ cm$^{-2}$, and the H$_2$ gas density n$_{\rm H_2}=(3-9)\times10^5$ cm$^{-3}$ for an assumed kinetic temperature 100-200K. The SiO clumps are interpreted as highly embedded protostellar outflows, signifying an early stage of massive star formation near Sgr A* in the last $10^4-10^5$ years. Support for this interpretation is provided by the SiO (5-4) line luminosities and velocity widths which lie in the range measured for protostellar outflows in star forming regions in the Galaxy. Furthermore, SED modeling of stellar sources shows two YSO candidates near SiO clumps, supporting in-situ star formation near Sgr A*. We discuss the nature of star formation where the gravitational potential of the black hole dominates. In particular, we suggest that external radiative pressure exerted on self-shielded molecular clouds enhances the gas density, before the gas cloud become gravitationally unstable near Sgr A*. Alternatively, collisions between clumps in the ring may trigger gravitational collapse. △ Less

Submitted 14 March, 2013; originally announced March 2013.

Comments: 14 pages, 4 figures, ApJL (in press)