Source author record

Paola Caselli

Paola Caselli appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.SR astro-ph.EP astro-ph.CO astro-ph.IM astro-ph.HE physics.chem-ph astro-ph

Catalog footprint

What is connected

79works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2023arXiv

Infall of material onto the filaments in Barnard 5

Aims. We aim to study the structure and kinematics of the two filaments inside the subsonic core Barnard 5 in Perseus using high-resolution ($\approx$ 2400 au) NH3 data and a multi-component fit analysis. Methods. We used observations of NH3 (1,1) and (2,2) inversion transitions using the Very Large Array (VLA) and the Green Bank Telescope (GBT). We smoothed the data to a beam of 8'' to reliably fit multiple velocity components towards the two filamentary structures identified in B5. Results. Along with the core and cloud components, which dominate the flux in the line of sight, we detected two components towards the two filaments showing signs of infall. We also detected two additional components that can possibly trace new material falling into the subsonic core of B5. Conclusions. Following comparison with previous simulations of filament formation scenarios in planar geometry, we conclude that either the formation of the B5 filaments is likely to be rather cylindrically symmetrical or the filaments are magnetically supported. We also estimate infall rates of $1.6\times10^{-4}\,M_\odot\,yr^{-1}$ and $1.8\times10^{-4}\,M_\odot\,yr^{-1}$ (upper limits) for the material being accreted onto the two filaments. At these rates, the filament masses can change significantly during the core lifetime. We also estimate an upper limit of $3.5\times10^{-5}\,M_\odot\,yr^{-1}$ for the rate of possible infall onto the core itself. Accretion of new material onto cores indicates the need for a significant update to current core evolution models, where cores are assumed to evolve in isolation.

preprint2022arXiv

A train of shocks at 3000 au scale? Exploring the clash of an expanding bubble into the NGC 1333 IRAS 4 region. SOLIS XIV

There is evidence that the star formation process is linked to the intricate net of filaments in molecular clouds, which may be also due to gas compression from external triggers. We studied the southern region of the Perseus NGC 1333 molecular cloud, known to be heavily shaped by similar external triggers, to shed light on the process that perturbed the filament where the Class 0 IRAS4 protostars lie. We use new IRAM-NOEMA observations of SiO and CH3OH, both known to trace violent events as shocks, toward IRAS 4A as part of the Large Program Seeds Of Life in Space (SOLIS). We detected three parallel elongated ($>$6000 au) structures, called fingers, with narrow line profiles (~1.5 $km s^{-1}$) peaked at the cloud systemic velocity, tracing gas with high density (5-20 $10^5 cm^{-3}$) and high temperature (80-160 K). They are chemically different, with the northern finger traced by both SiO and CH3OH ([CH3OH]/[SiO]~160-300), while the other two only by SiO ([CH3OH]/[SiO]$<$ 40). Among various possibilities, a train of three shocks, distanced by $>$5000 yr, would be consistent with the observations if a substantial fraction of silicon, frozen onto the grain mantles, is released by the shocks.We suggest that the shock train is due to an expanding gas bubble, coming behind NGC 1333 from the southwest and clashing against the filament, where IRAS 4A lies. Finally, we propose a solution to the two-decades long debate on the nature and origin of the widespread narrow SiO emission observed in the south part of NGC 1333, namely that it is due to unresolved trains of shocks.

preprint2022arXiv

An Interferometric View of H-MM1. I. Direct Observation of NH3 Depletion

Spectral lines of ammonia, NH$_3$, are useful probes of the physical conditions in dense molecular cloud cores. In addition to advantages in spectroscopy, ammonia has also been suggested to be resistant to freezing onto grain surfaces, which should make it a superior tool for studying the interior parts of cold, dense cores. Here we present high-resolution NH$_3$ observations with the Very Large Array (VLA) and Green Bank Telescope (GBT) towards a prestellar core. These observations show an outer region with a fractional NH$_3$ abundance of X(NH$_3$) = (1.975$\pm$0.005)$\times 10^{-8}$ ($\pm 10\%$ systematic), but it also reveals that after all, the X(NH$_3$) starts to decrease above a H$_2$ column density of $\approx 2.6 \times 10^{22}$ cm$^{-2}$. We derive a density model for the core and find that the break-point in the fractional abundance occurs at the density n(H$_2$) $\sim 2\times10^5$ cm$^{-3}$, and beyond this point the fractional abundance decreases with increasing density, following the power law $n^{-1.1}$. This power-law behavior is well reproduced by chemical models where adsorption onto grains dominates the removal of ammonia and related species from the gas at high densities. We suggest that the break-point density changes from core to core depending on the temperature and the grain properties, but that the depletion power law is anyway likely to be close to $n^{-1}$ owing to the dominance of accretion in the central parts of starless cores.

preprint2022arXiv

Astrochemical modelling of infrared dark clouds

Infrared dark clouds (IRDCs) are cold, dense regions of the interstellar medium (ISM) that are likely to represent the initial conditions for massive star formation. It is thus important to study the physical and chemical conditions of IRDCs to provide constraints and inputs for theoretical models of these processes. We aim to determine the astrochemical conditions, especially cosmic ray ionisation rate (CRIR) and chemical age, in different regions of the massive IRDC G28.37+00.07 by comparing observed abundances of multiple molecules and molecular ions with the predictions of astrochemical models. We have computed a series of single-zone astrochemical models with a gas-grain network that systematically explores the parameter space of density, temperature, CRIR, and visual extinction. We have also investigated the effects of choices of CO ice binding energy and temperatures achieved in transient heating of grains when struck by cosmic rays. We selected 10 positions across the IRDC that are known to have a variety of star formation activity. We utilised mid-infrared (MIR) extinction maps and sub-mm emission maps to measure the mass surface densities of these regions, needed for abundance and volume density estimates. The sub-mm emission maps were also used to measure temperatures. We then used IRAM-30m observations of various tracers to estimate column densities and thus abundances. Using estimates of the abundances of CO, HCO$^+$ and N$_2$H$^+$ we find consistency with astrochemical models that have relatively low CRIRs of $ζ\sim10^{-18}$ to $\sim10^{-17}\:{\rm s}^{-1}$, with no evidence for systematic variation with the level of star formation activity. Astrochemical ages are found to be < 1 Myr. We discuss potential sources of systematic uncertainties in these results and the overall implications for IRDC evolutionary history and astrochemical models.(abridged for arXiv)

preprint2022arXiv

CH$_3$CN deuteration in the SVS13-A Class I hot-corino. SOLIS XV

We studied the line emission from CH3CN and its deuterated isotopologue CH$_2$DCN towards the prototypical Class I object SVS13-A, where the deuteration of a large number of species has already been reported. Our goal is to measure the CH$_3$CN deuteration in a Class I protostar, for the first time, in order to constrain the CH$_3$CN formation pathways and the chemical evolution from the early prestellar core and Class 0 to the evolved Class I stages. We imaged CH2DCN towards SVS13-A using the IRAM NOEMA interferometer at 3mm in the context of the Large Program SOLIS (with a spatial resolution of 1.8"x1.2"). The NOEMA images have been complemented by the CH$_3$CN and CH$_2$DCN spectra collected by the IRAM-30m Large Program ASAI, that provided an unbiased spectral survey at 3mm, 2mm, and 1.3mm. The observed line emission has been analysed using LTE and non-LTE LVG approaches. The NOEMA/SOLIS images of CH2DCN show that this species emits in an unresolved area centered towards the SVS13-A continuum emission peak, suggesting that methyl cyanide and its isotopologues are associated with the hot corino of SVS13-A, previously imaged via other iCOMs. In addition, we detected 41 and 11 ASAI transitions of CH$_3$CN and CH2DCN, respectively, which cover upper level energies (Eup) from 13 to 442 K and from 18 K to 200 K, respectively. The derived [CH2DCN]/[CH3CN] ratio is $\sim$9\%. This value is consistent with those measured towards prestellar cores and a factor 2-3 higher than those measured in Class 0 protostars. Contrarily to what expected for other molecular species, the CH3CN deuteration does not show a decrease in SVS13-A with respect to measurements in younger prestellar cores and Class 0 protostars. Finally, we discuss why our new results suggest that CH3CN was likely synthesised via gas-phase reactions and frozen onto the dust grain mantles during the cold prestellar phase.

preprint2022arXiv

FAUST III. Misaligned rotations of the envelope, outflow, and disks in the multiple protostellar system of VLA 1623$-$2417

We report a study of the low-mass Class-0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H$^{13}$CO$^+$ ($J=3-2$), CS ($J=5-4$), and CCH ($N=3-2$) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the rotation of the circum-binary VLA 1623A disk as well as the VLA 1623B disk. We found that the minor axis of the circum-binary disk of VLA 1623A is misaligned by about 12 degrees with respect to the large-scale outflow and the rotation axis of the envelope. In contrast, the minor axis of the circum-binary disk is parallel to the large-scale magnetic field according to previous dust polarization observations, suggesting that the misalignment may be caused by the different directions of the envelope rotation and the magnetic field. If the velocity gradient of the outflow is caused by rotation, the outflow has a constant angular momentum and the launching radius is estimated to be $5-16$ au, although it cannot be ruled out that the velocity gradient is driven by entrainments of the two high-velocity outflows. Furthermore, we detected for the first time a velocity gradient associated with rotation toward the VLA 16293B disk. The velocity gradient is opposite to the one from the large-scale envelope, outflow, and circum-binary disk. The origin of its opposite gradient is also discussed.

preprint2022arXiv

First detection of CHD$_{2}$OH towards pre-stellar cores

The inheritance of material across the star and planet formation process is traced by deuterium fractionation. We report here the first detection of doubly deuterated methanol towards pre-stellar cores. We study the deuterium fractionation of methanol, CH$_{3}$OH, towards two starless and two pre-stellar cores. We derive a D/H ratio of 0.8-1.9$\%$ with CH$_{2}$DOH in pre-stellar cores H-MM1 and L694-2, consistent with measurements in more evolved Class 0/I objects and comet 67P/Churyumov-Gerasimenko, suggesting a direct chemical link arising in the pre-stellar stage. Furthermore, the column density ratios of CHD$_{2}$OH/CH$_{2}$DOH are $\sim$50-80$\%$, consistently high as that towards Class 0/I objects, indicating an efficient formation mechanism of CHD$_{2}$OH, possibly through H atom additions to D$_{2}$CO. The CH$_{2}$DOH/CH$_{3}$OH and CHD$_{2}$OH/CH$_{3}$OH column density ratios in the two pre-stellar cores are larger than that in the two starless cores B68 and L1521E, representing an evolutionary trend of methanol deuteration in early-stage cores.

preprint2022arXiv

Nitrogen fractionation towards a pre-stellar core traces isotope-selective photodissociation

Isotopologue abundance ratios are important to understand the evolution of astrophysical objects and ultimately the origins of a planetary system like our own. Being nitrogen a fundamental ingredient of pre-biotic material, understanding its chemistry and inheritance is of fundamental importance to understand the formation of the building blocks of life. We present here single-dish observations of the ground state rotational transitions of the $^{13}$C and $^{15}$N isotopologues of HCN, HNC and CN with the IRAM 30m telescope. We analyse their column densities and compute the $^{14}$N/$^{15}$N ratio map for HCN. The $^{15}$N-fractionation of CN and HNC is computed towards different offsets across L1544. The $^{15}$N-fractionation map of HCN shows a clear decrease of the $^{14}$N/$^{15}$N ratio towards the southern edge of L1544, where carbon chain molecules present a peak, strongly suggesting that isotope-selective photodissociation has a strong effect on the fractionation of nitrogen across pre-stellar cores. The $^{14}$N/$^{15}$N ratio in CN measured towards four positions across the core also shows a decrease towards the South-East of the core, while HNC shows opposite behaviour. The uneven illumination of the pre-stellar core L1544 provides clear evidence that $^{15}$N-fractionation of HCN and CN is enhanced toward the region more exposed to the interstellar radiation field. Isotope-selective photodissociation of N$_2$ is then a crucial process to understand $^{15}$N fractionation, as already found in protoplanetary disks. Therefore, the $^{15}$N-fractionation in pre-stellar material is expected to change depending on the environment within which pre-stellar cores are embedded. The $^{12}$CN/$^{13}$CN ratio also varies across the core, but its variation does not affect our conclusions on the effect of the environment on the fractionation of nitrogen.

preprint2022arXiv

Seeds of Life in Space (SOLIS). XIII. Nitrogen fractionation towards the protocluster OMC-2 FIR4

Isotopic fractionation is an important tool to investigate the chemical history of our Solar System (SS). In particular, the isotopic fraction of nitrogen (14N/15N) is lower in comets and other pristine SS bodies with respect to the value measured for the protosolar nebula, suggesting a local chemical enrichment of 15N during the SS formation. Therefore, interferometric studies of nitrogen fractionation in SS precursors are imperative for us to obtain clues about our astrochemical origins. In this work, we investigated the variation of the 14N/15N ratio in one of the closest analogues of the environment in which the SS was born: the protocluster OMC-2 FIR4. We present the first comparison at high angular resolution between HCN and N2H+ using interferometric data. We analysed observations of the HCN isotopologues H13CN and HC15N in the OMC-2 FIR4 protocluster, specifically the transitions H13CN (1-0) and HC15N (1-0), from NOEMA within the context of the IRAM Seeds Of Life In Space Large Program. We combined our results with analysis of archival data obtained with ALMA of N2H+ and its 15N isotopologues. Our results show a small regional variation in the ratio for HCN from ~250 to 500. The ratios in the central regions of FIR4, where the candidate protostars are located, are largely consistent (~300). They also show little variation from the part of the protocluster known to harbour a high cosmic-ray ionisation rate, to the portion with lower rate. We also found a small variation in the ratio of N2H+ across different regions from ~200 to ~400. These results suggest that local changes in the physical parameters occurring on small linear scales probed by our observations do not seem to affect the 14N/15N ratio in either HCN or N2H+ and hence that this is independent of the molecule used. Moreover, the high level of irradiation due to cosmic rays does not affect the N-fractionation either.

preprint2022arXiv

The Central 1000 au of a Pre-stellar Core Revealed with ALMA. II. Almost Complete Freeze-out

