Showing 1–39 of 39 results for author: Montillaud, J

The magnetic field in colliding filaments G202.3+2.5

Authors: Qi-Lao Gu, Tie Liu, Pak Shing Li, Zhi-Qiang Shen, Xunchuan Liu, Junhao Liu, Xing Lu, Julien Montillaud, Sihan Jiao, Mika Juvela, Mark G. Rawlings, Qizhou Zhang, Patrick Koch, Isabelle Ristorcelli, Jean-Sébastien Carriere, David Eden, Zhiyuan Ren, Ken'ichi Tatematsu, Naomi Hirano, Qiu-yi Luo, Xiaofeng Mai, Namitha Issac

Abstract: We observe the magnetic field morphology towards a nearby star-forming filamentary cloud, G202.3+2.5, by the JCMT/POL-2 850 μm thermal dust polarization observation with an angular resolution of 14.4" (~0.053 pc). The average magnetic field orientation is found to be perpendicular to the filaments while showing different behaviors in the four subregions, suggesting various effects from filaments'… ▽ More We observe the magnetic field morphology towards a nearby star-forming filamentary cloud, G202.3+2.5, by the JCMT/POL-2 850 μm thermal dust polarization observation with an angular resolution of 14.4" (~0.053 pc). The average magnetic field orientation is found to be perpendicular to the filaments while showing different behaviors in the four subregions, suggesting various effects from filaments' collision in these subregions. With the kinematics obtained by N2H+ observation by IRAM, we estimate the plane-of-sky (POS) magnetic field strength by two methods, the classical Davis-Chandrasekhar-Fermi (DCF) method and the angular dispersion function (ADF) method, B_{pos,dcf} and B_{pos,adf} are ~90 μG and ~53 μG. We study the relative importance between the gravity (G), magnetic field (B) and turbulence (T) in the four subregions, find G > T > B, G >= T > B, G ~ T > B and T > G > B in the north tail, west trunk, south root and east wing, respectively. In addition, we investigate the projection effect on the DCF and ADF methods based on a similar simulation case and find the 3D magnetic field strength may be underestimated by a factor of ~3 if applying the widely-used statistical B_{pos}-to-B_{3D} factor when using DCF or ADF method, which may further underestimate/overestimate related parameters. △ Less

Submitted 12 January, 2024; v1 submitted 10 January, 2024; originally announced January 2024.

Comments: Accepted by ApJ. 20 pages, 9 figures

A self-consistent dynamical model of the Milky Way disc adjusted to Gaia data

Authors: A. C. Robin, O. Bienaymé, J. B. Salomon, C. Reylé, N. Lagarde, F. Figueras, R. Mor, J. G. Fernández-Trincado, J. Montillaud

Abstract: This paper shows how a self-consistent dynamical model can be obtained by fitting the gravitational potential of the Milky Way to the stellar kinematics and densities from Gaia data. Using the Besancon Galaxy Model we derive a potential and the disc stellar distribution functions are computed based on three integrals of motion to model stationary stellar discs. The gravitational potential and the… ▽ More This paper shows how a self-consistent dynamical model can be obtained by fitting the gravitational potential of the Milky Way to the stellar kinematics and densities from Gaia data. Using the Besancon Galaxy Model we derive a potential and the disc stellar distribution functions are computed based on three integrals of motion to model stationary stellar discs. The gravitational potential and the stellar distribution functions are built self-consistently, and then adjusted to be in agreement with the kinematics and the density distributions obtained from Gaia observations. A Markov chain Monte Carlo (MCMC) is used to fit the free parameters of the dynamical model to Gaia parallax and proper motion distributions. The fit is done on several sets of Gaia eDR3 data, widely spread in longitudes and latitudes. We are able to determine the velocity dispersion ellipsoid and its tilt for sub-components of different ages, both varying with R and z. The density laws and their radial scale lengths, for the thin and thick disc populations are also obtained self-consistently. This new model has some interesting characteristics, such as a flaring thin disc. The thick disc is found to present very distinctive characteristics from the old thin disc, both in density and kinematics. This well supports the idea that thin and thick discs were formed in distinct scenarios as the density and kinematics transition between them is found to be abrupt. The dark matter halo is shown to be nearly spherical. We also derive the Solar motion to be (10.79 $\pm$ 0.56, 11.06 $\pm$ 0.94, 7.66 $\pm$ 0.43) km/s, in good agreement with recent studies. The resulting fully self-consistent gravitational potential, still axisymmetric, is a good approximation of a smooth mass distribution in the Milky Way and can be used for further studies, including to compute orbits for real stars in our Galaxy (abridged). △ Less

Submitted 13 September, 2022; v1 submitted 29 August, 2022; originally announced August 2022.

Comments: Final version, 26 pages, 14 figures. A&A in press

Journal ref: A&A 667, A98 (2022)

Chemical exploration of Galactic cold cores

Authors: Chenlin Zhou, Charlotte Vastel, Julien Montillaud, Cecilia Ceccarelli, Karine Demyk, Jorma Harju, Mika Juvela, Isabelle Ristorcelli, Tie Liu

Abstract: A solar-type system starts from an initial molecular core that acquires organic complexity as it evolves. The so-called prestellar cores that can be studied are rare, which has hampered our understanding of how organic chemistry sets in and grows. Aims. We selected the best prestellar core targets from the cold core catalogue that represent a diversity in terms of their environment to explore thei… ▽ More A solar-type system starts from an initial molecular core that acquires organic complexity as it evolves. The so-called prestellar cores that can be studied are rare, which has hampered our understanding of how organic chemistry sets in and grows. Aims. We selected the best prestellar core targets from the cold core catalogue that represent a diversity in terms of their environment to explore their chemical complexity: 1390 (in the compressed shell of Lambda Ori), 869 (in the MBM12 cloud), and 4149 (in the California nebula). We obtained a spectral survey with the IRAM 30 m telescope in order to explore the molecular complexity of the cores. We carried out a radiative transfer analysis of the detected transitions in order to place some constraints on the physical conditions of the cores and on the molecular column densities. We also used the molecular ions in the survey to estimate the cosmic-ray ionisation rate and the S/H initial elemental abundance using a gas-phase chemical model to reproduce their abundances. We found large differences in the molecular complexity (deuteration, complex organic molecules, sulphur, carbon chains, and ions) and compared their chemical properties with a cold core and two prestellar cores. The chemical diversity we found in the three cores seems to be correlated with their chemical evolution: two of them are prestellar (1390 and 4149), and one is in an earlier stage (869). The influence of the environment is likely limited because cold cores are strongly shielded from their surroundings. The high extinction prevents interstellar UV radiation from penetrating deeply into the cores. Higher spatial resolution observations of the cores are therefore needed to constrain the physical structure of the cores, as well as a larger-scale distribution of molecular ions to understand the influence of the environment on their molecular complexity. △ Less

Submitted 4 February, 2022; originally announced February 2022.

Comments: 26 pages, 13 figures, accepted in A&A in January 2022

Journal ref: A&A 658, A131 (2022)

3D extinction mapping of the Milky Way using Convolutional Neural Networks: Presentation of the method and demonstration in the Carina Arm region

Authors: D. Cornu, J. Montillaud, D. J. Marshall, A. C. Robin, L. Cambrésy

Abstract: Context. Several methods have been proposed to build 3D extinction maps of the Milky Way (MW), most often based on Bayesian approaches. Although some studies employed machine learning (ML) methods in part of their procedure, or to specific targets, no 3D extinction map of a large volume of the MW solely based on a Neural Network method has been reported so far. Aims. We aim to apply deep learning… ▽ More Context. Several methods have been proposed to build 3D extinction maps of the Milky Way (MW), most often based on Bayesian approaches. Although some studies employed machine learning (ML) methods in part of their procedure, or to specific targets, no 3D extinction map of a large volume of the MW solely based on a Neural Network method has been reported so far. Aims. We aim to apply deep learning as a solution to build 3D extinction maps of the MW. Methods. We built a convolutional neural network (CNN) using the CIANNA framework, and trained it with synthetic 2MASS data. We used the Besançon Galaxy model to generate mock star catalogs, and 1D Gaussian random fields to simulate the extinction profiles. From these data we computed color-magnitude diagrams (CMDs) to train the network, using the corresponding extinction profiles as targets. A forward pass with observed 2MASS CMDs provided extinction profile estimates for a grid of lines of sight. Results. We trained our network with data simulating lines of sight in the area of the Carina spiral arm tangent and obtained a 3D extinction map for a large sector in this region ($l = 257 - 303$ deg, $|b| \le 5$ deg), with distance and angular resolutions of $100$ pc and $30$ arcmin, respectively, and reaching up to $\sim 10$ kpc. Although each sightline is computed independently in the forward phase, the so-called fingers-of-God artifacts are weaker than in many other 3D extinction maps. We found that our CNN was efficient in taking advantage of redundancy across lines of sight, enabling us to train it with only 9 sightlines simultaneously to build the whole map. Conclusions. We found deep learning to be a reliable approach to produce 3D extinction maps from large surveys. With this methodology, we expect to easily combine heterogeneous surveys without cross-matching, and therefore to exploit several surveys in a complementary fashion. △ Less

Submitted 14 January, 2022; originally announced January 2022.

Comments: 7 pages

Deciphering the evolution of the Milky Way discs: Gaia APOGEE Kepler giant stars and the Besançon Galaxy Model

Authors: N. Lagarde, C. Reylé, C. Chiappini, R. Mor, F. Anders, F. Figueras, A. Miglio, M. Romero-Gómez, T. Antoja, N. Cabral, J. -B. Salomon, A. C. Robin, O. Bienaymé, C. Soubiran, D. Cornu, J. Montillaud

