Source author record

Aina Palau

Aina Palau 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

50works
4topics
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

50 published item(s)

preprint2026arXiv

Synthetic polarization observations of magnetized pillars in HII regions: The limitations of the Davis-Chandrasekhar-Fermi method

We investigated the morphology and strength of magnetic fields in pillar-shaped structures at the boundaries of HII regions by combining three-dimensional radiation-magnetohydrodynamic, R-MHD, simulations with synthetic polarimetric observations. Our analysis focuses on the first pillar formed self-consistently in the simulation and is used as a proof of concept to test the applicability of the Davis-Chandrasekhar-Fermi (DCF) method under conditions dominated by external agents. The pillar arises as the ionization front compresses a dense clump, producing a magnetically aligned, elongated structure whose morphology and field configuration closely resemble those observed in real systems such as the pillars of M16. Synthetic dust-polarization maps at 850 $μ$m reproduce the large-scale magnetic morphology of the simulated pillar, confirming that polarimetry is a reliable tracer of magnetic field geometry. However, magnetic field strengths inferred using DCF-type methods systematically overestimate intrinsic values from the simulations by factors $\gtrsim 2$. We attribute this discrepancy to the fact that field alignment and amplification are primarily driven by external gas pressure from the expanding HII region rather than internal turbulence, thus violating key assumptions of the DCF method. Our results highlight the need for caution when applying classical DCF-based analysis to pillars or other structures shaped by external compressions.

preprint2026arXiv

Unveiling the First O-Type Bloated Star Candidate through ALMA and EVLA Observations

We investigate the circumstellar environment of the O-type bloated star candidate IRAS 19520+2759 (I19520) using high-resolution observations from the Atacama Large Millimeter/submillimeter Array (ALMA) and the Expanded Very Large Array (EVLA). Radio continuum emission traced by the EVLA (C, K, and Q bands) exhibits a spectral index of 0.5, consistent with a thermal jet. ALMA 1.3 mm continuum map reveals a compact source coincident with the optical counterpart of I19520, likely tracing the dense core hosting the central massive young stellar object. A prominent molecular outflow in the east-west direction, along with a possible secondary outflow oriented northeast-southwest, is identified in the $^{13}\mathrm{CO}$ emission. A hot molecular core and a Keplerian disk are detected in several $\mathrm{SO}_2$ transitions. Assuming an edge-on disk geometry, the dynamical mass of the central object is estimated to be in the range of $10$-$15~M_\odot$.

preprint2022arXiv

Different molecular filament widths as tracers of accretion onto filaments

We explore how dense filament widths, when measured using different molecular species, may change as a consequence of gas accretion toward the filament. As a gas parcel falls into the filament, it will experience different density, temperature, and extinction values. The rate at which this environment changes will affect differently the abundance of different molecules. So, a molecule that forms quickly will better reflect the local physical conditions a gas parcel experiences than a slower-forming molecule. Since these differences depend on how the respective timescales compare, the different molecular distributions should reflect how rapidly the environment changes, i.e., the accretion rate toward the filament. We find that the filament widths measured from time-dependent abundances for C2H, CO, CN, CS, and C3H2, are sensitive the most to this effect, being those molecules the ones presenting also the wider filament widths. On the contrary, molecules such as N2H+, NH3, H2CO, HNC or CH3OH are not so sensitive to accretion and present the narrowest filament widths. We propose that ratios of filament widths for different tracers could be a useful tool to estimate the accretion rate onto the filament.

preprint2021arXiv

ALMA observations of NGC 6334S. II. Subsonic and Transonic Narrow Filaments in a High-mass Star Formation Cloud

