Collapse Of Prestellar Cores Research Articles

Mixing is easy: New insights for cosmochemical evolution from pre-stellar core collapse

Signposts of early planet formation are ubiquitous in substructured young discs. Dense, hot, and high-pressure regions that formed during the gravitational collapse process, integral to star formation, facilitate dynamical mixing of dust within the protostellar disc. This provides an incentive to constrain the role of gas and dust interaction and resolve potential zones of dust concentration during star and disc formation stages. We explore whether the thermal and dynamical conditions that developed during protostellar disc formation can generate gas flows that efficiently mix and transport the well-coupled gas and dust components. We simulated the collapse of dusty molecular cloud cores with the hydrodynamics code PLUTO augmented with radiation transport and self-gravity. We used a two-dimensional axisymmetric geometry and followed the azimuthal component of the velocity. The dust was treated as Lagrangian particles that are subject to drag from the gas, whose motion is computed on a Eulerian grid. We considered 1, 10, and 100 micron-sized neutral, spherical dust grains. Importantly, the equation of state accurately includes molecular hydrogen dissociation. We focus on molecular cloud core masses of 1 and 3 $M_ odot $ and explore the effects of different initial rotation rates and cloud core sizes. Our study underlines mechanisms for the early transport of dust from the inner hot disc regions via the occurrence of two transient gas motions, namely meridional flow and outflow. The vortical flow fosters dynamical mixing and retention of dust, while the thermal pressure driven outflow replenishes dust in the outer disc. Notably, these phenomena occur regardless of the initial cloud core mass, size, and rotation rate. Young dynamical precursors to planet-forming discs exhibit regions with complex hydrodynamical gas features and high-temperature structures. These can play a crucial role in concentrating dust for subsequent growth into protoplanets. Dust transport, especially, from sub-au scales surrounding the protostar to the outer relatively cooler parts, offers an efficient pathway for thermal reprocessing during pre-stellar core collapse.

Exploring the dust grain size and polarization mechanism in the hot and massive Class 0 disk IRAS 16293-2422 B

Context. Multiwavelength dust continuum and polarization observations arising from self-scattering have been used to investigate grain sizes in young disks. However, the likelihood of self-scattering being the polarization mechanism in embedded disks decreases for very highly optically thick disks and makes us reconsider some of the size constraints from polarization, particularly for younger and more massive disks. The 1.3 mm polarized emission detected toward the hot (≳400 K) Class 0 disk IRAS 16293-2422 B has been attributed to self-scattering, with predictions of bare grain sizes between 200 and 2000 µm. Aims. We aim to investigate the effects of changing the maximum grain sizes in the resultant continuum and continuum polarization fractions from self-scattering for a hot and massive Class 0 disk extracted from numerical simulations of prestellar core collapse and to compare them with IRAS 16293 B observations. Methods. We compared new and archival dust continuum and polarization observations at high resolution between 1.3 and 18 mm to a set of synthetic models. We developed a new publicly available tool to automate this process called Synthesizer. This tool is an easy-to-use program for generating synthetic observations from numerical simulations. Results. Optical depths are in the range of 130 to 2 from 1.3 to 18 mm, respectively. Predictions of significant grain growth populations, including amax = 1000 µm, are comparable to the observations from IRAS 16293 B at all observed wavelengths. The polarization fraction produced by self-scattering reaches a maximum of approximately 0.1% at 1.3 mm for a maximum grain size of 100 µm, which is an order of magnitude lower than the grain size observed toward IRAS 16293 B. Conclusions. From comparison of the Stokes I fluxes, we conclude that significant grain growth could be present in the young Class 0 disk IRAS 16293 B, particularly in the inner hot region (< 10 au, T > 300 K) where refractory organics evaporate. The polarization produced by self-scattering in our model is not high enough to explain the observations at 1.3 and 7 mm, and such effects as dichroic extinction and polarization reversal of elongated aligned grains remain other possible but untested scenarios.

Massive pre-stellar cores in radiation-magneto-turbulent simulations of molecular clouds

