Gravitational Instability Research Articles

Large-scale variability in style and magnitude of footwall rift-related unconformities, Northern Carnarvon Basin, offshore NW Australia

Fault-scarp degradation complexes record rift-related erosion of normal fault scarps. At the scale of an individual fault segment, the magnitude and distribution of erosion are related to the total and along-strike variability in fault throw. However, previous studies on degradation complexes focused on one type of rift-related erosional unconformity, with the factors controlling the development, magnitude, and style of footwall erosion at a far larger scale (i.e., several fault blocks) being poorly constrained. This study uses high-resolution subsurface data to constrain the timing, magnitude, and style of footwall erosion across the NW Shelf of Australia. By doing so, we test conceptual and numerical model predictions for rift-related unconformities development. We show that footwall erosion complexes development is dependent on the occurrence and duration of footwall exposure and that fault throw and footwall uplift control the magnitude of erosion. Furthermore, the style of footwall erosion depends on the driving process(es), with gravitational instabilities resulting in cuspate footwall crest erosional surfaces, whereas peneplain-style surfaces develop in response to wave-driven erosion dictated by the base level. The results of this study can help refine rift-related faults evolution and the effects of sea level changes in tectono-stratigraphic models of rift basins.

The behaviour of eccentric sub-pc massive black hole binaries embedded in massive discs

Using the 3D smoothed-particle hydrodynamics code PHANTOM, we investigated the evolution of the orbital properties of massive black hole binaries embedded in massive discs where gravitational instabilities (GIs) triggered by the disc’s self-gravity are the only significant source of angular momentum transport. In particular, we investigated the evolution of binaries with different initial eccentricities of e0 = 0.05, 0.5, 0.8 and mass ratios of q = 0.1, 0.3, 0.9. Previous studies indicate an equilibrium eccentricity of around 0.6 < e < 0.8, where the binary tidal torques and disc self-gravity torque are in dynamical balance. Our simulations suggest that there might not be a universal value of critical eccentricity. We find that binaries that are initially more eccentric attain higher asymptotic eccentricity than more circular ones do. This implies that there is a range of critical eccentricity values that depend on the initial condition and microphysics (initial eccentricity, mass ratio, cooling) of the system. In particular, we find the width of this range to be narrower for more unequal binaries. We furthermore measured preferential accretion on one of the binary components, only finding more accretion onto the primary for the mass ratio q = 0.3 and eccentricity e = 0.8. We discuss how this might have implications for the amplitude of the gravitational wave background detected by pulsar timing array (PTA) experiments. We finally measure the corresponding value of the viscosity parameter α due to GIs in our simulations and discuss how this depends on the binary properties.

The Stability of a Dense Crust Situated on Small Bodies

Small planetary bodies in the solar system, including Io, Ganymede, and Callisto, may have a crust denser than their underlying mantle. Despite the inherent gravitational instability of such structures, we show that the growth timescale of the Rayleigh–Taylor (RT) instability can be as long as the age of the solar system, owing to the strong temperature dependence of viscosity. Even in cases where the instability timescale is shorter, the instability is confined to a thin layer at the base of the crust, making the foundering of the entire crust improbable in many scenarios. This study delineates the onset and aftermath of the RT instability, applying a quantitative framework to assess the stability of (i) rock-contaminated crust on icy satellites, and (ii) silicate crust floating on top of a subsurface magma ocean on Io. Notably, for Io the RT instability peels off only 10–100 m from the crust’s base, and thermal diffusion rapidly recovers the crustal thickness through solidification of a magma ocean. Despite recurrent delamination of the crustal base, the initial crustal thickness is maintained by thermal diffusion, virtually stabilizing a floating dense crust. Cracking of the crust also is unlikely to result in the foundering of the crust. A dense crust on a small body is therefore difficult to be overturned, suggesting the potential ubiquity of dense surface layers throughout the solar system.

Fragmentation in Gravitationally Unstable Collapsar Disks and Subsolar Neutron Star Mergers

