Source author record

Mika Juvela

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

41works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

41 published item(s)

preprint2022arXiv

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): A Hot Corino Survey toward Protostellar Cores in the Orion Cloud

The presence of complex organic molecules (COMs) in the interstellar medium (ISM) is of great interest since it may link to the origin and prevalence of life in the universe. Aiming to investigate the occurrence of COMs and their possible origins, we conducted a chemical census toward a sample of protostellar cores as part of the ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP) project. We report the detection of 11 hot corino sources, which exhibit compact emissions from warm and abundant COMs, among 56 Class 0/I protostellar cores. All the hot corino sources discovered are likely Class 0 and their sizes of the warm region ($>$ 100 K) are comparable to 100 au. The luminosity of the hot corino sources exhibits positive correlations with the total number of methanol and the extent of its emissions. Such correlations are consistent with the thermal desorption picture for the presence of hot corino and suggest that the lower luminosity (Class 0) sources likely have a smaller region with COMs emissions. With the same sample selection method and detection criteria being applied, the detection rates of the warm methanol in the Orion cloud (15/37) and the Perseus cloud (28/50) are statistically similar when the cloud distances and the limited sample size are considered. Observing the same set of COM transitions will bring a more informative comparison between the cloud properties.

preprint2022arXiv

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): Evidence for a Molecular Jet Launched at an Unprecedented Early Phase of Protostellar evolution

Protostellar outflows and jets play a vital role in star formation as they carry away excess angular momentum from the inner disk surface, allowing the material to be transferred toward the central protostar. Theoretically, low velocity and poorly collimated outflows appear from the beginning of the collapse, at the first hydrostatic core (FHSC) stage. With growing protostellar core mass, high-density jets are launched which entrain an outflow from the infalling envelope. Until now, molecular jets have been observed at high velocity ($\gtrsim$ 100 km/s) in early Class\,0 protostars. We, for the first time, detect a dense molecular jet in SiO emission with small-velocity ($\sim$ 4.2 km\,s$^{-1}$, deprojected $\sim$ 24 km\,s$^{-1}$) from source G208.89-20.04Walma (hereafter, G208Walma) using ALMA Band\,6 observations. This object has some characteristics of FHSCs, such as a small outflow/jet velocity, extended 1.3\,mm continuum emission, and N$_2$D$^+$ line emission. Additional characteristics, however, are typical of early protostars: collimated outflow and SiO jet. The full extent of the outflow corresponds to a dynamical time scale of $\sim$ 930$^{+200}_{-100}$ years. The spectral energy distribution also suggests a very young source having an upper limit of T$_{bol}$ $\sim$ 31 K and L$_{bol}$ $\sim$ 0.8 L$_\sun$. We conclude that G208Walma is likely in the transition phase from FHSC to protostar, and the molecular jet has been launched within a few hundred years of initial collapse. Therefore, G208Walma may be the earliest object discovered in the protostellar phase with a molecular jet.

preprint2022arXiv

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): How do dense core properties affect the multiplicity of protostars?

During the transition phase from a prestellar to a protostellar cloud core, one or several protostars can form within a single gas core. The detailed physical processes of this transition, however, still remain unclear. We present 1.3 mm dust continuum and molecular line observations with the Atacama Large Millimeter/submillimeter Array (ALMA) toward 43 protostellar cores in the Orion Molecular Cloud Complex ($λ$ Orionis, Orion B, and Orion A) with an angular resolution of $\sim$ 0.35" ($\sim$ 140 au). In total, we detect 13 binary/multiple systems. We derive an overall multiplicity frequency (MF) of 28$\%$ $\pm$ 4$\%$ and a companion star fraction (CSF) of 51$\%$ $\pm$ 6$\%$, over a separation range of 300-8900 au. The median separation of companions is about 2100 au. The occurrence of stellar multiplicity may depend on the physical characteristics of the dense cores. Notably, those containing binary/multiple systems tend to show higher gas density and Mach number than cores forming single stars. The integral-shaped filament (ISF) of Orion A giant molecular cloud (GMC), which has the highest gas density and hosts high-mass star formation in its central region (the Orion Nebula cluster), shows the highest MF and CSF among the Orion GMCs. In contrast, the $λ$ Orionis Giant Molecular Cloud (GMC) has a lower MF and CSF than the Orion B and Orion A GMCs, indicating that feedback from HII regions may suppress the formation of multiple systems. We also find that the protostars comprising a binary/multiple system are usually at different evolutionary stages.

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

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- VII. A catalogue of SiO clumps from ACA observations

To understand the nature of SiO emission, we conducted ACA observations of the SiO (2-1) lines toward 146 massive star-forming regions, as part of the ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS) survey. We detected SiO emission in 128 (87.7$\%$) sources and identified 171 SiO clumps, 105 of which are spatially separated from 3 mm continuum emission. A large amount of the SiO line profiles (60$\%$) are non-Gaussian. The velocity dispersion of the SiO lines ranges from 0.3 to 5.43 km s$^{-1}$. In 63 sources the SiO clumps are associated with H$_\rm{II}$ regions characterized by H40$α$ emission. We find that 68$\%$ (116) of the SiO clumps are associated with strong outflows. The median velocity dispersion of the SiO line for outflow sources and non-outflow sources is 1.91 km s$^{-1}$ and 0.99 km s$^{-1}$, respectively. These results indicate that outflow activities could be connected to strongly shocked gas. The velocity dispersion and [SiO]/[H$^{13}$CO$^+$] intensity ratio do not show any correlation with the dust temperature and particle number density of clumps. We find a positive correlation between the SiO line luminosity and the bolometric luminosity, implying stronger shock activities are associated with more luminous proto-clusters. The SiO clumps in associations with H$_\rm{II}$ regions were found to show a steeper feature in $L_\rm{sio}$/$L_\rm{bol}$. The SiO line luminosity and the fraction of shocked gas have no apparent evidence of correlation with the evolutionary stages traced by luminosity to mass ratio ($L_\rm{bol}/M$).

preprint2022arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- VIII. A search for hot cores by using C$_2$H$_5$CN, CH$_3$OCHO and CH$_3$OH lines

