Search

Most investigations of star-forming regions (SFRs) focus on their cores, where disks, jets, outflows and nebular clouds are commonplace. However, searches around the peripheries of stellar nurseries are almost non-existent, yet ought to contain additional young members, informing us about cloud dynamics and the dissolution of SFRs. “It’s a fishing trip” -- who’s going to offer that kind of observing time, right? Welcome to HYPERS*, a HST project which obtained 100 cycles of mostly optical photometry around the outskirts of SFRs by using the vacant second observing facility whilst the primary mission targeted a known nebulous region. Observations, completed in February 2022 consist of ~30 SFRs at distances 100-1000pc, size-scales 5-50 pc and populations of a few to a few thousand. Early results show: (1) dozens of emission disk candidates and a few potential jets; (2) many resolved multiple systems as close as ~0.05" that could not be detected without HST's powerful angular resolution; (3) large, previously unreported, cavities in the molecular clouds around the Eagle nebula. Our cross-match with several astronomical surveys indicates that we have lots of new stuff in there!

Brown dwarfs and planetary-mass objects have photometric and spectroscopic properties similar to exoplanet orbiting stars but they are in isolation, making their characterisation easier. We present a solid spectral sequence of substellar late-M and L-type dwarf members with effective temperatures below 2000 Kelvins of five young regions and open clusters: sigma Orionis (3-8 Myr; 350 pc), Upper Scorpius (5-10 Myr; 145 pc), the Pleiades (125 Myr; 135 pc), Coma Berenices (500 Myr; 87 pc), and the Hyades (650 Myr; 47 pc). We compare their photometric and spectroscopic properties to members of young moving groups and older field L dwarfs., Most of the diagrams and results presented in this poster have been published in refereed astronomical journals (or submitted): Sigma Orionis: Zapatero Osorio et al. (2017, ApJ, 842, 65) Upper Scorpius: Lodieu et al. (2018, MNRAS, 473, 2020) Pleiades: Zapatero Osorio et al. (2018, MNRAS, 475, 139) Hyades: Perez-Garrido et al. (2017, A&A, 599, 78); Martin et al. (2018, ApJ, 856, 40); Lodieu et al. (2018, A&A, 615, 12) Coma Berenices: Olivares et al. (2023, A&A, submitted), {"references":["Zapatero Osorio et al. (2017, ApJ, 842, 65)","Lodieu et al. (2018, MNRAS, 473, 2020)","Zapatero Osorio et al. (2018, MNRAS, 475, 139)","Pérez-Garrido et al. (2017, A&A, 599, 78)","Martin et al. (2018, ApJ, 856, 40)","Lodieu et al. (2018, A&A, 615, 12)","Olivares et al. (2023, A&A, submitted)"]}

We present high resolution Karl G. Jansky Very Large Array (VLA) observations of the protostar L1527 IRS at 7 mm, 1.3 cm, and 2 cm wavelengths. We detect the edge-on dust disk at all three wavelengths and find that it is asymmetric, with the southern side of the disk brighter than the northern side. We confirm this asymmetry through analytic modeling and also find that the disk is flared at 7 mm. We test the data against models including gap features in the intensity profile, and though we cannot rule such models out, they do not provide a statistically significant improvement in the quality of fit to the data. From these fits, we can however place constraints on allowed properties of any gaps that could be present in the true, underlying intensity profile. The physical nature of the asymmetry is difficult to associate with physical features due to the edge-on nature of the disk, but could be related to spiral arms or asymmetries seen in other imaging of more face-on disks., 20 pages, 5 figures, 3 tables; accepted for publication in the Astrophysical Journal

Chemistry along the star- and planet-formation sequence regulates how prebiotic building blocks -- carriers of the elements CHNOPS -- are incorporated into nascent planetesimals and planets. Spectral line observations across the electromagnetic spectrum are needed to fully characterize interstellar CHNOPS chemistry, yet to date there are only limited astrochemical constraints at THz frequencies. Here, we highlight advances to the study of CHNOPS astrochemistry that will be possible with the Orbiting Astronomical Satellite for Investigating Stellar Systems (OASIS). OASIS is a NASA mission concept for a space-based observatory that will utilize an inflatable 14-m reflector along with a heterodyne receiver system to observe at THz frequencies with unprecedented sensitivity and angular resolution. As part of a survey of H2O and HD towards ~100 protostellar and protoplanetary disk systems, OASIS will also obtain statistical constraints on the inventories of light hydrides including NH3 and H2S towards protoplanetary disks, as well as complex organics in protostellar hot corinos and envelopes. Line surveys of additional star-forming regions, including high-mass hot cores, protostellar outflow shocks, and prestellar cores, will also leverage the unique capabilities of OASIS to probe high-excitation organics and small hydrides, as is needed to fully understand the chemistry of these objects., Accepted to Frontiers in Astronomy and Space Sciences

Tao-Chung Ching et al., We present 850 μm dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the subfilaments, the magnetic fields are mainly parallel to the filamentary structures and smoothly connect to the magnetic fields of the main filament. We compare the POL-2 and Planck dust polarization observations to study the magnetic field structures of the DR21 filament on 0.1–10 pc scales. The magnetic fields revealed in the Planck data are well-aligned with those of the POL-2 data, indicating a smooth variation of magnetic fields from large to small scales. The plane-of-sky magnetic field strengths derived from angular dispersion functions of dust polarization are 0.6–1.0 mG in the DR21 filament and ∼0.1 mG in the surrounding ambient gas. The mass-to-flux ratios are found to be magnetically supercritical in the filament and slightly subcritical to nearly critical in the ambient gas. The alignment between column density structures and magnetic fields changes from random alignment in the low-density ambient gas probed by Planck to mostly perpendicular in the high-density main filament probed by James Clerk Maxwell Telescope. The magnetic field structures of the DR21 filament are in agreement with MHD simulations of a strongly magnetized medium, suggesting that magnetic fields play an important role in shaping the DR21 main filament and subfilaments., This work is supported by National Natural Science Foundation of China (NSFC) grant Nos. 11988101, U1931117, 11725313, and 12073061 and the CAS International Partnership Program of Chinese Academy of Sciences grant No. 114A11KYSB20160008. T.-C.C. is funded by Chinese Academy of Sciences Taiwan Young Talent Program grant No. 2018TW2JB0002. T.-C.C. and C.E. were supported by Special Funding for Advanced Users, budgeted and administrated by Center for Astronomical Mega-Science (CAMS), Chinese Academy of Sciences. K.P. is a Royal Society University Research Fellow, supported by grant No. URF\R1\211322. D.J. is supported by the National Research Council of Canada and by a Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant. P.M.K. is supported by the Ministry of Science and Technology (MoST) in Taiwan through grants 109-2112-M-001-022 and 110-2112-M-001-057. C.W.L. is supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (NRF-2019R1A2C1010851), and by the Korea Astronomy and Space Science Institute grant funded by the Korea government (MSIT; Project No. 2022-1-840-05). T.H. is supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) through the Mid-career Research Program (2019R1A2C1087045). W.K. was supported by the NRF grant funded by the MSIT (2021R1F1A1061794). C.E. acknowledges the financial support from grant RJF/2020/000071 as a part of Ramanujan Fellowship awarded by the Science and Engineering Research Board (SERB), Department of Science and Technology (DST), Govt. of India. F.P. acknowledges support from the Spanish State Research Agency (AEI) under grant No. PID2019-105552RB-C43. M.T. is supported by JSPS KAKENHI grant Nos. 18H05442, 15H02063, and 22000005. J.K. is supported by JSPS KAKENHI grant No. 19K14775. L.F. and F.K. acknowledge the support by the MoST in Taiwan through grant 107-2119-M-001-031-MY3 and Academia Sinica through grant AS-IA-106-M03. L.F. acknowledges the support by the MoST in Taiwan through grants 111-2811-M-005-007 and 109-2112-M-005-003-MY3. C.L.H.H. acknowledges the support of the NAOJ Fellowship and JSPS KAKENHI grants 18K13586 and 20K14527. F.K. is supported by the Spanish program Unidad de Excelencia María de Maeztu CEX2020-001058-M, financed by MCIN/AEI/10.13039/501100011033. K.Q. is partially supported by National Key R&D Program of China No. 2022YFA1603100, and acknowledges the National Natural Science Foundation of China (NSFC) grant U1731237. S.P.L. acknowledges grants from the Ministry of Science and Technology of Taiwan 106-2119-M-007-021-MY3 and 109-2112-M-007-010-MY3. Y.D. acknowledges the support of JSPS KAKENHI grants 25247016 and 18H01250. Y.S.D. is supported by the National Key R&D Program of China for grant No. 2022YFA1605300, and NSFC grants Nos. 12273051, 11933003.

