Young Moving Groups (YMGs) are close (<100pc), coherent collections of young (<100Myr) stars that appear to have formed in the same star-forming molecular cloud. As such we would expect their individual initial mass functions (IMFs) to be similar to other star-forming regions, and by extension the Galactic field. Their close proximity to the Sun and their young ages means that YMGs are promising locations to search for young forming exoplanets. However, due to their low numbers of stars, stochastic sampling of the IMF means their stellar populations could vary significantly. We determine the range of planet-hosting stars (spectral types A, G and M) possible from sampling the IMF multiple times, and find that some YMGs appear deficient in M-dwarfs. We then use these data to show that the expected probability of detecting terrestrial magma ocean planets is highly dependent on the exact numbers of stars produced through stochastic sampling of the IMF.

Galactic Wolf-Rayet (WR) star membership of star forming regions can be used to constrain the formation environments of massive stars. Here, we utilise $\textit{Gaia}$ DR2 parallaxes and proper motions to reconsider WR star membership of clusters and associations in the Galactic disk, supplemented by recent near-IR studies of young massive clusters. We find that only 18$-$36% of 553 WR stars external to the Galactic Centre region are located in clusters, OB associations or obscured star-forming regions, such that at least 64% of the known disk WR population are isolated, in contrast with only 13% of O stars from the Galactic O star Catalogue. The fraction located in clusters, OB associations or star-forming regions rises to 25$-$41% from a global census of 663 WR stars including the Galactic Centre region. We use simulations to explore the formation processes of isolated WR stars. Neither runaways, nor low mass clusters, are numerous enough to account for the low cluster membership fraction. Rapid cluster dissolution is excluded as mass segregation ensures WR stars remain in dense, well populated environments. Only low density environments consistently produce WR stars that appeared to be isolated during the WR phase. We therefore conclude that a significant fraction of WR progenitors originate in low density association-like surroundings which expand over time. We provide distance estimates to clusters and associations host to WR stars, and estimate cluster ages from isochrone fitting.

Theory predicts that we should find fast, ejected (runaway) stars of all masses around dense, young star-forming regions. $N$-body simulations show that the number and distribution of these ejected stars could be used to constrain the initial spatial and kinematic substructure of the regions. We search for runaway and slower walkaway stars within 100 pc of the Orion Nebula Cluster (ONC) using $Gaia$ DR2 astrometry and photometry. We compare our findings to predictions for the number and velocity distributions of runaway stars from simulations that we run for 4 Myr with initial conditions tailored to the ONC. In $Gaia$ DR2, we find 31 runaway and 54 walkaway candidates based on proper motion, but not all of these are viable candidates in three dimensions. About 40 per cent are missing radial velocities, but we can trace back 9 3D-runaways and 24 3D-walkaways to the ONC, all of which are low/intermediate-mass (<8 M$_{\odot}$). Our simulations show that the number of runaways within 100 pc decreases the older a region is (as they quickly travel beyond this boundary), whereas the number of walkaways increases up to 3 Myr. We find fewer walkaways in $Gaia$ DR2 than the maximum suggested from our simulations, which may be due to observational incompleteness. However, the number of $Gaia$ DR2 runaways agrees with the number from our simulations during an age of $\sim$1.3-2.4 Myr, allowing us to confirm existing age estimates for the ONC (and potentially other star-forming regions) using runaway stars.

The Q-parameter is used extensively to quantify the spatial distributions of stars and gas in star-forming regions as well as older clusters and associations. It quantifies the amount of structure using the ratio of the average length of a minimum spanning tree, mbar, to the average length within the complete graph, sbar. The interpretation of the Q-parameter often relies on comparing observed values of Q, mbar and sbar to idealised synthetic geometries, where there is little or no match between the observed star-forming regions and the synthetic regions. We measure Q, mbar, and sbar over 10 Myr in N-body simulations which are compared to IC 348, NGC 1333, and the ONC. For each star-forming region we set up simulations that approximate their initial conditions for a combination of different virial rations and fractal dimensions. We find that dynamical evolution of idealised fractal geometries can account for the observed Q, mbar, and sbar values in nearby star-forming regions. In general, an initially fractal star-forming region will tend to evolve to become more smooth and centrally concentrated. However, we show that initial conditions, as well as where the edge of the region is defined, can cause significant differences in the path that a star-forming region takes across the mbar-sbar plot as it evolves. We caution that the observed Q-parameter should not be directly compared to idealised geometries. Instead, it should be used to determine the degree to which a star-forming region is either spatially substructured or smooth and centrally concentrated.

Observations of pre-/proto-stellar cores in young star-forming regions show them to be mass segregated, i.e. the most massive cores are centrally concentrated, whereas pre-main sequence stars in the same star-forming regions (and older regions) are not. We test whether this apparent contradiction can be explained by the massive cores fragmenting into stars of much lower mass, thereby washing out any signature of mass segregation in pre-main sequence stars. Whilst our fragmentation model can reproduce the stellar initial mass function, we find that the resultant distribution of pre-main sequence stars is mass segregated to an even higher degree than that of the cores, because massive cores still produce massive stars if the number of fragments is reasonably low (between one and five). We therefore suggest that the reason cores are observed to be mass segregated and stars are not is likely due to dynamical evolution of the stars, which can move significant distances in star-forming regions after their formation.

We examine the mass-dependence of the velocity dispersion of stars in the young cluster NGC 6530 to better understand how it formed. Using a large sample of members we find that the proper motion velocity dispersion increases with stellar mass. While this trend is the opposite to that predicted if the cluster were developing energy equipartition, it is in agreement with recent N-body simulations that find such a trend develops because of the Spitzer instability. In these simulations the massive stars sink to the centre of the cluster and form a self-gravitating system with a higher velocity dispersion. If the cluster has formed by the cool collapse of an initially substructured distribution then this occurs within 1-2 Myr, in agreement with our observations of NGC 6530. We therefore conclude that NGC 6530 formed from much more extended initial conditions and has since collapsed to form the cluster we see now. This cluster formation model is inconsistent with the idea that all stars form in dense, compact clusters and provides the first dynamical evidence that star clusters can form by hierarchical mergers between subclusters.

We use $N$-body simulations to probe the early phases of the dynamical evolution of star-forming regions and focus on mass and velocity distributions of unbound stars. In this parameter space study, we vary the initial virial ratio and degree of spatial and kinematic substructure and analyse the fraction of stars that become unbound in two different mass classes (above and below 8 M$_{\odot}$). We find that the fraction of unbound stars differs depending on the initial conditions. After 10 Myr, in initially highly subvirial, substructured simulations, the high-mass and lower-mass unbound fractions are similar at $\sim$23 per cent. In initially virialised, substructured simulations, we find only $\sim$16 per cent of all high-mass stars are unbound, whereas $\sim$37 per cent of all lower-mass stars are. The velocity distributions of unbound stars only show differences for extremely different initial conditions. The distributions are dominated by large numbers of lower-mass stars becoming unbound just above the escape velocity of $\sim$3 km s$^{-1}$ with unbound high-mass stars moving faster on average than lower-mass unbound stars. We see no high-mass runaway stars (velocity > 30 km s$^{-1}$) from any of our initial conditions and only an occasional lower-mass runaway star from initially subvirial/substructured simulations. In our simulations, we find a small number of lower-mass walkaway stars (with velocity 5-30 km s$^{-1}$) from all of our initial conditions. These walkaway stars should be observable around many nearby star-forming regions with Gaia.

Understanding how the birthplace of stars affects planet-forming discs is important for a comprehensive theory of planet formation. Most stars are born in dense star-forming regions where the external influence of other stars, particularly the most massive stars, will affect the survival and enrichment of their planet-forming discs. Simulations suggest that stellar dynamics play a central role in regulating how external feedback affects discs, but comparing models to observations requires an estimate of the initial stellar density in star-forming regions. Structural analyses constrain the amount of dynamical evolution a star-forming region has experienced; regions that maintain substructure and do not show mass segregation are likely dynamically young, and therefore close to their birth density. In this paper, we present a structural analysis of two clusters in the Carina Nebula, Tr14 and Tr16. We show that neither cluster shows evidence for mass segregation or a centrally concentrated morphology, suggesting that both regions are dynamically young. This allows us to compare to simulations from Nicholson et al. (2019) who predict disc survival rates in star-forming regions of different initial densities. The surviving disc fractions in Tr14 and Tr16 are consistent with their predictions (both are $\sim 10$%), supporting a growing body of evidence that the star-forming environment plays an important role in the survival and enrichment of protoplanetary discs.