Hot cores characterized by rich lines of complex organic molecules are considered as ideal sites for investigating the physical and chemical environments of massive star formation. We present a search for hot cores by using typical nitrogen- and oxygen-bearing complex organic molecules (C$_2$H$_5$CN, CH$_3$OCHO and CH$_3$OH), based on ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS). The angular resolutions and line sensitivities of the ALMA observations are better than 2 arcsec and 10 mJy/beam, respectively. A total of 60 hot cores are identified with 45 being newly detected, in which the complex organic molecules have high gas temperatures ($>$ 100 K) and small source sizes ($<$ 0.1 pc). So far this is the largest sample of hot cores observed with similar angular resolution and spectral coverage. The observations have also shown nitrogen and oxygen differentiation in both line emission and gas distribution in 29 hot cores. Column densities of CH$_3$OH and CH$_3$OCHO increase as rotation temperatures rise. The column density of CH$_3$OCHO correlates tightly with that of CH$_3$OH. The pathways for production of different species are discussed. Based on the spatial position difference between hot cores and UC~H{\sc ii} regions, we conclude that 24 hot cores are externally heated while the other hot cores are internally heated. The observations presented here will potentially help establish a hot core template for studying massive star formation and astrochemistry.

preprint2022arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- XI. From inflow to infall in hub-filament systems

We investigate the presence of hub-filament systems in a large sample of 146 active proto-clusters, using H$^{13}$CO$^{+}$ J=1-0 molecular line data obtained from the ATOMS survey. We find that filaments are ubiquitous in proto-clusters, and hub-filament systems are very common from dense core scales ($\sim$0.1 pc) to clump/cloud scales ($\sim$1-10 pc). The proportion of proto-clusters containing hub-filament systems decreases with increasing dust temperature ($T_d$) and luminosity-to-mass ratios ($L/M$) of clumps, indicating that stellar feedback from H{\sc ii} regions gradually destroys the hub-filament systems as proto-clusters evolve. Clear velocity gradients are seen along the longest filaments with a mean velocity gradient of 8.71 km s$^{-1}$pc$^{-1}$ and a median velocity gradient of 5.54 km s$^{-1}$pc$^{-1}$. We find that velocity gradients are small for filament lengths larger than $\sim$1~pc, probably hinting at the existence of inertial inflows, although we cannot determine whether the latter are driven by large-scale turbulence or large-scale gravitational contraction. In contrast, velocity gradients below $\sim$1~pc dramatically increase as filament lengths decrease, indicating that the gravity of the hubs or cores starts to dominate gas infall at small scales. We suggest that self-similar hub-filament systems and filamentary accretion at all scales may play a key role in high-mass star formation.

preprint2022arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- XII: Fragmentation and multi-scale gas kinematics in protoclusters G12.42+0.50 and G19.88-0.53

We present new continuum and molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS) survey for the two protoclusters, G12.42+0.50 and G19.88-0.53. The 3 mm continuum maps reveal seven cores in each of the two globally contracting protoclusters. These cores satisfy the radius-mass relation and the surface mass density criteria for high-mass star formation. Similar to their natal clumps, the virial analysis of the cores suggests that they are undergoing gravitational collapse ($\rm α_{vir} << 2$). The clump to core scale fragmentation is investigated and the derived core masses and separations are found to be consistent with thermal Jeans fragmentation. We detect large-scale filamentary structures with velocity gradients and multiple outflows in both regions. Dendrogram analysis of the H$^{13}$CO$^{+}$ map identifies several branch and leaf structures with sizes $\sim$ 0.1 and 0.03 pc, respectively. The supersonic gas motion displayed by the branch structures is in agreement with the Larson power-law indicating that the gas kinematics at this spatial scale is driven by turbulence. The transition to transonic/subsonic gas motion is seen to occur at spatial scales of $\sim$0.1 pc indicating the dissipation of turbulence. In agreement with this, the leaf structures reveal gas motions that deviate from the slope of Larson's law. From the large-scale converging filaments to the collapsing cores, the gas dynamics in G12.42+0.50 and G19.88-0.53 show scale-dependent dominance of turbulence and gravity and the combination of these two driving mechanisms needs to be invoked to explain massive star formation in the protoclusters.

preprint2022arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions-IX. A pilot study towards IRDC G034.43+00.24 on multi-scale structures and gas kinematics

We present a comprehensive study of the gas kinematics associated with density structures at different spatial scales in the filamentary infrared dark cloud, G034.43+00.24 (G34). This study makes use of the H13CO+ (1-0) molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS) survey, which has spatial and velocity resolution of 0.04 pc and 0.2 km/s, respectively. Several tens of dendrogram structures have been extracted in the position-position-velocity space of H13CO+, which include 21 small-scale leaves and 20 larger-scale branches. Overall, their gas motions are supersonic but they exhibit the interesting behavior where leaves tend to be less dynamically supersonic than the branches. For the larger-scale, branch structures, the observed velocity-size relation (i.e., velocity variation/dispersion versus size) are seen to follow the Larson scaling exponent while the smaller-scale, leaf structures show a systematic deviation and display a steeper slope. We argue that the origin of the observed kinematics of the branch structures is likely to be a combination of turbulence and gravity-driven ordered gas flows. In comparison, gravity-driven chaotic gas motion is likely at the level of small-scale leaf structures. The results presented in our previous paper and this current follow-up study suggest that the main driving mechanism for mass accretion/inflow observed in G34 varies at different spatial scales. We therefore conclude that a scale-dependent combined effect of turbulence and gravity is essential to explain the star-formation processes in G34.

preprint2022arXiv

Chemical exploration of Galactic cold cores

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

preprint2022arXiv

Nobeyama Survey of Inward Motions toward Cores in Orion Identified by SCUBA-2