Although stable neutron stars (NSs) can in principle exist down to masses M ns ≈ 0.1 M ⊙, standard models of stellar core-collapse predict a robust lower limit M ns ≳ 1.2 M ⊙, roughly commensurate with the Chandrasekhar mass M Ch of the progenitor’s iron core (electron fraction Y e ≈ 0.5). However, this limit may be circumvented in sufficiently dense neutron-rich environments (Y e < 0.5) for which MCh∝Ye2 is reduced to ≲1 M ⊙. Such physical conditions could arise in the black hole accretion disks formed from the collapse of rapidly rotating stars (“collapsars”), as a result of gravitational instabilities and cooling-induced fragmentation, similar to models for planet formation in protostellar disks. We confirm that the conditions to form subsolar-mass NS (ssNS) may be marginally satisfied in the outer regions of massive neutrino-cooled collapsar disks. If the disk fragments into multiple ssNSs, their subsequent coalescence offers a channel for precipitating subsolar mass LIGO/Virgo gravitational-wave mergers that does not implicate primordial black holes. The model makes several additional predictions: (1) ∼Hz frequency Doppler modulation of the ssNS-merger gravitational-wave signals due to the binary’s orbital motion in the disk; (2) at least one additional gravitational-wave event (coincident within ≲hours), from the coalescence of the ssNS-merger remnant(s) with the central black hole; (3) an associated gamma-ray burst and supernova counterpart, the latter boosted in energy and enriched with r-process elements from the NS merger(s) embedded within the exploding stellar envelope (“kilonovae inside a supernova”).

Orthopyroxene-dominated upper mantle melting built the early crust of the Moon

The paradigm of lunar crust formation has been widely applied to other terrestrial bodies, but the nature of early crust building on the Moon remains enigmatic. Here we report non-Apollo-like highland clasts from the Chang’e-5 mission and find high-alumina melts enclosed in a noritic anorthosite. Geochemistry and phase equilibria modeling suggest that the melt is compositionally parental to lunar magnesian-suite rocks, and was sourced from a plagioclase-bearing, orthopyroxene-dominated upper mantle ( ~ 4.5 kbar and 1225°C). It was formed as a direct consequence of upper mantle melting at the onset of gravitational instability. We propose a continuous early crust formation on the Moon, started from multiple anorthositic cumulate flotations, to upper mantle melting caused by small-scale, in-situ overturn, and eventually ended up by decompression melting of lower mantle cumulates following large-scale, global overturn.

Why does the Milky Way have a metallicity floor?

ABSTRACT The prevalence of light element enhancement in the most metal-poor stars is potentially an indication that the Milky Way has a metallicity floor for star formation around $\sim 10^{-3.5}$ Z$_{\odot }$. We propose that this metallicity floor has its origins in metal-enriched star formation in the minihaloes present during the Galaxy’s initial formation. To arrive at this conclusion, we analyse a cosmological radiation hydrodynamics simulation that follows the concurrent evolution of multiple Population III star-forming minihaloes. The main driver for the central gas within minihaloes is the steady increase in hydrostatic pressure as the haloes grow. We incorporate this insight into a hybrid one-zone model that switches between pressure-confined and modified free-fall modes to evolve the gas density with time according to the ratio of the free-fall and sound-crossing time-scales. This model is able to accurately reproduce the density and chemo-thermal evolution of the gas in each of the simulated minihaloes up to the point of runaway collapse. We then use this model to investigate how the gas responds to the absence of H$_{2}$. Without metals, the central gas becomes increasingly stable against collapse as it grows to the atomic cooling limit. When metals are present in the halo at a level of $\sim 10^{-3.7}$ Z$_{\odot }$, however, the gas is able to achieve gravitational instability while still in the minihalo regime. Thus, we conclude that the Galaxy’s metallicity floor is set by the balance within minihaloes of gas-phase metal cooling and the radiation background associated with its early formation environment.