We present a study of narrow filaments toward a massive infrared dark cloud, NGC 6334S, using the Atacama Large Millimeter/submillimeter Array (ALMA). Thirteen gas filaments are identified using the H$^{13}$CO$^{+}$ line, while a single continuum filament is revealed by the continuum emission. The filaments present a compact radial distribution with a median filament width of $\sim$0.04 pc narrower than the previously proposed `quasi-universal' 0.1~pc filament width. The higher spatial resolution observations and higher-density gas tracer tend to identify even narrower and lower mass filaments. The filament widths are roughly twice the size of embedded cores. The gas filaments are largely supported by thermal motions. The nonthermal motions are predominantly subsonic and transonic in both identified gas filaments and embedded cores, which may imply that stars are likely born in environments of low turbulence. A fraction of embedded objects show a narrower velocity dispersion compared with their corresponding natal filaments, which may indicate that the turbulent dissipation is taking place in these embedded cores. The physical properties (mass, mass per unit length, gas kinematics, and width) of gas filaments are analogous to those of narrow filaments found in low- to high-mass star-forming regions. The more evolved sources are found to be farther away from the filaments, a situation that may have resulted from the relative motions between the YSOs and their natal filaments.

preprint2021arXiv

Density profile evolution during prestellar core collapse: Collapse starts at the large scale

We study the gravitationally-dominated, accretion-driven evolution of a prestellar core. In our model, as the core's density increases, it remains immersed in a constant-density environment and so it accretes from this environment, increasing its mass and reducing its Jeans length. Assuming a power-law density profile $ρ\propto r^{-p}$, we compute the rate of change of the slope $p$, and show that the value $p=2$ is stationary, and furthermore, an attractor. The radial profile of the Jeans length scales as $r^{p/2}$, implying that, for $p<2$, there is a radius below which the region is smaller than its Jeans length, thus appearing gravitationally stable and in need of pressure confinement, while, in reality, it is part of a larger-scale collapse and is undergoing compression by the infalling material. In this region, the infall speed decreases towards the center, eventually becoming subsonic, thus appearing "coherent", without the need for turbulence dissipation. We present a compilation of observational determinations of density profiles in dense cores and show that the distribution of their slopes peaks at $p \sim 1.7$--1.9, supporting the notion that the profile steepens over time. Finally, we discuss the case of magnetic support in a core in which the field scales as $B \propto ρ^β$. For the expected value of $β= 2/3$, this implies that the mass to magnetic flux ratio also decreases towards the central parts of the cores, making them appear magnetically supported, while in reality they may be part of larger collapsing supercritical region. We conclude that local signatures of either thermal or magnetic support are not conclusive evidence of stability, that the gravitational instability of a region must be established at the large scales, and that the prestellar stage of collapse is dynamic rather than quasistatic.

preprint2020arXiv

ALMA observations of NGC 6334S $-$ I: Forming massive stars and cluster in subsonic and transonic filamentary clouds

We present Atacama Large Millimeter/submillimeter Array (ALMA) and Karl G. Jansky Very Large Array (JVLA) observations of the massive infrared dark cloud NGC 6334S (also known as IRDC G350.56+0.44), located at the southwestern end of the NGC 6334 molecular cloud complex. The H$^{13}$CO$^{+}$ and the NH$_{2}$D lines covered by the ALMA observations at a $\sim$3$^{\prime\prime}$ angular resolution ($\sim$0.02 pc) reveal that the spatially unresolved non-thermal motions are predominantly subsonic and transonic, a condition analogous to that found in low-mass star-forming molecular clouds. The observed supersonic non-thermal velocity dispersions in massive star forming regions, often reported in the literature, might be significantly biased by poor spatial resolutions that broaden the observed line widths due to unresolved motions within the telescope beam. Our 3~mm continuum image resolves 49 dense cores, whose masses range from 0.17 to 14 $M_{\odot}$. The majority of them are resolved with multiple velocity components. Our analyses of these gas velocity components find an anti-correlation between the gas mass and the virial parameter. This implies that the more massive structures tend to be more gravitationally unstable. Finally, we find that the external pressure in the NGC 6334S cloud is important in confining these dense structures, and may play a role in the formation of dense cores, and subsequently, the embedded young stars.

preprint2020arXiv

Tidal Interaction between the UX Tauri Disk A/C System Revealed by ALMA

We present sensitive and high angular resolution ($\sim$0.2-0.3$''$) (sub)millimeter (230 and 345 GHz) continuum and CO(2$-$1)/CO(3$-$2) line archive observations of the disk star system in UX Tauri carried out with ALMA (The Atacama Large Millimeter/Submillimeter Array). These observations reveal the gas and dusty disk surrounding the young star UX Tauri A with a large signal-to-noise ratio ($>$400 in the continuum and $>$50 in the line), and for the first time is detected the molecular gas emission associated with the disk of UX Tauri C (with a size for the disk of $<$56 au). No (sub)millimeter continuum emission is detected at 5$σ$-level (0.2 mJy at 0.85 mm) associated with UX Tauri C. For the component UX Tauri C, we estimate a dust disk mass of $\leq$ 0.05 M$_\oplus$. Additionally, we report a strong tidal disk interaction between both disks UX Tauri A/C, separated 360 au in projected distance. The CO line observations reveal marked spiral arms in the disk of UX Tauri A and an extended redshifted stream of gas associated with the UX Tauri C disk. No spiral arms are observed in the dust continuum emission of UX Tauri A. Assuming a Keplerian rotation we estimate the enclosed masses (disk$+$star) from their radial velocities in 1.4 $\pm$ 0.6 M$_\odot$ for UX Tauri A, and 70 $\pm$ 30 / $\sin i$ Jupiter masses for UX Tauri C (the latter coincides with the mass upper limit value for a brown dwarf). The observational evidence presented here lead us to propose that UX Tauri C is having a close approach of a possible wide, evolving and eccentric orbit around the disk of UX Tauri A causing the formation of spiral arms and the stream of molecular gas falling towards UX Tauri C.

preprint2020arXiv

VLA ammonia observations of L1287: analysis of the "Guitar" core and two filaments

The present work aims at studying the dense gas of the molecular cloud LDN 1287 (L1287), which harbors a double FU Ori system, an energetic molecular outflow and a still-forming cluster of deeply embedded low-mass, young stellar objects, showing a high level of fragmentation. We present optical Halpha and [SII], and VLA NH$_3$ (1,1) and (2,2) observations with an angular resolution of ~3.5''. The observed NH$_3$ spectra have been analyzed with the HfS tool, fitting simultaneously three different velocity components. The NH$_3$ emission from L1287 comes from four different structures: a core associated with RNO 1, a guitar-shaped core (the "Guitar") and two interlaced filaments (the Blue and Red Filaments) roughly centered towards the binary FU Ori system RNO 1B/C and its associated cluster. Regarding the Guitar Core, there are clear signatures of gas infall onto a central mass that has been estimated to be ~2.1 $M_\odot$. Regarding the two filaments, they have radii ~0.03 pc, masses per unit length ~50 $M_\odot$ pc$^{-1}$, and are near isothermal equilibrium. A central cavity, probably related with the outflow, and also traced by the Halpha and [SII] emission, is identified, with several young stellar objects near its inner walls. Both filaments show clear signs of perturbation by the high-velocity gas of the outflows driven by one or several young stellar objects of the cluster. The Blue and Red filaments are coherent in velocity and have nearly subsonic gas motions, except at the position of the embedded sources. Velocity gradients across the Blue Filament can be interpreted either as infalling material onto the filament or rotation. Velocity gradients along the filaments are interpreted as infall motions towards a gravitational well at the intersection of both filaments.

preprint2016arXiv

A correlation between chemistry, polarization and dust properties in the Pipe Nebula starless core FeSt 1-457

Pre-stellar cores within molecular clouds provide the very initial conditions in which stars are formed. We use the IRAM 30m telescope and the PdBI to study the chemical and physical properties of the starless core FeSt 1-457 (Core 109), in the Pipe nebula. We fit the hyperfine structure of the N2H+(1-0) IRAM 30m data. This allow us to measure with high precision the velocity field, line widths and opacity and derive the excitation temperature and column density in the core. We use a modified Bonnor-Ebert sphere model adding a temperature gradient towards the center to fit the 1.2 mm continuum emission and visual extinction maps. Using this model, we estimate the abundances of the N2H+ and the rest of molecular lines detected in the 30 GHz wide line survey performed at 3 mm with IRAM 30m using ARTIST software. The core presents a rich chemistry with emission from early (C3H2, HCN, CS) and late-time molecules (e.g., N2H+), with a clear chemical spatial differentiation for nitrogen, oxygen and sulphurated molecules. For most of the molecules detected (HCN, HCO+, CH3OH, CS, SO, 13CO and C18O), abundances are best fitted with three values, presenting a clear decrease of abundance of at least 1 or 2 orders of magnitude towards the center of the core. The Bonnor-Ebert analysis indicates the core is gravitationally unstable and the magnetic field is not strong enough to avoid the collapse. Depletion of molecules onto the dust grains occurs at the interior of the core, where dust grain growth and dust depolarization also occurs. This suggests that these properties may be related. On the other hand, some molecules exhibit asymmetries in their integrated emission maps, which appear to be correlated with a previously reported submillimetre polarization asymmetry. These asymmetries could be due to a stronger interstellar radiation field in the western side of the core.

preprint2016arXiv

Dense core properties in the Infrared Dark cloud G14.225-0.506 revealed by ALMA

We have performed a dense core survey toward the Infrared Dark Cloud G14.225-0.506 at 3 mm continuum emission with the Atacama Large Millimeter/Submillimeter Array (ALMA). This survey covers the two hub-filament systems with an angular resolution of $\sim3$\arcsec ($\sim0.03$ pc). We identified 48 dense cores. Twenty out of the 48 cores are protostellar due to their association with young stellar objects (YSOs) and/or X-ray point-sources, while the other 28 cores are likely prestellar and unrelated with known IR or X-ray emission. Using APEX 870 $μ$m continuum emission, we also identified the 18 clumps hosting these cores. Through virial analysis using the ALMA N$_2$H$^+$ and VLA/Effelsberg NH$_3$ molecular line data, we found a decreasing trend in the virial parameter with decreasing scales from filaments to clumps, and then to cores. The virial parameters of $0.1-1.3$ in cores, indicate that cores are likely undergoing dynamical collapse. The cumulative Core Mass Function (CMF) for the prestellar cores candidates has a power law index of $α=1.6$, with masses ranging from 1.5 to 22 $M_\odot$. We find no massive prestellar or protostellar cores. Previous studies suggest that massive O-tpye stars have not been produced yet in this region. Therefore, high-mass stars should be formed in the prestellar cores by accreting a significant amount of gas from the surrounding medium. Another possibility is that low-mass YSOs become massive by accreting from their parent cores that are fed by filaments. These two possibilities might be consistent with the scenario of global hierarchical collapse.

preprint2016arXiv

Phosphorus-bearing molecules in massive dense cores

Phosphorus is a crucial element for the development of life, but so far P-bearing molecules have been detected only in a few astrophysical objects, hence its interstellar chemistry is almost totally unknown. Here we show new detections of phosphorus nitride in a sample of dense cores in different evolutionary stages of the intermediate- and high-mass star formation process: starless, with protostellar objects, and with ultracompact HII regions. All detected PN line widths are smaller than ~5 km/s , and they arise from regions associated with kinetic temperatures smaller than 100 K. Because the few previous detections reported in the literature are associated with warmer and more turbulent sources, the results of this work show that PN can arise from relatively quiescent and cold gas. This information is challenging for theoretical models that invoke either high desorption temperatures or grain sputtering from shocks to release phosphorus into the gas phase. Derived column densities are of the order of 10^{11-12} cm^{-2}, marginally lower than the values derived in the few high-mass star forming regions detected so far. New constraints on the column density upper limit of phosphorus monoxide, the fundamental unit of biologically relevant molecules, are also given.

preprint2015arXiv

A submillimeter search for pre- and proto-brown dwarfs in Chamaeleon II

Context. Chamaeleon II molecular cloud is an active star forming region that offers an excellent opportunity for studying the formation of brown dwarfs in the southern hemisphere. Aims. Our aims are to identify a population of pre- and proto- brown dwarfs (5 sigma mass limit threshold of ~0.015 Msun) and provide information on the formation mechanisms of substellar objects. Methods. We performed high sensitivity observations at 870 microns using the LABOCA bolometer at the APEX telescope towards an active star forming region in Chamaeleon II. The data are complemented with an extensive multiwavelength catalogue of sources from the optical to the far-infrared to study the nature of the LABOCA detections. Results. We detect fifteen cores at 870 microns, and eleven of them show masses in the substellar regime. The most intense objects in the surveyed field correspond to the submillimeter counterparts of the well known young stellar objects DK Cha and IRAS 12500-7658. We identify a possible proto-brown dwarf candidate (ChaII-APEX-L) with IRAC emission at 3.6 and 4.5 microns. Conclusions. Our analysis indicates that most of the spatially resolved cores are transient, and that the point-like starless cores in the sub-stellar regime (with masses between 0.016 Msun and 0.066 Msun) could be pre-brown dwarfs cores gravitationally unstable if they have radii smaller than 220 AU to 907 AU (1.2" to 5" at 178 pc) respectively for different masses. ALMA observations will be the key to reveal the energetic state of these pre-brown dwarfs candidates.

preprint2015arXiv

Gravity or turbulence? -III. Evidence of pure thermal Jeans fragmentation at ~0.1 pc scale

We combine previously published interferometric and single-dish data of relatively nearby massive dense cores that are actively forming stars to test whether their `fragmentation level' is controlled by turbulent or thermal support. We find no clear correlation between the fragmentation level and velocity dispersion, nor between the observed number of fragments and the number of fragments expected when the gravitationally unstable mass is calculated including various prescriptions for `turbulent support'. On the other hand, the best correlation is found for the case of pure thermal Jeans fragmentation, for which we infer a core formation efficiency around 13 per cent, consistent with previous works. We conclude that the dominant factor determining the fragmentation level of star-forming massive dense cores at 0.1 pc scale seems to be thermal Jeans fragmentation.

preprint2015arXiv

Ongoing star formation in the proto-cluster IRAS 22134+5834