ABSTRACT We simulate the formation and collapse of pre-stellar cores at few-au resolution in a set of radiation-magnetohydrodynamic simulations of giant molecular clouds (GMCs) using the grid-based code RAMSES-RT. We adopt, for the first time to our best knowledge, realistic initial/boundary conditions by zooming in on to individual massive pre-stellar cores within the GMC. We identify two distinct modes of fragmentation: ‘quasi-spherical’ and ‘filamentary’. In both modes, the fragments eventually become embedded in a quasi-steady accretion disc or toroid with radii ∼500–5000 au and opening angles H/R ∼ 0.5 − 1. The discs/toroids are Toomre stable but the accreted pre-existing fragments are found orbiting the outer disc, appearing as disc fragmentation. Each core converts nearly 100 per cent of the gas mass into a few massive stars forming near the disc centre. Large and massive discs around high-mass stars are supported by magnetic pressure in the outer disc, at radii >200–1000 au, and turbulent pressure in the inner disc. The most massive core accretes several times more mass than its initial mass, forming a cluster of 8 massive (proto)stars enshrouded by a toroid, suggesting a competitive accretion scenario for the formation of stars above ∼30 M⊙. We also find that the H ii regions produced by a single massive star remain trapped in the dense circumstellar discs for a few hundred kiloyears, while the dynamic motions of massive stars in wide binaries or multiple systems displace the stars from the densest parts of the disc, allowing UV radiation to escape producing steady or pulsating bipolar H ii regions.

A Massive Star Is Born: How Feedback from Stellar Winds, Radiation Pressure, and Collimated Outflows Limits Accretion onto Massive Stars

Massive protostars attain high luminosities as they are actively accreting and the radiation pressure exerted on the gas in the star’s atmosphere may launch isotropic high-velocity winds. These winds will collide with the surrounding gas producing shock-heated (T ∼ 107 K) tenuous gas that adiabatically expands and pushes on the dense gas that may otherwise be accreted. We present a suite of 3D radiation-magnetohydrodynamic simulations of the collapse of massive prestellar cores and include radiative feedback from the stellar and dust-reprocessed radiation fields, collimated outflows, and, for the first time, isotropic stellar winds to model how these processes affect the formation of massive stars. We find that winds are initially launched when the massive protostar is still accreting and its wind properties evolve as the protostar contracts to the main sequence. Wind feedback drives asymmetric adiabatic wind bubbles that have a bipolar morphology because the dense circumstellar material pinches the expansion of the hot shock-heated gas. We term this the “wind tunnel effect.” If the core is magnetized, wind feedback is less efficient at driving adiabatic wind bubbles initially because magnetic tension delays their growth. We find that wind feedback eventually quenches accretion onto ∼30 M ⊙ protostars that form from the collapse of the isolated cores simulated here. Hence, our results suggest that ≳30 M ⊙ stars likely require larger-scale dynamical inflows from their host cloud to overcome wind feedback. Additionally, we discuss the implications of observing adiabatic wind bubbles with Chandra while the massive protostars are still highly embedded.

The burst mode of accretion in massive star formation with stellar inertia

ABSTRACT The burst mode of accretion in massive star formation is a scenario linking the initial gravitational collapse of parent pre-stellar cores to the properties of their gravitationally unstable discs and of their accretion-driven bursts. In this study, we present a series of high-resolution 3D radiation-hydrodynamics numerical simulations for young massive stars formed out of collapsing $100{\rm M}_{\odot }$ molecular cores, spinning with several values of the ratio of rotational-to-gravitational energies $\beta =5{{-}9\ per cent}$. The models include the indirect gravitational potential caused by disc asymmetries. We find that this modifies the barycentre of the disc, causing significant excursions of the central star position, which we term stellar wobbling. The stellar wobbling slows down and protracts the development of gravitational instability in the disc, reducing the number and magnitude of the accretion-driven bursts undergone by the young massive stars, whose properties are in good agreement with that of the burst monitored from the massive protostar M17 MIR. Including stellar wobbling is therefore important for accurate modelling disc structures. Synthetic alma interferometric images in the millimetre waveband show that the outcomes of efficient gravitational instability such as spiral arms and gaseous clumps can be detected for as long as the disc is old enough and has already entered the burst mode of accretion.

Sulfur gas-phase abundance in dense cores

