More Like this

1. Simulating diverse instabilities of dust in magnetized gas

   https://doi.org/10.1093/mnras/staa1046  

   Hopkins, Philip F ; Squire, Jonathan ; Seligman, Darryl ( August 2020 , Monthly Notices of the Royal Astronomical Society)

   null (Ed.)

   ABSTRACT Recently, Squire & Hopkins showed that charged dust grains moving through magnetized gas under the influence of a uniform external force (such as radiation pressure or gravity) are subject to a spectrum of instabilities. Qualitatively distinct instability families are associated with different Alfvén or magnetosonic waves and drift or gyro motion. We present a suite of simulations exploring these instabilities, for grains in a homogeneous medium subject to an external acceleration. We vary parameters such as the ratio of Lorentz-to-drag forces on dust, plasma β, size scale, and acceleration. All regimes studied drive turbulent motions and dust-to-gas fluctuations in the saturated state, rapidly amplify magnetic fields into equipartition with velocity fluctuations, and produce instabilities that persist indefinitely (despite random grain motions). Different parameters produce diverse morphologies and qualitatively different features in dust, but the saturated gas state can be broadly characterized as anisotropic magnetosonic or Alfvénic turbulence. Quasi-linear theory can qualitatively predict the gas turbulent properties. Turbulence grows from small to large scales, and larger scale modes usually drive more vigorous gas turbulence, but dust velocity and density fluctuations are more complicated. In many regimes, dust forms structures (clumps, filaments, sheets) that reach extreme overdensities (up to ≫109 times mean), and exhibit substantial substructure even in nearly incompressible gas. These can be even more prominent at lower dust-to-gas ratios. In other regimes, dust self-excites scattering via magnetic fluctuations that isotropize and amplify dust velocities, producing fast, diffusive dust motions. 

   more » « less
2. The acoustic resonant drag instability with a spectrum of grain sizes

   https://doi.org/10.1093/mnras/stab3377  

   Squire, Jonathan ; Moroianu, Stefania ; Hopkins, Philip F. ( November 2021 , Monthly Notices of the Royal Astronomical Society)

   We study the linear growth and non-linear saturation of the ‘acoustic Resonant Drag Instability’ (RDI) when the dust grains, which drive the instability, have a wide, continuous spectrum of different sizes. This physics is generally applicable to dusty winds driven by radiation pressure, such as occurs around red-giant stars, star-forming regions, or active galactic nuclei. Depending on the physical size of the grains compared to the wavelength of the radiation field that drives the wind, two qualitatively different regimes emerge. In the case of grains that are larger than the radiation’s wavelength – termed the constant-drift regime – the grain’s equilibrium drift velocity through the gas is approximately independent of grain size, leading to strong correlations between differently sized grains that persist well into the saturated non-linear turbulence. For grains that are smaller than the radiation’s wavelength – termed the non-constant-drift regime – the linear instability grows more slowly than the single-grain-size RDI and only the larger grains exhibit RDI-like behaviour in the saturated state. A detailed study of grain clumping and grain–grain collisions shows that outflows in the constant-drift regime may be effective sites for grain growth through collisions, with large collision rates but low collision velocities.

    

   more » « less
3. Streaming instability with multiple dust species – II. Turbulence and dust–gas dynamics at non-linear saturation

   https://doi.org/10.1093/mnras/stab2959  

   Yang, Chao-Chin ; Zhu, Zhaohuan ( October 2021 , Monthly Notices of the Royal Astronomical Society)

   ABSTRACT The streaming instability is a fundamental process that can drive dust–gas dynamics and ultimately planetesimal formation in protoplanetary discs. As a linear instability, it has been shown that its growth with a distribution of dust sizes can be classified into two distinct regimes, fast- and slow-growth, depending on the dust-size distribution and the total dust-to-gas density ratio ϵ. Using numerical simulations of an unstratified disc, we bring three cases in different regimes into non-linear saturation. We find that the saturation states of the two fast-growth cases are similar to its single-species counterparts. The one with maximum dimensionless stopping time τs,max = 0.1 and ϵ = 2 drives turbulent vertical dust–gas vortices, while the other with τs,max = 2 and ϵ = 0.2 leads to radial traffic jams and filamentary structures of dust particles. The dust density distribution for the former is flat in low densities, while the one for the latter has a low-end cut-off. By contrast, the one slow-growth case results in a virtually quiescent state. Moreover, we find that in the fast-growth regime, significant dust segregation by size occurs, with large particles moving towards dense regions while small particles remain in the diffuse regions, and the mean radial drift of each dust species is appreciably altered from the (initial) drag-force equilibrium. The former effect may skew the spectral index derived from multiwavelength observations and change the initial size distribution of a pebble cloud for planetesimal formation. The latter along with turbulent diffusion may influence the radial transport and mixing of solid materials in young protoplanetary discs. 

   more » « less
