Search

(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., This paper is published in the Journal of Astronomical Instrumentation, 12(1), 2340006 (2023) https://doi.org/10.1142/S2251171723400068

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.

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.

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., Comment: 44 pages, 41 figures

Pre-stellar cores represent the initial conditions in the process of star and planet formation. Their low temperatures ($, Comment: accepted for publication in The Astrophysical Journal

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., Comment: 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

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., 41 pages, 25 figures, accepted for publication in ApJ

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., Comment: Accepted MNRAS; 8 pages, 5 figures

Aims. The main goal of the following analysis is to assess the potential role of magnetic fields in regulating the envelope rotation, the formation of disks and the fragmentation of Class 0 protostars in multiple systems. Methods. We use the Submillimeter Array to carry out observations of the dust polarized emission at 0.87 mm, 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 to characterize the orientation of the main component of the organized magnetic field at the envelope scales in these embedded protostars. This direction is compared to that of the protostellar outflow in order to study the relation between their misalignment and the kinematics of the circumstellar gas. The latter is traced via velocity gradient observed in the molecular line emission (mainly N2H+) of the gas at intermediate envelope scales. Results. 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., This project has received funding from the European Research Council (ERC) under the European Union Horizon 2020 research and innovation programme (MagneticYSOs project, grant agreement No. 679937, PI: Maury). J.M.G. is supported by the grant AYA2017-84390- C2-R (AEI/FEDER, UE).

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., 34 pages, 27 figures

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., Comment: Accepted for Publication in ApJS

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., To be published in MNRAS

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., Comment: To be published in MNRAS

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., Comment: 13 pages, 13 figures, 1 appendix, accepted for publication in Montly Notices of the Royal Astronomical Society

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., Comment: accepted for publication in The Astrophysical Journal

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., Accepted for publication in MNRAS

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., 15 pages, 15 figures, accepted for publication in MNRAS

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., Comment: 12 pages, 7 figures. Accepted for publication in MNRAS

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., 22 pages, 13 figures, 3 tables. Accepted for publication in MNRAS

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., 19 pages, 10 figures, 3 tables. Accepted for publication in MNRAS

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., Comment: 23 pages, 26 figures, accepted for publication in Astronomy and Astrophysics

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., Accepted to ApJ, 35 pages, 15 figures, 4 Tables. Intensity-corrected FITS images and temperature/optical depths maps are available with the electronic version

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., Comment: 22 pages, 9 figures, accepted for publication in A&A

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.

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., submitted to MNRAS

The newly discovered Massive Molecular Filament (MMF) G32.02+0.05 (~ 70 pc long, 105 M⊙) has been shaped and compressed by older generations of massive stars. The similarity of this filament in physical structure (density profile, temperature) to much smaller star-forming filaments, suggests that the mechanism to form such filaments may be a universal process. The densest portion of the filament, apparent as an Infrared Dark Cloud (IRDC) shows a range of massive star formation signatures throughout. We investigate the kinematics in this filament and find widespread inverse P cygni asymmetric line profiles. These line asymmetries are interpreted as a signature of large-scale radial collapse. Using line asymmetries observed with optically thick HCO+ (1-0) and optically thin H13CO+ (1-0) across a range of massive star forming regions in the filament, we estimate the global radial infall rate of the filament to range from a few 100 to a few 1000 M⊙ Myr−1 pc−1. At its current infall rate the densest portions of the cloud will more than double their current mass within a Myr.

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\deg$ (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., Comment: Paper accepted for publications in A&A

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., 24 pages, 6 figures, 3 tables, accepted for publication in ApJ

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., 39 pages, 17 figures, accepted for publication in PASJ (as a review paper)

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., 6 pages, 4 figures, accepted for publication in ApJL

Ultracompact and hypercompact H II regions appear when a star with a mass larger than about 15 Mstarts to ionize its own environment. Recent observations of time variability in these objects are one of the pieces of evidence that suggest that at least some of them harbour stars that are still accreting from an infalling neutral accretion flow that becomes ionized in its innermost part. We present an analysis of the properties of the H II regions formed in the three-dimensional radiation-hydrodynamic simulations presented by Peters et al. as a function of time. Flickering of the H II regions is a natural outcome of this model. The radio-continuum fluxes of the simulated H II regions as well as their flux and size variations are in agreement with the available observations. From the simulations, we estimate that a small but non-negligible fraction (∼10 per cent) of observed H II regions should have detectable flux variations (larger than 10 per cent) on time-scales of ∼10 yr, with positive variations being more likely to happen than negative variations. A novel result of these simulations is that negative flux changes do happen, in contrast to the simple expectation of ever growing H II regions. We also explore the temporal correlations between properties that are directly observed (flux and size) and other quantities like density and ionization rates.

This is the first paper of our series of high resolution (1") studies of the massive star forming region G10.6--0.4. We present the emission line observations of the hot core type tracers (O$^{13}$CS, OCS, SO$_{2}$) with $\sim$0$"$.5 resolution. By comparing the results to the high--resolution NH$_{3}$ absorption line observation, we confirm for the first time the rotationally flattened hot toroid in the central $, Comment: 26 pages, 5 figures, Accepted by ApJ

In this paper we discuss two approximate methods previously suggested for modeling hyperfine spectral line emission for molecules whose collisional transitions rates between hyperfine levels are unknown. Hyperfine structure is seen in the rotational spectra of many commonly observed molecules such as HCN, HNC, NH3, N2H+, and C17O. The intensities of these spectral lines can be modeled by numerical techniques such as Lambda-iteration that alternately solve the equations of statistical equilibrium and the equation of radiative transfer. However, these calculations require knowledge of both the radiative and collisional rates for all transitions. For most commonly observed radio frequency spectral lines, only the net collisional rates between rotational levels are known. For such cases, two approximate methods have been suggested. The first method, hyperfine statistical equilibrium (HSE), distributes the hyperfine level populations according to their statistical weight, but allows the population of the rotational states to depart from local thermodynamic equilibrium (LTE). The second method, the proportional method approximates the collision rates between the hyperfine levels as fractions of the net rotational rate apportioned according to the statistical degeneracy of the final hyperfine levels. The second method is able to model non-LTE hyperfine emission. We compare simulations of N2H+ hyperfine lines made with approximate and more exact rates and find that satisfactory results are obtained., Comment: 34 pages. Pages 22-34 are data tables. For ApJ.

We address the question of whether the formation of high-mass stars is similar to or differs from that of solar-mass stars through new molecular line observations and modeling of the accretion flow around the massive protostar IRAS20126+4104. We combine new observations of NH3(1,1) and (2,2) made at the Very Large Array, new observations of CHCN(13-12) made at the Submillimeter Array, previous VLA observations of NH(3,3), NH(4,4), and previous Plateau de Bure observations of C34S(2-1), C34S(5-4), and CHCN(12-11) to obtain a data set of molecular lines covering 15 to 419 K in excitation energy. We compare these observations against simulated molecular line spectra predicted from a model for high-mass star formation based on a scaled-up version of the standard disk-envelope paradigm developed for accretion flows around low-mass stars. We find that in accord with the standard paradigm, the observations require both a warm, dense, rapidly-rotating disk and a cold, diffuse infalling envelope. This study suggests that accretion processes around 10 M stars are similar to those of solar mass stars.

An example of a cold massive core, JCMT 18354-0649S, a possible high mass analogue to a low mass star forming core is studied. Line and continuum observations from JCMT, Mopra Telescope and Spitzer are presented and modelled in detail using a 3D molecular line radiative transfer code. In almost every way JCMT 18354-0649S is a scaled-up version of a typical low mass core with similar temperatures, chemical abundances and densities. The difference is that both the infall velocity and the turbulent width of the line profiles are an order of magnitude larger. While the higher infall velocity is expected due to the large mass of JCMT 18354-0649S, we suggest that the dissipation of this highly supersonic turbulence may lead to the creation of dense clumps of gas that surround the high mass core., 12 Pages, 17 Figures, 2 Tables

We present Australia Telescope Compact Array observations towards 6 massive star formation regions which, from their strong 24 GHz continuum emission but no compact 8 GHz continuum emission, appeared good candidates for hyper-compact HII regions. However, the properties of the ionised gas derived from the 19 to 93 GHz continuum emission and H70 alpha + H57 alpha radio recombination line data show the majority of these sources are, in fact, regions of spatially-extended, optically-thin free-free emission. These extended sources were missed in the previous 8 GHz observations due to a combination of spatial-filtering, poor surface brightness sensitivity and primary beam attenuation. We consider the implications that a significant number of these extended HII regions may have been missed by previous surveys of massive star formation regions. If the original sample of 21 sources is representative of the population as a whole, the fact that 6 contain previously undetected extended free-free emission suggests a large number of regions have been mis-classified. Rather than being very young objects prior to UCHII region formation, they are, in fact, associated with extended HII regions and thus significantly older. In addition, inadvertently ignoring a potentially substantial flux contribution (up to ~0.5Jy) from free-free emission has implications for dust masses derived from sub-mm flux densities. The large spatial scales probed by single-dish telescopes, which do not suffer from spatial filtering, are particularly susceptible and dust masses may be overestimated by up to a factor of ~2., Comment: Accepted MNRAS. 22 pages, 13 figures, 5 tables

In a previous paper we demonstrated that nonradial hydrodynamic oscillations of a thermally supported (Bonnor-Ebert) sphere embedded in a low-density, high-temperature medium persist for many periods. The predicted column density variations and molecular spectral line profiles are similar to those observed in the Bok globule B68, suggesting that the motions in some starless cores may be oscillating perturbations on a thermally supported equilibrium structure. Such oscillations can produce molecular line maps which mimic rotation, collapse, or expansion and, thus, could make determining the dynamical state from such observations alone difficult. However, while B68 is embedded in a very hot, low-density medium, many starless cores are not, having interior/exterior density contrasts closer to unity. In this paper we investigate the oscillation damping rate as a function of the exterior density. For concreteness we use the same interior model employed by Broderick et al., with varying models for the exterior gas. We also develop a simple analytical formalism, based on the linear perturbation analysis of the oscillations, which predicts the contribution to the damping rates due to the excitation of sound waves in the external medium. We find that the damping rate of oscillations on globules in dense molecular environments is always many periods, corresponding to hundreds of thousands of years and persisting over the inferred lifetimes of the globules.

We present a model in which the supersonic motions observed in molecular clouds are driven by gravitational energy released as large structures fragment into smaller ones. The fragmentation process begins in large molecular clouds, and continues down to fragments of a critical mass defined as the mass at which gravitational confinement may be replaced by pressure confinement. The power laws that describe the scaling of density, mass, and number spectra of the fragments are given in terms of the observed velocity dispersion of the fragments. The results agree with observations over the range from several to about a third of a million solar masses.

In a previous paper, we modeled the oscillations of a thermally-supported (Bonnor-Ebert) sphere as non-radial, linear perturbations following a standard analysis developed for stellar pulsations. The predicted column density variations and molecular spectral line profiles are similar to those observed in the Bok globule B68 suggesting that the motions in some starless cores may be oscillating perturbations on a thermally supported equilibrium structure. However, the linear analysis is unable to address several questions, among them the stability, and lifetime of the perturbations. In this paper we simulate the oscillations using a three-dimensional numerical hydrodynamic code. We find that the oscillations are damped predominantly by non-linear mode-coupling, and the damping time scale is typically many oscillation periods, corresponding to a few million years, and persisting over the inferred lifetime of gobules., 7 pages, 7 figures, accepted by ApJ

The fine-structure lines of singly ([Ne II] 12.8 micron) and doubly ([Ne III] 15.6 micron) ionized neon are among the most prominent features in the mid-infrared spectra of star-forming regions, and have the potential to be a powerful new indicator of the star formation rate in galaxies. Using a sample of star-forming galaxies with measurements of the fine-structure lines available from the literature, we show that the sum of the [Ne II] and [Ne III] luminosities obeys a tight, linear correlation with the total infrared luminosity, over 5 orders of magnitude in luminosity. We discuss the formation of the lines and their relation with the Lyman continuum luminosity. A simple calibration between star formation rate and the [Ne II]+[Ne III] luminosity is presented., To appear in ApJ. 8 pages

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., Accepted for publication in MNRAS, 25 pages, 22 figures

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$\sigma$) 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$\sigma$) 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., Comment: 20 pages, 5 figures, Accepted for publication in The Astrophysical Journal

We test the hypothesis that the starless cores may be gravitationally bound clouds supported largely by thermal pressure by comparing observed molecular line spectra to theoretical spectra produced by a simulation that includes hydrodynamics, radiative cooling, variable molecular abundance, and radiative transfer in a simple one-dimensional model. The results suggest that the starless cores can be divided into two categories: stable starless cores that are in approximate equilibrium and will not evolve to form protostars, and unstable pre-stellar cores that are proceeding toward gravitational collapse and the formation of protostars. The starless cores might be formed from the interstellar medium as objects at the lower end of the inertial cascade of interstellar turbulence. Additionally, we identify a thermal instability in the starless cores. Under par ticular conditions of density and mass, a core may be unstable to expansion if the density is just above the critical density for the collisional coupling of the gas and dust so that as the core expands the gas-dust coupling that cools the gas is reduced and the gas warms, further driving the expansion., Submitted to ApJ

We observed the nearby starburst galaxy M82 in CO in the higher frequency (2--1) transition to achieve an angular resolution below 1 arc second or 17 pc at the target. We resolved the molecular gas into a large number of compact clouds, with masses ranging from about 2x10^3 to 2x10^6 solar masses. The mass spectrum scales as N(M) ~ M^-1.5, similar to the mass spectra of young massive star clusters suggesting that individual molecular clouds are transformed in the starburst into individual star clusters. The larger clouds are surrounded by supernovae and HII regions suggesting that star formation proceeds from the outside of the clouds and progresses inward consistent with triggering by a sudden increase in external pressure. The clouds with internal star formation have velocity gradients and inverse P-Cygni spectral line profiles indicating inward motions of 35 kms consistent with shock driven compression. Diffuse free-free radio emission and X-ray emission around the clouds provides evidence for superheated ionized gas sufficient to drive the compression. Clouds with spectral lines indicating expansion show little internal star formation suggesting that the dynamics precedes and is responsible for the star formation rather than the inverse. M82 is known to be in interaction with neighboring M81. The overall picture is consistent with the formation of massive star clusters from individual giant molecular clouds crushed by a sudden galactic scale increase in external pressure generated by the changing dynamics that result from a near-collision with a neighboring galaxy. Present day globular clusters may have formed in a similar fashion in primordial galaxies., Submitted to ApJ

We present the first 865 μm continuum image with subarcsecond resolution obtained with the Submillimeter Array. These data resolve the Orion KL region into the hot core, the nearby radio source I, the submillimeter counterpart to the infrared source n (radio source L), and new submillimeter continuum sources. The radio to submillimeter emission from source I can be modeled as either the result of proton-electron free-free emission that is optically thick to ~100 GHz plus dust emission that accounts for the majority of the submillimeter flux, or H- free-free emission that gives rise to a power-law spectrum with a power-law index of ~1.6. The latter model would indicate similar physical conditions as found in the inner circumstellar environment of Mira variable stars. Future subarcsecond resolution observations at shorter submillimeter wavelengths should easily discriminate between these two possibilities. The submillimeter continuum emission toward source n can be interpreted in the framework of emission from an accretion disk.

We develop a method of analyzing radio frequency spectral line observations to derive data on the temperature, density, velocity, and molecular abundance of the emitting gas. The method incorporates a radiative transfer code with a new technique for handling overlapping hyperfine emission lines within the accelerated lambda iteration algorithm and a heuristic search algorithm based on simulated annnealing. We apply this method to new observations of N_2H^+ in three Lynds clouds thought to be starless cores in the first stages of star formation and determine their density structure. A comparison of the gas densities derived from the molecular line emission and the millimeter dust emission suggests that the required dust mass opacity is about kappa_{1.3mm}=0.04 cm^2/g, consistent with models of dust grains that have opacities enhanced by ice mantles and fluffy aggregrates., Comment: 42 pages, 17 figures, to appear in ApJ