Context. The chemistry of sulphur-bearing species in the interstellar medium remains poorly understood, but might play a key role in the chemical evolution of star-forming regions. Aims. Coupling laboratory experiments to observations of sulphur-bearing species in different parts of star-forming regions, we aim to understand the chemical behavior of the sulphur species in cold and warm regions of protostars, and we ultimately hope to connect them. Methods. We performed laboratory experiments in which we tested the reactivity of hydrogen sulfide (H2S) on a cold substrate with hydrogen and/or carbon monoxide (CO) under different physical conditions that allowed us to determine the products from sulphur reactions using a quadrupole mass spectrometer. The laboratory experiments were complemented by observations. We observed two luminous binary sources in the Cygnus-X star-forming complex, Cygnus X-N30 and N12, covering a frequency range of 329–361 GHz at a spatial resolution of 1′′5 with the SubMillimeter Array (SMA). This study was complemented by a 3 mm line survey of Cygnus X-N12 covering specific frequency windows in the frequency ranges 72.0–79.8 GHz at a spatial resolution of 34′′0–30′′0 and 84.2–115.5 GHz at a spatial resolution of 29′′0–21′′0, with the IRAM-30 m single-dish telescope. Column densities and excitation temperatures were derived under the local thermodynamic equilibrium approximation. Results. We find that OCS is a direct product from H2S reacting with CO and H under cold temperatures (T < 100 K) from laboratory experiments. OCS is therefore found to be an important solid-state S-reservoir. We identify several S-species in the cold envelope of Cyg X-N12, principally organo-sulphurs (H2CS, CS, OCS, CCS, C3S, CH3SH, and HSCN). For the hot cores of Cyg X-N12 and N30, only OCS, CS and H2CS were detected. We found a difference in the S-diversity between the hot core and the cold envelope of N12, which is likely due to the sensitivity of the observations toward the hot core of N12. Moreover, based on the hot core analysis of N30, the difference in S-diversity is likely driven by chemical processes rather than the low sensitivity of the observations. Furthermore, we found that the column density ratio of NCS/NSO is also an indicator of the warm (NCS/NSO > 1), cold (NCS/NSO < 1) chemistries within the same source. The line survey and molecular abundances inferred for the sulphur species are similar for protostars N30 and N12 and depends on the protostellar component targeted (i.e., envelope or hot core) rather than on the source itself. However, the spatial distribution of emission toward Cyg X-N30 shows differences compared to N12: toward N12, all molecular emission peaks on the two continuum sources, whereas emission is spatially distributed and shows variations within molecular families (N, O, and C families) toward N30. Moreover, this spatial distribution of all the identified S-species is offset from the N30 continuum peaks. The sulphur-bearing molecules are therefore good tracers to connect the hot and cold chemistry and to provide insight into the type of object that is observed.

We aim to characterize a diverse selection of dense, potentially star-forming cores, clumps, and clouds within the Milky Way in terms of their dust emission and SF activity. We studied 53 fields that have been observed in the JCMT SCUBA-2 continuum survey SCOPE and have been mapped with Herschel. We estimated dust properties by fitting Herschel observations with modified blackbody functions, studied the relationship between dust temperature and dust opacity spectral index $\beta$, and estimated column densities. We extracted clumps from the SCUBA-2 850 $\mu$m maps with the FellWalker algorithm and examined their masses and sizes. Clumps are associated with young stellar objects found in several catalogs. We estimated the gravitational stability of the clumps with virial analysis. The clumps are categorized as unbound starless, prestellar, or protostellar. We find 529 dense clumps, typically with high column densities from (0.3-4.8)$\times 10^{22}$ cm$^{-2}$, with a mean of (1.5$\pm$0.04)$\times10^{22}$ cm$^{-2}$, low temperatures ($T\sim $10-20 K), and estimated submillimeter $\beta$ =1.7$\pm$0.1. We detect a slight increase in opacity spectral index toward millimeter wavelengths. Masses of the sources range from 0.04 $M_\odot$ to 4259 $M_\odot$. Mass, linear size, and temperature are correlated with distance. Furthermore, the estimated gravitational stability is dependent on distance, and more distant clumps appear more virially bound. Finally, we present a catalog of properties of the clumps.Our sources present a large array of SF regions, from high-latitude, nearby diffuse clouds to large SF complexes near the Galactic center. Analysis of these regions will continue with the addition of molecular line data, which will allow us to study the densest regions of the clumps in more detail., Comment: 60 pages (article 18 pages), 73 figures To be published in A&A

A full understanding of high-mass star formation requires the study of one of the most elusive components of the energy balance in the interstellar medium: magnetic fields. We report Atacama Large Millimeter/submillimeter Array (ALMA) 1.2 mm, high-resolution (700 au) dust polarization and molecular line observations of the rotating hot molecular core embedded in the high-mass star-forming region IRAS 18089-1732. The dust continuum emission and magnetic field morphology present spiral-like features resembling a whirlpool. The velocity field traced by the H13CO+ (J = 3-2) transition line reveals a complex structure with spiral filaments that are likely infalling and rotating, dragging the field with them. We have modeled the magnetic field and find that the best model corresponds to a weakly magnetized core with a mass-to-magnetic-flux ratio (λ) of 8.38. The modeled magnetic field is dominated by a poloidal component, but with an important contribution from the toroidal component that has a magnitude of 30% of the poloidal component. Using the Davis-Chandrasekhar-Fermi method, we estimate a magnetic field strength of 3.5 mG. At the spatial scales accessible to ALMA, an analysis of the energy balance of the system indicates that gravity overwhelms turbulence, rotation, and the magnetic field. We show that high-mass star formation can occur in weakly magnetized environments, with gravity taking the dominant role., P.S. and B.W. were partially supported by a Grant-in-Aid for Scientific Research (KAKENHI number 18H01259) of the Japan Society for the Promotion of Science (JSPS). J.M.G. is supported by the Spanish grant AYA2017-84390-C2-R (AEI/FEDER, UE). C.L.H.H. acknowledges the support of the NAOJ Fellowship and JSPS KAKENHI grants 18K13586 and 20K14527. J.M.J.'s research was conducted in part at the SOFIA Science Center, which is operated by the Universities Space Research Association under contract NNA17BF53C with the National Aeronautics and Space Administration. K.T. was supported by JSPS KAKENHI grant No. 20H05645. Data analysis was in part carried out on the Multi-wavelength Data Analysis System operated by the Astronomy Data Center (ADC), National Astronomical Observatory of Japan. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2017.1.00101.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ.

We unveil the stellar wind driven shell of the luminous massive star-forming region of RCW 49 using SOFIA FEEDBACK observations of the [CII] 158 $\mu$m line. The complementary dataset of the $^{12}$CO and $^{13}$CO J = 3 - 2 transitions is observed by the APEX telescope and probes the dense gas toward RCW 49. Using the spatial and spectral resolution provided by the SOFIA and APEX telescopes, we disentangle the shell from a complex set of individual components of gas centered around RCW 49. We find that the shell of radius ~ 6 pc is expanding at a velocity of 13 km s$^{-1}$ toward the observer. Comparing our observed data with the ancillary data at X-Ray, infrared, sub-millimeter and radio wavelengths, we investigate the morphology of the region. The shell has a well defined eastern arc, while the western side is blown open and is venting plasma further into the west. Though the stellar cluster, which is ~ 2 Myr old gave rise to the shell, it only gained momentum relatively recently as we calculate the shell's expansion lifetime ~ 0.27 Myr, making the Wolf-Rayet star WR20a a likely candidate responsible for the shell's re-acceleration., Comment: 31 pages, 17 figures