In this study, 36 cores (30 starless and 6 protostellar) identified in Orion were surveyed to search for inward motions. We used the Nobeyama 45 m radio telescope, and mapped the cores in the $J = 1\rightarrow0$ transitions of HCO$^+$, H$^{13}$CO$^+$, N$_2$H$^+$, HNC, and HN$^{13}$C. The asymmetry parameter $δV$, which was the ratio of the difference between the HCO$^+$ and H$^{13}$CO$^+$ peak velocities to the H$^{13}$CO$^+$ line width, was biased toward negative values, suggesting that inward motions were more dominant than outward motions. Three starless cores (10% of all starless cores surveyed) were identified as cores with blue-skewed line profiles (asymmetric profiles with more intense blue-shifted emission), and another two starless cores (7%) were identified as candidate blue-skewed line profiles. The peak velocity difference between HCO$^+$ and H$^{13}$CO$^+$ of them was up to 0.9 km s$^{-1}$, suggesting that some inward motions exceeded the speed of sound for the quiescent gas ($\sim10-17$ K). The mean of $δV$ of the five aforementioned starless cores was derived to be $-$0.5$\pm$0.3. One core, G211.16$-$19.33North3, observed using the ALMA ACA in DCO$^+$ $J = 3\rightarrow2$ exhibited blue-skewed features. Velocity offset in the blue-skewed line profile with a dip in the DCO$^+$ $J = 3\rightarrow2$ line was larger ($\sim 0.5$ km s$^{-1}$) than that in HCO$^+$ $J = 1\rightarrow0$ ($\sim 0.2$ km s$^{-1}$), which may represent gravitational acceleration of inward motions. It seems that this core is at the last stage in the starless phase, judging from the chemical evolution factor version 2.0 (CEF2.0).

preprint2021arXiv

Molecular Cloud Cores with High Deuterium Fractions: Nobeyama Mapping Survey

We present the results of on-the-fly mapping observations of 44 fields containing 107 SCUBA-2 cores in the emission lines of molecules, N$_2$H$^+$, HC$_3$N, and CCS at 82$-$94 GHz using the Nobeyama 45-m telescope. This study aimed at investigating the physical properties of cores that show high deuterium fractions and might be close to the onset of star formation. We found that the distributions of the N$_2$H$^+$ and HC$_3$N line emissions are approximately similar to that of 850-$μ$m dust continuum emission, whereas the CCS line emission is often undetected or is distributed in a clumpy structure surrounding the peak position of the 850-$μ$m dust continuum emission. Occasionally (12%), we observe the CCS emission which is an early-type gas tracer toward the young stellar object, probably due to local high excitation. Evolution toward star formation does not immediately affect nonthermal velocity dispersion.

preprint2020arXiv

ALMA ACA and Nobeyama observations of two Orion cores in deuterated molecular lines

We mapped two molecular cloud cores in the Orion A cloud with the ALMA ACA 7-m Array and with the Nobeyama 45-m radio telescope. These cores have bright N$_2$D$^+$ emission in single-pointing observations with the Nobeyama 45-m radio telescope, have relatively high deuterium fraction, and are thought to be close to the onset of star formation. One is a star-forming core, and the other is starless. These cores are located along filaments observed in N$_2$H$^+$, and show narrow linewidths of 0.41 km s$^{-1}$ and 0.45 km s$^{-1}$ in N$_2$D$^+$, respectively, with the Nobeyama 45-m telescope. Both cores were detected with the ALMA ACA 7m Array in the continuum and molecular lines at Band 6. The starless core G211 shows clumpy structure with several sub-cores, which in turn show chemical differences. Also, the sub-cores in G211 have internal motions that are almost purely thermal. The starless sub-core G211D, in particular, shows a hint of the inverse P Cygni profile, suggesting infall motion. The star-forming core G210 shows an interesting spatial feature of two N$_2$D$^+$ peaks of similar intensity and radial velocity located symmetrically with respect to the single dust continuum peak. One interpretation is that the two N$_2$D$^+$ peaks represent an edge-on pseudo-disk. The CO outflow lobes, however, are not directed perpendicular to the line connecting both N$_2$D$^+$ peaks.

preprint2020arXiv

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): I. Detection of New Hot Corinos with ACA

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

preprint2020arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- I. Survey description and a first look at G9.62+0.19

The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation.

preprint2020arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- II. Compact objects in ACA observations and star formation scaling relations

We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales.

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

Molecular Cloud Cores with High Deuterium Fraction: Nobeyama Single-Pointing Survey

We present the results of a single-pointing survey of 207 dense cores embedded in Planck Galactic Cold Clumps distributed in five different environments ($λ$ Orionis, Orion A, B, Galactic plane, and high latitudes) to identify dense cores on the verge of star formation for the study of the initial conditions of star formation. We observed these cores in eight molecular lines at 76-94 GHz using the Nobeyama 45-m telescope. We find that early-type molecules (e.g., CCS) have low detection rates and that late-type molecules (e.g., N$_2$H$^+$, c-C$_3$H$_2$) and deuterated molecules (e.g., N$_2$D$^+$, DNC) have high detection rates, suggesting that most of the cores are chemically evolved. The deuterium fraction (D/H) is found to decrease with increasing distance, indicating that it suffers from differential beam dilution between the D/H pair of lines for distant cores ($>$1 kpc). For $λ$ Orionis, Orion A, and B located at similar distances, D/H is not significantly different, suggesting that there is no systematic difference in the observed chemical properties among these three regions. We identify at least eight high D/H cores in the Orion region and two at high latitudes, which are most likely to be close to the onset of star formation. There is no clear evidence of the evolutionary change in turbulence during the starless phase, suggesting that the dissipation of turbulence is not a major mechanism for the beginning of star formation as judged from observations with a beam size of 0.04 pc.

preprint2020arXiv

The Origin of Massive Stars: The Inertial-Inflow Model

We address the problem of the origin of massive stars, namely the origin, path and timescale of the mass flows that create them. Based on extensive numerical simulations, we propose a scenario where massive stars are assembled by large-scale, converging, inertial flows that naturally occur in supersonic turbulence. We refer to this scenario of massive-star formation as the "Inertial-Inflow Model". This model stems directly from the idea that the mass distribution of stars is primarily the result of turbulent fragmentation. Under this hypothesis, the statistical properties of the turbulence determine the formation timescale and mass of prestellar cores, posing definite constraints on the formation mechanism of massive stars. We quantify such constraints by the analysis of a simulation of supernova-driven turbulence in a 250-pc region of the interstellar medium, describing the formation of hundreds of massive stars over a time of approximately 30 Myr. Due to the large size of our statistical sample, we can say with full confidence that massive stars in general do not form from the collapse of massive cores, nor from competitive accretion, as both models are incompatible with the numerical results. We also compute synthetic continuum observables in Herschel and ALMA bands. We find that, depending on the distance of the observed regions, estimates of core mass based on commonly-used methods may exceed the actual core masses by up to two orders of magnitude, and that there is essentially no correlation between estimated and real core masses.

preprint2016arXiv