Tests of subgrid models for star formation using simulations of isolated disc galaxies

ABSTRACT We use smoothed particle hydrodynamics simulations of isolated Milky Way-mass disc galaxies that include cold, interstellar gas to test subgrid prescriptions for star formation (SF). Our fiducial model combines a Schmidt law with a gravitational instability criterion, but we also test density thresholds and temperature ceilings. While SF histories are insensitive to the prescription for SF, the Kennicutt–Schmidt (KS) relations between SF rate and gas surface density can discriminate between models. We show that our fiducial model, with an SF efficiency per free-fall time of 1 per cent, agrees with spatially resolved and azimuthally averaged observed KS relations for neutral, atomic, and molecular gas. Density thresholds do not perform as well. While temperature ceilings selecting cold, molecular gas can match the data for galaxies with solar metallicity, they are unsuitable for very low-metallicity gas and hence for cosmological simulations. We argue that SF criteria should be applied at the resolution limit rather than at a fixed physical scale, which means that we should aim for numerical convergence of observables rather than of the properties of gas labelled as star-forming. Our fiducial model yields good convergence when the mass resolution is varied by nearly 4 orders of magnitude, with the exception of the spatially resolved molecular KS relation at low surface densities. For the gravitational instability criterion, we quantify the impact on the KS relations of gravitational softening, the SF efficiency, and the strength of supernova feedback, as well as of observable parameters such as the inclusion of ionized gas, the averaging scale, and the metallicity.

Unveiling protoplanetary structure equations: Semi-analytical solutions via the homotopy analysis method

This paper revisits the distribution of thermodynamic variables within initial protoplanets formed via gravitational instability (GI) across a broad mass spectrum ranging from 0.3MJ to 10MJ (where 1MJ denotes 1 Jupiter mass, equal to 1.8986×1030 g), using the Homotopy Analysis Method (HAM), a novel approach in this context. Concerning heat transfer within the protoplanets, consideration is given to the convective mode. Our findings reveal a noteworthy alignment between the results obtained via the HAM, utilizing only the first four terms (third approximation), and numerical outcomes. The HAM is found to demonstrate rapid convergence towards the exact solution, showcasing its effectiveness in obtaining such solutions to nonlinear problems. Comparative graphical illustrations of approximate series solutions obtained via HAM and the Adomian decomposition method highlight the superior accuracy and analytical depth of HAM within this framework. This establishes HAM as a powerful and insightful tool for studying astrophysical phenomena. The method offers significant advantages in accuracy and analytical depth over traditional methods, demonstrating its potential for broader applications in astrophysical research and providing robust solutions where conventional approaches may fall short.

How does dark matter stabilize disc galaxies?

ABSTRACT The study presents a theoretical framework for understanding the role of dark matter on the stability of the galactic disc. We model the galaxy as a two-component system consisting of stars and gas in equilibrium with an external dark matter halo. We derive the equations governing the growth of perturbations and obtain a stability criterion that connects the potential of the dark matter halo and the gas fraction with the stability levels of the galaxy. We find that a two-component disc is more susceptible to the growth of gravitational instabilities than individual components, particularly as gas fractions increase. However, the external field, due to the dark matter halo, acts as a stabilizing agent and increases the net stability levels even in the presence of a cold gas component. We apply the stability criterion to models of the Milky Way, low surface brightness galaxies, and baryon-dominated cold rotating disc galaxies observed in the early universe. Our results show that the potential due to the dark matter halo plays a significant role in stabilizing nearby galaxies, such as the Milky Way, and low surface brightness galaxies, which would otherwise be prone to local gravitational instabilities. However, we find that the baryon-dominated cold disc galaxies observed in the early universe remain susceptible to the growth of local gravitational instabilities despite the stabilizing effect of the dark matter halo.

Generation of deposit-derived pyroclastic density currents by repeated crater rim failures at Stromboli Volcano (Italy)

