We conduct a census of the high-mass protostellar population of the $\sim70,000\:M_\odot$ Infrared Dark Cloud (IRDC) G028.37+00.07, identifying 35 sources based on their $70\:\mu$m emission, as reported in the {\it Herschel} Hi-GAL catalog of Molinari et al. (2016). We perform aperture photometry to construct spectral energy distributions (SEDs), which are then fit with the massive protostar models of Zhang & Tan (2018). We find that the sources span a range of isotropic luminosities from $\sim$20 to 4,500$\:L_\odot$. The most luminous sources are predicted to have current protostellar masses of $m_{*}\sim10\:M_\odot$ forming from cores of mass $M_{c}\sim40$ to $400\:M_\odot$. The least luminous sources in our sample are predicted to be protostars with masses as low as $\sim 0.5\:M_\odot$ forming from cores with $M_{c}\sim10\:M_\odot$, which are the minimum values explored in the protostellar model grid. The detected protostellar population has a total estimated protostellar mass of $M_{*}\sim 100\:M_\odot$. Allowing for completeness corrections, which are constrained by comparison with an ALMA study in part of the cloud, we estimate a star formation efficiency per free-fall time of $\sim3\%$ in the IRDC. Finally, analyzing the spatial distribution of the sources, we find relatively low degrees of central concentration of the protostars. The protostars, including the most massive ones, do not appear to be especially centrally concentrated in the protocluster as defined by the IRDC boundary.

We use Gaia DR2 to hunt for runaway and walkaway stars from the Orion Nebula Cluster (ONC). We search a large area of the sky extending 45$^\circ$ in radius around the ONC and out to a distance of 1 kpc to find sources that overlapped in angular position with the cluster in the last ~10 Myr, also accounting for the ONC's own proper motion. We construct a main sample of ~17,000 runaway/walkaway candidates selected by this 2D traceback condition. However, most of these are expected to be contaminants, e.g., caused by Galactic streaming motions of stars at different distances. We thus examine six further tests or "flags" to help identify real runaways, namely: (1) possessing young stellar object (YSO) colors and magnitudes based on Gaia optical photometry; (2) having IR excess consistent with YSOs based on 2MASS and WISE photometry; (3) having a high degree of optical variability; (4) having closest approach distances to the ONC core well constrained within the half-mass radius of the cluster; (5) having direction of ejection from the ONC that avoids the main contamination zone due to Galactic streaming; and (6) having a required radial velocity (RV) for 3D overlap with the ONC that is of reasonable magnitude (or, for the small fraction, 7%, of candidate stars with Gaia measured RVs, satisfying 3D traceback). Thirteen sources, not previously noted as Orion members, pass all of these tests, while another twelve are similarly promising, except that they are located in the zone most contaminated by Galactic streaming. Among these 25 ejection candidates, there are ten with measured RVs that pass the most restrictive 3D traceback condition. We present full lists of runaway/walkaway candidate stars, estimate the high-velocity population ejected from the ONC and discuss its implications for star cluster formation theories via comparison with numerical simulations.

We investigate the dynamical evolution of star clusters during their formation, assuming that they are born from a turbulent starless clump of a given mass that is embedded within a parent self-gravitating molecular cloud characterized by a particular mass surface density. In contrast to the standard practice of most $N$-body studies, we do not assume that all stars are formed at once. Rather, we explore the effects of different star formation rates on the global structure and evolution of young embedded star clusters, also considering various primordial binary fractions and mass segregation levels. Our fiducial clumps studied in this paper have initial masses of $M_{\rm cl} = 3000\,M_\odot$, are embedded in ambient cloud environments of $\Sigma_{\rm cloud} = 0.1$ and 1 g cm$^{-2}$, and gradually form stars with an overall efficiency of 50% until the gas is exhausted. We investigate star formation efficiencies per free-fall time in the range $\epsilon_{\rm ff}=0.01$ to 1, and also compare to the instantaneous case ($\epsilon_{\rm ff}=\infty$) of Paper I. We show that most of the interesting dynamical processes that determine the future of the cluster, happen during the early formation phase. In particular, the ejected stellar population is sensitive to the duration of star cluster formation: for example, clusters with longer formation times produce more runaway stars, since these clusters remain in a dense state for longer, thus favouring occurrence of dynamical ejections. We also show that the presence of radial age gradients in star clusters depends sensitively on the star formation efficiency per free fall time, with observed values being matched best by our slowest forming clusters with $\epsilon_{\rm ff}\lesssim0.03$.

