Turbulent Cores Research Articles

SHAPING THE BROWN DWARF DESERT: PREDICTING THE PRIMORDIAL BROWN DWARF BINARY DISTRIBUTIONS FROM TURBULENT FRAGMENTATION

The formation of brown dwarfs (BDs) poses a key challenge to star formation theory. The observed dearth of nearby ($\leq 5$ AU) brown dwarf companions to solar-mass stars, known as the brown dwarf desert, as well as the tendency for low-mass binary systems to be more tightly-bound than stellar binaries, have been cited as evidence for distinct formation mechanisms for brown dwarfs and stars. In this paper, we explore the implications of the minimal hypothesis that brown dwarfs in binary systems originate via the same fundamental fragmentation mechanism as stars, within isolated, turbulent giant molecular cloud cores. We demonstrate analytically that the scaling of specific angular momentum with turbulent core mass naturally gives rise to the brown dwarf desert, as well as wide brown-dwarf binary systems. Further, we show that the turbulent core fragmentation model also naturally predicts that very low-mass (VLM) binary and BD/BD systems are more tightly-bound than stellar systems. In addition, in order to capture the stochastic variation intrinsic to turbulence, we generate $10^4$ model turbulent cores with synthetic turbulent velocity fields to show that the turbulent fragmentation model accommodates a small fraction of binary brown dwarfs with wide separations, similar to observations. Indeed, the picture which emerges from the turbulent fragmentation model is that a single fragmentation mechanism may largely shape both stellar and brown dwarf binary distributions during formation.

Disc formation in turbulent cloud cores: is magnetic flux loss necessary to stop the magnetic braking catastrophe or not?

Recent numerical analysis of Keplerian disc formation in turbulent, magnetized cloud cores by Santos-Lima et al. demonstrated that reconnection diffusion is an efficient process to remove the magnetic flux excess during the buildup of a rotationally supported disc. This process is induced by fast reconnection of the magnetic fields in a turbulent flow. In a similar numerical study, Seifried et al. concluded that reconnection diffusion or any other non-ideal magnetohydrodynamic effects would not be necessary and turbulence shear alone would provide a natural way to build up a rotating disc without requiring magnetic flux loss. Their conclusion was based on the fact that the mean mass-to-flux ratio (μ) evaluated over a spherical region with a radius much larger than the disc is nearly constant in their models. In this paper, we compare the two sets of simulations and show that this averaging over large scales can mask significant real increases of μ in the inner regions where the disc is built up. We demonstrate that turbulence-induced reconnection diffusion of the magnetic field happens in the initial stages of the disc formation in the turbulent envelope material that is accreting. Our analysis is suggestive that reconnection diffusion is present in both sets of simulations and provides a simple solution for the ‘magnetic braking catastrophe’ which is discussed in the literature in relation to the formation of protostellar accretion discs.

Formation and early evolution of circumstellar discs in turbulent molecular cloud cores

We investigate the formation and evolution of circumstellar disks in turbulent cloud cores until several 104 years after protostar formation using smoothed particle hydrodynamics (SPH) calculations. The formation and evolution process of circumstellar disk in turbulent cloud cores differs substantially from that in rigidly rotating cloud cores. In turbulent cloud cores, a filamentary structure appears before the protostar formation and the protostar forms in the filament. If the turbulence is initially sufficiently strong, the remaining filament twists around the protostar and directly becomes a rotation-supported disk. Upon formation, the disk orientation is generally misaligned with the angular momentum of its host cloud core and it dynamically varies during the main accretion phase, even though the turbulence is weak. This is because the angular momentum of the entire cloud core is mainly determined by the large scale velocity field whose wavelength is comparable to the cloud scale, whereas the angular momentum of the disk is determined by the local velocity field where the protostar forms and these two velocity fields do not correlate with each other. In the case of disk evolution in a binary or multiple stars, the disks are misaligned with each other at least during the main accretion phase, because there is no correlation between the velocity fields around the position where each protostar forms. In addition, each disk is also misaligned with the binary orbital plane. Such misalignment can explain the recent observations of misaligned disks and misaligned protostellar outflows.

Disc formation in turbulent massive cores: circumventing the magnetic braking catastrophe

