Resolution requirements for smoothed particle hydrodynamics simulations of self-gravitating accretion discs

Stimulated by recent results by Meru & Bate, we revisit the issue of resolution requirements for simulating self-gravitating accretion discs with smoothed particle hydrodynamics (SPH). We show that all the results by Meru & Bate are actually consistent if they are interpreted as driven by resolution effects, therefore implying that the resolution criterion for cooling gaseous discs is a function of the imposed cooling rate. We discuss two possible numerical origins of such dependence, which are both consistent with the limited number of available data. Our results tentatively indicate that convergence for current simulations is being reached for a number of SPH particles approaching 10 million (for a disc mass of the order of 10 per cent of the central object mass), which would set the critical cooling time for fragmentation at about 15Ω ―1 , roughly a factor of 2 larger than previously thought. More in general, we discuss the extent to which the large number of recent numerical results are reliable or not. We argue that those results that pertain to the dynamics associated with gravitational instabilities (such as the locality of angular momentum transport, and the relationship between density perturbation and induced stress) are robust, while those pertaining to the thermodynamics of the system (such as the determination of the critical cooling time for fragmentation) can be affected by poor resolution.

[1]  G. Lodato,et al.  The effects of opacity on gravitational stability in protoplanetary discs , 2009, 0910.0850.

[2]  P. J. Armitage,et al.  Investigating fragmentation conditions in self-gravitating accretion discs , 2005 .

[3]  Cambridge,et al.  Testing the locality of transport in self-gravitating accretion discs , 2004 .

[4]  Charles F. Gammie,et al.  Nonlinear Outcome of Gravitational Instability in Disks with Realistic Cooling , 2003, astro-ph/0312507.

[5]  G. Lodato,et al.  The response of self‐gravitating protostellar discs to slow reduction in cooling time‐scale: the fragmentation boundary revisited , 2007, 0708.0742.

[6]  Numerical requirements for simulations of self-gravitating and non-self-gravitating discs , 2006, astro-ph/0609493.

[7]  Richard H. Durisen,et al.  The Thermal Regulation of Gravitational Instabilities in Protoplanetary Disks. II. Extended Simulations with Varied Cooling Rates , 2005 .

[8]  A. Boss.  Convective Cooling of Protoplanetary Disks and Rapid Giant Planet Formation , 2004 .

[9]  M. Krumholz,et al.  ON THE DYNAMICS AND EVOLUTION OF GRAVITATIONAL INSTABILITY-DOMINATED DISKS , 2010, 1003.4513.

[10]  Joachim Stadel,et al.  The Evolution of Gravitationally Unstable Protoplanetary Disks: Fragmentation and Possible Giant Planet Formation , 2003, astro-ph/0310771.

[11]  J.C.B. Papaloizou,et al.  On the Dynamical Foundations of α Disks , 1999 .

[12]  Joachim Stadel,et al.  Gravitational instability in binary protoplanetary discs: new constraints on giant planet formation , 2004, astro-ph/0405502.

[13]  R. Rafikov,et al.  PROPERTIES OF GRAVITOTURBULENT ACCRETION DISKS , 2009, 0901.4739.

[14]  Cambridge,et al.  Testing the locality of transport in self-gravitating accretion discs — II. The massive disc case , 2005 .

[15]  P. Armitage,et al.  The role of disc self-gravity in the formation of protostars and protostellar discs , 2009, 0911.1202.

[16]  S. Michael,et al.  The Effects of Metallicity and Grain Size on Gravitational Instabilities in Protoplanetary Disks , 2005, astro-ph/0508354.

[17]  Charles F. Gammie,et al.  Nonlinear Outcome of Gravitational Instability in Cooling, Gaseous Disks , 2001, astro-ph/0101501.

[18]  Sergei Nayakshin.  Massive stars in subparsec rings around galactic centres , 2005 .

[19]  A. C. Boley,et al.  The Thermal Regulation of Gravitational Instabilities in Protoplanetary Disks. III. Simulations with Radiative Cooling and Realistic Opacities , 2006 .

[20]  Sergei Nayakshin,et al.  Simulations of star formation in a gaseous disc around Sgr A* - A failed active galactic nucleus , 2007, astro-ph/0701141.

[21]  J. E. Pringle,et al.  Accretion Discs in Astrophysics , 1981 .

[22]  C. Clarke,et al.  Pseudo-viscous modelling of self-gravitating discs and the formation of low mass ratio binaries , 2009, 0904.3549.

[23]  J. Monaghan.  Smoothed particle hydrodynamics , 2005 .

[24]  B. Macintosh,et al.  Direct Imaging of Multiple Planets Orbiting the Star HR 8799 , 2008, Science.

[25]  Giuseppe Lodato,et al.  Characterizing the gravitational instability in cooling accretion discs , 2008, 0811.3629.

[26]  James R. Murray,et al.  SPH simulations of tidally unstable accretion discs in cataclysmic variables , 1995, astro-ph/9511031.

[27]  S. Lubow,et al.  Mass Flow through Gaps in Circumbinary Disks , 1996 .

[28]  D. Balsara.  von Neumann stability analysis of smoothed particle hydrodynamics—suggestions for optimal algorithms , 1995 .

[29]  Farzana Meru,et al.  On the fragmentation criteria of self-gravitating protoplanetary discs , 2010, 1008.0465.

[30]  W. Tscharnuter,et al.  A Fluid Dynamical Treatment of the Common Action of Self-Gravitation, Collisions, and Rotation in Saturn's B-Ring , 1995 .

[31]  Richard I. Klein,et al.  The Jeans Condition: A New Constraint on Spatial Resolution in Simulations of Isothermal Self-Gravitational Hydrodynamics , 1997 .

[32]  James Wadsley,et al.  Fragmentation of Gravitationally Unstable Gaseous Protoplanetary Disks with Radiative Transfer , 2007 .

[33]  M. Bate,et al.  Resolution requirements for smoothed particle hydrodynamics calculations with self-gravity , 1997 .

[34]  A. P. Boss.  Gas giant protoplanets formed by disk instability in binary star systems , 2005 .

[35]  University of Exeter,et al.  On the diffusive propagation of warps in thin accretion discs , 2010, 1002.2973.

[36]  Mark Clampin,et al.  Optical Images of an Exosolar Planet 25 Light-Years from Earth , 2008, Science.