Prestellar cores are self-gravitating dense and cold structures within molecular clouds where future stars are born. They are expected, at the stage of transitioning to the protostellar phase, to harbor centrally concentrated dense (sub)structures that will seed the formation of a new star or the binary/multiple stellar systems. Characterizing this critical stage of evolution is key to our understanding of star formation. In this work, we report the detection of high-density (sub)structures on the thousand-astronomical-unit (au) scale in a sample of dense prestellar cores. Through our recent ALMA observations toward the Orion Planck Galactic Cold Clumps, we have found five extremely dense prestellar cores, which have centrally concentrated regions of ∼2000 au in size, and several 107 cm-3 in average density. Masses of these centrally dense regions are in the range of 0.30 to 6.89M⊙. For the first time, our higher resolution observations (0.8' ∼ 320 au) further reveal that one of the cores shows clear signatures of fragmentation; such individual substructures/fragments have sizes of 800-1700 au, masses of 0.08 to 0.84M⊙, densities of 2 - 8 × 107 cm-3, and separations of ∼1200 au. The substructures are massive enough (≳0.1M⊙) to form young stellar objects and are likely examples of the earliest stage of stellar embryos that can lead to widely (∼1200 au) separated multiple systems., L. acknowledges the support from the international partnership program of the Chinese Academy of Sciences through grant No.114231KYSB20200009, National Natural Science Foundation of China (NSFC) through grant NSFC No.12073061, and Shanghai Pujiang Program 20PJ1415500. N.H. acknowledges MoST 108-2112-M-001-017 and MoST 109-2112-M-001-023 grants. G.G. acknowledges support from ANID project AFB 170002. L.B. acknowledges support from ANID project AFB-170002. S.L.Q. is supported by the National Natural Science Foundation of China under grant No. U1631237. D.J. is supported by NRC Canada and by an NSERC Discovery Grant. V.M.P. acknowledges support by the Spanish MINECO under project AYA2017-88754-P, and financial support from the State Agency for Research of the Spanish Ministry of Science and Innovation through the “Unit of Excellence María de Maeztu 2020-2023” award to the Institute of Cosmos Sciences (CEX2019-000918-M). C.W.L. is supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (NRF-2019R1A2C1010851). A.S. acknowledges financial support from the NSF through grant AST-1715876. The research was carried out in part at the Jet Propulsion Laboratory, which is operated for NASA by the California Institute of Technology. D. L. acknowledges support from NSFC No. 11911530226 and 11725313. K.T. was supported by JSPS KAKENHI grant No. 20H05645. J.H. thanks the NSFC grant No. 11873086 and Yunnan Province of China (No. 2017HC018).This work is sponsored (in part) by the Chinese Academy of Sciences (CAS), through a grant to the CAS South America Center for Astronomy (CASSACA) in Santiago, Chile.

We analyze surveys of molecular cloud structures defined by tracers ranging from CO $J = 1-0$ through $^{13}$CO $J = 1-0$ to dust emission together with NH$_3$ data. The mean value of the virial parameter and the fraction of mass in bound structures depends on the method used to identify structures. Generally, the virial parameter decreases and the fraction of mass in bound structures increases with the effective density of the tracer, the surface density and mass of the structures, and the distance from the center of a galaxy. For the most complete surveys of structures in the Galaxy defined by CO $J = 1-0$, the fraction of mass that is in bound structures is 0.19. For catalogs of other galaxies based on CO $J = 2-1$, the fraction is 0.35. These results offer substantial alleviation of the fundamental problem of slow star formation. If only clouds found to be bound are counted and they are assumed to collapse in a free-fall time at their mean cloud density, the sum over all clouds in a complete survey of the Galaxy yields a predicted star formation rate of 46 solar masses per year, a factor of 6.5 less than if all clouds are bound., Comment: 26 pages, 14 figures, accepted to ApJ

Measurement of magnetic field strengths in a molecular cloud is essential for determining the criticality of magnetic support against gravitational collapse. In this paper, as part of the JCMT BISTRO survey, we suggest a new application of the Davis-Chandrasekhar-Fermi (DCF) method to estimate the distribution of magnetic field strengths in the OMC-1 region. We use observations of dust polarization emission at 450 $\mu$m and 850 $\mu$m, and C$^{18}$O (3-2) spectral line data obtained with the JCMT. We estimate the volume density, the velocity dispersion and the polarization angle dispersion in a box, 40$''$ $\times$ 40$''$ (5$\times$5 pixels), which moves over the OMC-1 region. By substituting three quantities in each box to the DCF method, we get magnetic field strengths over the OMC-1 region. We note that there are very large uncertainties in inferred field strengths, as discussed in detail in this paper. The field strengths vary from 0.8 to 26.4 mG and their mean value is about 6 mG. Additionally, we obtain maps of the mass-to-flux ratio in units of a critical value and the Alfv$\acute{e}$n mach number. The central parts of the BN-KL and South (S) clumps in the OMC-1 region are magnetically supercritical, so the magnetic field cannot support the clumps against gravitational collapse. However, the outer parts of the region are magnetically subcritical. The mean Alfv$\acute{e}$n mach number is about 0.4 over the region, which implies that the magnetic pressure exceeds the turbulent pressure in the OMC 1 region., Comment: 28 pages

Non-thermal particles and high-energy radiation can play a role in the dynamical processes in star-forming regions and provide an important piece of the multiwavelength observational picture of their structure and components. Powerful stellar winds and supernovae in compact clusters of massive stars and OB associations are known to be favourable sites of high-energy particle acceleration and sources of non-thermal radiation and neutrinos. Namely, young massive stellar clusters are likely sources of the PeV (petaelectronvolt) regime cosmic rays (CRs). They can also be responsible for the cosmic ray composition, e.g., 22Ne/20Ne anomalous isotopic ratio in CRs. Efficient particle acceleration can be accompanied by super-adiabatic amplification of the fluctuating magnetic fields in the systems converting a part of kinetic power of the winds and supernovae into the magnetic energy through the CR-driven instabilities. The escape and CR propagation in the vicinity of the sources are affected by the non-linear CR feedback. These effects are expected to be important in starburst galaxies, which produce high-energy neutrinos and gamma-rays. We give a brief review of the theoretical models and observational data on high-energy particle acceleration and their radiation in star-forming regions with young stellar population., 41 pages, Space Science Reviews v.216, id 42, topical collection Star Formation

H$_2$O maser emission associated with the massive star formation region W49N were observed with the Space-VLBI mission RadioAstron. The procedure for processing of the maser spectral line data obtained in the RadioAstron observations is described. Ultra-fine spatial structures in the maser emission were detected on space-ground baselines of up to 9.6 Earth diameters. The correlated flux densities of these features range from 0.1% to 0.6% of the total flux density. These low values of correlated flux density are probably due to turbulence either in the maser itself or in the interstellar medium., Comment: Accepted for publication in Advances in Space Research

Recent high angular resolution (≃40 mas) ALMA observations at 1.14 mm resolve a compact (R ≃ 200 au), flattened dust structure perpendicular to the HH 80-81 jet emanating from the GGD 27-MM1 high-mass protostar, making it a robust candidate for a true accretion disk. The jet-disk system (HH 80-81/GGD 27-MM1) resembles those found in association with low- A nd intermediate-mass protostars. We present radiative transfer models that fit the 1.14 mm ALMA dust image of this disk, which allow us to obtain its physical parameters and predict its density and temperature structure. Our results indicate that this accretion disk is compact (R disk ≃ 170 au) and massive (≃5 M o), at about 20% of the stellar mass of ≃20 M o. We estimate the total dynamical mass of the star-disk system from the molecular line emission, finding a range between 21 and 30 M o, which is consistent with our model. We fit the density and temperature structures found by our model with power-law functions. These results suggest that accretion disks around massive stars are more massive and hotter than their low-mass siblings, but they still are quite stable. We also compare the temperature distribution in the GGD 27-MM1 disk with that found in low- A nd intermediate-mass stars and discuss possible implications for the water snow line. We have also carried out a study of the distance based on Gaia DR2 data and the population of young stellar objects in this region and from the extinction maps. We conclude that the source distance is within 1.2 and 1.4 kpc, closer than what was derived in previous studies (1.7 kpc). © 2020. The American Astronomical Society. All rights reserved., We are grateful to Nuria Calvet for making D'Alessio's code available to us and for very useful discussions. We would also like to deeply thank Guillem Anglada for his critical comments on this work that have undoubtedly improved it very substantially. We thank our anonymous referee for her or his very helpful comments and suggestions to improve our manuscript. We also thank Rosine Lallement for the support on STILISM. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2015.1.00480.S. ALMA is a partnership of ESO (representing its member states), NSF (USA), and NINS (Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan) and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO, and NAOJ. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC; https://www.cosmos.esa.int/web/gaia/dpac/consortium).Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. N.A.-L., M.O., J.M.G, G.B., R.E., and J.M.T. are supported by the Spanish grant AYA2017-84390-C2-R (AEI/FEDER, UE). R.E. is also supported by MDM-2014-0369 of ICCUB (Unidad de Excelencia "Maria de Maeztu"). M.O. acknowledges financial support from the State Agency for Research of the Spanish MCIU through the "Center of Excellence Severo Ochoa" award for the Instituto de Astrofisica de Andalucia (SEV-2017-0709). S.C. acknowledges support from DGAPA, UNAM, and CONACyT, Mexico. M.F.-L. acknowledges support from the LACEGAL project, which has received funding from the European Union's Horizon 2020 Research and Innovation Program under the Marie SklodowskaCurie grant agreement No. 734374. C.C.-G. and R.G.M. acknowledges support from UNAM-DGAPA PAPIIT project IN 108218 and IN 104319. J.K. is supported by MEXT KAKENHI grant No. 19K14755.

We present 3, 15, and 33 GHz imaging towards galaxy nuclei and extranuclear star-forming regions using the Karl G. Jansky Very Large Array as part of the Star Formation in Radio Survey. With $3-33$ GHz radio spectra, we measured the spectral indices and corresponding thermal (free-free) emission fractions for a sample of 335 discrete regions having significant detections in at least two radio bands. After removing 14 likely background galaxies, we find that the median thermal fraction at 33 GHz is $92 \pm 0.8\%$ with a median absolute deviation of $11\%$, when a two-component power-law model is adopted to fit the radio spectrum. Limiting the sample to 238 sources that are confidently identified as star-forming regions, and not affected by potential AGN contamination (i.e., having galactocentric radii $r_{\rm G} \geq 250$ pc), results in a median thermal fraction of $93 \pm 0.8 \%$ with a median absolute deviation of $10\%$. We further measure the thermal fraction at 33 GHz for 163 regions identified at 7" resolution to be $94 \pm 0.8 \%$ with a median absolute deviation of $8\%$. Together, these results confirm that free-free emission dominates the radio spectra of star-forming regions on scales up to $\sim$500 pc in normal star-forming galaxies. We additionally find a factor of $\sim$1.6 increase in the scatter of the measured spectral index and thermal fraction distributions as a function of decreasing galactocentric radius. This trend is likely reflective of the continuous star-formation activity occurring in the galaxy centers, resulting a larger contribution of diffuse nonthermal emission relative to star-forming regions in the disk., Comment: 57 pages, 23 figures, 7 tables, accepted for publication in ApJS

We present the first spatially resolved study of molecular gas in the vicinity of a Gamma Ray Burst, using CO(2-1) emission line observations with the Atacama Large Millimetre Array (ALMA) at ~50 pc scales. The host galaxy of GRB 980425 contains a ring of high column density HI gas which is likely to have formed due to a collision between the GRB host and its companion galaxy, within which the GRB is located. We detect eleven molecular gas clumps in the galaxy, seven of which are within the gas ring. The clump closest to the GRB position is at a projected separation of ~280 pc. Although it is plausible that the GRB progenitor was ejected from clusters formed in this clump, we argue that the in situ formation of the GRB progenitor is the most likely scenario. We measure the molecular gas masses of the clumps and find them to be sufficient for forming massive star clusters. The molecular gas depletion times of the clumps show a variation of ~2 dex, comparable with the large variation in depletion times found in starburst galaxies in the nearby Universe. This demonstrates the presence of starburst modes of star formation on local scales in the galaxy, even while the galaxy as a whole cannot be categorised as a starburst based on its global properties. Our findings suggest that the progenitor of GRB 9802425 was originated in a young massive star cluster formed in the starburst mode of star formation., Comment: ApJ (in press), 12 pages, 4 figures, 1 table

The physical evolution of Young Stellar Objects (YSOs) is accompanied by an enrichment of the molecular complexity, mainly triggered by the heating and energetic processing of the astrophysical ices. In this paper, a study of how the ice column density varies across the protostellar evolution has been performed. Tabulated data of H$_2$O, CO$_2$, CH$_3$OH, HCOOH observed by ground- and space-based telescopes toward 27 early-stage YSOs were taken from the literature. The observational data shows that ice column density and spectral index ($\alpha$), used to classify the evolutionary stage, are well correlated. A 2D continuum radiative transfer simulation containing bare and grains covered by ices at different levels of cosmic-ray processing were used to calculate the Spectral Energy Distributions (SEDs) in different angle inclinations between face-on and edge-on configuration. The H$_2$O:CO$_2$ ice mixture was used to address the H$_2$O and CO$_2$ column density variation whereas the CH$_3$OH and HCOOH are a byproduct of the virgin ice after the energetic processing. The simulated spectra were used to calculate the ice column densities of YSOs in an evolutionary sequence. As a result, the models show that the ice column density variation of HCOOH with $\alpha$ can be justified by the envelope dissipation and ice energetic processing. On the other hand, the ice column densities are mostly overestimated in the cases of H$_2$O, CO$_2$ and CH$_3$OH, even though the physical and cosmic-ray processing effects are taken into account., Comment: Accepted for publication in The Astrophysical Journal. 29 pages, 16 figures

We present velocity resolved SOFIA/upGREAT observations of [OI] and [CII] lines toward a Class I protostar, L1551 IRS 5, and its outflows. The SOFIA observations detect [OI] emission toward only the protostar and [CII] emission toward the protostar and the red-shifted outflow. The [OI] emission has a width of $\sim$100 km s$^{-1}$ only in the blue-shifted velocity, suggesting an origin in shocked gas. The [CII] lines are narrow, consistent with an origin in a photodissociation region. Differential dust extinction from the envelope due to the inclination of the outflows is the most likely cause of the missing red-shifted [OI] emission. Fitting the [OI] line profile with two Gaussian components, we find one component at the source velocity with a width of $\sim$20 km s$^{-1}$ and another extremely broad component at -30 km s$^{-1}$ with a width of 87.5 km s$^{-1}$, the latter of which has not been seen in L1551 IRS 5. The kinematics of these two components resemble cavity shocks in molecular outflows and spot shocks in jets. Radiative transfer calculations of the [OI], high-J CO, and H$_2$O lines in the cavity shocks indicate that [OI] dominates the oxygen budget, making up more than 70% of the total gaseous oxygen abundance and suggesting [O]/[H] of $\sim$1.5$\times$10$^{-4}$. Attributing the extremely broad [OI] component to atomic winds, we estimate the intrinsic mass loss rate of (1.3$\pm$0.8)$\times$ 10$^{-6}$ M$_{\odot}$ yr$^{-1}$. The intrinsic mass loss rates derived from low-J CO, [OI], and HI are similar, supporting the model of momentum-conserving outflows, where the atomic wind carries most momentum and drives the molecular outflows., 19 pages, 14 figures, accepted for publication in The Astrophysical Journal

Context: Ultraviolet radiation (UV) influences the physics and chemistry of star-forming regions, but its properties and significance in the immediate surroundings of low-mass protostars are still poorly understood. Aims: We aim to extend the use of the CN/HCN ratio, already established for high-mass protostars, to the low-mass regime to trace and characterize the UV field around low-mass protostars on $\sim 0.6\times0.6$ pc scales. Methods: We present $5'\times5'$ maps of the Serpens Main Cloud encompassing 10 protostars observed with the EMIR receiver at the IRAM 30 m telescope in CN 1-0, HCN 1-0, CS 3-2, and some of their isotopologues. The radiative-transfer code RADEX and the chemical model Nahoon are used to determine column densities of molecules, gas temperature and density, and the UV field strength, $G_\mathrm{0}$. Results: The spatial distribution of HCN and CS are well-correlated with CO 6-5 emission that traces outflows. The CN emission is extended from the central protostars to their immediate surroundings also tracing outflows, likely as a product of HCN photodissociation. The ratio of CN to HCN total column densities ranges from $\sim$1 to 12 corresponding to G$_0$ $\approx$ $10^{1}-10^{3}$ for gas densities and temperatures typical for outflows of low-mass protostars. Conclusions: UV radiation associated with protostars and their outflows is indirectly identified in a significant part of the Serpens Main low-mass star-forming region. Its strength is consistent with the values obtained from the OH and H$_2$O ratios observed with Herschel and compared with models of UV-illuminated shocks. From a chemical viewpoint, the CN to HCN ratio is an excellent tracer of UV fields around low- and intermediate-mass star-forming regions., Comment: 32 pages, 25 figures, accepted by A\&A

The rarity of young massive stars combined with the fact that they are often deeply embedded has limited the understanding of their formation. Ground based mid-infrared (IR) interferometry is one way of securing the spatial resolution required to study massive young stellar objects (MYSOs) and as the spatial-frequency coverage of such observations is often incomplete, direct-imaging can be supplementary to such a dataset. By consolidating these observations with modelling, the features of a massive protostellar environment can be constrained. This work simultaneously fits the aforementioned observations and a spectral energy distribution (SED) with a 2.5D radiative transfer model, providing an extensive view of the physical characteristics of the accreting regions of the MYSO G305.20+0.21. The high-resolution observations were obtained using the Very Large Telescope's MIDI and VISIR instruments, producing visibilities in the N-band and near-diffraction-limited imaging in the Q-band respectively. A model including a central protostar with a luminosity of $\sim$5$\times$10$^{4}$L$_{\odot}$ surrounded by a low-density bipolar cavity, a flared 1M$_{\odot}$ disk and an envelope provides a sufficient fit all three types of observation. The need to include a disk in the model implies that this MYSO follows a scaled-up version of the low-mass star formation process. The weak silicate absorption feature within the SED requires low-density envelope cavities to be successfully fit and is an atypical characteristic in comparison to previously studied MYSOs. Additionally, the inner radius of the disk must be three times the dust sublimation radius to satisfy the MIDI visibilities. The low density, low extinction environment implies the object is a more evolved MYSO and this combined with large inner radius of the disk suggests that it could be an example of a transitional disk around an MYSO., Comment: 16 pages, 8 figures

We report the detection of four new hot corino sources, G211.47-19.27S, G208.68-19.20N1, G210.49-19.79W and G192.12-11.10 from a survey study of Planck Galactic Cold Clumps in the Orion Molecular Cloud Complex with the Atacama Compact Array (ACA). Three sources had been identified as low mass Class 0 protostars in the Herschel Orion Protostar Survey (HOPS). One source in the lambda Orionis region is firstly reported as a protostellar core. We have observed abundant complex organic molecules (COMs), primarily methanol but also other oxygen-bearing COMs (in G211.47-19.27S and G208.68-19.20N1) and the molecule of prebiotic interest NH2CHO (in G211.47-19.27S), signifying the presence of hot corinos. While our spatial resolution is not sufficient for resolving most of the molecular emission structure, the large linewidth and high rotational temperature of COMs suggest that they likely reside in the hotter and innermost region immediately surrounding the protostar. In G211.47-19.27S, the D/H ratio of methanol ([CH2DOH]/[CH3OH]) and the 12C/13C ratio of methanol ([CH3OH]/[13CH3OH]) are comparable to those of other hot corinos. Hydrocarbons and long carbon-chain molecules such as c-C3H2 and HCCCN are also detected in the four sources, likely tracing the outer and cooler molecular envelopes., 37 pages, 51 figures, to be published in ApJ

In the high-mass star-forming region G35.20-0.74N, small scale (about 800 AU) chemical segregation has been observed in which complex organic molecules containing the CN group are located in a small location. We aim to determine the physical origin of the large abundance difference (about 4 orders of magnitude) in complex cyanides within G35.20-0.74 B, and we explore variations in age, gas and dust temperature, and gas density. We performed gas-grain astrochemical modeling experiments with exponentially increasing (coupled) gas and dust temperature rising from 10 to 500 K at constant H$_2$ densities of 10$^7$, 10$^8$, and 10$^9$ cm$^{-3}$. We tested the effect of varying the initial ice composition, cosmic-ray ionization rate, warm-up time (over 50, 200, and 1000 kyr), and initial (10, 15, and 25 K) and final temperatures (300 and 500 K). Varying the initial ice compositions within the observed and expected ranges does not noticeably affect the modeled abundances indicating that the chemical make-up of hot cores is determined in the warm-up stage. Complex cyanides vinyl and ethyl cyanide (CH$_2$CHCN and C$_2$H$_5$CN, respectively) cannot be produced in abundances (versus H$_2$) greater than 5x10$^{-10}$ for CH$_2$CHCN and 2x10$^{-10}$ for C$_2$H$_5$CN with a fast warm-up time (52 kyr), while the lower limit for the observed abundance of C$_2$H$_5$CN toward source B3 is 3.4x10$^{-10}$. Complex cyanide abundances are reduced at higher initial temperatures and increased at higher cosmic-ray ionization rates. Reproducing the observed abundances toward G35.20-0.74 Core B3 requires a fast warm-up at a high cosmic-ray ionization rate (1x10$^{-16}$ s$^{-1}$) at a high gas density (>10$^9$ cm$^{-3}$). G35.20-0.74 source B3 only needs to be about 2000 years older than B1/B2 for the observed chemical difference to be present. (This abstract has been shortened), 14 pages (plus 11 pages appendices), 12 figures main text

It is not well known what drives the chemistry of a protostellar envelope, in particular the role of the stellar mass and the outflows on its chemical enrichment. We study the chemical structure of NGC 1333 IRAS 4A in order to (i) investigate the influence of the outflows on the chemistry, (ii) constrain the age of our object, (iii) compare it with a typical high-mass protostellar envelope. In our analysis we use JCMT line mapping and HIFI pointed spectra. To study the influence of the outflow on the degree of deuteration, we compare JCMT maps of HCO+ and DCO+ with non-LTE (RADEX) models in a region that spatially covers the outflow activity of IRAS 4A. To study the envelope chemistry, we derive empirical molecular abundance profiles for the observed species using the radiative transfer code (RATRAN) and adopting a 1D dust density/temperature profile from the literature. We compare our best-fit observed abundance profiles with the predictions from the time dependent gas grain chemical code (ALCHEMIC). The CO, HCN, HNC and CN abundance require an enhanced UV field which points towards an outflow cavity. The abundances (wrt H2) are 1 to 2 orders of magnitude lower than those observed in the high mass protostellar envelope (AFGL 2591), while they are found to be similar within factors of a few with respect to CO. Differences in UV radiation may be responsible for such chemical differentiation, but temperature differences seem a more plausible explanation. The CH3OH modeled abundance profile points towards an age of > 4x10^4 yrs for IRAS 4A. The spatial distribution of H2D+ differs from that of other deuterated species, indicating an origin from a foreground colder layer (, 26 pages, 23 Figures, 6 tables, accepted for publication in A&A

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., Comment: accepted for publication in A&A Letters

Context. Circumstellar discs around massive stars could mediate the accretion onto the star from the infalling envelope, and could minimize the effects of radiation pressure. Despite such a crucial role, only a few convincing candidates have been provided for discs around deeply embedded O-type (proto)stars. Aims: In order to establish whether disc-mediated accretion is the formation mechanism for the most massive stars, we have searched for circumstellar, rotating discs around a limited sample of six luminous (>105L☉) young stellar objects. These objects were selected on the basis of their IR and radio properties in order to maximize the likelihood of association with disc+jet systems. Methods: We used ALMA with 0.̋2 resolution to observe a large number of molecular lines typical of hot molecular cores. In this paper we limit our analysis to two disc tracers (methyl cyanide, CH3CN, and its isotopologue, 13CH3CN), and an outflow tracer (silicon monoxide, SiO). Results: We reveal many cores, although their number depends dramatically on the target. We focus on the cores that present prominent molecular line emission. In six of these a velocity gradient is seen across the core,three of which show evidence of Keplerian-like rotation. The SiO data reveal clear but poorly collimated bipolar outflow signatures towards two objects only. This can be explained if real jets are rare (perhaps short-lived) in very massive objects and/or if stellar multiplicity significantly affects the outflow structure.For all cores with velocity gradients, the velocity field is analysed through position-velocity plots to establish whether the gas is undergoing rotation with νrot ∝ R- α, as expected for Keplerian-like discs. Conclusions: Our results suggest that in three objects we are observing rotation in circumstellar discs, with three more tentative cases, and one core where no evidence for rotation is found. In all cases but one, we find that the gas mass is less than the mass of any embedded O-type star, consistent with the (putative) discs undergoing Keplerian-like rotation. With the caveat of low number statistics, we conclude that the disc detection rate could be sensitive to the evolutionary stage of the young stellar object. In young, deeply embedded sources, the evidence for discs could be weak because of confusion with the surrounding envelope, while in the most evolved sources the molecular component of the disc could have already been dispersed. Only in those objects that are at an intermediate stage of the evolution would the molecular disc be sufficiently prominent and relatively less embedded to be detectable by mm/submm observations.

We present results from a comprehensive submillimeter spectral survey toward the source Orion South, based on data obtained with the Heterodyne Instrument for the Far-Infrared instrument on board the Herschel Space Observatory, covering the frequency range of 480 to 1900 GHz. We detect 685 spectral lines with signal-to-noise ratios (S/Ns) > 3 sigma, originating from 52 different molecular and atomic species. We model each of the detected species assuming conditions of Local Thermodynamic Equilibrium. This analysis provides an estimate of the physical conditions of Orion South (column density, temperature, source size, and VLSR). We find evidence for three different cloud components: a cool (T-ex similar to 20-40 K), spatially extended (> 60 ''), and quiescent (Delta V-FWHM similar to 4 km s(-1)) component; a warmer (T-ex similar to 80-100 K), less spatially extended (similar to 30 ''), and dynamic (Delta V-FWHM similar to 8 km s(-1)) component, which is likely affected by embedded outflows; and a kinematically distinct region (T-ex > 100 K; V-LSR similar to 8 km s(-1)), dominated by emission from species that trace ultraviolet irradiation, likely at the surface of the cloud. We find little evidence for the existence of a chemically distinct "hot-core" component, likely due to the small filling factor of the hot core or hot cores within the Herschel beam. We find that the chemical composition of the gas in the cooler, quiescent component of Orion South more closely resembles that of the quiescent ridge in Orion-KL. The gas in the warmer, dynamic component, however, more closely resembles that of the Compact Ridge and Plateau regions of Orion-KL, suggesting that higher temperatures and shocks also have an influence on the overall chemistry of Orion South.

The era of interferometric observations leads to the need of a more and more precise description of physical structure and dynamics of star-forming regions, from pre-stellar cores to proto-planetary disks. The molecular and dust continuum emission can be described with multiple physical components. To compare with the observations, a precise and complex radiative transfer modelling of these regions is required. I have developed during this thesis a standalone application called GASS (Generator of Astrophysical Sources Structures, Quénard et al., submitted) for this purpose. Thanks to its interface, GASS allows to create, manipulate, and mix several different physical components such as spherical sources, disks, and outflows. In this thesis, I used GASS to work on different astrophysical cases and, among them, I studied in details the water and deuterated water emission in the pre-stellar core L1544 (Quénard et al., 2016) and the emission of ions in the low-mass proto-star IRAS16293-2422 (Quénard et al., submitted).; L'ère des observations interférométriques mène à la nécessité d'une description plus précise de la structure physique et de la dynamique des régions de formation d'étoiles, des coeurs pré-stellaires et des disques proto-planétaires. L'émission moléculaire et du continuum de la poussière peuvent être décrites par de multiples composantes physiques. Pour comparer avec les observations, un modèle de transfert radiatif précis et complexe de ces régions est nécessaire. J'ai développé au cours de cette thèse une application autonome appelée GASS (Generator of Astrophysical Sources Structures, Quénard et al., soumis) à cette fin. Grâce à son interface, GASS permet de créer, de manipuler et de mélanger différents composants physiques tels que des sources sphériques, des disques et des outflows. Dans cette thèse, j'ai utilisé GASS pour travailler sur différents cas astrophysiques et, entre autres, j'ai étudié en détail l'eau et l'émission de l'eau deutérée dans le coeur pré-stellaire L1544 (Quénard et al., 2016) ainsi que l'émission des ions dans la proto-étoile de faible masse IRAS16293-2422 (Quénard et al., soumis).

The era of interferometric observations leads to the need of a more and more precise description of physical structure and dynamics of star-forming regions, from pre-stellar cores to proto-planetary disks. The molecular and dust continuum emission can be described with multiple physical components. To compare with the observations, a precise and complex radiative transfer modelling of these regions is required. I have developed during this thesis a standalone application called GASS (Generator of Astrophysical Sources Structures, Quénard et al., submitted) for this purpose. Thanks to its interface, GASS allows to create, manipulate, and mix several different physical components such as spherical sources, disks, and outflows. In this thesis, I used GASS to work on different astrophysical cases and, among them, I studied in details the water and deuterated water emission in the pre-stellar core L1544 (Quénard et al., 2016) and the emission of ions in the low-mass proto-star IRAS16293-2422 (Quénard et al., submitted).; L'ère des observations interférométriques mène à la nécessité d'une description plus précise de la structure physique et de la dynamique des régions de formation d'étoiles, des coeurs pré-stellaires et des disques proto-planétaires. L'émission moléculaire et du continuum de la poussière peuvent être décrites par de multiples composantes physiques. Pour comparer avec les observations, un modèle de transfert radiatif précis et complexe de ces régions est nécessaire. J'ai développé au cours de cette thèse une application autonome appelée GASS (Generator of Astrophysical Sources Structures, Quénard et al., soumis) à cette fin. Grâce à son interface, GASS permet de créer, de manipuler et de mélanger différents composants physiques tels que des sources sphériques, des disques et des outflows. Dans cette thèse, j'ai utilisé GASS pour travailler sur différents cas astrophysiques et, entre autres, j'ai étudié en détail l'eau et l'émission de l'eau deutérée dans le coeur pré-stellaire L1544 (Quénard et al., 2016) ainsi que l'émission des ions dans la proto-étoile de faible masse IRAS16293-2422 (Quénard et al., soumis).

The era of interferometric observations leads to the need of a more and more precise description of physical structure and dynamics of star-forming regions, from pre-stellar cores to proto-planetary disks. The molecular and dust continuum emission can be described with multiple physical components. To compare with the observations, a precise and complex radiative transfer modelling of these regions is required. I have developed during this thesis a standalone application called GASS (Generator of Astrophysical Sources Structures, Quénard et al., submitted) for this purpose. Thanks to its interface, GASS allows to create, manipulate, and mix several different physical components such as spherical sources, disks, and outflows. In this thesis, I used GASS to work on different astrophysical cases and, among them, I studied in details the water and deuterated water emission in the pre-stellar core L1544 (Quénard et al., 2016) and the emission of ions in the low-mass proto-star IRAS16293-2422 (Quénard et al., submitted).; L'ère des observations interférométriques mène à la nécessité d'une description plus précise de la structure physique et de la dynamique des régions de formation d'étoiles, des coeurs pré-stellaires et des disques proto-planétaires. L'émission moléculaire et du continuum de la poussière peuvent être décrites par de multiples composantes physiques. Pour comparer avec les observations, un modèle de transfert radiatif précis et complexe de ces régions est nécessaire. J'ai développé au cours de cette thèse une application autonome appelée GASS (Generator of Astrophysical Sources Structures, Quénard et al., soumis) à cette fin. Grâce à son interface, GASS permet de créer, de manipuler et de mélanger différents composants physiques tels que des sources sphériques, des disques et des outflows. Dans cette thèse, j'ai utilisé GASS pour travailler sur différents cas astrophysiques et, entre autres, j'ai étudié en détail l'eau et l'émission de l'eau deutérée dans le coeur pré-stellaire L1544 (Quénard et al., 2016) ainsi que l'émission des ions dans la proto-étoile de faible masse IRAS16293-2422 (Quénard et al., soumis).

We have used the Odin satellite to obtain strip maps of the ground-state rotational transitions of ortho-water and ortho-ammonia, as well as CO(5-4) and 13CO(5-4) across the PDR, and H218O in the central position. A physi-chemical inhomogeneous PDR model was used to compute the temperature and abundance distributions for water, ammonia and CO. A multi-zone escape probability method then calculated the level populations and intensity distributions. These results are compared to a homogeneous model computed with an enhanced version of the RADEX code. H2O, NH3 and 13CO show emission from an extended PDR with a narrow line width of ~3 kms. Like CO, the water line profile is dominated by outflow emission, however, mainly in the red wing. The PDR model suggests that the water emission mainly arises from the surfaces of optically thick, high density clumps with n(H2)>10^6 cm^-3 and a clump water abundance, with respect to H2, of 5x10^-8. The mean water abundance in the PDR is 5x10^-9, and between ~2x10^-8 -- 2x10^-7 in the outflow derived from a simple two-level approximation. Ammonia is also observed in the extended clumpy PDR, likely from the same high density and warm clumps as water. The average ammonia abundance is about the same as for water: 4x10^-9 and 8x10^-9 given by the PDR model and RADEX, respectively. The similarity of water and ammonia PDR emission is also seen in the almost identical line profiles observed close to the bright rim. Around the central position, ammonia also shows some outflow emission although weaker than water in the red wing. Predictions of the H2O(110-101) and (111-000) antenna temperatures across the PDR are estimated with our PDR model for the forthcoming observations with the Herschel Space Observatory., 13 pages, 14 figures, 10 tables. Accepted for publication in Astronomy & Astrophysics 14 November 2008

The Delta-variance analysis is an efficient tool for measuring the structural scaling behaviour of interstellar turbulence in astronomical maps. In paper I we proposed essential improvements to the Delta-variance analysis. In this paper we apply the improved Delta-variance analysis to i) a hydrodynamic turbulence simulation with prominent density and velocity structures, ii) an observed intensity map of rho Oph with irregular boundaries and variable uncertainties of the different data points, and iii) a map of the turbulent velocity structure in the Polaris Flare affected by the intensity dependence on the centroid velocity determination. The tests confirm the extended capabilities of the improved Delta-variance analysis. Prominent spatial scales were accurately identified and artifacts from a variable reliability of the data were removed. The analysis of the hydrodynamic simulations showed that the injection of a turbulent velocity structure creates the most prominent density structures are produced on a scale somewhat below the injection scale. The new analysis of a rho Oph continuum map reveals an intermediate stage in the molecular cloud evolution showing both signatures of the typical molecular cloud scaling behaviour and the formation of condensed cores. When analysing the velocity structure of the Polaris Flare we show that a universal power law connects scales from 0.03 pc to 3 pc. However, a plateau in the Delta-variance spectrum around 5 pc indicates that the visible large-scale velocity gradient is not converted directly into a turbulent cascade., Comment: Accepted for publication in A&A, Section 6