IRAS 22134+5834 was observed in the centimeter with (E)VLA, 3~mm with CARMA, 2~mm with PdBI, and 1.3~mm with SMA, to study the continuum emission as well as the molecular lines, that trace different physical conditions of the gas to study the influence of massive YSOs on nearby starless cores, and the possible implications in the clustered star formation process. The multi-wavelength centimeter continuum observations revealed two radio sources within the cluster, VLA1 and VLA2. VLA1 is considered to be an optically thin UCHII region with a size of 0.01~pc and sits at the edge of the near-infrared (NIR) cluster. The flux of ionizing photons of the VLA1 corresponds to a B1 ZAMS star. VLA2 is associated with an infrared point source and has a negative spectral index. We resolved six millimeter continuum cores at 2~mm, MM2 is associated with the UCHII region VLA1, and other dense cores are distributed around the UCH{\sc ii} region. Two high-mass starless clumps (HMSC), HMSC-E (east) and HMSC-W (west), are detected around the NIR cluster with N$_2$H$^+$(1--0) and NH$_3$ emission, and show different physical and chemical properties. Two N$_2$D$^+$ cores are detected on an NH$_3$ filament close to the UCHII region, with a projected separation of $\sim$8000~AU at the assumed distance of 2.6~kpc. The kinematic properties of the molecular line emission confirm the expansion of the UCHII region and that the molecular cloud around the near infrared (NIR) cluster is also expanding. Our multi-wavelength study has revealed different generations of star formation in IRAS 22134+5834. The formed intermediate- to massive stars show strong impact on nearby starless clumps. We propose that while the stellar wind from the UCHII region and the NIR cluster drives the large scale bubble, the starless clumps and HMPOs formed at the edge of the cluster.

preprint2015arXiv

Origin and kinematics of the eruptive flow from XZ Tau revealed by ALMA