The abundance of volatile sulfur in dense clouds is long-standing problem in studies of the physics and chemistry of star-forming regions. Sulfur is an important species because its low ionization potential may possibly make it an important charge carrier. The observed sulfur-bearing species in the gas-phase of dense clouds represent only a minor fraction of the cosmic sulfur abundance, which has been interpreted as a signature of sulfur depletion into ices at the surface of dust grains. However, atomic sulfur, which could be the main gas-phase carrier, cannot be observed directly in cold cores. We present measurements of the nitrogen sulfide (NS) radical toward four dense cores performed with the IRAM-30m telescope. Analytical chemical considerations and chemical models over a wide parameter space show that the NS:N2H+ abundance ratio provides a direct constraint on the abundance of gas-phase atomic sulfur. Toward early-type cores, we find that n(S)∕nH is close, or even equal, to the cosmic abundance of sulfur, 14 × 10−6, demonstrating that sulfur is not depleted and is atomic, which is in agreement with chemical models. More chemically evolved cores show sulfur depletion by factors up to 100 in their densest parts. In L1544, atomic sulfur depletion is shown to increase with increasing density. Future observations are needed to discover the solid-phase carrier of sulfur. The initial steps of the collapse of pre-stellar cores in the high sulfur abundance regime also need to be explored from their chemical and dynamical perspectives.

Magnetohydrodynamic effect on first star formation: pre-stellar core collapse and protostar formation

ABSTRACT Recent theoretical studies have suggested that a magnetic field may play a crucial role in the first star formation in the universe. However, the influence of the magnetic field on the first star formation has yet to be understood well. In this study, we perform three-dimensional magnetohydrodynamic simulations taking into account all the relevant cooling processes and non-equilibrium chemical reactions up to the protostar density, in order to study the collapse of magnetized primordial gas cores with self-consistent thermal evolution. Our results show that the thermal evolution of the central core is hardly affected by a magnetic field because magnetic forces do not prevent the contraction along the field lines. We also find that the magnetic braking extracts the angular momentum from the core and suppresses fragmentation depending on the initial strength of the magnetic field. The angular momentum transport by the magnetic outflows is less effective than that by the magnetic braking because the outflows are launched only in a late phase of the collapse. Our results indicate that the magnetic effects become important for the field strength $B\gt 10^{-8}(n_{\rm H}/1\ \rm cm^{-3})^{2/3}\ \rm G$, where nH is the number density, during the collapse phase. Finally, we compare our results with simulations using a barotropic approximation and confirm that this approximation is reasonable at least for the collapse phase. Nevertheless, self-consistent treatment of the thermal and chemical processes is essential for extending simulations to the accretion phase, in which radiative feedback by protostars plays a crucial role.

Protoplanetary disk formation from the collapse of a prestellar core

Context.Between the two research communities that study star formation and protoplanetary disk evolution, only a few efforts have been made to understand and bridge the gap between studies of a collapsing prestellar core and a developed disk. While it has generally been accepted for about a decade that the magnetic field and its nonideal effects play important roles during the stellar formation, simple models of pure hydrodynamics and angular momentum conservation are still widely employed in the studies of disk assemblage in the framework of the so-called alpha-disk model because these models are simple.Aims.We revisit the assemblage phase of the protoplanetary disk and employ current knowledge of the prestellar core collapse.Methods.We performed 3D magnetohydrodynamic (MHD) simulations with ambipolar diffusion and full radiative transfer to follow the formation of the protoplanetary disk within a collapsing prestellar core. The global evolution of the disk and its internal properties were analyzed to understand how the infalling envelope regulates the buildup and evolution of the disk. We followed the global evolution of the protoplanetary disk from the prestellar core collapse during 100 kyr with a reasonable resolution of AU. Two snapshots from this reference run were extracted and rerun with significantly increased resolution to resolve the interior of the disk.Results.The disk that formed under our simulation setup is more realistic and agrees with recent observations of disks around class 0 young stellar objects. The source function of the mass flux that arrives at the disk and the radial mass accretion rate within the disk are measured and compared to analytical self-similar models based on angular momentum conservation. The source function is very centrally peaked compared to classical hydrodynamical models, implying that most of the mass falling onto the star does not transit through the midplane of the disk. We also found that the disk midplane is almost dead to turbulence, whereas upper layers and the disk outer edge are highly turbulent, and this is where the accretion occurs. The snow line, located at about 5–10 AU during the infall phase, is significantly farther away from the center than in a passive disk. This result might be of numerical origin.Conclusions.We studied self-consistent protoplanetary disk formation from prestellar core collapse, taking nonideal MHD effects into account. We developed a zoomed rerun technique to quickly obtain a reasonable disk that is highly stratified, weakly magnetized inside, and strongly magnetized outside. During the class 0 phase of protoplanetary disk formation, the interaction between the disk and the infalling envelope is important and ought not be neglected. We measured the complex flow pattern and compared it to the classical models of pure hydrodynamical infall. Accretion onto the star is found to mostly depend on dynamics at large scales, that is, the collapsing envelope, and not on the details of the disk structure.

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

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