Abstract: We investigate the properties of the double sequences of the Milky Way discs visible in the [$α$/Fe] vs [Fe/H] diagram. In the framework of Galactic formation and evolution, we discuss the complex relationships between age, metallicity, [$α$/Fe], and the velocity components. We study stars with measured chemical, seismic and astrometric properties from the APOGEE survey, the Kepler and Gaia satell… ▽ More We investigate the properties of the double sequences of the Milky Way discs visible in the [$α$/Fe] vs [Fe/H] diagram. In the framework of Galactic formation and evolution, we discuss the complex relationships between age, metallicity, [$α$/Fe], and the velocity components. We study stars with measured chemical, seismic and astrometric properties from the APOGEE survey, the Kepler and Gaia satellites, respectively. We separate the [$α$/Fe]-[Fe/H] diagram into 3 stellar populations: the thin disc, the high-$α$ metal-poor thick disc and the high-$α$ metal-rich thick disc and characterise each of these in the age-chemo-kinematics parameter space. We compare results obtained from different APOGEE data releases and using two recent age determinations. We use the Besançon Galaxy model (BGM) to highlight selection biases and mechanisms not included in the model. The thin disc exhibits a flat age-metallicity relation while [$α$/Fe] increases with stellar age. We confirm no correlation between radial and vertical velocities with [Fe/H], [$α$/Fe] and age for each stellar population. Considering both samples, V$_\varphi$ decreases with age for the thin disc, while it increases with age for the h$α$mp thick disc. Although the age distribution of the h$α$mr thick disc is very close to that of the h$α$mp thick disc between 7 and 14 Gyr, its kinematics seems to follow that of the thin disc. This feature, not predicted by the hypotheses included in the BGM, suggests a different origin and history for this population. Finally, we show that there is a maximum dispersion of the vertical velocity, $σ_Z$, with age for the h$α$mp thick disc around 8 Gyr. The comparisons with the BGM simulations suggest a more complex chemo-dynamical scheme to explain this feature, most likely including mergers and radial migration effects △ Less

Submitted 27 July, 2021; v1 submitted 29 June, 2021; originally announced June 2021.

Comments: Accepted for publication in Aastronomy & Astrophysics on 28/06/2021, 19 pages, 16 figures

Journal ref: A&A 654, A13 (2021)

Characterization of dense Planck clumps observed with Herschel and SCUBA-2

Authors: E. Mannfors, M. Juvela, L. Bronfman, D. J. Eden, Jinhua He, Gwanjeon Kim, Kee-Tae Kim, H. Kirppu, T. Liu, J. Montillaud, H. Parsons, Patricio Sanhueza, Hsien Shang, A. Soam, K. Tatematsu, A. Traficante, M. S. Väisälä, Chang Won Lee

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

Submitted 18 June, 2021; originally announced June 2021.

Comments: 60 pages (article 18 pages), 73 figures To be published in A&A

Journal ref: A&A 654, A123 (2021)

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): Detection of extremely high density compact structure of prestellar cores and multiple substructures within

Authors: Dipen Sahu, Sheng-Yuan Liu, Tie Liu, Neal J. Evans II, Naomi Hirano, Ken'ichi Tatematsu, Chin-Fei Lee, Kee-Tae Kim, Somnath Dutta, Dana Alina, Leonardo Bronfman, Maria Cunningham, David J. Eden, Guido Garay, Paul F. Goldsmith, Jinhua He, Shih-Ying Hsu, Kai-Syun Jhan, Doug Johnstone, Mika Juvela, Gwanjeong Kim, Yi-Jehng Kuan, Woojin Kwon, Chang Won Lee, Jeong-Eun Lee , et al. (16 additional authors not shown)

Abstract: Prestellar cores are self-gravitating dense and cold structures within molecular clouds where future stars are born. They are expected, at the stage of transitioning to the protostellar phase, to harbor centrally concentrated dense (sub)structures that will seed the formation of a new star or the binary/multiple stellar systems. Characterizing this critical stage of evolution is key to our underst… ▽ More Prestellar cores are self-gravitating dense and cold structures within molecular clouds where future stars are born. They are expected, at the stage of transitioning to the protostellar phase, to harbor centrally concentrated dense (sub)structures that will seed the formation of a new star or the binary/multiple stellar systems. Characterizing this critical stage of evolution is key to our understanding of star formation. In this work, we report the detection of high density (sub)structures on the thousand-au scale in a sample of dense prestellar cores. Through our recent ALMA observations towards the Orion molecular cloud, we have found five extremely dense prestellar cores, which have centrally concentrated regions $\sim$ 2000 au in size, and several $10^7$ $cm^{-3}$ in average density. Masses of these centrally dense regions are in the range of 0.30 to 6.89 M$_\odot$. {\it For the first time}, our higher resolution observations (0.8$'' \sim $ 320 au) further reveal that one of the cores shows clear signatures of fragmentation; such individual substructures/fragments have sizes of 800 -1700 au, masses of 0.08 to 0.84 M$_\odot$, densities of $2 - 8\times 10^7$ $cm^{-3}$ and separations of $\sim 1200$ au. The substructures are massive enough ($\gtrsim 0.1~M_\odot$) to form young stellar objects and are likely examples of the earliest stage of stellar embryos which can lead to widely ($\sim$ 1200 au) separated multiple systems. △ Less

Submitted 12 February, 2023; v1 submitted 15 December, 2020; originally announced December 2020.

Comments: Published in The Astrophysical Journal Letters (ApJL)

A neural network-based methodology to select young stellar object candidates from IR surveys

Authors: David Cornu, Julien Montillaud

Abstract: Observed Young Stellar Objects (YSOs) are used to study star formation and characterize star forming regions. For this purpose, YSO candidate catalogs are compiled from various surveys, especially in the infrared (IR), and simple selection schemes in colour-magnitude diagrams (CMDs) are often used to identify and classify YSOs. We propose a methodology for YSO classification through Machine Learni… ▽ More Observed Young Stellar Objects (YSOs) are used to study star formation and characterize star forming regions. For this purpose, YSO candidate catalogs are compiled from various surveys, especially in the infrared (IR), and simple selection schemes in colour-magnitude diagrams (CMDs) are often used to identify and classify YSOs. We propose a methodology for YSO classification through Machine Learning (ML) using Spitzer IR data. We detail our approach in order to ensure reproducibility and provide an in-depth example on how to efficiently apply ML to an astrophysical classification. We used feedforward Artificial Neural Networks (ANNs) that use the four IRAC bands ($3.6, 4.5, 5.8$ and $8 μm$) and the $24\ μm$ MIPS band from Spitzer to classify point source objects into CI and CII YSO candidates or as contaminants. We found that ANNs can efficiently be applied to YSO classification with a contained number of neurons ($\sim$ 25). Knowledge gathered on one star-forming region has shown to be partly efficient for prediction in new regions. Best generalization capacity was achieved using a combination of several star-forming regions to train the network. Carefully rebalancing the training proportions was necessary to achieve good results. We achieved above 90% and 97% recovery rate for CI and CII YSOs, respectively, with precision above 80% and 90% for our most general result. We took advantage of the ANN great flexibility to define, for each object, an effective membership probability to each output class. Using a threshold in this probability was found to efficiently improve the classification results at a reasonable cost of object exclusion. With this selection, we reached 90% precision on CI YSOs, for more than half of them. Our catalog of YSO candidates in Orion (365 CI, 2381 CII) and NGC 2264 (101 CI, 469 CII) predicted by our final ANN is publicly available at CDS. △ Less

Submitted 28 January, 2021; v1 submitted 4 October, 2020; originally announced October 2020.

Journal ref: A&A 647, A116 (2021)

Large-scale magnetic field in the Monoceros OB-1 East molecular cloud

Authors: D. Alina, J. Montillaud, Y. Hu, A. Lazarian, I. Ristorcelli, E. Abdikamalov, S. Sagynbayeva, M. Juvela, T. Liu, J. -S. Carrière

Abstract: We study the large-scale magnetic field structure and its interplay with the gas dynamics in the Monoceros OB1 East molecular cloud. We combine observations of dust polarised emission from the Planck telescope and CO molecular line emission observations from the Taeduk Radio Astronomy Observatory 14-metre telescope. We calculate the strength of the plane-of-the-sky magnetic field using a modified… ▽ More We study the large-scale magnetic field structure and its interplay with the gas dynamics in the Monoceros OB1 East molecular cloud. We combine observations of dust polarised emission from the Planck telescope and CO molecular line emission observations from the Taeduk Radio Astronomy Observatory 14-metre telescope. We calculate the strength of the plane-of-the-sky magnetic field using a modified Chandrasekhar-Fermi method and estimate mass over flux ratios in different regions of the cloud. We use the comparison of the velocity and intensity gradients of the molecular line observations with the polarimetric observations to trace dynamically active regions. The molecular complex shows an ordered large-scale plane-of-the-sky magnetic field structure. In the Northern part, it is mostly orientated along the filamentary structures while the Southern part shows at least two regions with distinct magnetic field orientations. We find that in the Northern filaments the magnetic field is unlikely to provide support against fragmentation at large scales. Our analysis reveals a shock region in the Northern part of the complex right in-between two filamentary clouds which were previously suggested to be in collision. Moreover, the shock seems to extend farther towards the Western part of the complex. In the Southern part, we find that either the magnetic field guides the accretion of interstellar matter towards the cloud or it was dragged by the matter towards the densest regions. The large-scale magnetic field in Monoceros OB-1 East molecular clouds is tightly connected to the global structure of the complex and, in the Northern part, it seems to be dominated by gravity and turbulence, while in the Southern part it influences the structuring of matter. △ Less

Submitted 6 October, 2021; v1 submitted 30 July, 2020; originally announced July 2020.

Comments: 20 pages, 21 figures Submitted to Astronomy & Astrophysics

Journal ref: A&A 658, A90 (2022)

Space Project for Astrophysical and Cosmological Exploration (SPACE), an ESA stand-alone mission and a possible contribution to the Origins Space Telescope

Authors: Denis Burgarella, Andrew Bunker, Rychard Bouwens, Laurent Pagani, Jose Afonso, Hakim Atek, Marc Audard, Sylvie Cabrit, Karina Caputi, Laure Ciesla, Christopher Conselice, Asantha Cooray, Giovanni Cresci, Mirko Curti, Jose Miguel Rodriguez Espinosa, Marc Ferrari, Chiaki Kobayashi, Nadege Lagarde, Jesus Gallego Maestro, Roberto Maiolino, Katarzyna Malek, Filippo Mannucci, Julien Montillaud, Pascal Oesch, Chris Pearson , et al. (8 additional authors not shown)

Abstract: We propose a new mission called Space Project for Astrophysical and Cosmological Exploration (SPACE) as part on the ESA long term planning Voyage 2050 programme. SPACE will study galaxy evolution at the earliest times, with the key goals of charting the formation of the heavy elements, measuring the evolution of the galaxy luminosity function, tracing the build-up of stellar mass in galaxies over… ▽ More We propose a new mission called Space Project for Astrophysical and Cosmological Exploration (SPACE) as part on the ESA long term planning Voyage 2050 programme. SPACE will study galaxy evolution at the earliest times, with the key goals of charting the formation of the heavy elements, measuring the evolution of the galaxy luminosity function, tracing the build-up of stellar mass in galaxies over cosmic time, and finding the first super-massive black holes (SMBHs) to form. The mission will exploit a unique region of the parameter space, between the narrow ultra-deep surveys with HST and JWST, and shallow wide-field surveys such as Roman Space Telescope and EUCLID, and should yield by far the largest sample of any current or planned mission of very high redshift galaxies at z > 10 which are sufficiently bright for detailed follow-up spectroscopy. Crucially, we propose a wide-field spectroscopic near-IR + mid-IR capability which will greatly enhance our understanding of the first galaxies by detecting and identifying a statistical sample of the first galaxies and the first SMBH, and to chart the metal enrichment history of galaxies in the early Universe - potentially finding signatures of the very first stars to form from metal-free primordial gas. The wide-field and wavelength range of SPACE will also provide us a unique opportunity to study star formation by performing a wide survey of the Milky Way in the near-IR + mid-IR. This science project can be enabled either by a stand-alone ESA-led M mission or by an instrument for an L mission (with ESA and/or NASA, JAXA and other international space agencies) with a wide-field (sub-)millimetre capability at wavelength > 500 microns. △ Less

