Showing 1–50 of 109 results for author: Keto, E

Velocity Analysis of Moving Objects in Earth Observation Satellite Images Using Multi-Spectral Push Broom Scanning

Authors: Eric Keto, Wesley Andres Watters

Abstract: In this study, we present a method for detecting and analyzing the velocities of moving objects in Earth observation satellite images, specifically using data from Planet Labs' push broom scanning satellites. By exploiting the sequential acquisition of multi-spectral images, we estimate the relative differences in acquisition times between spectral bands. This allows us to determine the velocities… ▽ More In this study, we present a method for detecting and analyzing the velocities of moving objects in Earth observation satellite images, specifically using data from Planet Labs' push broom scanning satellites. By exploiting the sequential acquisition of multi-spectral images, we estimate the relative differences in acquisition times between spectral bands. This allows us to determine the velocities of moving objects, such as aircraft, even without precise timestamp information from the image archive. We validate our method by comparing the velocities of aircraft observed in satellite images with those reported by onboard ADS-B transponders. The results demonstrate the potential, despite challenges posed by proprietary data limitations, of a new, useful application of commercial satellite data originally intended as an ongoing, once-daily survey of single images covering the entire land-area of the Earth. Our approach extends the applicability of satellite survey imagery for dynamic object tracking and contributes to the broader use of commercial satellite data in scientific research. △ Less

Submitted 22 June, 2024; v1 submitted 19 June, 2024; originally announced June 2024.

Comments: Submitted to Remote Sensing Letters

Detection of Moving Objects in Earth Observation Satellite Images: Verification

Authors: Eric Keto, Wesley Andres Watters

Abstract: In multi-spectral images made by Earth observation satellites that use push-broom scanning, such as those operated by Planet Labs Corp., moving objects can be identified by the appearance of the object at a different locations in each spectral band. The apparent velocity can be measured if the relative acquisition time between images in different spectral bands is known to millisecond accuracy. Th… ▽ More In multi-spectral images made by Earth observation satellites that use push-broom scanning, such as those operated by Planet Labs Corp., moving objects can be identified by the appearance of the object at a different locations in each spectral band. The apparent velocity can be measured if the relative acquisition time between images in different spectral bands is known to millisecond accuracy. The images in the Planet Labs archive are mosaics of individual exposures acquired at different times. Thus there is not a unique acquisition time for each spectral band. In an earlier paper, we proposed a method to determine the relative acquisition times from the information in the images themselves. High altitude balloons provide excellent targets to test our proposed method because of their high apparent velocity due to the orbital velocity of the satellite and geometric parallax in images aligned to the level of the ground. We use images of the Chinese balloon that crossed the US in February, 2024 as well as images of an identical balloon over Colombia to test our method. Our proposed method appears to be successful and allows the measurement of the apparent velocity of moving objects from the information available in the archive. △ Less

Submitted 19 June, 2024; originally announced June 2024.

Comments: Submitted to the Journal of Astronomical Instrumentation

Detection of Moving Objects in Earth Observation Satellite Images

Authors: Eric Keto, Wesley Andres Watters

Abstract: Moving objects have characteristic signatures in multi-spectral images made by Earth observation satellites that use push broom scanning. While the general concept is applicable to all satellites of this type, each satellite design has its own unique imaging system and requires unique methods to analyze the characteristic signatures. We assess the feasibility of detecting moving objects and measur… ▽ More Moving objects have characteristic signatures in multi-spectral images made by Earth observation satellites that use push broom scanning. While the general concept is applicable to all satellites of this type, each satellite design has its own unique imaging system and requires unique methods to analyze the characteristic signatures. We assess the feasibility of detecting moving objects and measuring their velocities in one particular archive of satellite images made by Planet Labs Corporation with their constellation of SuperDove satellites. Planet Labs data presents a particular challenge in that the images in the archive are mosaics of individual exposures and therefore do not have unique time stamps. We explain how the timing information can be restored indirectly. Our results indicate that the movement of common transportation vehicles, airplanes, cars, and boats, can be detected and measured. △ Less

Submitted 18 May, 2024; originally announced June 2024.

Journal ref: Journal of Astronomical Instrumentation, 2023, Volume 12, Issue 1, id. 2340007-515

Quasi-equilibrium chemical evolution in starless cores

Authors: Jonathan Rawlings, Eric Keto, Paola Caselli

Abstract: The chemistry of H2O, CO and other small molecular species in an isolated pre-stellar core, L1544, has been assessed in the context of a comprehensive gas-grain chemical model, coupled to an empirically constrained physical/dynamical model. Our main findings are (i) that the chemical network remains in near equilibrium as the core evolves towards star formation and the molecular abundances change… ▽ More The chemistry of H2O, CO and other small molecular species in an isolated pre-stellar core, L1544, has been assessed in the context of a comprehensive gas-grain chemical model, coupled to an empirically constrained physical/dynamical model. Our main findings are (i) that the chemical network remains in near equilibrium as the core evolves towards star formation and the molecular abundances change in response to the evolving physical conditions. The gas-phase abundances at any time can be calculated accurately with equilibrium chemistry, and the concept of chemical clocks is meaningless in molecular clouds with similar conditions and dynamical time scales, and (ii) A comparison of the results of complex and simple chemical networks indicates that the abundances of the dominant oxygen and carbon species, H2O, CO, C, and C+ are reasonably approximated by simple networks. In chemical equilibrium, the time-dependent differential terms vanish and a simple network reduces to a few algebraic equations. This allows rapid calculation of the abundances most responsible for spectral line radiative cooling in molecular clouds with long dynamical time scales. The dust ice mantles are highly structured and the ice layers retain a memory of the gas-phase abundances at the time of their deposition. A complex (gas-phase and gas-grain) chemical structure therefore exists, with cosmic-ray induced processes dominating in the inner regions. The inferred H2O abundance profiles for L1544 require that the outer parts of the core and also any medium exterior to the core are essentially transparent to the interstellar radiation field. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 19 pages, 6 figures

Scales of Stability and Turbulence in the Molecular ISM

Authors: Eric Keto

Abstract: From the archival data of the BU-FCRAO $^{13}$ CO GRS, we measure the radial profiles of column density and turbulent velocity dispersion from the centers of molecular clouds outward into the surrounding diffuse molecular ISM. Averaged across spatial scales, the profiles are consistent with turbulence that is on average in hydrostatic equilibrium. We measure the turbulent kinetic energy (KE) and t… ▽ More From the archival data of the BU-FCRAO $^{13}$ CO GRS, we measure the radial profiles of column density and turbulent velocity dispersion from the centers of molecular clouds outward into the surrounding diffuse molecular ISM. Averaged across spatial scales, the profiles are consistent with turbulence that is on average in hydrostatic equilibrium. We measure the turbulent kinetic energy (KE) and the gravitational potential energy (GE) within and the pressure energy (PE) outside the clouds. The distribution of the sum 2KE-|GE|-PE has a mean near zero indicating approximate virial equilibrium. The average pressure energy is consistent with estimates of the pressure of the multiphase ISM at the Galactic mid-plane. If the dynamical time scale of the turbulence scales with its crossing time, this apparent equilibrium may result from the rapid virialization of the KE and GE with respect to changes in PE. In a snapshot, the clouds appear to be in virial equilibrium within a confining pressure. However, the duration of the equilibrium is approximately the crossing time. An analysis of Larson's scaling relationships of line width and column density finds no correlation. Tthe inference of constant column density with cloud size is a misinterpretation of this lack of correlation. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: Submitted to Astronomische Nachrichten

The Scientific Investigation of Unidentified Aerial Phenomena (UAP) Using Multimodal Ground-Based Observatories

Authors: Wesley Andrés Watters, Abraham Loeb, Frank Laukien, Richard Cloete, Alex Delacroix, Sergei Dobroshinsky, Benjamin Horvath, Ezra Kelderman, Sarah Little, Eric Masson, Andrew Mead, Mitch Randall, Forrest Schultz, Matthew Szenher, Foteini Vervelidou, Abigail White, Angelique Ahlström, Carol Cleland, Spencer Dockal, Natasha Donahue, Mark Elowitz, Carson Ezell, Alex Gersznowicz, Nicholas Gold, Michael G. Hercz , et al. (13 additional authors not shown)

Abstract: (Abridged) Unidentified Aerial Phenomena (UAP) have resisted explanation and have received little formal scientific attention for 75 years. A primary objective of the Galileo Project is to build an integrated software and instrumentation system designed to conduct a multimodal census of aerial phenomena and to recognize anomalies. Here we present key motivations for the study of UAP and address hi… ▽ More (Abridged) Unidentified Aerial Phenomena (UAP) have resisted explanation and have received little formal scientific attention for 75 years. A primary objective of the Galileo Project is to build an integrated software and instrumentation system designed to conduct a multimodal census of aerial phenomena and to recognize anomalies. Here we present key motivations for the study of UAP and address historical objections to this research. We describe an approach for highlighting outlier events in the high-dimensional parameter space of our census measurements. We provide a detailed roadmap for deciding measurement requirements, as well as a science traceability matrix (STM) for connecting sought-after physical parameters to observables and instrument requirements. We also discuss potential strategies for deciding where to locate instruments for development, testing, and final deployment. Our instrument package is multimodal and multispectral, consisting of (1) wide-field cameras in multiple bands for targeting and tracking of aerial objects and deriving their positions and kinematics using triangulation; (2) narrow-field instruments including cameras for characterizing morphology, spectra, polarimetry, and photometry; (3) passive multistatic arrays of antennas and receivers for radar-derived range and kinematics; (4) radio spectrum analyzers to measure radio and microwave emissions; (5) microphones for sampling acoustic emissions in the infrasonic through ultrasonic frequency bands; and (6) environmental sensors for characterizing ambient conditions (temperature, pressure, humidity, and wind velocity), as well as quasistatic electric and magnetic fields, and energetic particles. The use of multispectral instruments and multiple sensor modalities will help to ensure that artifacts are recognized and that true detections are corroborated and verifiable. △ Less

Submitted 31 May, 2023; v1 submitted 29 May, 2023; originally announced May 2023.

Comments: This paper is published in the Journal of Astronomical Instrumentation, 12(1), 2340006 (2023) https://doi.org/10.1142/S2251171723400068

Journal ref: Journal of Astronomical Instrumentation, 12(1), 2340006 (2023)

A Hardware and Software Platform for Aerial Object Localization

Authors: Matthew Szenher, Alex Delacroix, Eric Keto, Sarah Little, Mitch Randall, Wesley Andrés Watters, Eric Masson, Richard Cloete

Abstract: To date, there are little reliable data on the position, velocity and acceleration characteristics of Unidentified Aerial Phenomena (UAP). The dual hardware and software system described in this document provides a means to address this gap. We describe a weatherized multi-camera system which can capture images in the visible, infrared and near infrared wavelengths. We then describe the software w… ▽ More To date, there are little reliable data on the position, velocity and acceleration characteristics of Unidentified Aerial Phenomena (UAP). The dual hardware and software system described in this document provides a means to address this gap. We describe a weatherized multi-camera system which can capture images in the visible, infrared and near infrared wavelengths. We then describe the software we will use to calibrate the cameras and to robustly localize objects-of-interest in three dimensions. We show how object localizations captured over time will be used to compute the velocity and acceleration of airborne objects. △ Less

Submitted 29 May, 2023; originally announced May 2023.

Journal ref: Journal of Astronomical Instrumentation, 12(1), 2340002 (2023)

Kinematics and stability of high-mass protostellar disk candidates at sub-arcsecond resolution -- Insights from the IRAM NOEMA large program CORE

Authors: Aida Ahmadi, H. Beuther, F. Bosco, C. Gieser, S. Suri, J. C. Mottram, R. Kuiper, Th. Henning, Á. Sánchez-Monge, H. Linz, R. E. Pudritz, D. Semenov, J. M. Winters, T. Möller, M. T. Beltrán, T. Csengeri, R. Galván-Madrid, K. G. Johnston, E. Keto, P. D. Klaassen, S. Leurini, S. N. Longmore, S. L. Lumsden, L. T. Maud, L. Moscadelli , et al. (6 additional authors not shown)

