Convergence of the critical cooling rate for protoplanetary disk fragmentation achieved; the key role of numerical dissipation of angular momentum

Abstract

We carry out simulations of gravitationally unstable disks using smoothed particle hydrodynamics(SPH) and the novel Lagrangian meshless finite mass (MFM) scheme in the GIZMO code (Hopkins 2015). Our aim is to understand the cause of the non-convergence of the cooling boundary for fragmentation reported in the literature. We run SPH simulations with two different artificial viscosity implementations, and compare them with MFM, which does not employ any artificial viscosity. With MFM we demonstrate convergence of the critical cooling time scale for fragmentation at β_{crit} =3.. Non-convergence persists in SPH codes, although it is significantly mitigated with schemes having reduced artificial viscosity such as inviscid SPH (ISPH) (Cullen & Dehnen 2010). We show how the non-convergence problem is caused by artificial fragmentation triggered by excessive dissipation of angular momentum in domains with large velocity derivatives. With increased resolution such domains become more prominent. Vorticity lags behind density due to numerical viscous dissipation in these regions, promoting collapse with longer cooling times. Such effect is shown to be dominant over the competing tendency of artificial viscosity to diminish with increasing resolution. When the initial conditions are first relaxed for several orbits, the flow is more regular, with lower shear and vorticity in non-axisymmetric regions, aiding convergence. Yet MFM is the only method that converges exactly. Our findings are of general interest as numerical dissipation via artificial viscosity or advection errors can also occur in grid-based codes. Indeed for the FARGO code values of β_{crit} significantly higher than our converged estimate have been reported in the literature. Finally, we discuss implications for giant planet formation via disk instability.

please check the ApJ version

References in corpus (13)

• Modelling discontinuities and Kelvin-Helmholtz instabilities in SPH
• Inviscid SPH
• The Two Modes of Gas Giant Planet Formation
• Characterising the Gravitational Instability in Cooling Accretion Discs
• A Superbubble Feedback Model for Galaxy Simulations
• Hydrodynamics of Embedded Planets' First Atmospheres. II. A Rapid Recycling of Atmospheric Gas
• Properties of gravitoturbulent accretion disks
• Numerical requirements for simulations of self gravitating and non-self gravitating disks
• On the gap-opening criterion of migrating planets in protoplanetary disks
• Hydrodynamics of Embedded Planets' First Atmospheres. I. A Centrifugal Growth Barrier for 2D Flows
• Circumplanetary disks around young giant planets: a comparison between core-accretion and disk instability
• The response of self-graviting protostellar discs to slow reduction in cooling timescale: the fragmentation boundary revisited
• Fragmentation of protoplanetary disks around M-dwarfs

Cited by in corpus (35)

• On the diversity and statistical properties of protostellar discs
• STARFORGE: Toward a comprehensive numerical model of star cluster formation and feedback
• The statistical properties of stars and their dependence on metallicity
• Towards a population synthesis model of self-gravitating disc fragmentation and tidal downsizing II: The effect of fragment-fragment interactions
• In-situ extreme mass ratio inspirals via sub-parsec formation and migration of stars in thin, gravitationally unstable AGN discs
• Formation of intermediate-mass planets via magnetically-controlled disk fragmentation
• Rapid Formation of Jupiter and Wide-Orbit Exoplanets in Disks with Pressure Bumps
• Chaotic Type I Migration in Turbulent Discs
• Characterizing gravito-turbulence in 3D: turbulent properties and stability against fragmentation
• The influence of infall on the properties of protoplanetary discs
• Chemulator: Fast, accurate thermochemistry for dynamical models through emulation
• Local simulations of MRI turbulence with meshless methods
• Giant planets and brown dwarfs on wide orbits: a code comparison project
• Direct Formation of Planetary Embryos in Self-Gravitating Disks
• Parametric instability in a free evolving warped protoplanetary disc
• A two-step gravitational cascade for the fragmentation of self-gravitating discs
• The role of the drag force in the gravitational stability of dusty planet-forming disc -- II. Numerical simulations
• Fragmentation favoured in discs around higher mass stars
• Thermal evolution of protoplanetary disks: from $\beta$-cooling to decoupled gas and dust temperatures
• Characterizing fragmentation and sub-Jovian clump properties in magnetized young protoplanetary disks
• The link between infall location, early disc size, and the fraction of self-gravitationally fragmenting discs
• Gravito-turbulence in local disk simulations with an adaptive moving mesh
• Turbulent Transport of Dust Particles in Protostellar Disks: The Effect of Upstream Diffusion
• Formation of free-floating planetary mass objects via circumstellar disk encounters
• Parametric instability in warped astrophysical discs: growth, saturation and feedback
• The Effect of the Approach to Gas Disk Gravitational Instability on the Rapid Formation of Gas Giant Planets. II. Quadrupled Spatial Resolution
• Hypothetical hyperbolic encounters between Venus and proto-Mercury that partially \ stripped away proto-Mercury's mantle
• Standing solitary waves as transitions to spiral structures in gravitationally unstable accretion disks
• Angular momentum transport via gravitational instability in the Elias 2-27 disc
• DIPSY: A new Disc Instability Population SYnthesis, I. Modeling, evolution of individual systems, and tests
• Collisions of young disc galaxies in the early universe
• Dust growth and planet formation by disc fragmentation
• A novel sub-grid model for super-Eddington accretion of spinning black holes in galaxy-scale simulations
• Gravitational instability in planet-forming discs
• Nonlinear behaviour of warped discs around a central object with a quadrupole moment