The combination of precise radial velocities from multi-object spectroscopy and highly accurate proper motions from Gaia DR2 opens up the possibility for detailed 3D kinematic studies of young star forming regions and clusters. Here, we perform such an analysis by combining Gaia-ESO Survey spectroscopy with Gaia astrometry for ~900 members of the Lagoon Nebula cluster, NGC 6530. We measure the 3D velocity dispersion of the region to be $5.35^{+0.39}_{-0.34}$~km~s$^{-1}$, which is large enough to suggest the region is gravitationally unbound. The velocity ellipsoid is anisotropic, implying that the region is not sufficiently dynamically evolved to achieve isotropy, though the central part of NGC 6530 does exhibit velocity isotropy that suggests sufficient mixing has occurred in this denser part. We find strong evidence that the stellar population is expanding, though this is preferentially occurring in the declination direction and there is very little evidence for expansion in the right ascension direction. This argues against a simple radial expansion pattern, as predicted by models of residual gas expulsion. We discuss these findings in the context of cluster formation, evolution and disruption theories.

Habitable zones are regions around stars where large bodies of liquid water can be sustained on a planet or satellite. As many stars form in binary systems with non-zero eccentricity, the habitable zones around the component stars of the binary can overlap and be enlarged when the two stars are at periastron (and less often when the stars are at apastron). We perform N-body simulations of the evolution of dense star-forming regions and show that binary systems where the component stars originally have distinct habitable zones can undergo interactions that push the stars closer together, causing the habitable zones to merge and become enlarged. Occasionally, overlapping habitable zones can occur if the component stars move further apart, but the binary becomes more eccentric. Enlargement of habitable zones happens to 1-2 binaries from an average initial total of 352 in each simulated star-forming region, and demonstrates that dense star-forming regions are not always hostile environments for planet formation and evolution.

We analyse N-body simulations of star-forming regions to investigate the effects of external far and extreme ultra-violet photoevaporation from massive stars on protoplanetary discs. By varying the initial conditions of simulated star-forming regions, such as the spatial distribution, net bulk motion (virial ratio), and density, we investigate which parameters most affect the rate at which discs are dispersed due to external photoevaporation. We find that disc dispersal due to external photoevaporation is faster in highly substructured star-forming regions than in smooth and centrally concentrated regions. Sub-virial star-forming regions undergoing collapse also show higher rates of disc dispersal than regions that are in virial equilibrium or are expanding. In moderately dense ($\sim$100 M$_{\odot}$ pc$^{-3}$) regions, half of all protoplanetary discs with radii $\geq$ 100 AU are photoevaporated within 1 Myr, three times faster than is currently suggested by observational studies. Discs in lower-density star-forming regions ($\sim$10 M$_{\odot}$ pc$^{-3}$) survive for longer, but half are still dispersed on short timescales ($\sim$2 Myr). This demonstrates that the initial conditions of the star forming regions will greatly impact the evolution and lifetime of protoplanetary discs. These results also imply that either gas giant planet formation is extremely rapid and occurs before the gas component of discs is evaporated, or gas giants only form in low-density star-forming regions where no massive stars are present to photoevaporate gas from protoplanetary discs.

Here we confirm Melnick 34, an X-ray bright star in the 30 Doradus region of the Large Magellanic Cloud, as an SB2 binary comprising WN5h+WN5h components. We present orbital solutions using 26 epochs of VLT/UVES spectra and 22 epochs of archival Gemini/GMOS spectra. Radial-velocity monitoring and automated template fitting methods both reveal a similar high eccentricity system with a mass ratio close to unity, and an orbital period in agreement with the 155.1 +/- 1 day X-ray light curve period previously derived by Pollock et al. Our favoured solution derived an eccentricity of 0.68 +/- 0.02 and mass ratio of 0.92 +/- 0.07, giving minimum masses of Ma_sin^{3}(i) = 65 +/- 7 Msun and Mb_sin^{3}(i) = 60 +/- 7 Msun. Spectral modelling using WN5h templates with CMFGEN reveals temperatures of T ~53 kK for each component and luminosities of log(La/Lsun) = 6.43 +/- 0.08 and log(Lb/Lsun) = 6.37 +/- 0.08, from which BONNSAI evolutionary modelling gives masses of Ma = 139 (+21,-18) Msun and Mb = 127 (+17,-17) Msun and ages of ~0.6 Myrs. Spectroscopic and dynamic masses would agree if Mk34 has an inclination of i ~50{\deg}, making Mk34 the most massive binary known and an excellent candidate for investigating the properties of colliding wind binaries. Within 2-3 Myrs, both components of Mk34 are expected to evolve to stellar mass black holes which, assuming the binary system survives, would make Mk34 a potential binary black hole merger progenitor and gravitational wave source.

The Friends of Friends algorithm identifies groups of objects with similar spatial and kinematic properties, and has recently been used extensively to quantify the distributions of gas and stars in young star-forming regions. We apply the Friends of Friends algorithm to $N$-body simulations of the dynamical evolution of subvirial (collapsing) and supervirial (expanding) star-forming regions. We find that the algorithm picks out a wide range of groups (1 -- 25) for statistically identical initial conditions, and cannot distinguish between subvirial and supervirial regions in that we obtain similar mode and median values for the number of groups it identifies. We find no correlation between the number of groups identified initially and either the initial or subsequent spatial and kinematic tracers of the regions' evolution, such as the amount of spatial substructure, dynamical mass segregation, or velocity dispersion. We therefore urge caution in using the Friends of Friends algorithm to quantify the initial conditions of star formation.

We quantify the spatial distributions of dense cores in three spatially distinct areas of the Orion B star-forming region. For L1622, NGC2068/NGC2071 and NGC2023/NGC2024 we measure the amount of spatial substructure using the $\mathcal{Q}$-parameter and find all three regions to be spatially substructured ($\mathcal{Q} < 0.8$). We quantify the amount of mass segregation using $\Lambda_{\rm MSR}$ and find that the most massive cores are mildly mass segregated in NGC2068/NGC2071 ($\Lambda_{\rm MSR} \sim 2$), and very mass segregated in NGC2023/NGC2024 ($\Lambda_{\rm MSR} = 28^{+13}_{-10}$ for the four most massive cores). Whereas the most massive cores in L1622 are not in areas of relatively high surface density, or deeper gravitational potentials, the massive cores in NGC2068/NGC2071 and NGC2023/NGC2024 are significantly so. Given the low density (10 cores pc$^{-2}$) and spatial substructure of cores in Orion B, the mass segregation cannot be dynamical. Our results are also inconsistent with simulations in which the most massive stars form via competitive accretion, and instead hint that magnetic fields may be important in influencing the primordial spatial distributions of gas and stars in star-forming regions.

We quantify the structure of a very large number of Galactic open clusters and look for evidence of mass segregation for the most massive stars in the clusters. We characterise the structure and mass segregation ratios of 1276 clusters in the Milky Way Stellar Cluster (MWSC) catalogue containing each at least 40 stars and that are located at a distance of up to $\approx 2$ kpc from the Sun. We use an approach based on the calculation of the minimum spanning tree of the clusters, and for each one of them, we calculate the structure parameter \Q\ and the mass segregation ratio $\Lambda_{\rm MSR}$. Our findings indicate that most clusters possess a \Q\ parameter that falls in the range 0.7-0.8 and are thus neither strongly concentrated nor do they show significant substructure. Only 27\% can be considered centrally concentrated with \Q\ values $> 0.8$. Of the 1276 clusters, only 14\% show indication of significant mass segregation ($\Lambda_{\rm MSR} > 1.5$). Furthermore, no correlation is found between the structure of the clusters or the degree of mass segregation with their position in the Galaxy. A comparison of the measured \Q\ values for the young open clusters in the MWSC to N-body numerical simulations that follow the evolution of the \Q\ parameter over the first 10 Myrs of the clusters life suggests that the young clusters found in the MWSC catalogue initially possessed local mean volume densities of $\rho_{*} \approx 10-100$ M$_{\odot}$ pc$^{-3}$.