We present high angular resolution ($\sim$0.94$"$) $^{12}$CO(1-0) Atacama Large Millimeter/Submillimeter Array (ALMA) observations obtained during the 2014 long baseline campaign from the eruptive bipolar flow from the multiple XZ Tau stellar system discovered by the Hubble Space Telescope (HST). These observations reveal, for the first time, the kinematics of the molecular flow. The kinematics of the different ejections close to XZ Tau reveal a rotating and expanding structure with a southeast-northwest velocity gradient. The youngest eruptive bubbles unveiled in the optical HST images are inside of this molecular expanding structure. Additionally, we report a very compact and collimated bipolar outflow emanating from XZ Tau A, which indicates that the eruptive outflow is indeed originating from this object. The mass (3 $\times$ 10$^{-7}$ $M_\odot$) and energetics (E$_{kin}$ $=$ 3 $\times$ 10$^{37}$ ergs) for the collimated outflow are comparable with those found in molecular outflows associated with young brown dwarfs.

preprint2015arXiv

Pre- and Post-burst Radio Observations of the Class 0 Protostar HOPS 383 in Orion

There is increasing evidence that episodic accretion is a common phenomenon in Young Stellar Objects (YSOs). Recently, the source HOPS 383 in Orion was reported to have a $\times 35$ mid-infrared -- and bolometric -- luminosity increase between 2004 and 2008, constituting the first clear example of a class 0 YSO (a protostar) with a large accretion burst. The usual assumption that in YSOs accretion and ejection follow each other in time needs to be tested. Radio jets at centimeter wavelengths are often the only way of tracing the jets from embedded protostars. We searched the Very Large Array archive for the available observations of the radio counterpart of HOPS 383. The data show that the radio flux of HOPS 383 varies only mildly from January 1998 to December 2014, staying at the level of $\sim 200$ to 300 $μ$Jy in the X band ($\sim 9$ GHz), with a typical uncertainty of 10 to 20 $μ$Jy in each measurement. We interpret the absence of a radio burst as suggesting that accretion and ejection enhancements do not follow each other in time, at least not within timescales shorter than a few years. Time monitoring of more objects and specific predictions from simulations are needed to clarify the details of the connection betwen accretion and jets/winds in YSOs.

preprint2015arXiv

Radio monitoring of the periodically variable IR source LRLL 54361: No direct correlation between the radio and IR emissions

LRLL 54361 is an infrared source located in the star forming region IC 348 SW. Remarkably, its infrared luminosity increases by a factor of 10 during roughly one week every 25.34 days. To understand the origin of these remarkable periodic variations, we obtained sensitive 3.3 cm JVLA radio continuum observations of LRLL 54361 and its surroundings in six different epochs: three of them during the IR-on state and three during the IR-off state. The radio source associated with LRLL 54361 remained steady and did not show a correlation with the IR variations. We suggest that the IR is tracing the results of fast (with a timescale of days) pulsed accretion from an unseen binary companion, while the radio traces an ionized outflow with an extent of $\sim$100 AU that smooths out the variability over a period of order a year. The average flux density measured in these 2014 observations, 27$\pm$5 $μ$Jy, is about a factor of two less than that measured about 1.5 years before, $53\pm$11 $μ$Jy, suggesting that variability in the radio is present, but over larger timescales than in the IR. We discuss other sources in the field, in particular two infrared/X-ray stars that show rapidly varying gyrosynchrotron emission.

preprint2014arXiv

A CO survey in planet-forming disks: characterizing the gas content in the epoch of planet formation

We carried out a 12CO(3-2) survey of 52 southern stars with a wide range of IR excesses (LIR/L*) using the single dish telescopes APEX and ASTE. The main aims were (1) to characterize the evolution of molecular gas in circumstellar disks using LIR/L* values as a proxy of disk dust evolution, and (2) to identify new gas-rich disk systems suitable for detailed study with ALMA. About 60% of the sample (31 systems) have LIR/L* > 0.01 typical of T-Tauri or Herbig AeBe stars, and the rest (21 systems) have LIR/L* < 0.01 typical of debris disks. We detect CO(3-2) emission from 20 systems, and 18 (90%) of these have LIR/L* > 0.01. However, the spectra of only four of the newly detected systems appear free of contamination from background or foreground emission from molecular clouds. These include the early-type stars HD 104237 (A4/5V, 116 pc) and HD 98922 (A2 III, 507 pc, as determined in this work), where our observations reveal the presence of CO-rich circumstellar disks for the first time. Of the other detected sources, many could harbor gaseous circumstellar disks, but our data are inconclusive. For these two newly discovered gas-rich disks, we present radiative transfer models that simultaneously reproduce their spectral energy distributions and the 12CO(3-2) line profiles. For both of these systems, the data are fit well by geometrically flat disks, placing them in the small class of non-flaring disks with significant molecular gas reservoirs.

preprint2014arXiv

A Luminous Blue Variable Star Interacting with a Nearby Infrared Dark Cloud

G79.29+0.46 is a nebula created by a Luminous Blue Variable (LBV) star candidate characterized by two almost circular concentric shells. In order to investigate whether the shells are interacting with the infrared dark cloud (IRDC) G79.3+0.3 located at the southwestern border of the inner shell, we conducted Jansky Very Large Array observations of NH3(1,1), (2,2) and c-C3H2, and combined them with previous Effelsberg data. The overall NH3 emission consists of one main clump, named G79A, elongated following the shape of the IRDC, plus two fainter and smaller cores to the north, which spatially match the inner infrared shell. We analysed the NH3 spectra at each position with detected emission and inferred linewidth, rotational temperature, column density and abundance maps, and find that: i) the linewidth of NH3(1,1) in the northern cores is 0.5 km/s, slightly larger than in their surroundings; ii) the NH3 abundance is enhanced by almost one order of magnitude towards the northwestern side of G79A; iii) there is one `hot slab' at the interface between the inner infrared shell and the NH3 peak of G79A; iv) the western and southern edges of G79A present chemical differentiation, with c-C3H2 tracing more external layers than NH3, similar to what is found in PDRs. Overall, the kinematics and physical conditions of G79A are consistent with both shock-induced and UV radiation-induced chemistry driven by the LBV star. Therefore, the IRDC is not likely associated with the star-forming region DR15, but located farther away, near G79.29+0.46 at 1.4 kpc.

preprint2014arXiv

A spider-like outflow in Barnard 5 - IRS 1: The transition from a collimated jet to a wide-angle outflow?

We present line and continuum observations made with the Submillimeter Array (SMA) of the young stellar object Barnard 5 - IRS1 located in the Perseus molecular cloud. Our $^{12}$CO(2-1) line observations resolve the high-velocity bipolar northeast-southwest outflow associated with this source. We find that the outflowing gas shows different structures at three different velocity regimes, in both lobes, resulting in a spider-like morphology. In addition to the low-velocity, cone-like (wide-angle) lobes that have previously been observed, we report the presence of intermediate-velocity parabolic shells emerging very close to the Class I protostar, as well as high velocity molecular bullets that appear to be associated to the optical/IR jet emanating from this source. These compact high-velocity features reach radial velocities of about 50 km s$^{-1}$ away from the cloud velocity. We interpret the peculiar spider-like morphology is a result of the molecular material being entrained by a wind with both a collimated jet-like component and a wide-angle component. We suggest the outflow is in a transitional evolutionary phase between a mostly jet-driven flow and an outflow in which the entrainment is dominated by the wide-angle wind component. We also detect 1300 $μ$m continuum emission at the position of the protostar, which likely arises from the dusty envelope and disk surrounding the protostar. Finally, we report the detection of $^{13}$CO(2-1) and SO(6$_5$-5$_4$) emission arising from the outflow and the location of the young stellar object.

preprint2014arXiv

ALMA reveals a candidate hot and compact disk around the O-type protostar IRAS 16547$-$4247

We present high angular resolution ($\sim$ 0.3$"$) submillimeter continuum (0.85 mm) and line observations of the O-type protostar IRAS 16547$-$4247 carried out with the Atacama Large Millimeter/Submillimeter Array (ALMA). In the 0.85 mm continuum band, the observations revealed two compact sources (with a separation of 2$"$), one of them associated with IRAS 16547$-$4247, and the other one to the west. Both sources are well resolved angularly, revealing a clumpy structure. On the other hand, the line observations revealed a rich variety of molecular species related to both continuum sources. In particular, we found a large number of S-bearing molecules, such as the rare molecule methyl mercaptan (CH$_3$SH). At scales larger than 10,000 AU, molecules (e.g., SO$_2$ or OCS) mostly with low excitation temperatures in the upper states (E$_k$ $\lesssim$ 300 K) are present in both millimeter continuum sources, and show a southeast-northwest velocity gradient of 7 km s$^{-1}$ over 3$"$ (165 km s$^{-1}$ pc$^{-1}$). We suggest that this gradient probably is produced by the thermal (free-free) jet emerging from this object with a similar orientation at the base. At much smaller scales (about 1000 AU), molecules with high excitation temperatures (E$_k$ $\gtrsim$ 500 K) are tracing a rotating structure elongated perpendicular to the orientation of the thermal jet, which we interpret as a candidate disk surrounding IRAS 16547$-$4247. The dynamical mass corresponding to the velocity gradient of the candidate to disk is about 20 M$_\odot$, which is consistent with the bolometric luminosity of IRAS 16547$-$4247.

preprint2014arXiv

Deuteration and evolution in the massive star formation process: the role of surface chemistry

An ever growing number of observational and theoretical evidence suggests that the deuterated fraction (column density ratio between a species containing D and its hydrogenated counterpart, Dfrac) is an evolutionary indicator both in the low- and the high-mass star formation process. However, the role of surface chemistry in these studies has not been quantified from an observational point of view. In order to compare how the deuterated fractions of species formed only in the gas and partially or uniquely on grain surfaces evolve with time, we observed rotational transitions of CH3OH, 13CH3OH, CH2DOH, CH3OD at 3 and 1.3~mm, and of NH2D at 3~mm with the IRAM-30m telescope, and the inversion transitions (1,1) and (2,2) of NH3 with the GBT, towards most of the cores already observed by Fontani et al.~(2011, 2014) in N2H+, N2D+, HNC, DNC. NH2D is detected in all but two cores, regardless of the evolutionary stage. Dfrac(NH3) is on average above 0.1, and does not change significantly from the earliest to the most evolved phases, although the highest average value is found in the protostellar phase (~0.3). Few lines of CH2DOH and CH3OD are clearly detected, and only towards protostellar cores or externally heated starless cores. This work clearly confirms an expected different evolutionary trend of the species formed exclusively in the gas (N2D+ and N2H+) and those formed partially (NH2D and NH3) or totally (CH2DOH and CH3OH) on grain mantles. The study also reinforces the idea that Dfrac(N2H+) is the best tracer of massive starless cores, while high values of Dfrac(CH3OH) seem rather good tracers of the early protostellar phases, at which the evaporation/sputtering of the grain mantles is most efficient.

preprint2014arXiv

Extended warm gas in Orion KL as probed by methyl cyanide

In order to study the temperature distribution of the extended gas within the Orion Kleinmann-Low nebula, we have mapped the emission by methyl cyanide (CH3CN) in its J=6_K-5_K, J=12_K-11_K, J=13_K-12_K, and J=14_K-13_K transitions at an average angular resolution of ~10 arcsec (22 arcsec for the 6_K-5_K lines), as part of a new 2D line survey of this region using the IRAM 30m telescope. These fully sampled maps show extended emission from warm gas to the northeast of IRc2 and the distinct kinematic signatures of the hot core and compact ridge source components. We have constructed population diagrams for the four sets of K-ladder emission lines at each position in the maps and have derived rotational excitation temperatures and total beam-averaged column densities from the fitted slopes. In addition, we have fitted LVG model spectra to the observations to determine best-fit physical parameters at each map position, yielding the distribution of kinetic temperatures across the region. The resulting temperature maps reveal a region of hot (T > 350 K) material surrounding the northeastern edge of the hot core, whereas the column density distribution is more uniform and peaks near the position of IRc2. We attribute this region of hot gas to shock heating caused by the impact of outflowing material from active star formation in the region, as indicated by the presence of broad CH3CN lines. This scenario is consistent with predictions from C-shock chemical models that suggest that gas-phase methyl cyanide survives in the post-shock gas and can be somewhat enhanced due to sputtering of grain mantles in the passing shock front.

preprint2014arXiv

Fragmentation of massive dense cores down to ~1000 AU: Relation between fragmentation and density structure

In order to shed light on the main physical processes controlling fragmentation of massive dense cores, we present a uniform study of the density structure of 19 massive dense cores, selected to be at similar evolutionary stages, for which their relative fragmentation level was assessed in a previous work. We inferred the density structure of the 19 cores through a simultaneous fit of the radial intensity profiles at 450 and 850 micron (or 1.2 mm in two cases) and the Spectral Energy Distribution, assuming spherical symmetry and that the density and temperature of the cores decrease with radius following power-laws. We find a weak (inverse) trend of fragmentation level and density power-law index, with steeper density profiles tending to show lower fragmentation, and vice versa. In addition, we find a trend of fragmentation increasing with density within a given radius, which arises from a combination of flat density profile and high central density and is consistent with Jeans fragmentation. We considered the effects of rotational-to-gravitational energy ratio, non-thermal velocity dispersion, and turbulence mode on the density structure of the cores, and found that compressive turbulence seems to yield higher central densities. Finally, a possible explanation for the origin of cores with concentrated density profiles, which are the cores showing no fragmentation, could be related with a strong magnetic field, consistent with the outcome of radiation magnetohydrodynamic simulations.

preprint2014arXiv

IC348-SMM2E: a Class 0 proto-brown dwarf candidate forming as a scaled-down version of low-mass stars

We report on Submillimeter Array observations of the 870 micron continuum and CO(3-2), 13CO(2-1) and C18O(2-1) line emission of a faint object, SMM2E, near the driving source of the HH797 outflow in the IC348 cluster. The continuum emission shows an unresolved source for which we estimate a mass of gas and dust of 30 Mjup, and the CO(3-2) line reveals a compact bipolar outflow centred on SMM2E, and barely seen also in 13CO(2-1). In addition, C18O(2-1) emission reveals hints of a possible rotating envelope/disk perpendicular to the outflow, for which we infer a dynamical mass of ~16 Mjup. In order to further constrain the accreted mass of the object, we gathered data from Spitzer, Herschel, and new and archive submillimetre observations, and built the Spectral Energy Distribution (SED). The SED can be fitted with one single modified black-body from 70 micron down to 2.1 cm, using a dust temperature of ~24 K, a dust emissivity index of 0.8, and an envelope mass of ~35 Mjup. The bolometric luminosity is 0.10 Lsun, and the bolometric temperature is 35 K. Thus, SMM2E is comparable to the known Class 0 objects in the stellar domain. An estimate of the final mass indicates that SMM2E will most likely remain substellar, and the SMM2E outflow force matches the trend with luminosity known for young stellar objects. Thus, SMM2E constitutes an excellent example of a Class 0 proto-brown dwarf candidate which forms as a scaled-down version of low-mass stars. Finally, SMM2E seems to be part of a wide (~2400 AU) multiple system of Class 0 sources.

preprint2014arXiv

The hot core towards the intermediate mass protostar NGC7129 FIRS 2: Chemical similarities with Orion KL

NGC 7129 FIRS 2 (hereafter FIRS 2) is an intermediate-mass (2 to 8 Msun) protostar located at a distance of 1250 pc. High spatial resolution observations are required to resolve the hot core at its center. We present a molecular survey from 218200 MHz to 221800 MHz carried out with the IRAM Plateau de Bure Interferometer. These observations were complemented with a long integration single-dish spectrum taken with the IRAM 30m telescope. We used a Local Thermodynamic Equilibrium (LTE) single temperature code to model the whole dataset. The interferometric spectrum is crowded with a total of ~300 lines from which a few dozens remain unidentified yet. The spectrum has been modeled with a total of 20 species and their isomers, isotopologues and deuterated compounds. Complex molecules like methyl formate (CH3OCHO), ethanol (CH3CH2OH), glycolaldehyde (CH2OHCHO), acetone (CH3COCH3), dimethyl ether (CH3OCH3), ethyl cyanide (CH3CH2CN) and the aGg' conformer of ethylene glycol (aGg'-(CH2OH)_2) are among the detected species. The detection of vibrationally excited lines of CH3CN, CH3OCHO, CH3OH, OCS, HC3N and CH3CHO proves the existence of gas and dust at high temperatures. In fact, the gas kinetic temperature estimated from the vibrational lines of CH3CN, ~405 K, is similar to that measured in massive hot cores. Our data allow an extensive comparison of the chemistry in FIRS~2 and the Orion hot core. We find a quite similar chemistry in FIRS 2 and Orion. Most of the studied fractional molecular abundances agree within a factor of 5. Larger differences are only found for the deuterated compounds D2CO and CH2DOH and a few molecules (CH3CH2CN, SO2, HNCO and CH3CHO). Since the physical conditions are similar in both hot cores, only different initial conditions (warmer pre-collapse phase in the case of Orion) and/or different crossing time of the gas in the hot core can explain this behavior.

preprint2013arXiv

A combined IRAM and Herschel/HIFI study of cyano(di)acetylene in Orion KL: tentative detection of DC3N

We present a study of cyanoacetylene (HC3N) and cyanodiacetylene (HC5N) in Orion KL, through observations from two line surveys performed with the IRAM 30m telescope and the HIFI instrument on board the Herschel telescope. The frequency ranges covered are 80-280 GHz and 480-1906 GHz. We model the observed lines of HC3N, HC5N, their isotopologues (including DC3N), and vibrational modes, using a non-LTE radiative transfer code. To investigate the chemical origin of HC3N and DC3N in Orion KL, we use a time-dependent chemical model. We detect 40 lines of the ground state of HC3N and 68 lines of its 13C isotopologues. We also detect 297 lines of six vibrational modes of this molecule (nu_7, 2nu_7, 3nu_7, nu_6, nu_5, and nu_6+nu_7) and 35 rotational lines of the ground state of HC5N. We report the first tentative detection of DC3N in a giant molecular cloud with a DC3N/HC3N abundance ratio of 0.015. We provide column densities and isotopic and molecular abundances. We also perform a 2x2" map around Orion IRc2 and we present maps of HC3N lines and maps of lines of the HC3N vibrational modes nu_6 and nu_7. In addition, a comparison of our results for HC3N with those in other clouds allows us to derive correlations between the column density, the FWHM, the mass, and the luminosity of the clouds. The high column densities of HC3N obtained in the hot core, make this molecule an excellent tracer of hot and dense gas. In addition, the large frequency range covered reveals the need to consider a temperature and density gradient in the hot core in order to obtain better line fits. The high D/H ratio (comparable to that obtained in cold clouds) that we derive suggests a deuterium enrichment. Our chemical models indicate that the possible deuterated HC3N present in Orion KL is formed during the gas-phase. This fact provides new hints concerning the processes leading to deuteration.

preprint2013arXiv

A line confusion-limited millimeter survey of Orion KL. III. Sulfur oxide species

We present a study of the sulfur-bearing species detected in a line confusion-limited survey towards Orion KL performed with the IRAM 30m telescope in the range 80-281 GHz. The study is part of an analysis of the line survey divided into families of molecules. Our aim is to derive accurate physical conditions and molecular abundances in the different components of Orion KL from observed SO and SO2 lines. First we assumed LTE conditions obtain rotational temperatures. We then used a radiative transfer model, assuming either LVG or LTE excitation to derive column densities of these molecules in the different components of Orion KL. We have detected 68 lines of SO, 34SO, 33SO, and S18O and 653 lines of SO2, 34SO2, 33SO2, SO18O and SO2 v2=1. We provide column densities for all of them and also upper limits for the column densities of S17O, 36SO, 34S18O, SO17O and 34SO2 v2=1 and for several undetected sulfur-bearing species. In addition, we present 2'x2' maps around Orion IRc2 of SO2 transitions with energies from 19 to 131 K and also maps with four transitions of SO, 34SO and 34SO2. We observe an elongation of the gas along the NE-SW direction. An unexpected emission peak appears at 20.5 km/s in most lines of SO and SO2. A study of the spatial distribution of this emission feature shows that it is a new component ~5" in diameter, which lies ~4" west of IRc2. We suggest the emission from this feature is related to shocks associated to the BN object. The highest column densities for SO and SO2 are found in the high-velocity plateau (a region dominated by shocks) and in the hot core. These values are up to three orders of magnitude higher than the results for the ridge components. We also find high column densities for their isotopologues in both components. Therefore, we conclude that SO and SO2 are good tracers, not only of regions affected by shocks, but also of regions with warm dense gas.

preprint2013arXiv

Deciphering the ionized gas content in the massive star forming complex G75.78+0.34

We present sub-arcsecond observations toward the massive star forming region G75.78+0.34. We used the Very Large Array to study the centimeter continuum and H2O and CH3OH maser emission, and the Owens Valley Radio Observatory and Submillimeter Array to study the millimeter continuum and recombination lines (H40alpha and H30alpha). We found radio continuum emission at all wavelengths, coming from three components: (1) a cometary ultracompact (UC) HII region with an electron density 3.7x10^4 cm^{-3}, excited by a B0 type star, and with no associated dust emission; (2) an almost unresolved UCHII region (EAST), located 6 arcsec to the east of the cometary UCHII region, with an electron density 1.3x10^5 cm^{-3}, and associated with a compact dust clump detected at millimeter and mid-infrared wavelengths; and (3) a compact source (CORE), located 2 arcsec to the southwest of the cometary arc, with a flux density increasing with frequency, and embedded in a dust condensation of 30 Msun. The CORE source is resolved into two compact and unresolved sources which can be well-fit by two homogeneous hypercompact HII regions each one photo-ionized by a B0.5 ZAMS star, or by free-free radiation from shock-ionized gas resulting from the interaction of a jet/outflow system with the surrounding environment. The spatial distribution and kinematics of water masers close to the CORE-N and S sources, together with excess emission at 4.5 mum and the detected dust emission, suggest that the CORE source is a massive protostar driving a jet/outflow.

preprint2013arXiv

On the origin of the molecular outflows in IRAS 16293-2422

We present CO 3-2, SiO 8-7, C34S 7-6, and 878 mum dust continuum subarcsecond angular resolution observations with the SMA toward IRAS 16293-2422 (I16293). The C34S emission traces the 878 mum dust continuum well, and clearly shows a smooth velocity gradient along the major axis of component I16293A. The CO shows emission at moderate high velocities arising from two bipolar outflows, which appear to be perpendicular with respect to each other. The high sensitivity and higher angular resolution of these observations allows us to pinpoint well the origin of these two outflows at the center of component I16293A. Interestingly, the most compact outflow appears to point toward I16293B. Our data show that the previously reported monopolar blueshifted CO outflow associated with component I16293B seems to be part of the compact outflow arising from component I16293A. In addition, the SiO emission is also tracing this compact outflow: on one hand, the SiO emission appears to have a jet-like morphology along the southern redshifted lobe; on the other hand, the SiO emission associated with the blueshifted northern lobe traces a well defined arc on the border of component I16293B facing I16293A. The blueshifted CO lobe of the compact outflow splits into two lobes around the position of this SiO arc. All these results lead us to propose that the compact outflow from component I16293A is impacting on the circumstellar gas around component I16293B, possibly being diverged as a consequence of the interaction.

preprint2013arXiv

Properties of dense cores in clustered massive star-forming regions at high angular resolution

We aim at characterising dense cores in the clustered environments associated with massive star-forming regions. For this, we present an uniform analysis of VLA NH3(1,1) and (2,2) observations towards a sample of 15 massive star-forming regions, where we identify a total of 73 cores, classify them as protostellar, quiescent starless, or perturbed starless, and derive some physical properties. The average sizes and ammonia column densities are 0.06 pc and 10^15 cm^-2, respectively, with no significant differences between the starless and protostellar cores, while the linewidth and rotational temperature of quiescent starless cores are smaller, 1.0 km/s and 16 K, than those of protostellar (1.8 km/s, 21 K), and perturbed starless (1.4 km/s, 19 K) cores. Such linewidths and temperatures for these quiescent starless cores in the surroundings of massive stars are still significantly larger than the typical values measured in starless cores of low-mass star-forming regions, implying an important non-thermal component. We confirm at high angular resolutions the correlations previously found with single-dish telescopes between the linewidth, the temperature of the cores, and the bolometric luminosity. In addition, we find a correlation between the temperature of each core and the incident flux from the most massive star in the cluster, suggesting that the large temperatures measured in the starless cores of our sample could be due to heating from the nearby massive star. A simple virial equilibrium analysis seems to suggest a scenario of a self-similar, self-graviting, turbulent, virialised hierarchy of structures from clumps (0.1-10 pc) to cores (0.05 pc). A closer inspection of the dynamical state taking into account external pressure effects, reveal that relatively strong magnetic field support may be needed to stabilise the cores, or that they are unstable and thus on the verge of collapse.

preprint2013arXiv

The outburst decay of the low magnetic field magnetar SGR 0418+5729

We report on the long term X-ray monitoring of the outburst decay of the low magnetic field magnetar SGR 0418+5729, using all the available X-ray data obtained with RXTE, SWIFT, Chandra, and XMM-Newton observations, from the discovery of the source in June 2009, up to August 2012. The timing analysis allowed us to obtain the first measurement of the period derivative of SGR 0418+5729: \dot{P}=4(1)x10^{-15} s/s, significant at ~3.5 sigma confidence level. This leads to a surface dipolar magnetic field of B_dip ~6x 10^{12} G. This measurement confirms SGR 0418+5729 as the lowest magnetic field magnetar. Following the flux and spectral evolution from the beginning of the outburst up to ~1200 days, we observe a gradual cooling of the tiny hot spot responsible for the X-ray emission, from a temperature of ~0.9 to 0.3 keV. Simultaneously, the X-ray flux decreased by about 3 orders of magnitude: from about 1.4x10^{-11} to 1.2x10^{-14} erg/s/cm^2 . Deep radio, millimeter, optical and gamma-ray observations did not detect the source counterpart, implying stringent limits on its multi-band emission, as well as constraints on the presence of a fossil disk. By modeling the magneto-thermal secular evolution of SGR 0418+5729, we infer a realistic age of ~550 kyr, and a dipolar magnetic field at birth of ~10^{14} G. The outburst characteristics suggest the presence of a thin twisted bundle with a small heated spot at its base. The bundle untwisted in the first few months following the outburst, while the hot spot decreases in temperature and size. We estimate the outburst rate of low magnetic field magnetars to be about one per year per galaxy, and we briefly discuss the consequences of such result in several other astrophysical contexts.

preprint2012arXiv

A search for pre-substellar cores and proto-brown dwarf candidates in Taurus: multiwavelength analysis in the B213-L1495 clouds

In an attempt to study whether the formation of brown dwarfs (BDs) takes place as a scaled-down version of low-mass stars, we conducted IRAM30m/MAMBO-II observations at 1.2 mm in a sample of 12 proto-BD candidates selected from Spitzer/IRAC data in the B213-L1495 clouds in Taurus. Subsequent observations with the CSO at 350 micron, VLA at 3.6 and 6 cm, and IRAM30m/EMIR in the 12CO(1-0), 13CO(1-0), and N2H+(1-0) transitions were carried out toward the two most promising Spitzer/IRAC source(s), J042118 and J041757. J042118 is associated with a compact (<10 arcsec or <1400 AU) and faint source at 350 micron, while J041757 is associated with a partially resolved (~16 arcsec or ~2000 AU) and stronger source emitting at centimetre wavelengths with a flat spectral index. The corresponding masses of the dust condensations are ~1 and ~5 Mjup for J042118 and J041757, respectively. In addition, about 40 arcsec to the northeast of J041757 we detect a strong and extended submillimetre source, J041757-NE, which is not associated with NIR/FIR emission down to our detection limits, but is clearly detected in 13CO and N2H+ at ~7 km/s, and for which we estimated a total mass of ~100 Mjup, close to the mass required to be gravitationally bound. In summary, our observational strategy has allowed us to find in B213-L1495 two proto-BD candidates and one pre-substellar core candidate, whose properties seem to be consistent with a scaled-down version of low-mass stars.

preprint2012arXiv

Dense gas in IRAS 20343+4129: an ultracompact HII region caught in the act of creating a cavity

The intermediate- to high-mass star-forming region IRAS 20343+4129 is an excellent laboratory to study the influence of high- and intermediate-mass young stellar objects on nearby starless dense cores, and investigate for possible implications in the clustered star formation process. We present 3 mm observations of continuum and rotational transitions of several molecular species (C2H, c-C3H2, N2H+, NH2D) obtained with the Combined Array for Research in Millimetre-wave Astronomy, as well as 1.3 cm continuum and NH3 observations carried out with the Very Large Array, to reveal the properties of the dense gas. We confirm undoubtedly previous claims of an expanding cavity created by an ultracompact HII region associated with a young B2 zero-age main sequence (ZAMS) star. The dense gas surrounding the cavity is distributed in a filament that seems squeezed in between the cavity and a collimated outflow associated with an intermediate-mass protostar. We have identified 5 millimeter continuum condensations in the filament. All of them show column densities consistent with potentially being the birthplace of intermediate- to high-mass objects. These cores appear different from those observed in low-mass clustered environments in sereval observational aspects (kinematics, temperature, chemical gradients), indicating a strong influence of the most massive and evolved members of the protocluster. We suggest a possible scenario in which the B2 ZAMS star driving the cavity has compressed the surrounding gas, perturbed its properties and induced the star formation in its immediate surroundings.

preprint2012arXiv

Early stages of cluster formation: fragmentation of massive dense cores down to ~1000 AU

In order to study the fragmentation of massive dense cores, which constitute the cluster cradles, we observed with the PdBI in the most extended configuration the continuum at 1.3 mm and the CO(2-1) emission of four massive cores. We detect dust condensations down to ~0.3 Msun and separate millimeter sources down to 0.4" or ~1000 AU, comparable to the sensitivities and separations reached in optical/infrared studies of clusters. The CO(2-1) high angular resolution images reveal high-velocity knots usually aligned with previously known outflow directions. This, in combination with additional cores from the literature observed at similar mass sensitivity and spatial resolution, allowed us to build a sample of 18 protoclusters with luminosities spanning 3 orders of magnitude. Among the 18 regions, ~30% show no signs of fragmentation, while 50% split up into ~4 millimeter sources. We compiled a list of properties for the 18 massive dense cores, such as bolometric luminosity, total mass, and mean density, and found no correlation of any of these parameters with the fragmentation level. In order to investigate the combined effects of magnetic field, radiative feedback and turbulence in the fragmentation process, we compared our observations to radiation magneto-hydrodynamic simulations, and found that the low-fragmented regions are well reproduced in the magnetized core case, while the highly-fragmented regions are consistent with cores where turbulence dominates over the magnetic field. Overall, our study suggests that the fragmentation in massive dense cores could be determined by the initial magnetic field/turbulence balance in each particular core.

preprint2012arXiv

From dusty filaments to massive stars: The case of NGC 7538 S

We report on high-sensitivity and high-angular resolution archival Submillimeter Array (SMA) observations of the large ($\sim$15000 AU) putative circumstellar disk associated with the O-type protostar NGC 7538 S. Observations of the continuum resolve this putative circumstellar disk into five compact sources, with sizes $\sim$ 3000 AU and masses $\sim10 M_\odot$. This confirm the results of recent millimeter observations made with CARMA/BIMA towards this object. However, we find that from most of these compact sources eject collimated bipolar outflows, revealed by our silicon monoxide (SiO {\it J}=5-4) observations and confirm that these sources have a (proto)stellar nature. All outflows are perpendicular to the large and rotating dusty structure. We propose therefore that, rather than being a single massive circumstellar disk, NGC 7538 S could be instead a large and massive contracting or rotating filament that is fragmenting at scales of 0.1 to 0.01 pc to form several B-type stars, via the standard process involving outflows and disks. As in recent high spatial resolution studies of dusty filaments, our observations also suggest that thermal pressure does not seem to be sufficient to support the filament, so that either additional support needs to be invoked, or else the filament must be in the process of collapsing. An SPH numerical simulation of the formation of a molecular cloud by converging warm neutral medium flows produces contracting filaments whose dimensions and spacings between the stars forming within them, as well as their column densities, strongly resemble those observed in the filament reported here.

preprint2012arXiv

IRAS 19520+2759: a 10^5 Lsun massive young stellar object driving a collimated outflow

The theory of massive star formation currently suffers from a scarce observational base of massive young stellar objects to compare with. In this paper, we present OVRO CO(1-0), 13CO(1-0), C18O(1-0), and 2.6 mm continuum images of the infrared source IRAS 19520+2759 together with complementary single-dish observations of CS(1-0), obtained with the 34m antenna DSS-54 at the Madrid Deep Space Communications Complex, as well as archive images at different wavelengths. As a result from our work, IRAS 19520+2759, with a controversial nature in the past, is firmly established as a massive young stellar object associated with a strong and compact millimetre source and driving a collimated outflow. In addition, a second fainter millimetre source is discovered about 4 arcsec to the south, which is also driving an outflow. Furthermore, the two millimetre sources are associated with C18O clumps elongated perpendicularly to the outflows, which may be related to rotating toroids. The masses of gas and dust of the millimetre sources are estimated to be around 100 and 50 Msun. MM1, the dominant source at all wavelengths, with a total luminosity of (1-2)x10^5 Lsun at 9 kpc, is however not associated with 6 cm emission down to a rms noise level of 0.1 mJy. We propose that IRAS 19520+2759 could be an example of the recent theoretical prediction of `bloated' or `swollen' star, i.e., a massive young stellar object whose radius has increased due to effects of accretion at a high-mass accretion rate.