ABSTRACT We present collapse simulations of 100 M⊙ turbulent cloud cores threaded by a strong magnetic field. During the initial collapse phase, filaments are generated which fragment quickly and form several protostars. Around these protostars Keplerian discs with typical sizes of up to 100 au build up in contrast to previous simulations neglecting turbulence. We examine three mechanisms potentially responsible for lowering the magnetic braking efficiency and therefore allowing for the formation of Keplerian discs. Analysing the condensations in which the discs form, we show that the build-up of Keplerian discs is neither caused by magnetic flux loss due to turbulent reconnection nor by the misalignment of the magnetic field and the angular momentum. It is rather a consequence of the turbulent surroundings of the disc which exhibit no coherent rotation structure while strong local shear flows carry large amounts of angular momentum. We suggest that the ‘magnetic braking catastrophe’, i.e. the formation of sub-Keplerian discs only, is an artefact of the idealized non-turbulent initial conditions and that turbulence provides a natural mechanism to circumvent this problem.

Statistical analysis of the mass-to-flux ratio in turbulent cores: effects of magnetic field reversals and dynamo amplification

We study the mass-to-flux ratio (M/\Phi) of clumps and cores in simulations of supersonic, magnetohydrodynamical turbulence for different initial magnetic field strengths. We investigate whether the (M/\Phi)-ratio of core and envelope, R = (M/\Phi)_{core}/(M/\Phi)_{envelope} can be used to distinguish between theories of ambipolar diffusion and turbulence-regulated star formation. We analyse R for different Lines-of-Sight (LoS) in various sub-cubes of our simulation box. We find that, 1) the average and median values of |R| for different times and initial magnetic field strengths are typically greater, but close to unity, 2) the average and median values of |R| saturate at average values of |R| ~ 1 for smaller magnetic fields, 3) values of |R| < 1 for small magnetic fields in the envelope are caused by field reversals when turbulence twists the field lines such that field components in different directions average out. Finally, we propose two mechanisms for generating values |R| ~< 1 for the weak and strong magnetic field limit in the context of a turbulent model. First, in the weak field limit, the small-scale turbulent dynamo leads to a significantly increased flux in the core and we find |R| ~< 1. Second, in the strong field limit, field reversals in the envelope also lead to values |R| ~< 1. These reversals are less likely to occur in the core region where the velocity field is more coherent and the internal velocity dispersion is typically subsonic.

Importance of the initial conditions for star formation - II. Fragmentation-induced starvation and accretion shielding

We investigate the impact of different initial conditions for the initial density profile and the initial turbulence on the formation process of protostellar clusters. We study the collapse of dense molecula r cloud cores with three-dimensional adaptive mesh refinement simulations. We focus our discussion on the distribution of the gas among the protostellar objects in the turbulent dynamical cores. Despite the large variations in the initial configurations and the resulting gas and cluster morphology we find that all stellar clusters follow a very similar gas accretion behaviour. Once secondary protostars begin to form, the central region of a cluster is efficiently shielded from further accretion. Hence, objects located close to the centre are starved of material, as indicated by a strong decrease of the central acc retion rate. This Fragmentation Induced Starvation occurs not only in rotationally supported discs and filaments, but also in more spherically symmetric clusters with complex chaotic motions.

The influence of the turbulent perturbation scale on pre-stellar core fragmentation and disc formation

The collapse of weakly turbulent pre-stellar cores is a critical stage in the process of star formation. Being highly non-linear and stochastic, the outcome of collapse can only be explored theoretically by performing large ensembles of numerical simulations. Standard practice is to quantify the initial turbulent velocity field in a core in terms of the amount of turbulent energy (or some equivalent) and the exponent in the power spectrum (n ≡−dlog Pk/dlog k). In this paper, we present a numerical study of the influence of the details of the turbulent velocity field on the collapse of an isolated, weakly turbulent, low-mass pre-stellar core. We show that, as long as n ≳ 3 (as is usually assumed), a more critical parameter than n is the maximum wavelength in the turbulent velocity field, λMAX. This is because λMAX carries most of the turbulent energy, and thereby influences both the amount and the spatial coherence of the angular momentum in the core. We show that the formation of dense filaments during collapse depends critically on λMAX, and we explain this finding using a force balance analysis. We also show that the core has only a high probability of fragmenting if λMAX > RCORE/2 (where RCORE is the core radius), the dominant mode of fragmentation involves the formation and break-up of filaments and although small protostellar discs (with radius RDISC≲ 20 au) form routinely, more extended discs are rare. In turbulent, low-mass cores of the type we simulate here, the formation of large, fragmenting protostellar discs is suppressed by early fragmentation in the filaments.

