• The Astropy Collaboration, A. M. Price-Whelan, B. M. Sipőcz, H. M. Günther, P. L. Lim, S. M. Crawford, S. Conseil, D. L. Shupe, M. W. Craig, N. Dencheva, A. Ginsburg, J. T. VanderPlas, L. D. Bradley, D. Pérez-Suárez, M. de Val-Borro, T. L. Aldcroft, K. L. Cruz, T. P. Robitaille, E. J. Tollerud, C. Ardelean, T. Babej, M. Bachetti, A. V. Bakanov, S. P. Bamford, G. Barentsen, P. Barmby, A. Baumbach, K. L. Berry, F. Biscani, M. Boquien, K. A. Bostroem, L. G. Bouma, G. B. Brammer, E. M. Bray, H. Breytenbach, H. Buddelmeijer, D. J. Burke, G. Calderone, J. L. Cano Rodríguez, M. Cara, J. V. M. Cardoso, S. Cheedella, Y. Copin, D. Crichton, D. DÁvella, C. Deil, É. Depagne, J. P. Dietrich, A. Donath, M. Droettboom, N. Earl, T. Erben, S. Fabbro, L. A. Ferreira, T. Finethy, R. T. Fox, L. H. Garrison, S. L. J. Gibbons, D. A. Goldstein, R. Gommers, J. P. Greco, P. Greenfield, A. M. Groener, F. Grollier, A. Hagen, P. Hirst, D. Homeier, A. J. Horton, G. Hosseinzadeh, L. Hu, J. S. Hunkeler, Ž. Ivezić, A. Jain, T. Jenness, G. Kanarek, S. Kendrew, N. S. Kern, W. E. Kerzendorf, A. Khvalko, J. King, D. Kirkby, A. M. Kulkarni, A. Kumar, A. Lee, D. Lenz, S. P. Littlefair, Z. Ma, D. M. Macleod, M. Mastropietro, C. McCully, S. Montagnac, B. M. Morris, M. Mueller, S. J. Mumford, D. Muna, N. A. Murphy, S. Nelson, G. H. Nguyen, J. P. Ninan, M. Nöthe, S. Ogaz, S. Oh, J. K. Parejko, N. Parley, S. Pascual, R. Patil, A. A. Patil, A. L. Plunkett, J. X. Prochaska, T. Rastogi, V. Reddy Janga, J. Sabater, P. Sakurikar, M. Seifert, L. E. Sherbert, H. Sherwood-Taylor, A. Y. Shih, J. Sick, M. T. Silbiger, S. Singanamalla, L. P. Singer, P. H. Sladen, K. A. Sooley, S. Sornarajah, O. Streicher, P. Teuben, S. W. Thomas, G. R. Tremblay, J. E. H. Turner, V. Terrón, M. H. van Kerkwijk, A. de la Vega, L. L. Watkins, B. A. Weaver, J. B. Whitmore, J. Woillez, V. Zabalza
    Jan. 16, 2018 astro-ph.IM
    The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.
  • Constraining the temperature and density structure of dense molecular cloud cores is fundamental for understanding the initial conditions of star formation. We use Herschel observations of the thermal FIR dust emission from nearby isolated molecular cloud cores and combine them with ground-based submillimeter continuum data to derive observational constraints on their temperature and density structure. The aim of this study is to verify the validity of a ray-tracing inversion technique developed to derive the dust temperature and density structure of isolated starless cores directly from the dust emission maps and to test if the resulting temperature and density profiles are consistent with physical models. Using this ray-tracing inversion technique, we derive the dust temperature and density structure of six isolated starless cloud cores. We employ self-consistent radiative transfer modeling to the derived density profiles, treating the ISRF as the only heating source. The best-fit values of local strength of the ISRF and the extinction by the outer envelope are derived by comparing the self-consistently calculated temperature profiles with those derived by the ray-tracing method. We find that all starless cores are significantly colder inside than outside, with the core temperatures showing a strong negative correlation with peak column density. This suggests that their thermal structure is dominated by external heating from the ISRF and shielding by dusty envelopes. The temperature profiles derived with the ray-tracing inversion method can be well-reproduced with self-consistent radiative transfer models.
  • We have updated our publicly available dust radiative transfer code (HOCHUNK3D) to include new emission processes and various 3-D geometries appropriate for forming stars. The 3-D geometries include warps and spirals in disks, accretion hotspots on the central star, fractal clumping density enhancements, and misaligned inner disks. Additional axisymmetric (2-D) features include gaps in disks and envelopes, "puffed-up inner rims" in disks, multiple bipolar cavity walls, and iteration of disk vertical structure assuming hydrostatic equilibrium. We include the option for simple power-law envelope geometry, which combined with fractal clumping, and bipolar cavities, can be used to model evolved stars as well as protostars. We include non-thermal emission from PAHs and very small grains, and external illumination from the interstellar radiation field. The grid structure was modified to allow multiple dust species in each cell; based on this, a simple prescription is implemented to model dust stratification. We describe these features in detail, and show example calculations of each. Some of the more interesting results include the following: 1) Outflow cavities may be more clumpy than infalling envelopes. 2) PAH emission in high-mass stars may be a better indicator of evolutionary stage than the broadband SED slope; and related to this, 3) externally illuminated clumps and high-mass stars in optically thin clouds can masquerade as YSOs. 4) Our hydrostatic equilibrium models suggest that dust settling is likely ubiquitous in T Tauri disks, in agreement with previous observations.
  • This paper aims to investigate the hypothesis that the embedded luminous star AFGL2591-VLA3 (2.3E5Lsun at 3.33kpc) is forming according to a scaled-up version of a low-mass star formation scenario. We present multi-configuration VLA 3.6cm and 7mm, as well as CARMA C18O and 3mm continuum observations to investigate the morphology and kinematics of the ionized gas, dust, and molecular gas around AFGL2591. We also compare our results to ancillary near-IR images, and model the SED and 2MASS image profiles of AFGL2591 using a dust continuum radiative transfer code. The observed 3.6cm images uncover for the first time that the central powering source AFGL2591-VLA3 has a compact core plus collimated jet morphology, extending 4000AU eastward from the central source with an opening angle of <10deg at this radius. However, at 7mm VLA3 does not show a jet morphology, but instead compact (<500AU) emission, some of which (<0.57 mJy of 2.9mJy) is estimated to be from dust. We determine that the momentum rate of the jet is not sufficient to ionize itself via only shocks, and thus a significant portion of the emission is instead likely created in a photoionized wind. The C18O emission uncovers dense entrained material in the outflow(s) from the young stars in the region. The main features of the SED and 2MASS images of AFGL2591-VLA3 are also reproduced by our model dust geometry of a rotationally flattened envelope with and without a disk. The above results are consistent with a picture of massive star formation similar to that seen for low-mass protostars. However, within its envelope, AFGL2591-VLA3 contains at least four other young stars, constituting a small cluster. Therefore it appears that AFGL2591-VLA3 may be able to source its accreting material from a shared gas reservoir while still exhibiting the phenomena expected during the formation of low-mass stars. (Abridged)
  • We present the results of a Nobeyama 45-m water maser and ammonia survey of all 94 northern GLIMPSE Extended Green Objects (EGOs), a sample of massive young stellar objects (MYSOs) identified based on their extended 4.5 micron emission. We observed the ammonia (1,1), (2,2), and (3,3) inversion lines, and detect emission towards 97%, 63%, and 46% of our sample, respectively (median rms ~50 mK). The water maser detection rate is 68% (median rms ~0.11 Jy). The derived water maser and clump-scale gas properties are consistent with the identification of EGOs as young MYSOs. To explore the degree of variation among EGOs, we analyze subsamples defined based on MIR properties or maser associations. Water masers and warm dense gas, as indicated by emission in the higher-excitation ammonia transitions, are most frequently detected towards EGOs also associated with both Class I and II methanol masers. 95% (81%) of such EGOs are detected in water (ammonia(3,3)), compared to only 33% (7%) of EGOs without either methanol maser type. As populations, EGOs associated with Class I and/or II methanol masers have significantly higher ammonia linewidths, column densities, and kinetic temperatures than EGOs undetected in methanol maser surveys. However, we find no evidence for statistically significant differences in water maser properties (such as maser luminosity) among any EGO subsamples. Combining our data with the 1.1 mm continuum Bolocam Galactic Plane Survey, we find no correlation between isotropic water maser luminosity and clump number density. Water maser luminosity is weakly correlated with clump (gas) temperature and clump mass.
  • The ATLASGAL 870mum continuum survey conducted with the APEX telescope is the first survey covering the whole inner Galactic plane (60deg>l>-60deg & b<1.5deg) in submm continuum emission tracing the cold dust of dense and young star-forming regions. Here, we present the overall distribution of sources within our Galactic disk. The submm continuum emission is confined to a narrow range around the galactic plane, but shifted on average by ~0.07deg below the plane. Source number counts show strong enhancements toward the Galactic center, the spiral arms and toward prominent star-forming regions. Comparing the distribution of ATLASGAL dust continuum emission to that of young intermediate- to high-mass young stellar objects (YSOs) derived from Spitzer data, we find similarities as well as differences. In particular, the distribution of submm dust continuum emission is significantly more confined to the plane than the YSO distribution (FWHM of 0.7 and 1.1deg, corresponding to mean physical scale heights of approximately 46 and 80pc, respectively). While this difference may partly be caused by the large extinction from the dense submm cores, gradual dispersal of stellar distributions after their birth could also contribute to this effect. Compared to other tracers of Galactic structure, the ATLASGAL data are strongly confined to a narrow latitude strip around the Galactic plane.
  • The spectral energy distributions (SEDs), spanning the mid-infrared to millimeter wavelengths, of a sample of 13 high-mass protostellar objects (HMPOs) were studied using a large archive of 2-D axisymmetric radiative transfer models. Measurements from the Spitzer GLIMPSE and MIPSGAL surveys and the MSX survey were used in addition to our own surveys at millimeter and submillimeter wavelengths to construct the SEDs, which were then fit to the archive of models. These models assumed that stars of all masses form via accretion and allowed us to make estimates for the masses, luminosities and envelope accretion rates for the HMPOs. The models fit the observed SEDs well. The implied envelope accretion rates are high, $\approx 10^{-2.5} \msun/yr$, consistent with the accretion-based scenario of massive star formation. With the fitted accretion rates and with mass estimates of up to $\sim 20 \msun$ for these objects, it appears plausible that stars with stellar masses $M_{\ast} > 20 \msun$ can form via accretion.
  • We have observed the central region of the IR-dark cloud filament associated with IRAS 18507+0121 at millimeter wavelengths in CO(1-0), 13CO(1-0), and C18O(1-0) line emission and with the Spitzer Space Telescope at mid-IR wavelengths. Five massive outflows from two cloud cores were discovered. Three outflows are centered on or near an Ultracompact HII region (G34.4+0.23) while the remaining two outflows originate from the millimeter core G34.4+0.23 MM. Modeling of the SEDs of the mid-IR sources identified 31 young stellar objects in the filament with a combined stellar mass of ~127 +/- 27 Msun. An additional 22 sources were identified as probable cluster members based on the presence of strong 24 micron emission. The total star formation efficiency in the G34.4 cloud filament is estimated to be ~7% while the massive and intermediate mass star formation efficiency in the entire cloud filament is estimated to be roughly 2%. A comparison of the gravitational binding energy with the outflow kinetic energy suggests that the compact core containing G34.4+0.23 MM is being destroyed by its molecular outflows whereas the outflows associated with more massive core surrounding the G34.4 UC HII region are not likely to totally disrupt the cloud. Additionally, a qualitative evaluation of the region appears to suggest that stars in this region may have formed in two stages: first lower mass stars formed and then, a few Myrs later, the more massive stars began to form.
  • (Abridged) Spiral galaxies often have extended outflows that permeate beyond the region of the disk. Such outflows have been seen both in starburst galaxies, actively star forming galaxies and galaxies with an AGN. In the latter galaxies it is unknown whether the large-scale outflows are driven by star formation activity or purely by the active nucleus. The aim of our investigation is to study the frequency of extended minor-axis outflows in edge-on Seyfert galaxies to investigate the role of the AGN, the circumnuclear environment and star formation activity within the disk regions, and their importance for IGM enrichment on large scales. Narrowband imaging in two different ionizational stages (H-alpha and [OIII]) was performed to attempt a discrimination between processes associated with the active nucleus and those connected to star forming activity within the disk. The H-alpha morphology of the Seyfert galaxies is usually complex, but only in three out of 14 galaxies did we find evidence for minor axis disk outflows. At the sensitivity of our observations [OIII] emission is generally detected only in the nuclear region. Overall, our results show that extraplanar emission of similar brightness and extent as in the previously known cases of NGC3079 and NGC4388 is not common in Seyfert galaxies of otherwise similar properties. Comparison with our previous results shows that for nearby edge-on spiral galaxies star formation may be a more powerful mechanism for producing DIG than AGN activity. While in general AGN activity undoubtedly plays some role in driving minor-axis outflows, this probably requires higher AGN luminosities than are encountered in our small distance-limited sample.
  • We have identified 230 Tycho-2 Spectral Catalog stars that exhibit 8 micron mid-infrared extraphotospheric excesses in the MidCourse Space Experiment (MSX) and Spitzer Space Telescope Galactic Legacy MidPlane Survey Extraordinaire (GLIMPSE) surveys. Of these, 183 are either OB stars earlier than B8 in which the excess plausibly arises from a thermal bremsstrahlung component or evolved stars in which the excess may be explained by an atmospheric dust component. The remaining 47 stars have spectral classifications B8 or later and appear to be main sequence or late pre-main-sequence objects harboring circumstellar disks. Six of the 47 stars exhibit multiple signatures characteristic of pre-main-sequence circumstellar disks, including emission lines, near-infrared K-band excesses, and X-ray emission. Approximately one-third of the remaining 41 sources have emission lines suggesting relative youth. Of the 25 GLIMPSE stars with SST data at >24 microns, 20 also show an excess at 24 microns. Three additional objects have 24 micron upper limits consistent with possible excesses, and two objects have photospheric measurements at 24 microns. Six MSX sources had a measurement at wavelengths >8 microns. We modeled the excesses in 26 stars having two or more measurements in excess of the expected photospheres as single-component blackbodies. We determine probable disk temperatures and fractional infrared luminosities in the range 191 < T < 787 and 3.9x10^-4 < L_IR/L_* < 2.7x10^-1. We estimate a lower limit on the fraction of Tycho-2 Spectral Catalog main-sequence stars having mid-IR, but not near-IR, excesses to be 1.0+-0.3%.