Source author record

Leonardo Testi

Leonardo Testi 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

72works
6topics
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

72 published item(s)

preprint2026arXiv

Circumplanetary Disk Candidate in the Disk of HD 163296 Traced by Localized Emission from Simple Organics

Atacama Large Millimeter/submillimeter Array observations suggest that the disc of HD 163296 is being actively shaped by embedded, yet unseen protoplanets, as indicated by numerous gas and dust substructures consistent with planet-disc interaction models. We report the first detection of simple organic molecules, HCN and C2H, tracing a candidate circumplanetary disc (CPD) in the HD 163296 system, located at an orbital radius of $R=88\pm7$ au and azimuth $ϕ=46\pm3^\circ$ (or $R=0.75''$, $\rm{PA}=350^\circ$ in projected sky coordinates), and originating near the midplane of the circumstellar disc. The signature is localised but spectrally resolved, and it overlaps with a previously reported planet candidate, P94, identified through kinematic perturbations traced by CO lines. We propose a scenario in which the observed chemical anomalies arise from increased heating driven by the forming planet and ongoing accretion through its CPD, facilitating the thermal desorption of species that would otherwise remain frozen out in the disc midplane, and potentially triggering the activation barriers of chemical reactions that lead to enhanced molecular production. Based on a first-order dynamical analysis of the HCN spectrum from the CPD--isolated with a 7$σ$ significance--we infer an upper limit on the planet mass of 1.8 $M_{\rm Jup}$, consistent with predictions from CO kinematics and constraints from direct imaging studies. By comparing the CPD sizes derived from our models with theoretical expectations where the CPD radius corresponds to roughly one-third of the planet's Hill radius, we favor CPD gas temperatures $T > 150$ K, planet masses $M_{\rm p} < 1.0$ $M_{\rm Jup}$, and CPD radii $R_{\rm CPD} < 2$ au.

preprint2026arXiv

Water vapor emission at the warm cavity wall of the HD 100546 disk as revealed by ALMA

We present spatially resolved ALMA observations of the water line at 183 GHz in the disk around the Herbig star HD 100546. The water vapor emission peaks at the inner edge of the warm dust cavity, located ~15 au from the central star. We attribute this to thermal desorption at the water snowline, shifted outward at the dust cavity wall directly heated by the intense radiation. This represents the first spatially resolved image of the water snowline using ALMA observations of the main water isotopologue in a protoplanetary disk. The water emission morphology peaking inside the first dust ring is consistent with previous ALMA detections of oxygen-bearing complex organic molecules in the disk, including thermally desorbed methanol. These findings signal that warm cavities of transition disks provide ideal targets to directly reconstruct the spatial distribution of water vapor and the snowline location with ALMA, and directly connect water vapor emission to ice desorption of complex organic species.

preprint2022arXiv

A new planet candidate detected in a dust gap of the disc around HD 163296 through localised kinematic signatures. An observational validation of the Discminer

We report the robust detection of coherent, localised deviations from Keplerian rotation possibly associated with the presence of two giant planets embedded in the disc around HD 163296. The analysis is performed using the DISCMINER channel-map modelling framework on $^{12}$CO $J=2-1$ DSHARP data. Not only orbital radius, but also azimuth of the planets are retrieved by our technique. One of the candidate planets, detected at $R=94\pm6$ au, $ϕ=50\pm3^\circ$ (P94), is near the centre of one of the gaps in dust continuum emission, and is consistent with a planet mass of 1 $\rm M_{\rm Jup}$. The other planet, located at $R=261\pm4$ au, $ϕ=57\pm1^\circ$ (P261), is in the region where a velocity kink was previously observed in $^{12}$CO channel maps. Also, we provide a simultaneous description of the height and temperature of the upper and lower emitting surfaces of the disc, and propose the line width as a solid observable to track gas substructure. Using azimuthally averaged line width profiles we detect gas gaps at $R=38$ au, $R=88$ au, and $R=136$ au, closely matching the location of their dust and kinematical counterparts. Furthermore, we observe strong azimuthal asymmetries in line widths around the gas gap at $R=88$ au, possibly linked to turbulent motions driven by the P94 planet. Our results confirm that the DISCMINER is capable of finding localised, otherwise unseen velocity perturbations thanks to its robust statistical framework, but also that it is well suited for studies of the gas properties and vertical structure of protoplanetary discs.

preprint2022arXiv

Dust resurgence in protoplanetary disks due to planetesimal-planet interactions

Observational data on the dust content of circumstellar disks show that the median dust content in disks around pre-main sequence stars in nearby star forming regions seem to increase from about 1 Myr to about 2 Myr, and then decline with time. This behaviour challenges the models where the small dust grains steadily decline by accumulating into larger bodies and drifting inwards on a short timescale (less than about 1 Myr). In this Letter we explore the possibility to reconcile this discrepancy in the framework of a model where the early formation of planets dynamically stirs the nearby planetesimals and causes high energy impacts between them, resulting in the production of second-generation dust. We show that the observed dust evolution can be naturally explained by this process within a suite of representative disk-planet architectures.

preprint2022arXiv

External or internal companion exciting the spiral arms in CQ Tau?

We present new high-contrast images in near-infrared wavelengths ($λ$ = 1.04, 1.24, 1.62, 2.18 and 3.78$μ$m) of the young variable star CQ Tau, aiming to constrain the presence of companions in the protoplanetary disc. We reached a Ks-band contrast of 14 magnitudes with SPHERE/IRDIS at separations greater than 0."4 from the star. Our mass sensitivity curve rules out giant planets above 4 M$_{\rm Jup}$ immediately outside the spiral arms at $\sim$60 au and above 2-3 M$_{\rm Jup}$ beyond 100 au to 5$σ$ confidence assuming hot-start models. We do, however, detect four spiral arms, a double-arc and evidence for shadows in scattered light cast by a misaligned inner disc. Our observations may be explained by an unseen close-in companion on an inclined and eccentric orbit. Such a hypothesis would also account for the disc CO cavity and disturbed kinematics.

preprint2022arXiv

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

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

preprint2022arXiv

Observations of PAHs in the atmospheres of discs and exoplanets

Polycyclic aromatic hydrocarbons (PAHs) play a key role in the chemical and hydrodynamical evolution of the atmospheres of exoplanets and planet-forming discs. If they can survive the planet formation process, PAHs are likely to be involved in pre-biotic chemical reactions eventually leading to more complex molecules such as amino acids and nucleotides, which form the basis for life as we know it. However, the abundance and specific role of PAHs in these environments is largely unknown due to limitations in sensitivity and range of wavelength of current and previous space-borne facilities. Upcoming infrared space spectroscopy missions, such as Twinkle and Ariel, present a unique opportunity to detect PAHs in the atmospheres of exoplanets and planet-forming discs. In this work we present synthetic observations based on conservative numerical modeling of typical planet-forming discs and a transiting hot Saturnian planet around solar type star. Our models show that Twinkle and Ariel might both be able to detect the 3.3 micron PAH feature within reasonable observing time in discs and transiting planets, assuming that PAHs are present with an abundance of at least one tenth of the interstellar medium value.

preprint2022arXiv

On the time evolution of the $M_{\rm d} - M_\star$ and $\dot M - M_\star$ correlations for protoplanetary discs: the viscous timescale increases with stellar mass

Large surveys of star-forming regions have unveiled power-law correlations between the stellar mass and the disc parameters, such as the disc mass $M_{\mathrm{d}} \propto {M_{\star}}^{λ_{\mathrm{m}}}$ and the accretion rate $\dot M \propto {M_{\star}}^{λ_{\mathrm{acc}}}$. The observed slopes appear to be increasing with time, but the reason behind the establishment of these correlations and their subsequent evolution is still uncertain. We conduct a theoretical analysis of the impact of viscous evolution on power-law initial conditions for a population of protoplanetary discs. We find that, for evolved populations, viscous evolution enforces the two correlations to have the same slope, $λ_{\mathrm{m}}$ = $λ_{\mathrm{acc}}$, and that this limit is uniquely determined by the initial slopes $λ_{\mathrm{m}, 0}$ and $λ_{\mathrm{acc}, 0}$. We recover the increasing trend claimed from the observations when the difference in the initial values, $δ_0 = λ_{\mathrm{m}, 0} - λ_{\mathrm{acc}, 0}$, is larger than $1/2$; moreover, we find that this increasing trend is a consequence of a positive correlation between the viscous timescale and the stellar mass. We also present the results of disc population synthesis numerical simulations, that allow us to introduce a spread and analyse the effect of sampling, which show a good agreement with our analytical predictions. Finally, we perform a preliminary comparison of our numerical results with observational data, which allows us to constrain the parameter space of the initial conditions to $λ_{\mathrm{m}, 0} \in [1.2, 2.1]$, $λ_{\mathrm{acc}, 0} \in [0.7, 1.5]$.

preprint2022arXiv

The Galactic dynamics revealed by the filamentary structure in atomic hydrogen emission

We present a study of the filamentary structure in the atomic hydrogen (HI) emission at the 21 cm wavelength toward the Galactic plane using the observations in the HI4PI survey. Using the Hessian matrix method across radial velocity channels, we identified the filamentary structures and quantified their orientations using circular statistics. We found that the regions of the Milky Way's disk beyond 10 kpc and up to roughly 18 kpc from the Galactic center display HI filamentary structures predominantly parallel to the Galactic plane. For regions at lower Galactocentric radii, we found that the HI filaments are mostly perpendicular or do not have a preferred orientation with respect to the Galactic plane. We interpret these results as the imprint of supernova feedback in the inner Galaxy and Galactic rotation in the outer Milky Way. We found that the HI filamentary structures follow the Galactic warp and that they highlight some of the variations interpreted as the effect of the gravitational interaction with satellite galaxies. In addition, the mean scale height of the filamentary structures is lower than that sampled by the bulk of the HI emission, thus indicating that the cold and warm atomic hydrogen phases have different scale heights in the outer galaxy. Finally, we found that the fraction of the column density in HI filaments is almost constant up to approximately 18 kpc from the Galactic center. This is possibly a result of the roughly constant ratio between the cold and warm atomic hydrogen phases inferred from the HI absorption studies. Our results indicate that the HI filamentary structures provide insight into the dynamical processes shaping the Galactic disk. Their orientations record how and where the stellar energy input, the Galactic fountain process, the cosmic ray diffusion, and the gas accretion have molded the diffuse interstellar medium in the Galactic plane.

preprint2022arXiv

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

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

preprint2020arXiv

ALMA chemical survey of disk-outflow sources in Taurus (ALMA-DOT) V: Sample, overview, and demography of disk molecular emission