The gravitational instability of hot material deposited during eruptive activity can lead to the formation of glowing avalanches, commonly known as deposit-derived pyroclastic density currents (PDCs). These currents can travel hundreds of metres to several kilometres from the source at exceptionally high temperatures, posing a catastrophic hazard to areas surrounding steep-slope volcanoes. The occurrence of deposit-derived PDCs is often associated with crater rim failure, which can be triggered by various factors such as magma thrust from dike injection, magma fingering, bulging or less commonly, powerful explosions. Here, the in-depth study of data from the multi-parametric monitoring network operating on Stromboli (Italy), including video surveillance, seismicity and ground deformation data, complemented by remote topographic sensing data, has facilitated the understanding of the events leading to the crater rim collapse on 9 October and 4 December 2022. The failures resulted in the remobilisation of 6.4 ± 1.0 × 103 m3 and 88.9 ± 26.7 × 103 m3 of material for the 9 October and the 4 December 2022, respectively, which propagated as PDCs along the NW side of the volcano and reached the sea in a few tens of seconds. These events were characterised by a preparatory phase marked by an increase in magmatic pressure in the preceding weeks, which correlated with an increase in the displacement rate of the volcano’s summit. There was also an escalation in explosive degassing, evidenced by spattering accompanied by seismic tremors in the hours before the collapse.These events have been interpreted as an initial increase in magma vesicularity, followed by the release of gas once percolation threshold was reached. The degassing process induced densification of the magma, resulting in increased thrust on the conduit walls due to increased magmastatic pressure. This phase coincided with crater rim collapse, often followed or accompanied by the onset of lava overflow phases. A mechanism similar to the one proposed may shed light on similar phenomena observed at other volcanoes. The analysis performed in this study highlights the need for a multi-parametric and multi-platform approach to fully understand such complex phenomena. By integrating different data sources, including seismic, deformation and remote sensing data, it is possible to identify the phenomena associated with the different phases leading to crater rim collapse and the subsequent development of deposit-derived PDCs.

Effects of finite Larmor radius corrections and finite electron inertia on transverse Jean’s gravitational instability of radiative plasma with Hall current for star formation in ISM

We learn the consequences of finite ion Larmor radius (FLR) corrections, finite electron inertia, Hall current and radiative heat-loss function on the Jean’s-gravitational instability of an infinite homogeneous, viscous plasma integrating the impacts of finite electrical resistivity, thermal conductivity and permeability. A general dispersion relation is derived by means of the normal mode analysis method with the assistance of linearized perturbation equations of the problem. We find that the original Jeans criterion of gravitational instability is modified into the radiative instability criterion by the inclusion of thermal conductivity and radiative heat-loss function. The finite electrical resistivity eliminates the impact of the magnetic field and the viscosity of the medium eliminates the impact of FLR from radiative instability. The Hall parameter has no impact on the transverse mode of propagation. Numerical calculations demonstrate the stabilizing impact of temperature-dependent heat-loss function, FLR corrections and the destabilizing impact of density-dependent heat-loss function, finite electron inertia on the Jean’s-gravitational instability.

Gravitational instability in magnetized viscoelastic fluid with radiation pressure effect

This study examines the gravitational instability of a finitely conducting magnetized viscoelastic fluid, under the influence of radiation pressure within the framework of a generalized hydrodynamic fluid model. The general dispersion relation, valid for strongly and weakly coupled fluids under kinetic and hydrodynamic limits respectively, is derived using normal mode analysis for a transverse wave propagation mode. The obtained Jeans instability criteria, in terms of critical Jeans wavenumbers under both limits, are modified by the presence of radiative effects, viscoelastic compressional speed, and Alfvén wave velocity. The findings highlight that viscoelastic compressional speed and radiative pressure slow the Jeans instability's growth rate, exerting a stabilizing effect on the initiation of gravitational collapse. The present investigation provides valuable insights into the role of radiative mechanisms in the gravitational collapse within the strongly coupled region of molecular cloud clumps.

C/O Ratios and the Formation of Wide-separation Exoplanets