The presence of an unseen `Planet 9' on the outskirts of the Solar system has been invoked to explain the unexpected clustering of the orbits of several Edgeworth--Kuiper Belt Objects. We use $N$-body simulations to investigate the probability that Planet 9 was a free-floating planet (FFLOP) that was captured by the Sun in its birth star-formation environment. We find that only 1 - 6 per cent of FFLOPs are ensnared by stars, even with the most optimal initial conditions for capture in star-forming regions (one FFLOP per star, and highly correlated stellar velocities to facilitate capture). Depending on the initial conditions of the star-forming regions, only 5 - 10 of 10000 planets are captured onto orbits that lie within the constraints for Planet 9. When we apply an additional environmental constraint for Solar system formation - namely the injection of short-lived radioisotopes into the Sun's protoplanetary disc from supernovae - we find that the probability for the capture of Planet 9 to be almost zero.

We examine the level of substructure and mass segregation in the massive, young cluster Westerlund 1. We find that it is relatively smooth, with little or no mass segregation, but with the massive stars in regions of significantly higher than average surface density. While an expanding or bouncing-back scenario for the evolution of Westerlund 1 cannot be ruled out, we argue that the most natural model to explain these observations is one in which Westerlund 1 formed with no primordial mass segregation and at a similar or larger size than we now observe.

Approximately 10 per cent of star clusters are found in pairs, known as binary clusters. We propose a mechanism for binary cluster formation; we use N-body simulations to show that velocity substructure in a single (even fairly smooth) region can cause binary clusters to form. This process is highly stochastic and it is not obvious from a region's initial conditions whether a binary will form and, if it does, which stars will end up in which cluster. We find the probability that a region will divide is mainly determined by its virial ratio, and a virial ratio above 'equilibrium' is generally necessary for binary formation. We also find that the mass ratio of the two clusters is strongly influenced by the initial degree of spatial substructure in the region.

We analyse N-body and Smoothed Particle Hydrodynamic (SPH) simulations of young star-forming regions to search for differences in the spatial distributions of massive stars compared to lower-mass stars. The competitive accretion theory of massive star formation posits that the most massive stars should sit in deeper potential wells than lower-mass stars. This may be observable in the relative surface density or spatial concentration of the most massive stars compared to other, lower-mass stars. Massive stars in cool--collapse N-body models do end up in significantly deeper potentials, and are mass segregated. However, in models of warm (expanding) star-forming regions, whilst the massive stars do come to be in deeper potentials than average stars, they are not mass segregated. In the purely hydrodynamical SPH simulations, the massive stars do come to reside in deeper potentials, which is due to their runaway growth. However, when photoionisation and stellar winds are implemented in the simulations, these feedback mechanisms regulate the mass of the stars and disrupt the inflow of gas into the clouds' potential wells. This generally makes the potential wells shallower than in the control runs, and prevents the massive stars from occupying deeper potentials. This in turn results in the most massive stars having a very similar spatial concentration and surface density distribution to lower-mass stars. Whilst massive stars do form via competitive accretion in our simulations, this rarely translates to a different spatial distribution and so any lack of primordial mass segregation in an observed star-forming region does not preclude competitive accretion as a viable formation mechanism for massive stars.

We present analyses of the spatial distributions of stars in the young (1 - 3 Myr) star-forming regions IC348 and NGC1333 in the Perseus Giant Molecular Cloud. We quantify the spatial structure using the $\mathcal{Q}$-parameter and find that both IC348 and NGC1333 are smooth and centrally concentrated with $\mathcal{Q}$-parameters of 0.98 and 0.89 respectively. Neither region exhibits mass segregation ($\Lambda_{\rm MSR} = 1.1^{+0.2}_{-0.3}$ for IC348 and $\Lambda_{\rm MSR} = 1.2^{+0.4}_{-0.3}$ for NGC1333, where $\Lambda_{\rm MSR} \sim 1$ corresponds to no mass segregation), nor do the most massive stars reside in areas of enhanced stellar surface density compared to the average surface density, according to the $\Sigma_{\rm LDR}$ method. We then constrain the dynamical histories and hence initial conditions of both regions by comparing the observed values to $N$-body simulations at appropriate ages. Stars in both regions likely formed with sub-virial velocities which contributed to merging of substructure and the formation of smooth clusters. The initial stellar densities were no higher than $\rho \sim 100 - 500$M$_\odot$pc$^{-3}$ for IC348 and $\rho \sim 500 - 2000$M$_\odot$pc$^{-3}$ for NGC1333. These initial densities, in particular that of NGC1333, are high enough to facilitate dynamical interactions which would likely affect $\sim$10 per cent of protoplanetary discs and binary stars.

The young (~2 Myr) cluster Chamaeleon I is one of the closest laboratories to study the early stages of star cluster dynamics in a low-density environment. We studied its structural and kinematical properties combining parameters from the high-resolution spectroscopic survey Gaia-ESO with data from the literature. Our main result is the evidence of a large discrepancy between the velocity dispersion (sigma = 1.14 \pm 0.35 km s^{-1}) of the stellar population and the dispersion of the pre-stellar cores (~0.3 km s^{-1}) derived from submillimeter observations. The origin of this discrepancy, which has been observed in other young star clusters is not clear. It may be due to either the effect of the magnetic field on the protostars and the filaments, or to the dynamical evolution of stars driven by two-body interactions. Furthermore, the analysis of the kinematic properties of the stellar population put in evidence a significant velocity shift (~1 km s^{-1}) between the two sub-clusters located around the North and South main clouds. This result further supports a scenario, where clusters form from the evolution of multiple substructures rather than from a monolithic collapse. Using three independent spectroscopic indicators (the gravity indicator $\gamma$, the equivalent width of the Li line, and the H_alpha 10\% width), we performed a new membership selection. We found six new cluster members located in the outer region of the cluster. Starting from the positions and masses of the cluster members, we derived the level of substructure Q, the surface density \Sigma and the level of mass segregation $\Lambda_{MSR}$ of the cluster. The comparison between these structural properties and the results of N-body simulations suggests that the cluster formed in a low density environment, in virial equilibrium or supervirial, and highly substructured.

The presence and abundance of short lived radioisotopes (SLRs) $^{26}$Al and $^{60}$Fe in chondritic meteorites implies that the Sun formed in the vicinity of one or more massive stars that exploded as supernovae (SNe). Massive stars are more likely to form in massive star clusters ($>$1000 M$_{\odot}$) than lower mass clusters. However, photoevaporation of protoplanetary discs from massive stars and dynamical interactions with passing stars can inhibit planet formation in clusters with radii of $\sim$1 pc. We investigate whether low-mass (50 - 200 M$_{\odot}$) star clusters containing one or two massive stars are a more likely avenue for early Solar system enrichment as they are more dynamically quiescent. We analyse $N$-body simulations of the evolution of these low-mass clusters and find that a similar fraction of stars experience supernova enrichment than in high mass clusters, despite their lower densities. This is due to two-body relaxation, which causes a significant expansion before the first supernova even in clusters with relatively low (100 stars pc$^{-3}$) initial densities. However, because of the high number of low mass clusters containing one or two massive stars, the absolute number of enriched stars is the same, if not higher than for more populous clusters. Our results show that direct enrichment of protoplanetary discs from supernovae occurs as frequently in low mass clusters containing one or two massive stars (>20 M$_{\odot}$) as in more populous star clusters (1000 M$_\odot$). This relaxes the constraints on the direct enrichment scenario and therefore the birth environment of the Solar System.

In this chapter I review the effects of supernovae explosions on the dynamical evolution of (1) binary stars and (2) star clusters. (1) Supernovae in binaries can drastically alter the orbit of the system, sometimes disrupting it entirely, and are thought to be partially responsible for `runaway' massive stars - stars in the Galaxy with large peculiar velocities. The ejection of the lower-mass secondary component of a binary occurs often in the event of the more massive primary star exploding as a supernova. The orbital properties of binaries that contain massive stars mean that the observed velocities of runaway stars (10s - 100s km s$^{-1}$) are consistent with this scenario. (2) Star formation is an inherently inefficient process, and much of the potential in young star clusters remains in the form of gas. Supernovae can in principle expel this gas, which would drastically alter the dynamics of the cluster by unbinding the stars from the potential. However, recent numerical simulations, and observational evidence that gas-free clusters are observed to be bound, suggest that the effects of supernova explosions on the dynamics of star clusters are likely to be minimal.

We present 3.7 arcsec (~0.05 pc) resolution 3.2 mm dust continuum observations from the IRAM PdBI, with the aim of studying the structure and fragmentation of the filamentary Infrared Dark Cloud G035.39-00.33. The continuum emission is segmented into a series of 13 quasi-regularly spaced (~0.18pc) cores, following the major axis of the IRDC. We compare the spatial distribution of the cores with that predicted by theoretical work describing the fragmentation of hydrodynamic fluid cylinders, finding a significant (a factor of ~8) discrepancy between the two. Our observations are consistent with the picture emerging from kinematic studies of molecular clouds suggesting that the cores are harboured within a complex network of independent sub-filaments. This result emphasizes the importance of considering the underlying physical structure, and potentially, dynamically important magnetic fields, in any fragmentation analysis. The identified cores exhibit a range in (peak) beam-averaged column density ($3.6{\rm x}10^{23}{\rm cm}^{-2}<N_{H,c}<8.0{\rm x}10^{23}{\rm cm}^{-2}$), mass ($8.1M_{\odot}<M_{c}<26.1M_{\odot}$), and number density ($6.1{\rm x}10^{5}{\rm cm}^{-3}<n_{H, c, eq}<14.7{\rm x}10^{5}{\rm cm}^{-3}$). Two of these cores, dark in the mid-infrared, centrally-concentrated, monolithic (with no traceable substructure at our PdBI resolution), and with estimated masses of the order ~20-25$M_{\odot}$, are good candidates for the progenitors of intermediate-to-high-mass stars. Virial parameters span a range $0.2<\alpha_{\rm vir}<1.3$. Without additional support, possibly from dynamically important magnetic fields with strengths of the order 230$\mu$G<B<670$\mu$G, the cores are susceptible to gravitational collapse. These results may imply a multilayered fragmentation process, which incorporates the formation of sub-filaments, embedded cores, and the possibility of further fragmentation.