Pre-stellar cores represent the initial conditions in the process of star and planet formation. Their low temperatures ($<$10 K) allow the formation of thick icy dust mantles, which will be partially preserved in the future protoplanetary disks, ultimately affecting the chemical composition of planetary systems. Previous observations have shown that carbon- and oxygen-bearing species, in particular CO, are heavily depleted in pre-stellar cores due to the efficient molecular freeze-out onto the surface of cold dust grains. However, N-bearing species such as NH$_3$ and, in particular, its deuterated isotopologues, appear to maintain high abundances where CO molecules are mainly in solid phase. Thanks to ALMA, we present here the first clear observational evidence of NH$_2$D freeze-out toward the L1544 pre-stellar core, suggestive of the presence of a"complete-depletion zone" within a $\simeq$1800 au radius, in agreement with astrochemical pre-stellar core model predictions. Our state-of-the-art chemical model coupled with a non-LTE radiative transfer code demonstrates that NH$_2$D becomes mainly incorporated in icy mantles in the central 2000 au and starts freezing-out already at $\simeq$7000 au. Radiative transfer effects within the pre-stellar core cause the NH$_2$D(1$_{11}$-1$_{01}$) emission to appear centrally concentrated, with a flattened distribution within the central $\simeq$3000 au, unlike the 1.3 mm dust continuum emission which shows a clear peak within the central $\simeq$1800 au. This prevented NH$_2$D freeze-out to be detected in previous observations, where the central 1000 au cannot be spatially resolved.

preprint2022arXiv

Tracking the ice mantle history in the Solar-type Protostars of NGC 1333 IRAS 4