The gas and solid-state C/O ratios provide context to potentially link the atmospheric composition of planets to that of the natal disk. We provide a synthesis of extant estimates of the gaseous C/O and C/H ratios in planet-forming disks obtained primarily through analysis of Atacama Large Millimeter/submillimeter Array observations. These estimates are compared to atmospheric abundances of wide-separation (>10 au) gas giants. The resolved disk gas C/O ratios, from seven systems, generally exhibit C/O ≥ 1 with subsolar, or depleted, carbon content. In contrast, wide-separation gas giants have atmospheric C/O ratios that cluster near or slightly above the presumed stellar value with a range of elemental C/H. From the existing disk composition, we infer that the solid-state millimeter/centimeter-sized pebbles have a total C/O ratio (solid cores and ices) that is solar (stellar) in content. We explore simple models that reconstruct the exoplanet atmospheric composition from the disk, while accounting for silicate cloud formation in the planet atmosphere. If wide-separation planets formed via the core-accretion mechanism, they must acquire their metals from pebble or planetesimal accretion. Further, the dispersion in giant planet C/H content is best matched by a disk composition with modest and variable factors of carbon depletion. An origin of the wide-separation gas giants via gravitational instability cannot be ruled out, as stellar C/O ratios should natively form in this scenario. However, the variation in planet metallicity with a stellar C/O ratio potentially presents challenges to these models.

Formation of Giant Planets by Gas Disk Gravitational Instability on Wide Orbits around Protostars with Varied Masses. II. Quadrupled Spatial Resolution and Beta Cooling

Exoplanet demographics are sufficiently advanced to provide important constraints on theories of planet formation. While core and pebble accretion are preferred for rocky and icy planets, there appears to be a need for gas disk gravitational instability (GDGI) to play a role in the formation of M-dwarf gas giants and those orbiting at large distances. Here we present GDGI models that go beyond those presented by Boss (2011) dealing with the formation of wide-orbit gas giants. The new models use quadrupled spatial resolution, in both the radial and azimuthal directions, to reduce the effects of finite spatial resolution. The new models also employ the β cooling approximation, instead of the diffusion approximation used by Boss (2011), in order to push the models further in time. As in Boss (2011), the central protostars have masses of 0.1, 0.5, 1.0, 1.5, or 2.0 M ⊙, surrounded by disks with masses ranging from 0.019 M ⊙ to 0.21 M ⊙. For each case, two models are computed, one with an initial minimum Toomre Q stability value ranging from 1.1 to 1.7, and one with a higher initial disk temperature, resulting in the initial minimum Q ranging from 2.2 to 3.4. These new models continue to show that GDGI can explain the formation of gas giants at distances of ∼30 to ∼50 au on eccentric orbits (e less than ∼0.2), though the number formed drops to 0 as the protostar mass decreases to 0.1 M ⊙.

Cloud Dissipation and Disk Wind in the Late Phase of Star Formation

We perform a long-term simulation of star and disk formation using three-dimensional nonideal magnetohydrodynamics. The simulation starts from a prestellar cloud and proceeds through the long-term evolution of the circumstellar disk until ∼1.5 × 105 yr after protostar formation. The disk has size ≲50 au and little substructure in the main accretion phase because of the action of magnetic braking and the magnetically driven outflow to remove angular momentum. The main accretion phase ends when the outflow breaks out of the cloud, causing the envelope mass to decrease rapidly. The outflow subsequently weakens as the mass accretion rate also weakens. While the envelope-to-disk accretion continues, the disk grows gradually and develops transient spiral structures, due to gravitational instability. When the envelope-to-disk accretion ends, the disk becomes stable and reaches a size ≳300 au. In addition, about 30% of the initial cloud mass has been ejected by the outflow. A significant finding of this work is that after the envelope dissipates a revitalization of the wind occurs, and there is mass ejection from the disk surface that lasts until the end of the simulation. This mass ejection (or disk wind) is generated because the magnetic pressure significantly dominates both the ram pressure and thermal pressure above and below the disk at this stage. Using the angular momentum flux and mass-loss rate estimated from the disk wind, the disk dissipation timescale is estimated to be ∼106 yr.