Abstract: The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in… ▽ More The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in the 1.37 mm wavelength regime in both continuum and line emission, reaching $\sim$0.4" resolution (800 au at 2 kpc). Using the dense gas tracer CH$_3$CN, we find velocity gradients across 13 cores perpendicular to the directions of bipolar molecular outflows, making them excellent disk candidates. Specific angular momentum ($j$) radial profiles are on average $\sim10^{-3}$ km /s pc and follow $j \propto r^{1.7}$, consistent with a poorly resolved rotating and infalling envelope/disk model. Fitting the velocity profiles with a Keplerian model, we find protostellar masses in the range of $\sim 10-25$ $M_\odot$. Modelling the level population of CH$_3$CN lines, we present temperature maps and find median gas temperatures in the range $70-210$ K. We create Toomre $Q$ maps to study the stability of the disks and find almost all (11 of 13) disk candidates to be prone to fragmentation due to gravitational instabilities at the scales probed by our observations. In particular, disks with masses greater than $\sim10-20\%$ of the mass of their host (proto)stars are Toomre unstable, and more luminous protostellar objects tend to have disks that are more massive and hence more prone to fragmentation. Our finings show that most disks around high-mass protostars are prone to disk fragmentation early in their formation due to their high disk to stellar mass ratio. This impacts the accretion evolution of high-mass protostars which will have significant implications for the formation of the most massive stars. △ Less

Submitted 3 May, 2023; v1 submitted 28 April, 2023; originally announced May 2023.

Comments: 27 pages, 12 figures, 6 appendices - accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 677, A171 (2023)

CMZoom III: Spectral Line Data Release

Authors: Daniel Callanan, Steven N. Longmore, Cara Battersby, H. Perry Hatchfield, Daniel L. Walker, Jonathan Henshaw, Eric Keto, Ashley Barnes, Adam Ginsburg, Jens Kauffmann, Diederik Kruijssen, Xing Lu, Elisabeth A. C. Mills, Thushara Pillai, Qizhou Zhang, John Bally, Natalie Butterfield, Yanett A. Contreras, Luis C. Ho, Katharina Immer, Katharine G. Johnston, Juergen Ott, Nimesh Patel, Volker Tolls

Abstract: We present an overview and data release of the spectral line component of the SMA Large Program, \textit{CMZoom}. \textit{CMZoom} observed $^{12}$CO(2-1), $^{13}$CO(2-1) and C$^{18}$O(2-1), three transitions of H$_{2}$CO, several transitions of CH$_{3}$OH, two transitions of OCS and single transitions of SiO and SO, within gas above a column density of N(H$_2$)$\ge 10^{23}$\,cm$^{-2}$ in the Centr… ▽ More We present an overview and data release of the spectral line component of the SMA Large Program, \textit{CMZoom}. \textit{CMZoom} observed $^{12}$CO(2-1), $^{13}$CO(2-1) and C$^{18}$O(2-1), three transitions of H$_{2}$CO, several transitions of CH$_{3}$OH, two transitions of OCS and single transitions of SiO and SO, within gas above a column density of N(H$_2$)$\ge 10^{23}$\,cm$^{-2}$ in the Central Molecular Zone (CMZ; inner few hundred pc of the Galaxy). We extract spectra from all compact 1.3\,mm \emph{CMZoom} continuum sources and fit line profiles to the spectra. We use the fit results from the H$_{2}$CO 3(0,3)-2(0,2) transition to determine the source kinematic properties. We find $\sim 90$\% of the total mass of \emph{CMZoom} sources have reliable kinematics. Only four compact continuum sources are formally self-gravitating. The remainder are consistent with being in hydrostatic equilibrium assuming that they are confined by the high external pressure in the CMZ. Based on the mass and density of virially bound sources, and assuming star formation occurs within one free-fall time with a star formation efficiency of $10\% - 75\%$, we place a lower limit on the future embedded star-formation rate of $0.008 - 0.06$\,M$_{\odot}$\,yr$^{-1}$. We find only two convincing proto-stellar outflows, ruling out a previously undetected population of very massive, actively accreting YSOs with strong outflows. Finally, despite having sufficient sensitivity and resolution to detect high-velocity compact clouds (HVCCs), which have been claimed as evidence for intermediate mass black holes interacting with molecular gas clouds, we find no such objects across the large survey area. △ Less

Submitted 11 January, 2023; originally announced January 2023.

Comments: 44 pages, 41 figures

Clustered Formation of Massive Stars within an Ionized Rotating Disk

Authors: Roberto Galván-Madrid, Qizhou Zhang, Andrés Izquierdo, Charles J. Law, Thomas Peters, Eric Keto, Hauyu Baobab Liu, Paul T. P. Ho, Adam Ginsburg, Carlos Carrasco-González

Abstract: We present ALMA observations with a 800 au resolution and radiative-transfer modelling of the inner part ($r\approx6000$ au) of the ionized accretion flow around a compact star cluster in formation at the center of the luminous ultra-compact (UC) HII region G10.6-0.4. We modeled the flow with an ionized Keplerian disk with and without radial motions in its outer part, or with an external Ulrich en… ▽ More We present ALMA observations with a 800 au resolution and radiative-transfer modelling of the inner part ($r\approx6000$ au) of the ionized accretion flow around a compact star cluster in formation at the center of the luminous ultra-compact (UC) HII region G10.6-0.4. We modeled the flow with an ionized Keplerian disk with and without radial motions in its outer part, or with an external Ulrich envelope. The MCMC fits to the data give total stellar masses $M_\star$ from 120 to $200~M_\odot$, with much smaller ionized-gas masses $M_\mathrm{ion-gas} = 0.2$ to $0.25~M_\odot$. The stellar mass is distributed within the gravitational radius $R_g\approx 1000$ to 1500 au, where the ionized gas is bound. The viewing inclination angle from the face-on orientation is $i = 49$ to $56~°$. Radial motions at radii $r > R_g$ converge to $v_{r,0} \approx 8.7$ km/s, or about the speed of sound of ionized gas, indicating that this gas is marginally unbound at most. From additional constraints on the ionizing-photon rate and far-IR luminosity of the region, we conclude that the stellar cluster consists of a few massive stars with $M_\mathrm{star} = 32$ to $60~M_\odot$, or one star in this range of masses accompanied by a population of lower-mass stars. Any active accretion of ionized gas onto the massive (proto)stars is residual. The inferred cluster density is very large, comparable to that reported at similar scales in the Galactic Center. Stellar interactions are likely to occur within the next Myr. △ Less

Submitted 11 December, 2022; originally announced December 2022.

Comments: Accepted to ApJ Letters, December 7 2022

The Central 1000 au of a Pre-stellar Core Revealed with ALMA. II. Almost Complete Freeze-out

Authors: Paola Caselli, Jaime E. Pineda, Olli Sipilä, Bo Zhao, Elena Redaelli, Silvia Spezzano, Maria José Maureira, Felipe Alves, Luca Bizzocchi, Tyler L. Bourke, Ana Chacón-Tanarro, Rachel Friesen, Daniele Galli, Jorma Harju, Izaskun Jiménez-Serra, Eric Keto, Zhi-Yun Li, Marco Padovani, Anika Schmiedeke, Mario Tafalla, Charlotte Vastel

Abstract: Pre-stellar cores represent the initial conditions in the process of star and planet formation. Their low temperatures ($<$10 K) allow the formation of thick icy dust mantles, which will be partially preserved in the future protoplanetary disks, ultimately affecting the chemical composition of planetary systems. Previous observations have shown that carbon- and oxygen-bearing species, in particula… ▽ More Pre-stellar cores represent the initial conditions in the process of star and planet formation. Their low temperatures ($<$10 K) allow the formation of thick icy dust mantles, which will be partially preserved in the future protoplanetary disks, ultimately affecting the chemical composition of planetary systems. Previous observations have shown that carbon- and oxygen-bearing species, in particular CO, are heavily depleted in pre-stellar cores due to the efficient molecular freeze-out onto the surface of cold dust grains. However, N-bearing species such as NH$_3$ and, in particular, its deuterated isotopologues, appear to maintain high abundances where CO molecules are mainly in solid phase. Thanks to ALMA, we present here the first clear observational evidence of NH$_2$D freeze-out toward the L1544 pre-stellar core, suggestive of the presence of a"complete-depletion zone" within a $\simeq$1800 au radius, in agreement with astrochemical pre-stellar core model predictions. Our state-of-the-art chemical model coupled with a non-LTE radiative transfer code demonstrates that NH$_2$D becomes mainly incorporated in icy mantles in the central 2000 au and starts freezing-out already at $\simeq$7000 au. Radiative transfer effects within the pre-stellar core cause the NH$_2$D(1$_{11}$-1$_{01}$) emission to appear centrally concentrated, with a flattened distribution within the central $\simeq$3000 au, unlike the 1.3 mm dust continuum emission which shows a clear peak within the central $\simeq$1800 au. This prevented NH$_2$D freeze-out to be detected in previous observations, where the central 1000 au cannot be spatially resolved. △ Less

Submitted 27 February, 2022; originally announced February 2022.

Comments: accepted for publication in The Astrophysical Journal

Stability of an Ionization Front in Bondi Accretion

Authors: Eric Keto, Rolf Kuiper

Abstract: Spherical Bondi accretion is used in astrophysics as an approximation to investigate many types of accretion processes. Two-phase accretion flows that transition from neutral to ionized have observational support in high-mass star formation, and have application to accretion flows around any ionizing source, but the hydrodynamic stability of two-phase Bondi accretion is not understood. With both s… ▽ More Spherical Bondi accretion is used in astrophysics as an approximation to investigate many types of accretion processes. Two-phase accretion flows that transition from neutral to ionized have observational support in high-mass star formation, and have application to accretion flows around any ionizing source, but the hydrodynamic stability of two-phase Bondi accretion is not understood. With both semi-analytic and fully numerical methods we find that these flows may be stable, conditionally stable or unstable depending on the initial conditions. The transition from an R-type to a D-type ionization front plays a key role in conditionally stable and unstable flows. △ Less

Submitted 15 October, 2021; originally announced October 2021.

Comments: Accepted MNRAS; 8 pages, 5 figures

Sub-arcsecond Imaging of the Complex Organic Chemistry in Massive Star-forming Region G10.6-0.4

Authors: Charles J. Law, Qizhou Zhang, Karin I. Öberg, Roberto Galván-Madrid, Eric Keto, Hauyu Baobab Liu, Paul T. P. Ho

Abstract: Massive star-forming regions exhibit an extremely rich and diverse chemistry, which in principle provides a wealth of molecular probes, as well as laboratories for interstellar prebiotic chemistry. Since the chemical structure of these sources displays substantial spatial variation among species on small scales (${\lesssim}10^4$ au), high angular resolution observations are needed to connect chemi… ▽ More Massive star-forming regions exhibit an extremely rich and diverse chemistry, which in principle provides a wealth of molecular probes, as well as laboratories for interstellar prebiotic chemistry. Since the chemical structure of these sources displays substantial spatial variation among species on small scales (${\lesssim}10^4$ au), high angular resolution observations are needed to connect chemical structures to local environments and inform astrochemical models of massive star formation. To address this, we present ALMA 1.3 mm observations toward OB cluster-forming region G10.6-0.4 (hereafter "G10.6") at a resolution of 0.14$^{\prime\prime}$ (700 au). We find highly-structured emission from complex organic molecules (COMs) throughout the central 20,000 au, including two hot molecular cores and several shells or filaments. We present spatially-resolved rotational temperature and column density maps for a large sample of COMs and warm gas tracers. These maps reveal a range of gas substructure in both O- and N-bearing species. We identify several spatial correlations that can be explained by existing models of COM formation, including NH$_2$CHO/HNCO and CH$_3$OCHO/CH$_3$OCH$_3$, but also observe unexpected distributions and correlations which suggest that our current understanding of COM formation is far from complete. Importantly, complex chemistry is observed throughout G10.6, rather than being confined to hot cores. The COM composition appears to be different in the cores compared to the more extended structures, which illustrates the importance of high spatial resolution observations of molecular gas in elucidating the physical and chemical processes associated with massive star formation. △ Less

Submitted 19 January, 2021; originally announced January 2021.