Submitted 5 August, 2020; v1 submitted 29 July, 2020; originally announced July 2020.

Comments: Submitted to Experimental Astronomy as one of the proposals for the European Space Agency Voyage 2050 Programme

Multi-scale analysis of the Monoceros OB 1 star-forming region: II. Colliding filaments in the Monoceros OB1 molecular cloud

Authors: Julien Montillaud, Mika Juvela, Charlotte Vastel, J. H. He, Tie Liu, Isabelle Ristorcelli, David J. Eden, Sung-ju Kang, Kee-Tae Kim, Patrick M. Koch, Chang Won Lee, Mark G. Rawlings, Mika Saajasto, Patricio Sanhueza, Archana Soam, Sarolta Zahorecz, Dana Alina, Rebeka Bögner, David Cornu, Yasuo Doi, Johanna Malinen, Douglas Marshall, E. R. Micelotta, V. M. Pelkonen, L. V. Tóth , et al. (2 additional authors not shown)

Abstract: We started a multi-scale analysis of G202.3+2.5, an intertwined filamentary region of Monoceros OB1. In Paper I, we examined the distributions of dense cores and protostars and found enhanced star formation (SF) activity in the junction region of the filaments. In this second paper, we aim to unveil the connections between the core and filament evolutions, and between the filament dynamics and the… ▽ More We started a multi-scale analysis of G202.3+2.5, an intertwined filamentary region of Monoceros OB1. In Paper I, we examined the distributions of dense cores and protostars and found enhanced star formation (SF) activity in the junction region of the filaments. In this second paper, we aim to unveil the connections between the core and filament evolutions, and between the filament dynamics and the global evolution of the cloud. We characterise the gas dynamics and energy balance using Herschel and WISE observations and molecular tracers observed with the IRAM 30m and TRAO 14m telescopes. The velocity field of the cloud is examined and velocity-coherent structures are put in perspective with the cloud environment. Two main velocity components (VCs) are revealed, well separated in the north and merged around the location of intense N2H+ emission where Paper I found the peak of SF activity. The relative position of the two VCs along the sightline, and the velocity gradient in N2H+ emission imply that the VCs have been undergoing collision for ~10^5 yrs. The dense gas where N2H+ is detected is interpreted as the compressed region between the two filaments, which corresponds to a high mass inflow rate of ~1e-3 Msun/yr and possibly leads to an increase in its SF efficiency. We identify a protostar in the junction region that possibly powers two crossed intermittent outflows. We show that the HII region around the nearby cluster NCG 2264 is still expanding and its role in the collision is examined. However, we cannot rule out the idea that the collision arises mostly from the global collapse of the cloud. The (sub-)filament-scale observables examined in this paper reveal a collision between G202.3+2.5 sub-structures and its probable role in feeding the cores in the junction region. One must now characterise the cloud morphology, its fragmentation, and magnetic field, all at high resolution. △ Less

Submitted 13 September, 2019; originally announced September 2019.

Comments: A&A, in press; IRAM and TRAO data cubes will be available at CDS

Journal ref: A&A 631, A3 (2019)

Multi-scale analysis of the Monoceros OB 1 star-forming region: I. The dense core population

Authors: Julien Montillaud, Mika Juvela, Charlotte Vastel, J. H. He, Tie Liu, Isabelle Ristorcelli, David J. Eden, Sung-ju Kang, Kee-Tae Kim, Patrick M. Koch, Chang Won Lee, Mark G. Rawlings, Mika Saajasto, Patricio Sanhueza, Archana Soam, Sarolta Zahorecz, Dana Alina, Rebeka Bögner, David Cornu, Yasuo Doi, Johanna Malinen, Douglas Marshall, E. R. Micelotta, V. M. Pelkonen, L. V. Tóth , et al. (2 additional authors not shown)

Abstract: Current theories and models attempt to explain star formation globally, from core scales to giant molecular cloud scales. A multi-scale observational characterisation of an entire molecular complex is necessary to constrain them. We investigate star formation in G202.3+2.5, a ~10x3 pc sub-region of the Monoceros OB1 cloud with a complex morphology harbouring interconnected filamentary structures.… ▽ More Current theories and models attempt to explain star formation globally, from core scales to giant molecular cloud scales. A multi-scale observational characterisation of an entire molecular complex is necessary to constrain them. We investigate star formation in G202.3+2.5, a ~10x3 pc sub-region of the Monoceros OB1 cloud with a complex morphology harbouring interconnected filamentary structures. We aim to connect the evolution of cores and filaments in G202.3+2.5 with the global evolution of the cloud and to identify the engines of the cloud dynamics. In this first paper, the star formation activity is evaluated by surveying the distributions of dense cores and protostars, and their evolutionary state, as characterised using both infrared observations from the Herschel and WISE telescopes and molecular line observations with the IRAM 30-m telescope. We find ongoing star formation in the whole cloud, with a local peak in star formation activity around the centre of G202.3+2.5 where a chain of massive cores (10-50 Msun) forms a massive ridge (>150 Msun). All evolutionary stages from starless cores to Class II protostars are found in G202.3+2.5, including a possibly starless, large column density (8x10^{22} cm-2), and massive (52 Msun) core. All the core-scale observables examined in this paper point to an enhanced star formation activity centred on the junction between the three main branches of the ramified structure of G202.3+2.5. This suggests that the increased star-formation activity results from the convergence of these branches. To further investigate the origin of this enhancement, it is now necessary to extend the analysis to larger scales, in order to examine the relationship between cores, filaments and their environment. We address these points through the analysis of the dynamics of G202.3+2.5 in a joint paper. △ Less

Submitted 9 September, 2019; originally announced September 2019.

Journal ref: A&A 631, L1 (2019)

Galactic continuum surveys with the new 850 micron MKID camera at the EAO/JCMT 15-m telescope

Authors: Tie Liu, Di Li, David Eden, James Di Francesco, JinHua He, Sarah I. Sadavoy, Kenichi Tatematsu, Kee-Tae Kim, Naomi Hirano, Gary Fuller, Yuefang Wu, Sheng-Yuan Liu, Ke Wang, Mark Thompson, Mika Juvela, Isabelle Ristorcelli, L. Viktor Toth, Archana Soam, Patricio Sanhueza, Sung-ju Kang, Woojin Kwon, D. Ward-Thompson, Julien Montillaud, John Richer, Jinjin Xie , et al. (2 additional authors not shown)

Abstract: This white paper gives a brief summary of Galactic continuum surveys with the next generation 850 micron camera at the James Clerk Maxwell Telescope (JCMT) in the next decade. This new camera will have mapping speeds at least 10-20 times faster than the present SCUBA-2 camera, and will enable deep (<10 mJy/beam) and extensive continuum surveys such as a wider (-5 deg.<l<240 deg., |b|<2-5 deg.) Gal… ▽ More This white paper gives a brief summary of Galactic continuum surveys with the next generation 850 micron camera at the James Clerk Maxwell Telescope (JCMT) in the next decade. This new camera will have mapping speeds at least 10-20 times faster than the present SCUBA-2 camera, and will enable deep (<10 mJy/beam) and extensive continuum surveys such as a wider (-5 deg.<l<240 deg., |b|<2-5 deg.) Galactic Plane survey, a larger Gould Belt survey, and a follow-up survey of all Planck compact objects visible from the northern hemisphere. These surveys will provide a complete census of molecular clouds, filaments, and dense cores across the Galaxy, vital for studying star formation in various environments. △ Less

Submitted 19 July, 2019; originally announced July 2019.

Comments: white paper for Galactic continuum surveys with the next generation 850 micron camera at JCMT. Comments are very welcome

Submillimeter continuum variability in Planck Galactic cold clumps

Authors: Geumsook Park, Kee-Tae Kim, Doug Johnstone, Sung-ju Kang, Tie Liu, Steve Mairs, Minho Choi, Jeong-Eun Lee, Patricio Sanhueza, Mika Juvela, Miju Kang, David Eden, Archana Soam, Julien Montillaud, Gary Fuller, Patrick M. Koch, Chang Won Lee, Dimitris Stamatellos, Jonathan Rawlings, Gwanjeong Kim, Chuan-Peng Zhang, Woojin Kwon, Hyunju Yoo

Abstract: In the early stages of star formation, a protostar is deeply embedded in an optically thick envelope such that it is not directly observable. Variations in the protostellar accretion rate, however, will cause luminosity changes that are reprocessed by the surrounding envelope and are observable at submillimeter wavelengths. We searched for submillimeter flux variability toward 12 Planck Galactic C… ▽ More In the early stages of star formation, a protostar is deeply embedded in an optically thick envelope such that it is not directly observable. Variations in the protostellar accretion rate, however, will cause luminosity changes that are reprocessed by the surrounding envelope and are observable at submillimeter wavelengths. We searched for submillimeter flux variability toward 12 Planck Galactic Cold Clumps detected by the James Clerk Maxwell Telescope (JCMT)-SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE) survey. These observations were conducted at 850 um using the JCMT/SCUBA-2. Each field was observed three times over about 14 months between 2016 April and 2017 June. We applied a relative flux calibration and achieved a calibration uncertainty of ~ 3.6% on average. We identified 136 clumps across 12 fields and detected four sources with flux variations of ~ 30%. For three of these sources, the variations appear to be primarily due to large-scale contamination, leaving one plausible candidate. The flux change of the candidate may be associated with low- or intermediate-mass star formation assuming a distance of 1.5 kpc, although we cannot completely rule out the possibility that it is a random deviation. Further studies with dedicated monitoring would provide a better understanding of the detailed relationship between submillimeter flux and accretion rate variabilities while enhancing the search for variability in star-forming clumps farther away than the Gould Belt. △ Less

Submitted 14 June, 2019; v1 submitted 28 May, 2019; originally announced May 2019.

Comments: 14 pages, 7 figures, 3 tables; accepted for publication in ApJS, LaTeX; minor corrections (typos, acronyms, etc.), update of authors' information