Heating by short-lived radioisotopes (SLRs) such as aluminum-26 and iron-60 fundamentally shaped the thermal history and interior structure of Solar System planetesimals during the early stages of planetary formation. The subsequent thermo-mechanical evolution, such as internal differentiation or rapid volatile degassing, yields important implications for the final structure, composition and evolution of terrestrial planets. SLR-driven heating in the Solar System is sensitive to the absolute abundance and homogeneity of SLRs within the protoplanetary disk present during the condensation of the first solids. In order to explain the diverse compositions found for extrasolar planets, it is important to understand the distribution of SLRs in active planet formation regions (star clusters) during their first few Myr of evolution. By constraining the range of possible effects, we show how the imprint of SLRs can be extrapolated to exoplanetary systems and derive statistical predictions for the distribution of aluminum-26 and iron-60 based on N-body simulations of typical to large clusters (1000-10000 stars) with a range of initial conditions. We quantify the pollution of protoplanetary disks by supernova ejecta and show that the likelihood of enrichment levels similar to or higher than the Solar System can vary considerably, depending on the cluster morphology. Furthermore, many enriched systems show an excess in radiogenic heating compared to Solar System levels, which implies that the formation and evolution of planetesimals could vary significantly depending on the birth environment of their host stars.

Mass segregation in star clusters is often thought to indicate the onset of energy equipartition, where the most massive stars impart kinetic energy to the lower-mass stars and brown dwarfs/free floating planets. The predicted net result of this is that the centrally concentrated massive stars should have significantly lower velocities than fast-moving low-mass objects on the periphery of the cluster. We search for energy equipartition in initially spatially and kinematically substructured N-body simulations of star clusters with N = 1500 stars, evolved for 100 Myr. In clusters that show significant mass segregation we find no differences in the proper motions or radial velocities as a function of mass. The kinetic energies of all stars decrease as the clusters relax, but the kinetic energies of the most massive stars do not decrease faster than those of lower-mass stars. These results suggest that dynamical mass segregation -- which is observed in many star clusters -- is not a signature of energy equipartition from two-body relaxation.

We follow the dynamical evolution of young star-forming regions with a wide range of initial conditions and examine how the radial velocity dispersion, $\sigma$, evolves over time. We compare this velocity dispersion to the theoretically expected value for the velocity dispersion if a region were in virial equilibrium, $\sigma_{\rm vir}$ and thus assess the virial state ($\sigma / \sigma_{\rm vir}$) of these systems. We find that in regions that are initially subvirial, or in global virial equilibrium but subvirial on local scales, the system relaxes to virial equilibrium within several million years, or roughly 25 - 50 crossing times, according to the measured virial ratio. However, the measured velocity dispersion, $\sigma$, appears to be a bad diagnostic of the current virial state of these systems as it suggests that they become supervirial when compared to the velocity dispersion estimated from the virial mass, $\sigma_{\rm vir}$. We suggest that this discrepancy is caused by the fact that the regions are never fully relaxed, and that the early non-equilibrium evolution is imprinted in the one-dimensional velocity dispersion at these early epochs. If measured early enough ($<$2 Myr in our simulations, or $\sim$20 crossing times), the velocity dispersion can be used to determine whether a region was highly supervirial at birth without the risk of degeneracy. We show that combining $\sigma$, or the ratio of $\sigma$ to the interquartile range (IQR) dispersion, with measures of spatial structure, places stronger constraints on the dynamical history of a region than using the velocity dispersion in isolation.

The presence and abundance of the short-lived radioisotopes (SLRs) $^{26}$Al and $^{60}$Fe during the formation of the Solar System is difficult to explain unless the Sun formed in the vicinity of one or more massive star(s) that exploded as supernovae. Two different scenarios have been proposed to explain the delivery of SLRs to the protosolar nebula: (i) direct pollution of the protosolar disc by supernova ejecta and (ii) the formation of the Sun in a sequential star formation event in which supernovae shockwaves trigger further star formation which is enriched in SLRs. The sequentially triggered model has been suggested as being more astrophysically likely than the direct pollution scenario. In this paper we investigate this claim by analysing a combination of $N$-body and SPH simulations of star formation. We find that sequential star formation would result in large age spreads (or even bi-modal age distributions for spatially coincident events) due to the dynamical relaxation of the first star-formation event(s). Secondly, we discuss the probability of triggering spatially and temporally discrete populations of stars and find this to be only possible in very contrived situations. Taken together, these results suggest that the formation of the Solar System in a triggered star formation event is as improbable, if not more so, than the direct pollution of the protosolar disc by a supernova.

We compare the spatial distribution of stars which form in hydrodynamical simulations to the spatial distribution of the gas, using the $\mathcal{Q}$-parameter. The $\mathcal{Q}$-parameter enables a self-consistent comparison between the stars and gas because it uses a pixelated image of the gas as a distribution of points, in the same way that the stars (sink particles in the simulations) are a distribution of points. We find that, whereas the stars have a substructured, or hierarchical spatial distribution ($\mathcal{Q} \sim 0.4 - 0.7$), the gas is dominated by a smooth, concentrated component and typically has $\mathcal{Q} \sim 0.9$. We also find no statistical difference between the structure of the gas in simulations that form with feedback, and those that form without, despite these two processes producing visually different distributions. These results suggest that the link between the spatial distributions of gas, and the stars which form from them, is non-trivial.

We examine the performance of four different methods which are used to measure mass segregation in star-forming regions: the radial variation of the mass function $\mathcal{M}_{\rm MF}$; the minimum spanning tree-based $\Lambda_{\rm MSR}$ method; the local surface density $\Sigma_{\rm LDR}$ method; and the $\Omega_{\rm GSR}$ technique, which isolates groups of stars and determines whether the most massive star in each group is more centrally concentrated than the average star. All four methods have been proposed in the literature as techniques for quantifying mass segregation, yet they routinely produce contradictory results as they do not all measure the same thing. We apply each method to synthetic star-forming regions to determine when and why they have shortcomings. When a star-forming region is smooth and centrally concentrated, all four methods correctly identify mass segregation when it is present. However, if the region is spatially substructured, the $\Omega_{\rm GSR}$ method fails because it arbitrarily defines groups in the hierarchical distribution, and usually discards positional information for many of the most massive stars in the region. We also show that the $\Lambda_{\rm MSR}$ and $\Sigma_{\rm LDR}$ methods can sometimes produce apparently contradictory results, because they use different definitions of mass segregation. We conclude that only $\Lambda_{\rm MSR}$ measures mass segregation in the classical sense (without the need for defining the centre of the region), although $\Sigma_{\rm LDR}$ does place limits on the amount of previous dynamical evolution in a star-forming region.

We present a large-scale, volume-limited companion survey of 245 late-K to mid-M (K7-M6) dwarfs within 15 pc. Infrared adaptive optics (AO) data were analysed from the Very Large Telescope, Subaru Telescope, Canada-France-Hawaii Telescope, and MMT Observatory to detect close companions to the sample from $\sim$1 au to 100 au, while digitised wide-field archival plates were searched for wide companions from $\sim$100 au to 10,000 au. With sensitivity to the bottom of the main sequence over a separation range of 3 au to 10,000 au, multiple AO and wide-field epochs allow us to confirm candidates with common proper motions, minimize background contamination, and enable a measurement of comprehensive binary statistics. We detected 65 co-moving stellar companions and find a companion star fraction of $23.5 \pm 3.2$ per cent over the 3 au to 10,000 au separation range. The companion separation distribution is observed to rise to a higher frequency at smaller separations, peaking at closer separations than measured for more massive primaries. The mass ratio distribution across the $q = 0.2 - 1.0$ range is flat, similar to that of multiple systems with solar-type primaries. The characterisation of binary and multiple star frequency for low-mass field stars can provide crucial comparisons with star forming environments and hold implications for the frequency and evolutionary histories of their associated disks and planets.