We present an overview of the ALMA chemical survey of disk-outflow sources in Taurus (ALMA-DOT), a campaign devoted to the characterization of the molecular emission from partly embedded, young stars. The project aims at better understanding the gaseous products that are delivered to planets by means of high-resolution maps of assorted lines probing disks at the time of the planet formation (less than 1 Myr). Nine different molecules are surveyed by our observations of six Class I/flat-spectrum sources. A series of accompanying articles analyze specific targets and molecules. Here we describe the sample and provide a general overview of the results, focusing on the spatial distribution, column densities, and abundance ratios of H$_2$CO, CS, and CN. The results of this work are a first step toward the characterization of the disk chemical evolution that need to be complemented by further observations of less exceptional disks and customized thermo-chemical modeling.

preprint2020arXiv

ALMA chemical survey of disk-outflow sources in Taurus (ALMA-DOT): I. CO, CS, CN, and H2CO around DG Tau B

The chemical composition of planets is inherited by the distribution of the various molecular species in the protoplanetary disk at thetime of their formation. As of today, only a handful of disks has been imaged in multiple spectral lines with high spatial resolution. As part of a small campaign devoted to the chemical characterization of disk-outflow sources in Taurus, we report on new ALMA Band 6 observations with 20 au resolution toward the embedded young star DG Tau B. Images of the continuum emission reveals a dust disk with rings and, putatively, a leading spiral arm. The disk, as well as the prominent outflow cavities, are detected in CO, H2CO, CS, and CN while they remain undetected in SO2, HDO, and CH3OH. From the absorption of the back-side outflow, we inferred that the disk emission is optically thick in the inner 50 au. This morphology explains why no line emission is detected from this inner region and poses some limitations toward the calculation of the dust mass and the characterization of the inner gaseous disk. The H2CO and CS emission from the inner 200 au is mostly from the disk and their morphology is very similar. The CN emission significantly differs from the other two molecules as it is observed only beyond 150 au. This ring-like morphology is consistent with previous observations and the predictions of thermochemical disk models. Finally, we constrained the disk-integrated column density of all molecules. In particular, we found that the CH3OH/H2CO ratio must be smaller than 2, making the methanol non-detection still consistent with the only such a ratio available from the literature (1.27 in TW Hya).

preprint2020arXiv

CAB: Towards the RNA-world in the interstellar medium -- detection of urea, and search of 2-amino-oxazole and simple sugars

In the past decade, Astrochemistry has witnessed an impressive increase in the number of detections of complex organic molecules. Some of these species are of prebiotic interest such as glycolaldehyde, the simplest sugar, or amino acetonitrile, a possible precursor of glycine. Recently, we have reported the detection of two new nitrogen-bearing complex organics, glycolonitrile and Z-cyanomethanimine, known to be intermediate species in the formation process of ribonucleotides within theories of a primordial ribonucleic acid (RNA)-world for the origin of life. In this paper, we present deep and high-sensitivity observations toward two of the most chemically rich sources in the Galaxy: a Giant Molecular Cloud in the center of the Milky Way (G+0.693-0.027) and a proto-Sun (IRAS16293-2422 B). Our aim is to explore whether the key precursors considered to drive the primordial RNA-world chemistry, are also found in space. Our high-sensitivity observations reveal that urea is present in G+0.693-0.027 with an abundance of about 5x10-11. This is the first detection of this prebiotic species outside a star-forming region. Urea remains undetected toward the proto-Sun IRAS16293-2422 B (upper limit to its abundance of less than 2x10-11). Other precursors of the RNA-world chemical scheme such as glycolaldehyde or cyanamide are abundant in space, but key prebiotic species such as 2- amino-oxazole, glyceraldehyde or dihydroxyacetone are not detected in either source. Future more sensitive observations targeting the brightest transitions of these species will be needed to disentangle whether these large prebiotic organics are certainly present in space.

preprint2020arXiv

Dynamical stellar masses of pre-main sequence stars in Lupus and Taurus obtained with ALMA surveys in comparison with stellar evolutionary models

We analysed archival molecular line data of pre-main sequence (PMS) stars in the Lupus and Taurus star-forming regions obtained with ALMA surveys with an integration time of a few minutes per source. We stacked the data of $^{13}$CO and C$^{18}$O (J = 2-1 & 3-2) and CN (N = 3-2, J = 7/2-5/2) lines to enhance the signal-to-noise ratios, and measured the stellar masses of 45 out of 67 PMS stars from the Keplerian rotation in their circumstellar disks. The measured dynamical stellar masses were compared to the stellar masses estimated from the spectroscopic measurements with seven different stellar evolutionary models. We found that the magnetic model of Feiden (2016) provides the best estimate of the stellar masses in the mass range of $0.6~M_{\odot}\leq M_{\star} \leq 1.3~M_{\odot}$ with a deviation of $<$0.7$σ$ from the dynamical masses, while all the other models underestimate the stellar masses in this mass range by 20% to 40%. In the mass range of $<0.6~M_{\odot}$, the stellar masses estimated with the magnetic model of Feiden (2016) have a larger deviation ($>2σ$) from the dynamical masses, and other, non-magnetic stellar evolutionary models of Siess et al. (2000), Baraffe et al. (2015) and Feiden (2016) show better agreements with the dynamical masses with the deviations of 1.4$σ$ to 1.6$σ$. Our results show the mass dependence of the accuracy of these stellar evolutionary models.

preprint2020arXiv

Hot Corinos Chemical Diversity: Myth or Reality?

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

preprint2020arXiv

HST survey of the Orion Nebula Cluster in the H$_2$O 1.4 $μ$m absorption band: I. A census of substellar and planetary mass objects

In order to obtain a complete census of the stellar and sub-stellar population, down to a few M$_{Jup}$ in the $\sim1$ Myr old Orion Nebula Cluster, we used the infrared channel of the Wide Field Camera 3 of the Hubble Space Telescope with the F139M and F130N filters. These bandpasses correspond to the $1.4 μ$m H$_2$O absorption feature and an adjacent line-free continuum region. Out of $4,504$ detected sources, $3,352$ (about $75\%$) appear fainter than m$_{130}=14$ (Vega mag) in the F130N filter, a brightness corresponding to the hydrogen-burning limit mass (M$\simeq 0.072 M_\odot$) at $\sim 1$ Myr. Of these, however, only $742$ sources have a negative F130M-139N color index, indicative of the presence of H$_2$O vapor in absorption, and can therefore be classified as bona-fide M and L dwarfs, with effective temperatures T$\lesssim 2850$ K at an assumed $1$ Myr cluster age. On our color-magnitude diagram, this population of sources with H$_2$O absorption appears clearly distinct from the larger background population of highly reddened stars and galaxies with positive F130M-F139N color index, and can be traced down to the sensitivity limit of our survey, m$_{130}\simeq 21.5$, corresponding to a $1$ Myr old $\simeq 3 $M$_{Jup}$, planetary mass object under about 2 magnitudes of visual extinction. Theoretical models of the BT-Settl family predicting substellar isochrones of $1, 2$ and $3$ Myr (down to $\sim 1 $M$_{Jup}$) fail to reproduce the observed H$_2$O color index at M$\lesssim 20 $M$_{Jup}$. We perform a Bayesian analysis to determine extinction, mass and effective temperature of each sub-stellar member of our sample, together with its membership probability.

preprint2020arXiv

Near-Infrared Imaging of a Spiral in the CQ Tau Disk

We present $L^\prime$-band Keck/NIRC2 imaging and $H$-band Subaru/AO188+HiCIAO polarimetric observations of CQ Tau disk with a new spiral arm. Apart from the spiral feature our observations could not detect any companion candidates. We traced the spiral feature from the $r^2$-scaled HiCIAO polarimetric intensity image and the fitted result is used for forward modeling to reproduce the ADI-reduced NIRC2 image. We estimated the original surface brightness after throughput correction in $L^\prime$-band to be $\sim126$ mJy/arcsec$^2$ at most. We suggest that the grain temperature of the spiral may be heated up to $\sim$200 K in order to explain both of the $H$- and $L^{\prime}$-bands results. The $H$-band emission at the location of the spiral originates from the scattering from the disk surface while both scattering and thermal emission may contribute to the $L^{\prime}$-band emission. If the central star is only the light source of scattered light, the spiral emission at $L^\prime$-band should be thermal emission. If an inner disk also acts as the light source, the scattered light and the thermal emission may equally contribute to the $L^\prime$-band spiral structure.

preprint2020arXiv

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

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

preprint2020arXiv

Size and structures of disks around very low mass stars in the Taurus star-forming region

We aim to estimate if structures, such as cavities, rings, and gaps, are common in disks around VLMS and to test models of structure formation in these disks. We also aim to compare the radial extent of the gas and dust emission in disks around VLMS, which can give us insight about radial drift. We studied six disks around VLMS in the Taurus star-forming region using ALMA Band 7 ($\sim 340\,$GHz) at a resolution of $\sim0.1''$. The targets were selected because of their high disk dust content in their stellar mass regime. Our observations resolve the disk dust continuum in all disks. In addition, we detect the $^{12}$CO ($J=3-2$) emission line in all targets and $^{13}$CO ($J=3-2$) in five of the six sources. The angular resolution allows the detection of dust substructures in three out of the six disks, which we studied by using UV-modeling. Central cavities are observed in the disks around stars MHO\,6 (M5.0) and CIDA\,1 (M4.5), while we have a tentative detection of a multi-ringed disk around J0433. Single planets of masses $0.1\sim0.4\,M_{\rm{Jup}}$ would be required. The other three disks with no observed structures are the most compact and faintest in our sample. The emission of $^{12}$CO and $^{13}$CO is more extended than the dust continuum emission in all disks of our sample. When using the $^{12}$CO emission to determine the gas disk extension $R_{\rm{gas}}$, the ratio of $R_{\rm{gas}}/R_{\rm{dust}}$ in our sample varies from 2.3 to 6.0, which is consistent with models of radial drift being very efficient around VLMS in the absence of substructures. Our observations do not exclude giant planet formation on the substructures observed. A comparison of the size and luminosity of VLMS disks with their counterparts around higher mass stars shows that they follow a similar relation.

preprint2020arXiv

The evolution of dust-disk sizes from a homogeneous analysis of 1-10 Myr-old stars

We utilize ALMA archival data to estimate the dust disk size of 152 protoplanetary disks in Lupus (1-3 Myr), Chamaeleon I (2-3 Myr), and Upper-Sco (5-11 Myr). We combine our sample with 47 disks from Tau/Aur and Oph whose dust disk radii were estimated, as here, through fitting radial profile models to visibility data. We use these 199 homogeneously derived disk sizes to identify empirical disk-disk and disk-host property relations as well as to search for evolutionary trends. In agreement with previous studies, we find that dust disk sizes and millimeter luminosities are correlated, but show for the first time that the relationship is not universal between regions. We find that disks in the 2-3 Myr-old Cha I are not smaller than disks in other regions of similar age, and confirm the Barenfeld et al. (2017) finding that the 5-10 Myr USco disks are smaller than disks belonging to younger regions. Finally, we find that the outer edge of the Solar System, as defined by the Kuiper Belt, is consistent with a population of dust disk sizes which have not experienced significant truncation.

preprint2019arXiv

