ABSTRACT The formation of protoplanetary disks during the collapse of molecular dense cores is significantly influenced by angular momentum transport, notably by the magnetic torque. In turn, the ...evolution of the magnetic field is determined by dynamical processes and non-ideal MHD effects such as ambipolar diffusion. Considering simple relations between various timescales characteristic of the magnetized collapse, we derive an expression for the early disk radius, where M is the total disk plus protostar mass, is the ambipolar diffusion coefficient, and Bz is the magnetic field in the inner part of the core. This is significantly smaller than the disks that would form if angular momentum was conserved. The analytical predictions are confronted against a large sample of 3D, non-ideal MHD collapse calculations covering variations of a factor 100 in core mass, a factor 10 in the level of turbulence, a factor 5 in rotation, and magnetic mass-to-flux over critical mass-to-flux ratios 2 and 5. The disk radius estimates are found to agree with the numerical simulations within less than a factor 2. A striking prediction of our analysis is the weak dependence of circumstellar disk radii upon the various relevant quantities, suggesting weak variations among class-0 disk sizes. In some cases, we note the onset of large spiral arms beyond this radius.
We present high-resolution (∼35 au) ALMA Band 6 1.3 mm dust polarization observations of IRAS 16293. These observations spatially resolve the dust polarization across the two protostellar sources and ...toward the filamentary structures between them. The dust polarization and inferred magnetic field have complicated structures throughout the region. In particular, we find that the magnetic field is aligned parallel to three filamentary structures. We characterize the physical properties of the filamentary structure that bridges IRAS 16293A and IRAS 16293B and estimate a magnetic field strength of 23-78 mG using the Davis-Chandrasekhar-Fermi method. We construct a toy model for the bridge material assuming that the young stars dominate the mass and gravitational potential of the system. We find that the expected gas flow to each star is of comparable order to the Alfvén speed, which suggests that the field may be regulating the gas flow. We also find that the bridging material should be depleted in ∼103 yr. If the bridge is part of the natal filament that formed the stars, then it must have accreted new material. Alternatively, the bridge could be a transient structure. Finally, we show that the 1.3 mm polarization morphology of the optically thick IRAS 16293B system is qualitatively similar to dust self-scattering. Based on similar polarization measurements at 6.9 mm, we propose that IRAS 16293B has produced a substantial population of large dust grains with sizes between 200 and 2000 m.
We report the results of radiation-magnetohydrodynamics calculations in the context of high-mass star formation, using for the first time a self-consistent model for photon emission (i.e., via ...thermal emission and in radiative shocks) and with the high resolution necessary to properly resolve magnetic braking effects and radiative shocks on scales <100 AU. We investigate the combined effects of magnetic field, turbulence, and radiative transfer on the early phases of the collapse and the fragmentation of massive dense cores. We identify a new mechanism that inhibits initial fragmentation of massive dense cores where magnetic field and radiative transfer interplay. We show that this interplay becomes stronger as the magnetic field strength increases. Magnetic braking is transporting angular momentum outward and is lowering the rotational support and is thus increasing the infall velocity. This enhances the radiative feedback owing to the accretion shock on the first core. We speculate that highly magnetized massive dense cores are good candidates for isolated massive star formation while moderately magnetized massive dense cores are more appropriate forming OB associations or small star clusters.
Context. Star formation begins with the gravitational collapse of a dense core inside a molecular cloud. As the collapse progresses, the centre of the core begins to heat up as it becomes optically ...thick. The temperature and density in the centre eventually reach high enough values where fusion reactions can ignite, and the protostar is born. This sequence of events entails many physical processes, of which radiative transfer is of paramount importance. Simulated collapsing cores without radiative transfer rapidly become thermally supported before reaching high enough temperatures and densities, preventing the formation of stars. Aims. Many simulations of protostellar collapse make use of a grey treatment of radiative transfer coupled to the hydrodynamics. However, interstellar gas and dust opacities present large variations as a function of frequency, which can potentially be overlooked by grey models and lead to significantly different results. In this paper, we follow up on a previous paper on the collapse and formation of Larson’s first core using multigroup radiation hydrodynamics (Paper I) by extending the calculations to the second phase of the collapse and the formation of Larson’s second core. Methods. We have made the use of a non-ideal gas equation of state as well as an extensive set of spectral opacities in a spherically symmetric fully implicit Godunov code to model all the phases of the collapse of a 0.1, 1, and 10 M⊙ cloud cores. Results. We find that, for an identical central density, there are only small differences between the grey and multigroup simulations. The first core accretion shock remains supercritical while the shock at the second core border is found to be strongly subcritical with all the accreted energy being transfered to the core. The size of the first core was found to vary somewhat in the different simulations (more unstable clouds form smaller first cores) while the size, mass, and temperature of the second cores are independent of initial cloud mass, size, and temperature. Conclusions. Our simulations support the idea of a standard (universal) initial second core size of ~3 × 10-3 AU and mass ~1.4 × 10-3 M⊙. The grey approximation for radiative transfer appears to perform well in one-dimensional simulations of protostellar collapse, most probably because of the high optical thickness of the majority of the protostar-envelope system. A simple estimate of the characteristic timescale of the second core suggests that the effects of using multigroup radiative transfer may be more important in the long-term evolution of the protostar.
We present high-resolution (∼30 au) ALMA Band 6 dust polarization observations of VLA 1623. The VLA 1623 data resolve compact ∼40 au inner disks around the two protobinary sources, VLA 1623-A and VLA ...1623-B, and also an extended ∼180 au ring of dust around VLA 1623-A. This dust ring was previously identified as a large disk in lower-resolution observations. We detect highly structured dust polarization toward the inner disks and the extended ring with typical polarization fractions 1.7% and 2.4%, respectively. The two components also show distinct polarization morphologies. The inner disks have uniform polarization angles aligned with their minor axes. This morphology is consistent with expectations from dust scattering. By contrast, the extended dust ring has an azimuthal polarization morphology not previously seen in lower-resolution observations. We find that our observations are well-fit by a static, oblate spheroid model with a flux-frozen, poloidal magnetic field. We propose that the polarization traces magnetic grain alignment likely from flux freezing on large scales and magnetic diffusion on small scales. Alternatively, the azimuthal polarization may be attributed to grain alignment by the anisotropic radiation field. If the grains are radiatively aligned, then our observations indicate that large (∼100 m) dust grains grow quickly at large angular extents. Finally, we identify significant proper motion of VLA 1623 using our observations and those in the literature. This result indicates that the proper motion of nearby systems must be corrected for when combining ALMA data from different epochs.
We present subarcsecond 1.3 mm continuum ALMA observations toward the Orion Molecular Cloud 1 South (OMC-1S) region, down to a spatial resolution of 74 au, which reveal a total of 31 continuum ...sources. We also present subarcsecond 7 mm continuum VLA observations of the same region, which allow further study of fragmentation down to a spatial resolution of 40 au. By applying a method of "mean surface density of companions" we find a characteristic spatial scale at ∼560 au, and we use this spatial scale to define the boundary of 19 "cores" in OMC-1S as groupings of millimeter sources. We find an additional characteristic spatial scale at ∼2800 au, which is the typical scale of the filaments in OMC-1S, suggesting a two-level fragmentation process. We measured the fragmentation level within each core and find a higher fragmentation toward the southern filament. In addition, the cores of the southern filament are also the densest cores (within 1100 au) in OMC-1S. This is fully consistent with previous studies of fragmentation at spatial scales one order of magnitude larger, and suggests that fragmentation down to 40 au seems to be governed by thermal Jeans processes in OMC-1S.
Astrophysical plasmas are subject to a tight connection between magnetic fields and the diffusion of particles, which leads to an anisotropic transport of energy. Under the fluid assumption, this ...effect can be reduced to an advection-diffusion equation, thereby augmenting the equations of magnetohydrodynamics. We introduce a new method for solving the anisotropic diffusion equation using an implicit finite-volume method with adaptive mesh refinement and adaptive time-stepping in the ramses code. We apply this numerical solver to the diffusion of cosmic ray energy and diffusion of heat carried by electrons, which couple to the ion temperature. We test this new implementation against several numerical experiments and apply it to a simple supernova explosion with a uniform magnetic field.
Abstract
Through the magnetic braking and the launching of protostellar outflows, magnetic fields play a major role in the regulation of angular momentum in star formation, which directly impacts the ...formation and evolution of protoplanetary disks and binary systems. The aim of this paper is to quantify those phenomena in the presence of nonideal magnetohydrodynamics effects, namely, the ohmic and ambipolar diffusion. We perform three-dimensional simulations of protostellar collapses varying the mass of the prestellar dense core, the thermal support (the
α
ratio), and the dust grain size distribution. The mass mostly influences the magnetic braking in the pseudo-disk, while the thermal support impacts the accretion rate and hence the properties of the disk. Removing the grains smaller than 0.1
μ
m in the Mathis–Rumpl–Nordsieck distribution enhances the ambipolar diffusion coefficient. Similar to previous studies, we find that this change in the distribution reduces the magnetic braking with an impact on the disk. The outflow is also significantly weakened. In either case, the magnetic braking largely dominates the outflow as a process to remove the angular momentum from the disk. Finally, we report a large ionic precursor to the outflow with velocities of several km s
−1
, which may be observable.
Context. Cosmic rays propagate through the galactic scales down to the smaller scales at which stars form. Cosmic rays are close to energy equipartition with the other components of the interstellar ...medium and can provide a support against gravity if pressure gradients develop. Aims. We study the propagation of cosmic rays within the turbulent and magnetised bi-stable interstellar gas. The conditions necessary for cosmic-ray trapping and cosmic-ray pressure gradient development are investigated. Methods. We derived an analytical value of the critical diffusion coefficient for cosmic-ray trapping within a turbulent medium, which follows the observed scaling relations. We then presented a numerical study using 3D simulations of the evolution of a mixture of interstellar gas and cosmic rays, in which turbulence is driven at varying scales by stochastic forcing within a box of 40 pc. We explored a large parameter space in which the cosmic-ray diffusion coefficient, the magnetisation, the driving scale, and the amplitude of the turbulence forcing, as well as the initial cosmic-ray energy density, vary. Results. We identify a clear transition in the interstellar dynamics for cosmic-ray diffusion coefficients below a critical value deduced from observed scaling relations. This critical diffusion depends on the characteristic length scale L of Dcrit ≃ 3.1 × 1023 cm2 s−1(L/1 pc)q+1, where the exponent q relates the turbulent velocity dispersion σ to the length scale as σ ~ Lq. Hence, in our simulations this transition occurs around Dcrit ≃ 1024–1025 cm2 s−1. The transition is recovered in all cases of our parameter study and is in very good agreement with our simple analytical estimate. In the trapped cosmic-ray regime, the induced cosmic-ray pressure gradients can modify the gas flow and provide a support against the thermal instability development. We discuss possible mechanisms that can significantly reduce the cosmic-ray diffusion coefficients within the interstellar medium. Conclusions. Cosmic-ray pressure gradients can develop and modify the evolution of thermally bi-stable gas for diffusion coefficients D0 ≤ 1025 cm2 s−1 or in regions where the cosmic-ray pressure exceeds the thermal one by more than a factor of ten. This study provides the basis for further works including more realistic cosmic-ray diffusion coefficients, as well as local cosmic-ray sources.
In spite of decades of theoretical efforts, the physical origin of the stellar initial mass function (IMF) is still debated. Particularly crucial is the question of what sets the peak of the ...distribution. To investigate this issue, we perform high-resolution numerical simulations with radiative feedback exploring, in particular, the role of the stellar and accretion luminosities. We also perform simulations with a simple effective equation of state (EOS), and we investigate 1000 solar-mass clumps having, respectively, 0.1 and 0.4 pc of initial radii. We found that most runs, both with radiative transfer or an EOS, present similar mass spectra with a peak broadly located around 0.3-0.5 M and a power-law-like mass distribution at higher masses. However, when accretion luminosity is accounted for, the resulting mass spectrum of the most compact clump tends to be moderately top-heavy. The effect remains limited for the less compact one, which overall remains colder. Our results support the idea that rather than the radiative stellar feedback, this is the transition from the isothermal to the adiabatic regime, which occurs at a gas density of about 1010 cm−3, that is responsible for setting the peak of the IMF. This stems from (i) the fact that extremely compact clumps for which the accretion luminosity has a significant influence are very rare and (ii) the luminosity problem, which indicates that the effective accretion luminosity is likely weaker than expected.