skip to main content

Explore Research Products in the NSF-PAR

Title: Fast methods for tracking grain coagulation and ionization: III. Protostellar collapse with non-ideal MHD
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
Award ID(s):
1815461
NSF-PAR ID:
10421740
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Astronomy & Astrophysics
Volume:
670
ISSN:
0004-6361
Page Range / eLocation ID:
A61
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ( , Astronomy & Astrophysics)
    null (Ed.)
    Context. The degree of coupling between the gas and the magnetic field during the collapse of a core and the subsequent formation of a disk depends on the assumed dust size distribution. Aims. We study the impact of grain–grain coagulation on the evolution of magnetohydrodynamic (MHD) resistivities during the collapse of a prestellar core. Methods. We use a 1D model to follow the evolution of the dust size distribution, out-of-equilibrium ionisation state, and gas chemistry during the collapse of a prestellar core. To compute the grain–grain collisional rate, we consider models for both random and systematic, size-dependent, velocities. We include grain growth through grain–grain coagulation and ice accretion, but ignore grain fragmentation. Results. Starting with a Mathis-Rumpl-Nordsieck (MRN) size distribution (Mathis et al. 1977, ApJ, 217, 425), we find that coagulation in grain–grain collisions generated by hydrodynamical turbulence is not efficient at removing the smallest grains and, as a consequence, does not have a large effect on the evolution of the Hall and ambipolar diffusion MHD resistivities, which still drop significantly during the collapse like in models without coagulation. The inclusion of systematic velocities, possibly induced by the presence of ambipolar diffusion, increases the coagulation rate between small and large grains, removing small grains earlier in the collapse and therefore limiting the drop in the Hall and ambipolar diffusion resistivities. At intermediate densities ( n H ~ 10 8 cm −3 ), the Hall and ambipolar diffusion resistivities are found to be higher by 1 to 2 orders of magnitude in models with coagulation than in models where coagulation is ignored, and also higher than in a toy model without coagulation where all grains smaller than 0.1 μ m would have been removed in the parent cloud before the collapse. Conclusions. When grain drift velocities induced by ambipolar diffusion are included, dust coagulation happening during the collapse of a prestellar core starting from an initial MRN dust size distribution appears to be efficient enough to increase the MHD resistivities to the values necessary to strongly modify the magnetically regulated formation of a planet-forming disk. A consistent treatment of the competition between fragmentation and coagulation is, however, necessary before reaching firm conclusions. 
    more » « less
  2. ; ; ; ( , Astronomy & Astrophysics)
    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
  3. ; ; ; ; ; ; ; ( , The Astrophysical Journal)
    Abstract The alignment of dust grains with the ambient magnetic field produces polarization of starlight as well as thermal dust emission. Using the archival SOFIA/HAWC+ polarimetric data observed toward the ρ Ophiuchus (Oph) A cloud hosted by a B star at 89 and 154 μ m, we find that the fractional polarization of thermal dust emission first increases with the grain temperature and then decreases once the grain temperature exceeds ≃25–32 K. The latter trend differs from the prediction of the popular RAdiative Torques (RATs) alignment theory, which implies a monotonic increase of the polarization fraction with the grain temperature. We perform numerical modeling of polarized dust emission for the ρ Oph-A cloud and calculate the degree of dust polarization by simultaneously considering the dust grain alignment and rotational disruption by RATs. Our modeling results could successfully reproduce both the rising and declining trends of the observational data. Moreover, we show that the alignment of only silicate grains or a mixture of silicate–carbon grains within a composite structure can reproduce the observational trends, assuming that all dust grains follow a power-law size distribution. Although there are a number of simplifications and limitations to our modeling, our results suggest grains in the ρ Oph-A cloud have a composite structure, and the grain size distribution has a steeper slope than the standard size distribution for the interstellar medium. Combination of SOFIA/HAWC+ data with JCMT observations 450 and 850  μ m would be useful to test the proposed scenario based on grain alignment and disruption by RATs. 
    more » « less
  4. ; ; ; ; ; ; ; ; ; ; et al ( , Astronomy & Astrophysics)

    Context.T Tauri stars are low-mass young stars whose disks provide the setting for planet formation, which is one of the most fundamental processes in astronomy. Yet the mechanisms of this are still poorly understood. SU Aurigae is a widely studied T Tauri star and here we present original state-of-the-art interferometric observations with better uv and baseline coverage than previous studies.

    Aims.We aim to investigate the characteristics of the circumstellar material around SU Aur, and constrain the disk geometry, composition and inner dust rim structure.

    Methods.The MIRC-X instrument at CHARA is a six-telescope optical beam combiner offering baselines up to 331 m. We undertook image reconstruction for model-independent analysis, and fitted geometric models such as Gaussian and ring distributions. Additionally, the fitting of radiative transfer models constrained the physical parameters of the disk.

    Results.Image reconstruction reveals a highly inclined disk with a slight asymmetry consistent with inclination effects obscuring the inner disk rim through absorption of incident star light on the near side and thermal re-emission/scattering of the far side. Geometric models find that the underlying brightness distribution is best modelled as a Gaussian with a Full-Width Half-Maximum of 1.53 ± 0.01 mas at an inclination of 56.9 ± 0.4° and a minor axis position angle of 55.9 ± 0.5°. Radiative transfer modelling shows a flared disk with an inner radius at 0.16 au which implies a grain size of 0.14 μm assuming astronomical silicates and a scale height of 9.0 au at 100 au. In agreement with the literature, only the dusty disk wind successfully accounts for the near infrared excess by introducing dust above the mid-plane.

    Conclusions.Our results confirm and provide better constraints than previous inner disk studies of SU Aurigae. We confirm the presence of a dusty disk wind in the cicumstellar environment, the strength of which is enhanced by a late infall event which also causes very strong misalignments between the inner and outer disks.

     
    more » « less
  5. ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    ABSTRACT

    In dusty cool-star outflow or ejection events around asymptotic giant branch (AGB) or R Coronae Borealis or RCB-like stars, dust is accelerated by radiation from the star and coupled to the gas via collisional drag forces. It has recently been shown that such dust-gas mixtures are unstable to a super-class of instabilities called the resonant drag instabilities (RDIs), which promote dust clustering. We therefore consider idealized simulations of the RDIs operating on a spectrum of dust grain sizes subject to radiative acceleration (allowing for different grain optical properties), coupled to the gas with a realistic drag law, including or excluding the effects of magnetic fields and charged grains, and calculate for the first time how the RDIs could contribute to observed variability. We show that the RDIs naturally produce significant variations (spatially and temporally) ($\sim 10\!-\!20{{\ \rm per\ cent}}$ 1 σ-level) in the extinction, corresponding to $\sim 0.1\!-\!1\,$mag level in the stellar types above, on time-scales of order months to a year. The fluctuations are surprisingly robust to the assumed size of the source as they are dominated by large-scale modes, which also means their spatial structure could be resolved in some nearby systems. We also quantify how this produces variations in the line-of-sight grain size-distribution. All of these variations are similar to those observed, suggesting that the RDIs may play a key role driving observed spatial and temporal variability in dust extinction within dusty outflow/ejection events around cool stars. We further propose that the measured variations in grain sizes could directly be used to identify the presence of the RDIs in close by systems with observations.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email