Showing 1–50 of 191 results for author: Lis, D

Bias versus variance when fitting multi-species molecular lines with a non-LTE radiative transfer model

Authors: Antoine Roueff, Jérôme Pety, Maryvonne Gerin, Léontine Ségal, Javier Goicoechea, Harvey Liszt, Pierre Gratier, Ivana Bešlić, Lucas Einig, M. Gaudel, Jan Orkisz, Pierre Palud, Miriam Santa-Maria, Victor de Souza Magalhaes, Antoine Zakardjian, Sebastien Bardeau, Emeric E. Bron, Pierre Chainais, Simon Coudé, Karine Demyk, Viviana Guzman Veloso, Annie Hughes, David Languignon, François Levrier, Dariusz C Lis , et al. (6 additional authors not shown)

Abstract: Robust radiative transfer techniques are requisite for efficiently extracting the physical and chemical information from molecular rotational lines.We study several hypotheses that enable robust estimations of the column densities and physical conditions when fitting one or two transitions per molecular species. We study the extent to which simplifying assumptions aimed at reducing the complexity… ▽ More Robust radiative transfer techniques are requisite for efficiently extracting the physical and chemical information from molecular rotational lines.We study several hypotheses that enable robust estimations of the column densities and physical conditions when fitting one or two transitions per molecular species. We study the extent to which simplifying assumptions aimed at reducing the complexity of the problem introduce estimation biases and how to detect them.We focus on the CO and HCO+ isotopologues and analyze maps of a 50 square arcminutes field. We used the RADEX escape probability model to solve the statistical equilibrium equations and compute the emerging line profiles, assuming that all species coexist. Depending on the considered set of species, we also fixed the abundance ratio between some species and explored different values. We proposed a maximum likelihood estimator to infer the physical conditions and considered the effect of both the thermal noise and calibration uncertainty. We analyzed any potential biases induced by model misspecifications by comparing the results on the actual data for several sets of species and confirmed with Monte Carlo simulations. The variance of the estimations and the efficiency of the estimator were studied based on the Cram{é}r-Rao lower bound.Column densities can be estimated with 30% accuracy, while the best estimations of the volume density are found to be within a factor of two. Under the chosen model framework, the peak 12CO(1--0) is useful for constraining the kinetic temperature. The thermal pressure is better and more robustly estimated than the volume density and kinetic temperature separately. Analyzing CO and HCO+ isotopologues and fitting the full line profile are recommended practices with respect to detecting possible biases.Combining a non-local thermodynamic equilibrium model with a rigorous analysis of the accuracy allows us to obtain an efficient estimator and identify where the model is misspecified. We note that other combinations of molecular lines could be studied in the future. △ Less

Submitted 29 March, 2024; originally announced March 2024.

Comments: Astronomy and Astrophysics - A\&A, In press

The magnetic field in the Flame nebula

Authors: Ivana Bešlić, Simon Coudé, Dariusz C. Lis, Maryvonne Gerin, Paul F. Goldsmith, Jerome Pety, Antoine Roueff, Karine Demyk, Charles D. Dowell, Lucas Einig, Javier R. Goicoechea, Francois Levrier, Jan Orkisz, Nicolas Peretto, Miriam G. Santa-Maria, Nathalie Ysard, Antoine Zakardjian

Abstract: Star formation is essential in galaxy evolution and the cycling of matter. The support of interstellar clouds against gravitational collapse by magnetic (B-) fields has been proposed to explain the low observed star formation efficiency in galaxies and the Milky Way. Despite the Planck satellite providing a 5-15' all-sky map of the B-field geometry in the diffuse interstellar medium, higher spatia… ▽ More Star formation is essential in galaxy evolution and the cycling of matter. The support of interstellar clouds against gravitational collapse by magnetic (B-) fields has been proposed to explain the low observed star formation efficiency in galaxies and the Milky Way. Despite the Planck satellite providing a 5-15' all-sky map of the B-field geometry in the diffuse interstellar medium, higher spatial resolution observations are required to understand the transition from diffuse gas to gravitationally unstable filaments. NGC 2024, the Flame Nebula, in the nearby Orion B molecular cloud, contains a young, expanding HII region and a dense filament that harbors embedded protostellar objects. Therefore, NGC 2024 is an excellent opportunity to study the role of B-fields in the formation, evolution, and collapse of filaments, as well as the dynamics and effects of young HII regions on the surrounding molecular gas. We combine new 154 and 216 micron dust polarization measurements carried out using the HAWC+ instrument aboard SOFIA with molecular line observations of 12CN(1-0) and HCO+(1-0) from the IRAM 30-meter telescope to determine the B-field geometry and to estimate the plane of the sky magnetic field strength across the NGC 2024. The HAWC+ observations show an ordered B-field geometry in NGC 2024 that follows the morphology of the expanding HII region and the direction of the main filament. The derived plane of the sky B-field strength is moderate, ranging from 30 to 80 micro G. The strongest B-field is found at the northern-west edge of the HII region, characterized by the highest gas densities and molecular line widths. In contrast, the weakest field is found toward the filament in NGC 2024. The B-field has a non-negligible influence on the gas stability at the edges of the expanding HII shell (gas impacted by the stellar feedback) and the filament (site of the current star formation). △ Less

Submitted 7 February, 2024; v1 submitted 30 January, 2024; originally announced January 2024.

Comments: 36 pages, 26 figures Accepted for publication in Astronomy & Astrophysics

PRIMA General Observer Science Book

Authors: A. Moullet, T. Kataria, D. Lis, S. Unwin, Y. Hasegawa, E. Mills, C. Battersby, A. Roc, M. Meixner

Abstract: PRIMA (The PRobe for-Infrared Mission for Astrophysics) is a concept for a far-infrared (IR) observatory. PRIMA features a cryogenically cooled 1.8 m diameter telescope and is designed to carry two science instruments enabling ultra-high sensitivity imaging and spectroscopic studies in the 24 to 235 microns wavelength range. The resulting observatory is a powerful survey and discovery machine, wit… ▽ More PRIMA (The PRobe for-Infrared Mission for Astrophysics) is a concept for a far-infrared (IR) observatory. PRIMA features a cryogenically cooled 1.8 m diameter telescope and is designed to carry two science instruments enabling ultra-high sensitivity imaging and spectroscopic studies in the 24 to 235 microns wavelength range. The resulting observatory is a powerful survey and discovery machine, with mapping speeds better by 2 - 4 orders of magnitude with respect to its far-IR predecessors. The bulk of the observing time on PRIMA should be made available to the community through a General Observer (GO) program offering 75% of the mission time over 5 years. In March 2023, the international astronomy community was encouraged to prepare authored contributions articulating scientific cases that are enabled by the telescope massive sensitivity advance and broad spectral coverage, and that could be performed within the context of GO program. This document, the PRIMA General Observer Science Book, is the edited collection of the 76 received contributions. △ Less

Submitted 31 October, 2023; originally announced October 2023.

Comments: A. Moullet, T. Kataria, D. Lis, S. Unwin, Y. Hasegawa, E. Mills, C. Battersby, A. Roc, M. Meixner are the editors of the PRIMA General Observer Science Book. The book compiles 76 authored contributions. 399 pages

Retrievals of Protoplanetary Disk Parameters using Thermochemical Models: I. Disk Gas Mass from Hydrogen Deuteride Spectroscopy

Authors: Young Min Seo, Karen Willacy, Geoffrey Bryden, Dariusz C. Lis, Paul F. Goldsmith, Klaus M. Pontoppidan, Wing-Fai Thi

Abstract: We discuss statistical relationships between the mass of protoplanetary disks and the hydrogen deuteride (HD) line emission and the dust spectral energy distribution (SED) determined using 3000 ProDiMo disk models. The models have 15 free parameters describing disk physical properties, the central star, and the local radiation field. The sampling of physical parameters is done using a Monte Carlo… ▽ More We discuss statistical relationships between the mass of protoplanetary disks and the hydrogen deuteride (HD) line emission and the dust spectral energy distribution (SED) determined using 3000 ProDiMo disk models. The models have 15 free parameters describing disk physical properties, the central star, and the local radiation field. The sampling of physical parameters is done using a Monte Carlo approach to evaluate the probability density functions of observables as a function of physical parameters. We find that the HD fractional abundance is almost constant even though the UV flux varies by several orders of magnitude. Probing the statistical relation between the physical quantities and the HD flux, we find that low-mass (optically thin) disks display a tight correlation between the average disk gas temperature and HD line flux, while massive disks show no such correlation. We demonstrate that the central star luminosity, disk size, dust size distribution, and HD flux may be used to determine the disk gas mass to within a factor of three. We also find that the far-IR and sub-mm/mm SEDs and the HD flux may serve as strong constraints for determining the disk gas mass to within a factor of two. If the HD lines are fully spectrally resolved ($R\gtrsim 1.5\times10^6, Δv=0.2~\rm km\,s^{-1}$), the 56 $μ$m and 112 $μ$m HD line profiles alone may constrain the disk gas mass to within a factor of two. △ Less

Submitted 24 April, 2024; v1 submitted 24 October, 2023; originally announced October 2023.

Comments: 31 pages, 16 figures, ApJ accepted, preprint

Mapping Physical Conditions in Neighboring Hot Cores: NOEMA Studies of W3(H$_2$O) and W3(OH)

Authors: Morgan M. Giese, Will E. Thompson, Dariusz C. Lis, Susanna L. Widicus Weaver

Abstract: The complex chemistry that occurs in star-forming regions can provide insight into the formation of prebiotic molecules at various evolutionary stages of star formation. To study this process, we present millimeter-wave interferometric observations of the neighboring hot cores W3(H$_2$O) and W3(OH) carried out using the NOEMA interferometer. We have analyzed distributions of six molecules that acc… ▽ More The complex chemistry that occurs in star-forming regions can provide insight into the formation of prebiotic molecules at various evolutionary stages of star formation. To study this process, we present millimeter-wave interferometric observations of the neighboring hot cores W3(H$_2$O) and W3(OH) carried out using the NOEMA interferometer. We have analyzed distributions of six molecules that account for most observed lines across both cores and have constructed physical parameter maps for rotational temperature, column density, and velocity field with corresponding uncertainties. We discuss the derived spatial distributions of these parameters in the context of the physical structure of the source. We propose the use of HCOOCH$_3$ as a new temperature tracer in W3(H$_2$O) and W3(OH) in addition to the more commonly used CH$_3$CN. By analyzing the physically-derived parameters for each molecule across both W3(H$_2$O) and W3(OH), the work presented herein further demonstrates the impact of physical environment on hot cores at different evolutionary stages. △ Less

Submitted 16 November, 2023; v1 submitted 18 October, 2023; originally announced October 2023.

Comments: Accepted to The Astrophysical Journal

A JWST inventory of protoplanetary disk ices: The edge-on protoplanetary disk HH 48 NE, seen with the Ice Age ERS program

Authors: J. A. Sturm, M. K. McClure, T. L. Beck, D. Harsono, J. B. Bergner, E. Dartois, A. C. A. Boogert, J. E. Chiar, M. A. Cordiner, M. N. Drozdovskaya, S. Ioppolo, C. J. Law, H. Linnartz, D. C. Lis, G. J. Melnick, B. A. McGuire, J. A. Noble, K. I. Öberg, M. E. Palumbo, Y. J. Pendleton, G. Perotti, K. M. Pontoppidan, D. Qasim, W. R. M. Rocha, H. Terada , et al. (2 additional authors not shown)

Abstract: Ices are the main carriers of volatiles in protoplanetary disks and are crucial to our understanding of the chemistry that ultimately sets the organic composition of planets. The ERS program Ice Age on the JWST follows the ice evolution through all stages of star and planet formation. JWST/NIRSpec observations of the edge-on Class II protoplanetary disk HH~48~NE reveal spatially resolved absorptio… ▽ More Ices are the main carriers of volatiles in protoplanetary disks and are crucial to our understanding of the chemistry that ultimately sets the organic composition of planets. The ERS program Ice Age on the JWST follows the ice evolution through all stages of star and planet formation. JWST/NIRSpec observations of the edge-on Class II protoplanetary disk HH~48~NE reveal spatially resolved absorption features of the major ice components H$_2$O, CO$_2$, CO, and multiple weaker signatures from less abundant ices NH$_3$, OCN$^-$, and OCS. Isotopologue $^{13}$CO$_2$ ice has been detected for the first time in a protoplanetary disk. Since multiple complex light paths contribute to the observed flux, the ice absorption features are filled in by ice-free scattered light. The $^{12}$CO$_2$/$^{13}$CO$_2$ ratio of 14 implies that the $^{12}$CO$_2$ feature is saturated, without the flux approaching 0, indicative of a very high CO$_2$ column density on the line of sight, and a corresponding abundance with respect to hydrogen that is higher than ISM values by a factor of at least a few. Observations of rare isotopologues are crucial, as we show that the $^{13}$CO$_2$ observation allows us to determine the column density of CO$_2$ to be at an order of magnitude higher than the lower limit directly inferred from the observed optical depth. Radial variations in ice abundance, e.g., snowlines, are significantly modified since all observed photons have passed through the full radial extent of the disk. CO ice is observed at perplexing heights in the disk, extending to the top of the CO-emitting gas layer. We argue that the most likely interpretation is that we observe some CO ice at high temperatures, trapped in less volatile ices like H$_2$O and CO$_2$. Future radiative transfer models will be required to constrain the implications on our current understanding of disk physics and chemistry. △ Less

Submitted 14 September, 2023; originally announced September 2023.

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

HCN emission from translucent gas and UV-illuminated cloud edges revealed by wide-field IRAM 30m maps of Orion B GMC: Revisiting its role as tracer of the dense gas reservoir for star formation

Authors: M. G. Santa-Maria, J. R. Goicoechea, J. Pety, M. Gerin, J. H. Orkisz, F. Le Petit, L. Einig, P. Palud, V. de Souza Magalhaes, I. Bešlić, L. Segal, S. Bardeau, E. Bron, P. Chainais, J. Chanussot, P. Gratier, V. V. Guzmán, A. Hughes, D. Languignon, F. Levrier, D. C. Lis, H. S. Liszt, J. Le Bourlot, Y. Oya, K. Öberg , et al. (6 additional authors not shown)

