Observations of young star-forming regions suggest that star clusters are born completely mass segregated. These initial conditions are, however, gradually lost as the star cluster evolves dynamically. For star clusters with single stars only and a canonical initial mass function, it has been suggested that traces of these initial conditions vanish at a time $\tau_\mathrm{v}$ between 3 and $3.5\,t_\mathrm{rh}$ (initial half-mass relaxation times). Since a significant fraction of stars are observed in binary systems and it is widely accepted that most stars are born in binary systems, we aim to investigate what role a primordial binary population (even up to $100\,\%$ binaries) plays in the loss of primordial mass segregation of young star clusters. We used numerical $N$-body models similar in size to the Orion Nebula Cluster (ONC) -- a representative of young open clusters -- integrated over several relaxation times to draw conclusions on the evolution of its mass segregation. We also compared our models to the observed ONC. We found that $\tau_\mathrm{v}$ depends on the binary star fraction and the distribution of initial binary parameters that include a semi-major axis, eccentricity, and mass ratio. For instance, in the models with $50\,\%$ binaries, we find $\tau_\mathrm{v} = (2.7 \pm 0.8)\,t_\mathrm{rh}$, while for $100\,\%$ binary fraction, we find a lower value $\tau_\mathrm{v} = (2.1 \pm 0.6)\,t_\mathrm{rh}$. We also conclude that the initially completely mass segregated clusters, even with binaries, are more compatible with the present-day ONC than the non-segregated ones.

Observational results of young star-forming regions suggest that star clusters are completely mass segregated at birth. As a star cluster evolves dynamically, these initial conditions are gradually lost. For star clusters with single stars only and a canonical IMF, it has been suggested that traces of these initial conditions vanish at $\tau_{\rm v}$ between 3 and 3.5 half-mass relaxation times. By the means of numerical models, here we investigate the role of the primordial binary population on the loss of primordial mass segregation. We found that $\tau_{\rm v}$ does not seem to depend on the binary star distribution, yielding $3 < \tau_{\rm v} / t_{\rm rh} < 3.5$. We also conclude that the completely mass segregated clusters, even with binaries, are more compatible with the present-day ONC than the non-segregated ones.

ALMA observations of the Serpens South star-forming region suggest that stellar protoclusters may be completely mass segregated at birth. Independent observations also suggest that embedded clusters form segregated by mass. As the primordial mass segregation seems to be lost over time, we aim to study on which timescale an initially perfectly mass-segregated star cluster becomes indistinguishable from an initially not mass-segregated cluster. As an example, the Orion Nebula Cluster (ONC) is also discussed. We used $N$-body simulations of star clusters with various masses and two different degrees of primordial mass segregation. We analysed their energy redistribution through two-body relaxation to quantify the time when the models agree in terms of mass segregation, which sets in only dynamically in the models that are primordially not mass segregated. A comprehensive cross-matched catalogue combining optical, infrared, and X-ray surveys of ONC members was also compiled and made available. The models evolve to a similar radial distribution of high-mass stars after the core collapse (about half a median two-body relaxation time, $t_\mathrm{rh}$) and become observationally indistinguishable from the point of view of mass segregation at time $\tau_\mathrm{v} \approx 3.3\,t_\mathrm{rh}$. In the case of the ONC, using the distribution of high-mass stars, we may not rule out either evolutionary scenario (regardless of whether they are initially mass segregated). When we account for extinction and elongation of the ONC, as reported elsewhere, an initially perfectly mass-segregated state seems to be more consistent with the observed cluster.

We present an account of the state of the art in the fields explored by the research community invested in 'Modeling and Observing DEnse STellar systems'. For this purpose, we take as a basis the activities of the MODEST-17 conference, which was held at Charles University, Prague, in September 2017. Reviewed topics include recent advances in fundamental stellar dynamics, numerical methods for the solution of the gravitational N-body problem, formation and evolution of young and old star clusters and galactic nuclei, their elusive stellar populations, planetary systems, and exotic compact objects, with timely attention to black holes of different classes of mass and their role as sources of gravitational waves. Such a breadth of topics reflects the growing role played by collisional stellar dynamics in numerous areas of modern astrophysics. Indeed, in the next decade, many revolutionary instruments will enable the derivation of positions and velocities of individual stars in the Milky Way and its satellites and will detect signals from a range of astrophysical sources in different portions of the electromagnetic and gravitational spectrum, with an unprecedented sensitivity. On the one hand, this wealth of data will allow us to address a number of long-standing open questions in star cluster studies; on the other hand, many unexpected properties of these systems will come to light, stimulating further progress of our understanding of their formation and evolution.