preprint2012arXiv

Temperature and kinematics of protoclusters with intermediate and high-mass stars: the case of IRAS 05345+3157

We have mapped at small spatial scales the temperature and the velocity field in the protocluster associated with IRAS 05345+3157, which contains both intermediate-/high-mass protostellar candidates and starless condensations, and is thus an excellent location to investigate the role of massive protostars on protocluster evolution. We observed the ammonia (1,1) and (2,2) inversion transitions with the VLA. Ammonia is the best thermometer for dense and cold gas, and the observed transitions have critical densities able to trace the kinematics of the intracluster gaseous medium. The ammonia emission is extended and distributed in two filamentary structures. The starless condensations are colder than the star-forming cores, but the gas temperature across the whole protocluster is higher (by a factor of ~1.3-1.5) than that measured typically in both infrared dark clouds and low-mass protoclusters. The non-thermal contribution to the observed line broadening is at least a factor of 2 larger than the expected thermal broadening even in starless condensations, contrary to the close-to-thermal line widths measured in low-mass quiescent dense cores. The NH3-to-N2H+ abundance ratio is greatly enhanced (a factor of 10) in the pre--stellar core candidates, probably due to freeze-out of most molecular species heavier than He. The more massive and evolved objects likely play a dominant role in the physical properties and kinematics of the protocluster. The high level of turbulence and the fact that the measured core masses are larger than the expected thermal Jeans masses indicate that turbulence likely was an important factor in the initial fragmentation of the parental clump.