Astrochemical Properties of Planck Cold Clumps

We observed thirteen Planck cold clumps with the James Clerk Maxwell Telescope/SCUBA-2 and with the Nobeyama 45 m radio telescope. The N$_2$H$^+$ distribution obtained with the Nobeyama telescope is quite similar to SCUBA-2 dust distribution. The 82 GHz HC$_3$N, 82 GHz CCS, and 94 GHz CCS emission are often distributed differently with respect to the N$_2$H$^+$ emission. The CCS emission, which is known to be abundant in starless molecular cloud cores, is often very clumpy in the observed targets. We made deep single-pointing observations in DNC, HN$^{13}$C, N$_2$D$^+$, cyclic-C$_3$H$_2$ toward nine clumps. The detection rate of N$_2$D$^+$ is 50\%. Furthermore, we observed the NH$_3$ emission toward 15 Planck cold clumps to estimate the kinetic temperature, and confirmed that most of targets are cold ($\lesssim$ 20 K). In two of the starless clumps observe, the CCS emission is distributed as it surrounds the N$_2$H$^+$ core (chemically evolved gas), which resembles the case of L1544, a prestellar core showing collapse. In addition, we detected both DNC and N$_2$D$^+$. These two clumps are most likely on the verge of star formation. We introduce the Chemical Evolution Factor (CEF) for starless cores to describe the chemical evolutionary stage, and analyze the observed Planck cold clumps.

preprint2016arXiv

Composite biasing in Monte Carlo radiative transfer

Biasing or importance sampling is a powerful technique in Monte Carlo radiative transfer, and can be applied in different forms to increase the accuracy and efficiency of simulations. One of the drawbacks of the use of biasing is the potential introduction of large weight factors. We discuss a general strategy, composite biasing, to suppress the appearance of large weight factors. We use this composite biasing approach for two different problems faced by current state-of-the-art Monte Carlo radiative transfer codes: the generation of photon packages from multiple components, and the penetration of radiation through high optical depth barriers. In both cases, the implementation of the relevant algorithms is trivial and does not interfere with any other optimisation techniques. Through simple test models, we demonstrate the general applicability, accuracy and efficiency of the composite biasing approach. In particular, for the penetration of high optical depths, the gain in efficiency is spectacular for the specific problems that we consider: in simulations with composite path length stretching, high accuracy results are obtained even for simulations with modest numbers of photon packages, while simulations without biasing cannot reach convergence, even with a huge number of photon packages.

preprint2016arXiv

Supernova Driving. III. Synthetic Molecular Cloud Observations

We present a comparison of molecular clouds (MCs) from a simulation of supernova-driven interstellar medium (ISM) turbulence with real MCs from the Outer Galaxy Survey. The radiative transfer calculations to compute synthetic CO spectra are carried out assuming the CO relative abundance depends only on gas density, according to four different models. Synthetic MCs are selected above a threshold brightness temperature value, $T_{\rm B,min}=1.4$ K, of the $J=1-0$ $^{12}$CO line, generating 16 synthetic catalogs (four different spatial resolutions and four CO abundance models), each containing up to several thousands MCs. The comparison with the observations focuses on the mass and size distributions and on the velocity-size and mass-size Larson relations. The mass and size distributions are found to be consistent with the observations, with no significant variations with spatial resolution or chemical model, except in the case of the unrealistic model with constant CO abundance. The velocity-size relation is slightly too steep for some of the models, while the mass-size relation is a bit too shallow for all models only at a spatial resolution $dx\approx 1$ pc. The normalizations of the Larson relations show a clear dependence on spatial resolution, for both the synthetic and the real MCs. The comparison of the velocity-size normalization suggests that the SN rate in the Perseus arm is approximately 70\% or less of the rate adopted in the simulation. Overall, the realistic properties of the synthetic clouds confirm that supernova-driven turbulence can explain the origin and dynamics of MCs.

preprint2015arXiv

Benchmarking the Calculation of Stochastic Heating and Emissivity of Dust Grains in the Context of Radiative Transfer Simulations

We define an appropriate problem for benchmarking dust emissivity calculations in the context of radiative transfer (RT) simulations, specifically including the emission from stochastically heated dust grains. Our aim is to provide a self-contained guide for implementors of such functionality, and to offer insights in the effects of the various approximations and heuristics implemented by the participating codes to accelerate the calculations. The benchmark problem definition includes the optical and calorimetric material properties, and the grain size distributions, for a typical astronomical dust mixture with silicate, graphite and PAH components; a series of analytically defined radiation fields to which the dust population is to be exposed; and instructions for the desired output. We process this problem using six RT codes participating in this benchmark effort, and compare the results to a reference solution computed with the publicly available dust emission code DustEM. The participating codes implement different heuristics to keep the calculation time at an acceptable level. We study the effects of these mechanisms on the calculated solutions, and report on the level of (dis)agreement between the participating codes. For all but the most extreme input fields, we find agreement within 10% across the important wavelength range from 3 to 1000 micron. We conclude that the relevant modules in RT codes can and do produce fairly consistent results for the emissivity spectra of stochastically heated dust grains.

preprint2015arXiv

Can we trace very cold dust from its emission alone ?

Context. Dust is a good tracer of cold dark clouds but its column density is difficult to quantify. Aims. We want to check whether the far-infrared and submillimeter high-resolution data from Herschel SPIRE and PACS cameras combined with ground-based telescope bolometers allow us to retrieve the whole dust content of cold dark clouds. Methods. We compare far-infrared and submillimeter emission across L183 to the 8 $μ$m absorption map from Spitzer data and fit modified blackbody functions towards three different positions. Results. We find that none of the Herschel SPIRE channels follow the cold dust profile seen in absorption. Even the ground-based submillimeter telescope observations, although more closely following the absorption profile, cannot help to characterize the cold dust without external information such as the dust column density itself. The difference in dust opacity can reach up to a factor of 3 in prestellar cores of high extinction. Conclusions. In dark clouds, the amount of very cold dust cannot be measured from its emission alone. In particular, studies of dark clouds based only on Herschel data can miss a large fraction of the dust content. This has an impact on core and filament density profiles, masse and stability estimates.

preprint2015arXiv

Planck cold clumps in the $λ$ Orionis complex: I. Discovery of an extremely young Class 0 protostellar object and a proto-brown dwarf candidate in a bright rimmed clump PGCC G192.32-11.88