Massive-star Formation via the Collapse of Subvirial and Virialized Turbulent Massive Cores

Similar to their low-mass counterparts, massive stars likely form via the collapse of pre-stellar molecular cores. Recent observations suggest that most massive cores are subvirial (i.e., not supported by turbulence) and therefore are likely unstable to gravitational collapse. Here we perform radiation hydrodynamic simulations to follow the collapse of turbulent massive pre-stellar cores with subvirial and virialized initial conditions to explore how their dynamic state affects the formation of massive stars and core fragmentation into companion stars. We find that subvirial cores undergo rapid monolithic collapse resulting in higher accretion rates at early times as compared to the collapse of virialized cores that have the same physical properties. In contrast, we find that virialized cores undergo a slower, gradual collapse and significant turbulent fragmentation at early times resulting in numerous companion stars. In the absence of strong magnetic fields and protostellar outflows we find that the faster growth rate of massive stars that are born out of subvirial cores leads to an increase in the radiative heating of the core thereby further suppressing fragmentation at early times when turbulent fragmentation occurs for virialized cores. Regardless of initial condition, we find that the massive accretion disks that form around massive stars dominant the accretion flow onto the star at late times and eventually become gravitationally unstable and fragment to form companion stars at late times.

Gravitoviscous protoplanetary disks with a dust component

Aims. The central region of a circumstellar disk is difficult to resolve in global numerical simulations of collapsing cloud cores, but its effect on the evolution of the entire disk can be significant. Methods. We used numerical hydrodynamics simulations to model the long-term evolution of self-gravitating and viscous circumstellar disks in the thin-disk limit. Simulations start from the gravitational collapse of pre-stellar cores of 0.5–1.0 M⊙ and both gaseous and dusty subsystems were considered, including a model for dust growth. The inner unresolved 1.0 au of the disk is replaced with a central smart cell (CSC), a simplified model that simulates physical processes that may occur in this region. Results. We found that the mass transport rate through the CSC has an appreciable effect on the evolution of the entire disk. Models with slow mass transport form more massive and warmer disks, and are more susceptible to gravitational instability and fragmentation, including a newly identified episodic mode of disk fragmentation in the T Tauri phase of disk evolution. Models with slow mass transport through the CSC feature episodic accretion and luminosity bursts in the early evolution, while models with fast transport are characterized by a steadily declining accretion rate with low-amplitude flickering. Dust grows to a larger, decimeter size in the slow transport models and efficiently drifts in the CSC, where it accumulates and reaches the limit where a streaming instability becomes operational. We argue that gravitational instability, together with a streaming instability likely operating in the inner disk regions, constitute two concurrent planet-forming mechanisms, which may explain the observed diversity of exoplanetary orbits. Conclusions. We conclude that sophisticated models of the inner unresolved disk regions should be used when modeling the formation and evolution of gaseous and dusty protoplanetary disks.

On the ALMA observability of nascent massive multiple systems formed by gravitational instability

