Search

Protoplanetary disks are routinely described as finite mass reservoirs left over by the gravitational collapse of the protostar, an assumption that strongly constrains both disk evolution and planet formation models. We propose a different scenario where protoplanetary disks of pre-main sequence stars are assembled primarily by Bondi-Hoyle accretion from the parent gas cloud. We demonstrate that Bondi-Hoyle accretion can supply not only the mass, but also the angular momentum necessary to explain the observed size of protoplanetary disks, and we predict the dependence of the disk specific angular momentum on the stellar mass. Our results are based on an analytical derivation of the scaling of the angular momentum in a turbulent flow, which we also confirm with a numerical simulation of supersonic turbulence. This new scenario for disk formation and evolution may alleviate a number of observational problems as well as compel major revisions of disk and planet formation models., Comment: Submitted to Nature Astronomy

Context. While Class II stars have already accreted most of their mass, the continued inflow of fresh material via Bondi-Hoyle accretion acts as an additional mass reservoir for their circumstellar disks. This may explain the observed accretion rates of pre-main- sequence (PMS) stars, as well as observational inconsistencies in the mass and angular momentum balance of their disks. Aims. Using a new simulation that reproduces the stellar IMF, we want to quantify the role of Bondi-Hoyle accretion in the formation of Class II disks, as well address the prospect of its observational detection with the James Webb Space Telescope (JWST). Methods. We study the mass and angular momentum of the accreting gas using passively-advected tracer particles in the simulation, and we carry out radiative transfer calculations of near-infrared scattering to generate synthetic JWST observations of Bondi-Hoyle trails of PMS stars. Results. Gas accreting on Class II PMS stars approximately 1 Myr after their formation has enough mass and angular momentum to strongly affect the evolution of the preexisting disks. The accreted angular momentum is large enough to also explain the observed size of Class II disks. The orientation of the angular momentum vector can differ significantly from that of the previously accreted gas, which may result in a significant disk warping or misalignment. We also predict that JWST observations of Class II stars will be able to detect Bondi-Hoyle trails with a 90% success rate with only 2 min exposure time, if stars with accretion rates \dot{M} > 5e-10 Msol/yr and luminosity of L > 0.5 Lsol are selected., Comment: 12 pages, 16 figures, submitted to A&A

The kinematic information of the Radcliffe Wave (RW) is essential for determining its existence and gaining insights into its origin and evolution. In this work, we present an accurate measurement of the vertical velocity ($V_Z$) of RW by incorporating the radial velocity (RV) measures through two methods, which is crucial but was neglected previously. First, the velocities are measured towards young stars, using their RV measurements from APOGEE-2 and proper motion measurements from Gaia DR3. Second, we combine RV measurements toward clouds with proper motion measurements of associated Young Stellar Objects (YSOs) to determine the vertical velocities of the clouds. The results reveal that the oscillations in $V_Z$ are not synchronous with the vertical coordinate $Z$, which differs from the conclusions of previous studies. Instead, we find a 5 km$\cdot$s$^{-1}$$\cdot$kpc$^{-1}$ gradient in $V_Z$ along the RW, exhibiting a dipole-like pattern. Consequently, the kinematic arrangement does not show a corresponding coherence with the spatial arrangement, bringing the Radcliffe Wave model into question., Comment: 9 pages, 6 figures, Accepted for publication in ApJ