Abstract: We present 5 deg^2 (~250 pc^2) HCN, HNC, HCO+, and CO J=1-0 maps of the Orion B GMC, complemented with existing wide-field [CI] 492 GHz maps, as well as new pointed observations of rotationally excited HCN, HNC, H13CN, and HN13C lines. We detect anomalous HCN J=1-0 hyperfine structure line emission almost everywhere in the cloud. About 70% of the total HCN J=1-0 luminosity arises from gas at A_V <… ▽ More We present 5 deg^2 (~250 pc^2) HCN, HNC, HCO+, and CO J=1-0 maps of the Orion B GMC, complemented with existing wide-field [CI] 492 GHz maps, as well as new pointed observations of rotationally excited HCN, HNC, H13CN, and HN13C lines. We detect anomalous HCN J=1-0 hyperfine structure line emission almost everywhere in the cloud. About 70% of the total HCN J=1-0 luminosity arises from gas at A_V < 8 mag. The HCN/CO J=1-0 line intensity ratio shows a bimodal behavior with an inflection point at A_V < 3 mag typical of translucent gas and UV-illuminated cloud edges. We find that most of the HCN J=1-0 emission arises from extended gas with n(H2) ~< 10^4 cm^-3, even lower density gas if the ionization fraction is > 10^-5 and electron excitation dominates. This result explains the low-A_V branch of the HCN/CO J=1-0 intensity ratio distribution. Indeed, the highest HCN/CO ratios (~0.1) at A_V < 3 mag correspond to regions of high [CI] 492 GHz/CO J=1-0 intensity ratios (>1) characteristic of low-density PDRs. Enhanced FUV radiation favors the formation and excitation of HCN on large scales, not only in dense star-forming clumps. The low surface brightness HCN and HCO+ J=1-0 emission scale with I_FIR (a proxy of the stellar FUV radiation field) in a similar way. Together with CO J=1-0, these lines respond to increasing I_FIR up to G0~20. On the other hand, the bright HCN J=1-0 emission from dense gas in star-forming clumps weakly responds to I_FIR once the FUV radiation field becomes too intense (G0>1500). The different power law scalings (produced by different chemistries, densities, and line excitation regimes) in a single but spatially resolved GMC resemble the variety of Kennicutt-Schmidt law indexes found in galaxy averages. As a corollary for extragalactic studies, we conclude that high HCN/CO J=1-0 line intensity ratios do not always imply the presence of dense gas. △ Less

Submitted 18 September, 2023; v1 submitted 6 September, 2023; originally announced September 2023.

Comments: accepted for publication in A&A. 24 pages, 18 figures, plus Appendix. Abridged Abstract

Journal ref: A&A 679, A4 (2023)

Deep learning denoising by dimension reduction: Application to the ORION-B line cubes

Authors: Lucas Einig, Jérôme Pety, Antoine Roueff, Paul Vandame, Jocelyn Chanussot, Maryvonne Gerin, Jan H. Orkisz, Pierre Palud, Miriam Garcia Santa-Maria, Victor de Souza Magalhaes, Ivana Bešlić, Sébastien Bardeau, Emeric E. Bron, Pierre Chainais, Javier R Goicoechea, Pierre Gratier, Viviana Guzman Veloso, Annie Hughes, Jouni Kainulainen, David Languignon, Rosine Lallement, François Levrier, Dariuscz C. Lis, Harvey Liszt, Jacques Le Bourlot , et al. (7 additional authors not shown)