Massive young stellar object (MYSOs) form during the collapse of high-mass pre-stellar cores, where infalling molecular material is accreted through a centrifugally-balanced accretion disc that is subject to efficient gravitational instabilities. In the resulting fragmented accretion disc of the MYSO, gaseous clumps and low-mass stellar companions can form, which will influence the future evolution of massive protostars in the Hertzsprung-Russell diagram. We perform dust continuum radiative transfer calculations and compute synthetic images of disc structures modelled by the gravito-radiation-hydrodynamics simulation of a forming MYSO, in order to investigate the Atacama Large Millimeter/submillimeter Array (ALMA) observability of circumstellar gaseous clumps and forming multiple systems. Both spiral arms and gaseous clumps located at ~a few 100 au from the protostar can be resolved by interferometric ALMA Cycle 7 C43-8 and C43-10 observations at band 6 (1.2 mm), using a maximal 0.015" beam angular resolution and at least 10-30 min exposure time for sources at distances of 1-2 kpc. Our study shows that substructures are observable regardless of their viewing geometry or can be inferred in the case of an edge-viewed disc. The observation probability of the clumps increases with the gradually increasing efficiency of gravitational instability at work as the disc evolves. As a consequence, large discs around MYSOs close to the zero-age-main-sequence line exhibit more substructures than at the end of the gravitational collapse. Our results motivate further observational campaigns devoted to the close surroundings of the massive protostars S255IR-NIRS3 and NGC 6334I-MM1, whose recent outbursts are a probable signature of disc fragmentation and accretion variability.

The ortho-to-para ratio of water in interstellar clouds

The nuclear-spin chemistry of interstellar water is investigated using the University of Grenoble Alpes Astrochemical Network (UGAN). This network includes reactions involving the different nuclear-spin states of the hydrides of carbon, nitrogen, oxygen and sulphur, as well as their deuterated forms. Nuclear-spin selection rules are implemented within the scrambling hypothesis for reactions involving up to seven protons. The abundances and ortho-to-para ratios (OPRs) of gas-phase water and water ions (H$_2$O$^+$ and H$_3$O$^+$) are computed under the steady-state conditions representative of a dark molecular cloud and during the early phase of gravitational collapse of a prestellar core. The model incorporates the freezing of the molecules on to grains, simple grain surface chemistry and cosmic-ray induced and direct desorption of ices. The predicted OPRs are found to deviate significantly from both thermal and statistical values and to be independent of temperature below $\sim $30~K. The OPR of H$_2$O is shown to lie between 1.5 and 2.6, depending on the spin-state of H$_2$, in good agreement with values derived in translucent clouds with relatively high extinction. In the prestellar core collapse calculations, the OPR of H$_2$O is shown to reach the statistical value of 3 in regions with severe depletion ($n_{\rm H}> 10^7$~cm$^{-3}$). We conclude that a low water OPR ($\lesssim 2.5$) is consistent with gas-phase ion-neutral chemistry and reflects a gas with OPR(H$_2)\lesssim 1$. Available OPR measurements in protoplanetary disks and comets are finally discussed.

New constraints on the initial parameters of low-mass star formation from chemical modeling

The complexity of physico-chemical models of star formation is increasing, with models that take into account new processes and more realistic setups. These models allow astrochemists to compute the evolution of chemical species throughout star formation. Hence, comparing the outputs of such models to observations allows to bring new constraints on star formation. The work presented in this paper is based on the recent public release of a database of radiation hydrodynamical low-mass star formation models. We used this database as physical parameters to compute the time dependent chemical composition of collapsing cores with a 3-phase gas-grain model. The results are analyzed to find chemical tracers of the initial physical parameters of collapse such as the mass, radius, temperature, density, and free-fall time. They are also compared to observed molecular abundances of Class 0 protostars. We find numerous tracers of the initial parameters of collapse, except for the initial mass. More particularly, we find that gas phase CH3CN, NS and OCS trace the initial temperature, while H2CS trace the initial density and free-fall time of the parent cloud. The comparison of our results with a sample of 12 Class 0 low mass protostars allows us to constrain the initial parameters of collapse of low-mass prestellar cores. We find that low-mass protostars are preferentially formed within large cores with radii greater than 20000 au, masses between 2 and 4 Msol, temperatures lower or equal to 15 K, and densities between 6e4 and 2.5e5 part.cm-3, corresponding to free-fall times between 100 and 200 kyrs.

Jets and outflows of massive protostars