Characterizing prestellar cores in star-forming regions is an important step towards the validation of theoretical models of star formation. Thanks to their sub-arcsecond resolution, ALMA observations can potentially provide samples of prestellar cores up to distances of a few kpc, where regions of massive star formation can be targeted. However, the extraction of real cores from dust-continuum observations of turbulent star-forming clouds is affected by complex projection effects. In this work, we study the problem of core extraction both in the idealized case of column-density maps and in the more realistic case of synthetic 1.3\,mm ALMA observations. The analysis is carried out on 12 regions of high column density from our 250 pc simulation. We find that derived core masses are highly unreliable, with only {\em a weak correlation between the masses of cores selected in the synthetic ALMA maps and those of the corresponding three-dimensional cores}. The fraction of real three-dimensional cores detected in the synthetic maps increases monotonically with mass and remains always below 50\%. Above $\sim 1\,M_{\odot}$, the core mass function derived from the column-density maps is steeper than that of the three-dimensional cores, while the core mass function from the synthetic ALMA maps has a slope closer to that of the real three-dimensional cores. Because of the mass uncertainties, proper guidance from realistic simulations is essential if ALMA observations of protoclusters at kpc distances are to be used to test star-formation models., Comment: 20 pages, 21 figures, MNRAS, in press

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

The dynamical state of massive clumps is key to our understanding of the formation of massive stars. In this work, we study the kinematic properties of massive clumps using synthetic observations. We have previously compiled a very large catalog of synthetic dust-continuum compact sources from our 250 pc, SN-driven, star formation simulation. Here, we compute synthetic $\rm N_{2}H^{+}$ line profiles for a subsample of those sources and compare their properties with the observations and with those of the corresponding three-dimensional (3D) clumps in the simulation. We find that the velocity dispersion of the sources estimated from the $\rm N_{2}H^{+}$ line is a good estimate of that of the 3D clumps, although its correlation with the source size is weaker than the velocity-size correlation of the 3D clumps. The relation between the mass of the 3D clumps, $M_{\rm main}$, and that of the corresponding synthetic sources, $M_{\rm SED}$, has a large scatter and a slope of 0.5, $M_{\rm main} \propto M_{\rm SED}^{0.5}$, due to uncertainties arising from the observational band-merging procedure and from projection effects along the line of sight. As a result, the virial parameters of the 3D clumps are not correlated with the clump masses, even if a negative correlation is found for the compact sources, and the virial parameter of the most massive sources may significantly underestimate that of the associated clumps., Comment: 19 pages, 12 figures, 1 table. Accepted for publication in MNRAS

Systematic surveys of massive clumps have been carried out to study the conditions leading to the formation of massive stars. These clumps are typically at large distances and unresolved, so their physical properties cannot be reliably derived from the observations alone. Numerical simulations are needed to interpret the observations. To this end, we generate synthetic Herschel observations using our large-scale star-formation simulation, where massive stars explode as supernovae driving the interstellar-medium turbulence. From the synthetic observations, we compile a catalog of compact sources following the exact same procedure as for the Hi-GAL compact source catalog. We show that the sources from the simulation have observational properties with statistical distributions consistent with the observations. By relating the compact sources from the synthetic observations to their three-dimensional counterparts in the simulation, we find that the synthetic observations overestimate the clump masses by about an order of magnitude on average due to line-of-sight projection, and projection effects are likely to be even worse for Hi-GAL Inner Galaxy sources. We also find that a large fraction of sources classified as protostellar are likely to be starless, and propose a new method to partially discriminate between true and false protostellar sources., Comment: 17 pages, 15 figures, submitted to MNRAS

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