SCOPE: SCUBA-2 Continuum Observations of Pre-protostellar Evolution - Survey Description and Compact Source Catalogue

Authors: D. J. Eden, Tie Liu, Kee-Tae Kim, S. -Y. Liu, K. Tatematsu, J. Di Francesco, K. Wang, Y. Wu, M. A. Thompson, G. A. Fuller, Di Li, I. Ristorcelli, Sung-ju Kang, N. Hirano, D. Johnstone, Y. Lin, J. H. He, P. M. Koch, Patricio Sanhueza, S. -L. Qin, Q. Zhang, P. F. Goldsmith, N. J. Evans II, J. Yuan, C. -P. Zhang , et al. (136 additional authors not shown)

Abstract: We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14… ▽ More We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other. △ Less

Submitted 26 February, 2019; originally announced February 2019.

Comments: 16 pages, 13 figures, 1 table. Accepted for publication in MNRAS

Dust spectrum and polarisation at 850 um in the massive IRDC G035.39-00.33

Authors: Mika Juvela, Vincent Guillet, Tie Liu, Isabelle Ristorcelli, Veli-Matti Pelkonen, Dana Alina, Leonardo Bronfman, David J. Eden, Kee Tae Kim, Patrick M. Koch, Woojin Kwon, Chang Won Lee, Johanna Malinen, Elisabetta Micelotta, Julien Montillaud, Mark G. Rawlings, Patricio Sanhueza, Archana Soam, Alessio Traficante, Nathalie Ysard, Chuan-Peng Zhang

Abstract: The dust sub-millimetre polarisation of star-forming clouds carries information on dust and the role of magnetic fields in cloud evolution. With observations of a dense filamentary cloud G035.39-00.33, we aim to characterise the dust emission properties and the variations of the polarisation fraction. JCMT SCUBA-2/POL-2 data at 850um are combined with Planck 850um (353GHz) data to map polarisation… ▽ More The dust sub-millimetre polarisation of star-forming clouds carries information on dust and the role of magnetic fields in cloud evolution. With observations of a dense filamentary cloud G035.39-00.33, we aim to characterise the dust emission properties and the variations of the polarisation fraction. JCMT SCUBA-2/POL-2 data at 850um are combined with Planck 850um (353GHz) data to map polarisation fractions. With previous SCUBA-2 observations (450um and 850um) and Herschel data, the column densities are determined via modified blackbody fits and via radiative transfer modelling. Models are constructed to examine how the polarisation angles and fractions depend on potential magnetic field geometries and grain alignment. POL-2 data show clear changes in the magnetic field orientation. The filament has a peak column density of N(H2)~7 10^22 cm-2, a minimum dust temperature of T~12 K, and a mass of some 4300Msun for the area N(H2)> 5 10^21 cm-2. The estimated average value of the dust opacity spectral index is beta ~ 1.9. The ratio of sub-millimetre and J band optical depths is tau(250 um)/tau(J) ~ 2.5 10^-3, more than four times the typical values for diffuse medium. The polarisation fraction decreases as a function of column density to p ~ 1% in the central filament. Because of noise, the observed decrease of p(N) is significant only at N(H2)>2 10^22 cm-2. The observations suggest that the grain alignment is not constant. Although the data can be explained with a complete loss of alignment at densities above ~ 10^4 cm-3 or using the predictions of radiative torques alignment, the uncertainty of the field geometry and the spatial filtering of the SCUBA-2 data prevent strong conclusions. G035.39-00.33 shows strong signs of dust evolution and the low polarisation fraction is suggestive of a loss of polarised emission from its densest parts. △ Less

Submitted 6 September, 2018; v1 submitted 4 September, 2018; originally announced September 2018.

Comments: accepted for publication in Astronomy and Astrophysics

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

The properties of Planck Galactic cold clumps in the L1495 dark cloud

Authors: Mengyao Tang, Tie Liu, Sheng-Li Qin, Kee-Tae Kim, Yuefang Wu, Ken'ichi Tatematsu, Jinghua Yuan, Ke Wang, Harriet Parsons, Patrick M. Koch, Patricio Sanhueza, D. Ward-Thompson, L. Viktor Tóth, Archana Soam, Chang Won Lee, David Eden, James Di Francesco, Jonathan Rawlings, Mark G. Rawlings, Julien Montillaud, Chuan-Peng Zhang, M. R. Cunningham

Abstract: Planck Galactic Cold Clumps (PGCCs) possibly represent the early stages of star formation. To understand better the properties of PGCCs, we studied 16 PGCCs in the L1495 cloud with molecular lines and continuum data from Herschel, JCMT/SCUBA-2 and the PMO 13.7 m telescope. Thirty dense cores were identified in 16 PGCCs from 2-D Gaussian fitting. The dense cores have dust temperatures of… ▽ More Planck Galactic Cold Clumps (PGCCs) possibly represent the early stages of star formation. To understand better the properties of PGCCs, we studied 16 PGCCs in the L1495 cloud with molecular lines and continuum data from Herschel, JCMT/SCUBA-2 and the PMO 13.7 m telescope. Thirty dense cores were identified in 16 PGCCs from 2-D Gaussian fitting. The dense cores have dust temperatures of $T_{\rm d}$ = 11-14 K, and H$_{2}$ column densities of $N_{\rm H_{2}}$ = 0.36-2.5$\times10^{22}$ cm$^{-2}$. We found that not all PGCCs contain prestellar objects. In general, the dense cores in PGCCs are usually at their earliest evolutionary stages. All the dense cores have non-thermal velocity dispersions larger than the thermal velocity dispersions from molecular line data, suggesting that the dense cores may be turbulence-dominated. We have calculated the virial parameter $α$ and found that 14 of the dense cores have $α$ $<$ 2, while 16 of the dense cores have $α$ $>$ 2. This suggests that some of the dense cores are not bound in the absence of external pressure and magnetic fields. The column density profiles of dense cores were fitted. The sizes of the flat regions and core radii decrease with the evolution of dense cores. CO depletion was found to occur in all the dense cores, but is more significant in prestellar core candidates than in protostellar or starless cores. The protostellar cores inside the PGCCs are still at a very early evolutionary stage, sharing similar physical and chemical properties with the prestellar core candidates. △ Less

Submitted 14 February, 2018; originally announced February 2018.

Structure of photodissociation fronts in star-forming regions revealed by observations of high-J CO emission lines with Herschel

Authors: C. Joblin, E. Bron, C. Pinto, P. Pilleri, F. Le Petit, M. Gerin, J. Le Bourlot, A. Fuente, O. Berne, J. R. Goicoechea, E. Habart, M. Koehler, D. Teyssier, Z. Nagy, J. Montillaud, C. Vastel, J. Cernicharo, M. Roellig, V. Ossenkopf-Okada, E. A. Bergin

Abstract: In bright photodissociation regions (PDRs) associated to massive star formation, the presence of dense "clumps" that are immersed in a less dense interclump medium is often proposed to explain the difficulty of models to account for the observed gas emission in high-excitation lines. We aim at presenting a comprehensive view of the modeling of the CO rotational ladder in PDRs, including the high-J… ▽ More In bright photodissociation regions (PDRs) associated to massive star formation, the presence of dense "clumps" that are immersed in a less dense interclump medium is often proposed to explain the difficulty of models to account for the observed gas emission in high-excitation lines. We aim at presenting a comprehensive view of the modeling of the CO rotational ladder in PDRs, including the high-J lines that trace warm molecular gas at PDR interfaces. We observed the 12CO and 13CO ladders in two prototypical PDRs, the Orion Bar and NGC 7023 NW using the instruments onboard Herschel. We also considered line emission from key species in the gas cooling of PDRs (C+, O, H2) and other tracers of PDR edges such as OH and CH+. All the intensities are collected from Herschel observations, the literature and the Spitzer archive and are analyzed using the Meudon PDR code. A grid of models was run to explore the parameter space of only two parameters: thermal gas pressure and a global scaling factor that corrects for approximations in the assumed geometry. We conclude that the emission in the high-J CO lines, which were observed up to Jup=23 in the Orion Bar (Jup=19 in NGC7023), can only originate from small structures of typical thickness of a few 1e-3 pc and at high thermal pressures (Pth~1e8 K cm-3). Compiling data from the literature, we found that the gas thermal pressure increases with the intensity of the UV radiation field given by G0, following a trend in line with recent simulations of the photoevaporation of illuminated edges of molecular clouds. This relation can help rationalising the analysis of high-J CO emission in massive star formation and provides an observational constraint for models that study stellar feedback on molecular clouds. △ Less

Submitted 15 April, 2018; v1 submitted 11 January, 2018; originally announced January 2018.

Comments: Accepted for publication in A&A; 20 pages, 11 figures

Journal ref: A&A 615, A129 (2018)

Galactic cold cores IX. Column density structures and radiative transfer modelling

Authors: M. Juvela, J. Malinen, J. Montillaud, V. -M. Pelkonen, I. Ristorcelli, L. V. Toth

Abstract: The Galactic Cold Cores (GCC) project has made Herschel observations of interstellar clouds where Planck detected compact sources of cold dust emission. Our aim is to characterise the structure of the clumps and their parent clouds. We also examine the accuracy to which the structure of dense clumps can be determined from submillimetre data. We use standard statistical methods to characterise the… ▽ More The Galactic Cold Cores (GCC) project has made Herschel observations of interstellar clouds where Planck detected compact sources of cold dust emission. Our aim is to characterise the structure of the clumps and their parent clouds. We also examine the accuracy to which the structure of dense clumps can be determined from submillimetre data. We use standard statistical methods to characterise the GCC fields. Clumps are extracted using column density thresholding and we construct for each field a three-dimensional radiative transfer (RT) model. These are used to estimate the relative radiation field intensities, clump stability, and the uncertainty of column density estimates. We examine the radial column density profiles of the clumps. In the GCC fields, the structure noise follows the relations previously established at larger scales. The fractal dimension has no significant dependence on column density and the values D = 1.25 +- 0.07 are only slightly lower than in typical molecular clouds. The column density PDFs exhibit large variations, e.g. in the case of externally compressed clouds. At scales r>0.1 pc, the radial column density distributions of the clouds follow an average relation of N~r^{-1}. In spite of a great variety of clump morphology, clumps tend to follow a similar N~r^{-1} relation below r~0.1 pc. RT calculations indicate only factor of 2.5 variation in the local radiation field intensity. The fraction of gravitationally bound clumps increases significantly in regions with A_V > 5 mag but most bound objects appear to be pressure-confined. The GCC host clouds have statistical properties similar to general molecular clouds. The gravitational stability, peak column density, and clump orientation are connected to the cloud background while most other statistics (e.g. D and radial profiles) are insensitive to the environment. △ Less