A 3D view on the local gravitational instability of cold gas discs in star-forming galaxies at 0 ≲ z ≲ 5

Local gravitational instability (LGI) is considered crucial for regulating star formation and gas turbulence in galaxy discs, especially at high redshift. Instability criteria usually assume infinitesimally thin discs or rely on approximations to include the stabilising effect of the gas disc thickness. We test a new 3D instability criterion for rotating gas discs that are vertically stratified in an external potential. This criterion reads Q3D &lt; 1, where Q3D is the 3D analogue of the Toomre parameter Q. The advantage of Q3D is that it allows us to study LGI in and above the galaxy midplane in a rigorous and self-consistent way. We apply the criterion to a sample of 44 star-forming galaxies at 0 ≲ z ≲ 5 hosting rotating discs of cold gas. The sample is representative of galaxies on the main sequence at z ≈ 0 and includes massive star-forming and starburst galaxies at 1 ≲ z ≲ 5. For each galaxy, we first apply the Toomre criterion for infinitesimally thin discs, finding ten unstable systems. We then obtain maps of Q3D from a 3D model of the gas disc derived in the combined potential of dark matter, stars and the gas itself. According to the 3D criterion, two galaxies with Q &lt; 1 show no evidence of instability and the unstable regions that are 20% smaller than those where Q &lt; 1. No unstable disc is found at 0 ≲ z ≲ 1, while ≈60% of the systems at 2 ≲ z ≲ 5 are locally unstable. In these latter, a relatively small fraction of the total gas (≈30%) is potentially affected by the instability. Our results disfavour LGI as the main regulator of star formation and turbulence in moderately star-forming galaxies in the present-day Universe. LGI likely becomes important at high redshift, but the input by other mechanisms seems required in a significant portion of the disc. We also estimate the expected mass of clumps in the unstable regions, offering testable predictions for observations.

Weighing protoplanetary discs with kinematics: Physical model, method, and benchmark

The mass of protoplanetary discs determines the amount of material available for planet formation, the level of coupling between gas and dust, and possibly also sets gravitational instabilities. Measuring the mass of a disc is challenging, as it is not possible to directly detect H$_ $, and CO-based estimates remain poorly constrained. An alternative method has recently been proposed that does not rely on tracer-to-H$_ $ ratios. It allows dynamical measurement of the disc mass together with the star mass and the disc critical radius by looking at deviations from Keplerian rotation induced by the self-gravity of the disc. So far, this method has been used to weigh three protoplanetary discs: Elias 2-27, IM Lup, and GM Aurigae. odot $ with the phantom code and post-processed them with radiative transfer ( mcfost ) to obtain synthetic observations. We find that dynamical weighing allows us to retrieve the expected disc masses as long as the disc-to-star mass ratio is larger than $M_d/M_ We estimate an uncertainty for the disc mass measurement of $ 25&lt;!PCT!&gt;$.

Dynamical modelling and the origin of gas turbulence in z ∼ 4.5 galaxies

Context. In recent years, a growing number of regularly rotating galaxy discs have been found at z ≥ 4. Such systems provide us with the unique opportunity to study the properties of dark matter (DM) halos at these early epochs, the turbulence within the interstellar medium and the evolution of scaling relations. Aims. Here, we investigate the dynamics of four gas discs in galaxies at z ∼ 4.5 observed with ALMA in the [CII] 158 μm fine-structure line. We aim to derive the structural properties of the gas, stars and DM halos of the galaxies and to study the mechanisms driving the turbulence in high-z discs. Methods. We decomposed the rotation curves into baryonic and DM components within the extent of the [CII] discs, that is, 3 to 5 kpc. Furthermore, we used the gas velocity dispersion profiles as a diagnostic tool in investigating the mechanisms driving the turbulence in the discs. Results. We obtain total stellar, gas and DM masses in the ranges of log(M/M⊙) = 10.3 − 11.0, 9.8 − 11.3, and 11.2 − 13.3, respectively. We find dynamical evidence in all four galaxies for the presence of compact stellar components conceivably, stellar bulges. The turbulence present in the galaxies appears to be primarily driven by stellar feedback, negating the necessity for large-scale gravitational instabilities. Finally, we investigate the position of our galaxies in the context of local scaling relations, in particular the stellar-to-halo mass and Tully–Fisher analogue relations.