The gravitational instability model of planet/brown dwarf formation proposes that protostellar discs can fragment into objects with masses above a few Jupiter masses at large semimajor axis. Tidal downsizing may reduce both the object mass and semimajor axis. However, most studies of tidal downsizing end when the protostellar disc disperses, while the system is embedded in its parent star-forming region. To compare disc fragment descendants with exoplanet and brown dwarf observations, the subsequent dynamical evolution must be explored. We carry out N-Body integrations of fragment-fragment scattering in multi-object star systems, and star systems embedded in substructured clusters. In both cases, we use initial conditions generated by population synthesis models of tidal downsizing. The scattering simulations produce a wide range of eccentricities. The ejection rate is around 25%. The ejecta mass distribution is similar to that for all objects, with a velocity dispersion consistent with those produced by full hydrodynamic simulations. The semimajor axis distribution after scattering extends to parsec scales. In the cluster simulations, 13% of objects are ejected from their planetary system, and around 10% experience significant orbit modification. A small number of objects are recaptured on high eccentricity, high inclination orbits. The velocity distribution of ejecta is similar to that produced by fragment-fragment scattering. If fragment-fragment scattering and cluster stripping act together, then disc fragmentation should be efficient at producing free-floating substellar objects, and hence characterising the free-floating planet population will provide strong constraints on the frequency of disc fragmentation.

We take the end result of smoothed particle hydrodynamics (SPH) simulations of star formation which include feedback from photoionisation and stellar winds and evolve them for a further 10Myr using $N$-body simulations. We compare the evolution of each simulation to a control run without feedback, and to a run with photoionisation feedback only. In common with previous work, we find that the presence of feedback prevents the runaway growth of massive stars, and the resulting star-forming regions are less dense, and preserve their initial substructure for longer. The addition of stellar winds to the feedback produces only marginal differences compared to the simulations with just photoionisation feedback. We search for mass segregation at different stages in the simulations; before feedback is switched on in the SPH runs, at the end of the SPH runs (before $N$-body integration) and during the $N$-body evolution. Whether a simulation is primordially mass segregated (i.e. before dynamical evolution) depends extensively on how mass segregation is defined, and different methods for measuring mass segregation give apparently contradictory results. Primordial mass segregation is also less common in the simulations when star formation occurs under the influence of feedback. Further dynamical mass segregation can also take place during the subsequent (gas-free) dynamical evolution. Taken together, our results suggest that extreme caution should be exercised when interpreting the spatial distribution of massive stars relative to low-mass stars in simulations.

The initial density of individual star-forming regions (and by extension the birth environment of planetary systems) is difficult to constrain due to the "density degeneracy problem": an initially dense region expands faster than a more quiescent region due to two-body relaxation and so two regions with the same observed present-day density may have had very different initial densities. We constrain the initial densities of seven nearby star-forming regions by folding in information on their spatial structure from the $\mathcal{Q}$-parameter and comparing the structure and present-day density to the results of $N$-body simulations. This in turn places strong constraints on the possible effects of dynamical interactions and radiation fields from massive stars on multiple systems and protoplanetary discs. We apply our method to constrain the initial binary population in each of these seven regions and show that the populations in only three - the Orion Nebula Cluster, $\rho$ Oph and Corona Australis - are consistent with having evolved from the Kroupa universal initial period distribution and a binary fraction of unity.

Recent observations of binary stars in the Galactic field show that the binary fraction and the mean orbital separation both decrease as a function of decreasing primary mass. We present $N$-body simulations of the effects of dynamical evolution in star-forming regions on primordial binary stars to determine whether these observed trends can be explained by the dynamical processing of a common binary population. We find that dynamical processing of a binary population with an initial binary fraction of unity and an initial excess of intermediate/wide separation (100 - 10$^4$ au) binaries does not reproduce the observed properties in the field, even in initially dense ($\sim 10^3$M$_\odot$ pc$^{-3}$) star-forming regions. If instead we adopt a field-like population as the initial conditions, most brown dwarf and M-dwarf binaries are dynamically hard and their overall fractions and separation distributions are unaffected by dynamical evolution. G-dwarf and A-star binaries in the field are dynamically intermediate in our simulated dense regions and dynamical processing does destroy some systems with separations $>$100 au. However, the formation of wide binaries through the dissolution of supervirial regions is a strong function of primary mass, and the wide G-dwarf and A-star binaries that are destroyed by dynamical evolution in subvirial regions are replenished by the formation of binaries in supervirial regions. We therefore suggest that the binary population in the field is $indicative$ of the primordial binary population in star-forming regions, at least for systems with primary masses in the range 0.02 - 3.0M$_\odot$.

One of the fundamental questions in star formation is whether or not brown dwarfs form in the same way as stars, or more like giant planets. If their formation scenarios are different, we might expect brown dwarfs to have a different spatial distribution to stars in nearby star-forming regions. In this contribution, we discuss methods to look for differences in their spatial distributions and show that in the only nearby star-forming region with a significantly different spatial distribution (the Orion Nebula Cluster), this is likely due to dynamical evolution. We then present a method for unravelling the past dynamical history of a star-forming region, and show that in tandem with Gaia, we will be able to discern whether observed differences are due to distinct formation mechanisms for brown dwarfs compared to stars.

We use $N$-body simulations to compare the evolution of spatial distributions of stars and brown dwarfs in young star-forming regions. We use three different diagnostics; the ratio of stars to brown dwarfs as a function of distance from the region's centre, $\mathcal{R}_{\rm SSR}$, the local surface density of stars compared to brown dwarfs, $\Sigma_{\rm LDR}$, and we compare the global spatial distributions using the $\Lambda_{\rm MSR}$ method. From a suite of twenty initially statistically identical simulations, 6/20 attain $\mathcal{R}_{\rm SSR} << 1$ $and$ $\Sigma_{\rm LDR} << 1$ $and$ $\Lambda_{\rm MSR} << 1$, indicating that dynamical interactions could be responsible for observed differences in the spatial distributions of stars and brown dwarfs in star-forming regions. However, many simulations also display apparently contradictory results - for example, in some cases the brown dwarfs have much lower local densities than stars ($\Sigma_{\rm LDR} << 1$), but their global spatial distributions are indistinguishable ($\Lambda_{\rm MSR} = 1$) and the relative proportion of stars and brown dwarfs remains constant across the region ($\mathcal{R}_{\rm SSR} = 1$). Our results suggest that extreme caution should be exercised when interpreting any observed difference in the spatial distribution of stars and brown dwarfs, and that a much larger observational sample of regions/clusters (with complete mass functions) is necessary to investigate whether or not brown dwarfs form through similar mechanisms to stars.

We examine substructure and mass segregation in the massive OB association Cygnus OB2 to better understand its initial conditions. Using a well understood Chandra X-ray selected sample of young stars we find that Cyg OB2 exhibits considerable physical substructure and has no evidence for mass segregation, both indications that the association is not dynamically evolved. Combined with previous kinematical studies we conclude that Cyg OB2 is dynamically very young, and what we observe now is very close to its initial conditions: Cyg OB2 formed as a highly substructured, unbound association with a low volume density (< 100 stars/pc^3). This is inconsistent with the idea that all stars form in dense, compact clusters. The massive stars in Cyg OB2 show no evidence for having formed particularly close to one another, nor in regions of higher than average density. Since Cyg OB2 contains stars as massive as ~100 Mo this result suggests that very massive stars can be born in relatively low-density environments. This would imply that the massive stars in Cyg OB2 did not form by competitive accretion, or by mergers.

We model the dynamical evolution of star forming regions with a wide range of initial properties. We follow the evolution of the regions' substructure using the Q-parameter, we search for dynamical mass segregation using the Lambda_MSR technique, and we also quantify the evolution of local density around stars as a function of mass using the Sigma_LDR method. The amount of dynamical mass segregation measured by Lambda_MSR is generally only significant for subvirial and virialised, substructured regions - which usually evolve to form bound clusters. The Sigma_LDR method shows that massive stars attain higher local densities than the median value in all regions, even those that are supervirial and evolve to form (unbound) associations. We also introduce the Q-Sigma_LDR plot, which describes the evolution of spatial structure as a function of mass-weighted local density in a star forming region. Initially dense (>1000 stars pc^{-2}), bound regions always have Q >1, Sigma_LDR > 2 after 5Myr, whereas dense unbound regions always have Q < 1, Sigma_LDR > 2 after 5Myr. Less dense regions (<100 stars pc^{-2}) do not usually exhibit Sigma_LDR > 2 values, and if relatively high local density around massive stars arises purely from dynamics, then the Q-Sigma_LDR plot can be used to estimate the initial density of a star forming region.

We use N-body simulations of star cluster evolution to explore the hypothesis that short-lived radioactive isotopes found in meteorites, such as 26-Al, were delivered to the Sun's protoplanetary disc from a supernova at the epoch of Solar System formation. We cover a range of star cluster formation parameter space and model both clusters with primordial substructure, and those with smooth profiles. We also adopt different initial virial ratios - from cool, collapsing clusters to warm, expanding associations. In each cluster we place the same stellar population; the clusters each have 2100 stars, and contain one massive 25M_Sun star which is expected to explode as a supernova at about 6.6Myr. We determine the number of Solar (G)-type stars that are within 0.1 - 0.3pc of the 25M_Sun star at the time of the supernova, which is the distance required to enrich the protoplanetary disc with the 26-Al abundances found in meteorites. We then determine how many of these G-dwarfs are unperturbed `singletons'; stars which are never in close binaries, nor suffer sub-100au encounters, and which also do not suffer strong dynamical perturbations. The evolution of a suite of twenty initially identical clusters is highly stochastic, with the supernova enriching over 10 G-dwarfs in some clusters, and none at all in others. Typically only ~25 per cent of clusters contain enriched, unperturbed singletons, and usually only 1 - 2 per cluster (from a total of 96 G-dwarfs in each cluster). The initial conditions for star formation do not strongly affect the results, although a higher fraction of supervirial (expanding) clusters would contain enriched G-dwarfs if the supernova occurred earlier than 6.6Myr. If we sum together simulations with identical initial conditions, then ~1 per cent of all G-dwarfs in our simulations are enriched, unperturbed singletons.