We explore scenarios for the dynamical ejection of stars BN and x from source I in the Kleinmann-Low nebula of the Orion Nebula Cluster (ONC), which is important for being the closest region of massive star formation. This ejection would cause source I to become a close binary or a merger product of two stars. We thus consider binary-binary encounters as the mechanism to produce this event. By running a large suite of $N$-body simulations, we find that it is nearly impossible to match the observations when using the commonly adopted masses for the participants, especially a source I mass of $7\:{\rm{M}}_\odot$. The only way to recreate the event is if source I is more massive, i.e., $\sim20\:{\rm{M}}_\odot$. However, even in this case, the likelihood of reproducing the observed system is low. We discuss the implications of these results for understanding this important star-forming region.

We investigate the response of initially substructured, young, embedded star clusters to instantaneous gas expulsion of their natal gas. We introduce primordial substructure to the stars and the gas by simplistically modelling the star formation process so as to obtain a variety of substructure distributed within our modelled star forming regions. We show that, by measuring the virial ratio of the stars alone (disregarding the gas completely), we can estimate how much mass a star cluster will retain after gas expulsion to within 10% accuracy, no matter how complex the background structure of the gas is, and we present a simple analytical recipe describing this behaviour. We show that the evolution of the star cluster while still embedded in the natal gas, and the behavior of the gas before being expelled, are crucial processes that affect the timescale on which the cluster can evolve into a virialized spherical system. Embedded star clusters that have high levels of substructure are subvirial for longer times, enabling them to survive gas expulsion better than a virialized and spherical system. By using a more realistic treatment for the background gas than our previous studies, we find it very difficult to destroy the young clusters with instantaneous gas expulsion. We conclude that gas removal may not be the main culprit for the dissolution of young star clusters.

We investigate the evolution of mass segregation in initially sub-structured young embedded star clusters with two different background potentials mimicking the gas. Our clusters are initially in virial or sub-virial global states and have different initial distributions for the most massive stars: randomly placed, initially mass segregated or even inverse segregation. By means of N-body simulation we follow their evolution for 5 Myr. We measure the mass segregation using the minimum spanning tree method Lambda_MSR and an equivalent restricted method. Despite this variety of different initial conditions, we find that our stellar distributions almost always settle very fast into a mass segregated and more spherical configuration, suggesting that once we see a spherical or nearly spherical embedded star cluster, we can be sure it is mass segregated no matter what the real initial conditions were. We, furthermore, report under which circumstances this process can be more rapid or delayed, respectively.

We investigate the formation and early evolution of star clusters assuming that they form from a turbulent starless clump of given mass bounded inside a parent self-gravitating molecular cloud characterized by a particular mass surface density. As a first step we assume instantaneous star cluster formation and gas expulsion. We draw our initial conditions from observed properties of starless clumps. We follow the early evolution of the clusters up to 20 Myr, investigating effects of different star formation efficiencies, primordial binary fractions and eccentricities and primordial mass segregation levels. We investigate clumps with initial masses of $M_{\rm cl}=3000\:{\rm M}_\odot$ embedded in ambient cloud environments with mass surface densities, $\Sigma_{\rm cloud}=0.1$ and $1\:{\rm g\:cm^{-2}}$. We show that these models of fast star cluster formation result, in the fiducial case, in clusters that expand rapidly, even considering only the bound members. Clusters formed from higher $\Sigma_{\rm cloud}$ environments tend to expand more quickly, so are soon larger than clusters born from lower $\Sigma_{\rm cloud}$ conditions. To form a young cluster of a given age, stellar mass and mass surface density, these models need to assume a parent molecular clump that is many times denser, which is unrealistic compared to observed systems. We also show that in these models the initial binary properties are only slightly modified by interactions, meaning that binary properties, e.g., at 20 Myr, are very similar to those at birth. With this study we set up the basis of future work where we will investigate more realistic models of star formation compared to this instantaneous, baseline case.

We attempt to find a progenitor for the ultra-faint object Segue 1 under the assumption that it formed as a dark matter free star cluster in the past. We look for orbits, using the elongation of Segue 1 on the sky as a tracer of its path. Those orbits are followed backwards in time to find the starting points of our N-body simulations. The successful orbit, with which we can reproduce Segue 1 has a proper motion of mu_alpha = -0.19 mas/yr and mu_delta = -1.9 mas/yr, placing Segue 1 near its apo-galacticon today. Our best fitting model has an initial mass of 6224 Msun and an initial scale-length of 5.75 pc.