4. Numerical study of cosmic ray confinement through dust resonant drag instabilities

   https://doi.org/10.1093/mnras/stac895  

   Ji, Suoqing ; Squire, Jonathan ; Hopkins, Philip F ( April 2022 , Monthly Notices of the Royal Astronomical Society)

   ABSTRACT We investigate the possibility of cosmic ray (CR) confinement by charged dust grains through resonant drag instabilities (RDIs). We perform magnetohydrodynamic particle-in-cell simulations of magnetized gas mixed with charged dust and cosmic rays, with the gyro-radii of dust and GeV CRs on ∼au scales fully resolved. As a first study, we focus on one type of RDI wherein charged grains drift super-Alfvénically, with Lorentz forces strongly dominating over drag forces. Dust grains are unstable to the RDIs and form concentrated columns and sheets, whose scale grows until saturating at the simulation box size. Initially perfectly streaming CRs are strongly scattered by RDI-excited Alfvén waves, with the growth rate of the CR perpendicular velocity components equaling the growth rate of magnetic field perturbations. These rates are well-predicted by analytic linear theory. CRs finally become isotropized and drift at least at ∼vA by unidirectional Alfvén waves excited by the RDIs, with a uniform distribution of the pitch angle cosine μ and a flat profile of the CR pitch angle diffusion coefficient Dμμ around μ = 0, without the ‘90○ pitch angle problem.’ With CR feedback on the gas included, Dμμ decreases by a factor of a few, indicating a lower CR scattering rate, because the backreaction on the RDI from the CR pressure adds extra wave damping, leading to lower quasi-steady-state scattering rates. Our study demonstrates that the dust-induced CR confinement can be very important under certain conditions, e.g. the dusty circumgalactic medium around quasars or superluminous galaxies. 

   more » « less
5. Near-cancellation of up- and down-gradient momentum transport in forced magnetized shear-flow turbulence

   https://doi.org/arXiv:2208.03342  

   Tripathi, B. ; Fraser, A. E. ; Terry, P. W. ; Zweibel, E. G. ; Pueschel, M. J. ( October 2022 , Physics of plasmas)

   Visco-resistive magnetohydrodynamic turbulence, driven by a two-dimensional unstable shear layer that is maintained by an imposed body force, is examined by decomposing it into dissipationless linear eigenmodes of the initial profiles. The down-gradient momentum flux, as expected, originates from the large-scale instability. However, continual up-gradient momentum transport by large-scale linearly stable but nonlinearly excited eigenmodes is identified, and found to nearly cancel the down-gradient transport by unstable modes. The stable modes effectuate this by depleting the large-scale turbulent fluctuations via energy transfer to the mean flow. This establishes a physical mechanism underlying the long-known observation that coherent vortices formed from nonlinear saturation of the instability reduce turbulent transport and fluctuations, as such vortices are composed of both the stable and unstable modes, which are nearly equal in their amplitudes. The impact of magnetic fields on the nonlinearly excited stable modes is then quantified. Even when imposing a strong magnetic field that almost completely suppresses the instability, the up-gradient transport by the stable modes is at least two-thirds of the down-gradient transport by the unstable modes, whereas for weaker fields, this fraction reaches up to 98% . These effects are persistent with variations in magnetic Prandtl number and forcing strength. Finally, continuum modes are shown to be energetically less important, but essential for capturing the magnetic fluctuations and Maxwell stress. A simple analytical scaling law is derived for their saturated turbulent amplitudes. It predicts the fall-off rate as the inverse of the Fourier wavenumber, a property which is confirmed in numerical simulations. 

   more » « less