We determine the fraction of G-dwarf stars that could host stable planetary systems based on the observed properties of binaries in the Galactic field, and in various postulated primordial binary populations, which assume that the primordial binary fraction is higher than that in the field. We first consider the frequency of Solar System analogues - planetary systems that form either around a single G-dwarf star, or a binary containing a G-dwarf where the binary separation exceeds 100-300au. If the primordial binary fraction and period distribution is similar to that in the field, then up to 63 per cent of G-dwarf systems could potentially host a Solar System analogue. However, if the primordial binary fraction is higher, the fraction of G-dwarf systems that could host a planetary system like our own is lowered to 38 per cent. We extend our analysis to consider the fraction of G-dwarf systems (both single and binary) that can host either circumprimary planets (orbiting the primary star of the binary) or circumbinary planets (orbiting both stars in the binary) for fiducial planetary separations between 1 - 100au. Depending on the assumed binary population, in the circumprimary case between 65 and 95 per cent of systems can host a planet at 1au, decreasing to between 20 and 65 per cent of systems that can host a planet at 100au. In the circumbinary case, between 5 and 59 per cent of systems can host a planet at 1au, increasing to between 34 and 75 per cent of systems that can host a planet at 100au. Our results suggest that the assumed binary fraction is the key parameter in determining the fraction of potentially stable planetary systems in G-dwarf systems and that using the present-day value may lead to significant overestimates if the binary fraction was initially higher.

We explore the effects of dynamical evolution in dense clusters on the companion mass ratio distribution (CMRD) of binary stars. Binary systems are destroyed by interactions with other stars in the cluster, lowering the total binary fraction and significantly altering the initial semi-major axis distribution. However, the shape of the CMRD is unaffected by dynamics; an equal number of systems with high mass ratios are destroyed compared to systems with low mass ratios. We might expect a weak dependence of the survivability of a binary on its mass ratio because its binding energy is proportional to both the primary and secondary mass components of the system. However, binaries are broken up by interactions in which the perturbing star has a significantly higher energy (by a factor of >10, depending on the particular binary properties) than the binding energy of the binary, or through multiple interactions in the cluster. We therefore suggest that the shape of the observed binary CMRD is an outcome of the star formation process, and should be measured in preference to the distributions of orbital parameters, such as the semi-major axis distribution.

We present the results of N-body simulations in which we take the masses, positions and velocities of sink particles from five pairs of hydrodynamical simulations of star formation by Dale et al. (2012, 2013) and evolve them for a further 10Myr. We compare the dynamical evolution of star clusters that formed under the influence of mass-loss driven by photoionization feedback, to the evolution of clusters that formed without feedback. We remove any remaining gas and follow the evolution of structure in the clusters (measured by the Q-parameter), half-mass radius, central density, surface density and the fraction of bound stars. There is little discernible difference in the evolution of clusters that formed with feedback compared to those that formed without. The only clear trend is that all clusters which form without feedback in the hydrodynamical simulations lose any initial structure over 10Myr, whereas some of the clusters which form with feedback retain structure for the duration of the subsequent N-body simulation. This is due to lower initial densities (and hence longer relaxation times) in the clusters from Dale et al. (2012, 2013) which formed with feedback, which prevents dynamical mixing from erasing substructure. However, several other conditions (such as supervirial initial velocities) also preserve substructure, so at a given epoch one would require knowledge of the initial density and virial state of the cluster in order to determine whether star formation in a cluster has been strongly influenced by feedback.

Radial velocity measurements can be used to constrain the dynamical state of a stellar cluster. However, for clusters with velocity dispersions smaller than a few km/s the observed radial velocity distribution tends to be dominated by the orbital motions of binaries rather than the stellar motions through the potential well of the cluster. Our goal is to characterize the intrinsic velocity distribution of a cluster from a single epoch of radial velocity data, even for a cluster with a velocity dispersion of a fraction of a km/s, using a maximum likelihood procedure. Assuming a period, mass ratio, and eccentricity distribution for the binaries in the observed cluster this procedure fits a dynamical model describing the velocity distribution for the single stars and center of masses of the binaries, simultaneously with the radial velocities caused by binary orbital motions, using all the information available in the observed velocity distribution. We find that the fits to the intrinsic velocity distribution depend only weakly on the binary properties assumed, so the uncertainty in the fitted parameters tends to be dominated by statistical uncertainties. Based on Monte Carlo simulations we provide an estimate of how these statistical uncertainties vary with the velocity dispersion, binary fraction, and the number of observed stars, which can be used to estimate the sample size needed to reach a specific accuracy. Finally we test the method on the well-studied open cluster NGC 188, showing that it can reproduce a velocity dispersion of only 0.5 km/s using a single epoch of the multi-epoch radial velocity data. If the binary period, mass ratio, and eccentricity distribution of the observed stars are roughly known, this procedure can be used to correct for the effect of binary orbital motions on an observed velocity distribution. [Abridged]

We have measured the multiplicity fractions and separation distributions of seven young star-forming regions using a uniform sample of young binaries. Both the multiplicity fractions and separation distributions are similar in the different regions. A tentative decline in the multiplicity fraction with increasing stellar density is apparent, even for binary systems with separations too close (19-100au) to have been dynamically processed. The separation distributions in the different regions are statistically indistinguishable over most separation ranges, and the regions with higher densities do not exhibit a lower proportion of wide (300-620au) relative to close (62-300au) binaries as might be expected from the preferential destruction of wider pairs. Only the closest (19-100au) separation range, which would be unaffected by dynamical processing, shows a possible difference in separation distributions between different regions. The combined set of young binaries, however, shows a distinct difference when compared to field binaries, with a significant excess of close (19-100au) systems among the younger binaries. Based on both the similarities and differences between individual regions, and between all seven young regions and the field, especially over separation ranges too close to be modified by dynamical processing, we conclude that multiple star formation is not universal and, by extension, the star formation process is not universal.

We determine the distribution of stellar surface densities, \Sigma, from models of static and dynamically evolving star clusters with different morphologies, including both radially smooth and substructured clusters. We find that the \Sigma distribution is degenerate, in the sense that many different cluster morphologies (smooth or substructured) produce similar cumulative distributions. However, when used in tandem with a measure of structure, such as the Q-parameter, the current spatial and dynamical state of a star cluster can be inferred. The effect of cluster dynamics on the \Sigma distribution and the Q-parameter is investigated using N-body simulations and we find that, depending on the assumed initial conditions, the \Sigma distribution can rapidly evolve from high to low densities in less than 5Myr. This suggests that the \Sigma distribution can only be used to assess the current density of a star forming region, and provides little information on its initial density. However, if the \Sigma distribution is used together with the Q-parameter, then information on the amount of substructure can be used as a proxy to infer the amount of dynamical evolution that has taken place. Substructure is erased quickly through dynamics, which can disrupt binary star systems and planets, as well as facilitate dynamical mass segregation. Therefore, dynamical processing in young star forming regions could still be significant, even without currently observed high densities.

In this contribution I present a review of star formation in clusters. I begin by discussing the various definitions of what constitutes a star cluster, and then compare the outcome of star formation (IMF, multiplicity, mass segregation and structure and morphology) in different star-forming regions. I also review recent numerical models of star formation in clusters, before ending with a summary of the potential effects of dynamical evolution in star clusters.

We apply two different algorithms to search for mass segregation to a recent observational census of the rho Ophiuchi star forming region. Firstly, we apply the Lambda_MSR method, which compares the minimum spanning tree (MST) of a chosen subset of stars to MSTs of random subsets of stars in the cluster, and determine the mass segregation ratio, Lambda_MSR. Secondly, we apply the m-Sigma method, which calculates the local stellar surface density around each star and determines the statistical significance of the average surface density for a chosen mass bin, compared to the average surface density in the whole cluster. Using both methods, we find no indication of mass segregation (normal or inverse) in the spatial distribution of stars and brown dwarfs in rho Ophiuchi. Although rho Ophiuchi suffers from high visual extinction, we show that a significant mass segregation signature would be detectable, albeit slightly diluted, despite dust obscuration of centrally located massive stars.