Context. Massive stars live short but intense lives. While less numerous than low-mass stars, they enormously impact their surroundings by several feedback mechanisms. They form in opaque and far-away regions of the galaxy, such that one of these feedback mechanisms also becomes a record of their evolution: their bright large-scale jets and outflows. Aims. In a comprehensive convergence study, we investigate the computational conditions necessary to resolve (pseudo-) disk formation and jet-launching processes, and analyze possible caveats. We explore the magneto-hydrodynamic (MHD) processes of the collapse of massive prestellar cores in detail, including an analysis of the forces involved and their temporal evolution for up to two free-fall times. Methods. We conduct MHD simulations using the state-of-the-art code PLUTO, combining nonideal MHD, self-gravity, and very high resolutions as they have never been achieved before. Our setup includes a 100 M⊙ cloud core that collapses under its own self-gravity to self-consistently form a dense disk structure and launch tightly collimated magneto-centrifugal jets and wide-angle tower flows. Results. We show a comprehensive evolutionary picture of the collapse of a massive prestellar core with a detailed analysis of the physical processes involved and our high-resolution simulations can resolve a magneto-centrifugal jet and a magnetic pressure-driven outflow, separately. The nature of the outflows depends critically on spatial resolution. Only high-resolution simulations are able to differentiate a magneto-centrifugally launched, highly collimated jet from a slow wide-angle magnetic-pressure-driven tower flow. Of these two outflow components, the tower flow dominates angular-momentum transport. The mass outflow rate is dominated by the entrained material from the interaction of the jet with the stellar environment and only part of the ejected medium is directly launched from the accretion disk. A tower flow can only develop to its full extent when much of the original envelope has already dispersed. Taking into account both the mass launched from the surface of the disk and the entrained material from the envelope, we find an ejection-to-accretion efficiency of 10%. Nonideal MHD is required to form centrifugally supported accretion disks and the disk size is strongly dependent on spatial resolution. A converged result for disk and both outflow components requires a spatial resolution of Δx ≤ 0.17 au at 1 au and sink-cell sizes ≤3.1 au. Conclusions. Massive stars not only possess slow wide-angle tower flows, but also produce magneto-centrifugal jets, just as their low-mass counterparts. The actual difference between low-mass and high-mass star formation lies in the “embeddedness” of the high-mass star which implies that the jet and tower flow interact with the infalling large-scale stellar environment, potentially resulting in entrainment.

ALMA observations of the young protostellar system Barnard 1b: signatures of an incipient hot corino in B1b-S.

The Barnard 1b core shows signatures of being at the earliest stages of low-mass star formation, with two extremely young and deeply embedded protostellar objects. Hence, this core is an ideal target to study the structure and chemistry of the first objects formed in the collapse of prestellar cores. We present ALMA Band 6 spectral line observations at ~0.6″ of angular resolution towards Barnard 1b. We have extracted the spectra towards both protostars, and used a Local Thermodynamic Equilibrium (LTE) model to reproduce the observed line profiles. B1b-S shows rich and complex spectra, with emission from high energy transitions of complex molecules, such as CH3OCOH and CH3CHO, including vibrational level transitions. We have tentatively detected for the first time in this source emission from NH2CN, NH2CHO, CH3CH2OH, CH2OHCHO, CH3CH2OCOH and both aGg' and gGg' conformers of (CH2OH)2. This is the first detection of ethyl formate (CH3CH2OCOH) towards a low-mass star forming region. On the other hand, the spectra of the FHSC candidate B1b-N are free of COMs emission. In order to fit the observed line profiles in B1b-S, we used a source model with two components: an inner hot and compact component (200 K, 0.35″) and an outer and colder one (60 K, 0.6″). The resulting COM abundances in B1b-S range from 10-13 for NH2CN and NH2CHO, up to 10-9 for CH3OCOH. Our ALMA Band 6 observations reveal the presence of a compact and hot component in B1b-S, with moderate abundances of complex organics. These results indicate that a hot corino is being formed in this very young Class 0 source.

Effect of accretion on the pre-main-sequence evolution of low-mass stars and brown dwarfs