Importance of the initial conditions for star formation - I. Cloud evolution and morphology

We present a detailed parameter study of collapsing turbulent cloud cores, varying the initial density profile and the initial turbulent velocity field. We systematically investigate the influence of different initial conditions on the star formation process, mainly focusing on the fragmentation, the number of formed stars, and the resulting mass distributions. Our study compares four different density profiles (uniform, Bonnor-Ebert type, $\rho\propto r^{-1.5}$, and $\rho\propto r^{-2}$), combined with six different supersonic turbulent velocity fields (compressive, mixed, and solenoidal, initialised with two different random seeds each) in three-dimensional simulations using the adaptive-mesh refinement, hydrodynamics code FLASH. The simulations show that density profiles with flat cores produce hundreds of low-mass stars, either distributed throughout the entire cloud or found in subclusters, depending on the initial turbulence. Concentrated density profiles always lead to the formation of one high-mass star in the centre of the cloud and, if at all, low-mass stars surrounding the central one. In uniform and Bonnor-Ebert type density distributions, compressive initial turbulence leads to local collapse about 25% earlier than solenoidal turbulence. However, central collapse in the steep power-law profiles is too fast for the turbulence to have any significant influence. We conclude that (I) the initial density profile and turbulence mainly determine the cloud evolution and the formation of clusters, (II) the initial mass function (IMF) is not universal for all setups, and (III) that massive stars are much less likely to form in flat density distributions. The IMFs obtained in the uniform and Bonnor-Ebert type density profiles are more consistent with the observed IMF, but shifted to lower masses.

Collapse, outflows and fragmentation of massive, turbulent and magnetized prestellar barotropic cores

Stars and more particularly massive stars, have a drastic impact on galaxy evolution. Yet the conditions in which they form and collapse are still not fully understood. In particular, the influence of the magnetic field on the collapse of massive clumps is relatively unexplored, it is thus of great relevance in the context of the formation of massive stars to investigate its impact. We perform high resolution, MHD simulations of the collapse of hundred solar masses, turbulent and magnetized clouds, using the adaptive mesh refinement code RAMSES. We compute various quantities such as mass distribution, magnetic field and angular momentum within the collapsing core and study the episodic outflows and the fragmentation that occurs during the collapse. The magnetic field has a drastic impact on the cloud evolution. We find that magnetic braking is able to substantially reduce the angular momentum in the inner part of the collapsing cloud. Fast and episodic outflows are being launched with typical velocities of the order of 3-5 km s$^{-1}$ although the highest velocities can be as high as 30-40 km s$^{-1}$. The fragmentation in several objects, is reduced in substantially magnetized clouds with respect to hydrodynamical ones by a factor of the order of 1.5-2. We conclude that magnetic fields have a significant impact on the evolution of massive clumps. In combination with radiation, magnetic fields largely determine the outcome of massive core collapse. We stress that numerical convergence of MHD collapse is a challenging issue. In particular, numerical diffusion appears to be important at high density therefore possibly leading to an over-estimation of the number of fragments.

PROTOSTELLAR COLLAPSE OF MAGNETO-TURBULENT CLOUD CORES: SHAPE DURING COLLAPSE AND OUTFLOW FORMATION

We investigate protostellar collapse of molecular cloud cores by numerical simulations, taking into account turbulence and magnetic fields. By using the adaptive mesh refinement technique, the collapse is followed over a wide dynamic range from the scale of a turbulent cloud core to that of the first core. The cloud core is lumpy in the low density region owing to the turbulence, while it has a smooth density distribution in the dense region produced by the collapse. The shape of the dense region depends mainly on the mass of the cloud core; a massive cloud core tends to be prolate while a less massive cloud core tends to be oblate. In both cases, anisotropy of the dense region increases during the isothermal collapse. The minor axis of the dense region is always oriented parallel to the local magnetic field. All the models eventually yield spherical first cores supported mainly by the thermal pressure. Most of turbulent cloud cores exhibit protostellar outflows around the first cores. These outflows are classified into two types, bipolar and spiral flows, according to the morphology of the associated magnetic field. Bipolar flow often appears in the less massive cloud core. The rotation axis of the first core is oriented parallel to the local magnetic field for bipolar flow, while the orientation of the rotation axis from the global magnetic field depends on the magnetic field strength. In spiral flow, the rotation axis is not aligned with the local magnetic field.

