Dust grains play a major role in many astrophysical contexts. They affect the chemical, magnetic, dynamical, and optical properties of their environment, from galaxies down to the interstellar medium, star-forming regions, and protoplanetary disks. Their coagulation leads to shifts in their size distribution and ultimately to the formation of planets. However, although the coagulation process is reasonably uncomplicated to numerically implement by itself, it is difficult to couple it with multidimensional hydrodynamics numerical simulations because of its high computational cost. We propose here a simple method for tracking the coagulation of grains at far lower cost. Given an initial grain size distribution, the state of the distribution at time t is solely determined by the value of a single variable integrated along the trajectory, independently of the specific path taken by the grains. Although this method cannot account for processes other than coagulation, it is mathematically exact, fast, inexpensive, and can be used to evaluate the effect of grain coagulation in most astrophysical contexts. It is applicable to all coagulation kernels in which local physical conditions and grain properties can be separated. We also describe another method for calculating the average electric charge of grains and the density of ions and electrons in environments that are shielded from radiation fields, given the density and temperature of the gas, the cosmic-ray ionization rate, and the average mass of the ions. The equations we provide are fast to integrate numerically and can be used in multidimensional numerical simulations to self-consistently calculate on the fly the local resistivities that are required to model nonideal magnetohydrodynamics.
more »
« less
Fast methods for tracking grain coagulation and ionization: II. Extension to thermal ionization
Thermal ionization is a critical process at temperatures T > 10 3 K, particularly during star formation. An increase in ionization leads to a decrease in nonideal magnetohydrodynamics (MHD) resistivities, which has a significant impact on protoplanetary disks and protostar formation. We developed an extension of the fast computational ionization method presented in our recent paper to include thermal ionization. The model can be used to inexpensively calculate the density of ions and electrons and the electric charge of each size of grains for an arbitrary size distribution. This tool should be particularly useful for the self-consistent calculation of nonideal MHD resistivities in multidimensional simulations, especially of protostellar collapse and protoplanetary disks.
more »
« less
- Award ID(s):
- 1815461
- NSF-PAR ID:
- 10421732
- Date Published:
- Journal Name:
- Astronomy & Astrophysics
- Volume:
- 666
- ISSN:
- 0004-6361
- Page Range / eLocation ID:
- A27
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The water snowline in circumstellar disks is a crucial component in planet formation, but direct observational constraints on its location remain sparse owing to the difficulty of observing water in both young embedded and mature protoplanetary disks. Chemical imaging provides an alternative route to locate the snowline, and HCO + isotopologues have been shown to be good tracers in protostellar envelopes and Herbig disks. Here we present ∼0.″5 resolution (∼35 au radius) Atacama Large Millimeter/submillimeter Array (ALMA) observations of HCO + J = 4 − 3 and H 13 CO + J = 3 − 2 toward the young (Class 0/I) disk L1527 IRS. Using a source-specific physical model with the midplane snowline at 3.4 au and a small chemical network, we are able to reproduce the HCO + and H 13 CO + emission, but for HCO + only when the cosmic-ray ionization rate is lowered to 10 −18 s −1 . Even though the observations are not sensitive to the expected HCO + abundance drop across the snowline, the reduction in HCO + above the snow surface and the global temperature structure allow us to constrain a snowline location between 1.8 and 4.1 au. Deep observations are required to eliminate the envelope contribution to the emission and to derive more stringent constraints on the snowline location. Locating the snowline in young disks directly with observations of H 2 O isotopologues may therefore still be an alternative option. With a direct snowline measurement, HCO + will be able to provide constraints on the ionization rate.more » « less
-
Dust grains influence many aspects of star formation, including planet formation and the opacities for radiative transfer, chemistry, and the magnetic field via Ohmic, Hall, as well as ambipolar diffusion. The size distribution of the dust grains is the primary characteristic influencing all these aspects. Grain size increases by coagulation throughout the star formation process. In this work, we describe numerical simulations of protostellar collapse using methods described in earlier papers of this series. We compute the evolution of the grain size distribution from coagulation and the non-ideal magnetohydrodynamics effects self-consistently and at low numerical cost. We find that the coagulation efficiency is mostly affected by the time spent in high-density regions. Starting from sub-micron radii, grain sizes reach more than 100 µm in an inner protoplanetary disk that is only 1000 yr old. We also show that the growth of grains significantly affects the resistivities, while also having an indirect effect on the dynamics and angular momentum of the disk.more » « less
-
null (Ed.)ABSTRACT Non-ideal magnetohydrodynamic (MHD) effects have been shown recently as a robust mechanism of averting the magnetic braking ‘catastrophe’ and promoting protostellar disc formation. However, the magnetic diffusivities that determine the efficiency of non-ideal MHD effects are highly sensitive to microphysics. We carry out non-ideal MHD simulations to explore the role of microphysics on disc formation and the interplay between ambipolar diffusion (AD) and Hall effect during the protostellar collapse. We find that removing the smallest grain population (≲10 nm) from the standard MRN size distribution is sufficient for enabling disc formation. Further varying the grain sizes can result in either a Hall-dominated or an AD-dominated collapse; both form discs of tens of au in size regardless of the magnetic field polarity. The direction of disc rotation is bimodal in the Hall-dominated collapse but unimodal in the AD-dominated collapse. We also find that AD and Hall effect can operate either with or against each other in both radial and azimuthal directions, yet the combined effect of AD and Hall is to move the magnetic field radially outward relative to the infalling envelope matter. In addition, microphysics and magnetic field polarity can leave profound imprints both on observables (e.g. outflow morphology, disc to stellar mass ratio) and on the magnetic field characteristics of protoplanetary discs. Including Hall effect relaxes the requirements on microphysics for disc formation, so that prestellar cores with cosmic ray ionization rate of ≲2–3 × 10−16 s−1 can still form small discs of ≲10 au radius. We conclude that disc formation should be relatively common for typical prestellar core conditions, and that microphysics in the protostellar envelope is essential to not only disc formation, but also protoplanetary disc evolution.more » « less
-
Abstract Giant planets have been discovered at large separations from the central star. Moreover, a striking number of young circumstellar disks have gas and/or dust gaps at large orbital separations, potentially driven by embedded planetary objects. To form massive planets at large orbital separations through core accretion within the disk lifetime, however, an early solid body to seed pebble and gas accretion is desirable. Young protoplanetary disks are likely self-gravitating, and these gravitoturbulent disks may efficiently concentrate solid material at the midplane driven by spiral waves. We run 3D local hydrodynamical simulations of gravitoturbulent disks with Lagrangian dust particles to determine whether particle and gas self-gravity can lead to the formation of dense solid bodies, seeding later planet formation. When self-gravity between dust particles is included, solids of size St = 0.1–1 concentrate within the gravitoturbulent spiral features and collapse under their own self-gravity into dense clumps up to several M ⊕ in mass at wide orbits. Simulations with dust that drift most efficiently, St = 1, form the most massive clouds of particles, while simulations with smaller dust particles, St = 0.1, have clumps with masses an order of magnitude lower. When the effect of dust backreaction onto the gas is included, dust clumps become smaller by a factor of a few but more numerous. The existence of large solid bodies at an early stage of the disk can accelerate the planet formation process, particularly at wide orbital separations, and potentially explain planets distant from the central stars and young protoplanetary disks with substructures.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1051/0004-6361/202142551