Planck Galactic Cold Clumps (PGCCs) are contemplated to be the ideal targets to probe the early phases of star formation. We have conducted a survey of 72 young dense cores inside PGCCs in the Orion complex with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.3\,mm (band 6) using three different configurations (resolutions $\sim$ 0$\farcs$35, 1$\farcs$0, and 7$\farcs$0) to statistically investigate their evolutionary stages and sub-structures. We have obtained images of the 1.3\,mm continuum and molecular line emission ($^{12}$CO, and SiO) at an angular resolution of $\sim$ 0$\farcs$35 ($\sim$ 140\,au) with the combined arrays. We find 70 substructures within 48 detected dense cores with median dust-mass $\sim$ 0.093\,M$_{\sun}$ and deconvolved size $\sim$ 0$\farcs$27. Dense substructures are clearly detected within the central 1000\,au of four candidate prestellar cores. The sizes and masses of the substructures in continuum emission are found to be significantly reduced with protostellar evolution from Class\,0 to Class\,I. We also study the evolutionary change in the outflow characteristics through the course of protostellar mass accretion. A total of 37 sources exhibit CO outflows, and 20 ($>$50\%) show high-velocity jets in SiO. The CO velocity-extents ($\Delta$Vs) span from 4 to 110 km/s with outflow cavity opening angle width at 400\,au ranging from $[\Theta_{obs}]_{400}$ $\sim$ 0$\farcs$6 to 3$\farcs$9, which corresponds to 33$\fdg$4$-$125$\fdg$7. For the majority of the outflow sources, the $\Delta$Vs show a positive correlation with $[\Theta_{obs}]_{400}$, suggesting that as protostars undergo gravitational collapse, the cavity opening of a protostellar outflow widens and the protostars possibly generate more energetic outflows., Comment: 43 pages, 25 figures (14 in the main text, 9 in the appendix), accepted for publication in the Astrophysical Journal Supplement Series (ApJS)

We study the near-infrared (NIR) scattering in LDN 1642, its correlation with the cloud structure, and the ability of dust models to simultaneously explain sub-millimetre emission, NIR extinction, and NIR scattering. We use observations from the HAWK-I instrument to measure the NIR surface brightness and extinction. These are compared with Herschel data on dust emission and, with radiative transfer modelling, with predictions calculated for different dust models. We find an optical depth ratio $\tau(250\,\mu{\rm m})/\tau(J)\approx 10^{-3}$, confirming earlier findings of high sub-millimetre emissivity. The relationships between the column density derived from dust emission and the NIR colour excesses is linear and consistent with the standard NIR extinction curve. The extinction peaks at $A_J=2.6\,$mag, the NIR surface brightness remaining correlated with $N({\rm H}_2)$ without saturation. Radiative transfer models can fit the sub-millimetre data with any of the tested dust models. However, these predict a NIR extinction that is higher and a NIR surface brightness that is lower than in observations. If the dust sub-millimetre emissivity is rescaled to the observed value of $\tau(250\,\mu{\rm m})/\tau(J)$, dust models with high NIR albedo can reach the observed level of NIR surface brightness. The NIR extinction of the models tends to be higher than directly measured, which is reflected in the shape of the NIR surface brightness spectra. The combination of emission, extinction, and scattering measurements provides strong constraints on dust models. The observations of LDN 1642 indicate clear dust evolution, including a strong increase in the sub-millimetre emissivity, not yet fully explained by the current dust models., Comment: 16 pages, 16 figures, submitted to A&A

While the importance of supernova feedback in galaxies is well established, its role on the scale of molecular clouds is still debated. In this work, we focus on the impact of supernovae on individual clouds, using a high-resolution magneto-hydrodynamic simulation of a region of 250 pc where we resolve the formation of individual massive stars. The supernova feedback is implemented with real supernovae that are the natural evolution of the resolved massive stars, so their position and timing are self-consistent. We select a large sample of molecular clouds from the simulation to investigate the supernova energy injection and the resulting properties of molecular clouds. We find that molecular clouds have a lifetime of a few dynamical times, less then half of them contract to the point of becoming gravitationally bound, and the dispersal time of bound clouds, of order one dynamical time, is a factor of two shorter than that of unbound clouds. We stress the importance of internal supernovae, that is massive stars that explode inside their parent cloud, in setting the cloud dispersal time, and their huge overdensity compared to models where the supernovae are randomly distributed. We also quantify the energy injection efficiency of supernovae as a function of supernova distance to the clouds. We conclude that intermittent driving by supernovae can maintain molecular-cloud turbulence and may be the main process of cloud dispersal. The role of supernovae in the evolution of molecular clouds cannot be fully accounted for without a self-consistent implementation of their feedback., Comment: 33 pages, 23 figures, submitted to ApJ