Searching for Galactic hidden gas through interstellar scintillation: results from a test with the NTT-SOFI detector

Aims: Stars twinkle because their light propagates through the atmosphere. The same phenomenon is expected at a longer time scale when the light of remote stars crosses an interstellar molecular cloud, but it has never been observed at optical wavelength. In a favorable case, the light of a background star can be subject to stochastic fluctuations on the order of a few percent at a characteristic time scale of a few minutes. Our ultimate aim is to discover or exclude these scintillation effects to estimate the contribution of molecular hydrogen to the Galactic baryonic hidden mass. This feasibility study is a pathfinder toward an observational strategy to search for scintillation, probing the sensitivity of future surveys and estimating the background level. Methods: We searched for scintillation induced by molecular gas in visible dark nebulae as well as by hypothetical halo clumpuscules of cool molecular hydrogen ($\mathrm{H_2-He}$) during two nights. We took long series of 10s infrared exposures with the ESO-NTT telescope toward stellar populations located behind visible nebulae and toward the Small Magellanic Cloud (SMC). We therefore searched for stars exhibiting stochastic flux variations similar to what is expected from the scintillation effect. According to our simulations of the scintillation process, this search should allow one to detect (stochastic) transverse gradients of column density in cool Galactic molecular clouds of order of $\sim 3\times 10^{-5}\,\mathrm{g/cm^2/10\,000\,km}$. Results: We found one light-curve that is compatible with a strong scintillation effect through a turbulent structure characterized by a diffusion radius $R_{diff}<100\, km$ in the B68 nebula. Complementary observations are needed to clarify the status of this candidate, and no firm conclusion can be established from this single observation. We can also infer limits on the existence of turbulent dense cores (of number density $n>10^9\, cm^{-3}$) within the dark nebulae. Because no candidate is found toward the SMC, we are also able to establish upper limits on the contribution of gas clumpuscules to the Galactic halo mass. Conclusions: The limits set by this test do not seriously constrain the known models, but we show that the short time-scale monitoring for a few $10^6 star\times hour$ in the visible band with a $>4$ meter telescope and a fast readout camera should allow one to quantify the contribution of turbulent molecular gas to the Galactic halo. The LSST (Large Synoptic Survey Telescope) is perfectly suited for this search.

Implementing and comparing sink particles in AMR and SPH

AbstractWe implemented sink particles in the Adaptive Mesh Refinement (AMR) code FLASH to model the gravitational collapse and accretion in turbulent molecular clouds and cores. Sink particles are frequently used to measure properties of star formation in numerical simulations, such as the star formation rate and efficiency, and the mass distribution of stars. We show that only using a density threshold for sink particle creation is insufficient in case of supersonic flows, because the density can exceed the threshold in strong shocks that do not necessarily lead to local collapse. Additional physical collapse indicators have to be considered. We apply our AMR sink particle module to the formation of a star cluster, and compare it to a Smoothed Particle Hydrodynamics (SPH) code with sink particles. Our comparison shows encouraging agreement of gas and sink particle properties between the AMR and SPH code.

The Luminosity Problem: Testing Theories of Star Formation

AbstractLow-mass protostars are less luminous than expected. This luminosity problem is important because the observations appear to be inconsistent with some of the basic premises of star formation theory. Two possible solutions are that stars form slowly, which is supported by recent data, and/or that protostellar accretion is episodic; current data suggest that the latter accounts for less than half the missing luminosity. The solution to the luminosity problem bears directly on the fundamental problem of the time required to form a low-mass star. The protostellar mass and luminosity functions provide powerful tools both for addressing the luminosity problem and for testing theories of star formation. Results are presented for the collapse of singular isothermal spheres, for the collapse of turbulent cores, and for competitive accretion.

Protostellar discs formed from turbulent cores

