Abstract:Context. There is growing evidence that massive stars grow by disk accretion in a similar way to their low-mass counterparts. Early in evolution, these disks can achieve masses that are comparable to the current stellar mass, and therefore the forming disks are highly susceptible to gravitational fragmentation. Aims. We investigate the formation and early evolution of an accretion disk around a forming massive protostar, focussing on its fragmentation physics. To this end, we follow the collapse of a molecular cloud of gas and dust, the formation of a massive protostar, the formation of its circumstellar disk, and the formation and evolution of the disk fragments. Methods. We used a grid-based, self-gravity radiation hydrodynamics code including a sub-grid module for stellar evolution and dust evolution. We purposely do not use a sub-grid module for fragmentation such as sink particles to allow for all paths of fragment formation and destruction, but instead we keep the spatial grid resolution high enough to properly resolve the physical length scales of the problem, namely the pressure scale height and Jeans length of the disk. Simulations are performed on a grid in spherical coordinates with a logarithmic spacing of the grid cells in the radial direction and a cosine distribution of the grid cells in the polar direction, focusing the spatial resolution on the disk midplane. As a consequence, roughly 25% of the total number of grid cells, corresponding to ~26 million grid cells, are used to model the disk physics. These constitute the highest resolution simulations performed up to now on disk fragmentation around a forming massive star with the physics considered here. For a better understanding of the effects of spatial resolution and to compare our high-resolution results with previous lower resolution studies in the literature, we perform the same simulation at five different resolutions, each run differing in resolution from its predecessor by a factor of two. Results. The cloud collapses and a massive (proto)star is formed in its center surrounded by a fragmenting Keplerian-like accretion disk with spiral arms. The fragments have masses of ~1 M ⊙ , and their continuous interactions with the disk, spiral arms, and other fragments result in eccentric orbits. Fragments form hydrostatic cores surrounded by secondary disks with spiral arms that also produce new fragments. We identified several mechanisms of fragment formation, interaction, and destruction. Central temperatures of the fragments can reach the hydrogen dissociation limit, form second Larson cores, and evolve into companion stars. Based on this, we study the multiplicity predicted by the simulations and find approximately six companions at different distances from the primary: from possible spectroscopic multiples, to companions at distances between 1000 and 2000 au.

A new binary formation mechanism

Formation and evolution of protostellar accretion discs. I. Angular-momentum budget, gravitational self-regulation, and numerical convergence

The Formation Mass of a Binary System via Fragmentation of a Rotating Parent Core with Increasing Total Mass

A new formation scenario of a counter-rotating circumstellar disk: spiral-arm accretion from a circumbinary disk in a triple protostar system

Fragmentation of dense rotation-dominated structures fed by collapsing gravo-magneto-sheetlets and origin of misaligned 100 au-scale binaries and multiple systems

Turbulence-induced disc formation in strongly magnetised cloud cores

Disc formation in turbulent massive cores: Circumventing the magnetic braking catastrophe

The role of disc self-gravity in the formation of protostars and protostellar discs

Formation of massive close binaries I. Detached evolution

Gravitationally Coupled Scale-Free Discs

Structure of a Protobinary System: An Asymmetric Circumbinary Disk and Spiral Arms

Formation of Large Circumstellar Discs in Multi-scale, ideal-MHD Simulations of Magnetically Critical Pre-stellar Cores

Modeling disk fragmentation and multiplicity in massive star formation

Coevolution of Binaries and Gaseous Discs

Pathways for producing binary black holes with large misaligned spins in the isolated formation channel

Formation of discs around super-massive black hole binaries

Fragmentation of Dense Rotation-Dominated Structures Fed by Collapsing Gravomagneto-Sheetlets and Origin of Misaligned 100 au-Scale Binaries and Multiple Systems

Formation of multiple low mass stars, brown dwarfs and planemos via gravitational collapse

Time-dependent models of the structure and evolution of self-gravitating protoplanetary discs

Modeling the secular evolution of embedded protoplanetary discs