Observations of binaries in clusters tend to be of visual binaries with separations of 10s - 100s au. Such binaries are 'intermediates' and their destruction or survival depends on the exact details of their individual dynamical history. We investigate the stochasticity of the destruction of such binaries and the differences between the initial and processed populations using N-body simulations. We concentrate on Orion Nebula Cluster-like clusters, where the observed binary separation distribution ranges from 62 - 620 au. We find that, starting from the same initial binary population in statistically identical clusters, the number of intermediate binaries that are destroyed after 1 Myr can vary by a factor of >2, and that the resulting separation distributions can be statistically completely different in initially substructured clusters. We also find that the mass ratio distributions are altered (destroying more low mass ratio systems), but not as significantly as the binary fractions or separation distributions. We conclude that finding very different intermediate (visual) binary populations in different clusters does not provide conclusive evidence that the initial populations were different.

Whether massive stars can occasionally form in relative isolation or if they require a large cluster of lower-mass stars around them is a key test in the differentiation of star formation theories as well as how the initial mass function of stars is sampled. Previous attempts to find O-type stars that formed in isolation were hindered by the possibility that such stars are merely runaways from clusters, i.e., their current isolation does not reflect their birth conditions. We introduce a new method to find O-type stars that are not affected by such a degeneracy. Using the VLT-FLAMES Tarantula Survey and additional high resolution imaging we have identified stars that satisfy the following constraints: 1) they are O-type stars that are not detected to be part of a binary system based on RV time series analysis; 2) they are designated spectral type O7 or earlier ; 3) their velocities are within 1\sigma of the mean of OB-type stars in the 30 Doradus region, i.e. they are not runaways along our line-of-sight; 4) the projected surface density of stars does not increase within 3 pc towards the O-star (no evidence for clusters); 5) their sight lines are associated with gaseous and/or dusty filaments in the ISM, and 6) if a second candidate is found in the direction of the same filament with which the target is associated, both are required to have similar velocities. With these criteria, we have identified 15 stars in the 30 Doradus region, which are strong candidates for being high-mass stars that have formed in isolation. Additionally, we employed extensive MC stellar cluster simulations to confirm that our results rule out the presence of clusters around the candidates. Eleven of these are classified as Vz stars, possibly associated with the zero-age main sequence. We include a newly discovered W-R star as a candidate, although it does not meet all of the above criteria.

We have collated multiplicity data for five clusters (Taurus, Chamaeleon I, Ophiuchus, IC348, and the Orion Nebula Cluster). We have applied the same mass ratio (flux ratios of delta K <= 2.5) and primary mass cuts (~0.1-3.0 Msun) to each cluster and therefore have directly comparable binary statistics for all five clusters in the separation range 62-620 au, and for Taurus, Chamaeleon I, and Ophiuchus in the range 18-830 au. We find that the trend of decreasing binary fraction with cluster density is solely due to the high binary fraction of Taurus, the other clusters show no obvious trend over a factor of nearly 20 in density. With N-body simulations we attempt to find a set of initial conditions that are able to reproduce the density, morphology and binary fractions of all five clusters. Only an initially clumpy (fractal) distribution with an initial total binary fraction of 73 per cent (17 per cent in the range 62-620 au) is able to reproduce all of the observations (albeit not very satisfactorily). Therefore, if star formation is universal the initial conditions must be clumpy and with a high (but not 100 per cent) binary fraction. This could suggest that most stars, including M-dwarfs, form in binaries.

We present N-body simulations of young substructured star clusters undergoing various dynamical evolutionary scenarios and examine the direct effects of interactions in the cluster on planetary systems. We model clusters initially in cool collapse, in virial equilibrium and expanding, and place a 1 Jupiter-mass planet at either 5au or 30au from their host stars, with zero eccentricity. We find that after 10Myr 10% of planets initially orbiting at 30au have been liberated from their parent star and form a population of free-floating planets. A small number of these planets are captured by other stars. A further 10% have their orbital eccentricity significantly altered. The change in eccentricity is often accompanied by a change in orbital inclination which may lead to additional dynamical perturbations in planetary systems. The fraction of liberated and disrupted planetary systems is highest for subvirial clusters, but virial and supervirial clusters also dynamically process planetary systems, due to interactions in the substructure. Of the planets that become free-floating, those that remain observationally associated with the cluster (within two half-mass radii of the cluster centre) have a similar velocity distribution to the entire star cluster, irrespective of whether they were on a 5au or a 30au orbit, with median velocities typically 1km/s. Conversely, those planets that are no longer associated with the cluster have similar velocities to the non-associated stars if they were originally at 5au (9km/s), whereas planets originally at 30au have much lower velocities (3.8km/s) than the stars (10.8km/s). These findings highlight potential pitfalls of concluding that (a) planets with similar velocities to the cluster stars represent the very low-mass end of the IMF, and (b) planets on the periphery of a cluster with very different observed velocities form through different mechanisms.

A large population of fragile, wide (> 1000 AU) binary systems exists in the Galactic field and halo. These wide binary stars cannot be primordial because of the high stellar density in star forming regions, while formation by capture in the Galactic field is highly improbable. We propose that these binary systems were formed during the dissolution phase of star clusters (see Kouwenhoven et al. 2010, for details). Stars escaping from a dissolving star cluster can have very similar velocities, which can lead to the formation of a wide binary systems. We carry out N-body simulations to test this hypothesis. The results indicate that this mechanism explains the origin of wide binary systems in the Galaxy. The resulting wide binary fraction and semi-major axis distribution depend on the initial conditions of the dissolving star cluster, while the distributions in eccentricity and mass ratio are universal. Finally, since most stars are formed in (relatively tight) primordial binaries, we predict that a large fraction of the wide "binary stars" are in fact higher-order multiple systems.

Observations and theory suggest that star clusters can form in a subvirial (cool) state and are highly substructured. Such initial conditions have been proposed to explain the level of mass segregation in clusters through dynamics, and have also been successful in explaining the origin of trapezium-like systems. In this paper we investigate, using N-body simulations, whether such a dynamical scenario is consistent with the observed binary properties in the Orion Nebula Cluster (ONC). We find that several different primordial binary populations are consistent with the overall fraction and separation distribution of visual binaries in the ONC (in the range 67 - 670 au), and that these binary systems are heavily processed. The substructured, cool-collapse scenario requires a primordial binary fraction approaching 100 per cent. We find that the most important factor in processing the primordial binaries is the initial level of substructure; a highly substructured cluster processes up to 20 per cent more systems than a less substructured cluster because of localised pockets of high stellar density in the substructure. Binaries are processed in the substructure before the cluster reaches its densest phase, suggesting that even clusters remaining in virial equilibrium or undergoing supervirial expansion would dynamically alter their primordial binary population. Therefore even some expanding associations may not preserve their primordial binary population.

We use the new minimum spanning tree (MST) method to look for mass segregation in the Taurus association. The method computes the ratio of MST lengths of any chosen subset of objects, including the most massive stars and brown dwarfs, to the MST lengths of random sets of stars and brown dwarfs in the cluster. This mass segregation ratio (Lambda_MSR) enables a quantitative measure of the spatial distribution of high-mass and low-mass stars, and brown dwarfs to be made in Taurus. We find that the most massive stars in Taurus are inversely mass segregated, with Lambda_MSR = 0.70 +/- 0.10 (Lambda_MSR = 1 corresponds to no mass segregation), which differs from the strong mass segregation signatures found in more dense and massive clusters such as Orion. The brown dwarfs in Taurus are not mass segregated, although we find evidence that some low-mass stars are, with an Lambda_MSR = 1.25 +/- 0.15. Finally, we compare our results to previous measures of the spatial distribution of stars and brown dwarfs in Taurus, and briefly discuss their implications.

Very low-mass binaries (VLMBs), with system masses <0.2 Msun appear to have very different properties to stellar binaries. This has led to the suggestion that VLMBs form a distinct and different population. As most stars are born in clusters, dynamical evolution can significantly alter any initial binary population, preferentially destroying wide binaries. In this paper we examine the dynamical evolution of initially different VLMB distributions in clusters to investigate how different the initial and final distributions can be. We find that the majority of the observed VLMB systems, which have separations <20 au, cannot be destroyed in even the densest clusters. Therefore, the distribution of VLMBs with separations <20 au now must have been the birth population (although we note that the observations of this population may be very incomplete). Most VLMBs with separations >100 au can be destroyed in high-density clusters, but are mainly unaffected in low-density clusters. Therefore, the initial VLMB population must contain many more binaries with these separations than now, or such systems must be made by capture during cluster dissolution. M-dwarf binaries are processed in the same way as VLMBs and so the difference in the current field populations either points to fundamentally different birth populations, or significant observational incompleteness in one or both samples.