Aims. We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution.Methods. A number of benchmark models have been created, covering low and high gas densities n = 10(3), 10(5.5) cm(-3) and far ultraviolet intensities chi = 10, 10(5) in units of the Draine field (FUV: 6 Results. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes such as the effect of PAHs on the electron density or the temperature dependence of the dissociation of CO by cosmic ray induced secondary photons, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly. Nevertheless, the remaining spread in the computed observables such as the atomic fine-structure line intensities serves as a warning that there is still a considerable uncertainty when interpreting astronomical data with our models.

Aims. We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. Methods. A number of benchmark models have been created, covering low and high gas densities n = 10^3, 10^(5.5) cm^(-3) and far ultraviolet intensities χ = 10, 10^5 in units of the Draine field (FUV: 6 < h υ < 13.6 eV). The benchmark models were computed in two ways: one set assuming constant temperatures, thus testing the consistency of the chemical network and photo-processes, and a second set determining the temperature self consistently by solving the thermal balance, thus testing the modeling of the heating and cooling mechanisms accounting for the detailed energy balance throughout the clouds. Results. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes such as the effect of PAHs on the electron density or the temperature dependence of the dissociation of CO by cosmic ray induced secondary photons, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly. Nevertheless, the remaining spread in the computed observables such as the atomic fine-structure line intensities serves as a warning that there is still a considerable uncertainty when interpreting astronomical data with our models.

The ISO spectra of the bilobal planetary nebula Hb 5 are presented. These spectra are combined with the spectra in the visual wavelength region to obtain a complete, extinction corrected, spectrum. The chemical composition of the nebula is then calculated in several ways. First by directly calculating and adding individual ion abundances, assuming that all the ionic lines are formed in an ionized region surrounding the ionizing star. Secondly by building an "end-to-end model" nebula in which we have included a neutral region and a photodissociation region (PDR) beyond the ionized nebula. In this way we attempt to interpret the molecular hydrogen lines observed by ISO in a more self-consistent way. In the final analysis, the model is found to be basically heuristic, but gives new insights about the PDR and the PN. The implications of these are discussed.

We present a mid-infrared study of four H II regions based on the combination of Infrared Space Observatory (ISO) SWS observations with spatial information from SpectroCam-10 images and long-slit spectra. We report a prominent broad emission feature near 8.9 mum. The peak position and width of this feature do not resemble known polycyclic aromatic hydrocarbon ( PAH) or dust emission features. Its spatial distribution is analyzed for two H II regions and is clearly distinct from that of the PAH emission features, being more similar to that of the ionized gas and especially to the spatial distribution of the dust continuum at 15 mum. Thus, the carrier of this band likely resides within the H II region. Possible carriers of this feature are highlighted.

We present a mid-infrared study of four H II regions based on the combination of Infrared Space Observatory (ISO) SWS observations with spatial information from SpectroCam-10 images and long-slit spectra. We report a prominent broad emission feature near 8.9 mum. The peak position and width of this feature do not resemble known polycyclic aromatic hydrocarbon ( PAH) or dust emission features. Its spatial distribution is analyzed for two H II regions and is clearly distinct from that of the PAH emission features, being more similar to that of the ionized gas and especially to the spatial distribution of the dust continuum at 15 mum. Thus, the carrier of this band likely resides within the H II region. Possible carriers of this feature are highlighted.

We present near-IR speckle images of 21 massive Young Stellar Objects (YSOs) associated with outflows. The aim of this study is to search for sub-arcsecond reflection nebulae associated with the outflow cavity. We find that 6 of the massive YSOs show a conical nebula which can be interpreted in terms of reflected light from the dusty walls of the outflow cavity. In all cases, the small scale structures seen in our images are compared with outflow indicators found in the literature. No clear correlation is found between the presence of the reflection nebulosity and any property such as degree of embeddedness. We also note that 3 of the sources show close companions, one of them belonging also to the sample with conical nebula., Comment: 17 pages, 16 figures, accepted for publication in A&A

We present ground-based 2.1 to 4.2 microns observations of four low-mass protostars. We searched for the 4.1 microns OD stretch band, characteristic of solid HDO in grain mantles. We did not detect solid HDO in any of the four sources, but we derive 3-sigma upper limits from 0.5% to 2% for the HDO/H2O ratio depending on the source. These ratios provide strong constraints to solid-state deuteration models when compared to deuterium fractionation values observed in the gas phase. We discuss various scenarios that could lead to such a low water deuteration compared to the high formaldehyde and methanol deuteration observed in the gas phase., 8 pages, 6 figures Accepted for publication in A&A

While warm dense gas is prevalent around low-mass protostars, the presence of complex saturated molecules-the chemical inventory characteristic of hot cores-has remained elusive in such environments. Here we report the results of an IRAM 30 m study of the molecular composition associated with the low-mass protostar IRAS 16293-2422. Our observations highlight an extremely rich organic inventory in this source with abundant amounts of complex O- and N-bearing molecules such as formic acid, HCOOH, acetaldehyde, CH3CHO, methyl formate, CH3OCHO, dimethyl ether, CH3OCH3, acetic acid, CH3COOH, methyl cyanide, CH3CN, ethyl cyanide, C2H5CN, and propyne, CH3CCH. We compare the composition of the hot core around this low-mass young stellar object with those around massive protostars and address the chemical processes involved in molecular complexity in regions of star formation.

While warm dense gas is prevalent around low-mass protostars, the presence of complex saturated molecules-the chemical inventory characteristic of hot cores-has remained elusive in such environments. Here we report the results of an IRAM 30 m study of the molecular composition associated with the low-mass protostar IRAS 16293-2422. Our observations highlight an extremely rich organic inventory in this source with abundant amounts of complex O- and N-bearing molecules such as formic acid, HCOOH, acetaldehyde, CH3CHO, methyl formate, CH3OCHO, dimethyl ether, CH3OCH3, acetic acid, CH3COOH, methyl cyanide, CH3CN, ethyl cyanide, C2H5CN, and propyne, CH3CCH. We compare the composition of the hot core around this low-mass young stellar object with those around massive protostars and address the chemical processes involved in molecular complexity in regions of star formation.

While warm dense gas is prevalent around low-mass protostars, the presence of complex saturated molecules-the chemical inventory characteristic of hot cores-has remained elusive in such environments. Here we report the results of an IRAM 30 m study of the molecular composition associated with the low-mass protostar IRAS 16293-2422. Our observations highlight an extremely rich organic inventory in this source with abundant amounts of complex O- and N-bearing molecules such as formic acid, HCOOH, acetaldehyde, CH3CHO, methyl formate, CH3OCHO, dimethyl ether, CH3OCH3, acetic acid, CH3COOH, methyl cyanide, CH3CN, ethyl cyanide, C2H5CN, and propyne, CH3CCH. We compare the composition of the hot core around this low-mass young stellar object with those around massive protostars and address the chemical processes involved in molecular complexity in regions of star formation.

We aim to reveal the gas energetics in the circumstellar environment of the prototypical high-mass protostellar object AFGL2591 using space-based far-infrared observations of linear rotor molecules. Rotational spectral line signatures of CO, HCO+, CS, HCN and HNC from a 490-1240 GHz survey with Herschel/HIFI, complemented by ground-based JCMT and IRAM 30m spectra, cover transitions with E(up)/k between 5 and ~300 K (750K for 12C16O, using selected frequency settings up to 1850 GHz). The resolved spectral line profiles are used to separate and study various kinematic components. The line profiles show two emission components, the widest and bluest of which is attributed to an approaching outflow and the other to the envelope. We find evidence for progressively more redshifted and wider line profiles from the envelope gas with increasing energy level, qualitatively explained by residual outflow contribution picked up in the systematically decreasing beam size. Integrated line intensities for each species decrease as E(up)/k increases from 150 K for all species. This indicates that warm gas should be added to the model and that the model's geometry should provide low optical depth pathways for line emission from this warm gas to escape, for example in the form of UV heated outflow cavity walls viewed at a favorable inclination angle. Physical and chemical conditions derived for the outflow gas are similar to those in the protostellar envelope, possibly indicating that the modest velocity (, Accepted for publication in A&A; 16 pages, 8 figures, 4 tables

We present results of 1.3 mm dust polarization observations toward 16 nearby, low-mass protostars, mapped with ~2.5" resolution at CARMA. The results show that magnetic fields in protostellar cores on scales of ~1000 AU are not tightly aligned with outflows from the protostars. Rather, the data are consistent with scenarios where outflows and magnetic fields are preferentially misaligned (perpendicular), or where they are randomly aligned. If one assumes that outflows emerge along the rotation axes of circumstellar disks, and that the outflows have not disrupted the fields in the surrounding material, then our results imply that the disks are not aligned with the fields in the cores from which they formed., Comment: 6 pages, 2 figures, 1 table

We present the sub-millimeter spectra from 450 GHz to 1550 GHz of eleven nearby active galaxies observed with the SPIRE Fourier Transform Spectrometer (SPIRE/FTS) onboard Herschel. We detect CO transitions from J_up = 4 to 12, as well as the two [CI] fine structure lines at 492 and 809 GHz and the [NII] 461 GHz line. We used radiative transfer models to analyze the observed CO spectral line energy distributions (SLEDs). The FTS CO data were complemented with ground-based observations of the low-J CO lines. We found that the warm molecular gas traced by the mid-J CO transitions has similar physical conditions (n_H2 ~ 10^3.2 - 10^3.9 cm^-3 and T_kin ~ 300 - 800 K) in most of our galaxies. Furthermore, we found that this warm gas is likely producing the mid-IR rotational H2 emission. We could not determine the specific heating mechanism of the warm gas, however it is possibly related to the star-formation activity in these galaxies. Our modeling of the [CI] emission suggests that it is produced in cold (T_kin < 30 K) and dense (n_H2 > 10^3 cm^-3) molecular gas. Transitions of other molecules are often detected in our SPIRE/FTS spectra. The HF J=1-0 transition at 1232 GHz is detected in absorption in UGC05101 and in emission in NGC7130. In the latter, near-infrared pumping, chemical pumping, or collisional excitation with electrons are plausible excitation mechanisms likely related to the AGN of this galaxy. In some galaxies few H2O emission lines are present. Additionally, three OH+ lines at 909, 971, and 1033 GHz are identified in NGC7130., Accepted for publication in ApJ; 20 pages, 9 figures

The abundances of interstellar CH+ and SH+ are not well understood as their most likely formation channels are highly endothermic. Using data from Herschel, we study the formation of CH+ and SH+ in a typical high UV-illumination photon-dominated region (PDR), the Orion Bar. Herschel/HIFI provides velocity-resolved data of CH+ 1-0 and 2-1 and three hyperfine transitions of SH+. Herschel/PACS provides information on the excitation and spatial distribution of CH+ (up to J=6-5). The widths of the CH+ 2-1 and 1-0 transitions are of ~5 km/s, significantly broader than the typical width of dense gas tracers in the Orion Bar (2-3 km/s) and are comparable to the width of tracers of the interclump medium such as C+ and HF. The detected SH+ transitions are narrower compared to CH+ and have line widths of 3 km/s, indicating that SH+ emission mainly originates in denser condensations. Non-LTE radiative transfer models show that electron collisions affect the excitation of CH+ and SH+, and that reactive collisions need to be taken into account to calculate the excitation of CH+. Comparison to PDR models shows that CH+ and SH+ are tracers of the warm surface region (AV, Accepted for publication in A&A, 12 pages, 11 figures, abstract shortened

Based on astrochemical theory, the complex molecular composition around high-mass YSOs should evolve from the outer envelope in toward the central hot region as a sequence of temperature dependent chemical pathways are activated in ices and in the gas-phase. The resulting complex molecules have been broadly classified into three generations dependent on the temperature (25, and >100 K) required for formation. We combine IRAM 30m and Submillimeter Array observations to explore the spatial distribution of organic molecules around the high-mass young stellar object NGC 7538 IRS9, whose weak complex molecule emission previously escaped detection, quantifying the emission and abundance profiles of key organic molecules as a function of distance from the central protostar. We find that emission from N-bearing organics and saturated O-bearing organics present large increases in emission around 8000 AU and R, Accepted for publication in ApJ. 33 pages including 15 figures

This paper contains a summary of the results from the first years of observations with the HIFI instrument onboard ESA's Herschel space observatory. The paper starts by outlining the goals and possibilities of far-infrared and submillimeter astronomy, the limitations of the Earth's atmosphere, and the scientific scope of the Herschel-HIFI mission. The presentation of science results from the mission follows the life cycle of gas in galaxies as grouped into five themes: Structure of the interstellar medium, First steps in interstellar chemistry, Formation of stars and planets, Solar system results and Evolved stellar envelopes. The HIFI observations paint a picture where the interstellar medium in galaxies has a mixed, rather than a layered structure; the same conclusion may hold for protoplanetary disks. In addition, the HIFI data show that exchange of matter between comets and asteroids with planets and moons plays a large role. The paper concludes with an outlook to future instrumentation in the far-infrared and submillimeter wavelength ranges., Invited review paper to appear in Advances of Space Research; 32 pages, 9 figures. See http://www.sron.rug.nl/~vdtak/adspres.pdf for a 14-page journal-style version

The HDO/H2O ratio is a powerful diagnostic to understand the evolution of water from the first stages of star formation to the formation of planets and comets. Our aim is to determine precisely the abundance distribution of HDO towards the low-mass protostar IRAS16293-2422 and learn more about the water formation mechanisms by determining the HDO/H2O abundance ratio. A spectral survey of the source IRAS16293-2422 was carried out in the framework of the CHESS Herschel Key program with the HIFI instrument, allowing the detection of numerous HDO lines. Other transitions have been observed previously with ground-based telescopes. The spherical Monte Carlo radiative transfer code RATRAN was used to reproduce the observed line profiles of HDO by assuming an abundance jump. To determine the H2O abundance throughout the envelope, a similar study was made of the H2-18O observed lines, as the H2O main isotope lines are contaminated by the outflows. We derive an inner HDO abundance of 1.7e-7 and an outer HDO abundance of 8e-11. To reproduce the HDO absorption lines, it is necessary to add an absorbing layer in front of the envelope. It may correspond to a water-rich layer created by the photodesorption of the ices at the edges of the molecular cloud. The HDO/H2O ratio is ~1.4-5.8% in the hot corino whereas it is ~0.2-2.2% in the outer envelope. It is estimated at ~4.8% in the added absorbing layer. Although it is clearly higher than the cosmic D/H abundance, the HDO/H2O ratio remains lower than the D/H ratio derived for other deuterated molecules observed in the same source. The similarity of the ratios derived in the hot corino and in the added absorbing layer suggests that water formed before the gravitational collapse of the protostar, contrary to formaldehyde and methanol, which formed later once the CO molecules had depleted on the grains., Comment: 12 pages, 11 figures

We present Herschel/HIFI observations of fourteen water lines in W43-MM1, a massive protostellar object in the luminous star cluster-forming region W43. We analyze the gas dynamics from the line profiles using Herschel-HIFI observations (WISH-KP) of fourteen far-IR water lines (H2O, H217O, H218O), CS(11-10), and C18O(9-8) lines, and using our modeling of the continuum spectral energy distribution. As for lower mass protostellar objects, the molecular line profiles are a mix of emission and absorption, and can be decomposed into 'medium', and 'broad' velocity components. The broad component is the outflow associated with protostars of all masses. Our modeling shows that the remainder of the water profiles can be well fitted by an infalling and passively heated envelope, with highly supersonic turbulence varying from 2.2 km/s in the inner region to 3.5 km/s in the outer envelope. Also, W43-MM1 has a high accretion rate, between 4.0 x 10^{-4} and 4.0 x 10^{-2} \msun /yr, derived from the fast (0.4-2.9 km/s) infall observed. We estimate a lower mass limit of gaseous water of 0.11 \msun and total water luminosity of 1.5 \lsun (in the 14 lines presented here). The central hot core is detected with a water abundance of 1.4 x 10^{-4} while the water abundance for the outer envelope is 8 x10^{-8}. The latter value is higher than in other sources, most likely related to the high turbulence and the micro-shocks created by its dissipation. Examining water lines of various energies, we find that the turbulent velocity increases with the distance to the center. While not in clear disagreement with the competitive accretion scenario, this behavior is predicted by the turbulent core model. Moreover, the estimated accretion rate is high enough to overcome the expected radiation pressure., Accepted in A&A on April 2, 2012. 12 pages 7 figures