We are performing a series of observations with ground-based telescopes toward Planck Galactic cold clumps (PGCCs) in the $λ$ Orionis complex in order to systematically investigate the effects of stellar feedback. In the particular case of PGCC G192.32-11.88, we discovered an extremely young Class 0 protostellar object (G192N) and a proto-brown dwarf candidate (G192S). G192N and G192S are located in a gravitationally bound bright-rimmed clump. The velocity and temperature gradients seen in line emission of CO isotopologues indicate that PGCC G192.32-11.88 is externally heated and compressed. G192N probably has the lowest bolometric luminosity ($\sim0.8$ L$_{\sun}$) and accretion rate (6.3$\times10^{-7}$ M$_{\sun}$~yr$^{-1}$) when compared with other young Class 0 sources (e.g. PACS Bright Red sources (PBRs)) in the Orion complex. It has slightly larger internal luminosity ($0.21\pm0.01$ L$_{\sun}$) and outflow velocity ($\sim$14 km~s$^{-1}$) than the predictions of first hydrostatic cores (FHSCs). G192N might be among the youngest Class 0 sources, which are slightly more evolved than a FHSC. Considering its low internal luminosity ($0.08\pm0.01$ L$_{\odot}$) and accretion rate (2.8$\times10^{-8}$ M$_{\sun}$~yr$^{-1}$), G192S is an ideal proto-brown dwarf candidate. The star formation efficiency ($\sim$0.3\%-0.4\%) and core formation efficiency ($\sim$1\%) in PGCC G192.32-11.88 are significantly smaller than in other giant molecular clouds or filaments, indicating that the star formation therein is greatly suppressed due to stellar feedback.

preprint2014arXiv

Dust properties inside molecular clouds from coreshine modeling and observations

Context. Using observations to deduce dust properties, grain size distribution, and physical conditions in molecular clouds is a highly degenerate problem. Aims. The coreshine phenomenon, a scattering process at 3.6 and 4.5 $μ$m that dominates absorption, has revealed its ability to explore the densest parts of clouds. We want to use this effect to constrain the dust parameters. The goal is to investigate to what extent grain growth (at constant dust mass) inside molecular clouds is able to explain the coreshine observations. We aim to find dust models that can explain a sample of Spitzer coreshine data. We also look at the consistency with near-infrared data we obtained for a few clouds. Methods. We selected four regions with a very high occurrence of coreshine cases: Taurus-Perseus, Cepheus, Chameleon and L183/L134. We built a grid of dust models and investigated the key parameters to reproduce the general trend of surface bright- nesses and intensity ratios of both coreshine and near-infrared observations with the help of a 3D Monte-Carlo radiative transfer code. The grid parameters allow to investigate the effect of coagulation upon spherical grains up to 5 $μ$m in size derived from the DustEm diffuse interstellar medium grains. Fluffiness (porosity or fractal degree), ices, and a handful of classical grain size distributions were also tested. We used the near- and mostly mid-infrared intensity ratios as strong discriminants between dust models. Results. The determination of the background field intensity at each wavelength is a key issue. In particular, an especially strong background field explains why we do not see coreshine in the Galactic plane at 3.6 and 4.5 $μ$m. For starless cores, where detected, the observed 4.5 $μ$m / 3.6 $μ$m coreshine intensity ratio is always lower than $\sim$0.5 which is also what we find in the models for the Taurus-Perseus and L183 directions. Embedded sources can lead to higher fluxes (up to four times greater than the strongest starless core fluxes) and higher coreshine ratios (from 0.5 to 1.1 in our selected sample). Normal interstellar radiation field conditions are sufficient to find suitable grain models at all wavelengths for starless cores. The standard interstellar grains are not able to reproduce observations and, due to the multi-wavelength approach, only a few grain types meet the criteria set by the data. Porosity does not affect the flux ratios while the fractal dimension helps to explain coreshine ratios but does not seem able to reproduce near-infrared observations without a mix of other grain types. Conclusions. Combined near- and mid-infrared wavelengths confirm the potential to reveal the nature and size distribution of dust grains. Careful assessment of the environmental parameters (interstellar and background fields, embedded or nearby reddened sources) is required to validate this new diagnostic.

preprint2014arXiv

Follow-up observations toward Planck cold clumps with ground-based radio telescopes

The physical and chemical properties of prestellar cores, especially massive ones, are still far from being well understood due to the lack of a large sample. The low dust temperature ($<$14 K) of Planck cold clumps makes them promising candidates for prestellar objects or for sources at the very initial stages of protostellar collapse. We have been conducting a series of observations toward Planck cold clumps (PCCs) with ground-based radio telescopes. In general, when compared with other star forming samples (e.g. infrared dark clouds), PCCs are more quiescent, suggesting that most of them may be in the earliest phase of star formation. However, some PCCs are associated with protostars and molecular outflows, indicating that not all PCCs are in a prestellar phase. We have identified hundreds of starless dense clumps from the mapping survey with the Purple Mountain Observatory (PMO) 13.7-m telescope. Follow-up observations suggest that these dense clumps are ideal targets to search for prestellar objects.

preprint2013arXiv

Determination of the far-infrared dust opacity in a prestellar core

(abridged) We correlated near-infrared stellar H-Ks colour excesses of background stars from NTT/SOFI with the far-IR optical depth map, tauFIR, derived from Herschel 160, 250, 350, and 500 um data. The Herschel maps were also used to construct a model for the cloud to examine the effect of temperature gradients on the estimated optical depths and dust absorption cross-sections. A linear correlation is seen between the colour H-Ks and tauFIR up to high extinctions (AV ~ 25). The correlation translates to the average extinction ratio A250um/AJ = 0.0014 +/- 0.0002, assuming a standard near-infrared extinction law and a dust emissivity index beta=2. Using an empirical NH/AJ ratio we obtain an average absorption cross-section per H nucleus of sigmaH(250um) = (1.8 +/- 0.3) * 10^(-25) cm^2 / H-atom, corresponding to a cross-section per unit mass of gas kappaG(250 um) = 0.08 +/- 0.01 cm^2 / g. The cloud model however suggests that owing to the bias caused by temperature changes along the line-of-sight these values underestimate the true cross-sections by up to 40% near the centre of the core. Assuming that the model describes the effect of the temperature variation on tauFIR correctly, we find that the relationship between H-Ks and tauFIR agrees with the recently determined relationship between sigmaH and NH in Orion A. The derived far-IR cross-section agrees with previous determinations in molecular clouds with moderate column densities, and is not particularly large compared with some other cold cores. We suggest that this is connected to the core not beng very dense (the central density is likely to be ~10^5 cm^-3) and judging from previous molecular line data, it appears to be at an early stage of chemical evolution.

