Particle Dynamics in 3D Self-gravitating Disks I: Spirals

Abstract

Spiral arms are distinctive features of many circumstellar disks, observed in scattered light, which traces the disk surface, millimeter dust emission, which probes the disk midplane, as well as molecular emission. The two leading explanations for spirals are wakes generated by a massive planet and the density waves excited by disk self-gravity. We use stratified 3D hydrodynamic shearing-box simulations including dust particles and disk self-gravity to investigate how gas and dust spirals in a self-gravitating disk depend on the simulation size, the cooling efficiency, and the aerodynamics properties of particles. We find that opening angles of spirals are universal ($\sim 10^o$), and not significantly affected by the size of the computational domain, the cooling time, or the particle size. In simulations with the biggest domain, the spirals in the gaseous disk become slightly more open with a higher cooling efficiency. Small dust follows the gaseous spirals very well, while intermediate-sized dust with dimensionless stopping time $(\mathrm{St})$ close to 1 concentrates to the spirals more and shows stronger spirals. However, large dust with $\mathrm{St}>1$ also shows spirals, which is different from some previous simulations. We identify that this is due to the gravity from the gas to the dust component. We show that when $\mathrm{St}\gtrsim Q$, the gravitational force from the gaseous spirals to the dust particles becomes stronger than the particles' aerodynamic drag force, so that the gas significantly affects these large particles through gravitational interaction. This has important implications for both spiral observations and planetesimal formation/dynamics.

15 pages, 7 figures, Accepted for publication in ApJ

References in corpus (18)

• The NumPy array: a structure for efficient numerical computation
• Spiral Density Waves in a Young Protoplanetary Disk
• Asymmetric features in the protoplanetary disk MWC758
• Shadows and spirals in the protoplanetary disk HD 100453
• Characterising the Gravitational Instability in Cooling Accretion Discs
• The Disk Substructures at High Angular Resolution Project (DSHARP): III. Spiral Structures in the Millimeter Continuum of the Elias 27, IM Lup, and WaOph 6 Disks
• Numerical requirements for simulations of self gravitating and non-self gravitating disks
• Global magnetohydrodynamical models of turbulence in protoplanetary disks I. A cylindrical potential on a Cartesian grid and transport of solids
• Vortices and spirals in the HD135344B transition disk
• Spiral Structure in the Gas Disk of TW Hya
• Spiral arms and instability within the AFGL 4176 mm1 disc
• The fragmentation criteria in local vertically stratified self-gravitating disk simulations
• Large-scale CO spiral arms and complex kinematics associated with the T Tauri star RU Lup
• Characterizing gravito-turbulence in 3D: turbulent properties and stability against fragmentation
• Kinematic Analysis of a Protostellar Multiple System: Measuring the Protostar Masses and Assessing Gravitational Instability in the Disks of L1448 IRS3B and L1448 IRS3A
• A Tight Relation between Spiral Arm Pitch Angle and Protoplanetary Disk Mass
• The observational impact of dust trapping in self-gravitating discs
• Global Spiral Density Wave Modes in Protoplanetary Disks: Morphology of Spiral Arms

Cited by in corpus (4)

• Particle Dynamics in 3D Self-gravitating Disks II: Strong Gas Accretion and Thin Dust Disks
• A Simple Time-Dependent Method for Calculating Spirals: Applications to Eccentric Planets in Protoplanetary discs
• Distributions of gas and small and large grains in the LkH$\alpha 330$ disk trace a young planetary system
• Continuing to Hide Signatures of Gravitational Instability in Protoplanetary Discs with Planets