Comments: 41 pages, 25 figures, accepted for publication in ApJ

An observational correlation between magnetic field, angular momentum and fragmentation in the envelopes of Class 0 protostars?

Authors: Maud Galametz, Anaelle Maury, Josep M. Girart, Ramprasad Rao, Qizhou Zhang, Mathilde Gaudel, Valeska Valdivia, Patrick Hennebelle, Victoria Cabedo-Soto, Eric Keto, Shih-Ping Lai

Abstract: To assess the potential role of magnetic fields in regulating the envelope rotation and the fragmentation of Class 0 protostars, we carried out observations of the dust polarized emission at 0.87 mm with the SMA, in the envelopes of a large sample of 20 Class 0 protostars. We estimate the mean magnetic field orientation over the central 1000 au envelope scales and compared it to that of the protos… ▽ More To assess the potential role of magnetic fields in regulating the envelope rotation and the fragmentation of Class 0 protostars, we carried out observations of the dust polarized emission at 0.87 mm with the SMA, in the envelopes of a large sample of 20 Class 0 protostars. We estimate the mean magnetic field orientation over the central 1000 au envelope scales and compared it to that of the protostellar outflow in order to study the relation between their misalignment and the kinematics of the circumstellar gas. We discover a strong relationship between the misalignment of the magnetic field orientation with the outflow and the amount of angular momentum observed at similar scales in the protostellar envelope, revealing a potential link between the kinetic and the magnetic energy at envelope scales. The relation could be driven by favored B misalignments in more dynamical envelopes or a dependence of the envelope dynamics with the large-scale B initial configuration. Comparing the trend with the presence of fragmentation, we observe that single sources are mostly associated with conditions of low angular momentum in the inner envelope and good alignment of the magnetic field with protostellar outflows, at intermediate scales. Our results suggest that the properties of the magnetic field in protostellar envelopes bear a tight relationship with the rotating-infalling gas directly involved in the star and disk formation: we find that it may not only influence the fragmentation of protostellar cores into multiple stellar systems, but also set the conditions establishing the pristine properties of planet-forming disks. △ Less

Submitted 27 January, 2021; v1 submitted 23 October, 2020; originally announced October 2020.

Comments: 21 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 644, A47 (2020)

CMZoom II: Catalog of Compact Submillimeter Dust Continuum Sources in the Milky Way's Central Molecular Zone

Authors: H Perry Hatchfield, Cara Battersby, Eric Keto, Daniel Walker, Ashley Barnes, Daniel Callanan, Adam Ginsburg, Jonathan D. Henshaw, Jens Kauffmann, J. M. Diederik Kruijssen, Steve N. Longmore, Xing Lu, Elisabeth A. C. Mills, Thushara Pillai, Qizhou Zhang, John Bally, Natalie Butterfield, Yanett A. Contreras, Luis C. Ho, Jürgen Ott, Nimesh Patel, Volker Tolls

Abstract: In this paper we present the CMZoom Survey's catalog of compact sources (< 10'', ~0.4pc) within the Central Molecular Zone (CMZ). CMZoom is a Submillimeter Array (SMA) large program designed to provide a complete and unbiased map of all high column density gas (N(H$_2$) $\geq$ 10$^{23}$ cm$^{-2}$) of the innermost 500pc of the Galaxy in the 1.3mm dust continuum. We generate both a robust catalog d… ▽ More In this paper we present the CMZoom Survey's catalog of compact sources (< 10'', ~0.4pc) within the Central Molecular Zone (CMZ). CMZoom is a Submillimeter Array (SMA) large program designed to provide a complete and unbiased map of all high column density gas (N(H$_2$) $\geq$ 10$^{23}$ cm$^{-2}$) of the innermost 500pc of the Galaxy in the 1.3mm dust continuum. We generate both a robust catalog designed to reduce spurious source detections, and a second catalog with higher completeness, both generated using a pruned dendrogram. In the robust catalog, we report 285 compact sources, or 816 in the high completeness catalog. These sources have effective radii between 0.04-0.4 pc, and are the potential progenitors of star clusters. The masses for both catalogs are dominated by the Sagittarius B2 cloud complex, where masses are likely unreliable due to free-free contamination, uncertain dust temperatures, and line-of-sight confusion. Given the survey selection and completeness, we predict that our robust catalog accounts for more than ~99% of compact substructure capable of forming high mass stars in the CMZ. This catalog provides a crucial foundation for future studies of high-mass star formation in the Milky Way's Galactic Center. △ Less

Submitted 18 December, 2020; v1 submitted 10 September, 2020; originally announced September 2020.

Comments: 34 pages, 27 figures

Journal ref: ApJS 251 14 (2020)

CMZoom: Survey Overview and First Data Release

Authors: Cara Battersby, Eric Keto, Daniel Walker, Ashley Barnes, Daniel Callanan, Adam Ginsburg, H Perry Hatchfield, Jonathan Henshaw, Jens Kauffmann, J. M. Diederik Kruijssen, Steven N. Longmore, Xing Lu, Elisabeth A. C. Mills, Thushara Pillai, Qizhou Zhang, John Bally, Natalie Butterfield, Yanett A. Contreras, Luis C. Ho, Jurgen Ott, Nimesh Patel, Volker Tolls

Abstract: We present an overview of the CMZoom survey and its first data release. CMZoom is the first blind, high-resolution survey of the Central Molecular Zone (CMZ; the inner 500 pc of the Milky Way) at wavelengths sensitive to the pre-cursors of high-mass stars. CMZoom is a 500-hour Large Program on the Submillimeter Array (SMA) that mapped at 1.3 mm all of the gas and dust in the CMZ above a molecular… ▽ More We present an overview of the CMZoom survey and its first data release. CMZoom is the first blind, high-resolution survey of the Central Molecular Zone (CMZ; the inner 500 pc of the Milky Way) at wavelengths sensitive to the pre-cursors of high-mass stars. CMZoom is a 500-hour Large Program on the Submillimeter Array (SMA) that mapped at 1.3 mm all of the gas and dust in the CMZ above a molecular hydrogen column density of 10^23 cm^-2 at a resolution of ~3" (0.1 pc). In this paper, we focus on the 1.3 mm dust continuum and its data release, but also describe CMZoom spectral line data which will be released in a forthcoming publication. While CMZoom detected many regions with rich and complex substructure, its key result is an overall deficit in compact substructures on 0.1 - 2 pc scales (the compact dense gas fraction: CDGF). In comparison with clouds in the Galactic disk, the CDGF in the CMZ is substantially lower, despite having much higher average column densities. CMZ clouds with high CDGFs are well-known sites of active star formation. The inability of most gas in the CMZ to form compact substructures is likely responsible for the dearth of star formation in the CMZ, surprising considering its high density. The factors responsible for the low CDGF are not yet understood but are plausibly due to the extreme environment of the CMZ, having far-reaching ramifications for our understanding of the star formation process across the cosmos. △ Less

Submitted 26 August, 2020; v1 submitted 9 July, 2020; originally announced July 2020.

Comments: Accepted for Publication in ApJS

Journal ref: 2020 ApJS 249 35

Stability and Solution of the Time-Dependent Bondi-Parker Flow

Authors: Eric Keto

Abstract: Bondi (1952) and Parker (1958} derived a steady-state solution for Bernouilli's equation in spherical symmetry around a point mass for two cases, respectively, an inward accretion flow and an outward wind. Left unanswered were the stability of the steady-state solution, the solution itself of time-dependent flows, whether the time-dependent flows would evolve to the steady-state, and under what co… ▽ More Bondi (1952) and Parker (1958} derived a steady-state solution for Bernouilli's equation in spherical symmetry around a point mass for two cases, respectively, an inward accretion flow and an outward wind. Left unanswered were the stability of the steady-state solution, the solution itself of time-dependent flows, whether the time-dependent flows would evolve to the steady-state, and under what conditions a transonic flow would develop. In a Hamiltonian description, we find that the steady state solution is equivalent to the Lagrangian implying that time-dependent flows evolve to the steady state. We find that the second variation is definite in sign for isothermal and adiabatic flows, implying at least linear stability. We solve the partial differential equation for the time-dependent flow as an initial-value problem and find that a transonic flow develops under a wide range of realistic initial conditions. We present some examples of time-dependent solutions. △ Less

Submitted 20 February, 2020; originally announced February 2020.

Comments: To be published in MNRAS

A Turbulent-Entropic Instability and the Fragmentation of Star-Forming Clouds

Authors: Eric Keto, George B. Field, Eric G. Blackman

Abstract: The kinetic energy of supersonic turbulence within interstellar clouds is subject to cooling by dissipation in shocks and subsequent line radiation. The clouds are therefore susceptible to a condensation process controlled by the specific entropy. In a form analogous to the thermodynamic entropy, the entropy for supersonic turbulence is proportional to the log of the product of the mean turbulent… ▽ More The kinetic energy of supersonic turbulence within interstellar clouds is subject to cooling by dissipation in shocks and subsequent line radiation. The clouds are therefore susceptible to a condensation process controlled by the specific entropy. In a form analogous to the thermodynamic entropy, the entropy for supersonic turbulence is proportional to the log of the product of the mean turbulent velocity and the size scale. We derive a dispersion relation for the growth of entropic instabilities in a spherical self-gravitating cloud and find that there is a critical maximum dissipation time scale, about equal to the crossing time, that allows for fragmentation and subsequent star formation. However, the time scale for the loss of turbulent energy may be shorter or longer, for example with rapid thermal cooling or the injection of mechanical energy. Differences in the time scale for energy loss in different star-forming regions may result in differences in the outcome, for example, in the initial mass function. △ Less

Submitted 8 January, 2020; originally announced January 2020.

Comments: To be published in MNRAS

Massive star formation via torus accretion: the effect of photoionization feedback

Authors: N. S. Sartorio, B. Vandenbroucke, D. Falceta-Goncalves, K. Wood, E. Keto

Abstract: The formation of massive stars is a long standing problem. Although a number of theories of massive star formation exist, ideas appear to converge to a disk-mediated accretion scenario. Here we present radiative hydrodynamic simulations of a star accreting mass via a disk embedded in a torus. We use a Monte Carlo based radiation hydrodynamics code to investigate the impact that ionizing radiation… ▽ More The formation of massive stars is a long standing problem. Although a number of theories of massive star formation exist, ideas appear to converge to a disk-mediated accretion scenario. Here we present radiative hydrodynamic simulations of a star accreting mass via a disk embedded in a torus. We use a Monte Carlo based radiation hydrodynamics code to investigate the impact that ionizing radiation has on the torus. Ionized regions in the torus midplane are found to be either gravitationally trapped or in pressure driven expansion depending on whether or not the size of the ionized region exceeds a critical radius. Trapped Hii regions in the torus plane allow accretion to progress, while expanding Hii regions disrupt the accretion torus preventing the central star from aggregating more mass, thereby setting the star's final mass. We obtain constraints for the luminosities and torus densities that lead to both scenarios. △ Less

Submitted 26 April, 2019; originally announced April 2019.

Comments: 15 pages, 15 figures, accepted for publication in MNRAS

Radiation hydrodynamic simulations of massive star formation via gravitationally trapped HII regions - Spherically symmetric ionised accretion flows

Authors: Kristin Lund, Kenneth Wood, Diego Falceta-Gonçalves, Bert Vandenbroucke, Nina Sartorio, Ian Bonnell, Katharine Johnston, Eric Keto

Abstract: This paper investigates the gravitational trapping of HII regions predicted by steady-state analysis using radiation hydrodynamical simulations. We present idealised spherically symmetric radiation hydrodynamical simulations of the early evolution of HII regions including the gravity of the central source. As with analytic steady state solutions of spherically symmetric ionised Bondi accretion flo… ▽ More This paper investigates the gravitational trapping of HII regions predicted by steady-state analysis using radiation hydrodynamical simulations. We present idealised spherically symmetric radiation hydrodynamical simulations of the early evolution of HII regions including the gravity of the central source. As with analytic steady state solutions of spherically symmetric ionised Bondi accretion flows, we find gravitationally trapped HII regions with accretion through the ionisation front onto the source. We found that, for a constant ionising luminosity, fluctuations in the ionisation front are unstable. This instability only occurs in this spherically symmetric accretion geometry. In the context of massive star formation, the ionising luminosity increases with time as the source accretes mass. The maximum radius of the recurring HII region increases on the accretion timescale until it reaches the sonic radius, where the infall velocity equals the sound speed of the ionised gas, after which it enters a pressure-driven expansion phase. This expansion prevents accretion of gas through the ionisation front, the accretion rate onto the star decreases to zero, and it stops growing from accretion. Because of the time required for any significant change in stellar mass and luminosity through accretion our simulations keep both mass and luminosity constant and follow the evolution from trapped to expanding in a piecewise manner. Implications of this evolution of HII regions include a continuation of accretion of material onto forming stars for a period after the star starts to emit ionising radiation, and an extension of the lifetime of ultracompact HII regions. △ Less