preprint2013arXiv

Luminous Infrared Galaxies With the Submillimeter Array. IV. 12CO J=6-5 Observations of VV 114

We present high-resolution (~2.5") observations of 12CO J=6-5 towards the luminous infrared galaxy VV 114 using the Submillimeter Array. We detect 12CO J=6-5 emission from the eastern nucleus of VV 114 but do not detect the western nucleus or the central region. We combine the new 12CO J=6-5 observations with previously published or archival low-J CO observations, that include 13CO J=1-0 Atacama Large Millimeter/submillimeter Array cycle 0 observations, to analyze the beam-averaged physical conditions of the molecular gas in the eastern nucleus. We use the radiative transfer code RADEX and a Bayesian likelihood code to constrain the temperature (T_kin), density (n(H2)) and column density (N(12CO)) of the molecular gas. We find that the most probable scenario for the eastern nucleus is a cold (T_kin = 38 K), moderately dense (n(H2) = 10^2.89 cm^-3) molecular gas component. We find the most probable 12CO to 13CO abundance ratio ([12CO]/[13CO]) is 229, roughly three times higher than the Milky Way value. This high abundance ratio may explain the observed high 12CO/ 13CO line ratio (> 25). The unusual 13CO J=2-1/J=1-0 line ratio of 0.6 is produced by a combination of moderate 13CO optical depths (tau = 0.4 - 1.1) and extremely subthermal excitation temperatures. We measure the CO-to-H2 conversion factor, alpha_co to be 0.5 M_sol (K km s^-1 pc^2)^-1, which agrees with the widely used factor for ultra luminous infrared galaxies of Downes & Solomon (1998; alpha_co =0.8 M_sol (K km s^-1 pc^2)^-1).

preprint2013arXiv

Variation in dust properties in a dense filament of the Taurus molecular complex (L1506)

We observed the L1506 filament, which is located in the Taurus molecular complex, with the Herschel PACS and SPIRE instruments. Our aim is to prove the variation in grain properties along the entire length of the filament. In particular, we want to determine above which gas density this variation arises and what changes in the grain optical properties/size distribution are required. We use the 3D radiative transfer code CRT, coupled to the dust emission and extinction code DustEM, to model the emission and extinction of the dense filament. We test a range of optical properties and size distributions for the grains: dust of the diffuse interstellar medium (interstellar PAHs and amorphous carbons and silicates) and both compact and fluffy aggregates. We find that the grain opacity has to increase across the filament to fit simultaneously the near-IR extinction and Herschel emission profiles of L1506. We interpret this change to be a consequence of the coagulation of dust grains to form fluffy aggregates. Grains similar to those in the diffuse medium have to be present in the outer layers of the cloud, whereas aggregates must prevail above gas densities of a few 1000 H/cm3. This corresponds to line-of-sights with visual extinction in the V band of the order of 2 to 3. The dust opacity at 250 microns is raised by a factor of 1.8 to 2.2, while the grain average size is increased by a factor of 5. These exact numbers depend naturally on the dust model chosen to fit the data. Our findings agree with the constraints given by the study of the gas molecular lines. Using a simple approach, we show that the aggregates may have time to form inside the filament within the cloud lifetime. Our model also characterises the density structure of the filament, showing that the filament width is not constant along L1506 but instead varies by a factor of the order of 4.

preprint2012arXiv

Luminous Infrared Galaxies With the Submillimeter Array. III. The Dense Kiloparsec Molecular Concentrations of Arp 299

We have used high resolution (~2.3") observations of the local (D = 46 Mpc) luminous infrared galaxy Arp 299 to map out the physical properties of the molecular gas which provides the fuel for its extreme star formation activity. The 12CO J=3-2, 12CO J=2-1 and 13CO J=2-1 lines were observed with the Submillimeter Array and the short spacings of the 12CO J=2-1 and J=3-2 observations have been recovered using James Clerk Maxwell Telescope single dish observations. We use the radiative transfer code RADEX to estimate the physical properties (density, column density and temperature) of the different regions in this system. The RADEX solutions of the two galaxy nuclei, IC 694 and NGC 3690, are consistent with a wide range of gas components, from warm moderately dense gas with T_{kin} > 30 K and n(H_{2}) ~ 0.3 - 3 x 10^{3} cm^{-3} to cold dense gas with T_{kin} ~ 10-30 K and n(H_{2}) > 3 x 10^{3} cm^{-3}. The overlap region is shown to have a better constrained solution with T_{\rm{kin}}$ ~ 10-50 K and n(H_{2}) ~ 1-30 x 10^{3} cm^{-3}. We estimate the gas masses and star formation rates of each region in order to derive molecular gas depletion times. The depletion times of all regions (20-60 Myr) are found to be about 2 orders of magnitude lower than those of normal spiral galaxies. This rapid depletion time can probably be explained by a high fraction of dense gas on kiloparsec scales in Arp 299. We estimate the CO-to-H_{2} factor, α_{co} to be 0.4 \pm 0.3 (3 x 10^{-4}/ x_{CO}) M_{sol} (K km s^{-1} pc^{2})^{-1} for the overlap region. This value agrees well with values determined previously for more advanced merger systems.

preprint2011arXiv

A (sub)millimetre study of dense cores in Orion B9

We aim to further constrain the properties and evolutionary stages of dense cores in Orion B9. The central part of Orion B9 was mapped at 350 micron with APEX/SABOCA. A sample of nine cores in the region were observed in C17O(2-1), H13CO+(4-3) (towards 3 sources), DCO+(4-3), N2H+(3-2), and N2D+(3-2) with APEX/SHFI. These data are used in conjunction with our previous APEX/LABOCA 870-micron dust continuum data. Many of the LABOCA cores show evidence of substructure in the higher-resolution SABOCA image. In particular, we report on the discovery of multiple very low-mass condensations in the prestellar core SMM 6. Based on the 350-to-870 micron flux density ratios, we determine dust temperatures of ~7.9-10.8 K, and dust emissivity indices of ~0.5-1.8. The CO depletion factors are in the range ~1.6-10.8. The degree of deuteration in N2H+ is ~0.04-0.99, where the highest value (seen towards the prestellar core SMM 1) is, to our knowledge, the most extreme level of N2H+ deuteration reported so far. The level of HCO+ deuteration is about 1-2%. We also detected D2CO towards two sources. The detection of subcondensations within SMM 6 shows that core fragmentation can already take place during the prestellar phase. The origin of this substructure is likely caused by thermal Jeans fragmentation of the elongated parent core. A low depletion factor and the presence of gas-phase D2CO in SMM 1 suggest that the core chemistry is affected by the nearby outflow. The very high N2H+ deuteration in SMM 1 is likely to be remnant of the earlier CO-depleted phase.

preprint2011arXiv

CH abundance gradient in TMC-1

We observed the 9-cm Lambda-doubling lines of CH along the dense filament of TMC-1. The CH column densities were compared with the total H2 column densities derived using the 2MASS NIR data and previously published SCUBA maps and with OH column densities derived using previous observations with Effelsberg. We also modelled the chemical evolution of TMC-1 adopting physical conditions typical of dark clouds using the UMIST Database for Astrochemistry gas-phase reaction network to aid the interpretation of the observed OH/CH abundance ratios. The CH column density has a clear peak in the vicinity of the cyanopolyyne maximum of TMC-1. The fractional CH abundance relative to H2 increases steadily from the northwestern end of the filament where it lies around 1.0e-8, to the southeast where it reaches a value of 2.0e-8. The OH and CH column densities are well correlated, and we obtained OH/CH abundance ratios of ~ 16 - 20. These values are clearly larger than what has been measured recently in diffuse interstellar gas and is likely to be related to C to CO conversion at higher densities. The good correlation between CH and OH can be explained by similar production and destruction pathways. We suggest that the observed CH and OH abundance gradients are mainly due to enhanced abundances in a low-density envelope which becomes more prominent in the southeastern part and seems to continue beyond the dense filament. An extensive envelope probably signifies an early stage of dynamical evolution, and conforms with the detection of a large CH abundance in the southeastern part of the cloud. The implied presence of other simple forms of carbon in the gas phase provides a natural explanation for the observation of "early-type" molecules in this region.

preprint2011arXiv

Mapping the prestellar core Ophiuchus D (L1696A) in ammonia

The gas kinetic temperature in the centres of starless, high-density cores is predicted to fall as low as 5-6 K. The aim of this study was to determine the kinetic temperature distribution in the low-mass prestellar core Oph D where previous observations suggest a very low central temperature. The densest part of the Oph D core was mapped in the NH3(1,1) and (2,2) inversion lines using the Very Large Array (VLA). The physical quantities were derived from the observed spectra by fitting the hyperfine structure of the lines, and subsequently the temperature distribution of Oph D was calculated using the standard rotational temperature techniques. A physical model of the cores was constructed, and the simulated spectra after radiative transfer calculations with a 3D Monte Carlo code were compared with the observed spectra. Temperature, density, and ammonia abundance of the core model were tuned until a satisfactory match with the observation was obtained. The high resolution of the interferometric data reveals that the southern part of Oph D comprises of two small cores. The gas kinetic temperatures, as derived from ammonia towards the centres of the southern and northern core are 7.4 and 8.9 K, respectively. The observed masses of the cores are only 0.2 M_Sun. Their potential collapse could lead to formation of brown dwarfs or low-mass stars.

preprint2010arXiv

Physical properties of dense cores in Orion B9

We aim to determine the physical and chemical properties of dense cores in Orion B9. We observed the NH3(1,1) and (2,2), and the N2H+(3-2) lines towards the submm peak positions. These data are used in conjunction with our LABOCA 870 micron dust continuum data. The gas kinetic temperature in the cores is between ~9.4-13.9 K. The non-thermal velocity dispersion is subsonic in most of the cores. The non-thermal linewidth in protostellar cores appears to increase with increasing bolometric luminosity. The core masses are very likely drawn from the same parent distribution as the core masses in Orion B North. Starless cores in the region are likely to be gravitationally bound, and thus prestellar. Some of the cores have a lower radial velocity than the systemic velocity of the region, suggesting that they are members of the "low-velocity part" of Orion B. The observed core-separation distances deviate from the corresponding random-like model distributions. The distances between the nearest-neighbours are comparable to the thermal Jeans length. The fractional abundances of NH3 and N2H+ in the cores are ~1.5-9.8x10^{-8} and ~0.2-5.9x10^{-10}, respectively. The NH3 abundance appears to decrease with increasing H2 column and number densities. The NH3/N2H+ column density ratio is larger in starless cores than in cores with embedded protostars. The core population in Orion B9 is comparable in physical properties to those in nearby low-mass star-forming regions. It is unclear if the origin of cores could be explained by turbulent fragmentation. On the other hand, many of the core properties conform with the picture of dynamic core evolution. The Orion B9 region has probably been influenced by the feedback from the nearby Ori OB 1b group, and the fragmentation of the parental cloud into cores could be caused by gravitational instability.

preprint2009arXiv

A Corona Australis cloud filament seen in NIR scattered light II: Comparison with sub-millimeter data

We study a northern part of the Corona Australis molecular cloud that consists of a filament and a dense sub-millimetre core inside the filament. Our aim is to measure dust temperature and sub-mm emissivity within the region. We also look for confirmation that near-infrared (NIR) surface brightness can be used to study the structure of even very dense clouds. We extend our previous NIR mapping south of the filament. The dust colour temperatures are estimated using Spitzer 160um and APEX/Laboca 870um maps. The column densities derived based on the reddening of background stars, NIR surface brightness, and thermal sub-mm dust emission are compared. A three dimensional toy model of the filament is used to study the effect of anisotropic illumination on near-infrared surface brightness and the reliability of dust temperature determination. Relative to visual extinction, the estimated emissivity at 870um is kappa(870) = (1.3 +- 0.4) x 10^{-5} 1/mag. This is similar to the values found in diffuse medium. A significant increase in the sub-millimetre emissivity seems to be excluded. In spite of saturation, NIR surface brightness was able to accurately pinpoint, and better than measurements of the colour excesses of background stars, the exact location of the column density maximum. Both near- and far-infrared data show that the intensity of the radiation field is higher south of the filament.

preprint2009arXiv

A far-infrared survey at the North Galactic Pole I: Nearby star-forming galaxies and the effect of confused sources on source counts

We present follow-up observations of the far-infrared (FIR) sources at 90, 150 and 180 mu detected as part of the ISOPHOT EBL project, which has recently measured the absolute surface brightness of the cosmic infrared background radiation (CIRB) for the first time independently from COBE data. We have observed the fields at the North Galactic Pole region in the optical and near-IR, and complement these data with SDSS photometry, and spectroscopy where available, and present identifications of the 25 FIR sources which reach down to ~150 mJy in all three ISOPHOT bands. Identifications are done by means of full spectral energy density fitting to all sources in the FIR error circle areas. Approximately 80 per cent are identified as star-forming or star-bursting galaxies at z<0.3. We also find that more than half of the counterparts have disturbed morphologies, with signs of past or present interactions. However, only 20 per cent of all the sources are uniquely matched with a single galaxy -- 40 per cent are blends of two or more of these nearby star-forming galaxies, while another 20 per cent are likely blends of nearby and fainter galaxies. The final 20 per cent are likely to be more luminous IR galaxies at higher redshifts. The blended sources have an effect on the FIR source counts. In particular, taking into account realistic confusion or blending of sources, the differential FIR counts move down by a factor of ~1.5 and steepen in the 100 to 400 mJy range.

preprint2009arXiv

The Power Spectrum of Turbulence in NGC 1333: Outflows or Large-Scale Driving?

Is the turbulence in cluster-forming regions internally driven by stellar outflows or the consequence of a large-scale turbulent cascade? We address this question by studying the turbulent energy spectrum in NGC 1333. Using synthetic 13CO maps computed with a snapshot of a supersonic turbulence simulation, we show that the VCS method of Lazarian and Pogosyan provides an accurate estimate of the turbulent energy spectrum. We then apply this method to the 13CO map of NGC 1333 from the COMPLETE database. We find the turbulent energy spectrum is a power law, E(k) k^-beta, in the range of scales 0.06 pc < ell < 1.5 pc, with slope beta=1.85\pm 0.04. The estimated energy injection scale of stellar outflows in NGC 1333 is ell_inj 0.3 pc, well resolved by the observations. There is no evidence of the flattening of the energy spectrum above the scale ell_inj predicted by outflow-driven simulations and analytical models. The power spectrum of integrated intensity is also a nearly perfect power law in the range of scales 0.16 pc < ell < 7.9 pc, with no feature above ell_inj. We conclude that the observed turbulence in NGC 1333 does not appear to be driven primarily by stellar outflows.

preprint2009arXiv

The Super-Alfvénic Model of Molecular Clouds: Predictions for Mass-to-Flux and Turbulent-to-Magnetic Energy Ratios

Recent measurements of the Zeeman effect in dark-cloud cores provide important tests for theories of cloud dynamics and prestellar core formation. In this Letter we report results of simulated Zeeman measurements, based on radiative transfer calculations through a snapshot of a simulation of supersonic and super-Alfvénic turbulence. We have previously shown that the same simulation yields a relative mass-to-flux ratio (core versus envelope) in agreement with the observations (and in contradiction with the ambipolar-drift model of core formation). Here we show that the mass-to-flux and turbulent-to-magnetic-energy ratios in the simulated cores agree with observed values as well. The mean magnetic field strength in the simulation is very low, \bar{B}=0.34 \muG, presumably lower than the mean field in molecular clouds. Nonetheless, high magnetic field values are found in dense cores, in agreement with the observations (the rms field, amplified by the turbulence, is B_{rms}=3.05 \muG). We conclude that a strong large-scale mean magnetic field is not required by Zeeman effect measurements to date, although it is not ruled out by this work.

preprint2008arXiv

The Super-Alfvénic Model of Molecular Clouds: Predictions for Zeeman Splitting Measurements

We present synthetic OH Zeeman splitting measurements of a super-Alfvenic model of molecular clouds. We select dense cores from synthetic 13CO maps computed from the largest simulation to date of supersonic and super-Alfvenic turbulence. The synthetic Zeeman splitting measurements in the cores yield a relation between the magnetic field strength, B, and the column density, N, in good agreement with the observations. The large scatter in B at a fixed value of N is partly due to intrinsic variations in the magnetic field strength from core to core. We also compute the relative mass-to-flux ratio between the center of the cores and their envelopes, ${\cal R}_μ$, and show that super-Alfvenic turbulence produces a significant scatter also in ${\cal R}_μ$, including negative values (field reversal between core center and envelope). We find ${\cal R}_μ<1$ for 70% of the cores, and ${\cal R}_μ<0$ for 12%. Of the cores with $|B_{\rm LOS}|>10$ \muG, 81% have ${\cal R}_μ<1$. These predictions of the super-Alfvenic model are in stark contrast to the ambipolar drift model of core formation, where only ${\cal R}_μ>1$ is allowed.

preprint1999arXiv

Super-Sonic Turbulence in the Perseus Molecular Cloud

We compare the statistical properties of J=1-0 13CO spectra observed in the Perseus Molecular Cloud with synthetic J=1-0 13CO spectra, computed solving the non-LTE radiative transfer problem for a model cloud obtained as solutions of the three dimensional magneto-hydrodynamic (MHD) equations. The model cloud is a randomly forced super-Alfvenic and highly super-sonic turbulent isothermal flow. The purpose of the present work is to test if idealized turbulent flows, without self-gravity, stellar radiation, stellar outflows, or any other effect of star formation, are inconsistent or not with statistical properties of star forming molecular clouds. We present several statistical results that demonstrate remarkable similarity between real data and the synthetic cloud. Statistical properties of molecular clouds like Perseus are appropriately described by random super-sonic and super-Alfvenic MHD flows. Although the description of gravity and stellar radiation are essential to understand the formation of single protostars and the effects of star formation in the cloud dynamics, the overall description of the cloud and of the initial conditions for star formation can apparently be provided on intermediate scales without accounting for gravity, stellar radiation, and a detailed modeling of stellar outflows. We also show that the relation between equivalent line width and integrated antenna temperature indicates the presence of a relatively strong magnetic field in the core B1, in agreement with Zeeman splitting measurements.