Spectroscopic analyses of H-rich WN5-6 stars within the young star clusters NGC 3603 and R136 are presented, using archival HST & VLT spectroscopy, & high spatial resolution near-IR photometry. We derive high T* for the WN stars in NGC 3603 (T*~42+/-2 kK) & R136 (T*~53+/-3 kK) plus clumping-corrected dM/dt ~ 2-5x10^-5 Msun/yr which closely agree with theoretical predictions. These stars make a disproportionate contribution to the global budget of their host clusters. R136a1 alone supplies ~7% of N(LyC) of the entire 30 Dor region. Comparisons with stellar models calculated for the main-sequence evolution of 85-500 Msun suggest ages of ~1.5 Myr & M_init in the range 105 - 170 Msun for 3 systems in NGC 3603, plus 165-320 Msun for 4 stars in R136. Our high stellar masses are supported by dynamical mass determinations for the components of NGC 3603 A1. We consider the predicted L_X of the R136 stars if they were close, colliding wind binaries. R136c is consistent with a colliding wind binary system. However, short period, colliding wind systems are excluded for R136a WN stars if mass ratios are of order unity. Widely separated systems would have been expected to harden owing to early dynamical encounters with other massive stars in such a dense environment. From simulated star clusters, whose constituents are randomly sampled from the Kroupa IMF, both clusters are consistent with a tentative upper mass limit of ~300 Msun. The Arches cluster is either too old, exhibits a deficiency of very massive stars, or more likely stellar masses have been underestimated - M_init for the most luminous stars in the Arches cluster approach 200 Msun according to contemporary stellar & photometric results. The potential for stars greatly exceeding 150 Msun within metal-poor galaxies suggests that such pair-instability SNe could occur within the local universe, as has been claimed for SN 2007bi (abridged).

Observations and theory both suggest that star clusters form sub-virial (cool) with highly sub-structured distributions. We perform a large ensemble of N-body simulations of moderate-sized (N=1000) cool, fractal clusters to investigate their early dynamical evolution. We find that cool, clumpy clusters dynamically mass segregate on a short timescale, that Trapezium-like massive higher-order multiples are commonly formed, and that massive stars are often ejected from clusters with velocities > 10 km/s (c.f. the average escape velocity of 2.5 km/s). The properties of clusters also change rapidly on very short timescales. Young clusters may also undergo core collapse events, in which a dense core containing massive stars is hardened due to energy losses to a halo of lower-mass stars. Such events can blow young clusters apart with no need for gas expulsion. The warmer and less substructured a cluster is initially, the less extreme its evolution.

Over the past few decades, numerous wide (>1000 au) binaries in the Galactic field and halo have been discovered. Their existence cannot be explained by the process of star formation or by dynamical interactions in the field, and their origin has long been a mystery. We explain the origin of these wide binaries by formation during the dissolution phase of young star clusters: an initially unbound pair of stars may form a binary when their distance in phase-space is small. Using N-body simulations, we find that the resulting wide binary fraction in the semi-major axis range 1000 au - 0.1 pc for individual clusters is 1-30%, depending on the initial conditions. The existence of numerous wide binaries in the field is consistent with observational evidence that most clusters start out with a large degree of substructure. The wide binary fraction decreases strongly with increasing cluster mass, and the semi-major axis of the newly formed binaries is determined by the initial cluster size. The resulting eccentricity distribution is thermal, and the mass ratio distribution is consistent with gravitationally-focused random pairing. As a large fraction of the stars form in primordial binaries, we predict that a large number of the observed 'wide binaries' are in fact triple or quadruple systems. By integrating over the initial cluster mass distribution, we predict a binary fraction of a few per cent in the semi-major axis range 1000 au - 0.1 pc in the Galactic field, which is smaller than the observed wide binary fraction. However, this discrepancy may be solved when we consider a broad range of cluster morphologies.

The majority of stars in the Galactic field and halo are part of binary or multiple systems. A significant fraction of these systems have orbital separations in excess of thousands of astronomical units, and systems wider than a parsec have been identified in the Galactic halo. These binary systems cannot have formed through the 'normal' star-formation process, nor by capture processes in the Galactic field. We propose that these wide systems were formed during the dissolution phase of young star clusters. We test this hypothesis using N-body simulations of evolving star clusters and find wide binary fractions of 1-30%, depending on initial conditions. Moreover, given that most stars form as part of a binary system, our theory predicts that a large fraction of the known wide 'binaries' are, in fact, multiple systems.

Over the last decades, numerous wide (>1000 AU) binaries have been discovered in the Galactic field and halo. The origin of these wide binaries cannot be explained by star formation or by dynamical interactions in the Galactic field. We explain their existence by wide binary formation during the dissolution phase of young star clusters. In this scenario, two single stars that leave the dissolving cluster at the same time, in the same direction, and with similar velocities, form a new, very wide binary. Using N-body simulations we study how frequently this occurs, and how the orbital parameters of such binaries depend on the properties of the cluster from which they originate. The resulting wide binary fraction for individual star clusters is 1-30%, depending on the initial conditions. As most stars form as part of a binary or multiple system, we predict that a large fraction of these wide binaries are in fact wide triple and quadruple systems.

We discuss the observations and theory of star cluster formation to argue that clusters form dynamically cool (subvirial) and with substructure. We then perform an ensemble of simulations of cool, clumpy (fractal) clusters and show that they often dynamically mass segregate on timescales far shorter than expected from simple models. The mass segregation comes about through the production of a short-lived, but very dense core. This shows that in clusters like the Orion Nebula Cluster the stars >4 Msun can dynamically mass segregate within the current age of the cluster. Therefore, the observed mass segregation in apparently dynamically young clusters need not be primordial, but could be the result of rapid and violent early dynamical evolution.

We examine the effects of dynamical evolution in clusters on planetary systems or protoplanetary disks orbiting the components of binary stars. In particular, we look for evidence that the companions of host stars of planetary systems or disks could have their inclination angles raised from zero to between the threshold angles (39.23 degrees and 140.77 degrees) that can induce the Kozai mechanism. We find that up to 20 per cent of binary systems have their inclination angles increased to within the threshold range. Given that half of all extrasolar planets could be in binary systems, we suggest that up to 10 per cent of extrasolar planets could be affected by this mechanism.

We examine the dynamical destruction of binary systems in star clusters of different densities. We find that at high densities (10^4 - 10^5 Msun pc^-3) almost all binaries with separations > 10^3 AU are destroyed after a few crossing times. At low densities (order(10^2) Msun pc^-3) many binaries with separations > 10^3 AU are destroyed, and no binaries with separations > 10^4 AU survive after a few crossing times. Therefore the binary separations in clusters can be used as a tracer of the dynamical age and past density of a cluster. We argue that the central region of the Orion Nebula Cluster was around 100 times denser in the past with a half-mass radius of only 0.1 - 0.2 pc as (a) it is expanding, (b) it has very few binaries with separations > 10^3 AU, and (c) it is well-mixed and therefore dynamically old. We also examine the origin of the field binary population. Binaries with separations < 10^2 AU are not significantly modified in any cluster, therefore at these separations the field reflects the sum of all star formation. Binaries with separations in the range 10^2 - 10^4 AU are progressively more and more heavily affected by dynamical disruption in increasingly dense clusters. If most star formation is clustered, these binaries must be over-produced relative to the field. Finally, no binary with a separation > 10^4 AU can survive in any cluster and so must be produced by isolated star formation, but only if all isolated star formation produces extremely wide binaries.

We present a new method to detect and quantify mass segregation in star clusters. It compares the minimum spanning tree (MST) of massive stars with that of random stars. If mass segregation is present, the MST length of the most massive stars will be shorter than that of random stars. This difference can be quantified (with an associated significance) to measure the degree of mass segregation. We test the method on simulated clusters in both 2D and 3D and show that the method works as expected. We apply the method to the Orion Nebula Cluster (ONC) and show that the method is able to detect the mass segregation in the Trapezium with a `mass segregation ratio' \Lambda_{MSR}=8.0 \pm 3.5 (where \Lambda_{MSR}=1 is no mass segregation) down to 16 \Msun, and also that the ONC is mass segregated at a lower level (~2.0 \pm 0.5) down to 5 \Msun. Below 5 \Msun we find no evidence for any further mass segregation in the ONC.

Around 4% of O-stars are observed in apparent isolation, with no associated cluster, and no indication of having been ejected from a nearby cluster. We define an isolated O-star as a star > 17.5 M_\odot in a cluster with total mass <100 M_\odot which contains no other massive (>10 M_\odot) stars. We show that the fraction of apparently isolated O-stars is reproduced when stars are sampled (randomly) from a standard initial mass function and a standard cluster mass function of the form N(M) \propto M^-2. This result is difficult to reconcile with the idea that there is a fundamental relationship between the mass of a cluster and the mass of the most massive star in that cluster. We suggest that such a relationship is a typical result of star formation in clusters, and that `isolated O-stars' are low-mass clusters in which massive stars have been able to form.