Submitted 1 March, 2019; originally announced March 2019.

Comments: Accepted for publication in MNRAS

Testing the stability of supersonic ionized Bondi accretion flows with radiation hydrodynamics

Authors: Bert Vandenbroucke, Nina S. Sartorio, Kenneth Wood, Kristin Lund, Diego Falceta-Gonçalves, Thomas J. Haworth, Ian Bonnell, Eric Keto, Daniel Tootill

Abstract: We investigate the general stability of 1D spherically symmetric ionized Bondi accretion onto a massive object in the specific context of accretion onto a young stellar object. We first derive a new analytic expression for a steady state two temperature solution that predicts the existence of compact and hypercompact HII regions. We then show that this solution is only marginally stable if ionizat… ▽ More We investigate the general stability of 1D spherically symmetric ionized Bondi accretion onto a massive object in the specific context of accretion onto a young stellar object. We first derive a new analytic expression for a steady state two temperature solution that predicts the existence of compact and hypercompact HII regions. We then show that this solution is only marginally stable if ionization is treated self-consistently. This leads to a recurring collapse of the HII region over time. We derive a semi-analytic model to explain this instability, and test it using spatially converged 1D radiation hydrodynamical simulations. We discuss the implications of the 1D instability on 3D radiation hydrodynamics simulations of supersonic accreting flows. △ Less

Submitted 1 March, 2019; originally announced March 2019.

Comments: 13 pages, 13 figures, 1 appendix, accepted for publication in Montly Notices of the Royal Astronomical Society

The central 1000 AUえーゆー of a pre-stellar core revealed with ALMA. I. 1.3 mm continuum observations

Authors: Paola Caselli, Jaime E. Pineda, Bo Zhao, Malcolm C. Walmsley, Eric Keto, Mario Tafalla, Ana Chacon-Tanarro, Tyler L. Bourke, Rachel Friesen, Daniele Galli, Marco Padovani

Abstract: Stars like our Sun form in self-gravitating dense and cold structures within interstellar clouds, called pre-stellar cores. Although much is known about the physical structure of dense clouds just before and soon after the switch-on of a protostar, the central few thousand astronomical units (au) of pre-stellar cores are unexplored. It is within these central regions that stellar systems assemble… ▽ More Stars like our Sun form in self-gravitating dense and cold structures within interstellar clouds, called pre-stellar cores. Although much is known about the physical structure of dense clouds just before and soon after the switch-on of a protostar, the central few thousand astronomical units (au) of pre-stellar cores are unexplored. It is within these central regions that stellar systems assemble and fragmentation may take place, with the consequent formation of binaries and multiple systems. We present ALMA Band 6 observations (ACA and 12m array) of the dust continuum emission of the 8 Msun pre-stellar core L1544, with angular resolution of 2'' x 1.6'' (linear resolution 270 au x 216 au). Within the primary beam, a compact region of 0.1 Msun, which we call a "kernel", has been unveiled. The kernel is elongated, with a central flat zone with radius Rker ~ 10'' (~ 1400 au). The average number density within Rker is ~1 x 10^6 cm^{-3}, with possible local density enhancements. The region within Rker appears to have fragmented, but detailed analysis shows that similar substructure can be reproduced by synthetic interferometric observations of a smooth centrally concentrated dense core with a similar central flat zone. The presence of a smooth kernel within a dense core is in agreement with non-ideal magneto-hydro-dynamical simulations of a contracting cloud core with a peak number density of 1 x 10^7 cm^{-3}. Dense cores with lower central densities are completely filtered out when simulated 12m-array observations are carried out. These observations demonstrate that the kernel of dynamically evolved dense cores can be investigated at high angular resolution with ALMA. △ Less

Submitted 14 February, 2019; originally announced February 2019.

Comments: accepted for publication in The Astrophysical Journal

The dynamical evolution of molecular clouds near the Galactic Centre - II. Spatial structure and kinematics of simulated clouds

Authors: J. M. D. Kruijssen, J. E. Dale, S. N. Longmore, D. L. Walker, J. D. Henshaw, S. M. R. Jeffreson, M. A. Petkova, A. Ginsburg, A. T. Barnes, C. D. Battersby, K. Immer, J. M. Jackson, E. R. Keto, N. Krieger, E. A. C. Mills, Á. Sánchez-Monge, A. Schmiedeke, S. T. Suri, Q. Zhang

Abstract: The evolution of molecular clouds in galactic centres is thought to differ from that in galactic discs due to a significant influence of the external gravitational potential. We present a set of numerical simulations of molecular clouds orbiting on the 100-pc stream of the Central Molecular Zone (the central $\sim500$ pc of the Galaxy) and characterise their morphological and kinematic evolution i… ▽ More The evolution of molecular clouds in galactic centres is thought to differ from that in galactic discs due to a significant influence of the external gravitational potential. We present a set of numerical simulations of molecular clouds orbiting on the 100-pc stream of the Central Molecular Zone (the central $\sim500$ pc of the Galaxy) and characterise their morphological and kinematic evolution in response to the background potential and eccentric orbital motion. We find that the clouds are shaped by strong shear and torques, by tidal and geometric deformation, and by their passage through the orbital pericentre. Within our simulations, these mechanisms control cloud sizes, aspect ratios, position angles, filamentary structure, column densities, velocity dispersions, line-of-sight velocity gradients, spin angular momenta, and kinematic complexity. By comparing these predictions to observations of clouds on the Galactic Centre 'dust ridge', we find that our simulations naturally reproduce a broad range of key observed morphological and kinematic features, which can be explained in terms of well-understood physical mechanisms. We argue that the accretion of gas clouds onto the central regions of galaxies, where the rotation curve turns over and the tidal field is fully compressive, is accompanied by transformative dynamical changes to the clouds, leading to collapse and star formation. This can generate an evolutionary progression of cloud collapse with a common starting point, which either marks the time of accretion onto the tidally-compressive region or of the most recent pericentre passage. Together, these processes may naturally produce the synchronised starbursts observed in numerous (extra)galactic nuclei. △ Less

Submitted 7 February, 2019; v1 submitted 5 February, 2019; originally announced February 2019.

Comments: 21 pages, 8 figures, 3 tables; accepted by MNRAS (February 5, 2019); Figures 2, 3, and 4 show the main results of the paper. An animated version of Figure 2 is available as an ancillary file in the arXiv source and will be included as online Supporting Information with the MNRAS publication

Core fragmentation and Toomre stability analysis of W3(H2O): A case study of the IRAM NOEMA large program CORE

Authors: A. Ahmadi, H. Beuther, J. C. Mottram, F. Bosco, H. Linz, Th. Henning, J. M. Winters, R. Kuiper, R. Pudritz, Á. Sánchez-Monge, E. Keto, M. Beltran, S. Bontemps, R. Cesaroni, T. Csengeri, S. Feng, R. Galvan-Madrid, K. G. Johnston, P. Klaassen, S. Leurini, S. N. Longmore, S. Lumsden, L. T. Maud, K. M. Menten, L. Moscadelli , et al. (8 additional authors not shown)

Abstract: The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.3… ▽ More The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.37 mm wavelength regime in both line and dust continuum emission. We present the spectral line setup of the CORE survey and a case study for W3(H2O). At ~0.35" (700 AUえーゆー at 2 kpc) resolution, the W3(H2O) clump fragments into two cores (West and East), separated by ~2300 AUえーゆー. Velocity shifts of a few km/s are observed in the dense-gas tracer, CH3CN, across both cores, consistent with rotation and perpendicular to the directions of two bipolar outflows, one emanating from each core. The kinematics of the rotating structure about W3(H2O) W shows signs of differential rotation of material, possibly in a disk-like object. The observed rotational signature around W3(H2O) E may be due to a disk-like object, an unresolved binary (or multiple) system, or a combination of both. We fit the emission of CH3CN (12-11) K = 4-6 and derive a gas temperature map with a median temperature of ~165 K across W3(H2O). We create a Toomre Q map to study the stability of the rotating structures against gravitational instability. The rotating structures appear to be Toomre unstable close to their outer boundaries, with a possibility of further fragmentation in the differentially-rotating core W3(H2O) W. Rapid cooling in the Toomre-unstable regions supports the fragmentation scenario. Combining millimeter dust continuum and spectral line data toward the famous high-mass star-forming region W3(H2O), we identify core fragmentation on large scales, and indications for possible disk fragmentation on smaller spatial scales. △ Less

Submitted 1 August, 2018; originally announced August 2018.

Comments: 23 pages, 26 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 618, A46 (2018)

Fragmentation and disk formation during high-mass star formation: The IRAM NOEMA (Northern Extended Millimeter Array) large program CORE

Authors: H. Beuther, J. C. Mottram, A. Ahmadi, F. Bosco, H. Linz, Th. Henning, P. Klaassen, J. M. Winters, L. T. Maud, R. Kuiper, D. Semenov, C. Gieser, T. Peters, J. S. Urquhart, R. Pudritz, S. E. Ragan, S. Feng, E. Keto, S. Leurini, R. Cesaroni, M. Beltran, A. Palau, A. Sanchez-Monge, R. Galvan-Madrid, Q. Zhang , et al. (8 additional authors not shown)

Abstract: Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high… ▽ More Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AUえーゆー. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation. △ Less

Submitted 3 May, 2018; originally announced May 2018.

Comments: 29 pages, 17 figures, Astronomy & Astrophysics in press, for a higher-resolution version of the paper see http://www.mpia.de/homes/beuther/papers.html

Journal ref: A&A 617, A100 (2018)

Radiative transfer modelling of W33A MM1: 3-D structure and dynamics of a complex massive star forming region

Authors: Andrés F. Izquierdo, Roberto Galván-Madrid, Luke T. Maud, Melvin G. Hoare, Katharine G. Johnston, Eric R. Keto, Qizhou Zhang, Willem-Jan de Wit

Abstract: We present a composite model and radiative transfer simulations of the massive star forming core W33A MM1. The model was tailored to reproduce the complex features observed with ALMA at $\approx 0.2$ arcsec resolution in CH$_3$CN and dust emission. The MM1 core is fragmented into six compact sources coexisting within $\sim 1000$ au. In our models, three of these compact sources are better represen… ▽ More We present a composite model and radiative transfer simulations of the massive star forming core W33A MM1. The model was tailored to reproduce the complex features observed with ALMA at $\approx 0.2$ arcsec resolution in CH$_3$CN and dust emission. The MM1 core is fragmented into six compact sources coexisting within $\sim 1000$ au. In our models, three of these compact sources are better represented as disc-envelope systems around a central (proto)star, two as envelopes with a central object, and one as a pure envelope. The model of the most prominent object (Main) contains the most massive (proto)star ($M_\star\approx7~M_\odot$) and disc+envelope ($M_\mathrm{gas}\approx0.4~M_\odot$), and is the most luminous ($L_\mathrm{Main} \sim 10^4~L_\odot$). The model discs are small (a few hundred au) for all sources. The composite model shows that the elongated spiral-like feature converging to the MM1 core can be convincingly interpreted as a filamentary accretion flow that feeds the rising stellar system. The kinematics of this filament is reproduced by a parabolic trajectory with focus at the center of mass of the region. Radial collapse and fragmentation within this filament, as well as smaller filamentary flows between pairs of sources are proposed to exist. Our modelling supports an interpretation where what was once considered as a single massive star with a $\sim 10^3$ au disc and envelope, is instead a forming stellar association which appears to be virialized and to form several low-mass stars per high-mass object. △ Less