Core collapse is a prominent evolutionary stage of self-gravitating systems. In an idealised collisionless approximation, the region around the cluster core evolves in a self-similar way prior to the core collapse. Thus, its radial density profile outside the core can be described by a power law, $\rho \propto r^{-\alpha}$. We aim to find the characteristics of core collapse in $N$-body models. In such systems, a complete collapse is prevented by transferring the binding energy of the cluster to binary stars. The contraction is, therefore, more difficult to identify. We developed a method that identifies the core collapse in $N$-body models of star clusters based on the assumption of their homologous evolution. We analysed different models (equal- and multi-mass), most of which exhibit patterns of homologous evolution, yet with significantly different values of $\alpha$: the equal-mass models have $\alpha \approx 2.3$, which agrees with theoretical expectations, the multi-mass models have $\alpha \approx 1.5$ (yet with larger uncertainty). Furthermore, most models usually show sequences of separated homologous collapses with similar properties. Finally, we investigated a correlation between the time of core collapse and the time of formation of the first hard binary star. The binding energy of such a binary usually depends on the depth of the collapse in which it forms, for example from $100\,kT$ to $10^4\,kT$ in the smallest equal-mass to the largest multi-mass model, respectively. However, not all major hardenings of binaries happened during the core collapse. In the multi-mass models, we see large transfers of binding energy of $\sim 10^4\,kT$ to binaries that occur on the crossing timescale and outside of the periods of the homologous collapses.

Globular clusters (GCs) are the ideal environment for the formation of neutron stars (NSs) and millisecond pulsars (MSPs). NSs origin and evolution provide a useful information on stellar dynamics and evolution in star clusters, and are among the most interesting astrophysical objects, being precursors of several high-energy phenomena such as gravitational waves and gamma-ray bursts. Due to a large velocity kick that they receive at birth, most of the NSs escape the local field, affecting the evolution and dynamics of their parent cluster. In this paper, we study the origin and dynamical evolution of NSs within GCs with different initial masses, metallicities and primordial binary fractions. We find that the radial profile of NSs is shaped by the BH content of the cluster, which partially quenches the NS segregation until most of the BHs are ejected from the system. Independently on the cluster mass and initial configuration, the NSs map the average stellar population, as their average radial distance is $\approx 60-80\%$ of the cluster half-mass radius. Finally, by assuming a recycling fraction of $f_\mathrm{rec}=0.1$ and an average MSP gamma-ray emission of $L_\gamma=2\times 10^{33}$ erg s$^{-1}$, we show that the typical gamma-ray emission from our GCs agrees with observations and supports the MSP origin of the gamma-ray excess signal observed by the Fermi-LAT telescope in the Galactic Centre.

Recent research has been constraining the retention fraction of black holes (BHs) in globular clusters by comparing the degree of mass segregation with $N$-body simulations. They are consistent with an upper limit of the retention fraction being $50\,\%$ or less. In this work, we focus on direct simulations of the dynamics of BHs in star clusters. We aim to constrain the effective distribution of natal kicks that BHs receive during supernova (SN) explosions and to estimate the BH retention fraction. We used the collisional $N$-body code nbody6 to measure the retention fraction of BHs for a given set of parameters, which are: the initial mass of a star cluster, the initial half-mass radius, and $\sigma_\mathrm{BH}$, which sets the effective Maxwellian BH velocity kick distribution. We compare these direct $N$-body models with our analytic estimates and newest observational constraints. The numerical simulations show that for the one-dimensional (1D) velocity kick dispersion $\sigma_\mathrm{BH} < 50\,\mathrm{km\,s^{-1}}$, clusters with radii of 2 pc and that are initially more massive than $5 \times 10^3\,M_\odot$ retain more than $20\,\%$ of BHs within their half-mass radii. Our simple analytic model yields a number of retained BHs that is in good agreement with the $N$-body models. Furthermore, the analytic estimates show that ultra-compact dwarf galaxies (UCDs) should have retained more than $80\,\%$ of their BHs for $\sigma_\mathrm{BH} \leq 190\,\mathrm{km\,s^{-1}}$. Although our models do not contain primordial binaries, in the most compact clusters with $10^3$ stars, we have found evidence of delayed SN explosions producing a surplus of BHs compared to the IMF due to dynamically formed binary stars. These cases do not occur in the more populous or expanded clusters.