Convergence of simulations of self-gravitating accretion discs II: Sensitivity to the implementation of radiative cooling and artificial viscosity

Abstract

Recently it has been suggested that the fragmentation boundary in Smoothed Particle Hydrodynamic (SPH) and FARGO simulations of self-gravitating accretion discs with beta-cooling do not converge as resolution is increased. Furthermore, this recent work suggests that by carefully optimising the artificial viscosity parameters in these codes it can be shown that fragmentation may occur for much longer cooling times than earlier work suggests. If correct, this result is intriguing as it suggests that gas giant planets could form, via direct gravitational collapse, reasonably close to their parent stars. This result is, however, slightly surprising and there have been a number of recent studies suggesting that the result is likely an indication of a numerical problem with the simulations. One suggestion, in particular, is that the SPH results are influenced by the manner in which the cooling is implemented. We extend this work here and show that if the cooling is implemented in a manner that removes a known numerical artefact in the shock regions, the fragmentation boundary converges to a value consistent with earlier work and that fragmentation is unlikely for the long cooling times suggested by this recent work. We also investigate the optimisation of the artificial viscosity parameters and show that the values that appear optimal are likely introducing numerical problems in both the SPH and FARGO simulations. We therefore conclude that earlier predictions for the cooling times required for fragmentation are likely correct and that, as suggested by this earlier work, fragmentation cannot occur in the inner parts (r < 50 au) of typical protostellar discs.

11 pages, 7 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

References in corpus (6)

• Modelling discontinuities and Kelvin-Helmholtz instabilities in SPH
• MHD simulations of the magnetorotational instability in a shearing box with zero net flux. II. The effect of transport coefficients
• Characterising the Gravitational Instability in Cooling Accretion Discs
• Star Formation Around Super-Massive Black Holes
• Can giant planets form by gravitational fragmentation of discs?
• 3D Radiative Hydrodynamics for Disk Stability Simulations: A Proposed Testing Standard and New Results

Cited by in corpus (38)

• On the diversity and statistical properties of protostellar discs
• The SPHERE infrared survey for exoplanets (SHINE). III. The demographics of young giant exoplanets below 300 au with SPHERE
• The VLT/NaCo large program to probe the occurrence of exoplanets and brown dwarfs at wide orbits. IV. Gravitational instability rarely forms wide, giant planets
• 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
• Forming spectroscopic massive proto-binaries by disk fragmentation
• On the gap-opening criterion of migrating planets in protoplanetary disks
• Convergence of the critical cooling rate for protoplanetary disk fragmentation achieved; the key role of numerical dissipation of angular momentum
• A revised condition for self-gravitational fragmentation of protoplanetary disks
• The fragmentation criteria in local vertically stratified self-gravitating disk simulations
• Grand challenges in protoplanetary disc modelling
• Gravito-Turbulent Disks in 3D: Turbulent Velocities vs. Depth
• Viscosity prescription for gravitationally unstable accretion disks
• On shocks driven by high-mass planets in radiatively inefficient disks. II. Three-dimensional global disk simulations
• Smoothed particle hydrodynamics simulations of gas and dust mixtures
• Triggered fragmentation in self-gravitating discs: forming fragments at small radii
• Ejection of gaseous clumps from gravitationally unstable protostellar disks
• Gravito-turbulence and the excitation of small-scale parametric instability in astrophysical discs
• On the fragmentation boundary in magnetised self-gravitating discs
• Quantification of stochastic fragmentation of self-gravitating discs
• Local simulations of MRI turbulence with meshless methods
• The Effect of Protoplanetary Disk Cooling Times on the Formation of Gas Giant Planets by Gravitational Instability
• Temperature Fluctuations driven by Magnetorotational Instability in Protoplanetary Disks
• On shocks driven by high-mass planets in radiatively inefficient disks. I. Two-dimensional global disk simulations
• Gravitational stability and fragmentation condition for discs around accreting supermassive stars
• A two-step gravitational cascade for the fragmentation of self-gravitating discs
• Fragmentation favoured in discs around higher mass stars
• The evolution of self-gravitating accretion discs
• The Role of the Cooling Prescription for Disk Fragmentation: Numerical Convergence & Critical Cooling Parameter in Self-Gravitating Disks
• Fragmentation of protoplanetary disks around M-dwarfs
• Gravitoturbulence in magnetised protostellar discs
• On the origin of the central 1" hole in the stellar disk of Sgr A* and the Fermi gamma-ray bubbles
• The impact of numerical oversteepening on the fragmentation boundary in self-gravitating disks
• Circumstellar discs: What will be next?
• Structure of Protoplanetary Discs with Magnetically-driven Winds
• On fragmentation of turbulent self-gravitating discs in the long cooling time regime
• Tensor classification of structure in smoothed particle hydrodynamics density fields
• A theoretical perspective on the formation and fragmentation of protostellar discs