To understand the origin of the diversity observed in exoplanetary systems, it is crucial to characterize the early stages of their formation, represented by Solar-type protostars. Likely, the gaseous chemical content of these objects directly depends on the composition of the dust grain mantles formed before the collapse. Directly retrieving the ice mantle composition is challenging, but it can be done indirectly by observing the major components, such as NH3 and CH3OH at cm wavelengths, once they are released into the gas-phase during the warm protostellar stage. We observed several CH3OH and NH3 lines toward three Class 0 protostars in NGC1333 (IRAS 4A1, IRAS 4A2, and IRAS 4B), at high angular resolution (1"; ~300 au) with the VLA interferometer at 24-26 GHz. Using a non-LTE LVG analysis, we derived a similar NH3/CH3OH abundance ratio in the three protostars (<0.5, 0.015-0.5, and 0.003-0.3 for IRAS 4A1, 4A2, and 4B, respectively). Hence, we infer they were born from pre-collapse material with similar physical conditions. Comparing the observed abundance ratios with astrochemical model predictions, we constrained the dust temperature at the time of the mantle formation to be ~17 K, which coincides with the average temperature of the southern NGC 1333 diffuse cloud. We suggest that a brutal event started the collapse that eventually formed IRAS 4A1, 4A2 and 4B, which,therefore, did not experience the usual pre-stellar core phase. This event could be the clash of a bubble with NGC 1333 south, that has previously been evoked in the literature.

preprint2021arXiv

Deuterium Chemodynamics of Massive Pre-Stellar Cores

High levels of deuterium fractionation of $\rm N_2H^+$ (i.e., $\rm D_{frac}^{N_2H^+} \gtrsim 0.1$) are often observed in pre-stellar cores (PSCs) and detection of $\rm N_2D^+$ is a promising method to identify elusive massive PSCs. However, the physical and chemical conditions required to reach such high levels of deuteration are still uncertain, as is the diagnostic utility of $\rm N_2H^+$ and $\rm N_2D^+$ observations of PSCs. We perform 3D magnetohydrodynamics simulations of a massive, turbulent, magnetised PSC, coupled with a sophisticated deuteration astrochemical network. Although the core has some magnetic/turbulent support, it collapses under gravity in about one freefall time, which marks the end of the simulations. Our fiducial model achieves relatively low $\rm D_{frac}^{N_2H^+} \sim 0.002$ during this time. We then investigate effects of initial ortho-para ratio of $\rm H_2$ ($\rm OPR^{H_2}$), temperature, cosmic ray (CR) ionization rate, CO and N-species depletion factors and prior PSC chemical evolution. We find that high CR ionization rates and high depletion factors allow the simulated $\rm D_{frac}^{N_2H^+}$ and absolute abundances to match observational values within one freefall time. For $\rm OPR^{H_2}$, while a lower initial value helps the growth of $\rm D_{frac}^{N_2H^+}$, the spatial structure of deuteration is too widespread compared to observed systems. For an example model with elevated CR ionization rates and significant heavy element depletion, we then study the kinematic and dynamic properties of the core as traced by its $\rm N_2D^+$ emission. The core, undergoing quite rapid collapse, exhibits disturbed kinematics in its average velocity map. Still, because of magnetic support, the core often appears kinematically sub-virial based on its $\rm N_2D^+$ velocity dispersion.

preprint2021arXiv

Dissecting the super-critical filaments embedded in the 0.5 pc subsonic region of Barnard 5

We characterize in detail the two ~0.3 pc long filamentary structures found within the subsonic region of Barnard 5. We use combined GBT and VLA observations of the molecular lines NH$_3$(1,1) and (2,2) at a resolution of 1800 au, as well as JCMT continuum observations at 850 and 450 $μ$m at a resolution of 4400 au and 3000 au, respectively. We find that both filaments are highly super-critical with a mean mass per unit length, $M/L$, of ~80 M$_\odot$ pc$^{-1}$, after background subtraction, with local increases reaching values of ~150 M$_\odot$ pc$^{-1}$. This would require a magnetic field strength of ~500 $μ$G to be stable against radial collapse. We extract equidistant cuts perpendicular to the spine of the filament and fit a modified Plummer profile as well as a Gaussian to each of the cuts. The filament widths (deconvolved FWHM) range between 6500-7000 au (~0.03 pc) along the filaments. This equals ~2.0 times the radius of the flat inner region. We find an anti-correlation between the central density and this flattening radius, suggestive of contraction. Further, we also find a strong correlation between the power-law exponent at large radii and the flattening radius. We note that the measurements of these three parameters fall in a plane and derive their empirical relation. Our high-resolution observations provide direct constraints of the distribution of the dense gas within super-critical filaments showing pre- and protostellar activity.

preprint2021arXiv

Methanol Mapping in Cold Cores: Testing Model Predictions

Chemical models predict that in cold cores gas-phase methanol is expected to be abundant at the outer edge of the CO depletion zone, where CO is actively adsorbed. CO adsorption correlates with volume density in cold cores, and, in nearby molecular clouds, the catastrophic CO freeze-out happens at volume densities above 10$^4$ cm$^{-3}$. The methanol production rate is maximized there and its freeze-out rate does not overcome its production rate, while the molecules are shielded from UV destruction by gas and dust. Thus, in cold cores, methanol abundance should generally correlate with visual extinction that depends both on volume and column density. In this work, we test the most basic model prediction that maximum methanol abundance is associated with a local $A_V\simeq$4 mag in dense cores and constrain the model parameters with the observational data. With the IRAM 30 m antenna, we mapped the CH$_3$OH (2-1) and (3-2) transitions toward seven dense cores in the L1495 filament in Taurus to measure the methanol abundance. We use the Herschel/SPIRE maps to estimate visual extinction, and the C$^{18}$O(2-1) maps from Tafalla & Hacar (2015) to estimate CO depletion. We explored the observed and modeled correlations between the methanol abundances, CO depletion, and visual extinction varying the key model parameters. The modeling results show that hydrogen surface diffusion via tunneling is crucial to reproduce the observed methanol abundances, and the needed reactive desorption efficiency matches the one deduced from laboratory experiments.

preprint2020arXiv

A new proxy to estimate the cosmic-ray ionisation rate in dense cores

Cosmic rays are a global source of ionisation, and the ionisation fraction represents a fundamental parameter in the interstellar medium. Ions couple to magnetic fields, affect the chemistry, and the dynamics of star-forming regions as well as planetary atmospheres. However, the cosmic-ray ionisation rate represents one of the bottlenecks for astrochemical models, and its determination is one of the most puzzling problems in astrophysics. While for diffuse clouds reasonable values have been provided from H$_3^+$ observations, for dense clouds, due to the lack of rotational transitions, this is not possible, and estimates are strongly biased by the employed model. We present here an analytical expression, obtained from first principles, to estimate the cosmic-ray ionisation rate from observational quantities. The theoretical predictions are validated with high-resolution three-dimensional numerical simulations and applied to the well known core L1544; we obtained an estimate of $ζ_2 \sim 2-3 \times 10^{-17}$ s$^{-1}$. Our results and the analytical formulae provided represent the first model-independent, robust tool to probe the cosmic-ray ionisation rate in the densest part of star-forming regions (on spatial scales of $R \leq 0.05$ pc). An error analysis is presented to give statistical relevance to our study.

preprint2020arXiv

A protostellar system fed by a streamer of 10,500 au length

Binary formation is an important aspect of star formation. One possible route for close-in binary formation is disk fragmentation$^{[1,2,3]}$. Recent observations show small scale asymmetries (<300 au) around young protostars$^{[2,4]}$, although not always resolving the circumbinary disk, are linked to disk phenomena$^{[5,6]}$. In later stages, resolved circumbinary disk observations$^{[7]}$ (<200 au) show similar asymmetries, suggesting the origin of the asymmetries arises from binary-disk interactions$^{[8,9,10]}$. We observed one of the youngest systems to study the connection between disk and dense core. We find for the first time a bright and clear streamer in chemically fresh material (Carbon-chain species) that originates from outside the dense core (>10,500 au). This material connects the outer dense core with the region where asymmetries arise near disk scales. This new structure type, 10x larger than those seen near disk scales, suggests a different interpretation of previous observations: large-scale accretion flows funnel material down to disk scales. These results reveal the under-appreciated importance of the local environment on the formation and evolution of disks in early systems$^{[13,14]}$ and a possible initial condition for the formation of annular features in young disks$^{[15,16]}$.

preprint2020arXiv

Carbon Chain Chemistry in Hot-Core Regions around Three Massive Young Stellar Objects Associated with 6.7 GHz Methanol Masers

We have carried out observations of CCH ($N=1-0$), CH$_{3}$CN ($J=5-4$), and three $^{13}$C isotopologues of HC$_{3}$N ($J=10-9$) toward three massive young stellar objects (MYSOs), G12.89+0.49, G16.86--2.16, and G28.28--0.36, with the Nobeyama 45-m radio telescope. Combined with previous results on HC$_{5}$N, the column density ratios of $N$(CCH)/$N$(HC$_{5}$N), hereafter the CCH/HC$_{5}$N ratios, in the MYSOs are derived to be $\sim 15$. This value is lower than that in a low-mass warm carbon chain chemistry (WCCC) source by more than one order of magnitude. We compare the observed CCH/HC$_{5}$N ratios with hot-core model calculations (Taniguchi et al. 2019). The observed ratios in the MYSOs can be best reproduced by models when the gas temperature is $\sim 85$ K, which is higher than in L1527, a low-mass WCCC source ($\sim 35$ K). These results suggest that carbon-chain molecules detected around the MYSOs exist at least partially in higher temperature regions than those in low-mass WCCC sources. There is no significant difference in column density among the three $^{13}$C isotopologues of HC$_{3}$N in G12.89+0.49 and G16.86-2.16, while HCC$^{13}$CN is more abundant than the others in G28.28--0.36. We discuss carbon-chain chemistry around the three MYSOs based on the CCH/HC$_{5}$N ratio and the $^{13}$C isotopic fractionation of HC$_{3}$N.

preprint2020arXiv

Efficient methanol production on the dark side of a prestellar core

We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particularly efficient. Methanol and sulphur monoxide are most abundant in a narrow zone that follows the eastern side of the core. This side is sheltered from the stronger external radiation field coming from the west. We show that photodissociation on the illuminated side can give rise to an asymmetric methanol distribution, but that the stark contrast observed in H-MM1 is hard to explain without assuming enhanced desorption on the shaded side. The region of the brightest emission has a wavy structure that rolls up at one end. This is the signature of Kelvin-Helmholtz instability occurring in sheared flows. We suggest that in this zone, methanol and sulphur are released as a result of grain-grain collisions induced by shear vorticity.

preprint2020arXiv

Hall Effect in Protostellar Disc Formation and Evolution

The Hall effect is recently shown to be efficient in magnetized dense molecular cores, and could lead to a bimodal formation of rotationally supported discs (RSDs) in the first core phase. However, how such Hall dominated systems evolve in the protostellar accretion phase remains unclear. We carry out 2D axisymmetric simulations including Hall effect and Ohmic dissipation, with realistic magnetic diffusivities computed from our equilibrium chemical network. We find that Hall effect only becomes efficient when the large population of very small grains (VSGs: $\lesssim$10 nm) is removed from the standard MRN size distribution. With such an enhanced Hall effect, however, the bimodality of disc formation does not continue into the main accretion phase. The outer part of the initial $\sim$40 AU disc formed in the anti-aligned configuration (${\bf Ω\cdot B}<0$) flattens into a thin rotationally supported Hall current sheet as Hall effect moves the poloidal magnetic field radially inward relative to matter, leaving only the inner $\lesssim$10--20 AU RSD. In the aligned configuration (${\bf Ω\cdot B}>0$), disc formation is suppressed initially but a counter-rotating disc forms subsequently due to efficient azimuthal Hall drift. The counter-rotating disc first grows to $\sim$30 AU as Hall effect moves the magnetic field radially outward, but only the inner $\lesssim$10 AU RSD is long-lived like in the anti-aligned case. Besides removing VSGs, cosmic ray ionization rate should be below a few 10$^{-16}$ s$^{-1}$ for Hall effect to be efficient in disc formation. We conclude that Hall effect produces small $\lesssim$10--20 AU discs regardless of the polarity of the magnetic field, and that radially outward diffusion of magnetic fields remains crucial for disc formation and growth.

preprint2020arXiv

Hot Corinos Chemical Diversity: Myth or Reality?

After almost 20 years of hunting, only about a dozen hot corinos, hot regions enriched in interstellar complex organic molecules (iCOMs), are known. Of them, many are binary systems with the two components showing drastically different molecular spectra. Two obvious questions arise. Why are hot corinos so difficult to find and why do their binary components seem chemically different? The answer to both questions could be a high dust opacity that would hide the molecular lines. To test this hypothesis, we observed methanol lines at centimeter wavelengths, where dust opacity is negligible, using the Very Large Array interferometer. We targeted the NGC 1333 IRAS 4A binary system, for which one of the two components, 4A1, has a spectrum deprived of iCOMs lines when observed at millimeter wavelengths, while the other component, 4A2, is very rich in iCOMs. We found that centimeter methanol lines are similarly bright toward 4A1 and 4A2. Their non-LTE analysis indicates gas density and temperature ($\geq2\times10^6$ cm$^{-3}$ and 100--190 K), methanol column density ($\sim10^{19}$ cm$^{-2}$) and extent ($\sim$35 au in radius) similar in 4A1 and 4A2, proving that both are hot corinos. Furthermore, the comparison with previous methanol line millimeter observations allows us to estimate the optical depth of the dust in front of 4A1 and 4A2, respectively. The obtained values explain the absence of iCOMs line emission toward 4A1 at millimeter wavelengths and indicate that the abundances toward 4A2 are underestimated by $\sim$30\%. Therefore, centimeter observations are crucial for the correct study of hot corinos, their census, and their molecular abundances.

preprint2020arXiv

Impact of magneto-rotational instability on grain growth in protoplanetary disks: I. Relevant turbulence properties

Turbulence in the protoplanetary disks induces collisions between dust grains, and thus facilitates grain growth. We investigate the two fundamental assumptions of the turbulence in obtaining grain collisional velocities -- the kinetic energy spectrum and the turbulence autocorrelation time -- in the context of the turbulence generated by the magneto-rotational instability (MRI). We carry out numerical simulations of the MRI as well as driven turbulence, for a range of physical and numerical parameters. We find that the convergence of the turbulence $α$-parameter does not necessarily imply the convergence of the energy spectrum. The MRI turbulence is largely solenoidal, for which we observe a persistent kinetic energy spectrum of $k^{-4/3}$. The same is obtained for solenoidal driven turbulence with and without magnetic field, over more than 1 dex near the dissipation scale. This power-law slope appears to be converged in terms of numerical resolution, and to be due to the bottleneck effect. The kinetic energy in the MRI turbulence peaks at the fastest growing mode of the MRI. In contrast, the magnetic energy peaks at the dissipation scale. The magnetic energy spectrum in the MRI turbulence does not show a clear power-law range, and is almost constant over approximately 1 dex near the dissipation scale. The turbulence autocorrelation time is nearly constant at large scales, limited by the shearing timescale, and shows a power-law drop close to $k^{-1}$ at small scales, with a slope steeper than that of the eddy crossing time. The deviation from the standard picture of the Kolmogorov turbulence with the injection scale at the disk scale height can potentially have a significant impact on the grain collisional velocities.

preprint2020arXiv

Orbital and mass constraints of the young binary system IRAS 16293-2422 A

We present 3 mm ALMA continuum and line observations at resolutions of 6.5 au and 13 au respectively, toward the Class 0 system IRAS 16293-2422 A. The continuum observations reveal two compact sources towards IRAS 16293-2422 A, coinciding with compact ionized gas emission previously observed at radio wavelengths (A1 and A2), confirming the long-known radio sources as protostellar. The emission towards A2 is resolved and traces a dust disk with a FWHM size of ~12 au, while the emission towards A1 sets a limit to the FWHM size of the dust disk of ~4 au. We also detect spatially resolved molecular kinematic tracers near the protostellar disks. Several lines of the J=5-4 rotational transition of HNCO, NH2CHO and t-HCOOH are detected, with which we derived individual line-of-sight velocities. Using these together with the CS (J=2-1), we fit Keplerian profiles towards the individual compact sources and derive masses of the central protostars. The kinematic analysis indicates that A1 and A2 are a bound binary system. Using this new context for the previous 30 years of VLA observations, we fit orbital parameters to the relative motion between A1 and A2 and find the combined protostellar mass derived from the orbit is consistent with the masses derived from the gas kinematics. Both estimations indicate masses consistently higher (0.5< M1<M2<2 Msun) than previous estimations using lower resolution observations of the gas kinematics. The ALMA high-resolution data provides a unique insight into the gas kinematics and masses of a young deeply embedded bound binary system.

preprint2020arXiv

Probabilistic detection of spectral line components

Resolved kinematical information, such as from molecular gas in star forming regions, is obtained from spectral line observations. However, these observations often contain multiple line-of-sight components, making estimates harder to obtain and interpret. We present a fully automatic method that determines the number of components along the line of sight, or the spectral multiplicity, through Bayesian model selection. The underlying open-source framework, based on nested sampling and conventional spectral line modeling, is tested using the large area ammonia maps of NGC 1333 in Perseus molecular cloud obtained by the Green Bank Ammonia Survey (GAS). Compared to classic approaches, the presented method constrains velocities and velocity dispersions in a larger area. In addition, we find that the velocity dispersion distribution among multiple components did not change substantially from that of a single fit component analysis of the GAS data. These results showcase the power and relative ease of the fitting and model selection method, which makes it a unique tool to extract maximum information from complex spectral data.

preprint2020arXiv

Propargylimine in the laboratory and in space: millimetre-wave spectroscopy and first detection in the ISM

Small imines containing up to three carbon atoms are present in the interstellar medium. As alkynyl compounds are abundant in this medium, propargylimine thus represents a promising candidate for a new interstellar detection. The goal of the present work is to perform a comprehensive laboratory investigation of the rotational spectrum of propargylimine in its ground vibrational state in order to obtain a highly precise set of rest-frequencies and to search it in space. The rotational spectra of $E$ and $Z$ geometrical isomers of propargylimine have been recorded in laboratory in the 83-500 GHz frequency interval. The measurements have been performed using a source-modulation millimetre-wave spectrometer equipped with a pyrolysis system for the production of unstable species. High-level ab initio calculations were performed to assist the analysis and to obtain reliable estimates for an extended set of spectroscopic quantities. The improved spectral data allow us to perform a successful search for this new imine in the quiescent G+0.693-0.027 molecular cloud. Eighteen lines of $Z$-propargylimine have been detected at level $>2.5 σ$, resulting in a column density estimate of $N = (0.24\pm 0.02)\times 10^{14}$ cm$^{-2}$. An upper limit was retrieved for the higher-energy $E$ isomer, which was not detected in the data. The fractional abundance (w.r.t. H$_2$) derived for $Z$-propargylimine is $1.8\times 10^{-10}$. We discuss the possible formation routes by comparing the derived abundance with those measured in the source for possible chemical precursors.

preprint2020arXiv

Relative Alignment between Dense Molecular Cores and Ambient Magnetic Field: The Synergy of Numerical Models and Observations

The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc- to core-scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flow along the magnetic field toward dense cores. When comparing the observed cores identified from the GBT Ammonia Survey (GAS) and Planck polarization-inferred magnetic field orientations, we find that the relative core-field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core-field orientation could be used to probe the relative significance of the magnetic field within the cloud.

preprint2020arXiv

Seeds of Life in Space (SOLIS). VI. Chemical evolution of sulfuretted species along the outflows driven by the low-mass protostellar binary NGC1333-IRAS4A

Context. Low-mass protostars drive powerful molecular outflows that can be observed with mm and sub-mm telescopes. Various sulfuretted species are known to be bright in shocks and could be used to infer the physical and chemical conditions throughout the observed outflows. Aims. The evolution of sulfur chemistry is studied along the outflows driven by the NGC1333-IRAS4A protobinary system located in the Perseus cloud to constrain the physical and chemical processes at work in shocks. Methods. We observed various transitions from OCS, CS, SO, and SO$_2$ towards NGC1333-IRAS4A in the 1.3, 2, and 3mm bands using the IRAM NOEMA array and we interpreted the observations through the use of the Paris-Durham shock model. Results. The targeted species clearly show different spatial emission along the two outflows driven by IRAS4A. OCS is brighter on small and large scales along the south outflow driven by IRAS4A1, whereas SO$_2$ is detected rather along the outflow driven by IRAS4A2 that is extended along the north east - south west (NE-SW) direction. Column density ratio maps estimated from a rotational diagram analysis allowed us to confirm a clear gradient of the OCS/SO$_2$ column density ratio between the IRAS4A1 and IRAS4A2 outflows. SO is detected at extremely high radial velocity up to 25 km/s relative to the source velocity, clearly allowing us to distinguish the two outflows on small scales. Conclusions. The observed chemical differentiation between the two outflows of the IRAS4A system could be explained by a different chemical history. The outflow driven by IRAS4A1 is likely younger and more enriched in species initially formed in interstellar ices, such as OCS, and recently sputtered into the shock gas. In contrast, the longer and likely older outflow triggered by IRAS4A2 is more enriched in species that have a gas phase origin, such as SO$_2$.

preprint2020arXiv

Submillimetre and far-infrared spectroscopy of monodeuterated amidogen radical (NHD): improved rest-frequencies for astrophysical observations

Observations of ammonia in interstellar environments have revealed high levels of deuteration, and all its D-containing variants, including ND$_3$, have been detected in cold prestellar cores and around young protostars. The observation of these deuterated isotopologues is very useful to elucidate the chemical and physical processes taking place during the very early stages of star formation, as the abundance of deuterated molecules is highly enhanced in dense and cold gas. Nitrogen hydride radicals are key species lying at the very beginning of the reaction pathway leading to the formation of NH$_3$ and organic molecules of pre-biotic interest, but relatively little information is known about their D-bearing isotopologues. To date, only ND has been detected in the interstellar gas. To aid the identification of further deuterated nitrogen radicals, we have thoroughly re-investigated the rotational spectrum of NHD employing two different instruments: a frequency-modulation submillimetre spectrometer operating in the THz region and a synchrotron-based Fourier Transform infrared spectrometer operating in the 50-240 cm$^{-1}$ wavelength range. NHD was produced in a plasma of NH$_3$ and D$_2$. A wide range of rotational energy levels has been probed thanks to the observation of high $N$ (up to 15) and high $K_a$ (up to 9) transitions. A global analysis including our new data and those already available in the literature has provided a comprehensive set of very accurate spectroscopic parameters. A highly reliable line catalogue has been generated to assist archival data searches and future astronomical observations of NHD at submillimetre and THz regimes.

preprint2020arXiv

Ubiquitous $\rm NH_3$ supersonic component in L1688 coherent cores

Context : Star formation takes place in cold dense cores in molecular clouds. Earlier observations have found that dense cores exhibit subsonic non-thermal velocity dispersions. In contrast, CO observations show that the ambient large-scale cloud is warmer and has supersonic velocity dispersions. Aims : We aim to study the ammonia ($\rm NH_3$) molecular line profiles with exquisite sensitivity towards the coherent cores in L1688 in order to study their kinematical properties in unprecedented detail. Methods : We used $\rm NH_3$ (1,1) and (2,2) data from the first data release (DR1) in the Green Bank Ammonia Survey (GAS). We first smoothed the data to a larger beam of 1' to obtain substantially more extended maps of velocity dispersion and kinetic temperature, compared to the DR1 maps. We then identified the coherent cores in the cloud and analysed the averaged line profiles towards the cores. Results : For the first time, we detected a faint (mean $\rm NH_3$(1,1) peak brightness $<$0.25 K in $T_{MB}$), supersonic component towards all the coherent cores in L1688. We fitted two components, one broad and one narrow, and derived the kinetic temperature and velocity dispersion of each component. The broad components towards all cores have supersonic linewidths ($\mathcal{M}_S \ge 1$). This component biases the estimate of the narrow dense core component's velocity dispersion by $\approx$28% and the kinetic temperature by $\approx$10%, on average, as compared to the results from single-component fits. Conclusions : Neglecting this ubiquitous presence of a broad component towards all coherent cores causes the typical single-component fit to overestimate the temperature and velocity dispersion. This affects the derived detailed physical structure and stability of the cores estimated from $\rm NH_3$ observations.

preprint2020arXiv

Velocity-coherent Filaments in NGC 1333: Evidence for Accretion Flow?

Recent observations of global velocity gradients across and along molecular filaments have been interpreted as signs of gas accreting onto and along these filaments, potentially feeding star-forming cores and proto-clusters. The behavior of velocity gradients in filaments, however, has not been studied in detail, particularly on small scales (< 0.1 pc). In this paper, we present MUFASA, an efficient, robust, and automatic method to fit ammonia lines with multiple velocity components, generalizable to other molecular species. We also present CRISPy, a Python package to identify filament spines in 3D images (e.g., position-position-velocity cubes), along with a complementary technique to sort fitted velocity components into velocity-coherent filaments. In NGC 1333, we find a wealth of velocity gradient structures on a beam-resolved scale of ~0.05 pc. Interestingly, these local velocity gradients are not randomly oriented with respect to filament spines and their perpendicular, i.e., radial, component decreases in magnitude towards the spine for many filaments. Together with remarkably constant velocity gradients on larger scales along many filaments, these results suggest a scenario in which gas falling onto filaments is progressively damped and redirected to flow along these filaments.

preprint2019arXiv

Efficient Production of S$_8$ in Interstellar Ices: The effects of cosmic ray-driven radiation chemistry and non-diffusive bulk reactions

In this work, we reexamine sulfur chemistry occurring on and in the ice mantles of interstellar dust grains, and report the effects of two new modifications to standard astrochemical models; namely, (a) the incorporation of cosmic ray-driven radiation chemistry and (b) the assumption of fast, non-diffusive reactions for key radicals in the bulk. Results from our models of dense molecular clouds show that these changes can have a profound influence on the abundances of sulfur-bearing species in ice mantles, including a reduction in the abundance of solid-phase H$_2$S and HS, and a significant increase in the abundances of OCS, SO$_2$, as well as pure allotropes of sulfur, especially S$_8$. These pure-sulfur species - though nearly impossible to observe directly - have long been speculated to be potential sulfur reservoirs and our results represent possibly the most accurate estimates yet of their abundances in the dense ISM. Moreover, the results of these updated models are found to be in good agreement with available observational data. Finally, we examine the implications of our findings with regard to the as-yet-unknown sulfur reservoir thought to exist in dense interstellar environments.

preprint2019arXiv

The first steps of Interstellar Phosphorus Chemistry

Phosphorus-bearing species are an essential key to form life on Earth, however they have barely been detected in the interstellar medium. Since only PN and PO have been identified so far towards star-forming regions, the chemical formation pathways of P-bearing molecules are not easy to constrain and are thus highly debatable. An important factor still missing in the chemical models is the initial elemental abundance of phosphorus, i.e. the depletion level of P. In order to overcome this problem, we study P-bearing species in diffuse/translucent clouds. In these objects phosphorus is mainly in the gas phase and therefore the elemental initial abundance needed in our chemical simulations corresponds to the cosmic one and is thus well constrained. An advanced chemical model with an updated P-chemistry network has been used. Single-pointing observations were performed with the IRAM 30m telescope towards the line of sight to the blazar B0355+508 aiming for the (2-1) transitions of PN, PO, HCP and CP. This line of sight incorporates five diffuse/translucent clouds. The (2-1) transitions of the PN, PO, HCP and CP were not detected. We report detections of the (1-0) lines of $\mathrm{^{13}CO}$, HNC and CN along with a first detection of $\mathrm{C^{34}S}$. We have reproduced the observations of HNC, CN, CS and CO in every cloud by applying typical conditions for diffuse/translucent clouds. According to our best-fit model, the most abundant P-bearing species are HCP and CP ($\sim10^{-10}$), followed by PN, PO and $\mathrm{PH_3}$ ($\sim10^{-11}$). We show that the production of P-bearing species is favoured towards translucent rather than diffuse clouds, where the environment provides a stronger shielding from the interstellar radiation. Based on our improved model, the (1-0) transitions of HCP, CP, PN and PO are expected to be detectable with estimated intensities up to $\sim200$ mK.

preprint2016arXiv

An Ordered Bipolar Outflow from a Massive Early-Stage Core

We present ALMA follow-up observations of two massive, early-stage core candidates, C1-N & C1-S, in Infrared Dark Cloud (IRDC) G028.37+00.07, which were previously identified by their N2D+(3-2) emission and show high levels of deuteration of this species. The cores are also dark at far infrared wavelengths up to ~100 microns. We detect 12CO(2-1) from a narrow, highly-collimated bipolar outflow that is being launched from near the center of the C1-S core, which is also the location of the peak 1.3mm dust continuum emission. This protostar, C1-Sa, has associated dense gas traced by C18O(2-1) and DCN(3-2), from which we estimate it has a radial velocity that is near the center of the range exhibited by the C1-S massive core. A second outflow-driving source is also detected within the projected boundary of C1-S, but appears to be at a different radial velocity. After considering properties of the outflows, we conclude C1-Sa is a promising candidate for an early-stage massive protostar and as such it shows that these early phases of massive star formation can involve highly ordered outflow, and thus accretion, processes, similar to models developed to explain low-mass protostars.

preprint2016arXiv

Chemical differentiation in a prestellar core traces non-uniform illumination

Dense cloud cores present chemical differentiation due to the different distribution of C-bearing and N-bearing molecules, the latter being less affected by freeze-out onto dust grains. In this letter we show that two C-bearing molecules, CH$_3$OH and $c$-C$_3$H$_2$, present a strikingly different (complementary) morphology while showing the same kinematics toward the prestellar core L1544. After comparing their distribution with large scale H$_2$ column density N(H$_2$) map from the Herschel satellite, we find that these two molecules trace different environmental conditions in the surrounding of L1544: the $c$-C$_3$H$_2$ distribution peaks close to the southern part of the core, where the surrounding molecular cloud has a N(H$_2$) sharp edge, while CH$_3$OH mainly traces the northern part of the core, where N(H$_2$) presents a shallower tail. We conclude that this is evidence of chemical differentiation driven by different amount of illumination from the interstellar radiation field: in the South, photochemistry maintains more C atoms in the gas phase allowing carbon chain (such as $c$-C$_3$H$_2$) production; in the North, C is mainly locked in CO and methanol traces the zone where CO starts to freeze out significantly. During the process of cloud contraction, different gas and ice compositions are thus expected to mix toward the central regions of the core, where a potential Solar-type system will form. An alternative view on carbon-chain chemistry in star-forming regions is also provided.

preprint2016arXiv

Infall/Expansion Velocities in the Low-Mass Dense Cores L492, L694-2, and L1521F: Dependence on Position and Molecular Tracer

Although surveys of infall motions in dense cores have been carried out for years, few surveys have focused on mapping infall across cores using multiple spectral line observations. To fill this gap, we present IRAM 30-m Telescope maps of N2H+(1-0), DCO+(2-1), DCO+(3-2), and HCO+(3-2) emission towards two prestellar cores (L492 and L694-2) and one protostellar core (L1521F). We find that the measured infall velocity varies with position across each core and choice of molecular line, likely as a result of radial variations in core chemistry and dynamics. Line-of-sight infall speeds estimated from DCO+(2-1) line profiles can decrease by 40-50 m/s when observing at a radial offset >= 0.04 pc from the core's dust continuum emission peak. Median infall speeds calculated from all observed positions across a core can also vary by as much as 65 m/s depending on the transition. These results show that while single-pointing, single-transition surveys of core infall velocities may be good indicators of whether a core is either contracting or expanding, the magnitude of the velocities they measure are significantly impacted by the choice of molecular line, proximity to the core center, and core evolutionary state.

preprint2016arXiv

Protostellar Disk Formation Enabled by Removal of Small Dust Grains

It has been shown that a realistic level of magnetization of dense molecular cloud cores can suppress the formation of a rotationally supported disk (RSD) through catastrophic magnetic braking in the axisymmetric ideal MHD limit. In this study, we present conditions for the formation of RSDs through non-ideal MHD effects computed self-consistently from an equilibrium chemical network. We find that removing from the standard MRN distribution the large population of very small grains (VSGs) of ~10 $Å$ to few 100 $Å$ that dominate the coupling of the bulk neutral matter to the magnetic field increases the ambipolar diffusivity by ~1--2 orders of magnitude at densities below 10$^{10}$ cm$^{-3}$. The enhanced ambipolar diffusion (AD) in the envelope reduces the amount of magnetic flux dragged by the collapse into the circumstellar disk-forming region. Therefore, magnetic braking is weakened and more angular momentum can be retained. With continuous high angular momentum inflow, RSDs of tens of AU are able to form, survive, and even grow in size, depending on other parameters including cosmic-ray ionization rate, magnetic field strength, and rotation speed. Some disks become self-gravitating and evolve into rings in our 2D (axisymmetric) simulations, which have the potential to fragment into (close) multiple systems in 3D. We conclude that disk formation in magnetized cores is highly sensitive to chemistry, especially to grain sizes. A moderate grain coagulation/growth to remove the large population of VSGs, either in the prestellar phase or during free-fall collapse, can greatly promote AD and help formation of tens of AU RSDs.

preprint2016arXiv

Structure, Dynamics and Deuterium Fractionation of Massive Pre-Stellar Cores

High levels of deuterium fraction in N$_2$H$^+$ are observed in some pre-stellar cores. Single-zone chemical models find that the timescale required to reach observed values ($D_{\rm frac}^{{\rm N}_2{\rm H}^+} \equiv {\rm N}_2{\rm D}^+/{\rm N}_2{\rm H}^+ \gtrsim 0.1$) is longer than the free-fall time, possibly ten times longer. Here, we explore the deuteration of turbulent, magnetized cores with 3D magnetohydrodynamics simulations. We use an approximate chemical model to follow the growth in abundances of N$_2$H$^+$ and N$_2$D$^+$. We then examine the dynamics of the core using each tracer for comparison to observations. We find that the velocity dispersion of the core as traced by N$_2$D$^+$ appears slightly sub-virial compared to predictions of the Turbulent Core Model of McKee & Tan, except at late times just before the onset of protostar formation. By varying the initial mass surface density, the magnetic energy, the chemical age, and the ortho-to-para ratio of H$_2$, we also determine the physical and temporal properties required for high deuteration. We find that low initial ortho-to-para ratios ($\lesssim 0.01$) and/or multiple free-fall times ($\gtrsim 3$) of prior chemical evolution are necessary to reach the observed values of deuterium fraction in pre-stellar cores.

preprint2016arXiv

The Deuterium Fraction in Massive Starless Cores and Dynamical Implications

We study deuterium fractionation in two massive starless/early-stage cores C1-N and C1-S in Infrared Dark Cloud (IRDC) G028.37+00.07, first identified by Tan et al. (2013) with ALMA. Line emission from multiple transitions of $\rm N_2H^+$ and $\rm N_2D^+$ were observed with the ALMA, CARMA, SMA, JCMT, NRO 45m and IRAM 30m telescopes. By simultaneously fitting the spectra, we estimate the excitation conditions and deuterium fraction, $D_{\rm frac}^{\rm N_2H^+} \equiv [\rm N_2D^+]/[N_2H^+]$, with values of $D_{\rm frac}^{\rm N_2H^+} \simeq 0.2$--$0.7$, several orders of magnitude above the cosmic [D]/[H] ratio. Additional observations of o-H$_2$D$^+$ are also presented that help constrain the ortho-to-para ratio of $\rm H_2$, which is a key quantity affecting the degree of deuteration. We then present chemodynamical modeling of the two cores, exploring especially the implications for the collapse rate relative to free-fall, $α_{\rm ff}$. In order to reach the high level of observed deuteration of $\rm N_2H^+$, we find that the most likely evolutionary history of the cores involves collapse at a relatively slow rate, $\lesssim1/10$th of free-fall.

preprint2016arXiv

The Spatial Distribution of Complex Organic Molecules in the L1544 Pre-stellar Core

The detection of complex organic molecules (COMs) toward cold sources such as pre-stellar cores (with T<10 K), has challenged our understanding of the formation processes of COMs in the interstellar medium. Recent modelling on COM chemistry at low temperatures has provided new insight into these processes predicting that COM formation depends strongly on parameters such as visual extinction and the level of CO freeze out. We report deep observations of COMs toward two positions in the L1544 pre-stellar core: the dense, highly-extinguished continuum peak with Av>=30 mag within the inner 2700 au; and a low-density shell with average Av~7.5-8 mag located at 4000 au from the core's center and bright in CH3OH. Our observations show that CH3O, CH3OCH3 and CH3CHO are more abundant (by factors ~2-10) toward the low-density shell than toward the continuum peak. Other COMs such as CH3OCHO, c-C3H2O, HCCCHO, CH2CHCN and HCCNC show slight enhancements (by factors <=3) but the associated uncertainties are large. This suggests that COMs are actively formed and already present in the low-density shells of pre-stellar cores. The modelling of the chemistry of O-bearing COMs in L1544 indicates that these species are enhanced in this shell because i) CO starts freezing out onto dust grains driving an active surface chemistry; ii) the visual extinction is sufficiently high to prevent the UV photo-dissociation of COMs by the external interstellar radiation field; and iii) the density is still moderate to prevent severe depletion of COMs onto grains.

preprint2016arXiv

Theory of Gas Phase Scattering and Reactivity for Astrochemistry

Because of the very peculiar conditions of chemistry in many astrophysical gases (low densities, mostly low temperatures, kinetics-dominated chemical evolution), great efforts have been devoted to study molecular signatures and chemical evolution. While experiments are being performed in many laboratories, it appears that the efforts directed towards theoretical works are not as strong. This report deals with the present status of chemical physics/physical chemistry theory, for the qualitative and quantitative understanding of kinetics of molecular scattering, being it reactive or inelastic. By gathering several types of expertise, from applied mathematics to physical chemistry, dialog is made possible, as a step towards new and more adapted theoretical frameworks, capable of meeting the theoretical, methodological and numerical challenges of kinetics-dominated gas phase chemistry in astrophysical environments. A state of the art panorama is presented, alongside present-day strengths and shortcomings. However, coverage is not complete, being limited in this report to actual attendance of the workshop. Some paths towards relevant progress are proposed.

preprint2015arXiv

Detectability of deuterated water in prestellar cores

Water is an important molecule in the chemical and thermal balance of dense molecular gas, but knowing its history through-out the various stages of the star formation is a fundamental problem. Its molecular deuteration provides us with a crucial clue to its formation history. H$_2$O has recently been detected for the first time towards the prestellar core L1544 with the Herschel Space Observatory with a high spectral resolution (HIFI instrument). Prestellar cores provide the original reservoir of material from which future planetary systems are built, but few observational constraints exist on the formation of water and none on its deuteration before the collapse starts and a protostar forms at the centre. We report on new APEX observations of the ground state 1$_{0,1}$-0$_{0,0}$ HDO transition at 464 GHz towards the prestellar core L1544. The line is undetected, and we present an extensive study of the conditions for its detectability in cold and dense cloud cores. The water and deuterated water abundances have been estimated using an advanced chemical model simplified for the limited number of reactions or processes that are active in cold regions (< 15 K). We use the LIME radiative transfer code to compute the expected intensity and profile of both H$_2$O and HDO lines and compare them with the observations. We present several ad hoc profiles that best-fit the observations and compare the profiles with results from an astrochemical modelling, coupling gas phase and grain surface chemistry. Our comparison between observations, radiative transfer, and chemical modelling shows the limits of detectability for singly deuterated water, through the ground-state transitions 1$_{0,1}$-0$_{0,0}$ and 1$_{1,1}$-0$_{0,0}$ at 464.9 and 893.6 GHz, respectively, with both single-dish telescope and interferometric observations.

preprint2015arXiv

Interstellar dust charging in dense molecular clouds: cosmic ray effects

The local cosmic-ray (CR) spectra are calculated for typical characteristic regions of a cold dense molecular cloud, to investigate two so far neglected mechanisms of dust charging: collection of suprathermal CR electrons and protons by grains, and photoelectric emission from grains due to the UV radiation generated by CRs. The two mechanisms add to the conventional charging by ambient plasma, produced in the cloud by CRs. We show that the CR-induced photoemission can dramatically modify the charge distribution function for submicron grains. We demonstrate the importance of the obtained results for dust coagulation: While the charging by ambient plasma alone leads to a strong Coulomb repulsion between grains and inhibits their further coagulation, the combination with the photoemission provides optimum conditions for the growth of large dust aggregates in a certain region of the cloud, corresponding to the densities $n(\mathrm{H_2})$ between $\sim10^4$ cm$^{-3}$ and $\sim10^6$ cm$^{-3}$. The charging effect of CR is of generic nature, and therefore is expected to operate not only in dense molecular clouds but also in the upper layers and the outer parts of protoplanetary discs.

preprint2015arXiv

Mid-J CO observations of Perseus B1-East 5: evidence for turbulent dissipation via low-velocity shocks

Giant molecular clouds contain supersonic turbulence and magnetohydrodynamic simulations predict that this turbulence should decay rapidly. Such turbulent dissipation has the potential to create a warm (T ~100 K) gas component within a molecular cloud. We present observations of the CO J = 5-4 and 6-5 transitions, taken with the Herschel Space Observatory, towards the Perseus B1-East 5 region. We combine these new observations with archival measurements of lower rotational transitions and fit photodissociation region models to the data. We show that Perseus B1-E5 has an anomalously large CO J = 6-5 integrated intensity, consistent with a warm gas component existing within the region. This excess emission is consistent with predictions for shock heating due to the dissipation of turbulence in low velocity shocks with the shocks having a volume filling factor of 0.15 per cent. We find that B1-E has a turbulent energy dissipation rate of 3.5 x 10$^{32}$ erg / s and a dissipation time-scale that is only a factor of 3 larger than the flow crossing time-scale.

preprint2015arXiv

Next Generation Very Large Array Memo No. 6, Science Working Group 1: The Cradle of Life

This paper discusses compelling science cases for a future long-baseline interferometer operating at millimeter and centimeter wavelengths, like the proposed Next Generation Vary Large Array (ngVLA). We report on the activities of the Cradle of Life science working group, which focused on the formation of low- and high-mass stars, the formation of planets and evolution of protoplanetary disks, the physical and compositional study of Solar System bodies, and the possible detection of radio signals from extraterrestrial civilizations. We propose 19 scientific projects based on the current specification of the ngVLA. Five of them are highlighted as possible Key Science Projects: (1) Resolving the density structure and dynamics of the youngest HII regions and high-mass protostellar jets, (2) Unveiling binary/multiple protostars at higher resolution, (3) Mapping planet formation regions in nearby disks on scales down to 1 AU, (4) Studying the formation of complex molecules, and (5) Deep atmospheric mapping of giant planets in the Solar System. For each of these projects, we discuss the scientific importance and feasibility. The results presented here should be considered as the beginning of a more in-depth analysis of the science enabled by such a facility, and are by no means complete or exhaustive.

preprint2015arXiv

The Deuteration Clock for Massive Starless Cores

To understand massive star formation requires study of its initial conditions. Two massive starless core candidates, C1-N & C1-S, have been detected in IRDC G028.37+00.07 in $\rm N_2D^+$(3-2) with $ALMA$. From their line widths, either the cores are subvirial and are thus young structures on the verge of near free-fall collapse, or they are threaded by $\sim1$ mG $B$-fields that help support them in near virial equilibrium and potentially have older ages. We modeled the deuteration rate of $\rm N_2H^+$ to constrain collapse rates of the cores. First, to measure their current deuterium fraction, $D_{\rm frac}^{\rm N_2H^+}$ $\equiv [\rm N_2D^+]/[N_2H^+]$, we observed multiple transitions of $\rm N_2H^+$ and $\rm N_2D^+$ with $CARMA$, $SMA$, $JCMT$, $NRO~45m$ and $IRAM~30m$, to complement the $ALMA$ data. For both cores we derived $D_{\rm frac}^{\rm N_2H^+}\sim0.3$, several orders of magnitude above the cosmic [D]/[H] ratio. We then carried out chemodynamical modeling, exploring how collapse rate relative to free-fall, $α_{\rm ff}$, affects the level of $D_{\rm frac}^{\rm N_2H^+}$ that is achieved from a given initial condition. To reach the observed $D_{\rm frac}^{\rm N_2H^+}$, most models require slow collapse with $α_{\rm ff}\sim0.1$, i.e., $\sim1/10$th of free-fall. This makes it more likely that the cores have been able to reach a near virial equilibrium state and we predict that strong $B$-fields will eventually be detected. The methods developed here will be useful for measurement of the pre-stellar core mass function.

preprint2015arXiv

The Deuterium Fractionation Timescale in Dense Cloud Cores: A Parameter Space Exploration

The deuterium fraction [N$_2$D$^+$]/[N$_2$H$^+$], may provide information about the ages of dense, cold gas structures, important to compare with dynamical models of cloud core formation and evolution. Here we introduce a complete chemical network with species containing up to three atoms, with the exception of the Oxygen chemistry, where reactions involving H$_3$O$^+$ and its deuterated forms have been added, significantly improving the consistency with comprehensive chemical networks. Deuterium chemistry and spin states of H$_2$ and H$_3^+$ isotopologues are included in this primarily gas-phase chemical model. We investigate dependence of deuterium chemistry on model parameters: density ($n_{\rm H}$), temperature, cosmic ray ionization rate, and gas-phase depletion factor of heavy elements ($f_{\rm D}$). We also explore the effects of time-dependent freeze-out of gas-phase species and dynamical evolution of density at various rates relative to free-fall collapse. For a broad range of model parameters, the timescales to reach large values of $D_{\rm frac}^{\rm N_2H^+} \gtrsim 0.1$, observed in some low- and high-mass starless cores, are relatively long compared to the local free-fall timescale. These conclusions are unaffected by introducing time-dependent freeze-out and considering models with evolving density, unless the initial $f_{\rm D} \gtrsim$ 10. For fiducial model parameters, achieving $D_{\rm frac}^{\rm N_2H^+} \gtrsim 0.1$ requires collapse to be proceeding at rates at least several times slower than that of free-fall collapse, perhaps indicating a dynamically important role for magnetic fields in the support of starless cores and thus the regulation of star formation.

preprint2014arXiv

Chemistry and Radiative Transfer of Water in Cold, Dense Clouds

The Herschel Space Observatory's recent detections of water vapor in the cold, dense cloud L1544 allow a direct comparison between observations and chemical models for oxygen species in conditions just before star formation. We explain a chemical model for gas phase water, simplified for the limited number of reactions or processes that are active in extreme cold ($<$ 15 K). In this model, water is removed from the gas phase by freezing onto grains and by photodissociation. Water is formed as ice on the surface of dust grains from O and OH and released into the gas phase by photodesorption. The reactions are fast enough with respect to the slow dynamical evolution of L1544 that the gas phase water is in equilibrium for the local conditions thoughout the cloud. We explain the paradoxical radiative transfer of the H$_2$O ($1_{10}-1_{01}$) line. Despite discouragingly high optical depth caused by the large Einstein A coefficient, the subcritical excitation in the cold, rarefied H$_2$ causes the line brightness to scale linearly with column density. Thus the water line can provide information on the chemical and dynamical processes in the darkest region in the center of a cold, dense cloud. The inverse P-Cygni profile of the observed water line generally indicates a contracting cloud. This profile is reproduced with a dynamical model of slow contraction from unstable quasi-static hydrodynamic equilibrium (an unstable Bonnor-Ebert sphere).

preprint2014arXiv

Detectability of Glycine in Solar-type System Precursors

Glycine (NH2CH2COOH) is the simplest amino acid relevant for life. Its detection in the interstellar medium is key to understand the formation mechanisms of pre-biotic molecules and their subsequent delivery onto planetary systems. Glycine has extensively been searched for toward hot molecular cores, although these studies did not yield any firm detection. In contrast to hot cores, low-mass star forming regions, and in particular their earliest stages represented by cold pre-stellar cores, may be better suited for the detection of glycine as well as more relevant for the study of pre-biotic chemistry in young Solar System analogs. We present 1D spherically symmetric radiative transfer calculations of the glycine emission expected to arise from the low-mass pre-stellar core L1544. Water vapour has recently been reported toward this core, indicating that a small fraction of the grain mantles in L1544 (~0.5%) has been injected into the gas phase. Assuming that glycine is photo-desorbed together with water in L1544, and considering a solid abundance of glycine on ices of ~1E-4 with respect to water, our calculations reveal that several glycine lines between 67 GHz and 80 GHz have peak intensities larger than 10 mK. These results show for the first time that glycine could reach detectable levels in cold objects such as L1544. This opens up the possibility to detect glycine, and other pre-biotic species, at the coldest and earliest stages in the formation of Solar-type systems with near-future instrumentation such as the Band 2 receivers of ALMA.

preprint2014arXiv

Deuterated methanol in the pre-stellar core L1544

Using the IRAM 30m telescope, we mapped the methanol emission in the pre-stellar core L1544 and observed singly deuterated methanol (CH$_2$DOH and CH$_3$OD) towards the dust peak of L1544. Non-LTE radiative transfer modelling was performed on three CH$_3$OH emissions lines at 96.7 GHz, using a Bonnor-Ebert sphere as a model for the source. We have also assumed a centrally decreasing abundance profile to take the molecule freeze-out in the inner core into account. The column density of CH$_2$DOH was derived assuming LTE excitation and optically thin emission. The CH$_3$OH emission has a highly asymmetric morphology, resembling a non-uniform ring surrounding the dust peak, where CO is mainly frozen onto dust grains. The observations provide an accurate measure of methanol deuteration in the cold pre-stellar gas. The derived abundance ratio is [CH$_2$DOH]/[CH$_3$OH] $= 0.10\pm 0.03$, which is significantly smaller than the ones found in low-mass Class 0 protostars and smaller than the deuterium Fraction measured in other molecules towards L1544. The low deuterium fractionation observed in L1544 and the morphology of the CH$_3$OH emission suggest that we are mainly tracing the outer parts of the core, where CO just started to freeze-out onto dust grains.

preprint2014arXiv

Deuterium chemistry of dense gas in the vicinity of low-mass and massive star forming regions

The standard interstellar ratio of deuterium to hydrogen (D/H) atoms is $\sim 1.5 \times 10^{-5}$. However, the deuterium fractionation is in fact found to be enhanced, to different degrees, in cold, dark cores, hot cores around massive star forming regions, lukewarm cores, and warm cores ({\it hereafter}, hot corinos) around low-mass star forming regions. In this paper, we investigate the overall differences in the deuterium chemistry between hot cores and hot corinos. We have modelled the chemistry of dense gas around low-mass and massive star forming regions using a gas-grain chemical model. We investigate the influence of varying the core density, the depletion efficiency of gaseous species on to dust grains, the collapse mode and the final mass of the protostar on the chemical evolution of star forming regions. We find that the deuterium chemistry is, in general, most sensitive to variations of the depletion efficiency on to grain surfaces, in agreement with observations. In addition, the results showed that the chemistry is more sensitive to changes in the final density of the collapsing core in hot cores than in hot corinos. Finally, we find that ratios of deuterated sulphur bearing species in dense gas around hot cores and corinos may be good evolutionary indicators in a similar way as their non deuterated counterparts.

preprint2014arXiv

Deuterium Fractionation: the Ariadne's Thread from the Pre-collapse Phase to Meteorites and Comets today

The Solar System formed about 4.6 billion years ago from a condensation of matter inside a molecular cloud. Trying to reconstruct what happened is the goal of this chapter. For that, we put together our understanding of Galactic objects that will eventually form new suns and planetary systems, with our knowledge on comets, meteorites and small bodies of the Solar System today. Our specific tool is the molecular deuteration, namely the amount of deuterium with respect to hydrogen in molecules. This is the Ariadne's thread that helps us to find the way out from a labyrinth of possible histories of our Solar System. The chapter reviews the observations and theories of the deuterium fractionation in pre-stellar cores, protostars, protoplanetary disks, comets, interplanetary dust particles and meteorites and links them together trying to build up a coherent picture of the history of the Solar System formation. We emphasise the interdisciplinary nature of the chapter, which gathers together researchers from different communities with the common goal of understanding the Solar System history.

preprint2014arXiv

Herschel HIFI observations of O$_2$ toward Orion: special conditions for shock enhanced emission

We report observations of molecular oxygen (O$_2$) rotational transitions at 487 GHz, 774 GHz, and 1121 GHz toward Orion Peak A. The O2 lines at 487 GHz and 774 GHz are detected at velocities of 10-12 km/s with line widths 3 km/s; however, the transition at 1121 GHz is not detected. The observed line characteristics, combined with the results of earlier observations, suggest that the region responsible for the O$_2$ emission is 9" (6e16 cm) in size, and is located close to the H2 Peak 1position (where vibrationally-excited H$_2$ emission peaks), and not at Peak A, 23" away. The peak O2 column density is 1.1e18/cm2. The line velocity is close to that of 621 GHz water maser emission found in this portion of the Orion Molecular Cloud, and having a shock with velocity vector lying nearly in the plane of the sky is consistent with producing maximum maser gain along the line-of-sight. The enhanced O$_2$ abundance compared to that generally found in dense interstellar clouds can be explained by passage of a low-velocity C-shock through a clump with preshock density 2e4/cm3, if a reasonable flux of UV radiation is present. The postshock O$_2$ can explain the emission from the source if its line of sight dimension is ~10 times larger than its size on the plane of the sky. The special geometry and conditions required may explain why O$_2$ emission has not been detected in the cores of other massive star-forming molecular clouds.

preprint2014arXiv

Massive Star Formation

The enormous radiative and mechanical luminosities of massive stars impact a vast range of scales and processes, from the reionization of the universe, to the evolution of galaxies, to the regulation of the interstellar medium, to the formation of star clusters, and even to the formation of planets around stars in such clusters. Two main classes of massive star formation theory are under active study, Core Accretion and Competitive Accretion. In Core Accretion, the initial conditions are self-gravitating, centrally concentrated cores that condense with a range of masses from the surrounding, fragmenting clump environment. They then undergo relatively ordered collapse via a central disk to form a single star or a small-N multiple. In this case, the pre-stellar core mass function has a similar form to the stellar initial mass function. In Competitive Accretion, the material that forms a massive star is drawn more chaotically from a wider region of the clump without passing through a phase of being in a massive, coherent core. In this case, massive star formation must proceed hand in hand with star cluster formation. If stellar densities become very high near the cluster center, then collisions between stars may also help to form the most massive stars. We review recent theoretical and observational progress towards understanding massive star formation, considering physical and chemical processes, comparisons with low and intermediate-mass stars, and connections to star cluster formation.

preprint2014arXiv

The dynamical properties of dense filaments in the infrared dark cloud G035.39-00.33

Infrared Dark Clouds (IRDCs) are unique laboratories to study the initial conditions of high-mass star and star cluster formation. We present high-sensitivity and high-angular resolution IRAM PdBI observations of N2H+ (1-0) towards IRDC G035.39-00.33. It is found that G035.39-00.33 is a highly complex environment, consisting of several mildly supersonic filaments (sigma_NT/c_s ~1.5), separated in velocity by <1 km s^-1 . Where multiple spectral components are evident, moment analysis overestimates the non-thermal contribution to the line-width by a factor ~2. Large-scale velocity gradients evident in previous single-dish maps may be explained by the presence of substructure now evident in the interferometric maps. Whilst global velocity gradients are small (<0.7 km s^-1 pc^-1), there is evidence for dynamic processes on local scales (~1.5-2.5 km s^-1 pc^-1 ). Systematic trends in velocity gradient are observed towards several continuum peaks. This suggests that the kinematics are influenced by dense (and in some cases, starless) cores. These trends are interpreted as either infalling material, with accretion rates ~(7 \pm 4)x10^-5 M_sun yr^-1 , or expanding shells with momentum ~24 \pm 12 M_sun km s^-1 . These observations highlight the importance of high-sensitivity and high-spectral resolution data in disentangling the complex kinematic and physical structure of massive star forming regions.

preprint2014arXiv

The dynamics of collapsing cores and star formation

Low-mass stars are generally understood to form by the gravitational collapse of the dense molecular clouds known as starless cores. Continuum observations have not been able to distinguish among the several different hypotheses that describe the collapse because the predicted density distributions are the almost the same, as they are for all spherical self-gravitating clouds. However, the predicted contraction velocities are different enough that the models can be discriminated by comparing the velocities at large and small radii. This can be done by observing at least two different molecular line transitions that are excited at different densities. For example, the spectral lines of the H2O (110 - 101) and C18O (1-0) have critical densities for collisional de-excitation that differ by 5 orders of magnitude. We compare observations of these lines from the contracting starless core L1544 against the spectra predicted for several different hypothetical models of contraction including the Larson-Penston flow, the inside-out collapse of the singular isothermal sphere, the quasi-equilibrium contraction of an unstable Bonnor-Ebert sphere, and the non-equilibrium collapse of an over-dense Bonnor-Ebert sphere. Only the model of the unstable quasi-equilibrium Bonnor-Ebert sphere is able to produce the observed shapes of both spectral lines. This model allows us to interpret other observations of molecular lines in L1544 to find that the inward velocities seen in observations of CS(2-1) and N2H+ are located within the starless core itself, in particular in the region where the density profile follows an inverse square law. If these conclusions were to hold in the analysis of other starless cores, this would imply that the formation of hydrostatic clouds within the turbulent interstellar medium is not only possible but not exceptional and may be an evolutionary phase in low-mass star formation.

preprint2013arXiv

Correlating Infall with Deuterium Fractionation in Dense Cores

We present a survey of HCO+ (3-2) observations pointed towards dense cores with previous measurements of N(N2D+)/N(N2H+). Of the 26 cores in this survey, five show the spectroscopic signature of outward motion, nine exhibit neither inward nor outward motion, eleven appear to be infalling, and one is not detected. We compare the degree of deuterium fractionation with infall velocities calculated from the HCO+ spectra and find that those cores with [D]/[H] > 0.1 are more likely to have the signature of inward motions than cores with smaller [D]/[H] ratios. Infall motions are also much more common in cores with masses exceeding their thermal Jeans masses. The fastest infall velocity measured belongs to one of the two protostellar cores in our survey, L1521F, and the observed motions are typically on the order of the sound speed.

preprint2013arXiv

Deep observations of O2 toward a low-mass protostar with Herschel-HIFI

According to traditional gas-phase chemical models, O2 should be abundant in molecular clouds, but until recently, attempts to detect interstellar O2 line emission with ground- and space-based observatories have failed. Following the multi-line detections of O2 with low abundances in the Orion and rho Oph A molecular clouds with Herschel, it is important to investigate other environments, and we here quantify the O2 abundance near a solar-mass protostar. Observations of O2, at 487 GHz toward a deeply embedded low-mass Class 0 protostar, NGC 1333-IRAS 4A, are presented, using the HIFI instrument on the Herschel Space Observatory. Complementary data of the chemically related NO and CO molecules are obtained as well. The high spectral resolution data are analysed using radiative transfer models to infer column densities and abundances, and are tested directly against full gas-grain chemical models. The deep HIFI spectrum fails to show O2 at the velocity of the dense protostellar envelope, implying one of the lowest abundance upper limits of O2/H2 at <6x10^-9 (3 sigma). However, a tentative (4.5 sigma) detection of O2 is seen at the velocity of the surrounding NGC 1333 molecular cloud, shifted by 1 km/s relative to the protostar. For the protostellar envelope, pure gas-phase models and gas-grain chemical models require a long pre-collapse phase (~0.7-1x10^6 years), during which atomic and molecular oxygen are frozen out onto dust grains and fully converted to H2O, to avoid overproduction of O2 in the dense envelope. The same model also reproduces the limits on the chemically related NO molecule. The tentative detection of O2 in the surrounding cloud is consistent with a low-density PDR model with small changes in reaction rates. The low O2 abundance in the collapsing envelope around a low-mass protostar suggests that the gas and ice entering protoplanetary disks is very poor in O2.

preprint2013arXiv

Detection of (15)NNH+ in L1544: non-LTE modelling of dyazenilium hyperfine line emission and accurate (14)N/(15)N values

Samples of pristine Solar System material found in meteorites and interplanetary dust particles are highly enriched in (15)N. Conspicuous nitrogen isotopic anomalies have also been measured in comets, and the (14)N/(15)N abundance ratio of the Earth is itself larger than the recognised pre-solar value by almost a factor of two. Ion--molecules, low-temperature chemical reactions in the proto-solar nebula have been repeatedly indicated as responsible for these (15)N-enhancements. We have searched for (15)N variants of the N2H+ ion in L1544, a prototypical starless cloud core which is one of the best candidate sources for detection owing to its low central core temperature and high CO depletion. The goal is the evaluation of accurate and reliable (14)N/(15)N ratio values for this species in the interstellar gas. A deep integration of the (15)NNH+ (1-0) line at 90.4 GHz has been obtained with the IRAM 30 m telescope. Non-LTE radiative transfer modelling has been performed on the J=1-0 emissions of the parent and (15)N-containing dyazenilium ions, using a Bonnor--Ebert sphere as a model for the source. A high-quality fit of the N2H+ (1--0) hyperfine spectrum has allowed us to derive a revised value of the N2H+ column density in L1544. Analysis of the observed N(15)NH+ and (15)NNH+ spectra yielded an abundance ratio N[N(15)NH+]/N[(15)NNH+] = 1.1 +/- 0.3. The obtained (14)N/(15)N ratio is ~ 1000 +/- 200, suggestive of a sizeable (15)N depletion in this molecular ion. Such a result is not consistent with the prediction of present nitrogen chemical models. As chemical models predict large (15)N fractionation of N2H+, we suggest that (15)N(14)N, or (15)N in some other molecular form, is preferentially depleted onto dust grains.

preprint2013arXiv

The Dynamics of Massive Starless Cores with ALMA

How do stars that are more massive than the Sun form, and thus how is the stellar initial mass function (IMF) established? Such intermediate- and high-mass stars may be born from relatively massive pre-stellar gas cores, which are more massive than the thermal Jeans mass. The Turbulent Core Accretion model invokes such cores as being in approximate virial equilibrium and in approximate pressure equilibrium with their surrounding clump medium. Their internal pressure is provided by a combination of turbulence and magnetic fields. Alternatively, the Competitive Accretion model requires strongly sub-virial initial conditions that then lead to extensive fragmentation to the thermal Jeans scale, with intermediate- and high-mass stars later forming by competitive Bondi-Hoyle accretion. To test these models, we have identified four prime examples of massive (~100Msun) clumps from mid-infrared extinction mapping of infrared dark clouds (IRDCs). Fontani et al. found high deuteration fractions of N2H+ in these objects, which are consistent with them being starless. Here we present ALMA observations of these four clumps that probe the N2D+(3-2) line at 2.3" resolution. We find six N2D+ cores and determine their dynamical state. Their observed velocity dispersions and sizes are broadly consistent with the predictions of the Turbulent Core model of self-gravitating, magnetized (with Alfven Mach number m_A~1) and virialized cores that are bounded by the high pressures of their surrounding clumps. However, in the most massive cores, with masses up to ~60Msun, our results suggest that moderately enhanced magnetic fields (so that m_A~0.3) may be needed for the structures to be in virial and pressure equilibrium. Magnetically regulated core formation may thus be important in controlling the formation of massive cores, inhibiting their fragmentation, and thus helping to establish the stellar IMF.

preprint2012arXiv

A Virialized Filamentary Infrared Dark Cloud

The initial conditions of massive star and star cluster formation are expected to be cold, dense and high column density regions of the interstellar medium, which can reveal themselves via near, mid and even far-infrared absorption as Infrared Dark Clouds (IRDCs). Elucidating the dynamical state of IRDCs thus constrains theoretical models of these complex processes. In particular, it is important to assess whether IRDCs have reached virial equilibrium, where the internal pressure balances that due to the self-gravitating weight of the cloud plus the pressure of the external environmental. We study this question for the filamentary IRDC G035.39-00.33 by deriving mass from combined NIR & MIR extinction maps and velocity dispersion from C18O (1-0) & (2-1) line emission. In contrast to our previous moderately super-virial results based on 13CO emission and MIR-only extinction mapping, with improved mass measurements we now find that the filament is consistent with being in virial equilibrium, at least in its central parsec-wide region where ~1000 M_Sun snakes along several parsecs. This equilibrium state does not require large-scale net support or confinement by magnetic fields.

preprint2012arXiv

ALMA reveals a chemically evolved submillimeter galaxy at z=4.76

The chemical properties of high-z galaxies provide important information to constrain galaxy evolutionary scenarios. However, widely-used metallicity diagnostics based on rest-frame optical emission lines are not usable for heavily dust-enshrouded galaxies (such as Sub-Millimeter Galaxies; SMGs), especially at z>3. Here we focus on the flux ratio of the far-infrared fine-structure emission lines [NII]205um and [CII]158um to assess the metallicity of high-z SMGs. Through ALMA cycle 0 observations, we have detected the [NII]205um emission in a strongly [CII]-emitting SMG, LESS J033229.4-275619 at z=4.76. The velocity-integrated [NII]/[CII] flux ratio is 0.043 +/- 0.008. This is the first measurement of the [NII]/[CII] flux ratio in high-z galaxies, and the inferred flux ratio is similar to the ratio observed in the nearby universe (~0.02-0.07). The velocity-integrated flux ratio and photoionization models suggest that the metallicity in this SMG is consistent with solar, implying the chemical evolution has progressed very rapidly in this system at z=4.76. We also obtain a tight upper limit on the CO(12-11) transition, which translates into CO(12-11)/CO(2-1) <3.8 (3 sigma). This suggests that the molecular gas clouds in LESS J033229.4-275619 are not affected significantly by the radiation field emitted by the AGN in this system.

preprint2012arXiv

Complex, Quiescent Kinematics in a Highly Filamentary Infrared Dark Cloud

Infrared Dark Clouds (IRDCs) host the initial conditions under which massive stars and stellar clusters form. We have obtained high sensitivity and high spectral resolution observations with the IRAM 30m antenna, which allowed us to perform detailed analysis of the kinematics within one IRDC, G035.39-00.33. We focus on the 1-0 and 3-2 transitions of N2H+, C18O (1-0), and make comparison with SiO (2-1) observations and extinction mapping. Three interacting filaments of gas are found. We report large-scale velocity coherence throughout the cloud, evidenced through small velocity gradients and relatively narrow line widths. This suggests that the merging of these filaments is somewhat "gentle", possibly regulated by magnetic fields. This merging of filaments may be responsible for the weak parsec-scale SiO emission detected by Jimenez-Serra et al. 2010, via grain mantle vaporization. A systematic velocity shift between the N2H+ (1-0) and C18O (1-0) gas throughout the cloud of 0.18 +/- 0.04 kms^{-1} is also found, consistent with a scenario of collisions between filaments which is still ongoing. The N2H+ (1-0) is extended throughout the IRDC and it does not only trace dense cores, as found in nearby low-mass star-forming regions. The average H2 number density across the IRDC is ~ 5 x 10^4 cm^{-3}, at least one order of magnitude larger than in nearby molecular clouds where low-mass stars are forming. A temperature gradient perpendicular to the filament is found. From our study, we conclude that G035.39-00.33 (clearly seen in the extinction map and in N2H+) has been formed via the collision between two relatively quiescent filaments with average densities of ~ 5 x 10^3 cm^{-3}, moving with relative velocities of ~ 5 kms^{-1}. The accumulation of material at the merging points started > 1 Myr ago and it is still ongoing.

preprint2012arXiv

First detection of water vapor in a pre-stellar core

Water is a crucial molecule in molecular astrophysics as it controls much of the gas/grain chemistry, including the formation and evolution of more complex organic molecules in ices. Pre-stellar cores provide the original reservoir of material from which future planetary systems are built, but few observational constraints exist on the formation of water and its partitioning between gas and ice in the densest cores. Thanks to the high sensitivity of the Herschel Space Observatory, we report on the first detection of water vapor at high spectral resolution toward a dense cloud on the verge of star formation, the pre-stellar core L1544. The line shows an inverse P-Cygni profile, characteristic of gravitational contraction. To reproduce the observations, water vapor has to be present in the cold and dense central few thousand AU of L1544, where species heavier than Helium are expected to freeze-out onto dust grains, and the ortho:para H2 ratio has to be around 1:1 or larger. The observed amount of water vapor within the core (about 1.5x10^{-6} Msun) can be maintained by Far-UV photons locally produced by the impact of galactic cosmic rays with H2 molecules. Such FUV photons irradiate the icy mantles, liberating water wapor in the core center. Our Herschel data, combined with radiative transfer and chemical/dynamical models, shed light on the interplay between gas and solids in dense interstellar clouds and provide the first measurement of the water vapor abundance profile across the parent cloud of a future solar-type star and its potential planetary system.

preprint2012arXiv

H2D+ in the high mass star-forming Region Cygnus-X

H2D+ is a primary ion which dominates the gas-phase chemistry of cold dense gas. Therefore it is hailed as a unique tool in probing the earliest, prestellar phase of star formation. Observationally, its abundance and distribution is however just beginning to be understood in low-mass prestellar and cluster-forming cores. In high mass star forming regions, H2D+ has been detected only in two cores, and its spatial distribution remains unknown. Here we present the first map of the 372 GHz ortho-H2D+ and N2H+ 4-3 transition in the DR21 filament of Cygnus-X with the JCMT, and N2D+ 3--2 and dust continuum with the SMA. We have discovered five very extended (<= 34000 AU diameter) weak structures in H2D+ in the vicinity of, but distinctly offset from embedded protostars. More surprisingly, the H2D+ peak is not associated with either a dust continuum or N2D+ peak. We have therefore uncovered extended massive cold dense gas that was undetected with previous molecular line and dust continuum surveys of the region. This work also shows that our picture of the structure of cores is too simplistic for cluster forming cores and needs to be refined: neither dust continuum with existing capabilities, nor emission in tracers like N2D+ can provide a complete census of the total prestellar gas in such regions. Sensitive H2D+ mapping of the entire DR21 filament is likely to discover more of such cold quiescent gas reservoirs in an otherwise active high mass star-forming region.

preprint2012arXiv

Herschel-PACS observations of [OI]63um towards submillimetre galaxies at z~1

We present Herschel-PACS spectroscopy of the [OI]63um far-infrared cooling line from a sample of six unlensed and spectroscopically-confirmed 870um-selected submillimetre (submm) galaxies (SMGs) at 1.1<z<1.6 from the LABOCA Extended Chandra Deep Field South (ECDFS) Submm Survey (LESS). This is the first survey of [OI]63um, one of the main photodissociation region (PDR) cooling lines, in SMGs. New high-resolution ALMA interferometric 870um continuum imaging confirms that these six Herschel-targeted SMG counterparts are bona fide sources of submm emission. We detect [OI]63um in two SMGs with a SNR >3, tentatively detect [OI]63um in one SMG, and constrain the line flux for the non-detections. We also exploit the combination of submm continuum photometry from 250-870um and our new PACS continuum measurements to constrain the far-infrared (FIR) luminosity, L_FIR, in these SMGs to < 30%. We find that SMGs do not show a deficit in their [OI]63um-to-far-infrared continuum luminosity ratios (with ratios ranging from ~0.5-1.5%), similar to what was seen previously for the [CII]158um-to-FIR ratios in SMGs. These observed ratios are about an order of magnitude higher than what is seen typically for local ultra luminous infrared galaxies (ULIRGs), which adds to the growing body of evidence that SMGs are not simply `scaled up' versions of local ULIRGs. Rather, the PDR line-to-L_FIR ratios suggest that the star formation modes of SMGs are likely more akin to that of local normal (lower-luminosity) star-forming galaxies, with the bulk of the star formation occurring in extended regions, galaxy-scale (~kpc) in size. These observations represent the first step towards a census of the major PDR cooling lines in typical SMGs that will be attainable with ALMA, enabling detailed modelling to probe the global properties of the star formation and the evolutionary status of SMGs.

preprint2012arXiv

Our astrochemical heritage

Our Sun and planetary system were born about 4.5 billion years ago. How did this happen and what is our heritage from these early times? This review tries to address these questions from an astrochemical point of view. On the one hand, we have some crucial information from meteorites, comets and other small bodies of the Solar System. On the other hand, we have the results of studies on the formation process of Sun-like stars in our Galaxy. These results tell us that Sun-like stars form in dense regions of molecular clouds and that three major steps are involved before the planet formation period. They are represented by the pre-stellar core, protostellar envelope and protoplanetary disk phases. Simultaneously with the evolution from one phase to the other, the chemical composition gains increasing complexity. In this review, we first present the information on the chemical composition of meteorites, comets and other small bodies of the Solar System, which is potentially linked to the first phases of the Solar System's formation. Then we describe the observed chemical composition in the pre-stellar core, protostellar envelope and protoplanetary disk phases, including the processes that lead to them. Finally, we draw together pieces from the different objects and phases to understand whether and how much we inherited chemically from the time of the Sun's birth.

preprint2011arXiv

Enhanced [CII] emission in a z=4.76 submillimetre galaxy

We present the detection of bright [CII] emission in the z=4.76 submillimetre galaxy LESS J033229.4-275619 using the Atacama Pathfinder EXperiment. This represents the highest redshift [CII] detection in a submm selected, star-formation dominated system. The AGN contributions to the [CII] and far-infrared (FIR) luminosities are small. We find an atomic mass derived from [CII] comparable to the molecular mass derived from CO. The ratio of the [CII], CO and FIR luminosities imply a radiation field strength G_0~10^3 and a density ~10^4 cm^-3 in a kpc-scale starburst, as seen in local and high redshift starbursts. The high L_[CII]/L_FIR=2.4x10^-3 and the very high L_[CII]/L_CO(1-0) ~ 10^4 are reminiscent of low metallicity dwarf galaxies, suggesting that the highest redshift star-forming galaxies may also be characterised by lower metallicities. We discuss the implications of a reduced metallicity on studies of the gas reservoirs, and conclude that especially at very high redshift, [CII] may be a more powerful and reliable tracer of the interstellar matter than CO.

preprint2011arXiv

EVLA Observations of the Barnard 5 Star-Forming Core: Embedded Filaments Revealed

We present a ~6.5'x8' Expanded Very Large Array (EVLA) mosaic observations of the NH3 (1,1) emission in the Barnard 5 region in Perseus, with an angular resolution of 6". This map covers the coherent region, where the dense gas presents subsonic non-thermal motions (as seen from single dish observations with the Green Bank Telescope, GBT). The combined EVLA and GBT observations reveal, for the first time, a striking filamentary structure (20" wide or 5,000 AU at the distance of Perseus) in this low-mass star forming region. The integrated intensity profile of this structure is consistent with models of an isothermal filament in hydrostatic equilibrium. The observed separation between the B5-IRS1 young stellar object (YSO), in the central region of the core, and the northern starless condensation matches the Jeans length of the dense gas. This suggests that the dense gas in the coherent region is fragmenting. The region observed displays a narrow velocity dispersion, where most of the gas shows evidence for subsonic turbulence, and where little spatial variations are present. It is only close to the YSO where an increase in the velocity dispersion is found, but still displaying subsonic non-thermal motions

preprint2011arXiv

Herschel Measurements of Molecular Oxygen in Orion

We report observations of three rotational transitions of molecular oxygen (O2) in emission from the H2 Peak 1 position of vibrationally excited molecular hydrogen in Orion. We observed the 487 GHz, 774 GHz, and 1121 GHz lines using HIFI on the Herschel Space Observatory, having velocities of 11 km s-1 to 12 km s-1 and widths of 3 km s-1. The beam-averaged column density is N(O2) = 6.5\times1016 cm-2, and assuming that the source has an equal beam filling factor for all transitions (beam widths 44, 28, and 19"), the relative line intensities imply a kinetic temperature between 65 K and 120 K. The fractional abundance of O2 relative to H2 is 0.3 - 7.3\times10-6. The unusual velocity suggests an association with a ~ 5" diameter source, denoted Peak A, the Western Clump, or MF4. The mass of this source is ~ 10 M\odot and the dust temperature is \geq 150 K. Our preferred explanation of the enhanced O2 abundance is that dust grains in this region are sufficiently warm (T \geq 100 K) to desorb water ice and thus keep a significant fraction of elemental oxygen in the gas phase, with a significant fraction as O2. For this small source, the line ratios require a temperature \geq 180 K. The inferred O2 column density \simeq 5\times1018 cm-2 can be produced in Peak A, having N(H2) \simeq 4\times1024 cm-2. An alternative mechanism is a low-velocity (10 to 15 km s-1) C-shock, which can produce N(O2) up to 1017 cm-2.

preprint2011arXiv

Initial conditions for star formation in clusters: physical and kinematical structure of the starless core Oph A-N6

We present high spatial (<300 AU) and spectral (0.07 km/s) resolution Submillimeter Array observations of the dense starless cluster core Oph A-N6, in the 1 mm dust continuum and the 3-2 line of N2H+ and N2D+. The dust continuum observations reveal a compact source not seen in single-dish observations, of size ~1000 AU and mass 0.005-0.01 M\odot. The combined line and single-dish observations reveal a core of size 3000 \times 1400 AU elongated in a NW-SE direction, with almost no variation in either line width or line center velocity across the map, and very small non-thermal motions. The deuterium fraction has a peak value of ~0.15 and is >0.05 over much of the core. The N2H+ column density profile across the major axis of Oph A-N6 is well represented by an isothermal cylinder, with temperature 20 K, peak density 7.1 \times 10^6 cm^{-3}, and N2H+ abundance 2.7 \times 10^{-10}. The mass of Oph A-N6 is estimated to be 0.29 M\odot, compared to a value of 0.18 M\odot from the isothermal cylinder analysis, and 0.63 M\odot for the critical mass for fragmentation of an isothermal cylinder. Compared to isolated low-mass cores, Oph A-N6 shows similar narrow line widths and small velocity variation, with a deuterium fraction similar to "evolved" dense cores. It is significantly smaller than isolated cores, with larger peak column and volume density. The available evidence suggests Oph A-N6 has formed through the fragmentation of the Oph A filament and is the precursor to a low-mass star. The dust continuum emission suggests it may already have begun to form a star.

preprint2011arXiv

Mapping Large-Scale CO Depletion in a Filamentary Infrared Dark Cloud

Infrared Dark Clouds (IRDCs) are cold, high mass surface density and high density structures, likely to be representative of the initial conditions for massive star and star cluster formation. CO emission from IRDCs has the potential to be useful for tracing their dynamics, but may be affected by depleted gas phase abundances due to freeze-out onto dust grains. Here we analyze C18O J=1-0 and J=2-1 emission line data, taken with the IRAM 30m telescope, of the highly filamentary IRDC G035.39.-0033. We derive the excitation temperature as a function of position and velocity, with typical values of ~7K, and thus derive total mass surface densities, Sigma_C18O, assuming standard gas phase abundances and accounting for optical depth in the line, which can reach values of ~1. The mass surface densities reach values of ~0.07 g/cm^2. We compare these results to the mass surface densities derived from mid-infrared (MIR) extinction mapping, Sigma_SMF, by Butler & Tan, which are expected to be insensitive to the dust temperatures in the cloud. With a significance of >10sigma, we find Sigma_C18O/Sigma_SMF decreases by about a factor of 5 as Sigma increases from ~0.02 to ~0.2 g/cm^2, which we interpret as evidence for CO depletion. Several hundred solar masses are being affected, making this one of the most massive clouds in which CO depletion has been observed directly. We present a map of the depletion factor in the filament and discuss implications for the formation of the IRDC.

preprint2011arXiv

Multi-line spectral imaging of dense cores in the Lupus molecular cloud

The molecular clouds Lupus 1, 3 and 4 were mapped with the Mopra telescope at 3 and 12 mm. Emission lines from high density molecular tracers were detected, i.e. NH$_3$ (1,1), NH$_3$ (2,2), N$_2$H$^+$ (1-0), HC$_3$N (3-2), HC$_3$N (10-9), CS (2-1), CH$_3$OH (2$_0-1_0$)A$^+$ and CH$_3$OH (2$_{-1}-1_{-1}$)E. Velocity gradients of more than 1 km s$^{-1}$ are present in Lupus 1 and 3 and multiple gas components are present in these clouds along some lines of sight. Lupus 1 is the cloud richest in high density cores, 8 cores were detected in it, 5 cores were detected in Lupus 3 and only 2 in Lupus 4. The intensity of the three species HC$_3$N, NH$_3$ and N$_2$H$^+$ changes significantly in the various cores: cores that are brighter in HC$_3$N are fainter or undetected in NH$_3$ and N$_2$H$^+$ and vice versa. We found that the column density ratios HC$_3$N/N$_2$H$^+$ and HC$_3$N/NH$_3$ change by one order of magnitude between the cores, indicating that also the chemical abundance of these species is different. The time dependent chemical code that we used to model our cores shows that the HC$_3$N/N$_2$H$^+$ and HC$_3$N/NH$_3$ ratios decrease with time therefore the observed column density of these species can be used as an indicator of the chemical evolution of dense cores. On this base we classified 5 out of 8 cores in Lupus 1 and 1 out of 5 cores in Lupus 3 as very young protostars or prestellar cores. Comparing the millimetre cores population with the population of the more evolved young stellar objects identified in the Spitzer surveys, we conclude that in Lupus 3 the bulk of the star formation activity has already passed and only a moderate number of stars are still forming. On the contrary, in Lupus 1 star formation is on-going and several dense cores are still in the pre--/proto--stellar phase. Lupus 4 is at an intermediate stage, with a smaller number of individual objects.

preprint2011arXiv

Observational studies of pre-stellar cores and infrared dark clouds

Stars like our Sun and planets like our Earth form in dense regions within interstellar molecular clouds, called pre-stellar cores (PSCs). PSCs provide the initial conditions in the process of star and planet formation. In the past 15 years, detailed observations of (low-mass) PSCs in nearby molecular cloud complexes have allowed us to find that they are cold (T < 10 K) and quiescent (molecular line widths are close to thermal), with a chemistry profoundly affected by molecular freeze-out onto dust grains. In these conditions, deuterated molecules flourish, becoming the best tools to unveil the PSC physical and chemical structure. Despite their apparent simplicity, PSCs still offer puzzles to solve and they are far from being completely understood. For example, what is happening to the gas and dust in their nuclei (the future stellar cradles) is still a mystery that awaits for ALMA. Other important questions are: how do different environments and external conditions affect the PSC physical/chemical structure? Are PSCs in high-mass star forming regions similar to the well-known low-mass PSCs? Here I review observational and theoretical work on PSCs in nearby molecular cloud complexes and the ongoing search and study of massive PSCs embedded in infrared dark clouds (IRDCs), which host the initial conditions for stellar cluster and high-mass star formation.

preprint2011arXiv

The enigmatic core L1451-mm: a first hydrostatic core? or a hidden VeLLO?

We present the detection of a dust continuum source at 3-mm (CARMA) and 1.3-mm (SMA), and 12CO(2-1) emission (SMA) towards the L1451-mm dense core. These detections suggest a compact object and an outflow where no point source at mid-infrared wavelengths is detected using Spitzer. An upper limit for the dense core bolometric luminosity of 0.05 Lsun is obtained. By modeling the broadband SED and the continuum interferometric visibilities simultaneously, we confirm that a central source of heating is needed to explain the observations. This modeling also shows that the data can be well fitted by a dense core with a YSO and disk, or by a dense core with a central First Hydrostatic Core (FHSC). Unfortunately, we are not able to decide between these two models, which produce similar fits. We also detect 12CO(2-1) emission with red- and blue-shifted emission suggesting the presence of a slow and poorly collimated outflow, in opposition to what is usually found towards young stellar objects but in agreement with prediction from simulations of a FHSC. This presents the best candidate, so far, for a FHSC, an object that has been identified in simulations of collapsing dense cores. Whatever the true nature of the central object in L1451-mm, this core presents an excellent laboratory to study the earliest phases of low-mass star formation.

preprint2010arXiv

Detection of N15NH+ in L1544

Excess levels of 15N isotopes which have been detected in primitive solar system materials are explained as a remnant of interstellar chemistry which took place in regions of the protosolar nebula. Chemical models of nitrogen fractionation in cold clouds predict an enhancement in the gas-phase abundance of 15N-bearing molecules, thus we have searched for 15N variants of the N2H+ ion in L1544, which is one of the best candidate sources for detection owing to its low central core temperature and high CO depletion. With the IRAM 30m telescope we have obtained deep integrations of the N2H+(1-0) line at 91.2 GHz. The N2H+(1-0) line has been detected toward the dust emission peak of L1544. The 14N/15N abundance ratio in N2H+ resulted 446+/-71, very close to the protosolar value of ~450, higher than the terrestrial ratio of ~270, and significantly lower than the lower limit in L1544 found by Gerin et al. (2009, ApJ, 570, L101) in the same object using ammonia isotopologues.

preprint2010arXiv

Direct observation of a sharp transition to coherence in Dense Cores

We present NH3 observations of the B5 region in Perseus obtained with the Green Bank Telescope (GBT). The map covers a region large enough (~11'x14') that it contains the entire dense core observed in previous dust continuum surveys. The dense gas traced by NH3(1,1) covers a much larger area than the dust continuum features found in bolometer observations. The velocity dispersion in the central region of the core is small, presenting subsonic non-thermal motions which are independent of scale. However, it is thanks to the coverage and high sensitivity of the observations that we present the detection, for the first time, of the transition between the coherent core and the dense but more turbulent gas surrounding it. This transition is sharp, increasing the velocity dispersion by a factor of 2 in less than 0.04 pc (the 31" beam size at the distance of Perseus, ~250 pc). The change in velocity dispersion at the transition is ~3 km/s/pc. The existence of the transition provides a natural definition of dense core: the region with nearly-constant subsonic non-thermal velocity dispersion. From the analysis presented here we can not confirm nor rule out a corresponding sharp density transition.

preprint2009arXiv

Dense Cores in Perseus: The Influence of Stellar Content and Cluster Environment

We present the chemistry, temperature, and dynamical state of a sample of 193 dense cores or core candidates in the Perseus Molecular cloud and compare the properties of cores associated with young stars and clusters with those which are not. The combination of our NH3 and CCS observations with previous millimeter, sub-millimeter, and Spitzer data available for this cloud enable us both to determine core properties precisely and to accurately classify cores as starless or protostellar. The properties of cores in different cluster environments and before-and-after star formation provide important constraints on simulations of star-formation, particularly under the paradigm that the essence of star formation is set by the turbulent formation of prestellar cores. We separate the influence of stellar content from that of cluster environment and find that cores within clusters have (1) higher kinetic temperatures and (2) lower fractional abundances of CCS and NH3. Cores associated with protostars have (1) slightly higher kinetic temperatures (2) higher NH3 excitation temperatures), (3) are at higher column density, have (4) slightly more non-thermal/turbulent NH3 linewidths, have (5) higher masses and have (6) lower fractional abundance of CCS. We find that neither cluster environment nor protostellar content makes a significant difference to the dynamical state of cores as estimated by the virial parameter -- most cores in each category are gravitationally bound. Overall, cluster environment and protostellar content have a smaller influence on the properties of the cores than is typically assumed, and the variation within categories is larger than the differences between categories.

preprint2009arXiv

Dynamics and Depletion in Thermally Supercritical Starless Cores

In previous studies we identified two classes of starless cores, thermally subcritical and supercritical, distinguished by different dynamical behavior and internal structure. Here we study the evolution of the dynamically-unstable, thermally-supercritical cores by means of a numerical hydrodynamic simulation that includes radiative equilibrium and simple molecular chemistry. We use our non-LTE radiative transfer code MOLLIE to predict observable CO and N2H+ line spectra, including the non-LTE hyperfine ratios of N2H+, during the contraction. These are compared against observations of the starless core L1544.

preprint2009arXiv

The Dust Emissivity Spectral Index in the Starless Core TMC-1C

In this paper we present a dust emission map of the starless core TMC-1C taken at 2100 microns. Along with maps at 160, 450, 850 and 1200 microns, we study the dust emissivity spectral index from the (sub)millimeter spectral energy distribution, and find that it is close to the typically assumed value of beta = 2. We also map the dust temperature and column density in TMC-1C, and find that at the position of the dust peak (A_V ~ 50), the line-of-sight-averaged temperature is ~7 K. Employing simple Monte Carlo modeling, we show that the data are consistent with a constant value for the emissivity spectral index over the whole map of TMC-1C.

preprint2002arXiv

Dense Cores in Dark Clouds. XIV. N2H+(1-0) maps of dense cloud cores

We present results of an extensive mapping survey of N2H+(1-0) in about 60 low mass cloud cores already mapped in the NH3(1,1) inversion transition line. The survey has been carried out at the FCRAO antenna with an angular resolution about 1.5 times finer than the previous ammonia observations. Cores with stars typically have map sizes about a factor of two smaller for N2H+ than for NH3, indicating the presence of denser and more centrally concentrated gas compared to starless cores. Significant correlations are found between NH3 and N2H+ column densities and excitation temperatures in starless cores, but not in cores with stars, suggesting a different chemical evolution of the two species. Velocity gradients range between 0.5 and 6 km/s/pc, similar to what has been found with NH3 data. ``Local'' velocity gradients show significant variation in both magnitude and direction, suggesting the presence of complexmotions not interpretable as simple solid body rotation. Integrated intensity profiles of starless cores present a ``central flattening'' and are consistent with a spherically symmetric density law n ~ r^{-1.2} for r < ~0.03 pc and n ~ r^{-2} at larger r. Cores with stars are better modelled with single density power laws with n ~ r^{-2}. Line widths change across the core but we did not find a general trend. The deviation in line width correlates with the mean line width, suggesting that the line of sight contains ~ 10 coherence lengths. The corresponding value of the coherence length, ~ 0.01 pc, is similar to the expected cutoff wavelength for MHD waves. This similarity may account for the increased ``coherence'' of line widths on small scales. Despite of the finer angular resolution, the majority of N2H+ and NH3 maps show a similar ``simple'' structure, with single peaks and no elongation.

Paola Caselli

What is connected

Connect this record

See the researcher in context

Building this map preview

79 published item(s)

Infall of material onto the filaments in Barnard 5

A train of shocks at 3000 au scale? Exploring the clash of an expanding bubble into the NGC 1333 IRAS 4 region. SOLIS XIV

An Interferometric View of H-MM1. I. Direct Observation of NH3 Depletion

Astrochemical modelling of infrared dark clouds

CH$_3$CN deuteration in the SVS13-A Class I hot-corino. SOLIS XV

FAUST III. Misaligned rotations of the envelope, outflow, and disks in the multiple protostellar system of VLA 1623$-$2417

First detection of CHD$_{2}$OH towards pre-stellar cores

Nitrogen fractionation towards a pre-stellar core traces isotope-selective photodissociation

Seeds of Life in Space (SOLIS). XIII. Nitrogen fractionation towards the protocluster OMC-2 FIR4

The Central 1000 au of a Pre-stellar Core Revealed with ALMA. II. Almost Complete Freeze-out

Tracking the ice mantle history in the Solar-type Protostars of NGC 1333 IRAS 4

Deuterium Chemodynamics of Massive Pre-Stellar Cores

Dissecting the super-critical filaments embedded in the 0.5 pc subsonic region of Barnard 5

Methanol Mapping in Cold Cores: Testing Model Predictions

A new proxy to estimate the cosmic-ray ionisation rate in dense cores

A protostellar system fed by a streamer of 10,500 au length

Carbon Chain Chemistry in Hot-Core Regions around Three Massive Young Stellar Objects Associated with 6.7 GHz Methanol Masers

Efficient methanol production on the dark side of a prestellar core

Hall Effect in Protostellar Disc Formation and Evolution

Hot Corinos Chemical Diversity: Myth or Reality?

Impact of magneto-rotational instability on grain growth in protoplanetary disks: I. Relevant turbulence properties

Orbital and mass constraints of the young binary system IRAS 16293-2422 A

Probabilistic detection of spectral line components

Propargylimine in the laboratory and in space: millimetre-wave spectroscopy and first detection in the ISM

Relative Alignment between Dense Molecular Cores and Ambient Magnetic Field: The Synergy of Numerical Models and Observations

Seeds of Life in Space (SOLIS). VI. Chemical evolution of sulfuretted species along the outflows driven by the low-mass protostellar binary NGC1333-IRAS4A

Submillimetre and far-infrared spectroscopy of monodeuterated amidogen radical (NHD): improved rest-frequencies for astrophysical observations

Ubiquitous $\rm NH_3$ supersonic component in L1688 coherent cores

Velocity-coherent Filaments in NGC 1333: Evidence for Accretion Flow?

Efficient Production of S$_8$ in Interstellar Ices: The effects of cosmic ray-driven radiation chemistry and non-diffusive bulk reactions

The first steps of Interstellar Phosphorus Chemistry

An Ordered Bipolar Outflow from a Massive Early-Stage Core

Chemical differentiation in a prestellar core traces non-uniform illumination

Infall/Expansion Velocities in the Low-Mass Dense Cores L492, L694-2, and L1521F: Dependence on Position and Molecular Tracer

Protostellar Disk Formation Enabled by Removal of Small Dust Grains

Structure, Dynamics and Deuterium Fractionation of Massive Pre-Stellar Cores

The Deuterium Fraction in Massive Starless Cores and Dynamical Implications

The Spatial Distribution of Complex Organic Molecules in the L1544 Pre-stellar Core

Theory of Gas Phase Scattering and Reactivity for Astrochemistry

Detectability of deuterated water in prestellar cores

Interstellar dust charging in dense molecular clouds: cosmic ray effects

Mid-J CO observations of Perseus B1-East 5: evidence for turbulent dissipation via low-velocity shocks

Next Generation Very Large Array Memo No. 6, Science Working Group 1: The Cradle of Life

The Deuteration Clock for Massive Starless Cores

The Deuterium Fractionation Timescale in Dense Cloud Cores: A Parameter Space Exploration

Chemistry and Radiative Transfer of Water in Cold, Dense Clouds

Detectability of Glycine in Solar-type System Precursors

Deuterated methanol in the pre-stellar core L1544

Deuterium chemistry of dense gas in the vicinity of low-mass and massive star forming regions

Deuterium Fractionation: the Ariadne's Thread from the Pre-collapse Phase to Meteorites and Comets today

Herschel HIFI observations of O$_2$ toward Orion: special conditions for shock enhanced emission

Massive Star Formation

The dynamical properties of dense filaments in the infrared dark cloud G035.39-00.33

The dynamics of collapsing cores and star formation

Correlating Infall with Deuterium Fractionation in Dense Cores

Deep observations of O2 toward a low-mass protostar with Herschel-HIFI

Detection of (15)NNH+ in L1544: non-LTE modelling of dyazenilium hyperfine line emission and accurate (14)N/(15)N values

The Dynamics of Massive Starless Cores with ALMA

A Virialized Filamentary Infrared Dark Cloud

ALMA reveals a chemically evolved submillimeter galaxy at z=4.76

Complex, Quiescent Kinematics in a Highly Filamentary Infrared Dark Cloud

First detection of water vapor in a pre-stellar core

H2D+ in the high mass star-forming Region Cygnus-X

Herschel-PACS observations of [OI]63um towards submillimetre galaxies at z~1

Our astrochemical heritage

Enhanced [CII] emission in a z=4.76 submillimetre galaxy

EVLA Observations of the Barnard 5 Star-Forming Core: Embedded Filaments Revealed

Herschel Measurements of Molecular Oxygen in Orion

Initial conditions for star formation in clusters: physical and kinematical structure of the starless core Oph A-N6

Mapping Large-Scale CO Depletion in a Filamentary Infrared Dark Cloud

Multi-line spectral imaging of dense cores in the Lupus molecular cloud

Observational studies of pre-stellar cores and infrared dark clouds

The enigmatic core L1451-mm: a first hydrostatic core? or a hidden VeLLO?

Detection of N15NH+ in L1544