We investigate the collapse and fragmentation of low-mass, trans-sonically turbulent prestellar cores, using SPH simulations. The initial conditions are slightly supercritical Bonnor-Ebert spheres, all with the same density profile, the same mass (M_O=6.1 Msun) and the same radius (R_O=17,000 AU), but having different initial turbulent velocity fields. Four hundred turbulent velocity fields have been generated, all scaled so that the mean Mach number is M=1. Then a subset of these, having a range of net angular momenta, j, has been evolved. The evolution of these turbulent cores is not strongly correlated with j. Instead it is moderated by the formation of filamentary structures due to converging turbulent flows. A high fraction (~ 82%) of the protostars forming from turbulent cores are attended by protostellar accretion discs, but only a very small fraction (~16%) of these discs is sufficiently cool and extended to develop non-linear gravitational instabilities and fragment.

PANCHROMATIC SPECTRAL ENERGY DISTRIBUTIONS OF DUSTY GALAXIES WITH RADISHE. I. PREDICTIONS FORHERSCHEL: CORRELATING COLORS WITH GALACTIC ENERGY SOURCES

We present three-dimensional, self-consistent radiative transfer solutions with a new Monte Carlo radiative equilibrium code. The code RADISHE can be applied to calculate the emergent spectral energy distributions (SEDs) and broadband images from optical to millimeter wavelengths of arbitrary density geometries with distributed sources of radiation. One of the primary uses of this code has been to interface with hydrodynamical codes to calculate emergent SEDs along a simulation time sequence. We focus the applications of this code in this paper on infrared bright dusty galaxies, but RADISHE is also ideal for calculating emission from star clusters or protostellar environments. The primary methodological focus of this paper is on the radiative equilibrium temperature calculation. We find that an iterative calculation of the temperature, which takes the sum of photon flight paths as the Monte Carlo estimator for the mean free intensity, is significantly faster than relaxation temperature calculation methods, particularly when large numbers of grid cells are required (i.e., in modeling three-dimensional geometries such as the dust envelopes of turbulent massive protostellar cores or infrared-bright galaxies). An accurate long-wavelength SED and corresponding temperature calculation will be essential for analyzing upcoming Herschel Space Observatory observations. We present simulated color-color plots for infrared-bright galaxies at a range of redshifts, and unfold these plots as color versus the fractional active galactic nucleus (AGN) luminosity, to demonstrate that Herschel will be able to effectively discriminate between submillimeter galaxies, where the energy source is dominated by AGN, and those where star formation dominates. We demarcate in particular the Class II or energetically active AGN evolutionary phase in Herschel color-color plots.

Dense Molecular Gas in a Young Cluster around MWC 1080: Rule of the Massive Star

We present CS -->J = 2→ 1,13CO -->J = 1→ 0, and C18O -->J = 1→ 0, observations with the 10 element Berkeley Illinois Maryland Association (BIMA) array toward the young cluster around the Be star MWC 1080. These observations reveal a biconical outflow cavity with size ~0.3 and 0.05 pc for the semimajor and semiminor axis and ~45° position angle. These transitions trace the dense gas, which is likely the swept-up gas of the outflow cavity, rather than the remaining natal gas or the outflow gas. The gas is clumpy; 32 clumps are identified. The identified clumps are approximately gravitationally bound and consistent with a standard isothermal sphere density, which suggests that they are likely collapsing protostellar cores. The gas kinematics suggests that velocity gradients exist that imply effects from the inclination of the cavity and MWC 1080. The kinematics of dense gas has also been affected by either outflows or stellar winds from MWC 1080, and lower mass clumps are possibly under stronger effects from MWC 1080 than higher mass clumps. In addition, low-mass cluster members tend to be formed in the denser and more turbulent cores, compared to isolated low-mass star-forming cores. This results from the contributions of nearby forming massive stars, such as outflows or stellar winds. Therefore, we conclude that in clusters like the MWC 1080 system, effects from massive stars dominate the star-forming environment in both the kinematics and dynamics of the natal cloud and the formation of low-mass cluster members. This study provides insights into the effects of MWC 1080 on its natal cloud, and suggests a different low-mass star-forming environment in clusters compared to isolated star formation.

A deeply embedded young protoplanetary disk around L1489 IRS observed by the Submillimeter Array

