No Arabic abstract
AIMS: The aim of this work is to understand whether there is a difference in the dispersion of discs around stars in high-density young stellar clusters like the Orion Nebula Cluster (ONC) according to the mass of the star. METHODS: Two types of simulations were combined -- N-body simulations of the dynamics of the stars in the ONC and mass loss results from simulations of star-disc encounters, where the disc mass loss of all stars is determined as a function of time. RESULTS: We find that in the Trapezium, the discs around high-mass stars are dispersed much more quickly and to a larger degree by their gravitational interaction than for intermediate-mass stars. This is consistent with the very recent observations of IC 348, where a higher disc frequency was found around solar mass stars than for more massive stars, suggesting that this might be a general trend in large young stellar clusters.
We use a suite of SPH simulations to investigate the susceptibility of protoplanetary discs to the effects of self-gravity as a function of star-disc properties. We also include passive irradiation from the host star using different models for the stellar luminosities. The critical disc-to-star mass ratio for axisymmetry (for which we produce criteria) increases significantly for low-mass stars. This could have important consequences for increasing the potential mass reservoir in a proto Trappist-1 system, since even the efficient Ormel et al. (2017) formation model will be influenced by processes like external photoevaporation, which can rapidly and dramatically deplete the dust reservoir. The aforementioned scaling of the critical $M_d/M_*$ for axisymmetry occurs in part because the Toomre $Q$ parameter has a linear dependence on surface density (which promotes instability) and only an $M_*^{1/2}$ dependence on shear (which reduces instability), but also occurs because, for a given $M_d/M_*$, the thermal evolution depends on the host star mass. The early phase stellar irradiation of the disc (for which the luminosity is much higher than at the zero age main sequence, particularly at low stellar masses) can also play a key role in significantly reducing the role of self-gravity, meaning that even Solar mass stars could support axisymmetric discs a factor two higher in mass than usually considered possible. We apply our criteria to the DSHARP discs with spirals, finding that self-gravity can explain the observed spirals so long as the discs are optically thick to the host star irradiation.
In young star clusters, the density can be high enough and the velocity dispersion low enough for stars to collide and merge with a significant probability. This has been suggested as a possible way to build up the high-mass portion of the stellar mass function and as a mechanism leading to the formation of one or two very massive stars (M > 150 Msun) through a collisional runaway. I quickly review the standard theory of stellar collisions, covering both the stellar dynamics of dense clusters and the hydrodynamics of encounters between stars. The conditions for collisions to take place at a significant rate are relatively well understood for idealised spherical cluster models without initial mass segregation, devoid of gas and composed of main-sequence (MS) stars. In this simplified situation, 2-body relaxation drives core collapse through mass segregation and a collisional phase ensues if the core collapse time is shorter than the MS lifetime of the most massive stars initially present. The outcome of this phase is still highly uncertain. A more realistic situation is that of a cluster still containing large amounts of interstellar gas from which stars are accreting. As stellar masses increase, the central regions of the cluster contracts. This little-explored mechanism can potentially lead to very high stellar densities but it is likely that, except for very rich systems, the contraction is halted by few-body interactions before collisions set in. A complete picture, combining both scenarios, will need to address many uncertainties, including the role of cluster sub-structure, the dynamical effect of interstellar gas, non-MS stars and the structure and evolution of merged stars.
Star clusters appear to be the ideal environment for the assembly of neutron star-neutron star (NS-NS) and black hole-neutron star (BH-NS) binaries. These binaries are among the most interesting astrophysical objects, being potential sources of gravitational waves (GWs) and gamma-ray bursts. We use for the first time high-precision N-body simulations of young massive and open clusters to study the origin and dynamical evolution of NSs, within clusters with different initial masses, metallicities, primordial binary fractions, and prescriptions for the compact object natal kicks at birth. We find that the radial profile of NSs is shaped by the BH content of the cluster, which partially quenches the NS segregation due to the BH-burning process. This leaves most of the NSs out of the densest cluster regions, where NS-NS and BH-NS binaries could potentially form. Due to a large velocity kick that they receive at birth, most of the NSs escape the host clusters, with the bulk of their retained population made up of NSs of $sim 1.3$ M$_odot$ coming from the electron-capture supernova process. The details of the primordial binary fraction and pairing can smear out this trend. Finally, we find that a subset of our models produce NS-NS mergers, leading to a rate of $sim 0.01$--$0.1$ Gpc$^{-3}$ yr$^{-1}$ in the local Universe, and compute an upper limit of $sim 3times 10^{-2}$--$3times 10^{-3}$ Gpc$^{-3}$ yr$^{-1}$ for the BH-NS merger rate. Our estimates are several orders of magnitude smaller than the current empirical merger rate from LIGO/Virgo, in agreement with the recent rate estimates for old globular clusters.
In the last decade we have come to realize that the traditional classification of stellar clusters into open and globular clusters cannot be easily extended beyond the realm of the Milky Way, and that even for our Galaxy it is not fully valid. The main failure of the traditional classification is the existence of Massive Young Clusters (MYCs), which are massive like Globular Clusters (GCs) but also young like open clusters. We describe here the mass and age distributions of clusters in general with an emphasis on MYCs. We also discuss the issue of what constitutes a cluster and try to establish a general classification scheme.
We perform $N$-body simulations of the early phases of open cluster evolution including a large population of planetesimals, initially arranged in Kuiper-belt like discs around each star. Using a new, 4th-order and time-reversible $N$-body code on Graphics Processing Units (GPUs), we evolve the whole system under the stellar gravity, i.e. treating planetesimals as test particles, and consider two types of initial cluster models, similar to IC348 and the Hyades, respectively. In both cases, planetesimals can be dynamically excited, transferred between stars or liberated to become free-floating (such as A/2017 U1 or Oumuamua) during the early cluster evolution. We find that planetesimals captured from another star are not necessarily dynamically distinct from those native to a star. After an encounter both native and captured planetesimals can exhibit aligned periastrons, qualitatively similar to that seen in the Solar system and commonly thought to be the signature of Planet 9. We discuss the implications of our results for both our Solar system and exoplanetary systems.