Abstract: Context. The availability of large bandwidth receivers for millimeter radio telescopes allows the acquisition of position-position-frequency data cubes over a wide field of view and a broad frequency coverage. These cubes contain much information on the physical, chemical, and kinematical properties of the emitting gas. However, their large size coupled with inhomogenous signal-to-noise ratio (SNR… ▽ More Context. The availability of large bandwidth receivers for millimeter radio telescopes allows the acquisition of position-position-frequency data cubes over a wide field of view and a broad frequency coverage. These cubes contain much information on the physical, chemical, and kinematical properties of the emitting gas. However, their large size coupled with inhomogenous signal-to-noise ratio (SNR) are major challenges for consistent analysis and interpretation.Aims. We search for a denoising method of the low SNR regions of the studied data cubes that would allow to recover the low SNR emission without distorting the signals with high SNR.Methods. We perform an in-depth data analysis of the 13 CO and C 17 O (1 -- 0) data cubes obtained as part of the ORION-B large program performed at the IRAM 30m telescope. We analyse the statistical properties of the noise and the evolution of the correlation of the signal in a given frequency channel with that of the adjacent channels. This allows us to propose significant improvements of typical autoassociative neural networks, often used to denoise hyperspectral Earth remote sensing data. Applying this method to the 13 CO (1 -- 0) cube, we compare the denoised data with those derived with the multiple Gaussian fitting algorithm ROHSA, considered as the state of the art procedure for data line cubes.Results. The nature of astronomical spectral data cubes is distinct from that of the hyperspectral data usually studied in the Earth remote sensing literature because the observed intensities become statistically independent beyond a short channel separation. This lack of redundancy in data has led us to adapt the method, notably by taking into account the sparsity of the signal along the spectral axis. The application of the proposed algorithm leads to an increase of the SNR in voxels with weak signal, while preserving the spectral shape of the data in high SNR voxels.Conclusions. The proposed algorithm that combines a detailed analysis of the noise statistics with an innovative autoencoder architecture is a promising path to denoise radio-astronomy line data cubes. In the future, exploring whether a better use of the spatial correlations of the noise may further improve the denoising performances seems a promising avenue. In addition, △ Less

Submitted 24 July, 2023; originally announced July 2023.

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

Comparing Complex Chemistry in Neighboring Hot Cores: NOEMA Studies of W3(H$_{2}$O) and W3(OH)

Authors: Will E. Thompson, Morgan M. Giese, Dariusz C. Lis, Susanna L. Widicus Weaver

Abstract: Presented here are NOEMA interferometric observations of the neighboring hot cores W3(H$_{2}$O) and W3(OH). The presence of two star-forming cores at different evolutionary stages within the same parent cloud presents a unique opportunity to study how the physics of the source and its evolutionary stage impact the chemistry. Through spectral analysis and imaging, we identify over twenty molecules… ▽ More Presented here are NOEMA interferometric observations of the neighboring hot cores W3(H$_{2}$O) and W3(OH). The presence of two star-forming cores at different evolutionary stages within the same parent cloud presents a unique opportunity to study how the physics of the source and its evolutionary stage impact the chemistry. Through spectral analysis and imaging, we identify over twenty molecules in these cores. Most notably, we have detected HDO and CH$_{3}$CH$_{2}$CN in W3(OH), which were previously not detected in this core. We have imaged the molecular emission, revealing new structural features within these sources. W3(OH) shows absorption in a "dusty cocoon" surrounded by molecular emission. These observations also reveal extended emission that is potentially indicative of a low-velocity shock. From the information obtained herein, we have constructed column density and temperature maps for methanol and compared this information to the molecular images. By comparing the spatial distribution of molecules which may be destroyed at later stages of star formation, this work demonstrates the impact of physical environment on chemistry in star-forming regions at different evolutionary stages. △ Less

Submitted 18 July, 2023; originally announced July 2023.

Comments: Published in The Astrophysical Journal

Journal ref: ApJ 952 50 (2023)

Gas Sources from the Coma and Nucleus of Comet 46P/Wirtanen Observed Using ALMA

Authors: M. A. Cordiner, N. X. Roth, S. N. Milam, G. Villanueva, D. Bockelee-Morvan, A. J. Remijan, S. B. Charnley, N. Biver, D. C. Lis, C. Qi, B. Bonev, J. Crovisier, J. Boissier

Abstract: Gas-phase molecules in cometary atmospheres (comae) originate primarily from (1) outgassing by the nucleus, (2) sublimation of icy grains in the near-nucleus coma, and (3) coma (photo-)chemical processes. However, the majority of cometary gases observed at radio wavelengths have yet to be mapped, so their production/release mechanisms remain uncertain. Here we present observations of six molecular… ▽ More Gas-phase molecules in cometary atmospheres (comae) originate primarily from (1) outgassing by the nucleus, (2) sublimation of icy grains in the near-nucleus coma, and (3) coma (photo-)chemical processes. However, the majority of cometary gases observed at radio wavelengths have yet to be mapped, so their production/release mechanisms remain uncertain. Here we present observations of six molecular species towards comet 46P/Wirtanen, obtained using the Atacama Large Millimeter/submillimeter Array (ALMA) during the comet's unusually close (~0.1 au) approach to Earth in December 2018. Interferometric maps of HCN, CH3OH, CH3CN, H2CO, CS and HNC were obtained at an unprecedented sky-projected spatial resolution of up to 25 km, enabling the nucleus and coma sources of these molecules to be accurately quantified. The HCN, CH3OH and CH3CN spatial distributions are consistent with production by direct outgassing from (or very near to) the nucleus, with a significant proportion of the observed CH3OH originating from sublimation of icy grains in the near-nucleus coma (at a scale-length $L_p=36\pm7$ km). On the other hand, H2CO, CS and HNC originate primarily from distributed coma sources (with $L_p$ values in the range 550-16,000 km), the identities of which remain to be established. The HCN, CH3OH and HNC abundances in 46P are consistent with the average values previously observed in comets, whereas the H2CO, CH3CN and CS abundances are relatively low. △ Less

Submitted 19 June, 2023; v1 submitted 8 May, 2023; originally announced May 2023.

Comments: Accepted for publication in the Astrophysical Journal on 2023-06-17

Molecular Outgassing in Centaur 29P/Schwassmann-Wachmann 1 During Its Exceptional 2021 Outburst: Coordinated Multi-Wavelength Observations Using nFLASH at APEX and iSHELL at the NASA-IRTF

Authors: Nathan X. Roth, Stefanie N. Milam, Michael A. DiSanti, Geronimo L. Villanueva, Sara Faggi, Boncho P. Bonev, Martin A. Cordiner, Anthony J. Remijan, Dominique Bockelée-Morvan, Nicolas Biver, Jacques Crovisier, Dariusz C. Lis, Steven B. Charnley, Emmanuel Jehin, Eva. S. Wirström, Adam J. McKay

Abstract: The extraordinary 2021 September-October outburst of Centaur 29P/Schwassmann-Wachmann 1 afforded an opportunity to test the composition of primitive Kuiper disk material at high sensitivity. We conducted nearly simultaneous multi-wavelength spectroscopic observations of 29P/Schwassmann-Wachmann 1 using iSHELL at the NASA Infrared Telescope Facility and nFLASH at the Atacama Pathfinder EXperiment (… ▽ More The extraordinary 2021 September-October outburst of Centaur 29P/Schwassmann-Wachmann 1 afforded an opportunity to test the composition of primitive Kuiper disk material at high sensitivity. We conducted nearly simultaneous multi-wavelength spectroscopic observations of 29P/Schwassmann-Wachmann 1 using iSHELL at the NASA Infrared Telescope Facility and nFLASH at the Atacama Pathfinder EXperiment (APEX) on 2021 October 6, with follow-up APEX/nFLASH observations on 2021 October 7 and 2022 April 3. This coordinated campaign between near-infrared and radio wavelengths enabled us to sample molecular emission from a wealth of coma molecules and to perform measurements that cannot be accomplished with either wavelength alone. We securely detected CO emission on all dates with both facilities, including velocity-resolved spectra of the CO (J=2-1) transition with APEX/nFLASH and multiple CO (v=1-0) rovibrational transitions with IRTF/iSHELL. We report rotational temperatures, coma kinematics, and production rates for CO and stringent (3-sigma) upper limits on abundance ratios relative to CO for CH4, C2H6, CH3OH, H2CO, CS, and OCS. Our upper limits for CS/CO and OCS/CO represent their first values in the literature for this Centaur. Upper limits for CH4, C2H6, CH3OH, and H2CO are the most stringent reported to date, and are most similar to values found in ultra CO-rich Oort cloud comet C/2016 R2 (PanSTARRS), which may have implications for how ices are preserved in cometary nuclei. We demonstrate the superb synergy of coordinated radio and near-infrared measurements, and advocate for future small body studies that jointly leverage the capabilities of each wavelength. △ Less

Submitted 27 April, 2023; originally announced April 2023.

An Ice Age JWST inventory of dense molecular cloud ices

Authors: M. K. McClure, W. R. M. Rocha, K. M. Pontoppidan, N. Crouzet, L. E. U. Chu, E. Dartois, T. Lamberts, J. A. Noble, Y. J. Pendleton, G. Perotti, D. Qasim, M. G. Rachid, Z. L. Smith, Fengwu Sun, Tracy L Beck, A. C. A. Boogert, W. A. Brown, P. Caselli, S. B. Charnley, Herma M. Cuppen, H. Dickinson, M. N. Drozdovskaya, E. Egami, J. Erkal, H. Fraser , et al. (17 additional authors not shown)

Abstract: Icy grain mantles are the main reservoir of the volatile elements that link chemical processes in dark, interstellar clouds with the formation of planets and composition of their atmospheres. The initial ice composition is set in the cold, dense parts of molecular clouds, prior to the onset of star formation. With the exquisite sensitivity of JWST, this critical stage of ice evolution is now acces… ▽ More Icy grain mantles are the main reservoir of the volatile elements that link chemical processes in dark, interstellar clouds with the formation of planets and composition of their atmospheres. The initial ice composition is set in the cold, dense parts of molecular clouds, prior to the onset of star formation. With the exquisite sensitivity of JWST, this critical stage of ice evolution is now accessible for detailed study. Here we show the first results of the Early Release Science program "Ice Age" that reveal the rich composition of these dense cloud ices. Weak ices, including, $^{13}$CO$_2$, OCN$^-$, $^{13}$CO, OCS, and COMs functional groups are now detected along two pre-stellar lines of sight. The $^{12}$CO$_2$ ice profile indicates modest growth of the icy grains. Column densities of the major and minor ice species indicate that ices contribute between 2 and 19% of the bulk budgets of the key C, O, N, and S elements. Our results suggest that the formation of simple and complex molecules could begin early in a water-ice rich environment. △ Less

Submitted 22 January, 2023; originally announced January 2023.

Comments: To appear in Nature Astronomy on January 23rd, 2023. 33 pages, 16 figures, 3 tables; includes extended and supplemental data sections. Part of the JWST Ice Age Early Release Science program's science enabling products. Enhanced spectra downloadable on Zenodo at the following DOI: 10.5281/zenodo.7501239

Atomic oxygen abundance toward Sagittarius B2

Authors: Dariusz C. Lis, Paul F. Goldsmith, Rolf Güsten, Peter Schilke, Helmut Wiesemeyer, Youngmin Seo, Michael W. Werner

Abstract: A substantial fraction of oxygen in diffuse clouds is unaccounted for by observations and is postulated to be in an unknown refractory form, referred to as unidentified depleted oxygen (UDO), which, depending on the local gas density, may contribute up to 50% of the total oxygen content. Previous Infrared Space Observatory (ISO) observations suggest that a significant fraction of oxygen in even de… ▽ More A substantial fraction of oxygen in diffuse clouds is unaccounted for by observations and is postulated to be in an unknown refractory form, referred to as unidentified depleted oxygen (UDO), which, depending on the local gas density, may contribute up to 50% of the total oxygen content. Previous Infrared Space Observatory (ISO) observations suggest that a significant fraction of oxygen in even denser, translucent clouds may be in atomic form. We have analyzed velocity-resolved archival SOFIA observations of the 63 $μ$m fine-structure [O I] transition toward the high-mass star-forming region Sgr B2(M) in the Central Molecular Zone. The foreground spiral-arm clouds as well as the extended Sgr B2 envelope between the Sun and the background dust continuum source produce multiple [O i] absorption components, spectrally separated in velocity space. The gas-phase atomic oxygen column density in foreground clouds toward Sgr B2 is well correlated with the total hydrogen column density, with an average atomic oxygen abundance of $(2.51 \pm 0.69) \times 10^{-4}$ with respect to hydrogen nuclei. This value is in good agreement with the earlier ISO measurements on the same line of sight, and is about 35% lower than the total interstellar medium oxygen abundance in the low-density warm gas, as measured in the UV. We find no evidence that a significant fraction of the oxygen on the line of sight toward Sagittarius B2 is in the form of UDO. △ Less

Submitted 9 January, 2023; originally announced January 2023.

Comments: 7 pages, 7 figures, and 2 appendices. Astronomy & Astrophysics Letters, in press

Journal ref: A&A 669, L15 (2023)

Water, hydrogen cyanide, carbon monoxide, and dust production from distant comet 29P/Schwassmann-Wachmann 1

Authors: D. Bockelée-Morvan, N. Biver, C. A. Schambeau, J. Crovisier, C. Opitom, M. de Val Borro, E. Lellouch, P. Hartogh, B. Vandenbussche, E. Jehin, M. Kidger, M. Küppers, D. C. Lis, R. Moreno, S. Szutowicz, V. Zakharov

Abstract: 29P/Schwassmann-Wachmann 1 is a distant Centaur/comet, showing persistent CO-driven activity and frequent outbursts. We used the Herschel space observatory in 2010, 2011, and 2013 to observe H$_2$O and NH$_3$ and to image the dust coma. Observations with the IRAM 30 m were undertaken in 2007, 2010, 2011, and 2021 to monitor the CO production rate and to search for HCN. Modeling was performed to co… ▽ More 29P/Schwassmann-Wachmann 1 is a distant Centaur/comet, showing persistent CO-driven activity and frequent outbursts. We used the Herschel space observatory in 2010, 2011, and 2013 to observe H$_2$O and NH$_3$ and to image the dust coma. Observations with the IRAM 30 m were undertaken in 2007, 2010, 2011, and 2021 to monitor the CO production rate and to search for HCN. Modeling was performed to constrain the size of the sublimating icy grains and to derive the dust production rate. HCN is detected for the first time in comet 29P (at 5$σ$ in the line area). H$_2$O is detected as well, but not NH$_3$. H$_2$O and HCN line shapes differ strongly from the CO line shape, indicating that these two species are released from icy grains. CO production rates are in the range (2.9-5.6) $\times$ 10$^{28}$ s$^{-1}$ (1400--2600 kg s$^{-1}$). A correlation between the CO production rate and coma brightness is observed, as is a correlation between CO and H$_2$O production. The correlation obtained between the excess of CO production and excess of dust brightness with respect to the quiescent state is similar to that established for the continuous activity of comet Hale-Bopp. The measured $Q$(H$_2$O)/$Q$(CO) and $Q$(HCN)/$Q$(CO) production rate ratios are 10.0 $\pm$ 1.5 % and 0.12 $\pm$ 0.03 %, respectively, averaging the April-May 2010 measurements ($Q$(H$_2$O) = (4.1 $\pm$ 0.6) $\times$ 10$^{27}$ s$^{-1}$, $Q$(HCN) = (4.8 $\pm$ 1.1) $\times$ 10$^{25}$ s$^{-1}$). We derive three independent and similar values of the effective radius of the nucleus, $\sim$ 31 $\pm$ 3 km. The inferred dust mass-loss rates during quiescent phases are in the range 30-120 kg s$^{-1}$, indicating a dust-to-gas mass ratio $<$ 0.1 during quiescent activity. We conclude that strong local heterogeneities exist on the surface of 29P, with quenched dust activity from most of the surface, but not in outbursting regions. △ Less

Submitted 24 May, 2022; v1 submitted 23 May, 2022; originally announced May 2022.

Comments: 28 pages, 15 figures, A&A in press

HyGAL: Characterizing the Galactic ISM with observations of hydrides and other small molecules -- I. Survey description and a first look toward W3(OH), W3 IRS5 and NGC 7538 IRS1

Authors: A. M. Jacob, D. A. Neufeld, P. Schilke, H. Wiesemeyer, W. Kim, S. Bialy, M. Busch, D. Elia, E. Falgarone, M. Gerin, B. Godard, R. Higgins, P. Hennebelle, N. Indriolo, D. C. Lis, K. M. Menten, A. Sanchez-Monge, V. Ossenkopf-Okada, M. R. Rugel, D. Seifried, P. Sonnentrucker, S. Walch, M. Wolfire, F. Wyrowski, V. Valdivia

Abstract: The HyGAL SOFIA legacy program surveys six hydride molecules -- ArH+, OH+, H2O+, SH, OH, and CH -- and two atomic constituents -- C+ and O -- within the diffuse interstellar medium (ISM) by means of absorption-line spectroscopy toward 25 bright Galactic background continuum sources. This detailed spectroscopic study is designed to exploit the unique value of specific hydrides as tracers and probes… ▽ More The HyGAL SOFIA legacy program surveys six hydride molecules -- ArH+, OH+, H2O+, SH, OH, and CH -- and two atomic constituents -- C+ and O -- within the diffuse interstellar medium (ISM) by means of absorption-line spectroscopy toward 25 bright Galactic background continuum sources. This detailed spectroscopic study is designed to exploit the unique value of specific hydrides as tracers and probes of different phases of the ISM, as demonstrated by recent studies with the Herschel Space Observatory. The observations performed under the HyGAL program will allow us to address several questions related to the lifecycle of molecular material in the ISM and the physical processes that impact its phase transition, such as: (1) What is the distribution function of the H2 fraction in the ISM? (2) How does the ionization rate due to low-energy cosmic-rays vary within the Galaxy? (3) What is the nature of interstellar turbulence, and what mechanisms lead to its dissipation? This overview discusses the observing strategy, synergies with ancillary and archival observations, the data reduction and analysis schemes adopted; and presents the first results obtained toward three of the survey targets, W3(OH), W3IRS5 and NGC7538IRS1. Robust measurements of the column densities of these hydrides -- obtained through widespread observations of absorption lines-- help address the questions raised, and there is a timely synergy between these observations and the development of theoretical models, particularly pertaining to the formation of H2 within the turbulent ISM. The provision of enhanced HyGAL data products will therefore serve as a legacy for future ISM studies. △ Less

Submitted 10 February, 2022; originally announced February 2022.

Comments: 33 pages, 17 figures, 7 tables. Accepted for publication in the Astrophysical Journal

Chemical compositions in the vicinity of protostars in Ophiuchus

Authors: Kotomi Taniguchi, Liton Majumdar, Adele Plunkett, Shigehisa Takakuwa, Dariusz C. Lis, Paul F. Goldsmith, Fumitaka Nakamura, Masao Saito, Eric Herbst

Abstract: We have analyzed Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 4 Band 6 data toward two young stellar objects (YSOs), Oph-emb5 and Oph-emb9, in the Ophiuchus star-forming region. The YSO Oph-emb5 is located in a relatively quiescent region, whereas Oph-emb9 is irradiated by a nearby bright Herbig Be star. Molecular lines from $cyclic$-C$_{3}$H$_{2}$ ($c$-C$_{3}$H$_{2}$), H$_{2}$CO, CH… ▽ More We have analyzed Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 4 Band 6 data toward two young stellar objects (YSOs), Oph-emb5 and Oph-emb9, in the Ophiuchus star-forming region. The YSO Oph-emb5 is located in a relatively quiescent region, whereas Oph-emb9 is irradiated by a nearby bright Herbig Be star. Molecular lines from $cyclic$-C$_{3}$H$_{2}$ ($c$-C$_{3}$H$_{2}$), H$_{2}$CO, CH$_{3}$OH, $^{13}$CO, C$^{18}$O, and DCO$^{+}$ have been detected from both sources, while DCN is detected only in Oph-emb9. Around Oph-emb5, $c$-C$_{3}$H$_{2}$ is enhanced at the west side, relative to the IR source, whereas H$_{2}$CO and CH$_{3}$OH are abundant at the east side. In the field of Oph-emb9, moment 0 maps of the $c$-C$_{3}$H$_{2}$ lines show a peak at the eastern edge of the field of view, which is irradiated by the Herbig Be star. Moment 0 maps of CH$_{3}$OH and H$_{2}$CO show peaks farther from the bright star. We derive the $N$($c$-C$_{3}$H$_{2}$)/$N$(CH$_{3}$OH) column density ratios at the peak positions of $c$-C$_{3}$H$_{2}$ and CH$_{3}$OH near each YSO, which are identified based on their moment 0 maps. The $N$($c$-C$_{3}$H$_{2}$)/$N$(CH$_{3}$OH) ratio at the $c$-C$_{3}$H$_{2}$ peak is significantly higher than at the CH$_{3}$OH peak by a factor of $\sim 19$ in Oph-emb9, while the difference in this column density ratio between these two positions is a factor of $\sim2.6 $ in Oph-emb5. These differences are attributed to the efficiency of the photon-dominated region (PDR) chemistry in Oph-emb9. The higher DCO$^{+}$ column density and the detection of DCN in Oph-emb9 are also discussed in the context of UV irradiation flux. △ Less

Submitted 26 October, 2021; v1 submitted 24 August, 2021; originally announced August 2021.

Comments: 20 pages, 14 figures, 4 tables, Accepted by The Astrophysical Journal. We changed continuum images (Figure 2). We show the continuum images with only 12-m array data in this version. Since there are unknown issues when we combined 12-m and 7-m array data for the continuum image, the scales were wrong in the previous versions

Herschel observations of extraordinary sources: Full Herschel/HIFI molecular line survey of Sagittarius B2(M)

Authors: T. Möller, P. Schilke, A. Schmiedeke, E. A. Bergin, D. C. Lis, Á. Sánchez-Monge, A. Schwörer, C. Comito

Abstract: We present a full analysis of a broadband spectral line survey of Sagittarius B2 (Main), one of the most chemically rich regions in the Galaxy located within the giant molecular cloud complex Sgr B2 in the Central Molecular Zone. Our goal is to derive the molecular abundances and temperatures of the high-mass star-forming region Sgr B2(M) and thus its physical and astrochemical conditions. Sgr B2(… ▽ More We present a full analysis of a broadband spectral line survey of Sagittarius B2 (Main), one of the most chemically rich regions in the Galaxy located within the giant molecular cloud complex Sgr B2 in the Central Molecular Zone. Our goal is to derive the molecular abundances and temperatures of the high-mass star-forming region Sgr B2(M) and thus its physical and astrochemical conditions. Sgr B2(M) was observed using the Heterodyne Instrument for the Far-Infrared (HIFI) on board the Herschel Space Observatory in a spectral line survey from 480 to 1907 GHz at a spectral resolution of 1.1 MHz, which provides one of the largest spectral coverages ever obtained toward this high-mass star-forming region in the submillimeter with high spectral resolution and includes frequencies > 1 THz unobservable from the ground. We model the molecular emission from the submillimeter to the far-IR using the XCLASS program. For each molecule, a quantitative description was determined taking all emission and absorption features of that species across the entire spectral range into account. Additionally, we derive velocity resolved ortho / para ratios for those molecules for which ortho and para resolved molecular parameters are available. Finally, the temperature and velocity distributions are analyzed and the derived abundances are compared with those obtained for Sgr B2(N) from a similar HIFI survey. A total of 92 isotopologues were identified, arising from 49 different molecules, ranging from free ions to complex organic compounds and originating from a variety of environments from the cold envelope to hot and dense gas within the cores. Sulfur dioxide, methanol, and water are the dominant contributors. For the ortho / para ratios we find deviations from the high temperature values between 13 and 27 %. In total 14 % of all lines remain unidentified. △ Less

Submitted 27 April, 2021; originally announced April 2021.

Comments: 67 pages, 102 figures, submitted to A&A

Journal ref: A&A 651, A9 (2021)

Leveraging the ALMA Atacama Compact Array for Cometary Science: An Interferometric Survey of Comet C/2015 ER61 (PanSTARRS) and Evidence for a Distributed Source of Carbon Monosulfide

Authors: Nathan X. Roth, Stefanie N. Milam, Martin A. Cordiner, Dominique Bockelée-Morvan, Nicolas Biver, Jérémie Boissier, Dariusz C. Lis, Anthony J. Remijan, Steven B. Charnley

Abstract: We report the first survey of molecular emission from cometary volatiles using standalone Atacama Compact Array (ACA) observations of the Atacama Large Millimeter/Submillimeter Array (ALMA) toward comet C/2015 ER61 (PanSTARRS) carried out on UT 2017 April 11 and 15, shortly after its April 4 outburst. These measurements of HCN, CS, CH$_3$OH, H$_2$CO, and HNC (along with continuum emission from dus… ▽ More We report the first survey of molecular emission from cometary volatiles using standalone Atacama Compact Array (ACA) observations of the Atacama Large Millimeter/Submillimeter Array (ALMA) toward comet C/2015 ER61 (PanSTARRS) carried out on UT 2017 April 11 and 15, shortly after its April 4 outburst. These measurements of HCN, CS, CH$_3$OH, H$_2$CO, and HNC (along with continuum emission from dust) probed the inner coma of C/2015 ER61, revealing asymmetric outgassing and discerning parent from daughter/distributed source species. This work presents spectrally integrated flux maps, autocorrelation spectra, production rates, and parent scale lengths for each molecule, and a stringent upper limit for CO. HCN is consistent with direct nucleus release in C/2015 ER61, whereas CS, H$_2$CO, HNC, and potentially CH$_3$OH are associated with distributed sources in the coma. Adopting a Haser model, parent scale lengths determined for H$_2$CO (L$_p$ $\sim$ 2200 km) and HNC (L$_p$ $\sim$ 3300 km) are consistent with previous work in comets, whereas significant extended source production (L$_p$ $\sim$ 2000 km) is indicated for CS, suggesting production from an unknown parent in the coma. The continuum presents a point-source distribution, with a flux density implying an excessively large nucleus, inconsistent with other estimates of the nucleus size. It is best explained by the thermal emission of slowly-moving outburst ejectas, with total mass 5--8 $\times$ 10$^{10}$ kg. These results demonstrate the power of the ACA for revealing the abundances, spatial distributions, and locations of molecular production for volatiles in moderately bright comets such as C/2015 ER61. △ Less

Submitted 25 October, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

The Core Mass Function in the Orion Nebula Cluster Region: What Determines the Final Stellar Masses?

Authors: Hideaki Takemura, Fumitaka Nakamura, Shuo Kong, Héctor G. Arce, John M. Carpenter, Volker Ossenkopf-Okada, Ralf Klessen, Patricio Sanhueza, Yoshito Shimajiri, Takashi Tsukagoshi, Ryohei Kawabe, Shun Ishii, Kazuhito Dobashi, Tomomi Shimoikura, Paul F. Goldsmith, Álvaro Sánchez-Monge, Jens Kauffmann, Thushara Pillai, Paolo Padoan, Adam Ginsberg, Rowan J. Smith, John Bally, Steve Mairs, Jaime E. Pineda, Dariusz C. Lis , et al. (7 additional authors not shown)

Abstract: Applying dendrogram analysis to the CARMA-NRO C$^{18}$O ($J$=1--0) data having an angular resolution of $\sim$ 8", we identified 692 dense cores in the Orion Nebula Cluster (ONC) region. Using this core sample, we compare the core and initial stellar mass functions in the same area to quantify the step from cores to stars. About 22 \% of the identified cores are gravitationally bound. The derived… ▽ More Applying dendrogram analysis to the CARMA-NRO C$^{18}$O ($J$=1--0) data having an angular resolution of $\sim$ 8", we identified 692 dense cores in the Orion Nebula Cluster (ONC) region. Using this core sample, we compare the core and initial stellar mass functions in the same area to quantify the step from cores to stars. About 22 \% of the identified cores are gravitationally bound. The derived core mass function (CMF) for starless cores has a slope similar to Salpeter's stellar initial mass function (IMF) for the mass range above 1 $M_\odot$, consistent with previous studies. Our CMF has a peak at a subsolar mass of $\sim$ 0.1 $M_\odot$, which is comparable to the peak mass of the IMF derived in the same area. We also find that the current star formation rate is consistent with the picture in which stars are born only from self-gravitating starless cores. However, the cores must gain additional gas from the surroundings to reproduce the current IMF (e.g., its slope and peak mass), because the core mass cannot be accreted onto the star with a 100\% efficiency. Thus, the mass accretion from the surroundings may play a crucial role in determining the final stellar masses of stars. △ Less

Submitted 25 February, 2021; originally announced March 2021.

Comments: 13 pages, 5 figures, accepted by ApJL

High-resolution CARMA Observation of Molecular Gas in the North America and Pelican Nebulae

Authors: Shuo Kong, Héctor G. Arce, John M. Carpenter, John Bally, Volker Ossenkopf-Okada, Álvaro Sánchez-Monge, Anneila I. Sargent, Sümeyye Suri, Peregrine McGehee, Dariusz C. Lis, Ralf Klessen, Steve Mairs, Catherine Zucker, Rowan J. Smith, Fumitaka Nakamura, Thushara G. S. Pillai, Jens Kauffmann, Shaobo Zhang

Abstract: We present the first results from a CARMA high-resolution $^{12}$CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the North America and Pelican (NAP) Nebulae. CARMA observations have been combined with single-dish data from the Purple Mountain 13.7m telescope to add short spacings and produce high-dynamic-range images. We find that the molecular gas is predominantly shaped by t… ▽ More We present the first results from a CARMA high-resolution $^{12}$CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the North America and Pelican (NAP) Nebulae. CARMA observations have been combined with single-dish data from the Purple Mountain 13.7m telescope to add short spacings and produce high-dynamic-range images. We find that the molecular gas is predominantly shaped by the W80 HII bubble that is driven by an O star. Several bright rims are probably remnant molecular clouds heated and stripped by the massive star. Matching these rims in molecular lines and optical images, we construct a model of the three-dimensional structure of the NAP complex. Two groups of molecular clumps/filaments are on the near side of the bubble, one being pushed toward us, whereas the other is moving toward the bubble. Another group is on the far side of the bubble and moving away. The young stellar objects in the Gulf region reside in three different clusters, each hosted by a cloud from one of the three molecular clump groups. Although all gas content in the NAP is impacted by feedback from the central O star, some regions show no signs of star formation, while other areas clearly exhibit star formation activity. Other molecular gas being carved by feedback includes the cometary structures in the Pelican Head region and the boomerang features at the boundary of the Gulf region. The results show that the NAP complex is an ideal place for the study of feedback effects on star formation. △ Less

Submitted 7 March, 2021; originally announced March 2021.

Comments: accepted by AJ

Molecular composition of comet 46P/Wirtanen from millimetre-wave spectroscopy

Authors: N. Biver, D. Bockelée-Morvan, J. Boissier, R. Moreno, J. Crovisier, D. C. Lis, P. Colom, M. Cordiner, S. Milam, N. X. Roth, B. P. Bonev, N. Dello Russo, R. Vervack, M. A. DiSanti

Abstract: We present the results of a molecular survey of comet 46P/Wirtanen undertaken with the IRAM 30-m and NOEMA radio telescopes in December 2018. Observations at IRAM 30-m during the 12-18 Dec. period comprise a 2 mm spectral survey covering 25 GHz and a 1 mm survey covering 62 GHz. The gas outflow velocity and kinetic temperature have been accurately constrained by the observations. We derive abundan… ▽ More We present the results of a molecular survey of comet 46P/Wirtanen undertaken with the IRAM 30-m and NOEMA radio telescopes in December 2018. Observations at IRAM 30-m during the 12-18 Dec. period comprise a 2 mm spectral survey covering 25 GHz and a 1 mm survey covering 62 GHz. The gas outflow velocity and kinetic temperature have been accurately constrained by the observations. We derive abundances of 11 molecules, some being identified remotely for the first time in a Jupiter-family comet, including complex organic molecules such as formamide, ethylene glycol, acetaldehyde, or ethanol. Sensitive upper limits on the abundances of 24 other molecules are obtained. The comet is found to be relatively rich in methanol (3.4 percent relative to water), but relatively depleted in CO, CS, HNC, HNCO, and HCOOH. △ Less

Submitted 25 February, 2021; originally announced February 2021.

Comments: 35 pages, 16 figures, to be published in Astronomy & Astrophysics

Journal ref: A&A 648, A49 (2021)

Carbon-Chain Chemistry vs. Complex-Organic-Molecule Chemistry in Envelopes around Three Low-Mass Young Stellar Objects in the Perseus Region

Authors: Kotomi Taniguchi, Liton Majumdar, Shigehisa Takakuwa, Masao Saito, Dariusz C. Lis, Paul F. Goldsmith, Eric Herbst

Abstract: We have analyzed ALMA Cycle 5 data in Band 4 toward three low-mass young stellar objects (YSOs), IRAS 03235+3004 (hereafter IRAS 03235), IRAS 03245+3002 (IRAS 03245), and IRAS 03271+3013 (IRAS 03271), in the Perseus region. The HC$_{3}$N ($J=16-15$; $E_{\rm {up}}/k = 59.4$ K) line has been detected in all of the target sources, while four CH$_{3}$OH lines ($E_{\rm {up}}/k = 15.4-36.3$ K) have been… ▽ More We have analyzed ALMA Cycle 5 data in Band 4 toward three low-mass young stellar objects (YSOs), IRAS 03235+3004 (hereafter IRAS 03235), IRAS 03245+3002 (IRAS 03245), and IRAS 03271+3013 (IRAS 03271), in the Perseus region. The HC$_{3}$N ($J=16-15$; $E_{\rm {up}}/k = 59.4$ K) line has been detected in all of the target sources, while four CH$_{3}$OH lines ($E_{\rm {up}}/k = 15.4-36.3$ K) have been detected only in IRAS 03245. Sizes of the HC$_{3}$N distributions ($\sim 2930-3230$ au) in IRAS 03235 and IRAS 03245 are similar to those of the carbon-chain species in the warm carbon chain chemistry (WCCC) source L1527. The size of the CH$_{3}$OH emission in IRAS 03245 is $\sim 1760$ au, which is slightly smaller than that of HC$_{3}$N in this source. We compare the CH$_{3}$OH/HC$_{3}$N abundance ratios observed in these sources with predictions of chemical models. We confirm that the observed ratio in IRAS 03245 agrees with the modeled values at temperatures around 30--35 K, which supports the HC$_{3}$N formation by the WCCC mechanism. In this temperature range, CH$_{3}$OH does not thermally desorb from dust grains. Non-thermal desorption mechanisms or gas-phase formation of CH$_{3}$OH seem to work efficiently around IRAS 03245. The fact that IRAS 03245 has the highest bolometric luminosity among the target sources seems to support these mechanisms, in particular the non-thermal desorption mechanisms. △ Less

Submitted 19 February, 2021; originally announced February 2021.

Comments: Accepted by The Astrophysical Journal, 13 pages, 5 figures, 7 tables

4GREAT -- a four-color receiver for high-resolution airborne terahertz spectroscopy

Authors: Carlos A. Durán, Rolf Güsten, Christophe Risacher, Andrej Görlitz, Bernd Klein, Nicolas Reyes, Oliver Ricken, Hans-Joachim Wunsch, Urs U. Graf, Karl Jacobs, Cornelia E. Honingh, Jürgen Stutzki, Gert de Lange, Yan Delorme, Jean-Michel Krieg, Dariusz C. Lis

Abstract: 4GREAT is an extension of the German Receiver for Astronomy at Terahertz frequencies (GREAT) operated aboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). The spectrometer comprises four different detector bands and their associated subsystems for simultaneous and fully independent science operation. All detector beams are co-aligned on the sky. The frequency bands of 4GREAT cover… ▽ More 4GREAT is an extension of the German Receiver for Astronomy at Terahertz frequencies (GREAT) operated aboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). The spectrometer comprises four different detector bands and their associated subsystems for simultaneous and fully independent science operation. All detector beams are co-aligned on the sky. The frequency bands of 4GREAT cover 491-635, 890-1090, 1240-1525 and 2490-2590 GHz, respectively. This paper presents the design and characterization of the instrument, and its in-flight performance. 4GREAT saw first light in June 2018, and has been offered to the interested SOFIA communities starting with observing cycle 6. △ Less

Submitted 9 December, 2020; originally announced December 2020.

Comments: Accepted for publication in IEEE Transactions on Terahertz Science and Technology

The CARMA-NRO Orion Survey: Filament Formation via Collision-Induced Magnetic Reconnection -- The Stick in Orion A

Authors: Shuo Kong, Volker Ossenkopf-Okada, Héctor G. Arce, John Bally, Álvaro Sánchez-Monge, Peregrine McGehee, Sümeyye Suri, Ralf S. Klessen, John M. Carpenter, Dariusz C. Lis, Fumitaka Nakamura, Peter Schilke, Rowan J. Smith, Steve Mairs, Alyssa Goodman, María José Maureira

Abstract: A unique filament is identified in the {\it Herschel} maps of the Orion A giant molecular cloud. The filament, which, we name the Stick, is ruler-straight and at an early evolutionary stage. Transverse position-velocity diagrams show two velocity components closing in on the Stick. The filament shows consecutive rings/forks in C$^{18}$O(1-0) channel maps, which is reminiscent of structures generat… ▽ More A unique filament is identified in the {\it Herschel} maps of the Orion A giant molecular cloud. The filament, which, we name the Stick, is ruler-straight and at an early evolutionary stage. Transverse position-velocity diagrams show two velocity components closing in on the Stick. The filament shows consecutive rings/forks in C$^{18}$O(1-0) channel maps, which is reminiscent of structures generated by magnetic reconnection. We propose that the Stick formed via collision-induced magnetic reconnection (CMR). We use the magnetohydrodynamics (MHD) code Athena++ to simulate the collision between two diffuse molecular clumps, each carrying an anti-parallel magnetic field. The clump collision produces a narrow, straight, dense filament with a factor of $>$200 increase in density. The production of the dense gas is seven times faster than free-fall collapse. The dense filament shows ring/fork-like structures in radiative transfer maps. Cores in the filament are confined by surface magnetic pressure. CMR can be an important dense-gas-producing mechanism in the Galaxy and beyond. △ Less

Submitted 31 October, 2020; originally announced November 2020.

Comments: 39 pages, 31 figures, 2 tables, accepted by ApJ

Reconstructing EUV spectrum of star forming regions from millimeter recombination lines of HI, HeI, and HeII

Authors: Lena Murchikova, Eric J. Murphy, Dariusz C. Lis, Lee Armus, Selma de Mink, Kartik Sheth, Nadia Zakamska, Frank Tramper, Angela Bongiorno, Martin Elvis, Lisa Kewley, Hugues Sana

Abstract: The extreme ultraviolet (EUV) spectra of distant star-forming regions cannot be probed directly using either ground- or space-based telescopes due to the high cross-section for interaction of EUV photons with the interstellar medium. This makes EUV spectra poorly constrained. The mm/submm recombination lines of H and He, which can be observed from the ground, can serve as a reliable probe of the E… ▽ More The extreme ultraviolet (EUV) spectra of distant star-forming regions cannot be probed directly using either ground- or space-based telescopes due to the high cross-section for interaction of EUV photons with the interstellar medium. This makes EUV spectra poorly constrained. The mm/submm recombination lines of H and He, which can be observed from the ground, can serve as a reliable probe of the EUV. Here we present a study based on ALMA observations of three Galactic ultra-compact HII regions and the starburst region Sgr B2(M), in which we reconstruct the key parameters of the EUV spectra using mm recombination lines of HI, HeI and HeII. We find that in all cases the EUV spectra between 13.6 and 54.4 eV have similar frequency dependence: L_ν~ ν^{-4.5 +/- 0.4}. We compare the inferred values of the EUV spectral slopes with the values expected for a purely single stellar evolution model (Starburst99) and the Binary Population and Spectral Synthesis code (BPASS). We find that the observed spectral slope differs from the model predictions. This may imply that the fraction of interacting binaries in HII regions is substantially lower than assumed in BPASS. The technique demonstrated here allows one to deduce the EUV spectra of star forming regions providing critical insight into photon production rates at λ< 912 A and can serve as calibration to starburst synthesis models, improving our understanding of star formation in distant universe and the properties of ionizing flux during reionization. △ Less

Submitted 8 September, 2020; v1 submitted 26 June, 2020; originally announced June 2020.

Comments: v1: 18 pages, 5 figures, 1 table, 1 appendix v2: BPASS version 2.1 -> 2.2.1, SB99 tracks explicitly added to figures, evolution with metallicity added, few points clarified. 20 pages, 6 figures, 1 table. Accepted to ApJ

Journal ref: The Astrophysical Journal, Volume 903, Issue 1, id.29, 12 pp. (2020)

Unusually High CO Abundance of the First Active Interstellar Comet

Authors: M. A. Cordiner, S. N. Milam, N. Biver, D. Bockelée-Morvan, N. X. Roth, E. A. Bergin, E. Jehin, A. J. Remijan, S. B. Charnley, M. J. Mumma, J. Boissier, J. Crovisier, L. Paganini, Y. -J. Kuan, D. C Lis

Abstract: Comets spend most of their lives at large distances from any star, during which time their interior compositions remain relatively unaltered. Cometary observations can therefore provide direct insight into the chemistry that occurred during their birth at the time of planet formation. To-date, there have been no confirmed observations of parent volatiles (gases released directly from the nucleus)… ▽ More Comets spend most of their lives at large distances from any star, during which time their interior compositions remain relatively unaltered. Cometary observations can therefore provide direct insight into the chemistry that occurred during their birth at the time of planet formation. To-date, there have been no confirmed observations of parent volatiles (gases released directly from the nucleus) of a comet from any planetary system other than our own. Here we present high-resolution, interferometric observations of 2I/Borisov, the first confirmed interstellar comet, obtained using the Atacama Large Millimeter/submillimeter Array (ALMA) on 15th-16th December 2019. Our observations reveal emission from hydrogen cyanide (HCN), and carbon monoxide (CO), coincident with the expected position of 2I/Borisov's nucleus, with production rates Q(HCN)=$(7.0\pm1.1)\times10^{23}$ s$^{-1}$ and Q(CO)=$(4.4\pm0.7)\times10^{26}$ s$^{-1}$. While the HCN abundance relative to water (0.06-0.16%) appears similar to that of typical, previously observed comets in our Solar System, the abundance of CO (35-105%) is among the highest observed in any comet within 2 au of the Sun. This shows that 2I/Borisov must have formed in a relatively CO-rich environment - probably beyond the CO ice-line in the very cold, outer regions of a distant protoplanetary accretion disk, as part of a population of small, icy bodies analogous to our Solar System's own proto-Kuiper Belt. △ Less

Submitted 27 April, 2020; v1 submitted 20 April, 2020; originally announced April 2020.

New limits on the resonant absorption of solar axions obtained with a $^{169}$Tm-containing cryogenic detector

Authors: A. H. Abdelhameed, S. V. Bakhlanov, P. Bauer, A. Bento, E. Bertoldo, L. Canonica, A. V. Derbin, I. S. Drachnev, N. Ferreiro Iachellini, D. Fuchs, D. Hauff, M. Laubenstein, D. A. Lis, I. S. Lomskaya, M. Mancuso, V. N. Muratova, S. Nagorny, S. Nisi, F. Petricca, F. Proebst, J. Rothe, V. V. Ryabchenkov, S. E. Sarkisov, D. A. Semenov, K. A. Subbotin , et al. (3 additional authors not shown)

Abstract: A search for resonant absorption of solar axions by $^{169}$Tm nuclei was carried out. A newly developed approach involving low-background cryogenic bolometer based on Tm$_3$Al$_5$O$_{12}$ crystal was used that allowed for significant improvement of sensitivity in comparison with previous $^{169}$Tm based experiments. The measurements performed with $8.18$ g crystal during $6.6$ days exposure yiel… ▽ More A search for resonant absorption of solar axions by $^{169}$Tm nuclei was carried out. A newly developed approach involving low-background cryogenic bolometer based on Tm$_3$Al$_5$O$_{12}$ crystal was used that allowed for significant improvement of sensitivity in comparison with previous $^{169}$Tm based experiments. The measurements performed with $8.18$ g crystal during $6.6$ days exposure yielded the following limits on axion couplings: $|g_{Aγ} (g_{AN}^0 + g_{AN}^3) \leq 1.44 \times 10^{-14}$ GeV$^{-1}$ and $|g_{Ae} (g_{AN}^0 + g_{AN}^3) \leq 2.81 \times 10^{-16}$. △ Less

Submitted 17 April, 2020; originally announced April 2020.

Comments: 7 pages, 5 figures

Star Cluster Formation in Orion A

Authors: Wanggi Lim, Fumitaka Nakamura, Benjamin Wu, Thomas G. Bisbas, Jonathan C. Tan, Edward Chambers, John Bally, Shuo Kong, Peregrine McGehee, Dariusz C. Lis, Volker Ossenkopf-Okada, Álvaro Sánchez-Monge

Abstract: We introduce new analysis methods for studying the star cluster formation processes in Orion A, especially examining the scenario of a cloud-cloud collision. We utilize the CARMA-NRO Orion survey $^{13}$CO (1-0) data to compare molecular gas to the properties of YSOs from the SDSS III IN-SYNC survey. We show that the increase of $v_{\rm 13CO} - v_{\rm YSO}$ and $Σ$ scatter of older YSOs can be sig… ▽ More We introduce new analysis methods for studying the star cluster formation processes in Orion A, especially examining the scenario of a cloud-cloud collision. We utilize the CARMA-NRO Orion survey $^{13}$CO (1-0) data to compare molecular gas to the properties of YSOs from the SDSS III IN-SYNC survey. We show that the increase of $v_{\rm 13CO} - v_{\rm YSO}$ and $Σ$ scatter of older YSOs can be signals of cloud-cloud collision. SOFIA-upGREAT 158$μ$m [CII] archival data toward the northern part of Orion A are also compared to the $^{13}$CO data to test whether the position and velocity offsets between the emission from these two transitions resemble those predicted by a cloud-cloud collision model. We find that the northern part of Orion A, including regions ONC-OMC-1, OMC-2, OMC-3 and OMC-4, shows qualitative agreements with the cloud-cloud collision scenario, while in one of the southern regions, NGC1999, there is no indication of such a process in causing the birth of new stars. On the other hand, another southern cluster, L1641N, shows slight tendencies of cloud-cloud collision. Overall, our results support the cloud-cloud collision process as being an important mechanism for star cluster formation in Orion A. △ Less

Submitted 25 June, 2020; v1 submitted 7 April, 2020; originally announced April 2020.

Comments: 16 pages, 11 figures, 4 tables, published in PASJ

Distribution of Water Vapor in Molecular Clouds. II

Authors: Gary J. Melnick, Volker Tolls, Ronald L. Snell, Michael J. Kaufman, Edwin A. Bergin, Javier R. Goicoechea, Paul F. Goldsmith, Eduardo González-Alfonso, David J. Hollenbach, Dariusz C. Lis, David A. Neufeld

Abstract: The depth-dependent abundance of both gas-phase and solid-state water within dense, quiescent, molecular clouds is important to both the cloud chemistry and gas cooling. Where water is in the gas phase, it's free to participate in the network of ion-neutral reactions that lead to a host of oxygen-bearing molecules, and its many ortho and para energy levels make it an effective coolant for gas temp… ▽ More The depth-dependent abundance of both gas-phase and solid-state water within dense, quiescent, molecular clouds is important to both the cloud chemistry and gas cooling. Where water is in the gas phase, it's free to participate in the network of ion-neutral reactions that lead to a host of oxygen-bearing molecules, and its many ortho and para energy levels make it an effective coolant for gas temperatures greater than 20K. Where water is abundant as ice on grain surfaces, and unavailable to cool the gas, significant amounts of oxygen are removed from the gas phase, suppressing the gas-phase chemical reactions that lead to a number of oxygen-bearing species, including O2. Models of FUV-illuminated clouds predict that the gas-phase water abundance peaks in the range Av ~3 and 8mag of the cloud surface, depending on the gas density and FUV field strength. Deeper within such clouds, water is predicted to exist mainly as ice on grain surfaces. More broadly, these models are used to analyze a variety of other regions, including outflow cavities associated with young stellar objects and the surface layers of protoplanetary disks. In this paper, we report the results of observational tests of FUV-illuminated cloud models toward the Orion Molecular Ridge and Cepheus B using data obtained from the Herschel Space Observatory and the Five College Radio Astronomy Observatory. Toward Orion, 2220 spatial positions were observed along the face-on Orion Ridge in the H2O 110-101 557GHz and NH3 J,K=1,0-0,0 572GHz lines. Toward Cepheus B, two strip scans were made in the same lines across the edge-on ionization front. These new observations demonstrate that gas-phase water exists primarily within a few magnitudes of dense cloud surfaces, strengthening the conclusions of an earlier study based on a much smaller data set, and indirectly supports the prediction that water ice is quite abundant in dense clouds. △ Less

Submitted 10 March, 2020; originally announced March 2020.

Comments: 59 pages, including 23 figures, Accepted for publication in ApJ

The water line emission and ortho-to-para ratio in the Orion Bar photon-dominated region

Authors: T. Putaud, X. Michaut, F. Le Petit, E. Roueff, D. C. Lis

Abstract: A very low ortho-to-para ratio (OPR) of 0.1-0.5 was previously reported in the Orion Bar photon-dominated region (PDR), based on observations of two optically thin $\mathrm{H_2^{18}O}$ lines which were analyzed by using a single-slab large velocity gradient model. The corresponding spin temperature does not coincide with the kinetic temperature of the molecular gas in this UV-illuminated region. T… ▽ More A very low ortho-to-para ratio (OPR) of 0.1-0.5 was previously reported in the Orion Bar photon-dominated region (PDR), based on observations of two optically thin $\mathrm{H_2^{18}O}$ lines which were analyzed by using a single-slab large velocity gradient model. The corresponding spin temperature does not coincide with the kinetic temperature of the molecular gas in this UV-illuminated region. This was interpreted as an indication of water molecules being formed on cold icy grains which were subsequently released by UV photodesorption. A more complete set of water observations in the Orion Bar, including seven $\mathrm{H_2^{16}O}$ lines and one $\mathrm{H_2^{18}O}$ line, carried out using Herschel/HIFI instrument, was reanalyzed using the Meudon PDR code to derive gas-phase water abundance and the OPR, taking into account the steep density and temperature gradients present in the region. The model line intensities are in good agreement with the observations assuming that water molecules formed with an OPR corresponding to thermal equilibrium conditions at the local kinetic temperature of the gas and when solely considering gas-phase chemistry and water gas-grain exchanges through adsorption and desorption. Gas-phase water is predicted to arise from a region deep into the cloud, corresponding to a visual extinction of $A_{\mathrm{V}} \sim 9$, with a $\mathrm{H_2^{16}O}$ fractional abundance of $\sim 2\times 10^{-7}$ and column density of $(1.4 \pm 0.8) \times 10^{15}$ cm$^{-2}$ for a total cloud depth of $A_{\mathrm{V}}=15$. A line-of-sight average ortho-to-para ratio of $2.8 \pm 0.2$ is derived. The observational data are consistent with a nuclear spin isomer repartition corresponding to the thermal equilibrium at a temperature of $36 \pm 2$ K, much higher than the spin temperature previously reported for this region and close to the gas kinetic temperature in the water-emitting gas. △ Less

Submitted 7 October, 2019; v1 submitted 1 August, 2019; originally announced August 2019.

Journal ref: A&A 632, A8 (2019)

The Physical and chemical structure of Sagittarius B2 -- IV. Converging filaments in the high-mass cluster forming region Sgr B2(N)

Authors: A. Schwörer, Á. Sánchez-Monge, P. Schilke, T. Möller, A. Ginsburg, F. Meng, A. Schmiedeke, H. S. P. Müller, D. Lis, S. -L. Qin

Abstract: We have used an unbiased, spectral line-survey that covers the frequency range from 211 to 275 GHz and was obtained with ALMA (angular resolution of 0.4 arcsec) to study the small-scale structure of the dense gas in Sagittarius B2 (north). Eight filaments are found converging to the central hub and extending for about 0.1 pc. The spatial structure, together with the presence of the massive central… ▽ More We have used an unbiased, spectral line-survey that covers the frequency range from 211 to 275 GHz and was obtained with ALMA (angular resolution of 0.4 arcsec) to study the small-scale structure of the dense gas in Sagittarius B2 (north). Eight filaments are found converging to the central hub and extending for about 0.1 pc. The spatial structure, together with the presence of the massive central region, suggest that these filaments may be associated with accretion processes. In order to derive the kinematic properties of the gas in a chemically line-rich source like Sgr B2(N), we have developed a new tool that stacks all the detected transition lines of any molecular species. This permits to increase the signal-to-noise ratio of our observations and average out line blending effects, which are a common problem in line-rich regions. We derive velocity gradients along the filaments of about 20-100 km s$^{-1}$ pc$^{-1}$, which are 10-100 times larger than those typically found on larger scales (1 pc) in other star-forming regions. The mass accretion rates of individual filaments are about 0.05 M$_\odot$ yr$^{-1}$, which result in a total accretion rate of 0.16 M$_\odot$ yr$^{-1}$. Some filaments harbor dense cores that are likely forming stars and stellar clusters. The stellar content of these dense cores is on the order of 50% of the total mass. We conclude that the cores may merge in the center when already forming stellar clusters but still containing a significant amount of gas, resulting in a "damp" merger. The high density and mass of the central region, combined with the presence of converging filaments with high mass, high accretion rates and embedded dense cores already forming stars, suggest that Sgr B2(N) may have the potential to evolve into a super stellar cluster. △ Less

Submitted 26 June, 2019; originally announced June 2019.

Comments: Accepted for publication in Astronomy & Astrophysics. 10 pages with 8 figures and 2 tables, plus Appendices with extra figures and tables. (Abstract modified)

Journal ref: A&A 628, A6 (2019)

A Test of Bolometric Properties of Tm-containing Crystals as a Perspective Detector for the Solar Axion Search

Authors: E. Bertoldo, A. V. Derbin, I. S. Drachnev, M. Laubenstein, D. A. Lis, M. Mancuso, V. N. Muratova, S. Nagorny, S. Nisi, F. Petricca, V. V. Ryabchenkov, S. E. Sarkisov, D. A. Semenov, K. A. Subbotin, E. V. Unzhakov, E. V. Zharikov

Abstract: The $^{169}$Tm nuclide has first nuclear level at 8.41 keV with magnetic type transition to the ground state and, therefore, can be used as a target nucleus for the search of resonant absorption of solar axions. We plan to use a Tm-containing crystal of a garnet family Tm$_3$Al$_5$O$_{12}$ as a bolometric detector in order to search for the excitation of the first nuclear level of $^{169}$Tm via t… ▽ More The $^{169}$Tm nuclide has first nuclear level at 8.41 keV with magnetic type transition to the ground state and, therefore, can be used as a target nucleus for the search of resonant absorption of solar axions. We plan to use a Tm-containing crystal of a garnet family Tm$_3$Al$_5$O$_{12}$ as a bolometric detector in order to search for the excitation of the first nuclear level of $^{169}$Tm via the resonant absorption of solar axions. With this perspective in mind, a sample of the Tm$_3$Al$_5$O$_{12}$ crystal was grown and tested for its bolometric and optical properties. Measurements of chemical and/or radioactive contaminations were performed as well. In this paper we present the test results and estimate the requirements for a future low-background experimental setup. △ Less

Submitted 30 May, 2019; originally announced May 2019.

Comments: 9 pages, 6 figures, 3 tables

Terrestrial deuterium-to-hydrogen ratio in water in hyperactive comets

Authors: D. C. Lis, D. Bockelée-Morvan, R. Güsten, N. Biver, J. Stutzki, Yan Delorme, C. Durán, H. Wiesemeyer, Y. Okada

Abstract: The D/H ratio in cometary water has been shown to vary between 1 and 3 times the Earth's oceans value, in both Oort cloud comets and Jupiter-family comets originating from the Kuiper belt. We present new sensitive spectroscopic observations of water isotopologues in the Jupiter-family comet 46P/Wirtanen carried out using the GREAT spectrometer aboard the Stratospheric Observatory for Infrared Astr… ▽ More The D/H ratio in cometary water has been shown to vary between 1 and 3 times the Earth's oceans value, in both Oort cloud comets and Jupiter-family comets originating from the Kuiper belt. We present new sensitive spectroscopic observations of water isotopologues in the Jupiter-family comet 46P/Wirtanen carried out using the GREAT spectrometer aboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). The derived D/H ratio of $(1.61 \pm 0.65) \times 10^{-4}$ is the same as in the Earth's oceans. Although the statistics are limited, we show that interesting trends are already becoming apparent in the existing data. A clear anti-correlation is seen between the D/H ratio and the active fraction, defined as the ratio of the active surface area to the total nucleus surface. Comets with an active fraction above 0.5 typically have D/H ratios in water consistent with the terrestrial value. These hyperactive comets, such as 46P/Wirtanen, require an additional source of water vapor in their coma, explained by the presence of subliming icy grains expelled from the nucleus. The observed correlation may suggest that hyperactive comets belong to a population of ice-rich objects that formed just outside the snow line, or in the outermost regions of the solar nebula, from water thermally reprocessed in the inner disk that was transported outward during the early disk evolution. The observed anti-correlation between the active fraction and the nucleus size seems to argue against the first interpretation, as planetesimals near the snow line are expected to undergo rapid growth. Alternatively, isotopic properties of water outgassed from the nucleus and icy grains may be different due to fractionation effects at sublimation. In this case, all comets may share the same Earth-like D/H ratio in water, with profound implications for the early solar system and the origin of Earth's oceans. △ Less

Submitted 23 April, 2019; v1 submitted 19 April, 2019; originally announced April 2019.

Comments: 9 pages, 5 figures

Journal ref: A&A 625, L5 (2019)

The trail of water and the delivery of volatiles to habitable planets

Authors: Klaus M. Pontoppidan, Andrea Banzatti, Edwin Bergin, Geoffrey A. Blake, Sean Brittain, Maryvonne Gerin, Paul Goldsmith, Quentin Kral, David Leisawitz, Dariusz Lis, Melissa McClure, Stefanie Milam, Gary Melnick, Joan Najita, Karin Öberg, Matt Richter, Colette Salyk, Martina Wiedner, Ke Zhang

Abstract: Water is fundamental to our understanding of the evolution of planetary systems and the delivery of volatiles to the surfaces of potentially habitable planets. Yet, we currently have essentially no facilities capable of observing this key species comprehensively. With this white paper, we argue that we need a relatively large, cold space-based observatory equipped with a high-resolution spectromet… ▽ More Water is fundamental to our understanding of the evolution of planetary systems and the delivery of volatiles to the surfaces of potentially habitable planets. Yet, we currently have essentially no facilities capable of observing this key species comprehensively. With this white paper, we argue that we need a relatively large, cold space-based observatory equipped with a high-resolution spectrometer, in the mid- through far-infrared wavelength range (25-600~$μ$m) in order to answer basic questions about planet formation, such as where the Earth got its water, how giant planets and planetesimals grow, and whether water is generally available to planets forming in the habitable zone of their host stars. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

ALMA Autocorrelation Spectroscopy of Comets: The HCN/H^13CN ratio in C/2012 S1 (ISON)

Authors: M. A. Cordiner, M. Y. Palmer, M. de Val-Borro, S. B. Charnley, L. Paganini, G. Villanueva, D. Bockelée-Morvan, N. Biver, A. J. Remijan, Y. -J. Kuan, S. N. Milam, J. Crovisier, D. C. Lis, M. J. Mumma

Abstract: The Atacama Large Millimeter/submillimeter Array (ALMA) is a powerful tool for high-resolution mapping of comets, but the main interferometer (comprised of 50x12-m antennas) is insensitive to the largest coma scales due to a lack of very short baselines. In this work, we present a new technique employing ALMA autocorrelation data (obtained simultaneously with the interferometric observations), eff… ▽ More The Atacama Large Millimeter/submillimeter Array (ALMA) is a powerful tool for high-resolution mapping of comets, but the main interferometer (comprised of 50x12-m antennas) is insensitive to the largest coma scales due to a lack of very short baselines. In this work, we present a new technique employing ALMA autocorrelation data (obtained simultaneously with the interferometric observations), effectively treating the entire 12-m array as a collection of single-dish telescopes. Using combined autocorrelation spectra from 28 active antennas, we recovered extended HCN coma emission from comet C/2012 S1 (ISON), resulting in a fourteen-fold increase in detected line brightness compared with the interferometer. This resulted in the first detection of rotational emission from H^13CN in this comet. Using a detailed coma radiative transfer model accounting for optical depth and non-LTE excitation effects, we obtained an H^12CN/H^13CN ratio of 88+-18, which matches the terrestrial value of 89, consistent with a lack of isotopic fractionation in HCN during comet formation in the protosolar accretion disk. The possibility of future discoveries in extended sources using autocorrelation spectroscopy from the main ALMA array is thus demonstrated. △ Less

Submitted 24 January, 2019; originally announced January 2019.

The CARMA-NRO Orion Survey: The filamentary structure as seen in C$^{18}$O emission

Authors: S. T. Suri, A. Sanchez-Monge, P. Schilke, S. D. Clarke, R. J. Smith, V. Ossenkopf-Okada, R. Klessen, P. Padoan, P. Goldsmith, H. G. Arce, J. Bally, J. M. Carpenter, A. Ginsburg, D. Johnstone, J. Kauffmann, S. Kong, D. C. Lis, S. Mairs, T. Pillai, J. E. Pineda, A. Duarte-Cabral

Abstract: We present an initial overview of the filamentary structure in the Orion A molecular cloud utilizing a high angular and velocity resolution C$^{18}$O(1-0) emission map that was recently produced as part of the CARMA-NRO Orion Survey. The main goal of this study is to build a credible method to study varying widths of filaments which has previously been linked to star formation in molecular clouds.… ▽ More We present an initial overview of the filamentary structure in the Orion A molecular cloud utilizing a high angular and velocity resolution C$^{18}$O(1-0) emission map that was recently produced as part of the CARMA-NRO Orion Survey. The main goal of this study is to build a credible method to study varying widths of filaments which has previously been linked to star formation in molecular clouds. Due to the diverse star forming activities taking place throughout its $\sim$20 pc length, together with its proximity of 388 pc, the Orion A molecular cloud provides an excellent laboratory for such an experiment to be carried out with high resolution and high sensitivity. Using the widely-known structure identification algorithm, DisPerSE, on a 3-dimensional (PPV) C$^{18}$O cube, we identified 625 relatively short (the longest being 1.74 pc) filaments over the entire cloud. We study the distribution of filament widths using FilChaP, a python package that we have developed and made publicly available. We find that the filaments identified in a 2 square degree PPV cube do not overlap spatially, except for the complex OMC-4 region that shows distinct velocity components along the line of sight. The filament widths vary between 0.02 and 0.3 pc depending on the amount of substructure that a filament possesses. The more substructure a filament has, the larger is its width. We also find that despite this variation, the filament width shows no anticorrelation with the central column density which is in agreement with previous Herschel observations. △ Less

Submitted 1 January, 2019; originally announced January 2019.

Comments: 19 pages, 20 figures. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 623, A142 (2019)

ALMA observations of the young protostellar system Barnard 1b: signatures of an incipient hot corino in B1b-S

Authors: N. Marcelino, M. Gerin, J. Cernicharo, A. Fuente, H. A. Wootten, E. Chapillon, J. Pety, D. C. Lis, E. Roueff, B. Commerçon, A. Ciardi

Abstract: The Barnard 1b core shows signatures of being at the earliest stages of low-mass star formation, with two extremely young and deeply embedded protostellar objects. Hence, this core is an ideal target to study the structure and chemistry of the first objects formed in the collapse of prestellar cores. We present ALMA Band 6 spectral line observations at ~0.6'' of angular resolution towards Barnard… ▽ More The Barnard 1b core shows signatures of being at the earliest stages of low-mass star formation, with two extremely young and deeply embedded protostellar objects. Hence, this core is an ideal target to study the structure and chemistry of the first objects formed in the collapse of prestellar cores. We present ALMA Band 6 spectral line observations at ~0.6'' of angular resolution towards Barnard 1b. We have extracted the spectra towards both protostars, and used a Local Thermodynamic Equilibrium (LTE) model to reproduce the observed line profiles. B1b-S shows rich and complex spectra, with emission from high energy transitions of complex molecules, such as CH3OCOH and CH3CHO, including vibrational level transitions. We have tentatively detected for the first time in this source emission from NH2CN, NH2CHO, CH3CH2OH, CH2OHCHO, CH3CH2OCOH and both aGg' and gGg' conformers of (CH2OH)2. This is the first detection of ethyl formate (CH3CH2OCOH) towards a low-mass star forming region. On the other hand, the spectra of the FHSC candidate B1b-N are free of COMs emission. In order to fit the observed line profiles in B1b-S, we used a source model with two components: an inner hot and compact component (200 K, 0.35'') and an outer and colder one (60 K, 0.6''). The resulting COM abundances in B1b-S range from 1e-13 for NH2CN and NH2CHO, up to 1e-9 for CH3OCOH. Our ALMA Band 6 observations reveal the presence of a compact and hot component in B1b-S, with moderate abundances of complex organics. These results indicate that a hot corino is being formed in this very young Class 0 source. △ Less

Submitted 21 September, 2018; originally announced September 2018.

Comments: Accepted for publication in A&A

Journal ref: A&A 620, A80 (2018)

The CARMA-NRO Orion Survey

Authors: Shuo Kong, Héctor G. Arce, Jesse R. Feddersen, John M. Carpenter, Fumitaka Nakamura, Yoshito Shimajiri, Andrea Isella, Volker Ossenkopf-Okada, Anneila I. Sargent, Álvaro Sánchez-Monge, Sümeyye T. Suri, Jens Kauffmann, Thushara Pillai, Jaime E. Pineda, Jin Koda, John Bally, Dariusz C. Lis, Paolo Padoan, Ralf Klessen, Steve Mairs, Alyssa Goodman, Paul Goldsmith, Peregrine McGehee, Peter Schilke, Peter J. Teuben , et al. (13 additional authors not shown)

Abstract: We present the first results from a new, high resolution, $^{12}$CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the Orion A cloud, hereafter referred to as the CARMA-NRO Orion Survey. CARMA observations have been combined with single-dish data from the Nobeyama 45m telescope to provide extended images at about 0.01 pc resolution, with a dynamic range of approximately 1200 in… ▽ More We present the first results from a new, high resolution, $^{12}$CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the Orion A cloud, hereafter referred to as the CARMA-NRO Orion Survey. CARMA observations have been combined with single-dish data from the Nobeyama 45m telescope to provide extended images at about 0.01 pc resolution, with a dynamic range of approximately 1200 in spatial scale. Here we describe the practical details of the data combination in uv space, including flux scale matching, the conversion of single dish data to visibilities, and joint deconvolution of single dish and interferometric data. A $Δ$-variance analysis indicates that no artifacts are caused by combining data from the two instruments. Initial analysis of the data cubes, including moment maps, average spectra, channel maps, position-velocity diagrams, excitation temperature, column density, and line ratio maps provides evidence of complex and interesting structures such as filaments, bipolar outflows, shells, bubbles, and photo-eroded pillars. The implications for star formation processes are profound and follow-up scientific studies by the CARMA-NRO Orion team are now underway. We plan to make all the data products described here generally accessible; some are already available at https://dataverse.harvard.edu/dataverse/CARMA-NRO-Orion △ Less

Submitted 30 March, 2018; originally announced March 2018.

Comments: 46 pages, 28 figures, 2 tables, accepted by ApJS

French SKA White Book - The French Community towards the Square Kilometre Array

Authors: F. Acero, J. -T. Acquaviva, R. Adam, N. Aghanim, M. Allen, M. Alves, R. Ammanouil, R. Ansari, A. Araudo, E. Armengaud, B. Ascaso, E. Athanassoula, D. Aubert, S. Babak, A. Bacmann, A. Banday, K. Barriere, F. Bellossi, J. -P. Bernard, M. G. Bernardini, M. Béthermin, E. Blanc, L. Blanchet, J. Bobin, S. Boissier , et al. (153 additional authors not shown)

Abstract: The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive p… ▽ More The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration. △ Less

Submitted 28 March, 2018; v1 submitted 19 December, 2017; originally announced December 2017.

Comments: Editor in chief: C. Ferrari; Editors: M. Alves, S. Bosse, S. Corbel, A. Ferrari, K. Ferrière, S. Gauffre, E. Josselin, G. Lagache, S. Lambert, G. Marquette, J.-M. Martin, M.-A. Miville-Deschênes, L. Montier, B. Semelin, G. Theureau, S. Vergani, N. Vilmer, P. Zarka; Original file with high resolution figures at SKA-France link: https://ska-france.oca.eu/images/SKA-France-Media/FWB_051017.pdf

Evidence for disks at an early stage in class 0 protostars?

Authors: M. Gerin, J. Pety, B. Commercon, A. Fuente, J. Cernicharo, N. Marcelino, A. Ciardi, D. C. Lis, E. Roueff, H. A. Wootten, E. Chapillon

Abstract: The formation epoch of protostellar disks is debated because of the competing roles of rotation, turbulence, and magnetic fields in the early stages of low-mass star formation. Magnetohydrodynamics simulations of collapsing cores predict that rotationally supported disks may form in strongly magnetized cores through ambipolar diffusion or misalignment between the rotation axis and the magnetic fie… ▽ More The formation epoch of protostellar disks is debated because of the competing roles of rotation, turbulence, and magnetic fields in the early stages of low-mass star formation. Magnetohydrodynamics simulations of collapsing cores predict that rotationally supported disks may form in strongly magnetized cores through ambipolar diffusion or misalignment between the rotation axis and the magnetic field orientation. Detailed studies of individual sources are needed to cross check the theoretical predictions. We present 0.06-0.1" resolution images at 350 GHz toward B1b-N and B1b-S, which are young class 0 protostars, possibly first hydrostatic cores. The images have been obtained with ALMA, and we compare these data with magnetohydrodynamics simulations of a collapsing turbulent and magnetized core. The submillimeter continuum emission is spatially resolved by ALMA. Compact structures with optically thick 350 GHz emission are detected toward both B1b-N and B1b-S, with 0.2 and 0.35" radii (46 and 80 au at the Perseus distance of 230 pc), within a more extended envelope. The flux ratio between the compact structure and the envelope is lower in B1b-N than in B1b-S, in agreement with its earlier evolutionary status. The size and orientation of the compact structure are consistent with 0.2" resolution 32 GHz observations obtained with the Very Large Array as a part of the VANDAM survey, suggesting that grains have grown through coagulation. The morphology, temperature, and densities of the compact structures are consistent with those of disks formed in numerical simulations of collapsing cores. Moreover, the properties of B1b-N are consistent with those of a very young protostar, possibly a first hydrostatic core. These observations provide support for the early formation of disks around low-mass protostars. △ Less

Submitted 22 September, 2017; originally announced September 2017.

Comments: 8 pages, 10 PDF figures. Accepted for publication in A&A. Uses aa and natbib latex macros

Review: Far-Infrared Instrumentation and Technology Development for the Next Decade

Authors: Duncan Farrah, Kimberly Ennico Smith, David Ardila, Charles M. Bradford, Michael Dipirro, Carl Ferkinhoff, Jason Glenn, Paul Goldsmith, David Leisawitz, Thomas Nikola, Naseem Rangwala, Stephen A. Rinehart, Johannes Staguhn, Michael Zemcov, Jonas Zmuidzinas, James Bartlett, Sean Carey, William J. Fischer, Julia Kamenetzky, Jeyhan Kartaltepe, Mark Lacy, Dariusz C. Lis, Lisa Locke, Enrique Lopez-Rodriguez, Meredith MacGregor , et al. (11 additional authors not shown)

Abstract: Far-infrared astronomy has advanced rapidly since its inception in the late 1950's, driven by a maturing technology base and an expanding community of researchers. This advancement has shown that observations at far-infrared wavelengths are important in nearly all areas of astrophysics, from the search for habitable planets and the origin of life, to the earliest stages of galaxy assembly in the f… ▽ More Far-infrared astronomy has advanced rapidly since its inception in the late 1950's, driven by a maturing technology base and an expanding community of researchers. This advancement has shown that observations at far-infrared wavelengths are important in nearly all areas of astrophysics, from the search for habitable planets and the origin of life, to the earliest stages of galaxy assembly in the first few hundred million years of cosmic history. The combination of a still developing portfolio of technologies, particularly in the field of detectors, and a widening ensemble of platforms within which these technologies can be deployed, means that far-infrared astronomy holds the potential for paradigm-shifting advances over the next decade. In this review, we examine current and future far-infrared observing platforms, including ground-based, sub-orbital, and space-based facilities, and discuss the technology development pathways that will enable and enhance these platforms to best address the challenges facing far-infrared astronomy in the 21st century. △ Less

Submitted 8 January, 2019; v1 submitted 7 September, 2017; originally announced September 2017.

Comments: Invited review article, submitted to the Journal of Astronomical Telescopes, Instruments, and Systems. Updated in light of referee reports

Journal ref: J. Astron. Telesc. Instrum. Syst. 5(2), 020901 (2019)

Chemical segregation in the young protostars Barnard 1b-N and S: evidence of pseudo-disk rotation in Barnard 1b-S

Authors: Asunción Fuente, Maryvonne Gerin, Jerome Pety, Benoit Commercon, Marcelino Agúndez, José Cernicharo, Nuria Marcelino, Evelyne Roueff, Darek C. Lis, H. A. Wootten

Abstract: The extremely young Class 0 object B1b-S and the first hydrostatic core (FSHC) candidate, B1b-N, provide a unique opportunity to study the chemical changes produced in the elusive transition from the prestellar core to the protostellar phase. We present 40"x70" images of Barnard 1b in the 13CO 1->0, C18O 1->0, NH2D 1_{1,1}a->1_{0,1}s, and SO 3_2->2_1 lines obtained with the NOEMA interferometer. T… ▽ More The extremely young Class 0 object B1b-S and the first hydrostatic core (FSHC) candidate, B1b-N, provide a unique opportunity to study the chemical changes produced in the elusive transition from the prestellar core to the protostellar phase. We present 40"x70" images of Barnard 1b in the 13CO 1->0, C18O 1->0, NH2D 1_{1,1}a->1_{0,1}s, and SO 3_2->2_1 lines obtained with the NOEMA interferometer. The observed chemical segregation allows us to unveil the physical structure of this young protostellar system down to scales of ~500au. The two protostellar objects are embedded in an elongated condensation, with a velocity gradient of ~0.2-0.4 m s^{-1} au^{-1} in the east-west direction, reminiscent of an axial collapse. The NH2D data reveal cold and dense pseudo-disks (R~500-1000 au) around each protostar. Moreover, we observe evidence of pseudo-disk rotation around B1b-S. We do not see any signature of the bipolar outflows associated with B1b-N and B1b-S, which were previously detected in H2CO and CH3OH, in any of the imaged species. The non-detection of SO constrains the SO/CH3OH abundance ratio in the high-velocity gas. △ Less

Submitted 31 August, 2017; originally announced September 2017.

Comments: 9 pages, 4 figures

The onset of energetic particle irradiation in Class 0 protostars

Authors: C. Favre, A. Lopez-Sepulcre, C. Ceccarelli, C. Dominik, P. Caselli, E. Caux, A. Fuente, M. Kama, J. Le Bourlot, B. Lefloch, D. Lis, T. Montmerle, M. Padovani, C. Vastel

Abstract: The early stages of low-mass star formation are likely to be subject to intense ionization by protostellar energetic MeV particles. As a result, the surrounding gas is enriched in molecular ions, such as HCO$^{+}$ and N$_{2}$H$^{+}$. Nonetheless, this phenomenon remains poorly understood for Class 0 objects. Recently, based on Herschel observations taken as part of the key program Chemical HErsche… ▽ More The early stages of low-mass star formation are likely to be subject to intense ionization by protostellar energetic MeV particles. As a result, the surrounding gas is enriched in molecular ions, such as HCO$^{+}$ and N$_{2}$H$^{+}$. Nonetheless, this phenomenon remains poorly understood for Class 0 objects. Recently, based on Herschel observations taken as part of the key program Chemical HErschel Surveys of Star forming regions (CHESS), a very low HCO$^{+}$/N$_{2}$H$^{+}$ abundance ratio of about 3-4, has been reported toward the protocluster OMC-2 FIR4. This finding suggests a cosmic-ray ionization rate in excess of 10$^{-14}$ s$^{-1}$, much higher than the canonical value of $ζ$ = 3$\times$10$^{-17}$ s$^{-1}$ (value expected in quiescent dense clouds). To assess the specificity of OMC-2 FIR4, we have extended this study to a sample of sources in low- and intermediate mass. More specifically, we seek to measure the HCO$^{+}$/N$_2$H$^{+}$ abundance ratio from high energy lines (J $\ge$ 6) toward this source sample in order to infer the flux of energetic particles in the warm and dense gas surrounding the protostars. We use observations performed with the Heterodyne Instrument for the FarInfrared spectrometer on board the Herschel Space Observatory toward a sample of 9 protostars. We report HCO$^{+}$/N$_2$H$^{+}$ abundance ratios in the range of 5 up to 73 toward our source sample. The large error bars do not allow us to conclude whether OMC-2~FIR4 is a peculiar source. Nonetheless, an important result is that the measured HCO$^{+}$/N$_2$H$^{+}$ ratio does not vary with the source luminosity. At the present time, OMC-2 FIR4 remains the only source where a high flux of energetic particles is clearly evident. More sensitive and higher angular resolution observations are required to further investigate this process. △ Less

Submitted 28 August, 2017; originally announced August 2017.

Comments: Accepted for publication in A&A 18 pages, 12 figures, 15 tables

Journal ref: A&A 608, A82 (2017)

Survey of cold water lines in protoplanetary disks: indications of systematic volatile depletion

Authors: Fujun Du, Edwin Anthony Bergin, Michiel Hogerheijde, Ewine F. van Dishoeck, Geoff Blake, Simon Bruderer, Ilse Cleeves, Carsten Dominik, Davide Fedele, Dariusz C. Lis, Gary Melnick, David Neufeld, John Pearson, Umut Yildiz

Abstract: We performed very deep searches for 2 ground-state water transitions in 13 protoplanetary disks with the HIFI instrument on-board the Herschel Space Observatory, with integration times up to 12 hours per line. Two other water transitions that sample warmer gas were also searched for with shallower integrations. The detection rate is low, and the upper limits provided by the observations are genera… ▽ More We performed very deep searches for 2 ground-state water transitions in 13 protoplanetary disks with the HIFI instrument on-board the Herschel Space Observatory, with integration times up to 12 hours per line. Two other water transitions that sample warmer gas were also searched for with shallower integrations. The detection rate is low, and the upper limits provided by the observations are generally much lower than predictions of thermo-chemical models with canonical inputs. One ground-state transition is newly detected in the stacked spectrum of AA Tau, DM Tau, LkCa 15, and MWC 480. We run a grid of models to show that the abundance of gas-phase oxygen needs to be reduced by a factor of at least ~100 to be consistent with the observational upper limits (and positive detections) if a dust-to-gas mass ratio of 0.01 were to be assumed. As a continuation of previous ideas, we propose that the underlying reason for the depletion of oxygen (hence the low detection rate) is the freeze-out of volatiles such as water and CO onto dust grains followed by grain growth and settling/migration, which permanently removes these gas-phase molecules from the emissive upper layers of the outer disk. Such depletion of volatiles is likely ubiquitous among different disks, though not necessarily to the same degree. The volatiles might be returned back to the gas phase in the inner disk (within about 15 AU), which is consistent with current constraints. Comparison with studies on disk dispersal due to photoevaporation indicates that the timescale for volatile depletion is shorter than that of photoevaporation. △ Less

Submitted 2 May, 2017; originally announced May 2017.

Comments: Accepted for publication in ApJ

Nature of shocks revealed by SOFIA OI observations in the Cepheus E protostellar outflow

Authors: A. Gusdorf, S. Anderl, B. Lefloch, S. Leurini, H. Wiesemeyer, R. Guesten, M. Benedettini, C. Codella, B. Godard, A. I. Gomez-Ruiz, K. Jacobs, L. E. Kristensen, P. Lesaffre, G. Pineau des Forets, D. C. Lis

Abstract: Protostellar jets and outflows are key features of the star-formation process, and primary processes of the feedback of young stars on the interstellar medium. Understanding the underlying shocks is necessary to explain how jets and outflows are launched, and to quantify their chemical and energetic impacts on the surrounding medium. We performed a high-spectral resolution study of the [OI]… ▽ More Protostellar jets and outflows are key features of the star-formation process, and primary processes of the feedback of young stars on the interstellar medium. Understanding the underlying shocks is necessary to explain how jets and outflows are launched, and to quantify their chemical and energetic impacts on the surrounding medium. We performed a high-spectral resolution study of the [OI]$_{\rm 63 μm}$ emission in the outflow of the intermediate-mass Class 0 protostar Cep E-mm. We present observations of the OI $^3$P$_1 \rightarrow$ $^3$P$_2$, OH between $^2Π_{1/2}$ $J = 3/2$ and $J = 1/2$ at 1837.8 GHz, and CO (16-15) lines with SOFIA-GREAT at three positions in the Cep E outflow: mm (the driving protostar), BI (in the southern lobe), and BII (the terminal position in the southern lobe). The CO line is detected at all three positions. The OI line is detected in BI and BII, whereas the OH line is not detected. In BII, we identify three kinematical components in OI and CO, already detected in CO: the jet, the HH377 terminal bow-shock, and the outflow cavity. The OI column density is higher in the outflow cavity than in the jet, which itself is higher than in the terminal shock. The terminal shock is where the abundance ratio of OI to CO is the lowest (about 0.2), whereas the jet component is atomic (ratio $\sim$2.7). In the jet, we compare the OI observations with shock models that successfully fit the integrated intensity of 10 CO lines: these models do not fit the OI data. The high intensity of OI emission points towards the propagation of additional dissociative or alternative FUV-irradiated shocks, where the illumination comes from the shock itself. From the sample of low-to-high mass protostellar outflows where similar observations have been performed, the effects of illumination seem to increase with the mass of the protostar. △ Less

Submitted 12 April, 2017; originally announced April 2017.

Comments: 9 pages, 4 figures, 6 tables, Astronomy and astrophysics, in press

Journal ref: A&A 602, A8 (2017)

The physical and chemical structure of Sagittarius B2. II. Continuum millimeter emission of SgrB2(M) and SgrB2(N) with ALMA

Authors: A. Sanchez-Monge, P. Schilke, A. Schmiedeke, A. Ginsburg, R. Cesaroni, D. C. Lis, S. -L. Qin, H. S. P. Mueller, E. Bergin, C. Comito, Th. Moeller

Abstract: The high-mass star forming sites SgrB2(M) and SgrB2(N) have been the target of numerous studies, revealing e.g. a rich chemistry. We want to characterize their physical and chemical structure using ALMA high-angular resolution observations at mm wavelengths, reaching spatial scales of about 4000 au, and covering the whole band 6 (from 211 to 275 GHz). In order to determine the continuum emission i… ▽ More The high-mass star forming sites SgrB2(M) and SgrB2(N) have been the target of numerous studies, revealing e.g. a rich chemistry. We want to characterize their physical and chemical structure using ALMA high-angular resolution observations at mm wavelengths, reaching spatial scales of about 4000 au, and covering the whole band 6 (from 211 to 275 GHz). In order to determine the continuum emission in line-rich sources, we use a new statistical method: STATCONT. We detect 27 continuum sources in SgrB2(M) and 20 in SgrB2(N). We study the continuum emission across the ALMA band 6, and compare it with previous SMA 345 GHz and VLA 40 GHz observations, to study the nature of the sources detected. The brightest sources are dominated by (partially optically thick) dust emission, while there is an important degree of contamination from ionized gas free-free emission in weaker sources. While the total mass in SgrB2(M) is distributed in many fragments, most of the mass in SgrB2(N) arises from a single object, with filamentary-like structures converging towards the center. There seems to be a lack of low-mass dense cores in both regions. We determine H2 volume densities for the cores of about 10^5-10^7 Msun pc^-3, one to two orders of magnitude higher than the stellar densities of super star clusters. In general, SgrB2(N) is chemically richer than SgrB2(M). There seems to be a correlation between the chemical richness and the mass of the fragments, with more massive clumps being more chemically rich. Both SgrB2(N) and SgrB2(M) harbour a cluster of hot molecular cores. We compare the continuum images with predictions from a detailed 3D radiative transfer model that reproduces the structure of SgrB2 from 45 pc down to 100 au. This dataset, together with ongoing projects in the range 5 to 200 GHz, better constrain the 3D structure of SgrB2, and allow us to understand its physical and chemical structure. △ Less

Submitted 6 April, 2017; originally announced April 2017.

Comments: Accepted for publication in Astronomy & Astrophysics. 20 pages with 13 figures and 4 tables, plus Appendices with extra figures and tables. (Abstract modified)

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

ALMA mapping of rapid gas and dust variations in comet C/2012 S1 (ISON): new insights into the origin of cometary HNC

Authors: M. A. Cordiner, J. Boissier, S. B. Charnley, A. J. Remijan, M. J. Mumma, G. Villanueva, D. C. Lis, S. N. Milam, L. Paganini, J. Crovisier, D. Bockelee-Morvan, Y. -J. Kuan, N. Biver, I. M. Coulson

Abstract: Observations of the sungrazing comet C/2012 S1 (ISON) were carried out using the Atacama Large Millimeter/submillimeter Array (ALMA) at a heliocentric distance of 0.58-0.54 AU (pre-perihelion) on 2013 November 16-17. Temporally resolved measurements of the coma distributions of HNC, CH$_3$OH, H$_2$CO and dust were obtained over the course of about an hour on each day. During the period UT 10:10-11… ▽ More Observations of the sungrazing comet C/2012 S1 (ISON) were carried out using the Atacama Large Millimeter/submillimeter Array (ALMA) at a heliocentric distance of 0.58-0.54 AU (pre-perihelion) on 2013 November 16-17. Temporally resolved measurements of the coma distributions of HNC, CH$_3$OH, H$_2$CO and dust were obtained over the course of about an hour on each day. During the period UT 10:10-11:00 on Nov. 16, the comet displayed a remarkable drop in activity, manifested as a $>42$% decline in the molecular line and continuum fluxes. The H$_2$CO observations are consistent with an abrupt, $\approx50$% reduction in the cometary gas production rate soon after the start of our observations. On Nov. 17, the total observed fluxes remained relatively constant during a similar period, but strong variations in the morphology of the HNC distribution were detected as a function of time, indicative of a clumpy, intermittent outflow for this species. Our observations suggest that at least part of the detected HNC originated from degradation of nitrogen-rich organic refractory material, released intermittently from confined regions of the nucleus. By contrast, the distributions of CH$_3$OH and H$_2$CO during the Nov. 17 observations were relatively uniform, consistent with isotropic outflow and stable activity levels for these species. These results highlight a large degree of variability in the production of gas and dust from comet ISON during its pre-perihelion outburst, consistent with repeated disruption of the nucleus interspersed with periods of relative quiescence. △ Less

Submitted 10 February, 2017; originally announced February 2017.

Comments: 9 pages. Accepted for publication in ApJ, February 2017

Thermal physics of the inner coma: ALMA studies of the methanol distribution and excitation in comet C/2012 K1 (PanSTARRS)

Authors: M. A. Cordiner, N. Biver, J. Crovisier, D. Bockelee-Morvan, M. J. Mumma, S. B. Charnley, G. Villanueva, L. Paganini, D. C. Lis, S. N. Milam, A. J. Remijan, I. M. Coulson, Y. -J. Kuan, J. Boissier

Abstract: We present spatially and spectrally-resolved observations of CH$_3$OH emission from comet C/2012 K1 (PanSTARRS) using The Atacama Large Millimeter/submillimeter Array (ALMA) on 2014 June 28-29. Two-dimensional maps of the line-of-sight average rotational temperature ($T_{rot}$) were derived, covering spatial scales $0.3''-1.8''$ (corresponding to sky-projected distances $ρ\sim500$-2500 km). The CH… ▽ More We present spatially and spectrally-resolved observations of CH$_3$OH emission from comet C/2012 K1 (PanSTARRS) using The Atacama Large Millimeter/submillimeter Array (ALMA) on 2014 June 28-29. Two-dimensional maps of the line-of-sight average rotational temperature ($T_{rot}$) were derived, covering spatial scales $0.3''-1.8''$ (corresponding to sky-projected distances $ρ\sim500$-2500 km). The CH$_3$OH column density distributions are consistent with isotropic, uniform outflow from the nucleus, with no evidence for extended sources of CH$_3$OH in the coma. The $T_{rot}(ρ)$ radial profiles show a significant drop within a few thousand kilometers of the nucleus, falling from about 60 K to 20 K between $ρ=0$ and 2500 km on June 28, whereas on June 29, $T_{rot}$ fell from about 120 K to 40 K between $ρ=$ 0 km and 1000 km. The observed $T_{rot}$ behavior is interpreted primarily as a result of variations in the coma kinetic temperature due to adiabatic cooling of the outflowing gas, as well as radiative cooling of the CH$_3$OH rotational levels. Our excitation model shows that radiative cooling is more important for the $J=7-6$ transitions (at 338 GHz) than for the $K=3-2$ transitions (at 252 GHz), resulting in a strongly sub-thermal distribution of levels in the $J=7-6$ band at $ρ\gtrsim1000$ km. For both bands, the observed temperature drop with distance is less steep than predicted by standard coma theoretical models, which suggests the presence of a significant source of heating in addition to the photolytic heat sources usually considered. △ Less

Submitted 28 January, 2017; originally announced January 2017.

Comments: Accepted for publication in ApJ, January 2017

CO Spectral Line Energy Distributions in Galactic Sources: Empirical Interpretation of Extragalactic Observations

Authors: Nick Indriolo, E. A. Bergin, J. R. Goicoechea, J. Cernicharo, M. Gerin, A. Gusdorf, D. C. Lis, P. Schilke

Abstract: The relative populations in rotational transitions of CO can be useful for inferring gas conditions and excitation mechanisms at work in the interstellar medium. We present CO emission lines from rotational transitions observed with Herschel/HIFI in the star-forming cores Orion S, Orion KL, Sgr B2(M), and W49N. Integrated line fluxes from these observations are combined with those from Herschel/PA… ▽ More The relative populations in rotational transitions of CO can be useful for inferring gas conditions and excitation mechanisms at work in the interstellar medium. We present CO emission lines from rotational transitions observed with Herschel/HIFI in the star-forming cores Orion S, Orion KL, Sgr B2(M), and W49N. Integrated line fluxes from these observations are combined with those from Herschel/PACS observations of the same sources to construct CO spectral line energy distributions (SLEDs) from $5\leq J_u\leq48$. These CO SLEDs are compared to those reported in other galaxies, with the intention of empirically determining which mechanisms dominate excitation in such systems. We find that CO SLEDs in Galactic star-forming cores cannot be used to reproduce those observed in other galaxies, although the discrepancies arise primarily as a result of beam filling factors. The much larger regions sampled by the Herschel beams at distances of several Mpc contain significant amounts of cooler gas which dominate the extragalactic CO SLEDs, in contrast to observations of Galactic star-forming regions which are focused specifically on cores containing primarily hot molecular gas. △ Less

Submitted 13 January, 2017; originally announced January 2017.

Comments: Accepted for publication in the Astrophysical Journal; 15 pages, 5 figures, 2 tables

Journal ref: ApJ 836 (2017) 117

Analysis of the Herschel/HEXOS Spectral Survey Towards Orion South: A massive protostellar envelope with strong external irradiation

Authors: K. Tahani, R. Plume, E. A. Bergin, V. Tolls, T. G. Phillips, E. Caux, S. Cabrit, J. R. Goicoechea, P. F. Goldsmith, D. Johnstone, D. C. Lis, L. Pagani, K. M. Menten, H. S. P. Muller, V. Ossenkopf-Okada, J. C. Pearson, F. F. S. van der Tak

Abstract: We present results from a comprehensive submillimeter spectral survey toward the source Orion South, based on data obtained with the HIFI instrument aboard the \textit{Herschel Space Observatory}, covering the frequency range 480 to 1900 GHz. We detect 685 spectral lines with S/N $>$ 3$σ$, originating from 52 different molecular and atomic species. We model each of the detected species assuming co… ▽ More We present results from a comprehensive submillimeter spectral survey toward the source Orion South, based on data obtained with the HIFI instrument aboard the \textit{Herschel Space Observatory}, covering the frequency range 480 to 1900 GHz. We detect 685 spectral lines with S/N $>$ 3$σ$, 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, \& V$_{LSR}$). We find evidence for three different cloud components: a cool (T$_{ex} \sim 20-40$ K), spatially extended ($> 60"$), and quiescent ($ΔV_{FWHM} \sim 4$ km s $^{-1}$) component; a warmer (T$_{ex} \sim 80-100$ K), less spatially extended ($\sim 30"$), and dynamic ($ΔV_{FWHM} \sim 8$ km s $^{-1}$) component, which is likely affected by embedded outflows; and a kinematically distinct region (T$_{ex}$ $>$ 100 K; V$_{LSR}$ $\sim$ 8 km s $^{-1}$), dominated by emission from species which 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 \textit{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. △ Less

Submitted 22 August, 2016; originally announced August 2016.

Journal ref: ApJ 832, 12 (2016)