More and more observations indicate that young star clusters could retain imprints of their formation process. In particular, the degree of substructuring and rotation are possibly the direct result of the collapse of the parent molecular cloud from which these systems form. Such properties can, in principle, be washed-out, but they are also expected to have an impact on the relaxation of these systems. We ran and analyzed a set of ten hydrodynamical simulations of the formation of embedded star clusters through the collapse of turbulent massive molecular clouds. We systematically studied the fractality of our star clusters, showing that they are all extremely substructured (fractal dimension $D=1.0-1.8$). We also found that fractality is slowly reduced, with time, on small scales, while it persists on large scales on longer timescales. Signatures of rotation are found in different simulations at every time of the evolution, even for slightly supervirial substructures, proving that the parent molecular gas transfers part of its angular momentum to the new stellar systems.

Young star clusters are the most common birth-place of massive stars and are dynamically active environments. Here, we study the formation of black holes (BHs) and binary black holes (BBHs) in young star clusters, by means of 6000 N-body simulations coupled with binary population synthesis. We probe three different stellar metallicities (Z=0.02, 0.002 and 0.0002) and two initial density regimes (density at the half-mass radius $\rho_{\rm h}\ge{}3.4\times10^4$ and $\ge{1.5\times10^2}$ M$_\odot$ pc$^{-3}$ in dense and loose star clusters, respectively). Metal-poor clusters tend to form more massive BHs than metal-rich ones. We find $\sim{}6$, $\sim{}2$, and $<1$% of BHs with mass $m_{\rm BH}>60$ M$_\odot$ at Z=0.0002, 0.002 and 0.02, respectively. In metal-poor clusters, we form intermediate-mass BHs with mass up to $\sim{}320$ M$_\odot$. BBH mergers born via dynamical exchanges (exchanged BBHs) can be more massive than BBH mergers formed from binary evolution: the former (latter) reach total mass up to $\sim{}140$ M$_\odot$ ($\sim{}80$ M$_\odot$). The most massive BBH merger in our simulations has primary mass $\sim{}88$ M$_\odot$, inside the pair-instability mass gap, and a mass ratio of $\sim{}0.5$. Only BBHs born in young star clusters from metal-poor progenitors can match the masses of GW170729, the most massive event in O1 and O2, and those of GW190412, the first unequal-mass merger. We estimate a local BBH merger rate density $\sim{}110$ and $\sim{}55$ Gpc$^{-3}$ yr$^{-1}$, if we assume that all stars form in loose and dense star clusters, respectively.

In weakly collisional stellar systems such as some globular clusters, partial energy equipartition and mass segregation are expected to develop as a result of the cumulative effect of stellar encounters even in systems initially characterized by star-mass independent density and energy distributions. In parallel, numerical simulations have demonstrated that radially-biased pressure anisotropy slowly builds up in realistic models of globular clusters from initial isotropic conditions, leading to anisotropy profiles that, to some extent, mimic those resulting from incomplete violent relaxation known to be relevant to elliptical galaxies. In this paper we consider a set of realistic simulations realized by means of Monte Carlo methods and analyze them by means of self-consistent two-component models. For the purpose, we refer to an underlying distribution function, originally conceived to describe elliptical galaxies, that has been recently truncated and adapted to the context of globular clusters. The two components are meant to represent light stars (combining all main sequence stars) and heavy stars (giants, dark remnants, and binaries). We show that this conceptually simple family of two-component truncated models provides a reasonable description of simulated density, velocity dispersion, and anisotropy profiles, especially for the most relaxed systems, with the capability to express quantitatively the attained levels of energy equipartition and mass segregation. In contrast, two-component isotropic models based on the King distribution function do not offer a comparably satisfactory representation of the simulated globular clusters. With this work we provide a new reliable diagnostic tool applicable to nonrotating globular clusters that are characterized by significant gradients in the local value of the mass-to-light ratio, beyond the commonly used one-component dynamical models.

As a result of the slow action of two-body encounters, globular clusters develop mass segregation and attain a condition of only partial energy equipartition even in their central, most relaxed regions. Realistic numerical simulations show that, during the process, a radially-biased anisotropy profile slowly builds up, mimicking that resulting from incomplete violent relaxation. Commonly used dynamical models, such as the one-component King models, cannot describe these properties. Here we show that simple two-component models based on a distribution function originally conceived to describe elliptical galaxies, recently truncated and adapted to the context of globular clusters, can describe in detail what is observed in complex and realistic numerical simulations.