Submitted 24 April, 2018; originally announced April 2018.

Comments: 22 pages, 13 figures, 3 tables. Accepted for publication in MNRAS

SMA observations of the polarized dust emission in solar-type Class 0 protostars: the magnetic field properties at envelope scales

Authors: Maud Galametz, Anaelle Maury, Josep M. Girart, Ramprasad Rao, Qizhou Zhang, Mathilde Gaudel, Valeska Valdivia, Eric Keto, Shih-Ping Lai

Abstract: Although, from a theoretical point of view, magnetic fields are believed to have a significant role during the early stages of star formation, especially during the main accretion phase, the magnetic field strength and topology is poorly constrained in the youngest accreting Class 0 protostars that lead to the formation of solar-type stars. We carried out observations of the polarized dust continu… ▽ More Although, from a theoretical point of view, magnetic fields are believed to have a significant role during the early stages of star formation, especially during the main accretion phase, the magnetic field strength and topology is poorly constrained in the youngest accreting Class 0 protostars that lead to the formation of solar-type stars. We carried out observations of the polarized dust continuum emission with the SMA interferometer at 0.87mm, in order to probe the structure of the magnetic field in a sample of 12 low-mass Class 0 envelopes, including both single protostars and multiple systems, in nearby clouds. Our SMA observations probe the envelope emission at scales of 600-5000 au with a spatial resolution ranging from 600 to 1500 au depending on the source distance. We report the detection of linearly polarized dust continuum emission in all of our targets, with average polarization fractions ranging from 2% to 10% in these protostellar envelopes. The polarization fraction decreases with the continuum flux density, which translates into a decrease with the H2 column density within an individual envelope. Our analysis show that the envelope-scale magnetic field is preferentially observed either aligned or perpendicular to the outflow direction. Interestingly, our results suggest for the first time a relation between the orientation of the magnetic field and the rotational energy of envelopes, with a larger occurrence of misalignment in sources where strong rotational motions are detected at hundreds to thousands of au scales. We also show that the best agreement between the magnetic field and outflow orientation is found in sources showing no small-scale multiplicity and no large disks at about 100 au scales. △ Less

Submitted 16 April, 2018; originally announced April 2018.

Comments: 22 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 616, A139 (2018)

Magnetically regulated collapse in the B335 protostar ? I. ALMA observations of the polarized dust emission

Authors: Anaëlle J. Maury, Josep Miquel Girart, Qizhou Zhang, Patrick Hennebelle, Eric Keto, Ramprasad Rao, Shih-Ping Lai, Nagayoshi Ohashi, Maud Galametz

Abstract: The role of the magnetic field during protostellar collapse is poorly constrained from an observational point of view, although it could be significant if we believe state-of-the-art models of protostellar formation. We present polarimetric observations of the 233 GHz thermal dust continuum emission obtained with ALMA in the B335 Class 0 protostar. Linearly polarized dust emission arising from the… ▽ More The role of the magnetic field during protostellar collapse is poorly constrained from an observational point of view, although it could be significant if we believe state-of-the-art models of protostellar formation. We present polarimetric observations of the 233 GHz thermal dust continuum emission obtained with ALMA in the B335 Class 0 protostar. Linearly polarized dust emission arising from the circumstellar material in the envelope of B335 is detected at all scales probed by our observations, from radii of 50 to 1000 au. The magnetic field structure producing the dust polarization has a very ordered topology in the inner envelope, with a transition from a large-scale poloidal magnetic field, in the outflow direction, to strongly pinched in the equatorial direction. This is probably due to magnetic field lines being dragged along the dominating infall direction since B335 does not exhibit prominent rotation. Our data and their qualitative comparison to a family of magnetized protostellar collapse models show that, during the magnetized collapse in B335, the magnetic field is maintaining a high level of organization from scales 1000 au to 50 au: this suggests the field is dynamically relevant and capable of influencing the typical outcome of protostellar collapse, such as regulating the disk size in B335. △ Less

Submitted 28 February, 2018; originally announced March 2018.

Comments: Accepted for publication in MNRAS

Intensity-Corrected Herschel Observations of Nearby Isolated Low-Mass Clouds

Authors: Sarah I. Sadavoy, Eric Keto, Tyler L. Bourke, Michael M. Dunham, Philip C. Myers, Ian W. Stephens, James Di Francesco, Kristi Webb, Amelia Stutz, Ralf Launhardt, John Tobin

Abstract: We present intensity-corrected Herschel maps at 100 um, 160 um, 250 um, 350 um, and 500 um for 56 isolated low-mass clouds. We determine the zero-point corrections for Herschel PACS and SPIRE maps from the Herschel Science Archive (HSA) using Planck data. Since these HSA maps are small, we cannot correct them using typical methods. Here, we introduce a technique to measure the zero-point correctio… ▽ More We present intensity-corrected Herschel maps at 100 um, 160 um, 250 um, 350 um, and 500 um for 56 isolated low-mass clouds. We determine the zero-point corrections for Herschel PACS and SPIRE maps from the Herschel Science Archive (HSA) using Planck data. Since these HSA maps are small, we cannot correct them using typical methods. Here, we introduce a technique to measure the zero-point corrections for small Herschel maps. We use radial profiles to identify offsets between the observed HSA intensities and the expected intensities from Planck. Most clouds have reliable offset measurements with this technique. In addition, we find that roughly half of the clouds have underestimated HSA-SPIRE intensities in their outer envelopes relative to Planck, even though the HSA-SPIRE maps were previously zero-point corrected. Using our technique, we produce corrected Herschel intensity maps for all 56 clouds and determine their line-of-sight average dust temperatures and optical depths from modified black body fits. The clouds have typical temperatures of ~ 14-20 K and optical depths of ~ 1e-5 - 1e-3. Across the whole sample, we find an anti-correlation between temperature and optical depth. We also find lower temperatures than what was measured in previous Herschel studies, which subtracted out a background level from their intensity maps to circumvent the zero-point correction. Accurate Herschel observations of clouds are key to obtain accurate density and temperature profiles. To make such future analyses possible, intensity-corrected maps for all 56 clouds are publicly available in the electronic version. △ Less

Submitted 30 November, 2017; originally announced December 2017.

Comments: Accepted to ApJ, 35 pages, 15 figures, 4 Tables. Intensity-corrected FITS images and temperature/optical depths maps are available with the electronic version

Star formation in a high-pressure environment: An SMA view of the Galactic centre dust ridge

Authors: D. L. Walker, S. N. Longmore, Q. Zhang, C. Battersby, E. Keto, J. M. D. Kruijssen, A. Ginsburg, X. Lu, J. D. Henshaw, J. Kauffmann, T. Pillai, E. A. C. Mills, A. J. Walsh, J. Bally, L. C. Ho, K. Immer, K. G. Johnston

Abstract: The star formation rate in the Central Molecular Zone (CMZ) is an order of magnitude lower than predicted according to star formation relations that have been calibrated in the disc of our own and nearby galaxies. Understanding how and why star formation appears to be different in this region is crucial if we are to understand the environmental dependence of the star formation process. Here, we pr… ▽ More The star formation rate in the Central Molecular Zone (CMZ) is an order of magnitude lower than predicted according to star formation relations that have been calibrated in the disc of our own and nearby galaxies. Understanding how and why star formation appears to be different in this region is crucial if we are to understand the environmental dependence of the star formation process. Here, we present the detection of a sample of high-mass cores in the CMZ's "dust ridge" that have been discovered with the Submillimeter Array as part of the CMZoom survey. These cores range in mass from ~ 50 - 2150 Msun within radii of 0.1 - 0.25 pc. All appear to be young (pre-UCHII), meaning that they are prime candidates for representing the initial conditions of high-mass stars and sub-clusters. We report that at least two of these cores ('c1' and 'e1') contain young, high-mass protostars. We compare all of the detected cores with high-mass cores in the Galactic disc and find that they are broadly similar in terms of their masses and sizes, despite being subjected to external pressures that are several orders of magnitude greater - ~ 10^8 K/cm^3, as opposed to ~ 10^5 K/cm^3. The fact that > 80% of these cores do not show any signs of star-forming activity in such a high-pressure environment leads us to conclude that this is further evidence for an increased critical density threshold for star formation in the CMZ due to turbulence. △ Less

Submitted 2 November, 2017; originally announced November 2017.

Comments: 19 pages, 10 figures, 3 tables. Accepted for publication in MNRAS

NH_3(1_0-0_0) in the pre-stellar core L1544

Authors: P. Caselli, L. Bizzocchi, E. Keto, O. Sipila, M. Tafalla, L. Pagani, L. E. Kristensen, F. F. S. van der Tak, C. M. Walmsley, C. Codella, B. Nisini, Y. Aikawa, A. Faure, E. F. van Dishoeck

Abstract: Pre-stellar cores represent the initial conditions in the process of star and planet formation, therefore it is important to study their physical and chemical structure. Because of their volatility, nitrogen-bearing molecules are key to study the dense and cold gas present in pre-stellar cores. The NH_3 rotational transition detected with Herschel-HIFI provides a unique combination of sensitivity… ▽ More Pre-stellar cores represent the initial conditions in the process of star and planet formation, therefore it is important to study their physical and chemical structure. Because of their volatility, nitrogen-bearing molecules are key to study the dense and cold gas present in pre-stellar cores. The NH_3 rotational transition detected with Herschel-HIFI provides a unique combination of sensitivity and spectral resolution to further investigate physical and chemical processes in pre-stellar cores. Here we present the velocity-resolved Herschel-HIFI observations of the ortho-NH_3(1_0-0_0) line at 572 GHz and study the abundance profile of ammonia across the pre-stellar core L1544 to test current theories of its physical and chemical structure. Recently calculated collisional coefficients have been included in our non-LTE radiative transfer code to reproduce Herschel observations. A gas-grain chemical model, including spin-state chemistry and applied to the (static) physical structure of L1544 is also used to infer the abundance profile of ortho-NH_3 . The hyperfine structure of ortho-NH_3(1_0-0_0) is resolved for the first time in space. All the hyperfine components are strongly self-absorbed. The profile can be reproduced if the core is contracting in quasi-equilibrium, consistent with previous work, and if the NH_3 abundance is slightly rising toward the core centre, as deduced from previous interferometric observations of para-NH_3(1,1). The chemical model overestimates the NH_3 abundance at radii between ~ 4000 and 15000 AUえーゆー by about two orders of magnitude and underestimates the abundance toward the core centre by more than one order of magnitude. Our observations show that chemical models applied to static clouds have problems in reproducing NH_3 observations. △ Less

Submitted 9 June, 2017; originally announced June 2017.

Comments: accepted for publication in A&A Letters

The ALMA view of W33A: A Spiral Filament Feeding the Candidate Disc in MM1-Main

Authors: L. T. Maud, M. G. Hoare, R. Galván-Madrid, Q. Zhang, E. Keto, K. G. Johnston, J. E. Pineda