Circumstellar disks are expected to form early in the process that leads to the formation of a young star, during the collapse of the dense molecular cloud core. It is currently not well understood at what stage of the collapse the disk is formed or how it subsequently evolves. We aim to identify whether an embedded Keplerian protoplanetary disk resides in the L1489 IRS system. Given the amount of envelope material still present, such a disk would respresent a very young example of a protoplanetary disk. Using the Submillimeter Array (SMA) we have observed the HCO$^+$ $J=$ 3--2 line with a resolution of about 1$''$. At this resolution a protoplanetary disk with a radius of a few hundred AUs should be detectable, if present. Radiative transfer tools are used to model the emission from both continuum and line data. We find that these data are consistent with theoretical models of a collapsing envelope and Keplerian circumstellar disk. Models reproducing both the SED and the interferometric continuum observations reveal that the disk is inclined by 40$^\circ$ which is significantly different to the surrounding envelope (74$^\circ$). This misalignment of the angular momentum axes may be caused by a gradient within the angular momentum in the parental cloud or if L1489 IRS is a binary system rather than just a single star. In the latter case, future observations looking for variability at sub-arcsecond scales may be able to constrain these dynamical variations directly. However, if stars form from turbulent cores, the accreting material will not have a constant angular momentum axis (although the average is well defined and conserved) in which case it is more likely to have a misalignment of the angular momentum axes of the disk and the envelope.

Massive star and star cluster formation

AbstractI review the status of massive star formation theories: accretion from collapsing, massive, turbulent cores; competitive accretion; and stellar collisions. I conclude the observational and theoretical evidence favors the first of these models. I then discuss: the initial conditions of star cluster formation as traced by infrared dark clouds; the cluster formation timescale; and comparison of the initial cluster mass function in different galactic environments.

High-Resolution Studies of the Dense Molecular Cores toward Massive Star-Forming Regions

We present the results of the imaging of dense molecular cores in three massive star-forming regions (IRAS 02461$+$6147, IRAS 03035$+$5819, and IRAS 06058$+$2138), which are associated with luminous infrared sources with $1 \hbox{--} 5 \times 10^{4} \,{{{L}_{\odot}}}$. The C${{18\atop} \mathrm{O}}$ ($J = 1 \hbox{--} 0$) molecular emission from the dense cores and the thermal dust and/or free-free emission at 98 GHz and 110 GHz in the hot core candidates have been imaged with a resolution of $\sim {4{}^{\mathrm {\prime \prime }}}$ using the Nobeyama Millimeter Array. We identified 28 C${{18\atop} \mathrm{O}}$ dense cores, whose mass, effective radius, and line-width range from 2.1 to 29${M}_{\odot}$, from 0.013 to 0.108 pc, and from 0.7 to $2.7 \,\mathrm{km} \,\mathrm{s}^{-1}$, respectively, and identified four thermal dust millimeter continuum sources (MCSs). We divided the C${{18\atop} \mathrm{O}}$ cores into two types, a turbulent core and a non-turbulent core. The non-turbulent cores are similar to the typical cores in the low-mass star-forming regions. The turbulent cores have a higher average $\mathrm{H}_2$ density than those of the non-turbulent cores, and the external pressure of these cores is 100–1000 times higher than that of low-mass star-forming regions. Three of the turbulent cores are associated with massive protostar candidates and the intensity peak of MCS. Massive stars are formed from the turbulent cores which are gravitationally bound. In order to form such a turbulent core a molecular cloud has a large kinetic motion and a large mass.

Centrally condensed turbulent cores: massive stars or fragmentation?

We present numerical investigations into the formation of massive stars from centrally condensed turbulent cores. The results of five hydrodynamical simulations are described, following the collapse of the core, fragmentation and the formation of small clusters of protostars. We generate two different initial turbulent velocity fields corresponding to power-law spectra with power varying as k^{-4} and k^{-3.5}, and apply two different initial core radii. Calculations are included for both completely isothermal collapse, and a non-isothermal equation of state above a critical density 10^{-14}gcm^{-3}. Our calculations reveal the preference of fragmentation over monolithic star formation in turbulent cores. Fragmentation was prevalent in all the isothermal cases. Although disc fragmentation was largely suppressed in the non-isothermal runs due to the small dynamic range between the initial density and the critical density, our results show that some fragmentation still persisted. This is inconsistent with previous suggestions that turbulent cores result in the formation of a single massive star. We conclude that turbulence cannot be measured as an isotropic pressure term.