preprint2012arXiv

Unveiling a network of parallel filaments in the Infrared Dark Cloud G14.225-0.506

We present the results of combined NH3(1,1) and (2,2) line emission observed with the Very Large Array and the Effelsberg 100m telescope of the Infrared Dark Cloud G14.225-0.506. The NH3 emission reveals a network of filaments constituting two hub-filament systems. Hubs are associated with gas of rotational temperature Trot \sim 25 K, non-thermal velocity dispersion ~1.1 km/s, and exhibit signs of star formation, while filaments appear to be more quiescent (Trot \sim 11 K, non-thermal velocity dispersion ~0.6 km/s). Filaments are parallel in projection and distributed mainly along two directions, at PA \sim 10 deg and 60 deg, and appear to be coherent in velocity. The averaged projected separation between adjacent filaments is between 0.5 pc and 1pc, and the mean width of filaments is 0.12 pc. Cores within filaments are separated by ~0.33 pc, which is consistent with the predicted fragmentation of an isothermal gas cylinder due to the 'sausage'-type instability. The network of parallel filaments observed in G14.225-0.506 is consistent with the gravitational instability of a thin gas layer threaded by magnetic fields. Overall, our data suggest that magnetic fields might play an important role in the alignment of filaments, and polarization measurements in the entire cloud would lend further support to this scenario.