The pre-main-sequence evolution of low-mass stars and brown dwarfs is studied numerically starting from the formation of a protostellar/proto-brown dwarf seed and taking into account the mass accretion onto the central object during the initial several Myr of evolution. The stellar evolution was computed using the STELLAR evolution code developed by Yorke & Bodenheimer with recent modifications by Hosokawa et al. The mass accretion rates were taken from numerical hydrodynamics models of Vorobyov & Basu computing the circumstellar disk evolution starting from the gravitational collapse of pre-stellar cloud cores of various mass and angular momentum. The resulting stellar evolution tracks were compared with the isochrones and isomasses calculated using non-accreting models. We find that mass accretion in the initial several Myr of protostellar evolution can have a strong effect on the subsequent evolution of young stars and brown dwarfs. The disagreement between accreting and non-accreting models in terms of the total stellar luminosity L_st, stellar radius R_st and effective temperature T_eff depends on the thermal efficiency of accretion, i.e., on the fraction of accretion energy absorbed by the central object. The largest mismatch is found for the cold accretion case, in which essentially all accretion energy is radiated away. The relative deviations in L_st and R_st in this case can reach 50% for 1.0-Myr-old objects and remain notable even for 10-Myr-old objects. In the hot and hybrid accretion cases, in which a constant fraction of accretion energy is absorbed, the disagreement between accreting and non-accreting models becomes less pronounced, but still remains notable for 1.0-Myr-old objects. These disagreements may lead to the wrong age estimate for objects of (sub-)solar mass when using the isochrones based on non-accreting models (abridged).

An ALMA and MagAO Study of the Substellar Companion GQ Lup B∗

Abstract Multi-wavelength observations provide a complementary view of the formation of young, directly imaged planet-mass companions. We report the ALMA 1.3 mm and Magellan adaptive optics Hα, , , and Y S observations of the GQ Lup system, a classical T Tauri star with a substellar companion at ∼110 au projected separation. We estimate the accretion rates for both components from the observed Hα fluxes. In our ∼0.″05 resolution ALMA map, we resolve GQ Lup A’s disk in the dust continuum, but no signal is found from the companion. The disk is compact, with a radius of ∼22 au, a dust mass of ∼6 M ⊕, an inclination angle of ∼56°, and a very flat surface density profile indicative of a radial variation in dust grain sizes. No gaps or inner cavity are found in the disk, so there is unlikely a massive inner companion to scatter GQ Lup B outward. Thus, GQ Lup B might have formed in situ via disk fragmentation or prestellar core collapse. We also show that GQ Lup A’s disk is misaligned with its spin axis, and possibly with GQ Lup B’s orbit. Our analysis on the tidal truncation radius of GQ Lup A’s disk suggests that GQ Lup B’s orbit might have a low eccentricity.

An unstable truth: how massive stars get their mass

The pressure exerted by massive stars' radiation fields is an important mechanism regulating their formation. Detailed simulation of massive star formation therefore requires an accurate treatment of radiation. However, all published simulations have either used a diffusion approximation of limited validity; have only been able to simulate a single star fixed in space, thereby suppressing potentially-important instabilities; or did not provide adequate resolution at locations where instabilities may develop. To remedy this we have developed a new, highly accurate radiation algorithm that properly treats the absorption of the direct radiation field from stars and the re-emission and processing by interstellar dust. We use our new tool to perform three-dimensional radiation-hydrodynamic simulations of the collapse of massive pre-stellar cores with laminar and turbulent initial conditions and properly resolve regions where we expect instabilities to grow. We find that mass is channeled to the stellar system via gravitational and Rayleigh-Taylor (RT) instabilities, in agreement with previous results using stars capable of moving, but in disagreement with methods where the star is held fixed or with simulations that do not adequately resolve the development of RT instabilities. For laminar initial conditions, proper treatment of the direct radiation field produces later onset of instability, but does not suppress it entirely provided the edges of radiation-dominated bubbles are adequately resolved. Instabilities arise immediately for turbulent pre-stellar cores because the initial turbulence seeds the instabilities. Our results suggest that RT features are significant and should be present around accreting massive stars throughout their formation.

Chemical solver to compute molecule and grain abundances and non-ideal MHD resistivities in prestellar core-collapse calculations

We develop a detailed chemical network relevant to the conditions characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of Potassium, Sodium and Hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to n_H = 10^12 cm^-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process.