We examine the effects of gas-expulsion on initially substructured distributions of stars. We perform N-body simulations of the evolution of these distributions in a static background potential to mimic the gas. We remove the static potential instantaneously to model gas-expulsion. We find that the exact dynamical state of the cluster plays a very strong role in affecting a cluster's survival, especially at early times: they may be entirely destroyed or only weakly affected. We show that knowing both detailed dynamics and relative star-gas distributions can provide a good estimate of the post-gas expulsion state of the cluster, but even knowing these is not an absolute way of determining the survival or otherwise of the cluster.

Hercules is a dwarf spheroidal satellite of the Milky Way, found at a distance of about 138 kpc, and showing evidence of tidal disruption. It is very elongated and exhibits a velocity gradient of 16 +/- 3 km/s/kpc. Using this data a possible orbit of Hercules has previously been deduced in the literature. In this study we make use of a novel approach to find a best fit model that follows the published orbit. Instead of using trial and error, we use a systematic approach in order to find a model that fits multiple observables simultaneously. As such, we investigate a much wider parameter range of initial conditions and ensure we have found the best match possible. Using a dark matter free progenitor that undergoes tidal disruption, our best-fit model can simultaneously match the observed luminosity, central surface brightness, effective radius, velocity dispersion, and velocity gradient of Hercules. However, we find it is impossible to reproduce the observed elongation and the position angle of Hercules at the same time in our models. This failure persists even when we vary the duration of the simulation significantly, and consider a more cuspy density distribution for the progenitor. We discuss how this suggests that the published orbit of Hercules is very likely to be incorrect.

Recent deep photometry of the dwarf spheroidal Ursa Major II's morphology, and spectroscopy of individual stars, have provided a number of new constraints on its properties. With a velocity dispersion $\sim$6 km s$^{-1}$, and under the assumption that the galaxy is virialised, the mass-to-light ratio is found to be approaching $\sim$2000 - apparently heavily dark matter dominated. Using N-Body simulations, we demonstrate that the observed luminosity, ellipticity, irregular morphology, velocity gradient, and the velocity dispersion can be well reproduced through processes associated with tidal mass loss, and in the absence of dark matter. These results highlight the considerable uncertainty that exists in measurements of the dark matter content of Ursa Major II. The dynamics of the inner tidal tails, and tidal stream, causes the observed velocity dispersion of stars to be boosted to values of $>$5 km s$^{-1}$ ($>$20 km s$^{-1}$ at times). This effect is responsible for raising the velocity dispersion of our model to the observed values in UMaII. We test an iterative rejection technique for removing unbound stars from samples of UMaII stars whose positions on the sky, and line-of-sight velocities, are provided. We find this technique is very effective at providing an accurate bound mass from this information, and only fails when the galaxy has a bound mass less than 10$%$ of its initial mass. However when $<2%$ mass remains bound, mass overestimation by $>$3 orders of magnitude are seen. Additionally we find that mass measurements are sensitive to measurement uncertainty in line-of-sight velocities. Measurement uncertainties of 1-4 km s$^{-1}$ result in mass overestimates by a factor of $\sim$1.3-5.7.

The dynamics of globular cluster systems (GCSs) around galaxies are often used to assess the total enclosed mass, and even to constrain the dark matter distribution. The globular cluster system of a galaxy is typically assumed to be in dynamical equilibrium within the potential of the host galaxy. However cluster galaxies are subjected to a rapidly evolving and, at times, violently destructive tidal field. We investigate the impact of the harassment on the dynamics of GCs surrounding early type cluster dwarfs, using numerical simulations. We find that the dynamical behaviour of the GCS is strongly influenced by the fraction of bound dark matter f_{DM} remaining in the galaxy. Only when f_{DM} falls to ~15%, do stars and GCs begin to be stripped. Still the observed GC velocity dispersion can be used to measure the true enclosed mass to within a factor of 2, even when f_{DM} falls as low as ~3%. This is possible partly because unbound GCs quickly separate from the galaxy body. However even the distribution of {\it{bound}} GCs may spatially expand by a factor of 2-3. Once f_{DM} falls into the <3% regime, the galaxy is close to complete disruption, and GCS dynamics can no longer be used to reliably estimate the enclosed mass. In this regime, the remaining bound GCS may spatially expand by a factor of 4 to 8. It may be possible to test if a galaxy is in this regime by measuring the dynamics of the stellar disk. We demonstrate that if a stellar disk is rotationally supported, it is likely that a galaxy has sufficient dark matter, that the dynamics of the GCS can be used to reliably estimate the enclosed mass.