preprint2011arXiv

Deuteration as an evolutionary tracer in massive-star formation

Theory predicts, and observations confirm, that the column density ratio of a molecule containing D to its counterpart containing H can be used as an evolutionary tracer in the low-mass star formation process. Since it remains unclear if the high-mass star formation process is a scaled-up version of the low-mass one, we investigated whether the relation between deuteration and evolution can be applied to the high-mass regime. With the IRAM-30m telescope, we observed rotational transitions of N2D+ and N2H+ and derived the deuterated fraction in 27 cores within massive star-forming regions understood to represent different evolutionary stages of the massive-star formation process. Results. Our results clearly indicate that the abundance of N2D+ is higher at the pre-stellar/cluster stage, then drops during the formation of the protostellar object(s) as in the low-mass regime, remaining relatively constant during the ultra-compact HII region phase. The objects with the highest fractional abundance of N2D+ are starless cores with properties very similar to typical pre-stellar cores of lower mass. The abundance of N2D+ is lower in objects with higher gas temperatures as in the low-mass case but does not seem to depend on gas turbulence. Our results indicate that the N2D+-to-N2H+ column density ratio can be used as an evolutionary indicator in both low- and high-mass star formation, and that the physical conditions influencing the abundance of deuterated species likely evolve similarly during the processes that lead to the formation of both low- and high-mass stars.

preprint2011arXiv

Intermediate-mass hot cores at 500 AU: disks or outflows?

Observations with the Plateau de Bure Interferometer in the most extended configuration toward two intermediate-mass star-forming regions, IRAS22198+6336 and AFGL5142, reveal the presence of several complex organic molecules at ~500 AU scales, confirming the presence of hot cores in both regions. The hot cores are not rich in CN-bearing molecules, as often seen in massive hot cores, and are mainly traced by CH3CH2OH, (CH2OH)2, CH3COCH3, and CH3OH, with additionally CH3CHO, CH3OD and HCOOD for IRAS22198+6336, and C6H, and O13CS for AFGL5142. The emission of complex molecules is resolved down to sizes of ~300 and ~600 AU, for IRAS22198+6336 and AFGL5142, respectively, and most likely is tracing protostellar disks rather than flattened envelopes or toroids as usually found. This is specially clear for the case of IRAS22198+6336, where we detect a velocity gradient for all the mapped molecules perpendicular to the most chemically rich outflow of the region, yielding a dynamic mass >4 Msun. As for AFGL5142, the hot core emission is resolved into two elongated cores separated 1800 AU. A detailed comparison of the complex molecule peaks to the new CO(2-1) data and H2O maser data from literature suggests that also for AFGL5142 the complex molecules are mainly associated with disks, except for a faint and extended molecular emission found to the west, which is possibly produced in the interface between one of the outflows and the dense surrounding gas.

preprint2010arXiv

IRAS 22198+6336: Discovery of an Intermediate-Mass Hot Core

We present new SMA and PdBI observations of the intermediate-mass object IRAS 22198+6336 in the millimeter continuum and in several molecular line transitions. The millimeter continuum emission reveals a strong and compact source with a mass of around 5 Msun and with properties of Class 0 objects. CO emission shows an outflow with a quadrupolar morphology centered on the position of the dust condensation. The CO outflow emission seems to come from two distinct outflows, one of them associated with SiO outflow emission. A large set of molecular lines has been detected toward a compact dense core clearly coincident with the compact millimeter source, and showing a velocity gradient perpendicular to the outflow traced by CO and SiO. The chemically rich spectrum and the rotational temperatures derived from CH$_3$CN and CH$_3$OH (100-150 K) indicate that IRAS 22198+6336 is harbouring one the few intermediate-mass hot cores known at present.

preprint2010arXiv

Multiple shells around G79.29+0.46 revealed from near-IR to millimeter data

Aiming to perform a study of the warm dust and gas in the luminous blue variable star G79.29+0.46 and its associated nebula, we present infrared Spitzer imaging and spectroscopy, and new CO J=2-->1 and 4-->3 maps obtained with the IRAM 30m radio telescope and with the Submillimeter Telescope, respectively. We have analyzed the nebula detecting multiple shells of dust and gas connected to the star. Using Infrared Spectrograph-Spitzer spectra, we have compared the properties of the central object, the nebula, and their surroundings. These spectra show a rich variety of solid-state features (amorphous silicates, polycyclic aromatic hydrocarbons, and CO2 ices) and narrow emission lines, superimposed on a thermal continuum. We have also analyzed the physical conditions of the nebula, which point to the existence of a photo-dissociation region.

preprint2010arXiv

N2H+ depletion in the massive protostellar cluster AFGL 5142

We aim at investigating with high angular resolution the NH3/N2H+ ratio toward the high-mass star-forming region AFGL 5142 in order to study whether this ratio behaves similarly to the low-mass case, for which the ratio decreases from starless cores to cores associated with YSOs. CARMA was used to observe the 3.2 mm continuum and N2H+(1-0) emission. We used NH3(1,1) and (2,2), HCO+(1-0) and H13CO+(1-0) data from the literature and we performed a time-dependent chemical modeling of the region. The 3.2 mm continuum emission reveals a dust condensation of ~23 Msun associated with the massive YSOs, deeply embedded in the strongest NH3 core (hereafter central core). The N2H+ emission reveals two main cores, the western and eastern core, located to the west and to the east of the mm condensation, and surrounded by a more extended and complex structure of ~0.5 pc. Toward the central core the N2H+ emission drops significantly, indicating a clear chemical differentiation in the region. We found low values of the NH3/N2H+ ratio ~50-100 toward the western/eastern cores, and high values up to 1000 in the central core. The chemical model indicates that density, and in particular temperature, are key parameters in determining the NH3/N2H+ ratio. The high density and temperature reached in the central core allow molecules like CO to evaporate from grain mantles. The CO desorption causes a significant destruction of N2H+, favoring the formation of HCO+. This result is supported by our observations, which show that N2H+ and HCO+ are anticorrelated in the central core. The observed values of the NH3/N2H+ ratio in the central core can be reproduced by our model for times t~4.5-5.3x10^5 yr (central) and t~10^4-3x10^6 yr (western/eastern). The NH3/N2H+ ratio in AFGL 5142 does not follow the same trend as in regions of low-mass star formation mainly due to the high temperature reached in hot cores.