Submitted 15 January, 2018; v1 submitted 8 January, 2018; originally announced January 2018.

Comments: accepted to A&A, 21 pages (37 pages with appendices, this version including language corrections)

Journal ref: A&A 614, A83 (2018)

A CO survey on a sample of Herschel cold clumps

Authors: O. Fehér, M. Juvela, T. Lunttila, J. Montillaud, I. Ristorcelli, S. Zahorecz, L. V. Tóth

Abstract: The physical state of cold cloud clumps has a great impact on the process and efficiency of star formation and the masses of the forming stars inside these objects. The sub-millimetre survey of the Planck space observatory and the far-infrared follow-up mapping of the Herschel space telescope provide an unbiased, large sample of these cold objects. We have observed $^{12}$CO(1$-$0) and $^{13}$CO(1… ▽ More The physical state of cold cloud clumps has a great impact on the process and efficiency of star formation and the masses of the forming stars inside these objects. The sub-millimetre survey of the Planck space observatory and the far-infrared follow-up mapping of the Herschel space telescope provide an unbiased, large sample of these cold objects. We have observed $^{12}$CO(1$-$0) and $^{13}$CO(1$-$0) emission in 35 clumps in 26 Herschel fields sampling different environments in the Galaxy. Densities and temperatures were calculated from both the dust continuum and the molecular line data, kinematic distances were derived using $^{13}$CO line velocities and clump sizes and masses were calculated by fitting 2D Gaussian functions to the optical depth distribution maps. Clump masses and virial masses were estimated assuming an upper and lower limit on the kinetic temperatures and considering uncertainties due to distance limitations. The excitation temperatures are between 8.5$-$19.5 K, while the Herschel-derived dust colour temperatures are 12$-$16 K. The sizes (0.1$-$3 pc), $^{13}$CO column densities (0.5$-$44$\times$10$^{15}$ cm$^{-2}$) and masses (from less than 0.1 $M_{\odot}$ to more than 1500 $M_{\odot}$) of the objects span broad ranges. Eleven gravitationally unbound clumps were found, many of them smaller than 0.3 pc, but large, parsec-scale clouds with a few hundred solar masses appear as well. Colder clumps have generally high column densities but warmer objects appear at both low and higher column densities. The clump column densities derived from the line and dust observations correlate well, but are heavily affected by uncertainties of the dust properties, varying molecular abundances and optical depth effects. △ Less

Submitted 31 July, 2017; originally announced August 2017.

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

Journal ref: A&A 606, A102 (2017)

Correlation of gas dynamics and dust in the evolved filament G82.65-02.00

Authors: M. Saajasto, M. Juvela, K. Dobashi, T. Shimoikura, I. Ristorcelli, J. Montillaud, D. J. Marshall, J. Malinen, V. -M. Pelkonen, O. Fehér, A. Rivera-Ingraham, L. V. Toth, L. Montier, J-Ph. Bernard, T. Onishi

Abstract: The combination of line and continuum observations can provide vital insights to the formation and fragmentation of filaments and the initial conditions for star formation. We have carried out line observations to map the kinematics of an evolved, actively star forming filament G82.65-2.00. The filament was first The combination of line and continuum observations can provide vital insights to the… ▽ More The combination of line and continuum observations can provide vital insights to the formation and fragmentation of filaments and the initial conditions for star formation. We have carried out line observations to map the kinematics of an evolved, actively star forming filament G82.65-2.00. The filament was first The combination of line and continuum observations can provide vital insights to the formation and fragmentation of filaments and the initial conditions for star formation. We have carried out line observations to map the kinematics of an evolved, actively star forming filament G82.65-2.00. The filament was first identified from the Planck data as a region of particularly cold dust emission and was mapped at 100-500 $μ$m as a part of Herschel key program Galactic Cold Cores. The observations reveal several velocity components in the field, with strongest line emission concentrated to velocity range $\sim [3,5]$ km $\rm s^{-1}$. The column density of molecular hydrogen along the filament varies from 1.0 to 2.3 $\times 10^{22}$ $\rm cm^{2}$. We have examined six cold clumps from the central part of the filament. The clumps have masses in the range $10 - 20$ $M_{\odot}$ ($\sim70M_{\odot}$ in total) and are close to or above the virial mass. Furthermore, the main filament is heavily fragmented and most of the the substructures have a mass lower than or close to the virial mass, suggesting that the filament is dispersing as a whole. Position-velocity maps of $\rm ^{12}CO$ and $\rm ^{13}CO$ lines indicate that at least one of the striations is kinematically connected to two of the clumps, potentially indicating mass accretion from the striation onto the main filament. We tentatively estimate the accretion rate to be $\dot{M}$ = $2.23 \times 10^{-6}$ $ M_{\odot} / \rm year$. △ Less

Submitted 3 July, 2017; originally announced July 2017.

Comments: Accepted for publication in A&A

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

Galactic Cold Cores. VIII. Filament formation and evolution: Filament properties in context with evolutionary models

Authors: A. Rivera-Ingraham, I. Ristorcelli, M. Juvela, J. Montillaud, A. Men'shchikov, J. Malinen, V. -M. Pelkonen, A. Marston, P. G. Martin, L. Pagani, R. Paladini, D. Paradis, N. Ysard, D. Ward-Thompson, J. -P. Bernard, D. J. Marshall, L. Montier, L. V. Toth

Abstract: Filaments are key for star formation models. As part of the study carried out by the Herschel GCC Programme, here we study the filament properties presented in GCC.VII in context with theoretical models of filament formation and evolution. A conservative sample of filaments at a distance D<500pc was extracted with the Getfilaments algorithm. Their physical structure was quantified according to two… ▽ More Filaments are key for star formation models. As part of the study carried out by the Herschel GCC Programme, here we study the filament properties presented in GCC.VII in context with theoretical models of filament formation and evolution. A conservative sample of filaments at a distance D<500pc was extracted with the Getfilaments algorithm. Their physical structure was quantified according to two main components: the central (Gaussian) region (core component), and the power-law like region dominating the filament column density profile at larger radii (wing component). The properties and behaviour of these components relative to the total linear mass density of the filament and its environmental column density were compared with theoretical models describing the evolution of filaments under gravity-dominated conditions. The feasibility of a transition to supercritical state by accretion is dependent on the combined effect of filament intrinsic properties and environmental conditions. Reasonably self-gravitating (high Mline-core) filaments in dense environments (av\sim3mag) can become supercritical in timescales of t\sim1Myr by accreting mass at constant or decreasing width. The trend of increasing Mline-tot (Mline-core and Mline-wing), and ridge Av with background also indicates that the precursors of star-forming filaments evolve coevally with their environment. The simultaneous increase of environment and filament Av explains the association between dense environments and high Mline-core values, and argues against filaments remaining in constant single-pressure equilibrium states. The simultaneous growth of filament and background in locations with efficient mass assembly, predicted in numerical models of collapsing clouds, presents a suitable scenario for the fulfillment of the combined filament mass-environment criterium that is in quantitative agreement with Herschel observations. △ Less

Submitted 24 February, 2017; originally announced February 2017.

Comments: Accepted in A&A

Journal ref: A&A 601, A94 (2017)

Galactic Cold Cores VII: Filament Formation and Evolution - Methods & Observational Constraints

Authors: A. Rivera-Ingraham, I. Ristorcelli, M. Juvela, J. Montillaud, A. Men'shchikov, J. Malinen, V. -M. Pelkonen, A. Marston, P. G. Martin, L. Pagani, R. Paladini, D. Paradis, N. Ysard, D. Ward-Thompson, J. -P. Bernard, D. J. Marshall, L. Montier, V. Tóth

Abstract: The association of filaments with protostellar objects has made these structures a priority target in star formation studies. The datasets of the Herschel Galactic Cold Cores Key Programme allow for a statistical study of filaments with a wide range of intrinsic and environmental characteristics. Characterisation of this sample can be used to identify key physical parameters and quantify the role… ▽ More The association of filaments with protostellar objects has made these structures a priority target in star formation studies. The datasets of the Herschel Galactic Cold Cores Key Programme allow for a statistical study of filaments with a wide range of intrinsic and environmental characteristics. Characterisation of this sample can be used to identify key physical parameters and quantify the role of environment in the formation of supercritical filaments. Filaments were extracted from fields at D<500pc with the getfilaments algorithm and characterised according to their column density profiles and intrinsic properties. Each profile was fitted with a beam-convolved Plummer-like function and quantified based on the relative contributions from the filament 'core', represented by a Gaussian, and 'wing' component, dominated by the power-law of the Plummer-like function. These parameters were examined for populations associated with different background levels. We find that filaments increase their core (Mcore) and wing (Mwing) contributions while increasing their total linear mass density (Mtot). Both components appear to be linked to the local environment, with filaments in higher backgrounds having systematically more massive Mcore and Mwing. This dependence on the environment supports an accretion-based model for filament evolution in the local neighbourhood (D<500pc). Structures located in the highest backgrounds develop the highest central Av, Mcore, and Mwing as Mtot increases with time, favoured by the local availability of material and the enhanced gravitational potential. Our results indicate that filaments acquiring a significantly massive central region with Mcore>Mcrit/2 may become supercritical and form stars. This translates into a need for filaments to become at least moderately self-gravitating in order to undergo localised star formation or become star-forming filaments. △ Less

Submitted 30 March, 2016; originally announced March 2016.

Comments: Accepted for publication in A&A

Journal ref: A&A 591, A90 (2016)

Matching dust emission structures and magnetic field in high-latitude cloud L1642: comparing Herschel and Planck maps

Authors: J. Malinen, L. Montier, J. Montillaud, M. Juvela, I. Ristorcelli, S. E. Clark, O. Berné, J. -Ph. Bernard, V. -M. Pelkonen, D. C. Collins

Abstract: The nearby cloud L1642 is one of only two known very high latitude (|b| > 30 deg) clouds actively forming stars. It is a rare example of star formation in isolated conditions, and can reveal important details of star formation in general, e.g., of the effect of magnetic fields. We compare Herschel dust emission structures and magnetic field orientation revealed by Planck polarization maps in L1642… ▽ More The nearby cloud L1642 is one of only two known very high latitude (|b| > 30 deg) clouds actively forming stars. It is a rare example of star formation in isolated conditions, and can reveal important details of star formation in general, e.g., of the effect of magnetic fields. We compare Herschel dust emission structures and magnetic field orientation revealed by Planck polarization maps in L1642. The high-resolution ($\sim20"$) Herschel data reveal a complex structure including a dense, compressed central clump, and low density striations. The Planck polarization data (at 10$'$ resolution) reveal an ordered magnetic field pervading the cloud and aligned with the surrounding striations. There is a complex interplay between the cloud structure and large scale magnetic field. This suggests that the magnetic field is closely linked to the formation and evolution of the cloud. CO rotational emission confirms that the striations are connected with the main clumps and likely to contain material either falling into or flowing out of the clumps. There is a clear transition from aligned to perpendicular structures approximately at a column density of $N_{\rm{H}} = 1.6 \times 10^{21}\, {\rm cm}^{-2}$. Comparing the Herschel maps with the Planck polarization maps shows the close connection between the magnetic field and cloud structure even in the finest details of the cloud. △ Less

Submitted 3 May, 2016; v1 submitted 11 December, 2015; originally announced December 2015.

Comments: accepted to MNRAS; 14 pages; changes: moderate text and figure editing

Near-infrared extinction with discretised stellar colours

Authors: M. Juvela, J. Montillaud

Abstract: Several methods exist to convert near-infrared (NIR) stellar observations into extinction maps. We present a new method based on NIR multiband observations. The method uses a discretised version of the distribution of intrinsic stellar colours. A number of variations of the basic method are tested, and the results are compared to NICER calculations. When photometric errors are large, the results a… ▽ More Several methods exist to convert near-infrared (NIR) stellar observations into extinction maps. We present a new method based on NIR multiband observations. The method uses a discretised version of the distribution of intrinsic stellar colours. A number of variations of the basic method are tested, and the results are compared to NICER calculations. When photometric errors are large, the results are close to those of NICER method but some advantages can be seen when the distribution of intrinsic colours cannot be described well with a single covariance matrix. A priori information about relative column density variations at sub-beam scales can result in a significant increase in accuracy. The results may be further improved by considering the magnitude dependence of the intrinsic colours. Thus, the new methods are useful mostly when photometric errors are small, the distribution of intrinsic colours is well known, or one has prior knowledge of the small-scale structures. △ Less

Submitted 3 November, 2015; originally announced November 2015.

Comments: Accepted to A&A

Allsky NICER and NICEST extinction maps based on the 2MASS near-infrared survey

Authors: M. Juvela, J. Montillaud

Abstract: Extinction remains one of the most reliable methods of measuring column density of nearby Galactic interstellar clouds. The current and ongoing near-infrared surveys enable the mapping of extinction over large sky areas. We produce allsky extinction maps using the 2MASS near-infrared survey. We use the NICER and NICEST methods to convert the near-infrared colour excesses to extinction estimates. T… ▽ More Extinction remains one of the most reliable methods of measuring column density of nearby Galactic interstellar clouds. The current and ongoing near-infrared surveys enable the mapping of extinction over large sky areas. We produce allsky extinction maps using the 2MASS near-infrared survey. We use the NICER and NICEST methods to convert the near-infrared colour excesses to extinction estimates. The results are presented in Healpix format at the resolutions of 3.0, 4.5, and 12.0 arcmin. The main results of this study are the calculated J-band extinction maps. The comparison with earlier large-scale extinction mappings shows good correspondence but also demonstrates the presence of resolution-dependent bias. A large fraction of the bias can be corrected by using the NICEST method. For individual regions, best extinction estimates are obtained by careful analysis of the local stellar population and the use of the highest resolution afforded by the stellar density. However, the uniform allsky maps should still be useful for many global studies and as the first step into the investigation of individual clouds. △ Less

Submitted 3 November, 2015; originally announced November 2015.

Comments: Accepted to A&A

Journal ref: A&A 585, A38 (2016)

30 years of cosmic fullerenes

Authors: O. Berne, J. Montillaud, G. Mulas, C. Joblin

Abstract: In 1985, "During experiments aimed at understanding the mechanisms by which long-chain carbon molecules are formed in interstellar space and circumstellar shells", Harry Kroto and his collaborators serendipitously discovered a new form of carbon: fullerenes. The most emblematic fullerene (i.e. C$_{60}$ "buckminsterfullerene"), contains exactly 60 carbon atoms organized in a cage-like structure sim… ▽ More In 1985, "During experiments aimed at understanding the mechanisms by which long-chain carbon molecules are formed in interstellar space and circumstellar shells", Harry Kroto and his collaborators serendipitously discovered a new form of carbon: fullerenes. The most emblematic fullerene (i.e. C$_{60}$ "buckminsterfullerene"), contains exactly 60 carbon atoms organized in a cage-like structure similar to a soccer ball. Since their discovery impacted the field of nanotechnologies, Kroto and colleagues received the Nobel prize in 1996. The cage-like structure, common to all fullerene molecules, gives them unique properties, in particular an extraordinary stability. For this reason and since they were discovered in experiments aimed to reproduce conditions in space, fullerenes were sought after by astronomers for over two decades, and it is only recently that they have been firmly identified by spectroscopy, in evolved stars and in the interstellar medium. This identification offers the opportunity to study the molecular physics of fullerenes in the unique physical conditions provided by space, and to make the link with other large carbonaceous molecules thought to be present in space : polycyclic aromatic hydrocarbons. △ Less

Submitted 27 October, 2015; v1 submitted 6 October, 2015; originally announced October 2015.

Comments: To appear in the Proceedings of the annual meeting of the French society of astronomy and astrophysics (SF2A 2015)

Galactic cold cores VI. Dust opacity spectral index

Authors: M. Juvela, K. Demyk, Y. Doi, A. Hughes, C. Lefevre, D. J. Marshall, C. Meny, J. Montillaud, L. Pagani, D. Paradis, I. Ristorcelli, J. Malinen, L. A. Montier, R. Paladini, V. -M. Pelkonen, A. Rivera-Ingraham

Abstract: The Galactic Cold Cores project has made Herschel observations of 116 fields where the Planck survey has found signs of cold dust emission. The fields contain sources in different environments and different phases of star formation. The dust opacity spectral index beta and the dust colour temperature T are derived using Herschel and Planck data. The relation between beta and T is examined for the… ▽ More The Galactic Cold Cores project has made Herschel observations of 116 fields where the Planck survey has found signs of cold dust emission. The fields contain sources in different environments and different phases of star formation. The dust opacity spectral index beta and the dust colour temperature T are derived using Herschel and Planck data. The relation between beta and T is examined for the whole sample and inside individual fields. Based on IRAS and Planck data, the fields are characterised by a median colour temperature of 16.1 K and a median opacity spectral index of beta=1.84. We observe a clear T-beta anti-correlation. In Herschel observations, constrained at lower resolution by Planck data, the variations follow the column density structure and beta(FIR) can rise to ~2.2 in individual clumps. The Planck 217 GHz band shows a systematic excess that is consistent with a general flattening of the dust emission spectrum at millimetre wavelengths. When fitted separately below and above 700 um, the median spectral index values are beta(FIR) ~ 1.91 and beta(mm) ~ 1.66. The spectral index changes as a function of column density and wavelength. Beta variations are partly masked by temperature gradients and the changes in the intrinsic grain properties may be even greater. △ Less

Submitted 26 September, 2015; originally announced September 2015.

Comments: Accepted to A&A, 138 pages

Journal ref: A&A 584, A94 (2015)

Top-down formation of fullerenes in the interstellar medium

Authors: O. Berne, J. Montillaud, C. Joblin

Abstract: [Abridged] Fullerenes have been recently detected in various circumstellar and interstellar environments, raising the question of their formation pathway. It has been proposed that they can form by the photo-chemical processing of large polycyclic aromatic hydrocarbons (PAHs). Following our previous work on the evolution of PAHs in the NGC 7023 reflection nebula, we evaluate, using photochemical m… ▽ More [Abridged] Fullerenes have been recently detected in various circumstellar and interstellar environments, raising the question of their formation pathway. It has been proposed that they can form by the photo-chemical processing of large polycyclic aromatic hydrocarbons (PAHs). Following our previous work on the evolution of PAHs in the NGC 7023 reflection nebula, we evaluate, using photochemical modeling, the possibility that the PAH C$_{66}$H$_{20}$ (i.e. circumovalene) can lead to the formation of C$_{60}$ upon irradiation by ultraviolet photons. The chemical pathway involves full dehydrogenation, folding into a floppy closed cage and shrinking of the cage by loss of C$_2$ units until it reaches the symmetric C$_{60}$ molecule. At 10" from the illuminating star and with realistic molecular parameters, the model predicts that 100% of C$_{66}$H$_{20}$ is converted into C$_{60}$ in $\sim$ 10$^5$ years, a timescale comparable to the age of the nebula. Shrinking appears to be the kinetically limiting step of the whole process. Hence, PAHs larger than C$_{66}$H$_{20}$ are unlikely to contribute significantly to the formation of C$_{60}$, while PAHs containing between 60 and 66 C atoms should contribute to the formation of C$_{60}$ with shorter timescales, and PAHs containing less than 60 C atoms will be destroyed. Assuming a classical size distribution for the PAH precursors, our model predicts absolute abundances of C$_{60}$ are up to several $10^{-4}$ of the elemental carbon, i.e. less than a percent of the typical interstellar PAH abundance, which is consistent with observational studies. According to our model, once formed, C$_{60}$ can survive much longer than other fullerenes because of the remarkable stability of the C$_{60}$ molecule at high internal energies.Hence, a natural consequence is that C$_{60}$ is more abundant than other fullerenes in highly irradiated environments. △ Less

Submitted 25 January, 2016; v1 submitted 12 March, 2015; originally announced March 2015.

Comments: Accepted for publication in A&A. Latest version contains the corrected version of Fig. 2

Journal ref: A&A 577, A133 (2015)

Galactic cold cores V. Dust opacity

Authors: M. Juvela, I. Ristorcelli, D. J. Marshall, J. Montillaud, V. -M. Pelkonen, N. Ysard, P. McGehee, R. Paladini, L. Pagani, J. Malinen, A., Rivera-Ingraham, C. Lefevre, L. V. Toth, L. A. Montier, J. -P. Bernard, P. Martin

Abstract: The project Galactic Cold Cores has made Herschel observations of interstellar clouds where the Planck satellite survey has located cold and compact clumps. The sources range from starless clumps to protostellar cores. We examine 116 Herschel fields to estimate the submillimetre dust opacity and its variations. The submillimetre dust opacity was derived from Herschel data, and near-infrared observ… ▽ More The project Galactic Cold Cores has made Herschel observations of interstellar clouds where the Planck satellite survey has located cold and compact clumps. The sources range from starless clumps to protostellar cores. We examine 116 Herschel fields to estimate the submillimetre dust opacity and its variations. The submillimetre dust opacity was derived from Herschel data, and near-infrared observations of the reddening of background stars are converted into near-infrared optical depth. We studied the systematic errors affecting these parameters and used modelling to correct for the expected biases. The ratio of 250um and J band opacities is correlated with the cloud location and star formation activity. We find a median ratio of tau(250um)/tau(J)= (1.6+-0.2)*10^-3, which is more than three times the mean value in diffuse medium. Assuming a spectral index beta=1.8 instead of beta=2.0, the value would be lower by ~30%. No significant systematic variation is detected with Galactocentric distance or with Galactic height. The tau(250um)/tau(J) maps reveal six fields with clear increase of submillimetre opacity of up to tau(250um)/tau(J) ~ 4*10^-3. These are all nearby fields with spatially resolved clumps of high column density. We interpret the increase in the far-infrared opacity as a sign of grain growth in the densest and coldest regions of interstellar clouds. △ Less

Submitted 28 January, 2015; originally announced January 2015.

Comments: 28 pages, A&A, in press

Absolute evaporation rates of non-rotating neutral PAH clusters

Authors: Julien Montillaud, Christine Joblin

Abstract: Clusters of polycyclic aromatic hydrocarbons (PAHs) have been proposed as candidates for evaporating very small grains, which are thought to be precursors of free-flying PAHs. Evaporation rates have been calculated so far only for species containing up to a few 100 C atoms, whereas interstellar PAH clusters could contain up to ~1000 C atoms. We present a method that generalises the calculation of… ▽ More Clusters of polycyclic aromatic hydrocarbons (PAHs) have been proposed as candidates for evaporating very small grains, which are thought to be precursors of free-flying PAHs. Evaporation rates have been calculated so far only for species containing up to a few 100 C atoms, whereas interstellar PAH clusters could contain up to ~1000 C atoms. We present a method that generalises the calculation of the statistical evaporation rate of large PAH clusters and provides rates for species containing up to ~1000 C-atoms. The evaporation of non-rotating neutral homo-molecular PAH clusters containing up to 12 molecules from a family of highly symmetric compact PAHs is studied. Statistical calculations were performed and completed with molecular dynamics simulations at high internal energies to provide absolute values for the evaporation rate and distributions of kinetic energy released. The calculations used explicit atom-atom Lennard-Jones potentials in the rigid molecule approximation. A new method is proposed to take both inter- and intra-molecular vibrations into account. Without any parameter adjustment, the calculated evaporation rates agree well with available experimental data. We find that the non-rotation assumption has a limited impact on the evaporation rates. The photostability of PAH clusters increases dramatically with the size of molecules in the clusters, and to a lesser extent with the number of molecules in the clusters. For values of the UV radiation field that are typical of the regions where evaporating very small grains are observed, the smallest clusters in this study (~50 C-atoms) are found to be quickly photo-evaporated, whereas the largest clusters (~1000 C-atoms) are photostable. Our results support the idea that large PAH clusters are good candidates for evaporating very small grains. △ Less

Submitted 8 April, 2014; originally announced April 2014.

Comments: 13 pages, 10 figures

Journal ref: A&A 567, A45 (2014)

Multiwavelength study of the high-latitude cloud L1642: chain of star formation

Authors: J. Malinen, M. Juvela, S. Zahorecz, A. Rivera-Ingraham, J. Montillaud, K. Arimatsu, J. -Ph. Bernard, Y. Doi, L. Haikala, R. Kawabe, G. Marton, P. McGehee, V. -M. Pelkonen, I. Ristorcelli, Y. Shimajiri, S. Takita, L. V. Toth, T. Tsukagoshi, N. Ysard

Abstract: L1642 is one of the two high galactic latitude (|b| > 30deg) clouds confirmed to have active star formation. We examine the properties of this cloud, especially the large-scale structure, dust properties, and compact sources in different stages of star formation. We present high-resolution far-infrared and submm observations with the Herschel and AKARI satellites and mm observations with the AzTEC… ▽ More L1642 is one of the two high galactic latitude (|b| > 30deg) clouds confirmed to have active star formation. We examine the properties of this cloud, especially the large-scale structure, dust properties, and compact sources in different stages of star formation. We present high-resolution far-infrared and submm observations with the Herschel and AKARI satellites and mm observations with the AzTEC/ASTE telescope, which we combined with archive data from near- and mid-infrared (2MASS, WISE) to mm observations (Planck). The Herschel observations, combined with other data, show a sequence of objects from a cold clump to young stellar objects at different evolutionary stages. Source B-3 (2MASS J04351455-1414468) appears to be a YSO forming inside the L1642 cloud, instead of a foreground brown dwarf, as previously classified. Herschel data reveal striation in the diffuse dust emission around L1642. The western region shows striation towards NE and has a steeper column density gradient on its southern side. The densest central region has a bow-shock like structure showing compression from the west and a filamentary tail extending towards east. The differences suggest that these may be spatially distinct structures, aligned only in projection. We derive values of the dust emission cross-section per H nucleon for different regions of the cloud. Modified black-body fits to the spectral energy distribution of Herschel and Planck data give emissivity spectral index beta values 1.8-2.0 for the different regions. The compact sources have lower beta values and show an anticorrelation between T and beta. Markov chain Monte Carlo calculations demonstrate the strong anticorrelation between beta and T errors and the importance of mm Planck data in constraining the estimates. L1642 reveals a more complex structure and sequence of star formation than previously known. △ Less

Submitted 9 April, 2014; v1 submitted 11 February, 2014; originally announced February 2014.

Comments: 22 pages, 18 figures, accepted to Astronomy & Astrophysics; abstract shortened and figures reduced for astroph

Estimation of high-resolution dust column density maps: Empirical model fits

Authors: M. Juvela, J. Montillaud

Abstract: Sub-millimetre dust emission is an important tracer of density N of dense interstellar clouds. One has to combine surface brightness information at different spatial resolutions, and specific methods are needed to derive N at a resolution higher than the lowest resolution of the observations. Some methods have been discussed in the literature, including a method (in the following, method B) that c… ▽ More Sub-millimetre dust emission is an important tracer of density N of dense interstellar clouds. One has to combine surface brightness information at different spatial resolutions, and specific methods are needed to derive N at a resolution higher than the lowest resolution of the observations. Some methods have been discussed in the literature, including a method (in the following, method B) that constructs the N estimate in stages, where the smallest spatial scales being derived only use the shortest wavelength maps. We propose simple model fitting as a flexible way to estimate high-resolution column density maps. Our goal is to evaluate the accuracy of this procedure and to determine whether it is a viable alternative for making these maps. The new method consists of model maps of column density (or intensity at a reference wavelength) and colour temperature. The model is fitted using Markov chain Monte Carlo (MCMC) methods, comparing model predictions with observations at their native resolution. We analyse simulated surface brightness maps and compare its accuracy with method B and the results that would be obtained using high-resolution observations without noise. The new method is able to produce reliable column density estimates at a resolution significantly higher than the lowest resolution of the input maps. Compared to method B, it is relatively resilient against the effects of noise. The method is computationally more demanding, but is feasible even in the analysis of large Herschel maps. The proposed empirical modelling method E is demonstrated to be a good alternative for calculating high-resolution column density maps, even with considerable super-resolution. Both methods E and B include the potential for further improvements, e.g., in the form of better a priori constraints. △ Less

Submitted 27 July, 2013; originally announced July 2013.

Comments: Accepted to A&A

The degeneracy between dust colour temperature and spectral index. Comparison of methods for estimating the beta(T) relation

Authors: M. Juvela, J. Montillaud, N. Ysard, T. Lunttila

Abstract: Sub-millimetre dust emission provides information on the physics of interstellar clouds and dust. Noise can produce anticorrelation between the colour temperature T_C and the spectral index beta. This must be separated from the intrinsic beta(T) relation of dust. We compare methods for the analysis of the beta(T) relation. We examine sub-millimetre observations simulated as simple modified black b… ▽ More Sub-millimetre dust emission provides information on the physics of interstellar clouds and dust. Noise can produce anticorrelation between the colour temperature T_C and the spectral index beta. This must be separated from the intrinsic beta(T) relation of dust. We compare methods for the analysis of the beta(T) relation. We examine sub-millimetre observations simulated as simple modified black body emission or using 3D radiative transfer modelling. In addition to chi^2 fitting, we examine the results of the SIMEX method, basic Bayesian model, hierarchical models, and one method that explicitly assumes a functional form for beta(T). All methods exhibit some bias. Bayesian method shows significantly lower bias than direct chi^2 fits. The same is true for hierarchical models that also result in a smaller scatter in the temperature and spectral index values. However, significant bias was observed in cases with high noise levels. Beta and T estimates of the hierarchical model are biased towards the relation determined by the data with the highest S/N ratio. This can alter the recovered beta(T) function. With the method where we explicitly assume a functional form for the beta(T) relation, the bias is similar to the Bayesian method. In the case of an actual Herschel field, all methods agree, showing some degree of anticorrelation between T and beta. The Bayesian method and the hierarchical models can both reduce the noise-induced parameter correlations. However, all methods can exhibit non-negligible bias. This is particularly true for hierarchical models and observations of varying signal-to-noise ratios and must be taken into account when interpreting the results. △ Less

Submitted 9 May, 2013; originally announced May 2013.

Comments: Submitted to A&A, 18 pages

Spatial distribution of small hydrocarbons in the neighborhood of the Ultra Compact HII region Monoceros R2

Authors: P. Pilleri, S. Treviño-Morales, A. Fuente, C. Joblin, J. Cernicharo, M. Gerin, S. Viti, O. Berné, J. R. Goicoechea, J. Pety, M. Gonzalez-García, J. Montillaud, V. Ossenkopf, C. Kramer, S. García-Burillo, F. Le Petit, J. Le Bourlot

Abstract: We study the chemistry of small hydrocarbons in the photon-dominated regions (PDRs) associated with the ultra-compact HII region Mon R2. Our goal is to determine the variations of the abundance of small hydrocarbons in a high-UV irradiated PDR and investigate their chemistry. We present an observational study of CH, CCH and c-C$_3$H$_2$ in Mon R2 combining data obtained with the IRAM 30m telescope… ▽ More We study the chemistry of small hydrocarbons in the photon-dominated regions (PDRs) associated with the ultra-compact HII region Mon R2. Our goal is to determine the variations of the abundance of small hydrocarbons in a high-UV irradiated PDR and investigate their chemistry. We present an observational study of CH, CCH and c-C$_3$H$_2$ in Mon R2 combining data obtained with the IRAM 30m telescope and Herschel. We determine the column densities of these species, and compare their spatial distributions with that of polycyclic aromatic hydrocarbon (PAH). We compare the observational results with different chemical models to explore the relative importance of gas-phase, grain-surface and time-dependent chemistry in these environments. The emission of the small hydrocarbons show different patterns. The CCH emission is extended while CH and c-C$_3$H$_2$ are concentrated towards the more illuminated layers of the PDR. The ratio of the column densities of c-C$_3$H$_2$ and CCH shows spatial variations up to a factor of a few, increasing from $N(c-C$_3$H$_2$)/N(CCH)\approx0.004$ in the envelope to a maximum of $\sim0.015-0.029$ towards the 8$μ$m emission peak. Comparing these results with other galactic PDRs, we find that the abundance of CCH is quite constant over a wide range of G$_0$, whereas the abundance of c-C$_3$H$_2$ is higher in low-UV PDRs. In Mon R2, the gas-phase steady-state chemistry can account relatively well for the abundances of CH and CCH in the most exposed layers of the PDR, but falls short by a factor of 10 to reproduce c-C$_3$H$_2$. In the molecular envelope, time-dependent effects and grain surface chemistry play a dominant role in determining the hydrocarbons abundances. Our study shows that CCH and c-C$_3$H$_2$ present a complex chemistry in which UV photons, grain-surface chemistry and time dependent effects contribute to determine their abundances. △ Less

Submitted 20 March, 2013; originally announced March 2013.

Comments: 18 pages, 11 figures, 7 tables. Proposed for acceptance in A&A. Abstract abridged

Evolution of PAHs in photodissociation regions: Hydrogenation and charge states

Authors: J. Montillaud, C. Joblin, D. Toublanc

Abstract: Various studies have emphasized variations of the charge state and composition of the interstellar polycyclic aromatic hydrocarbon (PAH) population in photodissociation regions (PDRs). We aim to model the spatial evolution of the charge and hydrogenation states of PAHs in PDRs. We focus on the specific case of the north-west (NW) PDR of NGC 7023 and also discuss the case of the diffuse interstella… ▽ More Various studies have emphasized variations of the charge state and composition of the interstellar polycyclic aromatic hydrocarbon (PAH) population in photodissociation regions (PDRs). We aim to model the spatial evolution of the charge and hydrogenation states of PAHs in PDRs. We focus on the specific case of the north-west (NW) PDR of NGC 7023 and also discuss the case of the diffuse interstellar medium (ISM). The physical conditions in NGC 7023 NW are modelled using a state-of-the-art PDR code. We then use a new PAH chemical evolution model that includes recent experimental data on PAHs and describes multiphoton events. We consider a family of compact PAHs bearing up to 96 carbon atoms. The calculated ionization ratio is in good agreement with observations in NGC 7023 NW. Within the PDR, PAHs evolve into three major populations: medium-sized PAHs (50<Nc<90) are normally hydrogenated, larger PAHs (Nc>90) can be superhydrogenated, and smaller species (Nc<50) are fully dehydrogenated. In the cavity, where the fullerene C60 was recently detected, all the studied PAHs are found to be quickly fully dehydrogenated. PAH chemical evolution exhibits a complex non-linear behaviour as a function of the UV radiation field because of multiphoton events. Steady state for hydrogenation is reached on timescales ranging from less than a year for small PAHs, up to 10000 years for large PAHs at Av=1. We identified critical reactions that need more studies. Our new model allows us to rationalize the observational constraints without any fitting parameter. PAHs smaller than 50 carbon atoms are not expected to survive in the NGC 7023 NW PDR. A similar conclusion is obtained for the diffuse ISM. Carbon clusters turn out to be end products of PAH photodissociation, and the evolution of these clusters needs to be investigated further to evaluate their impact on the chemical and physical evolution of PDRs. △ Less

Submitted 28 January, 2013; originally announced January 2013.

Comments: 16 pages, 10 figures; Accepted for publication in A&A

Herschel / HIFI observations of CO, H2O and NH3 in Mon R2

Authors: P. Pilleri, A. Fuente, J. Cernicharo, V. Ossenkopf, O. Berné, M. Gerin, J. Pety, J. R. Goicoechea, J. R. Rizzo, J. Montillaud, M. González-García, C. Joblin, J. Le Bourlot, F. Le Petit, C. Kramer

Abstract: Context. Mon R2 is the only ultracompact HII region (UCHII) where the associated photon-dominated region (PDR) can be resolved with Herschel. Due to its brightness and proximity, it is the best source to investigate the chemistry and physics of highly UV-irradiated PDRs. Aims. Our goal is to estimate the abundance of H2O and NH3 in this region and investigate their origin. Methods. We present new… ▽ More Context. Mon R2 is the only ultracompact HII region (UCHII) where the associated photon-dominated region (PDR) can be resolved with Herschel. Due to its brightness and proximity, it is the best source to investigate the chemistry and physics of highly UV-irradiated PDRs. Aims. Our goal is to estimate the abundance of H2O and NH3 in this region and investigate their origin. Methods. We present new observations obtained with HIFI and the IRAM-30m telescope. Using a large velocity gradient approach, we model the line intensities and derive an average abundance of H2O and NH3 across the region. Finally, we model the line profiles with a non-local radiative transfer model and compare these results with the abundance predicted by the Meudon PDR code. Results. The variations of the line profiles and intensities indicate complex geometrical and kinematical patterns. The H2O lines present a strong absorption at the ambient velocity and emission in high velocity wings towards the HII region. The spatial distribution of the o-H2^18O line shows that the its emission arises in the PDR surrounding the HII region. By modeling the o-H2^18O emission we derive a mean abundance of o-H2O of ~10^-8 relative to H2. The ortho-H2O abundance is however larger, ~1x10^-7, in the high velocity wings. Possible explanations for this larger abundance include an expanding hot PDR and/or an outflow. Ammonia seems to be present only in the envelope with an average abundance of ~2x10^-9 relative to H2. Conclusions. The Meudon PDR code can account for the measured water abundance in the high velocity gas as long as we assume that it originates from a <1 mag hot expanding layer of the PDR, i.e. that the outflow has only a minor contribution to this emission. To explain the abundances in the rest of the cloud the molecular freeze out and grain surface chemistry would need to be included. △ Less

Submitted 6 July, 2012; v1 submitted 6 June, 2012; originally announced June 2012.

Comments: 12 pages, 7 figures, 3 tables. Accepted for publication in A&A. Abstract shortened. Updated references, language editing applied in v2

Evaporating Very Small Grains as tracers of the UV radiation field in Photo-dissociation Regions

Authors: P. Pilleri, J. Montillaud, O. Berné, C. Joblin

Abstract: Context. In photo-dissociation regions (PDRs), Polycyclic Aromatic Hydrocarbons (PAHs) could be produced by evaporation of Very Small Grains (VSGs) by the impinging UV radiation field from a nearby star. Aims. We investigate quantitatively the transition zone between evaporating Very Small Grains (eVSGs) and PAHs in several PDRs. Methods. We study the relative contribution of PAHs and eVSGs to the… ▽ More Context. In photo-dissociation regions (PDRs), Polycyclic Aromatic Hydrocarbons (PAHs) could be produced by evaporation of Very Small Grains (VSGs) by the impinging UV radiation field from a nearby star. Aims. We investigate quantitatively the transition zone between evaporating Very Small Grains (eVSGs) and PAHs in several PDRs. Methods. We study the relative contribution of PAHs and eVSGs to the mid-IR emission in a wide range of excitation conditions. We fit the observed mid-IR emission of PDRs by using a set of template band emission spectra of PAHs, eVSGs and gas lines. The fitting tool PAHTAT (PAH Toulouse Astronomical Templates) is made available to the community as an IDL routine. From the results of the fit, we derive the fraction of carbon f_eVSG locked in eVSGs and compare it to the intensity of the local UV radiation field. Results. We show a clear decrease of f_eVSG with increasing intensity of the local UV radiation field, which supports the scenario of photo-destruction of eVSGs. Conversely, this dependence can be used to quantify the intensity of the UV radiation field for different PDRs, including non resolved ones. Conclusions. PAHTAT can be used to trace the intensity of the local UV radiation field in regions where eVSGs evaporate, which correspond to relatively dense (nH = [100, 10^5 ] cm-3) and UV irradiated PDRs (G0 = [100, 5x10^4]) where H2 emits in rotational lines. △ Less

Submitted 26 April, 2012; v1 submitted 20 April, 2012; originally announced April 2012.

Comments: 13 pages, 11 figures. Accepted for publication in A&A. Typos corrected

Gas morphology and energetics at the surface of PDRs: new insights with Herschel observations of NGC 7023

Authors: C. Joblin, P. Pilleri, J. Montillaud, A. Fuente, M. Gerin, O. Berné, V. Ossenkopf, J. Le Bourlot, D. Teyssier, J. R. Goicoechea, F. Le Petit, M. Röllig, M. Akyilmaz, A. O. Benz, F. Boulanger, S. Bruderer, C. Dedes, K. France, R. Güsten, A. Harris, T. Klein, C. Kramer, S. D. Lord, P. G. Martin, J. Martin-Pintado , et al. (15 additional authors not shown)

Abstract: We investigate the physics and chemistry of the gas and dust in dense photon-dominated regions (PDRs), along with their dependence on the illuminating UV field. Using Herschel-HIFI observations, we study the gas energetics in NGC 7023 in relation to the morphology of this nebula. NGC 7023 is the prototype of a PDR illuminated by a B2V star and is one of the key targets of Herschel. Our approach co… ▽ More We investigate the physics and chemistry of the gas and dust in dense photon-dominated regions (PDRs), along with their dependence on the illuminating UV field. Using Herschel-HIFI observations, we study the gas energetics in NGC 7023 in relation to the morphology of this nebula. NGC 7023 is the prototype of a PDR illuminated by a B2V star and is one of the key targets of Herschel. Our approach consists in determining the energetics of the region by combining the information carried by the mid-IR spectrum (extinction by classical grains, emission from very small dust particles) with that of the main gas coolant lines. In this letter, we discuss more specifically the intensity and line profile of the 158 micron (1901 GHz) [CII] line measured by HIFI and provide information on the emitting gas. We show that both the [CII] emission and the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) arise from the regions located in the transition zone between atomic and molecular gas. Using the Meudon PDR code and a simple transfer model, we find good agreement between the calculated and observed [CII] intensities. HIFI observations of NGC 7023 provide the opportunity to constrain the energetics at the surface of PDRs. Future work will include analysis of the main coolant line [OI] and use of a new PDR model that includes PAH-related species. △ Less

Submitted 7 August, 2010; originally announced August 2010.

Comments: Accepted for publication in Astronomy and Astrophysics Letters (Herschel HIFI special issue), 5 pages, 5 figures