Abstract: We targeted the massive star forming region W33A using the Atacama Large Sub/Millimeter Array (ALMA) in band 6 (230 GHz) and 7 (345 GHz) to search for a sub-1000au disc around the central O-type massive young stellar object (MYSO) W33A MM1-Main. Our data achieves a resolution of ~0.2" (~500au) and resolves the central core, MM1, into multiple components and reveals complex and filamentary structur… ▽ More We targeted the massive star forming region W33A using the Atacama Large Sub/Millimeter Array (ALMA) in band 6 (230 GHz) and 7 (345 GHz) to search for a sub-1000au disc around the central O-type massive young stellar object (MYSO) W33A MM1-Main. Our data achieves a resolution of ~0.2" (~500au) and resolves the central core, MM1, into multiple components and reveals complex and filamentary structures. There is strong molecular line emission covering the entire MM1 region. The kinematic signatures are inconsistent with only Keplerian rotation although we propose that the shift in the emission line centroids within ~1000au of MM1-Main could hint at an underlying compact disc with Keplerian rotation. We cannot however rule out the possibility of an unresolved binary or multiple system. A putative smaller disc could be fed by the large scale spiral `feeding filament' we detect in both gas and dust emission. We also discuss the nature of the now-resolved continuum sources. △ Less

Submitted 24 January, 2017; originally announced January 2017.

A Brief Update on the CMZoom Survey

Authors: Cara Battersby, Eric Keto, Qizhou Zhang, Steven N. Longmore, J. M. Diederik Kruijssen, Thushara Pillai, Jens Kauffmann, Dan Walker, Xing Lu, Adam Ginsburg, John Bally, Elisabeth A. C. Mills, Jonathan D. Henshaw, Katharina Immer, Nimesh Patel, Volker Tolls, Andrew J. Walsh, Katharine Johnston, Luis C. Ho

Abstract: The inner few hundred parsecs of the Milky Way, the Central Molecular Zone (CMZ), is our closest laboratory for understanding star formation in the extreme environments (hot, dense, turbulent gas) that once dominated the universe. We present an update on the first large-area survey to expose the sites of star formation across the CMZ at high-resolution in submillimeter wavelengths: the CMZoom surv… ▽ More The inner few hundred parsecs of the Milky Way, the Central Molecular Zone (CMZ), is our closest laboratory for understanding star formation in the extreme environments (hot, dense, turbulent gas) that once dominated the universe. We present an update on the first large-area survey to expose the sites of star formation across the CMZ at high-resolution in submillimeter wavelengths: the CMZoom survey with the Submillimeter Array (SMA). We identify the locations of dense cores and search for signatures of embedded star formation. CMZoom is a three-year survey in its final year and is mapping out the highest column density regions of the CMZ in dust continuum and a variety of spectral lines around 1.3 mm. CMZoom combines SMA compact and subcompact configurations with single-dish data from BGPS and the APEX telescope, achieving an angular resolution of about 4" (0.2 pc) and good image fidelity up to large spatial scales. △ Less

Submitted 18 October, 2016; originally announced October 2016.

Comments: 4 pages, Astrophysics of the Galactic Centre, Proceedings IAUえーゆー Symposium No. 322, 2016; R. Crocker, S. Longmore & G. Bicknell. CMZoom website: https://www.cfa.harvard.edu/sma/LargeScale/CMZ/

Fourier-space combination of Planck and Herschel images

Authors: J. Abreu-Vicente, A. Stutz, Th. Henning, E. Keto, J. Ballesteros-Paredes, T. Robitaille

Abstract: Herschel has revolutionized our ability to measure column densities (N$_{\rm H}$) and temperatures (T) of molecular clouds thanks to its far infrared multiwavelength coverage. However, the lack of a well defined background intensity level in the Herschel data limits the accuracy of the N$_{\rm H}$ and T maps. We provide a method that corrects the missing Herschel background intensity levels using… ▽ More Herschel has revolutionized our ability to measure column densities (N$_{\rm H}$) and temperatures (T) of molecular clouds thanks to its far infrared multiwavelength coverage. However, the lack of a well defined background intensity level in the Herschel data limits the accuracy of the N$_{\rm H}$ and T maps. We provide a method that corrects the missing Herschel background intensity levels using the Planck model for foreground Galactic thermal dust emission. We present a Fourier method that combines the publicly available Planck model on large angular scales with the Herschel images on smaller angular scales. We apply our method to two regions spanning a range of Galactic environments: Perseus and the Galactic plane region around $l = 11°$ (HiGal--11). We post-process the combined dust continuum emission images to generate column density and temperature maps. We compare these to previously adopted constant--offset corrections. We find significant differences ($\gtrsim$20\%) over significant ($\sim$15\%) areas of the maps, at low column densities ($N_{\rm H}\lesssim10^{22}$\,cm$^{-2}$) and relatively high temperatures ($T\gtrsim20$\,K). We also apply our method to synthetic observations of a simulated molecular cloud to validate our method. Our method successfully corrects the Herschel images, including both the constant--offset intensity level and the scale-dependent background variations measured by Planck. Our method improves the previous constant--offset corrections, which did not account for variations in the background emission levels. △ Less

Submitted 29 June, 2017; v1 submitted 10 May, 2016; originally announced May 2016.

Comments: Paper accepted for publications in A&A

Journal ref: A&A 604, A65 (2017)

Radiative Transfer of HCN: Interpreting observations of hyperfine anomalies

Authors: A. M. Mullins, R. M. Loughnane, M. P. Redman, B. Wiles, N. Guegan, J. Barrett, E. R. Keto

Abstract: Molecules with hyperfine splitting of their rotational line spectra are useful probes of optical depth, via the relative line strengths of their hyperfine components.The hyperfine splitting is particularly advantageous in interpreting the physical conditions of the emitting gas because with a second rotational transition, both gas density and temperature can be derived. For HCN however, the relati… ▽ More Molecules with hyperfine splitting of their rotational line spectra are useful probes of optical depth, via the relative line strengths of their hyperfine components.The hyperfine splitting is particularly advantageous in interpreting the physical conditions of the emitting gas because with a second rotational transition, both gas density and temperature can be derived. For HCN however, the relative strengths of the hyperfine lines are anomalous. They appear in ratios which can vary significantly from source to source, and are inconsistent with local thermodynamic equilibrium. This is the HCN hyperfine anomaly, and it prevents the use of simple LTE models of HCN emission to derive reliable optical depths. In this paper we demonstrate how to model HCN hyperfine line emission, and derive accurate line ratios, spectral line shapes and optical depths. We show that by carrying out radiative transfer calculations over each hyperfine level individually, as opposed to summing them over each rotational level, the anomalous hyperfine emission emerges naturally. To do this requires not only accurate radiative rates between hyperfine states, but also accurate collisional rates. We investigate the effects of different sets of hyperfine collisional rates, derived via the 'proportional method' and through direct recoupling calculations. Through an extensive parameter sweep over typical low mass star forming conditions, we show the HCN line ratios to be highly variable to optical depth. We also reproduce an observed effect whereby the red-blue asymmetry of the hyperfine lines (an infall signature) switches sense within a single rotational transition. △ Less

Submitted 11 April, 2016; originally announced April 2016.

Comments: 12 pages, 7 figures. Accepted for publication in MNRAS

A Hot and Massive Accretion Disk around the High-Mass Protostar IRAS 20126+4104

Authors: Huei-Ru Vivien Chen, Eric Keto, Qizhou Zhang, T. K. Sridharan, Sheng-Yuan Liu, Yu-Nung Su

Abstract: We present new spectral line observations of the CH3CN molecule in the accretion disk around the massive protostar IRAS 20126+4104 with the Submillimeter Array that for the first time measure the disk density, temperature, and rotational velocity with sufficient resolution (0.37", equivalent to ~600 AUえーゆー) to assess the gravitational stability of the disk through the Toomre-Q parameter. Our observati… ▽ More We present new spectral line observations of the CH3CN molecule in the accretion disk around the massive protostar IRAS 20126+4104 with the Submillimeter Array that for the first time measure the disk density, temperature, and rotational velocity with sufficient resolution (0.37", equivalent to ~600 AUえーゆー) to assess the gravitational stability of the disk through the Toomre-Q parameter. Our observations resolve the central 2000 AUえーゆー region that shows steeper velocity gradients with increasing upper state energy, indicating an increase in the rotational velocity of the hotter gas nearer the star. Such spin-up motions are characteristics of an accretion flow in a rotationally supported disk. We compare the observed data with synthetic image cubes produced by three-dimensional radiative transfer models describing a thin flared disk in Keplerian motion enveloped within the centrifugal radius of an angular-momentum-conserving accretion flow. Given a luminosity of 1.3x10^4 Lsun, the optimized model gives a disk mass of 1.5 Msun and a radius of 858 AUえーゆー rotating about a 12.0 Msun protostar with a disk mass accretion rate of 3.9x10^{-5} Msun/yr. Our study finds that, in contrast to some theoretical expectations, the disk is hot and stable to fragmentation with Q > 2.8 at all radii which permits a smooth accretion flow. These results put forward the first constraints on gravitational instabilities in massive protostellar disks, which are closely connected to the formation of companion stars and planetary systems by fragmentation. △ Less

Submitted 2 April, 2016; originally announced April 2016.

Comments: 24 pages, 6 figures, 3 tables, accepted for publication in ApJ

Evidence of Short Timescale Flux Density Variations of UC HII regions in Sgr B2 Main and North

Authors: C. G. De Pree, T. Peters, M. -M. Mac Low, D. J. Wilner, W. M. Goss, R. Galván-Madrid, E. R. Keto, R. S. Klessen, A. Monsrud

Abstract: We have recently published observations of significant flux density variations at 1.3 cm in HII regions in the star forming regions Sgr B2 Main and North (De Pree et al. 2014). To further study these variations, we have made new 7 mm continuum and recombination line observations of Sgr B2 at the highest possible angular resolution of the Karl G. Jansky Very Large Array (VLA). We have observed Sgr… ▽ More We have recently published observations of significant flux density variations at 1.3 cm in HII regions in the star forming regions Sgr B2 Main and North (De Pree et al. 2014). To further study these variations, we have made new 7 mm continuum and recombination line observations of Sgr B2 at the highest possible angular resolution of the Karl G. Jansky Very Large Array (VLA). We have observed Sgr B2 Main and North at 42.9 GHz and at 45.4 GHz in the BnA configuration (Main) and the A configuration (North). We compare these new data to archival VLA 7 mm continuum data of Sgr B2 Main observed in 2003 and Sgr B2 North observed in 2001. We find that one of the 41 known ultracompact and hypercompact HII regions in Sgr B2 (K2-North) has decreased $\sim$27% in flux density from 142$\pm$14 mJy to 103$\pm$10 mJy (2.3$σしぐま$) between 2001 and 2012. A second source, F3c-Main has increased $\sim$30% in flux density from 82$\pm$8 mJy to 107 $\pm$11 mJy (1.8$σしぐま$) between 2003 and 2012. F3c-Main was previously observed to increase in flux density at 1.3 cm over a longer time period between 1989 and 2012 (De Pree et al. 2014). An observation of decreasing flux density, such as that observed in K2-North, is particularly significant since such a change is not predicted by the classical hypothesis of steady expansion of HII regions during massive star accretion. Our new observations at 7 mm, along with others in the literature, suggest that the formation of massive stars occurs through time-variable and violent accretion. △ Less

Submitted 16 November, 2015; originally announced November 2015.

Comments: 20 pages, 5 figures, Accepted for publication in The Astrophysical Journal

First-generation Science Cases for Ground-based Terahertz Telescopes

Authors: Hiroyuki Hirashita, Patrick M. Koch, Satoki Matsushita, Shigehisa Takakuwa, Masanori Nakamura, Keiichi Asada, Hauyu Baobab Liu, Yuji Urata, Ming-Jye Wang, Wei-Hao Wang, Satoko Takahashi, Ya-Wen Tang, Hsian-Hong Chang, Kuiyun Huang, Oscar Morata, Masaaki Otsuka, Kai-Yang Lin, An-Li Tsai, Yen-Ting Lin, Sundar Srinivasan, Pierre Martin-Cocher, Hung-Yi Pu, Francisca Kemper, Nimesh Patel, Paul Grimes , et al. (11 additional authors not shown)

Abstract: Ground-based observations at terahertz (THz) frequencies are a newly explorable area of astronomy for the next ten years. We discuss science cases for a first-generation 10-m class THz telescope, focusing on the Greenland Telescope as an example of such a facility. We propose science cases and provide quantitative estimates for each case. The largest advantage of ground-based THz telescopes is the… ▽ More Ground-based observations at terahertz (THz) frequencies are a newly explorable area of astronomy for the next ten years. We discuss science cases for a first-generation 10-m class THz telescope, focusing on the Greenland Telescope as an example of such a facility. We propose science cases and provide quantitative estimates for each case. The largest advantage of ground-based THz telescopes is their higher angular resolution (~ 4 arcsec for a 10-m dish), as compared to space or airborne THz telescopes. Thus, high-resolution mapping is an important scientific argument. In particular, we can isolate zones of interest for Galactic and extragalactic star-forming regions. The THz windows are suitable for observations of high-excitation CO lines and [N II] 205 um lines, which are scientifically relevant tracers of star formation and stellar feedback. Those lines are the brightest lines in the THz windows, so that they are suitable for the initiation of ground-based THz observations. THz polarization of star-forming regions can also be explored since it traces the dust population contributing to the THz spectral peak. For survey-type observations, we focus on ``sub-THz'' extragalactic surveys, whose uniqueness is to detect galaxies at redshifts z ~ 1--2, where the dust emission per comoving volume is the largest in the history of the Universe. Finally we explore possibilities of flexible time scheduling, which enables us to monitor active galactic nuclei, and to target gamma-ray burst afterglows. For these objects, THz and submillimeter wavelength ranges have not yet been explored. △ Less

Submitted 3 November, 2015; originally announced November 2015.

Comments: 39 pages, 17 figures, accepted for publication in PASJ (as a review paper)

SOFIA/FORCAST Observations of Warm Dust in S106: A Fragmented Environment

Authors: J. D. Adams, T. L. Herter, J. L. Hora, N. Schneider, R. M. Lau, J. G. Staughn, R. Simon, N. Smith, R. D. Gehrz, L. E. Allen, S. Bontemps, S. J. Carey, G. G. Fazio, R. A. Gutermuth, A. Guzman Fernandez, M. Hankins, T. Hill, E. Keto, X. P. Koenig, K. E. Kraemer, S. T. Megeath, D. R. Mizuno, F. Motte, P. C. Myers, H. A. Smith

Abstract: We present mid-IR (19 - 37 microns) imaging observations of S106 from SOFIA/FORCAST, complemented with IR observations from Spitzer/IRAC (3.6 - 8.0 microns), IRTF/MIRLIN (11.3 and 12.5 microns), and Herschel/PACS (70 and 160 microns). We use these observations, observations in the literature, and radiation transfer modeling to study the heating and composition of the warm (~ 100 K) dust in the reg… ▽ More We present mid-IR (19 - 37 microns) imaging observations of S106 from SOFIA/FORCAST, complemented with IR observations from Spitzer/IRAC (3.6 - 8.0 microns), IRTF/MIRLIN (11.3 and 12.5 microns), and Herschel/PACS (70 and 160 microns). We use these observations, observations in the literature, and radiation transfer modeling to study the heating and composition of the warm (~ 100 K) dust in the region. The dust is heated radiatively by the source S106 IR, with little contributions from grain-electron collisions and Ly-alpha radiation. The dust luminosity is >~ (9.02 +/- 1.01) x 10^4 L_sun, consistent with heating by a mid- to late-type O star. We find a temperature gradient (~ 75 - 107 K) in the lobes, which is consistent with a dusty equatorial geometry around S106 IR. Furthermore, the SOFIA observations resolve several cool (~ 65 - 70 K) lanes and pockets of warmer (~ 75 - 90 K) dust in the ionization shadow, indicating that the environment is fragmented. We model the dust mass as a composition of amorphous silicates, amorphous carbon, big grains, very small grains, and PAHs. We present the relative abundances of each grain component for several locations in S106. △ Less

Submitted 14 October, 2015; originally announced October 2015.

Comments: 18 pages, 11 figures, 4 tables, accepted by ApJ

A Keplerian-like disk around the forming O-type star AFGL 4176

Authors: Katharine G. Johnston, Thomas P. Robitaille, Henrik Beuther, Hendrik Linz, Paul Boley, Rolf Kuiper, Eric Keto, Melvin G. Hoare, Roy van Boekel

Abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) line and continuum observations at 1.2mm with ~0.3" resolution that uncover a Keplerian-like disk around the forming O-type star AFGL 4176. The continuum emission from the disk at 1.21 mm (source mm1) has a deconvolved size of 870+/-110 AUえーゆー x 330+/-300 AUえーゆー and arises from a structure ~8 M_sun in mass, calculated assuming a dust temperatu… ▽ More We present Atacama Large Millimeter/submillimeter Array (ALMA) line and continuum observations at 1.2mm with ~0.3" resolution that uncover a Keplerian-like disk around the forming O-type star AFGL 4176. The continuum emission from the disk at 1.21 mm (source mm1) has a deconvolved size of 870+/-110 AUえーゆー x 330+/-300 AUえーゆー and arises from a structure ~8 M_sun in mass, calculated assuming a dust temperature of 190 K. The first-moment maps, pixel-to-pixel line modeling, assuming local thermodynamic equilibrium (LTE), and position-velocity diagrams of the CH3CN J=13-12 K-line emission all show a velocity gradient along the major axis of the source, coupled with an increase in velocity at small radii, consistent with Keplerian-like rotation. The LTE line modeling shows that where CH3CN J=13-12 is excited, the temperatures in the disk range from ~70 to at least 300 K and that the H2 column density peaks at 2.8x10^24 cm^-2. In addition, we present Atacama Pathfinder Experiment (APEX) 12CO observations which show a large-scale outflow from AFGL 4176 perpendicular to the major axis of mm1, supporting the disk interpretation. Finally, we present a radiative transfer model of a Keplerian disk surrounding an O7 star, with a disk mass and radius of 12 M_sun and 2000 AUえーゆー, that reproduces the line and continuum data, further supporting our conclusion that our observations have uncovered a Keplerian disk around an O-type star. △ Less

Submitted 28 September, 2015; originally announced September 2015.

Comments: 6 pages, 4 figures, accepted for publication in ApJL

Combining radiative transfer and diffuse interstellar medium physics to model star formation

Authors: Matthew R. Bate, Eric R. Keto

Abstract: We present a method for modelling star-forming clouds that combines two different models of the thermal evolution of the interstellar medium (ISM). In the combined model, where the densities are low enough that at least some part of the spectrum is optically thin, a model of the thermodynamics of the diffuse ISM is more significant in setting the temperatures. Where the densities are high enough t… ▽ More We present a method for modelling star-forming clouds that combines two different models of the thermal evolution of the interstellar medium (ISM). In the combined model, where the densities are low enough that at least some part of the spectrum is optically thin, a model of the thermodynamics of the diffuse ISM is more significant in setting the temperatures. Where the densities are high enough to be optically thick across the spectrum, a model of flux limited diffusion is more appropriate. Previous methods either model the low-density interstellar medium and ignore the thermal behaviour at high densities (e.g. inside collapsing molecular cloud cores), or model the thermal behaviour near protostars but assume a fixed background temperature (e.g. approximately 10 K) on large-scales. Our new method treats both regimes. It also captures the different thermal evolution of the gas, dust, and radiation separately. We compare our results with those from the literature, and investigate the dependence of the thermal behaviour of the gas on the various model parameters. This new method should allow us to model the ISM across a wide range of densities and, thus, develop a more complete and consistent understanding of the role of thermodynamics in the star formation process. △ Less

Submitted 18 March, 2015; originally announced March 2015.

Comments: Accepted for publication in MNRAS, 25 pages, 22 figures

The dynamics of collapsing cores and star formation

Authors: Eric Keto, Paola Caselli, Jonathan Rawlings

Abstract: Low-mass stars are generally understood to form by the gravitational collapse of the dense molecular clouds known as starless cores. Continuum observations have not been able to distinguish among the several different hypotheses that describe the collapse because the predicted density distributions are the almost the same, as they are for all spherical self-gravitating clouds. However, the predict… ▽ More Low-mass stars are generally understood to form by the gravitational collapse of the dense molecular clouds known as starless cores. Continuum observations have not been able to distinguish among the several different hypotheses that describe the collapse because the predicted density distributions are the almost the same, as they are for all spherical self-gravitating clouds. However, the predicted contraction velocities are different enough that the models can be discriminated by comparing the velocities at large and small radii. This can be done by observing at least two different molecular line transitions that are excited at different densities. For example, the spectral lines of the H2O (110 - 101) and C18O (1-0) have critical densities for collisional de-excitation that differ by 5 orders of magnitude. We compare observations of these lines from the contracting starless core L1544 against the spectra predicted for several different hypothetical models of contraction including the Larson-Penston flow, the inside-out collapse of the singular isothermal sphere, the quasi-equilibrium contraction of an unstable Bonnor-Ebert sphere, and the non-equilibrium collapse of an over-dense Bonnor-Ebert sphere. Only the model of the unstable quasi-equilibrium Bonnor-Ebert sphere is able to produce the observed shapes of both spectral lines. This model allows us to interpret other observations of molecular lines in L1544 to find that the inward velocities seen in observations of CS(2-1) and N2H+ are located within the starless core itself, in particular in the region where the density profile follows an inverse square law. If these conclusions were to hold in the analysis of other starless cores, this would imply that the formation of hydrostatic clouds within the turbulent interstellar medium is not only possible but not exceptional and may be an evolutionary phase in low-mass star formation. △ Less

Submitted 21 October, 2014; originally announced October 2014.

Comments: 11 pages,7 figures

Magnetic Fields and Massive Star Formation

Authors: Qizhou Zhang, Keping Qiu, Josep M. Girart, Hauyu, Liu, Ya-Wen Tang, Patrick M. Koch, Zhi-Yun Li, Eric Keto, Paul T. P. Ho, Ramprasad Rao, Shih-Ping Lai, Tao-Chung Ching, Pau Frau, How-Huan Chen, Hua-Bai Li, Marco Padovani, Sylvain Bontemps, Timea Csengeri, Carmen Juarez

Abstract: Massive stars ($M > 8$ \msun) typically form in parsec-scale molecular clumps that collapse and fragment, leading to the birth of a cluster of stellar objects. We investigate the role of magnetic fields in this process through dust polarization at 870 $μみゅー$m obtained with the Submillimeter Array (SMA). The SMA observations reveal polarization at scales of $\lsim$ 0.1 pc. The polarization pattern in… ▽ More Massive stars ($M > 8$ \msun) typically form in parsec-scale molecular clumps that collapse and fragment, leading to the birth of a cluster of stellar objects. We investigate the role of magnetic fields in this process through dust polarization at 870 $μみゅー$m obtained with the Submillimeter Array (SMA). The SMA observations reveal polarization at scales of $\lsim$ 0.1 pc. The polarization pattern in these objects ranges from ordered hour-glass configurations to more chaotic distributions. By comparing the SMA data with the single dish data at parsec scales, we found that magnetic fields at dense core scales are either aligned within $40^\circ$ of or perpendicular to the parsec-scale magnetic fields. This finding indicates that magnetic fields play an important role during the collapse and fragmentation of massive molecular clumps and the formation of dense cores. We further compare magnetic fields in dense cores with the major axis of molecular outflows. Despite a limited number of outflows, we found that the outflow axis appears to be randomly oriented with respect to the magnetic field in the core. This result suggests that at the scale of accretion disks ($\lsim 10^3$ AUえーゆー), angular momentum and dynamic interactions possibly due to close binary or multiple systems dominate over magnetic fields. With this unprecedentedly large sample massive clumps, we argue on a statistical basis that magnetic fields play an important role during the formation of dense cores at spatial scale of 0.01 - 0.1 pc in the context of massive star and cluster star formation. △ Less

Submitted 15 July, 2014; originally announced July 2014.

Comments: Accepted for publication in Astrophysical Journal

Greenland Telescope Project --- Direct Confirmation of Black Hole with Sub-millimeter VLBI

Authors: M. Inoue, J. C. Algaba-Marcos, K. Asada, C. -C. Chang, M. -T. Chen, J. Han, H. Hirashita, P. T. P. Ho, S. -N. Hsieh, T. Huang, H. Jiang, P. M. Koch, D. Y. Kubo, C. -Y. Kuo, B. Liu, P. Martin-Cocher, S. Matsushita, Z. Meyer-Zhao, M. Nakamura, H. Nishioka, G. Nystrom, N. Pradel, H. -Y. Pu, P. A. Raffin, H. -Y. Shen , et al. (14 additional authors not shown)

Abstract: A 12-m diameter radio telescope will be deployed to the Summit Station in Greenland to provide direct confirmation of a Super Massive Black Hole (SMBH) by observing its shadow image in the active galaxy M87. The telescope (Greenland Telescope: GLT) is to become one of the Very Long Baseline Interferometry (VLBI) stations at sub-millimeter (submm) regime, providing the longest baseline > 9,000 km t… ▽ More A 12-m diameter radio telescope will be deployed to the Summit Station in Greenland to provide direct confirmation of a Super Massive Black Hole (SMBH) by observing its shadow image in the active galaxy M87. The telescope (Greenland Telescope: GLT) is to become one of the Very Long Baseline Interferometry (VLBI) stations at sub-millimeter (submm) regime, providing the longest baseline > 9,000 km to achieve an exceptional angular resolution of 20 micro arc sec at 350 GHz, which will enable us to resolve the shadow size of ~40 micro arc sec. The triangle with the longest baselines formed by the GLT, the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile, and the Submillimeter Array (SMA) in Hawaii will play a key role for the M87 observations. We have been working on the image simulations based on realistic conditions for a better understanding of the possible observed images. In parallel, retrofitting of the telescope and the site developments are in progress. Based on three years of opacity monitoring at 225 GHz, our measurements indicate that the site is excellent for submm observations, comparable to the ALMA site. The GLT is also expected to make single-dish observations up to 1.5 THz. △ Less

Submitted 9 July, 2014; originally announced July 2014.

Comments: 15 pages, 6 figures, accepted by Radio Science journal

Core and filament formation in magnetized, self-gravitating isothermal layers

Authors: Sven Van Loo, Eric Keto, Qizhou Zhang

Abstract: We examine the role of the gravitational instability in an isothermal, self-gravitating layer threaded by magnetic fields on the formation of filaments and dense cores. Using numerical simulation we follow the non-linear evolution of a perturbed equilibrium layer. The linear evolution of such a layer is described in the analytic work of Nagai et al (1998). We find that filaments and dense cores fo… ▽ More We examine the role of the gravitational instability in an isothermal, self-gravitating layer threaded by magnetic fields on the formation of filaments and dense cores. Using numerical simulation we follow the non-linear evolution of a perturbed equilibrium layer. The linear evolution of such a layer is described in the analytic work of Nagai et al (1998). We find that filaments and dense cores form simultaneously. Depending on the initial magnetic field, the resulting filaments form either a spiderweb-like network (for weak magnetic fields) or a network of parallel filaments aligned perpendicular to the magnetic field lines (for strong magnetic fields). Although the filaments are radially collapsing, the density profile of their central region (up to the thermal scale height) can be approximated by a hydrodynamical equilibrium density structure. Thus, the magnetic field does not play a significant role in setting the density distribution of the filaments. The density distribution outside of the central region deviates from the equilibrium. The radial column density distribution is then flatter than the expected power law of $r^{-4}$ and similar to filament profiles observed with Herschel. Our results does not explain the near constant filament width of $\sim 0.1$pc. However, our model does not include turbulent motions. It is expected that accretion-driven amplification of these turbulent motions provides additional support within the filaments against gravitational collapse. Finally, we interpret the filamentary network of the massive star forming complex G14.225-0.506 in terms of the gravitational instability model and find that the properties of the complex are consistent with being formed out of an unstable layer threaded by a strong, parallel magnetic field. △ Less

Submitted 5 May, 2014; originally announced May 2014.

Comments: Accepted by ApJ, 14 pages, 22 figures

Chemistry and Radiative Transfer of Water in Cold, Dense Clouds

Authors: Eric Keto, Jonathan Rawlings, Paola Caselli

Abstract: The Herschel Space Observatory's recent detections of water vapor in the cold, dense cloud L1544 allow a direct comparison between observations and chemical models for oxygen species in conditions just before star formation. We explain a chemical model for gas phase water, simplified for the limited number of reactions or processes that are active in extreme cold ($<$ 15 K). In this model, water i… ▽ More The Herschel Space Observatory's recent detections of water vapor in the cold, dense cloud L1544 allow a direct comparison between observations and chemical models for oxygen species in conditions just before star formation. We explain a chemical model for gas phase water, simplified for the limited number of reactions or processes that are active in extreme cold ($<$ 15 K). In this model, water is removed from the gas phase by freezing onto grains and by photodissociation. Water is formed as ice on the surface of dust grains from O and OH and released into the gas phase by photodesorption. The reactions are fast enough with respect to the slow dynamical evolution of L1544 that the gas phase water is in equilibrium for the local conditions thoughout the cloud. We explain the paradoxical radiative transfer of the H$_2$O ($1_{10}-1_{01}$) line. Despite discouragingly high optical depth caused by the large Einstein A coefficient, the subcritical excitation in the cold, rarefied H$_2$ causes the line brightness to scale linearly with column density. Thus the water line can provide information on the chemical and dynamical processes in the darkest region in the center of a cold, dense cloud. The inverse P-Cygni profile of the observed water line generally indicates a contracting cloud. This profile is reproduced with a dynamical model of slow contraction from unstable quasi-static hydrodynamic equilibrium (an unstable Bonnor-Ebert sphere). △ Less

Submitted 1 March, 2014; originally announced March 2014.

Comments: submitted to MNRAS

From Filaments to Oscillating Starless Cores

Authors: Eric Keto, Andreas Burkert

Abstract: Long wavelength sonic oscillations are observed or inferred in many of the small, dark molecular clouds, the starless cores, that are the precursors to protostars. The oscillations provide significant internal energy and the time scale for their dissipation may control the rate of star formation in the starless cores. Despite their importance, their origin is unknown. We explore one hypothesis tha… ▽ More Long wavelength sonic oscillations are observed or inferred in many of the small, dark molecular clouds, the starless cores, that are the precursors to protostars. The oscillations provide significant internal energy and the time scale for their dissipation may control the rate of star formation in the starless cores. Despite their importance, their origin is unknown. We explore one hypothesis that the oscillations develop as a starless core forms from a filament. We model this process with a numerical hydrodynamic simulation and generate synthetic molecular line observations with a radiative transfer simulation. Comparison with actual observations shows general agreement suggesting the proposed mechanism is viable. △ Less

Submitted 14 February, 2014; originally announced February 2014.

Comments: submitted to MNRAS

Flickering of 1.3 cm Sources in Sgr B2: Towards a Solution to the Ultracompact HII Region Lifetime Problem

Authors: C. G. De Pree, T. Peters, M. -M. Mac Low, D. J. Wilner, W. M. Goss, R. Galván-Madrid, E. R. Keto, R. S. Klessen, A. Monsrud

Abstract: Accretion flows onto massive stars must transfer mass so quickly that they are themselves gravitationally unstable, forming dense clumps and filaments. These density perturbations interact with young massive stars, emitting ionizing radiation, alternately exposing and confining their HII regions. As a result, the HII regions are predicted to flicker in flux density over periods of decades to centu… ▽ More Accretion flows onto massive stars must transfer mass so quickly that they are themselves gravitationally unstable, forming dense clumps and filaments. These density perturbations interact with young massive stars, emitting ionizing radiation, alternately exposing and confining their HII regions. As a result, the HII regions are predicted to flicker in flux density over periods of decades to centuries rather than increasing monotonically in size as predicted by simple Spitzer solutions. We have recently observed the Sgr B2 region at 1.3 cm with the VLA in its three hybrid configurations (DnC, CnB and BnA) at a resolution of 0.25''. These observations were made to compare in detail with matched continuum observations from 1989. At 0.25'' resolution, Sgr B2 contains 41 UC HII regions, 6 of which are hypercompact. The new observations of Sgr B2 allow comparison of relative peak flux densites for the HII regions in Sgr B2 over a 23 year time baseline (1989-2012) in one of the most source-rich massive star forming regions in the Milky Way. The new 1.3 cm continuum images indicate that four of the 41 UC HII regions exhibit significant changes in their peak flux density, with one source (K3) dropping in peak flux density, and the other 3 sources (F10.303, F1 and F3) increasing in peak flux density. The results are consistent with statistical predictions from simulations of high mass star formation, suggesting that they offer a solution to the lifetime problem for ultracompact HII regions. △ Less

Submitted 30 December, 2013; originally announced December 2013.

Comments: 12 pages, 3 figures, Accepted for publication in the Astrophysical Journal Letters

The Dynamics of Ultracompact HII Regions

Authors: Nathaniel Roth, Steven W. Stahler, Eric Keto

Abstract: Many ultracompact HII regions exhibit a cometary morphology in radio continuum emission. In such regions, a young massive star is probably ablating, through its ultraviolet radiation, the molecular cloud clump that spawned it. On one side of the star, the radiation drives an ionization front that stalls in dense molecular gas. On the other side, ionized gas streams outward into the more rarefied e… ▽ More Many ultracompact HII regions exhibit a cometary morphology in radio continuum emission. In such regions, a young massive star is probably ablating, through its ultraviolet radiation, the molecular cloud clump that spawned it. On one side of the star, the radiation drives an ionization front that stalls in dense molecular gas. On the other side, ionized gas streams outward into the more rarefied environment. This wind is underpressured with respect to the neutral gas. The difference in pressure draws in more cloud material, feeding the wind until the densest molecular gas is dissipated. Recent, time-dependent simulations of massive stars turning on within molecular gas show the system evolving in a direction similar to that just described. Here, we explore a semi-analytic model in which the wind is axisymmetric and has already achieved a steady state. Adoption of this simplified picture allows us to study the dependence of both the wind and its bounding ionization front on the stellar luminosity, the peak molecular density, and the displacement of the star from the center of the clump. For typical parameter values, the wind accelerates transonically to a speed of about 15 km / s, and transports mass outward at a rate of 10^-4 solar masses per year. Stellar radiation pressure acts to steepen the density gradient of the wind. △ Less

Submitted 22 November, 2013; originally announced November 2013.

Comments: Preprint format, including 21 figures at the end

MUSCLE W49 : A Multi-Scale Continuum and Line Exploration of the Most Luminous Star Formation Region in the Milky Way. I. Data and The Mass Structure of the Giant Molecular Cloud

Authors: R. Galván-Madrid, H. B. Liu, Z. -Y. Zhang, J. E. Pineda, T. -C. Peng, Q. Zhang, E. R. Keto, P. T. P. Ho, L. F. Rodríguez, L. Zapata, T. Peters, C. G. De Pree

Abstract: The Multi-Scale Continuum and Line Exploration of W49 (MUSCLE W49) is a comprehensive gas and dust survey of the giant molecular cloud (GMC) of W49A, the most luminous star-formation region in the Milky Way. It covers the entire GMC at different scales and angular resolutions. In this paper we present: 1) an all-configuration SMA mosaic in the 230-GHz band covering the central 3 arcmin (10 pc, kno… ▽ More The Multi-Scale Continuum and Line Exploration of W49 (MUSCLE W49) is a comprehensive gas and dust survey of the giant molecular cloud (GMC) of W49A, the most luminous star-formation region in the Milky Way. It covers the entire GMC at different scales and angular resolutions. In this paper we present: 1) an all-configuration SMA mosaic in the 230-GHz band covering the central 3 arcmin (10 pc, known as W49N), with most of the embedded massive stars; and 2) PMO 14m telescope observations in the 90-GHz band, covering the entire GMC with maps up to 35 arcmin in size, or 113 pc. We also make use of archival data from the VLA, JCMT-SCUBA, IRAM 30m, and the CSO BOLOCAM GPS. Our main findings are: 1) The W49 GMC is one of the most massive in the Galaxy, with a total mass ~1.1x10^6 Msun within a radius of 60 pc. Within a radius of 6 pc, the total gas mass is ~2x10^5 Msun. At these scales only 1% of the material is photoionized. The mass reservoir is sufficient to form several young massive clusters (YMCs) as massive as a globular cluster. 2) The mass of the GMC is distributed in a hierarchical network of filaments. At scales <10 pc, a triple, centrally condensed structure peaks toward the ring of HC HII regions in W49N. This structure extends to scales from ~10 to 100 pc. The W49A starburst most likely formed from global gravitational contraction with localized collapse in a "hub-filament" geometry. 3) Currently, feedback from the central YMCs (with a present mass Mcl > 5x10^4 Msun) is still not enough to entirely disrupt the GMC, but further stellar mass growth could be enough to allow radiation pressure to clear the cloud and halt star formation. 4) The resulting stellar content will probably remain as a gravitationally bound massive star cluster, or a small system of bound clusters. (ABRIDGED) △ Less

Submitted 9 January, 2014; v1 submitted 16 September, 2013; originally announced September 2013.

Comments: Author version of paper published in ApJ. For high-quality figures please see the Journal version

Journal ref: 2013, ApJ, 779, 121