Long lived dust rings around HD169142

Recent ALMA observations of the protoplanetary disc around HD~169142 reveal a peculiar structure made of concentric dusty rings: a main ring at $\sim$20 au, a triple system of rings at $\sim 55-75$ au in millimetric continuum emission and a perturbed gas surface density from the $^{12}$CO,$^{13}$CO and C$^{18}$O$(J=2-1)$ surface brightness profile. In this Letter, we perform three-dimensional numerical simulations and radiative transfer modeling exploring the possibility that two giant planets interacting with the disc and orbiting in resonant locking can be responsible for the origin of the observed dust inner rings structure. We find that in this configuration the dust structure is actually long lived while the gas mass of the disc is accreted onto the star and the giant planets, emptying the inner region. In addition, we also find that the innermost planet is located at the inner edge of the dust ring, and can accrete mass from the disc, generating a signature in the dust ring shape that can be observed in mm ALMA observations.

preprint2018arXiv

Gas density perturbations induced by forming planet(s) in the AS 209 protoplanetary disk as seen with ALMA

The formation of planets occurs within protoplanetary disks surrounding young stars, resulting in perturbation of the gas and dust surface densities. Here, we report the first evidence of spatially resolved gas surface density ($Σ_{g}$) perturbation towards the AS~209 protoplanetary disk from the optically thin C$^{18}$O ($J=2-1$) emission. The observations were carried out at 1.3~mm with ALMA at a spatial resolution of about 0.3$\arcsec$ $\times$ 0.2$\arcsec$ (corresponding to $\sim$ 38 $\times$ 25 au). The C$^{18}$O emission shows a compact ($\le$60~au), centrally peaked emission and an outer ring peaking at 140~au, consistent with that observed in the continuum emission and, its azimuthally averaged radial intensity profile presents a deficit that is spatially coincident with the previously reported dust map. This deficit can only be reproduced with our physico-thermochemical disk model by lowering $Σ_{gas}$ by nearly an order of magnitude in the dust gaps. Another salient result is that contrary to C$^{18}$O, the DCO$^{+}$ ($J=3-2$) emission peaks between the two dust gaps. We infer that the best scenario to explain our observations (C$^{18}$O deficit and DCO$^{+}$ enhancement) is a gas perturbation due to forming-planet(s), that is commensurate with previous continuum observations of the source along with hydrodynamical simulations. Our findings confirm that the previously observed dust gaps are very likely due to perturbation of the gas surface density that is induced by a planet of at least 0.2~M$\rm_{Jupiter}$ in formation. Finally, our observations also show the potential of using CO isotopologues to probe the presence of saturn mass planet(s).

preprint2018arXiv

Investigating fragmentation of gas structures in OB cluster-forming molecular clump G33.92+0.11 with 1000 AU resolution observations of ALMA

We report new, $\sim$1000 AU spatial resolution observations of 225 GHz dust continuum emission towards the OB cluster-forming molecular clump G33.92+0.11. On parsec scales, this molecular clump presents a morphology with several arm-like dense gas structures surrounding the two central massive ($\gtrsim$100 $M_{\odot}$) cores. From the new, higher resolution observations, we identified 28 localized, spatially compact dust continuum emission sources, which may be candidates of young stellar objects. Only one of them is not embedded within known arm-like (or elongated) dense gas structures. The spatial separations of these compact sources can be very well explained by Jeans lengths. We found that G33.92+0.11 may be consistently described by a marginally centrifugally supported, Toomre unstable accretion flow which is approximately in a face-on projection. The arm-like overdensities are natural consequence of the Toomre instability, which can fragment to form young stellar objects in shorter time scales than the timescale of the global clump contraction. On our resolved spatial scales, there is not yet evidence that the fragmentation is halted by turbulence, magnetic field, or stellar feedback.

preprint2016arXiv

A Census of Large-Scale ($\ge$ 10 pc), Velocity-Coherent, Dense Filaments in the Northern Galactic Plane: Automated Identification Using Minimum Spanning Tree

Large-scale gaseous filaments with length up to the order of 100 pc are on the upper end of the filamentary hierarchy of the Galactic interstellar medium. Their association with respect to the Galactic structure and their role in Galactic star formation are of great interest from both observational and theoretical point of view. Previous "by-eye" searches, combined together, have started to uncover the Galactic distribution of large filaments, yet inherent bias and small sample size limit conclusive statistical results to be drawn. Here, we present (1) a new, automated method to identify large-scale velocity-coherent dense filaments, and (2) the first statistics and the Galactic distribution of these filaments. We use a customized minimum spanning tree algorithm to identify filaments by connecting voxels in the position-position-velocity space, using the Bolocam Galactic Plane Survey spectroscopic catalog. In the range of $7.^{\circ}5 \le l \le 194^{\circ}$, we have identified 54 large-scale filaments and derived mass ($\sim 10^3 - 10^5 \, M_\odot$), length (10-276 pc), linear mass density (54-8625 $M_\odot \, \rm{pc}^{-1}$), aspect ratio, linearity, velocity gradient, temperature, fragmentation, Galactic location and orientation angle. The filaments concentrate along major spiral arms. They are widely distributed across the Galactic disk, with 50% located within $\pm$20 pc from the Galactic mid-plane and 27% run in the center of spiral arms (aka "bones"). An order of 1% of the molecular ISM is confined in large filaments. Massive star formation is more favorable in large filaments compared to elsewhere. This is the first comprehensive catalog of large filaments useful for a quantitative comparison with spiral structures and numerical simulations.

preprint2016arXiv

ALMA Cycle 0 Publication Statistics

The scientific impact of a facility is the most important measure of its success. Monitoring and analysing the scientific return can help to modify and optimise operations and adapt to the changing needs of scientific research. The methodology that we have developed to monitor the scientific productivity of the ALMA Observatory, as well as the first results, are described. We focus on the outcome of the first cycle (Cycle 0) of ALMA Early Science operations. Despite the fact that only two years have passed since the completion of Cycle 0 and operations have already changed substantially, this analysisconfirms the effectiveness of the underlying concepts. We find that ALMA is fulfilling its promise as a transformational facility for the observation of the Universe in the submillimetre.

preprint2016arXiv

Brown dwarf disks with Herschel: Linking far-infrared and (sub)-mm fluxes

Brown dwarf disks are excellent laboratories to test our understanding of disk physics in an extreme parameter regime. In this paper we investigate a sample of 29 well-characterized brown dwarfs and very low mass stars, for which Herschel far-infrared fluxes as well as (sub)-mm fluxes are available. We have measured new Herschel PACS fluxes for 11 objects and complement these with (sub)-mm data and Herschel fluxes from the literature. We analyze their spectral energy distributions in comparison with results from radiative transfer modeling. Fluxes in the far-infrared are strongly affected by the shape and temperature of the disk (and hence stellar luminosity), whereas the (sub)-mm fluxes mostly depend on disk mass. Nevertheless, there is a clear correlation between far-infrared and (sub)-mm fluxes. We argue that the link results from the combination of the stellar mass-luminosity relation and a scaling between disk mass and stellar mass. We find strong evidence of dust settling to the disk midplane. The spectral slopes between near- and far-infrared are mostly between $-0.5$ and $-1.2$ in our sample, comparable to more massive T Tauri stars, which may imply that the disk shapes are similar as well, though highly-flared disks are rare among brown dwarfs. We find that dust temperatures in the range of 7-15 K, calculated with $T\approx25\,(L/L_\odot)^{0.25}$ K, are appropriate for deriving disk masses from (sub)-mm fluxes for these low luminosity objects. About half of our sample hosts disks with at least one Jupiter mass, confirming that many brown dwarfs harbour sufficient material for the formation of Earth-mass planets in their midst.

preprint2016arXiv

Calibration of evolutionary diagnostics in high-mass star formation

The evolutionary classification of massive clumps that are candidate progenitors of high-mass young stars and clusters relies on a variety of independent diagnostics based on observables from the near-infrared to the radio. A promising evolutionary indicator for massive and dense cluster-progenitor clumps is the L/M ratio between the bolometric luminosity and the mass of the clumps. With the aim of providing a quantitative calibration for this indicator we used SEPIA/APEX to obtain CH3C2H(12-11) observations, that is an excellent thermometer molecule probing densities > 10^5 cm^-3 , toward 51 dense clumps with M>1000 solar masses, and uniformly spanning -2 < Log(L/M) < 2.3. We identify three distinct ranges of L/M that can be associated to three distinct phases of star formation in massive clumps. For L/M <1 no clump is detected in CH3C2H , suggesting an inner envelope temperature below 30K. For 1< L/M < 10 we detect 58% of the clumps, with a temperature between 30 and 35 K independently from the exact value of L/M; such clumps are building up luminosity due to the formation of stars, but no star is yet able to significantly heat the inner clump regions. For L/M> 10 we detect all the clumps, with a gas temperature rising with Log(L/M), marking the appearance of a qualitatively different heating source within the clumps; such values are found towards clumps with UCHII counterparts, suggesting that the quantitative difference in T - L/M behaviour above L/M >10 is due to the first appearance of ZAMS stars in the clumps.

preprint2016arXiv

CARMA Large Area Star Formation Survey: Dense Gas in the Young L1451 Region of Perseus

We present a 3 mm spectral line and continuum survey of L1451 in the Perseus Molecular Cloud. These observations are from the CARMA Large Area Star Formation Survey (CLASSy), which also imaged Barnard 1, NGC 1333, Serpens Main and Serpens South. L1451 is the survey region with the lowest level of star formation activity---it contains no confirmed protostars. HCO+, HCN, and N2H+ (J=1-0) are all detected throughout the region, with HCO+ the most spatially widespread, and molecular emission seen toward 90% of the area above N(H_2) column densities of 1.9x10^21 cm^-2. HCO+ has the broadest velocity dispersion, near 0.3 km/s on average, compared to ~0.15 km/s for the other molecules, thus representing a range from supersonic to subsonic gas motions. Our non-binary dendrogram analysis reveals that the dense gas traced by each molecule has similar hierarchical structure, and that gas surrounding the candidate first hydrostatic core (FHSC), L1451-mm, and other previously detected single-dish continuum clumps have similar hierarchical structure; this suggests that different sub-regions of L1451 are fragmenting on the pathway to forming young stars. We determined the three-dimensional morphology of the largest detectable dense gas structures to be relatively ellipsoidal compared to other CLASSy regions, which appeared more flattened at largest scales. A virial analysis shows the most centrally condensed dust structures are likely unstable against collapse. Additionally, we identify a new spherical, centrally condensed N2H+ feature that could be a new FHSC candidate. The overall results suggest L1451 is a young region starting to form its generation of stars within turbulent, hierarchical structures.

preprint2016arXiv

Cloud Structure of Galactic OB Cluster Forming Regions from Combining Ground and Space Based Bolometric Observations