preprint2010arXiv

The NH2D/NH3 ratio toward pre-protostellar cores around the UCHII region in IRAS 20293+3952

The deuterium fractionation, Dfrac, has been proposed as an evolutionary indicator in pre-protostellar and protostellar cores of low-mass star-forming regions. We investigate Dfrac, with high angular resolution, in the cluster environment surrounding the UCHII region IRAS 20293+3952. We performed high angular resolution observations with the IRAM Plateau de Bure Interferometer (PdBI) of the ortho-NH2D 1_{11}-1_{01} line at 85.926 GHz and compared them with previously reported VLA NH3 data. We detected strong NH2D emission toward the pre-protostellar cores identified in NH3 and dust emission, all located in the vicinity of the UCHII region IRAS 20293+3952. We found high values of Dfrac~0.1-0.8 in all the pre-protostellar cores and low values, Dfrac<0.1, associated with young stellar objects. The high values of Dfrac in pre-protostellar cores could be indicative of evolution, although outflow interactions and UV radiation could also play a role.

preprint2010arXiv

The reflection-symmetric wiggle of the young protostellar jet HH 211

HH 211 is a highly collimated jet originating from a nearby young Class 0 protostar. Here is a follow-up study of the jet with our previous observations at unprecedented resolution up to ~ 0.3" in SiO (J=8-7), CO (J=3-2), and SO (N_J=8_9-7_8). SiO, CO, and SO can all be a good tracer of the HH 211 jet, tracing the internal shocks in the jet. Although the emissions of these molecules show roughly the same morphology of the jet, there are detailed differences. In particular, the CO emission traces the jet closer to the source than the SiO and SO emissions. In addition, in the better resolved internal shocks, both the CO and SO emission are seen slightly ahead of the SiO emission. The jet is clearly seen on both sides of the source with more than one cycle of wiggle. The wiggle is reflection-symmetric about the source and can be reasonably fitted by an orbiting source jet model. The best-fit parameters suggest that the source itself could be a very low-mass protobinary with a total mass of ~ 60 M_Jup and a binary separation of ~ 4.6 AU. The abundances of SiO and SO in the gas phase are found to be highly enhanced in the jet as compared to the quiescent molecular clouds, even close to within 300 AU from the source where the dynamical time scale is <10 yrs. The abundance enhancements of these molecules are closely related to the internal shocks. The detected SiO is either the consequence of the release of Si-bearing material from dust grains or of its formation via gas chemistry in the shocks. The SO, on the other hand, seems to form via gas chemistry in the shocks.

preprint2009arXiv

A ring/disk/outflow system associated with W51 North: a very massive star in the making

Sensitive and high angular resolution ($\sim$ 0.4\arcsec) SO$_2$[22$_{2,20}$ $\to$ 22$_{1,21}$] and SiO[5$\to$4] line and 1.3 and 7 mm continuum observations made with the Submillimeter Array (SMA) and the Very Large Array (VLA) towards the young massive cluster W51 IRS2 are presented. We report the presence of a large (of about 3000 AU) and massive (40 M$_\odot$) dusty circumstellar disk and a hot gas molecular ring around a high-mass protostar or a compact small stellar system associated with W51 North. The simultaneous observations of the silicon monoxide molecule, an outflow gas tracer, further revealed a massive (200 M$_\odot$) and collimated ($\sim14^\circ$) outflow nearly perpendicular to the dusty and molecular structures suggesting thus the presence of a single very massive protostar with a bolometric luminosity of more than 10$^5$ L$_\odot$. A molecular hybrid LTE model of a Keplerian and infalling ring with an inner cavity and a central stellar mass of more than 60 M$_\odot$ agrees well with the SO$_2$[22$_{2,20}$ $\to$ 22$_{1,21}$] line observations. Finally, these results suggest that mechanisms, such as mergers of low- and intermediate- mass stars, might be not necessary for forming very massive stars.

preprint2009arXiv

Low-mass protostars and dense cores in different evolutionary stages in IRAS 00213+6530

We aim at studying with high angular resolution a dense core associated with a low-luminosity IRAS source, IRAS 00213+6530, in order to investigate whether low mass star formation is really taking place in isolation. We performed observations at 1.2mm with the IRAM 30m telescope, VLA observations at 6cm, 3.6cm, 1.3cm, 7mm, and H2O maser and NH3 lines, and observations with the NASA 70m antenna in CCS and H2O maser. The cm and mm continuum emission, together with the near infrared data from the 2MASS allowed us to identify 3 YSOs, IRS1, VLA8A, and VLA8B, with different radio and infrared properties, and which seem to be in different evolutionary stages. The NH3 emission consists of three clouds. Two of these, MM1 and MM2, are associated with dust emission, while the southern cloud is only detected in NH3. The YSOs are embedded in MM1, where we found evidence of line broadening and temperature enhancements. On the other hand, the southern cloud and MM2 appear to be quiescent and starless. We modeled the radial intensity profile at 1.2mm of MM1. The model fits reasonably well the data, but it underestimates the intensity at small projected distances from the 1.2mm peak, probably due to the presence of multiple YSOs embedded in the envelope. There is a differentiation in the relative NH3 abundance with low values, ~2x10^-8, toward MM1, and high values, up to 10^-6, toward the southern cloud and MM2, suggesting that these clouds could be in a young evolutionary stage. IRAS 00213+6530 is harboring a multiple system of low-mass protostars, indicating that star formation in this cloud is taking place in groups, rather than in isolation. The low-mass YSOs found in IRAS 00213+6530 are in different evolutionary stages suggesting that star formation is taking place in different episodes.

preprint2009arXiv

Rotation and Outflow motions in the very low-mass Class 0 protostellar system HH 211 at subarcsecond resolution

HH 211 is a nearby young protostellar system with a highly collimated jet. We have mapped it in 352 GHz continuum, SiO (J=8-7), and HCO+ (J=4-3) emission at up to ~ 0.2" resolution with the Submillimeter Array (SMA). The continuum source is now resolved into two sources, SMM1 and SMM2, with a separation of ~ 84 AU. SMM1 is seen at the center of the jet, probably tracing a (inner) dusty disk around the protostar driving the jet. SMM2 is seen to the southwest of SMM1 and may trace an envelope-disk around a small binary companion. A flattened envelope-disk is seen in HCO+ around SMM1 with a radius of ~ 80 AU perpendicular to the jet axis. Its velocity structure is consistent with a rotation motion and can be fitted with a Keplerian law that yields a mass of ~ 50+-15 Jupiter mass (a mass of a brown dwarf) for the protostar. Thus, the protostar could be the lowest mass source known to have a collimated jet and a rotating flattened envelope-disk. A small-scale (~ 200 AU) low-speed (~ 2 km/s) outflow is seen in HCOP+ around the jet axis extending from the envelope-disk. It seems to rotate in the same direction as the envelope-disk and may carry away part of the angular momentum from the envelope-disk. The jet is seen in SiO close to ~ 100 AU from SMM1. It is seen with a "C-shaped" bending. It has a transverse width of <~ 40 AU and a velocity of ~ 170+-60 km/s. A possible velocity gradient is seen consistently across its innermost pair of knots, with ~ 0.5 km/s at ~ 10 AU, consistent with the sense of rotation of the envelope-disk. If this gradient is an upper limit of the true rotational gradient of the jet, then the jet carries away a very small amount of angular momentum of ~ 5 AU km/s and thus must be launched from the very inner edge of the disk near the corotation radius.

preprint2009arXiv

Three intermediate-mass YSOs with different properties emerging from the same natal cloud in IRAS 00117+6412

We observed with the VLA, PdBI, and SMA the centimeter and millimeter continuum, N2H+(1-0), and CO(2-1) emission associated with a dusty cloud harboring a nascent cluster with intermediate-mass protostars. At centimeter wavelengths we found a strong source, tracing a UCHII region, at the eastern edge of the dusty cloud, with a shell-like structure, and with the near-infrared counterpart falling in the center of the shell. This is presumably the most massive source of the forming cluster. About 15'' to the west of the UCHII region and well embedded in the dusty cloud, we detected a strong millimeter source, MM1, associated with centimeter and near-infrared emission. MM1 seems to be driving a prominent high-velocity CO bipolar outflow, and is embedded in a ridge of dense gas traced by N2H+. We estimated that MM1 is an intermediate-mass source in the Class 0/I phase. About 15'' to the south of MM1, and still more deeply embedded in the dusty cloud, we detected a compact millimeter source, MM2, with neither centimeter nor near-infrared emission, but with water maser emission. MM2 is associated with a clump of N2H+, whose kinematics reveal a clear velocity gradient and additionally we found signposts of infall motions. MM2, being deeply embedded within the dusty cloud, with an associated water maser but no hints of CO outflow emission, is an intriguing object, presumably of intermediate mass. In conclusion, the UCHII region is found at the border of a dusty cloud which is currently undergoing active star formation. Two intermediate-mass protostars in the dusty cloud seem to have formed after the UCHII region and have different properties related to the outflow phenomenon.