M17 MIR: A Massive Star Is Forming via Episodic Mass Accretion

We analyzed the Atacama Large Millimeter/submillimeter Array (ALMA) band 6 data for the outbursting massive protostar M17 MIR. The ALMA CO J = 2–1 data reveal a collimated and bipolar north–south outflow from M17 MIR. The blueshifted outflow exhibits four CO knots (N1 to N4) along the outflow axis, while the redshifted outflow appears as a single knot (S1). The extremely high velocity (EHV) emissions of N1 and S1 are jetlike and contain subknots along the outflow axis. Assuming the nearest EHV subknots trace the ejecta from the accretion outbursts in the past decades, a tangential ejection velocity of ∼421 km s−1 is derived for M17 MIR. Assuming the same velocity, the dynamical times of the multiple ejecta, traced by the four blueshifted CO knots, range from 20 to 364 yr. The four blueshifted CO knots imply four clustered accretion outbursts with a duration of tens of years in the past few hundred years. The intervals between the four clustered accretion outbursts are also about tens of years. These properties of the four clustered accretion outbursts are in line with the disk gravitational instability and fragmentation model. The episodic accretion history of M17 MIR traced by episodic outflow suggests that a massive star can form from a lower-mass protostar via frequent episodic accretion events triggered by disk gravitational instability and fragmentation. The first detection of the knotty outflow from an outbursting massive protostar suggests that mass ejections accompanied with accretion events could serve as an effective diagnostic tool for the episodic accretion histories of massive protostars.

The backreaction of stellar wobbling on accretion discs of massive protostars

black In recent years, it has been demonstrated that massive stars see their infant circumstellar medium shaped into a large irradiated, gravitationally unstable accretion disc during their early formation phase. Such discs constitute the gas reservoir from which nascent high-mass stars gain a substantial fraction of their mass by episodic accretion of dense gaseous circumstellar clumps, simultaneously undergoing accretion-driven bursts and producing close-orbit spectroscopic companions of the young high-mass stellar object. We aim to evaluate the effects of stellar motion caused by the disc non-axisymmetric gravitational field on the disc evolution and its spatial morphology. In particular, we analyse the disc's propensity to gravitational instability and fragmentation and the disc's appearance in synthetic millimetre band images pertinent to the ALMA facility. We employed three-dimensional radiation-hydrodynamical simulations of the surroundings of a young massive star in the non-inertial spherical coordinate system, adopting the highest spatial resolution to date and including the indirect star-disc gravitational potential caused by the asymmetries in the circumstellar disc. The resulting disc configurations were post-processed with the radiation transfer tool RADMC-3D and CASA software to obtain synthetic images of the disc. We confirm that the early evolution of the accretion disc is notably different when stellar wobbling is taken into account. The redistribution of angular momentum in the system makes the disc smaller and rounder, reduces the number of circumstellar gaseous clumps formed via disc gravitational fragmentation, and prevents the ejection of gaseous clumps from the disc. The synthetic predictive images at millimetre wavelengths of the accretion disc that includes stellar wobbling are in better agreement with the observations of the surroundings of massive young stellar objects, namely black AFGL 4176 mm1, G17.64+0.16, and G353.273 than our simulations of numerical hydrodynamics that omit this physical mechanism. Our work confirms that stellar wobbling is an essential ingredient to account for in numerical simulations of accretion discs of massive protostars.