We have developed an iterative procedure to systematically combine the millimeter and submillimeter images of OB cluster-forming molecular clouds, which were taken by ground based (CSO, JCMT, APEX, IRAM-30m) and space telescopes (Herschel, Planck). For the seven luminous ($L$$>$10$^{6}$ $L_{\odot}$) Galactic OB cluster-forming molecular clouds selected for our analyses, namely W49A, W43-Main, W43-South, W33, G10.6-0.4, G10.2-0.3, G10.3-0.1, we have performed single-component, modified black-body fits to each pixel of the combined (sub)millimeter images, and the Herschel PACS and SPIRE images at shorter wavelengths. The $\sim$10$"$ resolution dust column density and temperature maps of these sources revealed dramatically different morphologies, indicating very different modes of OB cluster-formation, or parent molecular cloud structures in different evolutionary stages. The molecular clouds W49A, W33, and G10.6-0.4 show centrally concentrated massive molecular clumps that are connected with approximately radially orientated molecular gas filaments. The W43-Main and W43-South molecular cloud complexes, which are located at the intersection of the Galactic near 3-kpc (or Scutum) arm and the Galactic bar, show a widely scattered distribution of dense molecular clumps/cores over the observed $\sim$10 pc spatial scale. The relatively evolved sources G10.2-0.3 and G10.3-0.1 appear to be affected by stellar feedback, and show a complicated cloud morphology embedded with abundant dense molecular clumps/cores. We find that with the high angular resolution we achieved, our visual classification of cloud morphology can be linked to the systematically derived statistical quantities (i.e., the enclosed mass profile, the column density probability distribution function, the two-point correlation function of column density, and the probability distribution function of clump/core separations).

preprint2016arXiv

Probing the wind launching regions of the Herbig Be star HD 58647 with high spectral resolution interferometry