The Planck Catalogue of Galactic Cold Clumps (PGCC) contains over 13000 sources detected based on their cold dust signature. They are believed to consist of a mixture of quiescent, pre-stellar, and already star-forming objects. We extracted PGCC-type objects from cloud simulations and examined their physical and polarisation properties. The comparison with the PGCC catalogue helps to characterise the PGCC sample and provides valuable tests for numerical simulations of interstellar medium. We used several MHD snapshots to define the density field of our models. Sub-millimetre images of the surface brightness and polarisation were obtained with radiative transfer calculations. We examined the statistics of synthetic cold clump catalogues and examined the variations of the clump polarisation fraction p. The clump sizes, aspect ratios, and temperatures in the synthetic catalogue are similar to the PGCC. The fluxes and column densities are smaller by a factor of a few. Rather than with an increased dust opacity, this could be explained by increasing the average column density of the models by a factor of two to three, close to N(H2)= 10^22 cm-2. When the line of sight is parallel to the mean magnetic field, the polarisation fraction tends to increase towards the clump centres, contrary to observations. When the field is perpendicular, the polarisation fraction tends to decrease towards the clumps, but the drop in $p$ is small (e.g. from p~8% to p~7%). Magnetic field geometry reduces the polarisation fraction in the simulated clumps by only \Delta p~1% on average. The larger drop seen towards the actual PGCC clumps suggests some loss of grain alignment in the dense medium, such as predicted by the radiative torque mechanism. The statistical study is not able to quantify dust opacity changes at the scale of the PGCC clumps., Comment: Accepted to A&A (20 pages, 26 figures)

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

Analysis of all-sky Planck submillimetre observations and the IRAS 100um data has led to the detection of a population of Galactic cold clumps. The clumps can be used to study star formation and dust properties in a wide range of Galactic environments. Our aim is to measure dust spectral energy distribution (SED) variations as a function of the spatial scale and the wavelength. We examine the SEDs at large scales using IRAS, Planck, and Herschel data. At smaller scales, we compare with JCMT/SCUBA-2 850um maps with Herschel data that are filtered using the SCUBA-2 pipeline. Clumps are extracted using the Fellwalker method and their spectra are modelled as modified blackbody functions. According to IRAS and Planck data, most fields have dust colour temperatures T_C ~ 14-18K and opacity spectral index values of beta=1.5-1.9. The clumps/cores identified in SCUBA-2 maps have T~ 13K and similar beta values. There are some indications of the dust emission spectrum becoming flatter at wavelengths longer than 500um. In fits involving Planck data, the significance is limited by the uncertainty of the corrections for CO line contamination. The fits to the SPIRE data give a median beta value slightly above 1.8. In the joint SPIRE and SCUBA-2 850um fits the value decreases to beta ~1.6. Most of the observed T-beta anticorrelation can be explained by noise. The typical submillimetre opacity spectral index beta of cold clumps is found to be ~1.7. This is above the values of diffuse clouds but lower than in some previous studies of dense clumps. There is only tentative evidence of T-beta anticorrelation and beta decreasing at millimetre wavelengths., Comment: 37 pages, accepted to A&A, final version with language corrections

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

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

We discuss the possibility of mapping interstellar clouds at unprecedentedly high spatial resolution by means of near-IR imaging of their scattered light. We calculate the scattering of the interstellar radiation field by a cloud model obtained from the simulation of a supersonic turbulent flow. Synthetic maps of scattered light are computed in the J, H and K bands and are found to allow an accurate estimate of column density, in the range of visual extinction between 1 and 20 magnitudes. We provide a formalism to convert the intensity of scattered light at these near-IR bands into a total gas column density. We also show that this new method of mapping interstellar clouds is within the capability of existing near-IR facilities, which can achieve a spatial resolution of up to ~ 0.1 arcsec. This opens new perspectives in the study of interstellar dust and gas structure on very small scales. The validity of the method has been recently demonstrated by the extraordinary images of the Perseus region obtained by Foster & Goodman (2005), which motivated this investigation., Comment: Submitted to ApJL