We present a study of the wind launching region of the Herbig Be star HD 58647 using high angular (lambda/2B=0.003") and high spectral (R=12000) resolution interferometric VLTI-AMBER observations of the near-infrared hydrogen emission line, Br-gamma. The star displays double peaks in both Br-gamma line profile and wavelength-dependent visibilities. The wavelength-dependent differential phases show S-shaped variations around the line centre. The visibility level increases in the line (by ~0.1) at the longest projected baseline (88 m), indicating that the size of the line emission region is smaller than the size of the K-band continuum-emitting region, which is expected to arise near the dust sublimation radius of the accretion disc. The data have been analysed using radiative transfer models to probe the geometry, size and physical properties of the wind that is emitting Br-gamma. We find that a model with a small magnetosphere and a disc wind with its inner radius located just outside of the magnetosphere can well reproduce the observed Br-gamma profile, wavelength-dependent visibilities, differential and closure phases, simultaneously. The mass-accretion and mass-loss rates adopted for the model are Mdot_a = 3.5 x 10^{-7} Msun/yr and Mdot_dw = 4.5 x 10^{-8} Msun/yr, respectively (Mdot_dw/Mdot_a =0.13). Consequently, about 60 per cent of the angular momentum loss rate required for a steady accretion with the measured accretion rate is provide by the disc wind. The small magnetosphere in HD 58647 does not contribute to the Br-gamma line emission significantly.

preprint2016arXiv

Spiral Density Waves in a Young Protoplanetary Disk

Gravitational forces are expected to excite spiral density waves in protoplanetary disks, disks of gas and dust orbiting young stars. However, previous observations that showed spiral structure were not able to probe disk midplanes, where most of the mass is concentrated and where planet formation takes place. Using the Atacama Large Millimeter/submillimeter Array we detected a pair of trailing symmetric spiral arms in the protoplanetary disk surrounding the young star Elias 2-27. The arms extend to the disk outer regions and can be traced down to the midplane. These millimeter-wave observations also reveal an emission gap closer to the star than the spiral arms. We argue that the observed spirals trace shocks of spiral density waves in the midplane of this young disk.

preprint2016arXiv

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

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

preprint2016arXiv

The VLA view of the HL Tau Disk - Disk Mass, Grain Evolution, and Early Planet Formation

The first long-baseline ALMA campaign resolved the disk around the young star HL Tau into a number of axisymmetric bright and dark rings. Despite the very young age of HL Tau these structures have been interpreted as signatures for the presence of (proto)planets. The ALMA images triggered numerous theoretical studies based on disk-planet interactions, magnetically driven disk structures, and grain evolution. Of special interest are the inner parts of disks, where terrestrial planets are expected to form. However, the emission from these regions in HL Tau turned out to be optically thick at all ALMA wavelengths, preventing the derivation of surface density profiles and grain size distributions. Here, we present the most sensitive images of HL Tau obtained to date with the Karl G. Jansky Very Large Array at 7.0 mm wavelength with a spatial resolution comparable to the ALMA images. At this long wavelength the dust emission from HL Tau is optically thin, allowing a comprehensive study of the inner disk. We obtain a total disk dust mass of 0.001 - 0.003 Msun, depending on the assumed opacity and disk temperature. Our optically thin data also indicate fast grain growth, fragmentation, and formation of dense clumps in the inner densest parts of the disk. Our results suggest that the HL Tau disk may be actually in a very early stage of planetary formation, with planets not already formed in the gaps but in the process of future formation in the bright rings.

preprint2016arXiv

Using young massive star clusters to understand star formation and feedback in high-redshift-like environments

The formation environment of stars in massive stellar clusters is similar to the environment of stars forming in galaxies at a redshift of 1 - 3, at the peak star formation rate density of the Universe. As massive clusters are still forming at the present day at a fraction of the distance to high-redshift galaxies they offer an opportunity to understand the processes controlling star formation and feedback in conditions similar to those in which most stars in the Universe formed. Here we describe a system of massive clusters and their progenitor gas clouds in the centre of the Milky Way, and outline how detailed observations of this system may be able to: (i) help answer some of the fundamental open questions in star formation and (ii) quantify how stellar feedback couples to the surrounding interstellar medium in this high-pressure, high-redshift analogue environment.

preprint2015arXiv

Dense gas in the Galactic central molecular zone is warm and heated by turbulence

The Galactic center is the closest region in which we can study star formation under extreme physical conditions like those in high-redshift galaxies. We measure the temperature of the dense gas in the central molecular zone (CMZ) and examine what drives it. We mapped the inner 300 pc of the CMZ in the temperature-sensitive J = 3-2 para-formaldehyde (p-H$_2$CO) transitions. We used the $3_{2,1} - 2_{2,0} / 3_{0,3} - 2_{0,2}$ line ratio to determine the gas temperature in $n \sim 10^4 - 10^5 $cm$^{-3}$ gas. We have produced temperature maps and cubes with 30" and 1 km/s resolution and published all data in FITS form. Dense gas temperatures in the Galactic center range from ~60 K to > 100 K in selected regions. The highest gas temperatures T_G > 100 K are observed around the Sgr B2 cores, in the extended Sgr B2 cloud, the 20 km/s and 50 km/s clouds, and in "The Brick" (G0.253+0.016). We infer an upper limit on the cosmic ray ionization rate $ζ_{CR} < 10^{-14}$ 1/s. The dense molecular gas temperature of the region around our Galactic center is similar to values found in the central regions of other galaxies, in particular starburst systems. The gas temperature is uniformly higher than the dust temperature, confirming that dust is a coolant in the dense gas. Turbulent heating can readily explain the observed temperatures given the observed line widths. Cosmic rays cannot explain the observed variation in gas temperatures, so CMZ dense gas temperatures are not dominated by cosmic ray heating. The gas temperatures previously observed to be high in the inner ~75 pc are confirmed to be high in the entire CMZ.

preprint2015arXiv

Dust trapping by spiral arms in gravitationally unstable protostellar discs

In this paper we discuss the influence of gravitational instabilities in massive protostellar discs on the dynamics of dust grains. Starting from a Smoothed Particle Hydrodynamics (SPH) simulation, we have computed the evolution of the dust in a quasi-static gas density structure typical of self-gravitating disc. For different grain size distributions we have investigated the capability of spiral arms to trap particles. We have run 3D radiative transfer simulations in order to construct maps of the expected emission at (sub-)millimetre and near-infrared wavelengths. Finally, we have simulated realistic observations of our disc models at (sub-)millimetre and near-infrared wavelengths as they may appear with the Atacama Large Millimetre/sub-millimetre Array (ALMA) and the High-Contrast Coronographic Imager for Adaptive Optics (HiCIAO) in order to investigate whether there are observational signatures of the spiral structure. We find that the pressure inhomogeites induced by gravitational instabilities produce a non-negligible dynamical effect on centimetre sized particles leading to significant overdensities in spiral arms. We also find that the spiral structure is readily detectable by ALMA over a wide range of (sub-)millimetre wavelengths and by HiCIAO in near-infrared scattered light for non-face-on discs located in the Ophiucus star-forming region. In addition, we find clear spatial spectral index variations across the disc, revealing that the dust trapping produces a migration of large grains that can be potentially investigated through multi-wavelenghts observations in the (sub-)millimetric. Therefore, the spiral arms observed to date in protoplanetary disc might be interpreted as density waves induced by the development of gravitational instabilities.

preprint2015arXiv

Grain Growth in the Circumstellar Disks of the Young Stars CY Tau and DoAr 25

We present new results from the Disks@EVLA program for two young stars: CY Tau and DoAr 25. We trace continuum emission arising from their circusmtellar disks from spatially resolved observations, down to tens of AU scales, at λ = 0.9, 2.8, 8.0, and 9.8 mm for DoAr25 and at λ = 1.3, 2.8, and 7.1 mm for CY Tau. Additionally, we constrain the amount of emission whose origin is different from thermal dust emission from 5 cm observations. Directly from interferometric data, we find that observations at 7 mm and 1 cm trace emission from a compact disk while millimeter-wave observations trace an extended disk structure. From a physical disk model, where we characterize the disk structure of CY Tau and DoAr 25 at wavelengths shorter than 5 cm, we find that (1) dust continuum emission is optically thin at the observed wavelengths and over the spatial scales studied, (2) a constant value of the dust opacity is not warranted by our observations, and (3) a high-significance radial gradient of the dust opacity spectral index, β, is consistent with the observed dust emission in both disks, with low-β in the inner disk and high-β in the outer disk. Assuming that changes in dust properties arise solely due to changes in the maximum particle size (amax), we constrain radial variations of amax in both disks, from cm-sized particles in the inner disk (R < 40 AU) to millimeter sizes in the outer disk (R > 80 AU). These observational constraints agree with theoretical predictions of the radial-drift barrier, however, fragmentation of dust grains could explain our amax(R) constraints if these disks have lower turbulence and/or if dust can survive high-velocity collisions.

preprint2015arXiv

Large-scale filaments associated with Milky Way spiral arms

The ubiquity of filamentary structure at various scales through out the Galaxy has triggered a renewed interest in their formation, evolution, and role in star formation. The largest filaments can reach up to Galactic scale as part of the spiral arm structure. However, such large scale filaments are hard to identify systematically due to limitations in identifying methodology (i.e., as extinction features). We present a new approach to directly search for the largest, coldest, and densest filaments in the Galaxy, making use of sensitive Herschel Hi-GAL data complemented by spectral line cubes. We present a sample of the 9 most prominent Herschel filaments, including 6 identified from a pilot search field plus 3 from outside the field. These filaments measure 37-99 pc long and 0.6-3.0 pc wide with masses (0.5-8.3)$\times10^4 \, M_\odot$, and beam-averaged ($28"$, or 0.4-0.7 pc) peak H$_2$ column densities of (1.7-9.3)$\times 10^{22} \, \rm{cm^{-2}}$. The bulk of the filaments are relatively cold (17-21 K), while some local clumps have a dust temperature up to 25-47 K. All the filaments are located within <~60 pc from the Galactic mid-plane. Comparing the filaments to a recent spiral arm model incorporating the latest parallax measurements, we find that 7/9 of them reside within arms, but most are close to arm edges. These filaments are comparable in length to the Galactic scale height and therefore are not simply part of a grander turbulent cascade.

preprint2015arXiv

Large-scale latitude distortions of the inner Milky Way Disk from the Herschel/Hi-GAL Survey

We use the Herschel Hi-GAL survey data to study the spatial distribution in Galactic longitude and latitude of the interstellar medium and of dense, star-forming clumps in the inner Galaxy. The peak position and width of the latitude distribution of the dust column density as well as of number density of compact sources from the band-merged Hi-GAL photometric catalogues are analysed as a function of longitude. The width of the diffuse dust column density traced by the Hi-GAL 500 micron emission varies across the inner Galaxy, with a mean value of 1°.2-1°.3, similar to that of the 250um Hi-GAL sources. 70um Hi-GAL sources define a much thinner disk, with a mean FWHM of 0°.75, and an average latitude of b=0°.06, coincident with the results from ATLASGAL. The GLAT distribution as a function of GLON shows modulations, both for the diffuse emission and for the compact sources, with ~0°.2 displacements mostly toward negative latitudes at l~ +40°, +12°, -25° and -40°. No such modulations can be found in the MIPSGAL 24 or WISE 22 um data when the entire source samples are considered. The distortions revealed by Herschel are interpreted as large-scale bending modes of the Plane. The lack of similar distortions in tracers of more evolved YSOs or stars rules out gravitational instabilities or satellite-induced perturbations, as they should act on both the diffuse and stellar disk components. We propose that the observed bends are caused by incoming flows of extra-planar gas interacting with the gaseous disk. Stars decouple from the gaseous ISM and relax into the stellar disk potential. The time required for the disappearance of the distortions from the diffuse ISM to the relatively evolved YSO stages are compatible with star-formation timescales.

preprint2015arXiv

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

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

preprint2014arXiv

CARMA Large Area Star Formation Survey: Structure and Kinematics of Dense Gas in Serpens Main

We present observations of N2H+(1-0), HCO+(1-0), and HCN(1-0) toward the Serpens Main molecular cloud from the CARMA Large Area Star Formation Survey (CLASSy). We mapped 150 square arcminutes of Serpens Main with an angular resolution of 7 arcsecs. The gas emission is concentrated in two subclusters (the NW and SE subclusters). The SE subcluster has more prominent filamentary structures and more complicated kinematics compared to the NW subcluster. The majority of gas in the two subclusters has subsonic to sonic velocity dispersions. We applied a dendrogram technique with N2H+(1-0) to study the gas structures; the SE subcluster has a higher degree of hierarchy than the NW subcluster. Combining the dendrogram and line fitting analyses reveals two distinct relations: a flat relation between nonthermal velocity dispersion and size, and a positive correlation between variation in velocity centroids and size. The two relations imply a characteristic depth of 0.15 pc for the cloud. Furthermore, we have identified six filaments in the SE subcluster. These filaments have lengths of 0.2 pc and widths of 0.03 pc, which is smaller than a characteristic width of 0.1 pc suggested by Herschel observations. The filaments can be classified into two types based on their properties. The first type, located in the northeast of the SE subcluster, has larger velocity gradients, smaller masses, and nearly critical mass-per-unit-length ratios. The other type, located in the southwest of the SE subcluster, has the opposite properties. Several YSOs are formed along two filaments which have supercritical mass per unit length ratios, while filaments with nearly critical mass-per-unit-length ratios are not associated with YSOs, suggesting that stars are formed on gravitationally unstable filaments.

preprint2014arXiv

Constraints on photoevaporation models from (lack of) radio emission in the Corona Australis protoplanetary disks

Photoevaporation due to high-energy stellar photons is thought to be one of the main drivers of protoplanetary disk dispersal. The fully or partially ionized disk surface is expected to produce free-free continuum emission at centimeter (cm) wavelengths that can be routinely detected with interferometers such as the upgraded Very Large Array (VLA). We use deep (rms noise down to 8 $μ$Jy beam$^{-1}$ in the field of view center) 3.5 cm maps of the nearby (130 pc) Corona Australis (CrA) star formation (SF) region to constrain disk photoevaporation models. We find that the radio emission from disk sources in CrA is surprisingly faint. Only 3 out of 10 sources within the field of view are detected, with flux densities of order $10^2$ $μ$Jy. However, a significant fraction of their emission is non-thermal. Typical upper limits for non-detections are $3σ\sim 60~μ$Jy beam$^{-1}$. Assuming analytic expressions for the free-free emission from extreme-UV (EUV) irradiation, we derive stringent upper limits to the ionizing photon luminosity impinging on the disk surface $Φ_\mathrm{EUV}<1-4\times10^{41}$ s$^{-1}$. These limits constrain $Φ_\mathrm{EUV}$ to the low end of the values needed by EUV-driven photoevaporation models to clear protoplanetary disks in the observed few Myr timescale. Therefore, at least in CrA, EUV-driven photoevaporation is not likely to be the main agent of disk dispersal. We also compare the observed X-ray luminosities $L_X$ of disk sources with models in which photoevaporation is driven by such photons. Although predictions are less specific than for the EUV case, most of the observed fluxes (upper limits) are roughly consistent with the (scaled) predictions. Deeper observations, as well as predictions spanning a wider parameter space, are needed to properly test X-ray driven photoevaporation.

preprint2014arXiv

Detectability of Glycine in Solar-type System Precursors

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

preprint2014arXiv

Dust Evolution in Protoplanetary Disks

(abridged) In the core accretion scenario for the formation of planetary rocky cores, the first step toward planet formation is the growth of dust grains into larger and larger aggregates and eventually planetesimals. Although dust grains are thought to grow from the submicron sizes typical of interstellar dust to micron size particles in the dense regions of molecular clouds and cores, the growth from micron size particles to pebbles and kilometre size bodies must occur in protoplanetary disks. This step in the formation of planetary systems is the last stage of solids evolution that can be observed directly in young extrasolar systems. In this chapter we review the constraints on the physics of grain-grain collisions as they have emerged from laboratory experiments and numerical computations. We then review the current theoretical understanding of the global processes governing the evolution of solids in protoplanetary disks, including dust settling, growth, and radial transport. The predicted observational signatures are summarized. We discuss recent developments in the study of grain growth in molecular cloud cores and in collapsing envelopes of protostars as these provide the initial conditions for the dust in disks. We discuss the observational evidence for the growth of grains in young disks from mm surveys, as well as the recent evidence of radial variations of the dust properties in disks. We include a brief discussion of the constraints on the small end of the grain size distribution and on dust settling as derived from optical and IR observations. The observations are discussed in the context of global dust evolution models, in particular we focus on the emerging evidence for a very efficient early growth of grains and the radial distribution of grain sizes in disks. We also highlight the limits of current models, including the need to slow the radial drift of grains.

preprint2014arXiv

Hierarchical fragmentation and differential star formation in the Galactic "Snake": infrared dark cloud G11.11-0.12

We present Submillimeter Array (SMA) $λ=$ 0.88 and 1.3 mm broad band observations, and the Jansky Very Large Array (VLA) observations in $\rm{NH_3}$ $(J,K) = (1,1)$ up to $(5,5)$, $\rm{H_2O}$ and $\rm{CH_3OH}$ maser lines toward the two most massive molecular clumps in infrared dark cloud (IRDC) G11.11-0.12. Sensitive high-resolution images reveal hierarchical fragmentation in dense molecular gas from the $\sim 1$ pc clump scale down to $\sim 0.01$ pc condensation scale. At each scale, the mass of the fragments is orders of magnitude larger than the Jeans mass. This is common to all four IRDC clumps we studied, suggesting that turbulence plays an important role in the early stages of clustered star formation. Masers, shock heated $\rm{NH_3}$ gas, and outflows indicate intense ongoing star formation in some cores while no such signatures are found in others. Furthermore, chemical differentiation may reflect the difference in evolutionary stages among these star formation seeds. We find $\rm{NH_3}$ ortho/para ratios of $1.1\pm0.4$, $2.0\pm0.4$, and $3.0\pm0.7$ associated with three outflows, and the ratio tends to increase along the outflows downstream. Our combined SMA and VLA observations of several IRDC clumps present the most in depth view so far of the early stages prior to the hot core phase, revealing snapshots of physical and chemical properties at various stages along an apparent evolutionary sequence.

preprint2014arXiv

How to Detect the Signatures of Self-Gravitating Circumstellar Discs with the Atacama Large Millimetre/sub-millimetre Array

In this paper we present simulated Atacama Large Millimetre/sub-millimetre Array (ALMA) observations of self-gravitating circumstellar discs with different properties in size, mass and inclination, located in four of the most extensively studied and surveyed star-forming regions. Starting from a Smoothed Particle Hydrodynamics (SPH) simulation and representative dust opacities, we have initially constructed maps of the expected emission at sub-mm wavelengths of a large sample of discs with different properties. We have then simulated realistic observations of discs as they may appear with ALMA using the Common Astronomy Software Application ALMA simulator. We find that, with a proper combination of antenna configuration and integration time, the spiral structure characteristic of self-gravitating discs is readily detectable by ALMA over a wide range of wavelengths at distances comparable to TW Hydrae ($\sim 50 \,$pc), Taurus - Auriga and Ophiucus ($\sim 140 \,$pc) star-forming regions. However, for discs located in Orion complex ($\sim 400 \,$pc) only the largest discs in our sample (outer radius of 100 au) show a spatially resolved structure while the smaller ones (outer radius of 25 au) are characterized by a spiral structure that is not conclusively detectable with ALMA.

preprint2014arXiv

Resolved Multifrequency Radio Observations of GG Tau

We present sub-arcsecond resolution observations of continuum emission associated with the GG Tau quadruple star system at wavelengths of 1.3, 2.8, 7.3, and 50 mm. These data confirm that the GG Tau A binary is encircled by a circumbinary ring at a radius of 235 AU with a FWHM width of ~60 AU. We find no clear evidence for a radial gradient in the spectral shape of the ring, suggesting that the particle size distribution is spatially homogeneous on angular scales of ~0.1". A central point source, likely associated with the primary component (GG Tau Aa), exhibits a composite spectrum from dust and free-free emission. Faint emission at 7.3 mm is observed toward the low-mass star GG Tau Ba, although its origin remains uncertain. Using these measurements of the resolved, multifrequency emission structure of the GG Tau A system, models of the far-infrared to radio spectrum are developed to place constraints on the grain size distribution and dust mass in the circumbinary ring. The non-negligible curvature present in the ring spectrum implies a maximum particle size of 1-10 mm, although we are unable to place strong constraints on the distribution shape. The corresponding dust mass is 30-300 M_earth, at a temperature of 20-30 K. We discuss how this significant concentration of relatively large particles in a narrow ring at a large radius might be produced in a local region of higher gas pressures (i.e., a particle "trap") located near the inner edge of the circumbinary disk.

preprint2014arXiv

Temperaments of young stars: Rapid mass-accretion rate changes in T Tauri and Herbig Ae stars

Variability in emission lines is a characteristic feature in young stars and can be used as a tool to study the physics of the accretion process. Here we present a study of Hα variability in 15 T Tauri and Herbig Ae stars (K7-B2) over a wide range of time windows, from minutes, to hours, to days, and years. We assess the variability using linewidth measurements and the time series of line profiles. All objects show gradual, slow profile changes on time-scales of days. In addition, in three cases there is evidence for rapid variations in Hα with typical time-scales of 10 min, which occurs in 10% of the total covered observing time. The mean accretion-rate changes, inferred from the line fluxes,are 0.01-0.07 dex for time-scales of < 1 hour, 0.04-0.4 dex for time-scales of days, and 0.13-0.52 dex for time-scales of years. In Costigan et al. 2012 we derived an upper limit finding that the intermediate (days) variability dominated over longer (years) variability. Here our new results, based on much higher cadence observations, also provide a lower limit to accretion-rate variability on similar time-scales (days), thereby constraining the accretion rate variability physics in a much more definitive way. A plausible explanation for the gradual variations over days is an asymmetric accretion flow resulting in a rotational modulation of the accretion-related emission, although other interpretations are possible as well. In conjunction with our previous work, we find that the time-scales and the extent of the variability is similar for objects ranging in mass from ~ 0.1 to ~ 5 Msol. This confirms that a single mode of accretion is at work from T Tauri to Herbig Ae stars - across a wide range of stellar masses.

preprint2014arXiv

The Formation and Early Evolution of Young Massive Clusters

We review the formation and early evolution of the most massive and dense young stellar clusters, focusing on the role they can play in our understanding of star and planet formation as a whole. Young massive cluster (YMC) progenitor clouds in the Galactic Center can accumulate to a high enough density without forming stars that the initial protostellar densities are close to the final stellar density. For this to hold in the disk, the time scale to accumulate the gas to such high densities must be much shorter than the star formation timescale. Otherwise the gas begins forming stars while it is being accumulated to high density. The distinction between the formation regimes in the two environments is consistent with the predictions of environmentally-dependent density thresholds for star formation. This implies that stars in YMCs of similar total mass and radius can have formed at widely different initial protostellar densities. The fact that no systematic variations in fundamental properties are observed between YMCs in the disk and Galactic Center suggests stellar mass assembly is not strongly affected by the initial protostellar density. We review recent theoretical advances and summarize the debate on three key open questions: the initial (proto)stellar distribution, infant (im)mortality and age spreads within YMCs. We conclude: the initial protostellar distribution is likely hierarchical; YMCs likely experienced a formation history that was dominated by gas exhaustion rather than gas expulsion; YMCs are dynamically stable from a young age; and YMCs have age spreads much smaller than their mean age. Finally, we show that it is plausible that metal-rich globular clusters may have formed in a similar way to YMCs in nearby galaxies. In summary, the study of YMC formation bridges star/planet formation in the solar neighborhood to the oldest structures in the local Universe. [abridged]

preprint2014arXiv

The impact of the SKA on Galactic Radioastronomy: continuum observations

The SKA will be a state of the art radiotelescope optimized for both large area surveys as well as for deep pointed observations. In this paper we analyze the impact that the SKA will have on Galactic studies, starting from the immense legacy value of the all-sky survey proposed by the continuum SWG but also presenting some areas of Galactic Science that particularly benefit from SKA observations both surveys and pointed. The planned all-sky survey will be characterized by unique spatial resolution, sensitivity and survey speed, providing us with a wide-field atlas of the Galactic continuum emission. Synergies with existing, current and planned radio Galactic Plane surveys will be discussed. SKA will give the opportunity to create a sensitive catalog of discrete Galactic radio sources, most of them representing the interaction of stars at various stages of their evolution with the environment: complete census of all stage of HII regions evolution; complete census of late stages of stellar evolution such as PNe and SNRs; detection of stellar winds, thermal jets, Symbiotic systems, Chemically Peculiar and dMe stars, active binary systems in both flaring and quiescent states. Coherent emission events like Cyclotron Maser in the magnetospheres of different classes of stars can be detected. Pointed, deep observations will allow new insights into the physics of the coronae and plasma processes in active stellar systems and single stars, enabling the detection of flaring activity in larger stellar population for a better comprehension of the mechanism of energy release in the atmospheres of stars with different masses and age.

preprint2014arXiv

The Milky Way as a Star Formation Engine

The cycling of material from the interstellar medium (ISM) into stars and the return of stellar ejecta into the ISM is the engine that drives the "galactic ecology" in normal spirals, a cornerstone in the formation and evolution of galaxies through cosmic time. Major observational and theoretical challenges need to be addressed in determining the processes responsible for converting the low-density ISM into dense molecular clouds, forming dense filaments and clumps, fragmenting them into stars, OB associations and bound clusters, and characterizing the feedback that limits the rate and efficiency of star formation. This formidable task can be now effectively attacked thanks to the combination of new global-scale surveys of the Milky Way Galactic Plane from infrared to radio wavelengths, offering the possibility of bridging the gap between local and extragalactic star formation studies. The Herschel, Spitzer and WISE mid to far infrared continuum surveys, complemented by analogue surveys from ground-based facilities in the millimetre and radio wavelengths, enables us to measure the Galactic distribution and physical properties of dust on all scales and in all components of the ISM from diffuse clouds to filamentary complexes and tens of thousands of dense clumps. A complementary suite of spectroscopic surveys in various atomic and molecular tracers is providing the chemical fingerprinting of dense clumps and filaments, as well as essential kinematic information to derive distances and thus transform panoramic data into a 3D representation. The latest results emerging from these Galaxy-scale surveys are reviewed. New insights into cloud formation and evolution, filaments and their relationship to channeling gas onto gravitationally-bound clumps, the properties of these clumps, density thresholds for gravitational collapse, and star and cluster formation rates are discussed.

preprint2014arXiv

What controls star formation in the central 500 pc of the Galaxy?

The star formation rate (SFR) in the Central Molecular Zone (CMZ, i.e. the central 500 pc) of the Milky Way is lower by a factor of >10 than expected for the substantial amount of dense gas it contains, which challenges current star formation theories. In this paper, we quantify which physical mechanisms could be responsible. On scales larger than the disc scale height, the low SFR is found to be consistent with episodic star formation due to secular instabilities or possibly variations of the gas inflow along the Galactic bar. The CMZ is marginally Toomre-stable when including gas and stars, but highly Toomre-stable when only accounting for the gas, indicating a low condensation rate of self-gravitating clouds. On small scales, we find that the SFR in the CMZ may be caused by an elevated critical density for star formation due to the high turbulent pressure. The existence of a universal density threshold for star formation is ruled out. The HI-H$_2$ phase transition of hydrogen, the tidal field, a possible underproduction of massive stars due to a bottom-heavy initial mass function, magnetic fields, and cosmic ray or radiation pressure feedback also cannot individually explain the low SFR. We propose a self-consistent cycle of star formation in the CMZ, in which the effects of several different processes combine to inhibit star formation. The rate-limiting factor is the slow evolution of the gas towards collapse - once star formation is initiated it proceeds at a normal rate. The ubiquity of star formation inhibitors suggests that a lowered central SFR should be a common phenomenon in other galaxies. We discuss the implications for galactic-scale star formation and supermassive black hole growth, and relate our results to the star formation conditions in other extreme environments.

preprint2013arXiv

A UV-to-MIR monitoring of DR Tau: exploring how water vapor in the planet formation region of the disk is affected by stellar accretion variability

Young stars are known to show variability due to non-steady mass accretion rate from their circumstellar disks. Accretion flares can produce strong energetic irradiation and heating that may affect the disk in the planet formation region, close to the central star. During an extreme accretion outburst in the young star EX Lupi, the prototype of EXor variables, remarkable changes in molecular gas emission from $\sim1$ AU in the disk have recently been observed (Banzatti et al. 2012). Here, we focus on water vapor and explore how it is affected by variable accretion luminosity in T Tauri stars. We monitored a young highly variable solar-mass star, DR Tau, using simultaneously two high/medium-resolution ESO-VLT spectrographs: VISIR at 12.4 $μ$m to observe water lines from the disk, and X-shooter covering from 0.3 to 2.5 $μ$m to constrain the stellar accretion. Three epochs spanning timescales from several days to several weeks were obtained. Accretion luminosity was estimated to change within a factor $\sim2$, and no change in water emission was detected at a significant level. In comparison to EX Lupi and EXor outbursts, DR Tau suggests that the less long-lived and weaker variability phenomena typical of T Tauri stars may leave water at planet-forming radii in the disk mostly unaffected. We propose that these systems may provide evidence for two processes that act over different timescales: UV photochemistry in the disk atmosphere (faster) and heating of the disk deeper layers (slower).

preprint2012arXiv

Changes Of Dust Opacity With Density in the Orion A Molecular Cloud

We have studied the opacity of dust grains at submillimeter wavelengths by estimating the optical depth from imaging at 160, 250, 350, and 500 um from the Herschel Gould Belt Survey and comparing this to a column density obtained from the 2MASS-derived color excess E(J-Ks). Our main goal was to investigate the spatial variations of the opacity due to "big" grains over a variety of environmental conditions and thereby quantify how emission properties of the dust change with column (and volume) density. The central and southern areas of the Orion A molecular cloud examined here, with NH ranging from 1.5X10^21 cm^-2 to 50X10^21 cm^-2, are well suited to this approach. We fit the multi-frequency Herschel spectral energy distributions (SEDs) of each pixel with a modified blackbody to obtain the temperature, T, and optical depth, τ(1200), at a fiducial frequency of 1200 GHz (250 um). Using a calibration of NH/E(J-Ks)for the interstellar medium (ISM) we obtained the opacity (dust emission cross-section per H nucleon), σ_e(1200), for every pixel. From a value of ~ 1X10^-25 cm^2 H^-1 at the lowest column densities that is typical of the high latitude diffuse ISM, σ_e(1200) increases as NH^0.28 over the range studied. This is suggestive of grain evolution. Integrating the SEDs over frequency, we also calculated the specific power P (emission power per H) for the big grains. In low column density regions where dust clouds are optically thin to the interstellar radiation field (ISRF), P is typically 3.7 X 10^-31 W H^-1, again close to that in the high latitude diffuse ISM. However, we find evidence for a decrease of P in high column density regions, which would be a natural outcome of attenuation of the ISRF that heats the grains, and for localized increases for dust illuminated by nearby stars or embedded protostars.

preprint2012arXiv

Circumstellar disks and planets. Science cases for next-generation optical/infrared long-baseline interferometers

We present a review of the interplay between the evolution of circumstellar disks and the formation of planets, both from the perspective of theoretical models and dedicated observations. Based on this, we identify and discuss fundamental questions concerning the formation and evolution of circumstellar disks and planets which can be addressed in the near future with optical and infrared long-baseline interferometers. Furthermore, the importance of complementary observations with long-baseline (sub)millimeter interferometers and high-sensitivity infrared observatories is outlined.

preprint2012arXiv

Constraints on the Radial Variation of Grain Growth in the AS 209 Circumstellar Disk

We present dust continuum observations of the protoplanetary disk surrounding the pre-main sequence star AS 209, spanning more than an order of magnitude in wavelength from 0.88 to 9.8 mm. The disk was observed with sub-arcsecond angular resolution (0.2"-0.5") to investigate radial variations in its dust properties. At longer wavelengths, the disk emission structure is notably more compact, providing model-independent evidence for changes in the grain properties across the disk. We find that physical models which reproduce the disk emission require a radial dependence of the dust opacity κ_ν. Assuming that the observed wavelength-dependent structure can be attributed to radial variations in the dust opacity spectral index (β), we find that β(R) increases from β<0.5 at \sim20 AU to β>1.5 for R>80 AU, inconsistent with a constant value of β across the disk (at the 10σ level). Furthermore, if radial variations of κ_ν are caused by particle growth, we find that the maximum size of the particle-size distribution (a_{max}) increases from sub-millimeter-sized grains in the outer disk (R>70 AU) to millimeter and centimeter-sized grains in the inner disk regions (R< 70 AU). We compare our observational constraint on a_{max}(R) with predictions from physical models of dust evolution in proto-planetary disks. For the dust composition and particle-size distribution investigated here, our observational constraints on a_{max}(R) are consistent with models where the maximum grain size is limited by radial drift.

preprint2012arXiv

How to Find Young Massive Cluster Progenitors

We propose that bound, young massive stellar clusters form from dense clouds that have escape speeds greater than the sound speed in photo-ionized gas. In these clumps, radiative feedback in the form of gas ionization is bottled up, enabling star formation to proceed to sufficiently high efficiency so that the resulting star cluster remains bound even after gas removal. We estimate the observable properties of the massive proto-clusters (MPCs) for existing Galactic plane surveys and suggest how they may be sought in recent and upcoming extragalactic observations. These surveys will potentially provide a significant sample of MPC candidates that will allow us to better understand extreme star-formation and massive cluster formation in the Local Universe.

preprint2012arXiv

Modes of star formation from Herschel

We summarize some of the results obtained from Herschel surveys of the nearby star forming regions and the Galactic plane. We show that in the nearby star forming regions the starless core spatial surface density distribution is very similar to that of the young stellar objects. This, taken together with the similarity between the core mass function and the initial mass function for stars and the relationship between the amount of dense gas and star formation rate, suggest that the cloud fragmentation process defines the global outcome of star formation. This "simple" view of star formation may not hold on all scales. In particular dynamical interactions are expected to become important at the conditions required to form young massive clusters. We describe the successes of a simple criterion to identify young massive cluster precursors in our Galaxy based on (sub-)millimetre wide area surveys. We further show that in the location of our Galaxy where the best candidate for a precursor of a young massive cluster is found, the "simple" scaling relationship between dense gas and star formation rate appear to break down. We suggest that in regions where the conditions approach those of the central molecular zone of our Galaxy it may be necessary to revise the scaling laws for star formation.

preprint2012arXiv

The first ALMA view of IRAS 16293-2422: Direct detection of infall onto source B and high-resolution kinematics of source A

Aims: We focus on the kinematical properties of a proto-binary to study the infall and rotation of gas towards its two protostellar components. Methods: We present ALMA Science Verification observations with high-spectral resolution of IRAS 16293-2422 at 220.2 GHz. The wealth of molecular lines in this source and the very high spectral resolution offered by ALMA allow us to study the gas kinematics with unprecedented detail. Results: We present the first detection of an inverse P-Cygni profile towards source B in the three brightest lines. The line profiles are fitted with a simple two-layer model to derive an infall rate of 4.5x10^-5 Msun/yr. This infall detection would rule-out the previously suggested possibility that source B is a T Tauri star. A position velocity diagram for source A shows evidence for rotation with an axis close to the line-of-sight.

preprint2012arXiv

X-shooter spectroscopy of young stellar objects: I - Mass accretion rates of low-mass T Tauri stars in σOrionis

We present high-quality, medium resolution X-shooter/VLT spectra in the range 300-2500 nm for a sample of 12 very low-mass stars in the σOrionis cluster. The sample includes stars with masses ranging from 0.08 to 0.3 M$_\odot$. The aim of this first paper is to investigate the reliability of the many accretion tracers currently used to measure the mass accretion rate in low-mass, young stars. We use our spectra to measure the accretion luminosity from the continuum excess emission in the UV and visual; the derived mass accretion rates range from 10$^{-9}$ M$_{\odot}$ yr$^{-1}$ down to 5$\times10^{-11}$ M$_{\odot}$ yr$^{-1}$, allowing us to investigate the behavior of the accretion-driven emission lines in very-low mass accretion rate regimes. We compute the luminosity of ten accretion-driven emission lines, from the UV to the near-IR, obtained simultaneously. Most of the secondary tracers correlate well with the accretion luminosity derived from the continuum excess emission. We confirm the validity of the correlations between accretion luminosities and line luminosities given in the literature, with the possible exception of Hα. When looking at individual objects, we find that the Hydrogen recombination lines, from the UV to the near-IR, give good and consistent measurements of accretion luminosities, often in better agreement than the uncertainties introduced by the adopted correlations. The average accretion luminosity derived from several Hydrogen lines, measured simultaneously, have a much reduced error. This suggests that some of the spread in the literature correlations may be due to the use of non-simultaneous observations of lines and continuum. Three stars in our sample deviate from this behavior, and we discuss them individually.

preprint2011arXiv

Filaments and ridges in Vela C revealed by Herschel: from low-mass to high-mass star-forming sites

We present the first Herschel PACS and SPIRE results of the Vela C molecular complex in the far-infrared and submillimetre regimes at 70, 160, 250, 350, and 500 um, spanning the peak of emission of cold prestellar or protostellar cores. Column density and multi-resolution analysis (MRA) differentiates the Vela C complex into five distinct sub-regions. Each sub-region displays differences in their column density and temperature probability distribution functions (PDFs), in particular, the PDFs of the `Centre-Ridge' and `South-Nest' sub-regions appear in stark contrast to each other. The Centre-Ridge displays a bimodal temperature PDF representative of hot gas surrounding the HII region RCW 36 and the cold neighbouring filaments, whilst the South-Nest is dominated by cold filamentary structure. The column density PDF of the Centre-Ridge is flatter than the South-Nest, with a high column density tail, consistent with formation through large-scale flows, and regulation by self-gravity. At small to intermediate scales MRA indicates the Centre-Ridge to be twice as concentrated as the South-Nest, whilst on larger scales, a greater portion of the gas in the South-Nest is dominated by turbulence than in the Centre-Ridge. In Vela C, high-mass stars appear to be preferentially forming in ridges, i.e., dominant high column density filaments.

preprint2011arXiv

G0.253+0.016: a molecular cloud progenitor of an Arches-like cluster

Young massive clusters (YMCs) with stellar masses of 10^4 - 10^5 Msun and core stellar densities of 10^4 - 10^5 stars per cubic pc are thought to be the `missing link' between open clusters and extreme extragalactic super star clusters and globular clusters. As such, studying the initial conditions of YMCs offers an opportunity to test cluster formation models across the full cluster mass range. G0.253+0.016 is an excellent candidate YMC progenitor. We make use of existing multi-wavelength data including recently available far-IR continuum (Herschel/Hi-GAL) and mm spectral line (HOPS and MALT90) data and present new, deep, multiple-filter, near-IR (VLT/NACO) observations to study G0.253+0.016. These data show G0.253+0.016 is a high mass (1.3x10^5 Msun), low temperature (T_dust~20K), high volume and column density (n ~ 8x10^4 cm^-3; N_{H_2} ~ 4x10^23 cm^-2) molecular clump which is close to virial equilibrium (M_dust ~ M_virial) so is likely to be gravitationally-bound. It is almost devoid of star formation and, thus, has exactly the properties expected for the initial conditions of a clump that may form an Arches-like massive cluster. We compare the properties of G0.253+0.016 to typical Galactic cluster-forming molecular clumps and find it is extreme, and possibly unique in the Galaxy. This uniqueness makes detailed studies of G0.253+0.016 extremely important for testing massive cluster formation models.

preprint2011arXiv

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

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

preprint2010arXiv

A hot compact dust disk around a massive young stellar object

Circumstellar disks are an essential ingredient of the formation of low-mass stars. It is unclear, however, whether the accretion-disk paradigm can also account for the formation of stars more massive than about 10 solar masses, in which strong radiation pressure might halt mass infall. Massive stars may form by stellar merging, although more recent theoretical investigations suggest that the radiative-pressure limit may be overcome by considering more complex, nonspherical infall geometries. Clear observational evidence, such as the detection of compact dusty disks around massive young stellar objects, is needed to identify unambiguously the formation mode of the most massive stars. Here we report near-infrared interferometric observations that spatially resolve the astronomical unit-scale distribution of hot material around a high-mass (approx. 20 solar masses) young stellar object. The image shows an elongated structure with a size of about 13 x 19 astronomical units, consistent with a disk seen at an inclination angle of 45 degree. Using geometric and detailed physical models, we found a radial temperature gradient in the disk, with a dust-free region less than 9.5 astronomical units from the star, qualitatively and quantitatively similar to the disks observed in low-mass star formation. Perpendicular to the disk plane we observed a molecular outflow and two bow shocks, indicating that a bipolar outflow emanates from the inner regions of the system.

preprint2010arXiv

First results from VLTI near-infrared interferometry on high-mass young stellar objects

Due to the recent dramatic technological advances, infrared interferometry can now be applied to new classes of objects, resulting in exciting new science prospects, for instance, in the area of high-mass star formation. Although extensively studied at various wavelengths, the process through which massive stars form is still only poorly understood. For instance, it has been proposed that massive stars might form like low-mass stars by mass accretion through a circumstellar disk/envelope, or otherwise by coalescence in a dense stellar cluster. After discussing the technological challenges which result from the special properties of these objects, we present first near-infrared interferometric observations, which we obtained on the massive YSO IRAS 13481-6124 using VLTI/AMBER infrared long-baseline interferometry and NTT speckle interferometry. From our extensive data set, we reconstruct a model-independent aperture synthesis image which shows an elongated structure with a size of 13x19 AU, consistent with a disk seen under an inclination of 45 degree. The measured wavelength-dependent visibilities and closure phases allow us to derive the radial disk temperature gradient and to detect a dust-free region inside of 9.5 AU from the star, revealing qualitative and quantitative similarities with the disks observed in low-mass star formation. In complementary mid-infrared Spitzer and sub-millimeter APEX imaging observations we detect two bow shocks and a molecular out ow which are oriented perpendicular to the disk plane and indicate the presence of a bipolar outflow emanating from the inner regions of the system.

preprint2010arXiv

Millimeter imaging of MWC 758: probing the disk structure and kinematics

We investigate the structure and kinematics of the circumstellar disk around the Herbig Ae star MWC 758 using high resolution observations of the 12CO (3-2) and dust continuum emission at the wavelengths of 0.87 and 3.3 mm. We find that the dust emission peaks at an orbital radius of about 100 AU, while the CO intensity has a central peak coincident with the position of the star. The CO emission is in agreement with a disk in keplerian rotation around a 2.0 Msun star, confirming that MWC758 is indeed an intermediate mass star. By comparing the observation with theoretical disk models, we derive that the disk surface density Sigma(r) steeply increases from 40 to 100 AU, and decreases exponentially outward. Within 40 AU, the disk has to be optically thin in the continuum emission at millimeter wavelengths to explain the observed dust morphology, though our observations lack the angular resolution and sensitivity required to constrain the surface density on these spatial scales. The surface density distribution in MWC 758 disk is similar to that of ``transition'' disks, though no disk clearing has been previously inferred from the analysis of the spectral energy distribution (SED). Moreover, the asymmetries observed in the dust and CO emission suggest that the disk may be gravitationally perturbed by a low mass companion orbiting within a radius of 30 AU. Our results emphasize that SEDs alone do not provide a complete picture of disk structure and that high resolution millimeter-wave images are essential to reveal the structure of the cool disk mid plane.

preprint2010arXiv

Testing circumstellar disk lifetimes in young embedded clusters associated with the Vela Molecular Ridge

Context. The Vela Molecular Ridge hosts a number of young embedded star clusters in the same evolutionary stage. Aims. The main aim of the present work is testing whether the fraction of members with a circumstellar disk in a sample of clusters in the cloud D of the Vela Molecular Ridge, is consistent with relations derived for larger samples of star clusters with an age spread. Besides, we want to constrain the age of the young embedded star clusters associated with cloud D. Methods. We carried out L (3.78 microns) photometry on images of six young embedded star clusters associated with cloud D of the Vela Molecular Ridge, taken with ISAAC at the VLT. These data are complemented with the available HKs photometry. The 6 clusters are roughly of the same size and appear to be in the same evolutionary stage. The fraction of stars with a circumstellar disk was measured in each cluster by counting the fraction of sources displaying a NIR excess in colour-colour (HKsL) diagrams. Results. The L photometry allowed us to identify the NIR counterparts of the IRAS sources associated with the clusters. The fraction of stars with a circumstellar disk appears to be constant within errors for the 6 clusters. There is a hint that this is lower for the most massive stars. The age of the clusters is constrained to ~1-2 Myr. Conclusions. The fraction of stars with a circumstellar disk in the observed sample is consistent with the relations derived from larger samples of star clusters and with other age estimates for cloud D. The fraction may be lower for the most massive stars. Our results agree with a scenario where all intermediate and low-mass stars form with a disk, whose lifetime is shorter for higher mass stars.

preprint2010arXiv

The near-infrared reflected spectrum of source I in Orion-KL

Source I in the Orion-KL nebula is believed to be the nearest example of a massive star still in the main accretion phase. It is thus one of the best cases for studying the properties of massive protostars to constrain high-mass star formation theories. Near-infrared radiation from source I escapes through the cavity opened by the OMC1 outflow and is scattered by dust towards our line of sight. The reflected spectrum offers a unique possibility of observing the emission from the innermost regions of the system and probing the nature of source I and its immediate surroundings. We obtained moderately high spectral-resolution (R~9000) observations of the near-infrared diffuse emission in several locations around source I/Orion-KL. We observed a widespread rich absorption line spectrum that we compare with cool stellar photospheres and protostellar accretion disk models. The spectrum is broadly similar to strongly veiled, cool, low-gravity stellar photospheres in the range Teff~3500-4500 K, luminosity class I-III. An exact match explaining all features has not been found, and a plausible explanation is that a range of different temperatures contribute to the observed absorption spectrum. The 1D velocity dispersions implied by the absorption spectra, sigma~30 km/s, can be explained by the emission from a disk around a massive, mstar~10 Msun, protostar that is accreting at a high rate, mdot~3x10^{-3} Msun/yr. Our observations suggest that the near-infrared reflection spectrum observed in the Orion-KL region is produced close to source I and scattered to our line of sight in the OMC1 outflow cavity. The spectrum allows us to exclude source I being a very large, massive protostar rotating at breakup speed. We suggest that the absorption spectrum is produced in a disk surrounding a ~10 Msun protostar, accreting from its disk at a high rate of a few 10^{-3} Msun/yr.

preprint2009arXiv

Properties of Stellar Clusters around High-Mass Young Stars

[Abridged] Twenty-six high-luminosity IRAS sources believed to be collection of stars in the early phases of high-mass star formation have been observed in the NIR (J, H, K) to characterize the clustering properties of their young stellar population and gain insight into the initial conditions of star formation in these clusters (Initial Mass Function [IMF], Star Formation History [SFH]), and to deduce mean values for cluster ages. K luminosity functions (KLFs) are compared with simulated ones from a model that generates populations of synthetic clusters starting from assumptions on the IMF, the SFH, and the Pre-MS evolution, and using the average properties of the observed clusters as boundary conditions Twenty-two sources show evidence of clustering from a few up to several tens of objects, and a median cluster radius of 0.7 pc. A considerable number of cluster members present an infrared excess characteristic of young Pre-Main-Sequence objects. We find that the median stellar age ranges between 2.5 10^5 and 5 10^6 years, with evidence of an age spread of the same entity within each cluster. We also find evidence that older clusters tend to be smaller in size, in line with the fact that our clusters are on average larger than those around relatively older Herbig Ae/Be stars. The relationship of the mass of the most massive star in the cluster with both the clusters richness and their total stellar mass suggest that our modeled clusters may not be consistent with them resulting from random sampling of the IMF. Our results are consistent with a star formation which takes place continuously over a period of time which is longer than a typical crossing time.

preprint2008arXiv

The star formation in the L1615/L1616 cometary cloud

The present work aims at performing a comprehensive census and characterisation of the pre-main sequence (PMS) population in the cometary cloud L1615/L1616, in order to assess the significance of the triggered star formation scenario and investigate the impact of massive stars on its star formation history and mass spectrum. Our study is based on UBVRcIc and JHKs photometry, as well as optical multi-object spectroscopy. We performed a physical parametrisation of the young stellar population in L1615/L1616. We identified 25 new T Tauri stars mainly projected on the dense head of the cometary cloud, almost doubling the current number of known members. We studied the spatial distribution of the cloud members as a function of the age and H$α$ emission. The star formation efficiency in the cloud is about 7-8 %, as expected for molecular clouds in the vicinity of OB associations. The slope of the initial mass function (IMF), in the mass range 0.1<M<5.5 $M_{\odot}$, is consistent with that of other T and OB associations, providing further support of an universal IMF down to the hydrogen burning limit, regardless of environmental conditions. The cometary appearance, as well as the high star formation efficiency, can be explained in terms of triggered star formation induced by the strong UV radiation from OB stars or supernovae shockwaves. The age spread as well as both the spatial and age distribution of the PMS objects provide strong evidence of sequential, multiple events and possibly still ongoing star formation activity in the cloud.

preprint2005arXiv

Search for massive protostar candidates in the southern hemisphere: II. Dust continuum emission

In an ongoing effort to identify and study high-mass protostellar candidates we have observed in various tracers a sample of 235 sources selected from the IRAS Point Source Catalog, mostly with dec < -30 deg, with the SEST antenna at millimeter wavelengths. The sample contains 142 Low sources and 93 High, which are believed to be in different evolutionary stages. Both sub-samples have been studied in detail by comparing their physical properties and morphologies. Massive dust clumps have been detected in all but 8 regions, with usually more than one clump per region. The dust emission shows a variety of complex morphologies, sometimes with multiple clumps forming filaments or clusters. The mean clump has a linear size of ~0.5 pc, a mass of ~320 Msolar for a dust temperature Td=30 K, an H_2 density of 9.5E5 cm-3, and a surface density of 0.4 g cm-2. The median values are 0.4 pc, 102 Msolar, 4E4 cm-3, and 0.14 g cm-2, respectively. The mean value of the luminosity-to-mass ratio, L/M ~99 Lsolar/Msolar, suggests that the sources are in a young, pre-ultracompact HII phase. We have compared the millimeter continuum maps with images of the mid-IR MSX emission, and have discovered 95 massive millimeter clumps non-MSX emitters, either diffuse or point-like, that are potential prestellar or precluster cores. The physical properties of these clumps are similar to those of the others, apart from the mass that is ~3 times lower than for clumps with MSX counterpart. Such a difference could be due to the potential prestellar clumps having a lower dust temperature. The mass spectrum of the clumps with masses above M ~100 Msolar is best fitted with a power-law dN/dM proportional to M-alpha with alpha=2.1, consistent with the Salpeter (1955) stellar IMF, with alpha=2.35.

preprint2005arXiv

VITRUV - Science Cases

VITRUV is a second generation spectro-imager for the PRIMA enabled Very Large Telescope Interferometer. By combining simultaneously up to 8 telescopes VITRUV makes the VLTI up to 6 times more efficient. This operational gain allows two novel scientific methodologies: 1) massive surveys of sizes; 2) routine interferometric imaging. The science cases presented concentrate on the qualitatively new routine interferometric imaging methodology. The science cases are not exhaustive but complementary to the PRIMA reference mission. The focus is on: a) the close environment of young stars probing for the initial conditions of planet formation and disk evolution; b) the surfaces of stars tackling dynamos, activity, pulsation, mass-loss and evolution; c) revealing the origin of the extraordinary morphologies of Planetary Nebulae and related stars; d) studying the accretion-ejection structures of stellar black-holes (microquasars) in our galaxy; e) unveiling the different interacting components (torus, jets, BLRs) of Active Galactic Nuclei; and f) probing the environment of nearby supermassive black-holes and relativistic effects in the Galactic Center black-hole.