The kinematic structure of the Cygnus OB2 association is investigated. No evidence of expansion or contraction is found at any scale within the region. Stars that are within $\sim$ 0.5 parsecs of one another are found to have more similar velocities than would be expected by random chance, and so it is concluded that velocity substructure exists on these scales. At larger scales velocity substructure is not found. We suggest that bound substructures exist on scales of $\sim$ 0.5 parsecs, despite the region as a whole being unbound. We further suggest that any velocity substructure that existed on scales > 0.5 parsecs has been erased. The results of this study are then compared to those of other kinematic studies of Cygnus OB2.

We show that dynamical evolution in a strong (Galactic Centre-like) tidal field can create clusters that would appear to have very top-heavy IMFs. The tidal disruption of single star forming events can leave several bound `clusters' spread along 20 pc of the orbit within 1-2 Myr. These surviving (sub)clusters tend to contain an over-abundance of massive stars, with low-mass stars tending to be spread along the whole `tidal arm'. Therefore observing a cluster in a strong tidal field with a top-heavy IMF might well not mean the stars formed with a top-heavy IMF.

We present a detailed study of the Orion B clouds (d~400 pc), imaged with the PACS/SPIRE cameras at 70-500 $\mu$m by the Herschel Gould Belt survey (HGBS). We release new high-res. maps of column density and dust temperature. In the filamentary sub-regions NGC2023/2024, NGC2068/2071, and L1622, 1768 starless dense cores were identified, ~28-45% of which are self-gravitating prestellar cores. A total of 76 protostellar dense cores were also found. The typical lifetime of the prestellar cores was found to be $t_{\rm pre}=1.7_{-0.6}^{+0.8}$ Myr. The prestellar core mass function (CMF) peaks at ~0.5 $M_\odot$ and is consistent with a power law with log slope -1.27$\pm$0.24 at the high-mass end, compared to the Salpeter slope of -1.35. In this region, we confirm the existence of a transition in prestellar core formation efficiency (CFE) around a fiducial value A_V_bg~7 mag in background visual extinction, similar to the trend observed with Herschel in other clouds. This is not a sharp threshold, but a smooth transition between a regime with very low prestellar CFE at A_V_bg<5 and a regime with higher, roughly constant CFE at A_V_bg$\gtrsim$10. The total mass in the form of prestellar cores represents only ~20% of the dense molecular cloud gas at A_V_bg$\gtrsim$7 mag. About 60-80% of the prestellar cores are closely associated with filaments, and this fraction increases up to >90% when a more complete sample of filamentary structures is considered. Interestingly, the median separation between nearest core neighbors corresponds to the typical inner filament width of ~0.1 pc commonly observed in nearby molecular clouds. Analysis of the CMF observed as a function of background cloud column density shows that the most massive prestellar cores are spatially segregated in the highest column density areas, and suggests that both higher- and lower-mass prestellar cores may form in denser filaments.

We present the first search for binary companions of Massive Young Stellar Objects (MYSOs) using AO-assisted $K$ band observations, with NaCo at the VLT. We have surveyed 32 MYSOs from the RMS catalogue, probing the widest companions, with a physical separation range of 400 - 46,000 au, within the predictions of models and observations for multiplicity of MYSOs. Statistical methods are employed to discern whether these companions are physical rather than visual binaries. We find 18 physical companions around 10 target objects, amounting to a multiplicity fraction of 31$\pm$8\% and a companion fraction of 53$\pm$9\%. For similar separation and mass ratio ranges, MYSOs seem to have more companions than T Tauri or O stars, respectively. This suggests that multiplicity increases with mass and decreases with evolutionary stage. We compute very rough estimates for the mass ratios from the $K$ band magnitudes, and these appear to be generally larger than 0.5. This is inconsistent with randomly sampling the IMF, as predicted by the binary capture formation theory. Finally, we find that MYSOs with binaries do not show any different characteristics to the average MYSO in terms of luminosity, distance, outflow or disc presence.

We present a new method of analysing and quantifying velocity structure in star forming regions suitable for the rapidly increasing quantity and quality of stellar position-velocity data. The method can be applied to data in any number of dimensions, does not require the centre or characteristic size (e.g. radius) of the region to be determined, and can be applied to regions with any underlying density and velocity structure. We test the method on a variety of example datasets and show it is robust with realistic observational uncertainties and selection effects. This method identifies velocity structures/scales in a region, and allows a direct comparison to be made between regions.

We present a survey for the tightest visual binaries among 0.3-2 Msun members the Orion Nebula Cluster (ONC). Among 42 targets, we discovered 13 new 0.025-0.15" companions. Accounting for the Branch bias, we find a companion star fraction (CSF) in the 10-60 au range of 21+8/-5%, consistent with that observed in other star-forming regions (SFRs) and twice as high as among field stars; this excess is found with a high level of confidence. Since our sample is dominated by disk-bearing targets, this indicates that disk disruption by close binaries is inefficient, or has not yet taken place, in the ONC. The resulting separation distribution in the ONC drops sharply outside 60\,au. These findings are consistent with a scenario in which the initial multiplicity properties, set by the star formation process itself, are identical in the ONC and in other SFRs and subsequently altered by the cluster's dynamical evolution. This implies that the fragmentation process does not depend on the global properties of a molecular cloud, but on the local properties of prestellar cores, and that the latter are self-regulated to be nearly identical in a wide range of environments. These results, however, raise anew the question of the origin of field stars as the tight binaries we have discovered will not be destroyed as the ONC dissolves into the galactic field. It thus appears that most field stars formed in regions differ from well-studied SFRs in the Solar neighborhood, possibly due to changes in core fragmentation on Gyr timescales.

We examine the dynamical evolution of both Plummer sphere and substructured (fractal) star forming regions in Galactic Centre (GC) strong tidal fields to see what initial conditions could give rise to an Arches-like massive star cluster by $\sim 2$ Myr. We find that any initial distribution has to be contained within its initial tidal radius to survive, which sets a lower limit of the initial density of the Arches of $\sim$ 600 M$_\odot$ pc$^{-3}$ if the Arches is at 30 pc from the GC, or $\sim$ 200 M$_\odot$ pc$^{-3}$ if the Arches is at 100 pc from the GC. Plummer spheres that survive change little other than to dynamically mass segregate, but initially fractal distributions rapidly erase substructure, dynamically mass segregate and by 2 Myr look extremely similar to initial Plummer spheres, therefore it is almost impossible to determine the initial conditions of clusters in strong tidal fields.

Massive stars can be found in wide (hundreds to thousands AU) binaries with other massive stars. We use $N$-body simulations to show that any bound cluster should always have approximately one massive wide binary: one will probably form if none are present initially; and probably only one will survive if more than one are present initially. Therefore any region that contains many massive wide binaries must have been composed of many individual subregions. Observations of Cyg OB2 show that the massive wide binary fraction is at least a half (38/74) which suggests that Cyg OB2 had at least 30 distinct massive star formation sites. This is further evidence that Cyg OB2 has always been a large, low-density association. That Cyg OB2 has a normal high-mass IMF for its total mass suggests that however massive stars form they 'randomly sample' the IMF (as the massive stars did not 'know' about each other).

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.

The origin of very low-mass hydrogen-burning stars, brown dwarfs, and planetary-mass objects at the low-mass end of the initial mass function is not yet fully understood. Gravitational fragmentation of circumstellar discs provides a possible mechanism for the formation of such low-mass objects. The kinematic and binary properties of very low-mass objects formed through disc fragmentation at early times (< 10 Myr) were discussed in Li et al. (2015). In this paper we extend the analysis by following the long-term evolution of disc-fragmented systems, up to an age of 10 Gyr, covering the ages of the stellar and substellar population in the Galactic field. We find that the systems continue to decay, although the rates at which companions escape or collide with each other are substantially lower than during the first 10 Myr, and that dynamical evolution is limited beyond 1 Gyr. By t = 10 Gyr, about one third of the host stars is single, and more than half have only one companion left. Most of the other systems have two companions left that orbit their host star in widely separated orbits. A small fraction of companions have formed binaries that orbit the host star in a hierarchical triple configuration. The majority of such double companion systems have internal orbits that are retrograde with respect to their orbits around their host stars. Our simulations allow a comparison between the predicted outcomes of disc-fragmentation with the observed low-mass hydrogen-burning stars, brown dwarfs, and planetary-mass objects in the Solar neighborhood. Imaging and radial velocity surveys for faint binary companions among nearby stars are necessary for verification or rejection for the formation mechanism proposed in this paper.

We perform SPH simulations to study precession and changes in alignment between the circumprimary disc and the binary orbit in misaligned binary systems. We find that the precession process can be described by the rigid-disc approximation, where the disc is considered as a rigid body interacting with the binary companion only gravitationally. Precession also causes change in alignment between the rotational axis of the disc and the spin axis of the primary star. This type of alignment is of great important for explaining the origin of spin-orbit misaligned planetary systems. However, we find that the rigid-disc approximation fails to describe changes in alignment between the disc and the binary orbit. This is because the alignment process is a consequence of interactions that involve the fluidity of the disc, such as the tidal interaction and the encounter interaction. Furthermore, simulation results show that there are not only alignment processes, which bring the components towards alignment, but also anti-alignment processes, which tend to misalign the components. The alignment process dominates in systems with misalignment angle near 90 degrees, while the anti-alignment process dominates in systems with the misalignment angle near 0 or 180 degrees. This means that highly misaligned systems will become more aligned but slightly misaligned systems will become more misaligned.

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.

It is widely accepted that stars do not form in isolation but result from the fragmentation of molecular clouds, which in turn leads to star cluster formation. Over time, clusters dissolve or are destroyed by interactions with molecular clouds or tidal stripping, and their members become part of the general field population. Star clusters are thus among the basic building blocks of galaxies. In turn, star cluster populations, from young associations and open clusters to old globulars, are powerful tracers of the formation, assembly, and evolutionary history of their parent galaxies. Although their importance had been recognised for decades, major progress in this area has only become possible in recent years, both for Galactic and extragalactic cluster populations. Star clusters are the observational foundation for stellar astrophysics and evolution, provide essential tracers of galactic structure, and are unique stellar dynamical environments. Star formation, stellar structure, stellar evolution, and stellar nucleosynthesis continue to benefit and improve tremendously from the study of these systems. Additionally, fundamental quantities such as the initial mass function can be successfully derived from modelling either the H-R diagrams or the integrated velocity structures of, respectively, resolved and unresolved clusters and cluster populations. Star cluster studies thus span the fields of Galactic and extragalactic astrophysics, while heavily affecting our detailed understanding of the process of star formation in dense environments.This report highlights science results of the last decade in the major fields covered by IAU Commission 37: Star clusters and associations.

Theory and simulations suggest that it is possible to form low-mass hydrogen-burning stars, brown dwarfs and planetary-mass objects via disc fragmentation. As disc fragmentation results in the formation of several bodies at comparable distances to the host star, their orbits are generally unstable. Here, we study the dynamical evolution of these objects. We set up the initial conditions based on the outcomes of the SPH simulations of Stamatellos & Whitworth, and for comparison we also study the evolution of systems resulting from lower-mass fragmenting discs. We refer to these two sets of simulations as set 1 and set 2. At 10 Myr, approximately half of the host stars have one companion left, and approximately 22% (set 1) to 9.8% (set 2) of the host stars are single. Systems with multiple secondaries in relatively stable configurations are common (about 30% and 44%, respectively). The majority of the companions are ejected within 1 Myr with velocities mostly below 5 km/s, with some runaway escapers with velocities over 30 km/s. About 6% (set 1) and 2% (set 2) of the companions pair up into very low-mass binary systems. The majority of these pairs escape as very low-mass binaries, while others remain bound to the host star in hierarchical configurations (often with retrograde inner orbits). Physical collisions with the host star (0.43 and 0.18 events per host star for set 1 and set 2) and between companions (0.08 and 0.04 events per host star for set 1 and set 2) are relatively common and their frequency increases with increasing disc mass. Our study predicts observable properties of very low-mass binaries, low-mass hierarchical systems, the brown dwarf desert, and free-floating brown dwarfs and planetary-mass objects in and near young stellar groupings, which can be used to distinguish between different formation scenarios of very low-mass stars, brown dwarfs and planetary-mass objects.

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.

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.

A systematic investigation on the evolution of a prolate cloud at an Hii boundary is conducted using Smoothed Particle Hydrodynamics (SPH) in order to understand the mechanism for a variety of irregular morphological structures found at the boundaries of various Hii regions. The prolate molecular clouds in this investigation are set with their semi-major axes at inclinations between 0 and 90 degrees to a plane parallel ionizing radiation flux. A set of 4 parameters, the number density n, the ratio of major to minor axis gamma, the inclination angle phi and the incident flux F_EUV, are used to define the initial state of the simulated clouds. The dependence of the evolution of a prolate cloud under Radiation Driven Implosion (RDI) on each of the four parameters is investigated. It is found that: i) in addition to the well studied standard type A, B or C Bright Rimmed Clouds (BRCs), many other types such as asymmetrical BRCs, filamentary structures and irregular horse-head structures could also be developed at Hii boundaries with only simple initial conditions; ii) the final morphological structures are very sensitive to the 4 initial parameters, especially to the initial density and the inclination; iii) The previously defined ionizing radiation penetration depth can still be used as a good indicator of the final morphology. Based on the simulation results, the efficiency of the RDI triggered star formation from clouds of different initial conditions is also estimated. Finally a unified mechanism for the various morphological structures found in many different Hii boundaries is suggested.

We investigate how different stellar initial mass functions (IMFs) can affect the mass loss and survival of star clusters. We find that IMFs with radically different low-mass cut-offs (between 0.1 and 2 Msun) do not change cluster destruction time-scales as much as might be expected. Unsurprisingly, we find that clusters with more high-mass stars lose relatively more mass through stellar evolution, but the response to this mass loss is to expand and hence significantly slow their dynamical evolution. We also argue that it is very difficult, if not impossible, to have clusters with different IMFs that are initially "the same", since the mass, radius and relaxation times depend on each other and on the IMF in a complex way. We conclude that changing the IMF to be biased towards more massive stars does speed up mass loss and dissolution, but that it is not as dramatic as might be thought.

The evolution of a prolate cloud at an Hii boundary is investigated using Smoothed Particle Hydrodynamics (SPH). The prolate molecular clouds in our investigation are set with their semi-major axis perpendicular to the radiative direction of a plane parallel ionising Extreme Ultraviolet (EUV) flux. Simulations on three high mass prolate clouds reveal that EUV radiation can trigger distinctive high density core formation embedded in a final linear structure. This contrasts with results of the previous work in which only an isotropic Far Ultraviolet (FUV) interstellar background flux was applied. A systematic investigation on a group of prolate clouds of equal mass but different initial densities and geometric shapes finds that the distribution of the cores over the final linear structure changes with the initial conditions of the prolate cloud and the strength of the EUV radiation flux. These highly condensed cores may either scatter over the full length of the final linear structure or form two groups of high density cores at two foci, depending on the value of the ionising radiation penetration depth d_EUV, the ratio of the physical ionising radiation penetration depth to the minor axis of the cloud. Data anlysis on the total mass of the high density cores and the core formation time finds that the potential for EUV radiation triggered star formation efficiency is higher in prolate clouds with shallow ionisation penetration depth and intermediate major to minor axial ratio, for the physical environments investigated. Finally, it is suggested that the various fragment-core structures observed at Hii boundaries may result from the interaction between ionising radiation and pre-existing prolate clouds of different initial geometrical and physical conditions.

Observations from optical to centimeter wavelengths have demonstrated that multiple systems of two or more bodies is the norm at all stellar evolutionary stages. Multiple systems are widely agreed to result from the collapse and fragmentation of cloud cores, despite the inhibiting influence of magnetic fields. Surveys of Class 0 protostars with mm interferometers have revealed a very high multiplicity frequency of about 2/3, even though there are observational difficulties in resolving close protobinaries, thus supporting the possibility that all stars could be born in multiple systems. Near-infrared adaptive optics observations of Class I protostars show a lower binary frequency relative to the Class 0 phase, a declining trend that continues through the Class II/III stages to the field population. This loss of companions is a natural consequence of dynamical interplay in small multiple systems, leading to ejection of members. We discuss observational consequences of this dynamical evolution, and its influence on circumstellar disks, and we review the evolution of circumbinary disks and their role in defining binary mass ratios. Special attention is paid to eclipsing PMS binaries, which allow for observational tests of evolutionary models of early stellar evolution. Many stars are born in clusters and small groups, and we discuss how interactions in dense stellar environments can significantly alter the distribution of binary separations through dissolution of wider binaries. The binaries and multiples we find in the field are the survivors of these internal and external destructive processes, and we provide a detailed overview of the multiplicity statistics of the field, which form a boundary condition for all models of binary evolution. Finally we discuss various formation mechanisms for massive binaries, and the properties of massive trapezia.

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.

Globular clusters (GCs), once thought to be well approximated as simple stellar populations (i.e. all stars having the same age and chemical abundance), are now known to host a variety of anomalies, such as multiple discrete (or spreads in) populations in colour-magnitude diagrams and abundance variations in light elements (e.g., Na, O, Al). Multiple models have been put forward to explain the observed anomalies, although all have serious shortcomings (e.g., requiring a non-standard initial mass function of stars and GCs to have been initially 10-100 times more massive than observed today). These models also do not agree with observations of massive stellar clusters forming today, which do not display significant age spreads nor have gas/dust within the cluster. Here we present a model for the formation of GCs, where low mass pre-main sequence (PMS) stars accrete enriched material released from interacting massive binary and rapidly rotating stars onto their circumstellar discs, and ultimately onto the young stars. As was shown in previous studies, the accreted material matches the unusual abundances and patterns observed in GCs. The proposed model does not require multiple generations of star-formation, conforms to known properties of massive clusters forming today, and solves the "mass budget problem" without requiring GCs to have been significantly more massive at birth. Potential caveats to the model as well as model predictions are discussed.

Whether or not the rich star cluster population in the Large Magellanic Cloud (LMC) is affected by significant disruption during the first few x 10^8 yr of its evolution is an open question and the subject of significant current debate. Here, we revisit the problem, adopting a homogeneous data set of broad-band imaging observations. We base our analysis mainly on two sets of self-consistently determined LMC cluster ages and masses, one using standard modelling and one which takes into account the effects of stochasticity in the clusters' stellar mass functions. On their own, the results based on any of the three complementary analysis approaches applied here are merely indicative of the physical conditions governing the cluster population. However, the combination of our results from all three different diagnostics leaves little room for any conclusion other than that the optically selected LMC star cluster population exhibits no compelling evidence of significant disruption -- for clusters with masses, M_cl, of log(M_cl/M_sun) >= 3.0-3.5 -- between the age ranges of [3-10] Myr and [30-100] Myr, either "infant mortality" or otherwise. In fact, there is no evidence of any destruction beyond that expected from simple models just including stellar dynamics and stellar evolution for ages up to 1 Gyr. It seems, therefore, that the difference in environmental conditions in the Magellanic Clouds on the one hand and significantly more massive galaxies on the other may be the key to understanding the apparent variations in cluster disruption behaviour at early times.

Observations indicate that the central portions of the Present-Day Prestellar Core Mass Function (CMF) and the Stellar Initial Mass Function (IMF) both have approximately log-normal shapes, but that the CMF is displaced to higher mass than the IMF by a factor F = 4+/-1. This has lead to suggestions that the shape of the IMF is directly inherited from the shape of the CMF - and therefore, by implication, that there is a self-similar mapping from the CMF onto the IMF. If we assume a self-similar mapping, it follows (i) that F = N0/eta, where eta is the mean fraction of a core's mass that ends up in stars, and N0 is the mean number of stars spawned by a single core; and (ii) that the stars spawned by a single core must have an approximately log-normal distribution of relative masses, with universal standard deviation sigma0. Observations can be expected to deliver ever more accurate estimates of F, but this still leaves a degeneracy between eta and N0; and sigma0 is also unconstrained by observation. Here we show that these parameters can be estimated by invoking binary statistics. Specifically, if (a) each core spawns one long-lived binary system, and (b) the probability that a star of mass M is part of this long-lived binary is proportional to M^alpha, current observations of the binary frequency as a function of primary mass and the distribution of mass ratios strongly favour eta = 1.0 +/- 0.3, N0 = 4.3 +/- 0.4, sigma0 = 0.3 +/- 0.03 and alpha = 0.9 +/- 0.6. The mapping from CMF to IMF is not necessarily self-similar - there are many possible motivations for a non self-similar mapping - but if it is not, then the shape of the IMF cannot be inherited from the CMF. Given the limited observational constraints currently available and the ability of a self-similar mapping to satisfy them, the possibility that the shape of the IMF is inherited from the CMF cannot be ruled out at this juncture.

We compare the results for a set of hydrodynamical tests performed with the AMR finite volume code, MG and the SPH code, SEREN. The test suite includes shock tube tests, with and without cooling, the non-linear thin-shell instability and the Kelvin-Helmholtz instability. The main conclusions are : (i) the two methods converge in the limit of high resolution and accuracy in most cases. All tests show good agreement when numerical effects (e.g. discontinuities in SPH) are properly treated. (ii) Both methods can capture adiabatic shocks and well-resolved cooling shocks perfectly well with standard prescriptions. However, they both have problems when dealing with under-resolved cooling shocks, or strictly isothermal shocks, at high Mach numbers. The finite volume code only works well at 1st order and even then requires some additional artificial viscosity. SPH requires either a larger value of the artificial viscosity parameter, alpha_AV, or a modified form of the standard artificial viscosity term using the harmonic mean of the density, rather than the arithmetic mean. (iii) Some SPH simulations require larger kernels to increase neighbour number and reduce particle noise in order to achieve agreement with finite volume simulations. However, this is partly due to the need to reduce noise that can corrupt the growth of small-scale perturbations. In contrast, instabilities seeded from large-scale perturbations do not require more neighbours and hence work well with standard SPH formulations and converge with the finite volume simulations. (iv) For purely hydrodynamical problems, SPH simulations take an order of magnitude longer to run than finite volume simulations when running at equivalent resolutions, i.e. when they both resolve the underlying physics to the same degree. This requires about 2-3 times as many particles as the number of cells.

We present a new hybrid Smoothed Particle Hydrodynamics (SPH)/N-body method for modelling the collisional stellar dynamics of young clusters in a live gas background. By deriving the equations of motion from Lagrangian mechanics we obtain a formally conservative combined SPH/N-body scheme. The SPH gas particles are integrated with a 2nd order Leapfrog, and the stars with a 4th order Hermite scheme. Our new approach is intended to bridge the divide between the detailed, but expensive, full hydrodynamical simulations of star formation, and pure N-body simulations of gas-free star clusters. We have implemented this hybrid approach in the SPH code SEREN (Hubber et al. 2011) and perform a series of simple tests to demonstrate the fidelity of the algorithm and its conservation properties. We investigate and present resolution criteria to adequately resolve the density field and to prevent strong numerical scattering effects. Future developments will include a more sophisticated treatment of binaries.

We present deep Hubble Space Telescope/Wide Field and Planetary Camera 2 photometry of the young HD 97950 star cluster in the giant H {\sc ii} region NGC 3603. The data were obtained in 1997 and 2007 permitting us to derive membership based on proper motions of the stars. Our data are consistent with an age of 1 Myr for the HD 97950 cluster. A possible age spread, if present in the cluster, appears to be small. The global slope of the incompleteness-corrected mass function for member stars within 60$"$ is $\rm \Gamma=-0.88\pm0.15$, which is flatter than the value of a Salpeter slope of -1.35. The radially varying mass function shows pronounced mass segregation ranging from slopes of $-0.26 \pm 0.32$ in the inner $5"$ to $-0.94\pm 0.36$ in the outermost annulus ($40"$ -- $60"$). Stars more massive than 50 M$_{\odot}$ are found only in the cluster center. The $\Lambda$ minimum spanning tree technique confirms significant mass segregation down to 30 M$_{\odot}$. The dependence of $\Lambda$ on mass, i.e., that high-mass stars are more segregated than low mass stars, and the (weak) dependence of the velocity dispersion on stellar mass might imply that the mass segregation is dynamical in origin. While primordial segregation cannot be excluded, the properties of the mass segregation indicate that dynamical mass segregation may have been the dominant process for segregation of high-mass stars.

Binary properties are usually expressed (for good observational reasons) as a function of primary mass. It has been found that the distribution of companion masses -- the mass ratio distribution -- is different for different primary masses. We argue that system mass is the more fundamental physical parameter to use. We show that if system masses are drawn from a log-normal mass function, then the different observed mass ratio distributions as a function of primary mass, from M-dwarfs to A-stars, are all consistent with a universal, flat, system mass ratio distribution. We also show that the brown dwarf mass ratio distribution is not drawn from the same flat distribution, suggesting that the process which decides upon mass ratios is very different in brown dwarfs and stars.

We examine the effects of gas expulsion on initially sub-structured and out-of-equilibrium star clusters. We perform N-body simulations of the evolution of star clusters in a static background potential before adjusting that potential to model gas expulsion. We investigate the impact of varying the rate at which the gas is removed, and the instant at which gas removal begins. Reducing the rate at which the gas is expelled results in an increase in cluster survival. Quantitatively, this dependency is approximately in agreement with previous studies, despite their use of smooth, and virialised initial stellar distributions. However, the instant at which gas expulsion occurs is found to have a strong effect on cluster response to gas removal. We find if gas expulsion occurs prior to one crossing time, cluster response is poorly described by any global parameters. Furthermore in real clusters the instant of gas expulsion is poorly constrained. Therefore our results emphasis the highly stochastic and variable response of star clusters to gas expulsion.

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.

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.

Recent discoveries of strongly misaligned transiting exoplanets pose a challenge to the established planet formation theory which assumes planetary systems to form and evolve in isolation. However, the fact that the majority of stars actually do form in star clusters raises the question how isolated forming planetary systems really are. Besides radiative and tidal forces the presence of dense gas aggregates in star-forming regions are potential sources for perturbations to protoplanetary discs or systems. Here we show that subsequent capture of gas from large extended accretion envelopes onto a passing star with a typical circumstellar disc can tilt the disc plane to retrograde orientation, naturally explaining the formation of strongly inclined planetary systems. Furthermore, the inner disc regions may become denser, and thus more prone to speedy coagulation and planet formation. Pre-existing planetary systems are compressed by gas inflows leading to a natural occurrence of close-in misaligned hot Jupiters and short-period eccentric planets. The likelihood of such events mainly depends on the gas content of the cluster and is thus expected to be highest in the youngest star clusters.

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 have investigated the brown dwarf (BD) and stellar disc fractions in the Upper Scorpius OB Association (USco) and compared them with several other young regions. We have compiled the most complete sample of of all spectroscopically confirmed BDs in USco, and have made use of the WISE catalog to identify the disc candidates. We report on the discovery of 12 new BD discs in USco, with spectral type (SpT) between M6 and M8.5. The WISE colors for the new discs are similar to the primordial (transition) discs earlier detected in USco. Combining with previous surveys, we find the lowest inner disc fractions (~20-25%) for a wide range in stellar masses (~0.01-4.0 Msun) in the USco association. The low disc fractions for high-mass stars in USco (and the other clusters) are consistent with an evolutionary decline in inner disc frequency with age. However, BD disc fractions are higher than those for the stars in 1-3 Myr clusters, but very low in the ~5 Myr old USco. Also, primordial BD discs are still visible in the ~10 Myr TW Hydrae association, whereas the higher mass stars have all transitioned to the debris stage by this age. The disc frequencies for BDs do not show any dependence on the stellar density or the BD/star number ratio in a cluster. We suggest that the large differences in the observed BD disc fractions between regions may well be due to different BD formation mechanisms and therefore different initial disc fractions/properties. We also present a WISE SED classification scheme, based on the Ks and WISE bands of 3.4-12micron, to distinguish between the Class I/II and the Class III sequence. Our work includes a comparison of the sensitivities of WISE and Spitzer disc surveys. We estimate that WISE can be incomplete for discs at SpT later than M8 in distant clusters such as SOri. WISE should be able to recover the M8-M9 discs in the nearby young clusters.

We perform N-body simulations of young triple systems consisting of two low-mass companions orbiting around a significantly more massive primary. We find that, when the orbits of the companions are coplanar and not too widely separated, the chance of a collision between the two companions can be as high as 20 per cent. Collisions between one of the companions (always the less massive) and the primary can also occur in systems with unequal-mass companions. The chance of collisions is a few per cent in systems with more realistic initial conditions, such as with slightly non-coplanar orbits and unequal-mass companions. If the companions start widely separated then collision are very rare except in some cases when the total mass of the companions is large. We suggest that collisions between members of young multiple systems may explain some unusual young multiple systems such as apparently non-coeval companions.

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.

Stars form with a complex and highly structured distribution. For a smooth star cluster to form from these initial conditions, the star cluster must erase this substructure. We study how substructure is removed using N-body simulations that realistically handle two-body relaxation. In contrast to previous studies, we find that hierarchical cluster formation occurs chiefly as a result of scattering of stars out of clumps, and not through clump merging. Two-body relaxation, in particular within the body of a clump, can significantly increase the rate at which substructure is erased beyond that of clump-merging alone. Hence the relaxation time of individual clumps is a key parameter controlling the rate at which smooth, spherical star clusters can form. The initial virial ratio of the clumps is an additional key parameter controlling the formation rate of a cluster. Reducing the initial virial ratio causes a star cluster to lose its substructure more rapidly.

We investigate the formation and evolution of high-order massive star multiples similar to the Trapezium in the Orion Nebula Cluster. We perform ensembles of N-body simulations of the evolution of N=1000 Orion-like clusters with initial conditions ranging from cool and clumpy to relatively smooth and relaxed. We find that trapezium-like systems are frequently formed in the first 2 Myr in initially cool and clumpy clusters and can survive for significant amounts of time in such clusters. We also find that these systems are highly dynamical entities, constantly interacting with the surrounding cluster, changing their appearance and membership regularly. The eventual decay of trapezium-like systems can even destroy the host cluster. We argue that the current state of any trapezium-like system is transient and care should be taken when analysing and drawing conclusions from a single snapshot in the life of a highly dynamic object.

We examine the effects of gas expulsion on initially sub-structured and out-of-equilibrium star clusters. We perform $N$-body simulations of the evolution of star clusters in a static background potential before removing that potential to model gas expulsion. We find that the initial star formation efficiency is not a good measure of the survivability of star clusters. This is because the stellar distribution can change significantly, causing a large change in the relative importance of the stellar and gas potentials. We find that the initial stellar distribution and velocity dispersion are far more important parameters than the initial star formation efficiency, and that clusters with very low star formation efficiencies can survive gas expulsion. We suggest that it is variations in cluster initial conditions rather than in their star formation efficiencies that cause some clusters to be destroyed while a few survive.

We recently determined the mass of the most massive star known to the date, R136a1 with a mass at birth 320 times the mass of our sun, as well as the mass of several other stars that are more massive than 150 M. Such massive stars (~150-300 M) may end their life as pair-instability supernovae (PISN) if they retain enough mass until they die. We have calculated a grid of stellar evolution models in order to investigate the impact of mass loss and rotation on the evolution and fate of these very massive stars. As mass loss is very strong at solar metallicity, our models predict that most of the very massive stars will die as type Ic SNe. Only slowly and non-rotating stars at metallicities below that of the LMC might retain enough mass to produce a PISN. This would mean that the first stellar generations might have produced PISN although their chemical signature is not observed in extremely metal poor stars in the halo of our galaxy.

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.

We use high-resolution {\sl Hubble Space Telescope} imaging observations of the young ($\sim 15-25$ Myr-old) star cluster NGC 1818 in the Large Magellanic Cloud to derive an estimate for the binary fraction of F stars ($1.3 < M_\star/M_\odot < 1.6$). This study provides the strongest constraints yet on the binary fraction in a young star cluster in a low-metallicity environment (${[Fe/H]} \sim -0.4$ dex). Employing artificial-star tests, we develop a simple method that can efficiently measure the probabilities of stellar blends and superpositions from the observed stellar catalog. We create synthetic color-magnitude diagrams matching the fundamental parameters of NGC 1818, with different binary fractions and mass-ratio distributions. We find that this method is sensitive to binaries with mass ratios, $q \ga 0.4$. For binaries with F-star primaries and mass ratios $q > 0.4$, the binary fraction is $\sim 0.35$. This suggests a total binary fraction for F stars of 0.55 to unity, depending on assumptions about the form of the mass-ratio distribution at low $q$.

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).

Most stars are born in clusters and the resulting gravitational interactions between cluster members may significantly affect the evolution of circumstellar discs and therefore the formation of planets and brown dwarfs. Recent findings suggest that tidal perturbations of typical circumstellar discs due to close encounters may inhibit rather than trigger disc fragmentation and so would seem to rule out planet formation by external tidal stimuli. However, the disc models in these calculations were restricted to disc radii of 40 AU and disc masses below 0.1 M_sun. Here we show that even modest encounters can trigger fragmentation around 100 AU in the sorts of massive (~0.5 M_sun), extended (>=100 AU) discs that are observed around young stars. Tidal perturbation alone can do this, no disc-disc collision is required. We also show that very-low-mass binary systems can form through the interaction of objects in the disc. In our computations, otherwise non-fragmenting massive discs, once perturbed, fragment into several objects between about 0.01 and 0.1 M_sun, i.e., over the whole brown dwarf mass range. Typically these orbit on highly eccentric orbits or are even ejected. While probably not suitable for the formation of Jupiter- or Neptune-type planets, our scenario provides a possible formation mechanism for brown dwarfs and very massive planets which, interestingly, leads to a mass distribution consistent with the canonical substellar IMF. As a minor outcome, a possible explanation for the origin of misaligned extrasolar planetary systems is discussed.

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.

We present B and R band spectroastrometry of a sample of 45 Herbig Ae/Be stars in order to study their binary properties. All but one of the targets known to be binary systems with a separation of ~0.1-2.0 arcsec are detected by a distinctive spectroastrometric signature. Some objects in the sample exhibit spectroastrometric features that do not appear attributable to a binary system. We find that these may be due to light reflected from dusty halos or material entrained in winds. We present 8 new binary detections and 4 detections of an unknown component in previously discovered binary systems. The data confirm previous reports that Herbig Ae/Be stars have a high binary fraction, 74+/-6 per cent in the sample presented here. We use a spectroastrometric deconvolution technique to separate the spatially unresolved binary spectra into the individual constituent spectra. The separated spectra allow us to ascertain the spectral type of the individual binary components, which in turn allows the mass ratio of these systems to be determined. In addition, we appraise the method used and the effects of contaminant sources of flux. We find that the distribution of system mass ratios is inconsistent with random pairing from the Initial Mass Function, and that this appears robust despite a detection bias. Instead, the mass ratio distribution is broadly consistent with the scenario of binary formation via disk fragmentation.

Many, possibly most, stars form in binary and higher-order multiple systems. Therefore, the properties and frequency of binary systems provide strong clues to the star-formation process, and constraints on star-formation models. However, the majority of stars also form in star clusters in which the birth binary properties and frequency can be altered rapidly by dynamical processing. Thus, we almost never see the birth population, which makes it very difficult to know if star formation (as traced by binaries, at least) is universal, or if it depends on environment. In addition, the field population consists of a mixture of systems from different clusters which have all been processed in different ways.

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.

Obtaining accurate measurements of the initial mass function (IMF) is often considered to be the key to understanding star formation, and a universal IMF is often assumed to imply a universal star formation process. Here, we illustrate that different modes of star formation can result in the same IMF, and that, in order to truly understand star formation, a deeper understanding of the primordial binary population is necessary. Detailed knowledge on the binary fraction, mass ratio distribution, and other binary parameters, as a function of mass, is a requirement for recovering the star formation process from stellar population measurements.

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 show that the same initial mass function (IMF) can result from very different modes of star formation from very similar underlying core and/or system mass functions. In particular, we show that the canonical IMF can be recovered from very similar system mass functions, but with very different mass ratio distributions within those systems. This is a consequence of the basically log-normal shapes of all of the distributions. We also show that the relationships between the shapes of the core, system, and stellar mass functions may not be trivial. Therefore, different star formation in different regions could still result in the same IMF.

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.

We perform SPH simulations of the collapse and fragmentation of low-mass cores having different initial levels of turbulence (alpha_turb=0.05,0.10,0.25). We use a new treatment of the energy equation which captures the transport of cooling radiation against opacity due to both dust and gas (including the effects of dust sublimation, molecules, and H^- ions). We also perform comparison simulations using a standard barotropic equation of state. We find that -- when compared with the barotropic equation of state -- our more realistic treatment of the energy equation results in more protostellar objects being formed, and a higher proportion of brown dwarfs; the multiplicity frequency is essentially unchanged, but the multiple systems tend to have shorter periods (by a factor ~3), higher eccentricities, and higher mass ratios. The reason for this is that small fragments are able to cool more effectively with the new treatment, as compared with the barotropic equation of state. We find that the process of fragmentation is often bimodal. The first protostar to form is usually, at the end, the most massive, i.e. the primary. However, frequently a disc-like structure subsequently forms round this primary, and then, once it has accumulated sufficient mass, quickly fragments to produce several secondaries. We believe that this delayed fragmentation of a disc-like structure is likely to be an important source of very low-mass hydrogen-burning stars and brown dwarfs.

Knowledge of the binary population in stellar groupings provides important information about the outcome of the star forming process in different environments. Binarity is also a key ingredient in stellar population studies and is a prerequisite to calibrate the binary evolution channels. In these proceedings we present an overview of several commonly used methods to pair individual stars into binary systems, which we refer to as the pairing function. Many pairing functions are frequently used by observers and computational astronomers, either for the mathematical convenience, or because they roughly describe the expected outcome of the star forming process. We discuss the consequences of each pairing function for the interpretation of observations and numerical simulations. The binary fraction and mass ratio distribution generally depend strongly on the selection of the range in primary spectral type in a sample. These quantities, when derived from a binary survey with a mass-limited sample of target stars, are thus not representative for the population as a whole.

Knowledge of the binary population in stellar groupings provides important information about the outcome of the star forming process in different environments (see, e.g., Blaauw 1991, and references therein). Binarity is also a key ingredient in stellar population studies, and is a prerequisite to calibrate the binary evolution channels. In this paper we present an overview of several commonly used methods to pair individual stars into binary systems, which we refer to as pairing functions. These pairing functions are frequently used by observers and computational astronomers, either for their mathematical convenience, or because they roughly describe the expected outcome of the star forming process. We discuss the consequences of each pairing function for the interpretation of observations and numerical simulations. The binary fraction and mass ratio distribution generally depend strongly on the selection of the range in primary spectral type in a sample. The mass ratio distribution and binary fraction derived from a binarity survey among a mass-limited sample of targets is thus not representative for the population as a whole. Neither theory nor observations indicate that random pairing of binary components from the mass distribution, the simplest pairing function, is realistic. It is more likely that companion stars are formed in a disk around a star, or that a pre-binary core fragments into two binary components. The results of our analysis are important for (i) the interpretation of the observed mass ratio distribution and binary fraction for a sample of stars, (ii) a range of possible initial condition algorithms for star cluster simulations, and (iii) how to discriminate between the different star formation scenarios.

The diagnostic age versus mass-to-light ratio diagram is often used in attempts to constrain the shape of the stellar initial mass function (IMF), and the potential longevity of extragalactic young to intermediate-age massive star clusters. Here, we explore its potential for Galactic open clusters. On the basis of a small, homogenised cluster sample we provide useful constraints on the presence of significant binary fractions. Using the massive young Galactic cluster Westerlund 1 as a key example, we caution that stochasticity in the IMF introduces significant additional uncertainties. We conclude that for an open cluster to survive for any significant length of time, and in the absence of substantial external perturbations, it is a necessary but not a sufficient condition to be located close to or (in the presence of a significant binary population) somewhat BELOW the predicted photometric evolutionary sequences for "normal" simple stellar populations (although such a location may be dominated by a remaining "bound" cluster core and thus not adequately reflect the overall cluster dynamics).

The diagnostic age versus mass-to-light ratio diagram is often used in attempts to constrain the shape of the stellar initial mass function, and the stability and the potential longevity of extragalactic young to intermediate-age massive star clusters. Here, we explore the pitfalls associated with this approach and its potential for use with Galactic open clusters. We conclude that for an open cluster to survive for any significant fraction of a Hubble time (in the absence of substantial external perturbations), it is a necessary but not a sufficient condition to be located close to the predicted photometric evolutionary sequences for "normal" simple stellar populations.

The evolution of star clusters in the Magellanic Clouds has been the subject of significant recent controversy, particularly regarding the importance and length of the earliest, largely mass-independent disruption phase (referred to as "infant mortality"). Here, we take a fresh approach to the problem, using a large, independent, and homogeneous data set of UBVR imaging observations, from which we obtain the cluster age and mass distributions in both the Large and Small Magelanic Clouds (LMC, SMC) in a self-consistent manner. We conclude that the (optically selected) SMC star cluster population has undergone at most ~30% (1sigma) infant mortality between the age range from about 3-10 Myr, to that of approximately 40-160 Myr. We rule out a 90% cluster mortality rate per decade of age (for the full age range up to 10^9 yr) at a >6sigma level. Using a simple approach, we derive a "characteristic" cluster disruption time-scale for the cluster population in the LMC that implies that we are observing the INITIAL cluster mass function. Preliminary results suggest that the LMC cluster population may be affected by <10% infant mortality.

The observed properties of young star clusters, such as the core radius and luminosity profile, change rapidly during the early evolution of the clusters. Here we present observations of 6 young clusters in M51 where we derive their sizes using HST imaging and ages using deep Gemini-North spectroscopy. We find evidence for a rapid expansion of the cluster cores during the first 20 Myr of their evolution. We confirm this trend by including data from the literature of both Galactic and extra-galactic embedded and young clusters, and possible mechanisms (rapid gas removal, stellar evolutionary mass-loss, and internal dynamical heating) are discussed. We explore the implications of this result, focussing on the fact that clusters were more concentrated in the past, implying that their stellar densities were much higher and relaxation times correspondingly shorter. Thus, when estimating if a particular cluster is dynamically relaxed, (i.e. when determining if a cluster's mass segregation is due to primordial or dynamical processes), the current relaxation time is only an upper-limit, with the relaxation time likely being significantly shorter in the past.

The star formation efficiency (SFE) of a star cluster is thought to be the critical factor in determining if the cluster can survive for a significant (>50 Myr) time. There is an often quoted critical SFE of ~30 per cent for a cluster to survive gas expulsion. I reiterate that the SFE is not the critical factor, rather it is the dynamical state of the stars (as measured by their virial ratio) immediately before gas expulsion that is the critical factor. If the stars in a star cluster are born in an even slightly cold dynamical state then the survivability of a cluster can be greatly increased.

We explore how the expulsion of gas from star-cluster forming cloud-cores due to supernova explosions affects the shape of the initial cluster mass function, that is, the mass function of star clusters when effects of gas expulsion are over. We demonstrate that if the radii of cluster-forming gas cores are roughly constant over the core mass range, as supported by observations, then more massive cores undergo slower gas expulsion. Therefore, for a given star formation efficiency, more massive cores retain a larger fraction of stars after gas expulsion. The initial cluster mass function may thus differ from the core mass function substantially, with the final shape depending on the star formation efficiency. A mass-independent star formation efficiency of about 20 per cent turns a power-law core mass function into a bell-shaped initial cluster mass function, while mass-independent efficiencies of order 40 per cent preserve the shape of the core mass function.

Stars form from dense molecular cores, and the mass function of these cores (the CMF) is often found to be similar to the form of the stellar initial mass function (IMF). This suggests that the form of the IMF is the result of the form of the CMF. However, most stars are thought to form in binary and multiple systems, therefore the relationship between the IMF and the CMF cannot be trivial. We test two star formation scenarios - one in which all stars form as binary or triple systems, and one in which low-mass stars form in a predominantly single mode. We show that from a log-normal CMF, similar to those observed, and expected on theoretical grounds, the model in which all stars form as multiples gives a better fit to the IMF.

We have analysed near-infrared NTT/SofI observations of the starburst cluster Westerlund 1, which is among the most massive young clusters in the Milky Way. A comparison of colour-magnitude diagrams with theoretical main-sequence and pre-main sequence evolutionary tracks yields improved extinction and distance estimates of A_Ks = 1.13+-0.03 mag and d = 3.55+-0.17 kpc (DM = 12.75+-0.10 mag). The pre-main sequence population is best fit by a Palla & Stahler isochrone for an age of 3.2 Myr, while the main sequence population is in agreement with a cluster age of 3 to 5 Myr. An analysis of the structural parameters of the cluster yields that the half-mass radius of the cluster population increases towards lower mass, indicative of the presence of mass segregation. The cluster is clearly elongated with an eccentricity of 0.20 for stars with masses between 10 and 32 Msun, and 0.15 for stars with masses in the range 3 to 10 Msun. We derive the slope of the stellar mass function for stars with masses between 3.4 and 27 Msun. In an annulus with radii between 0.75 and 1.5 pc from the cluster centre, we obtain a slope of Gamma = -1.3. Closer in, the mass function of Westerlund 1 is shallower with Gamma = -0.6. The extrapolation of the mass function for stars with masses from 0.08 to 120 Msun yields an initial total stellar mass of ~52,000 Msun, and a present-day mass of 20,000 to 45,000 Msun (about 10 times the stellar mass of the Orion Nebula Cluster, and 2 to 4 times the mass of the NGC 3603 young cluster), indicating that Westerlund 1 is the most massive starburst cluster identified to date in the Milky Way.

The early evolution of star clusters in the Small Magellanic Cloud (SMC) has been the subject of significant recent controversy, particularly regarding the importance and length of the earliest, largely mass-independent disruption phase (referred to as "infant mortality"). Here, we take a fresh approach to the problem, using an independent, homogeneous data set of UBVR imaging observations, from which we obtain the SMC's cluster age and mass distributions in a self-consistent manner. We conclude that the (optically selected) SMC star cluster population has undergone at most ~30 per cent (1 sigma) infant mortality between the age range from about (3-10) Myr, to that of approximately (40-160) Myr. We rule out a 90 per cent cluster mortality rate per decade of age (for the full age range up to 10^9 yr) at a > 6 sigma level. We independently affirm this scenario based on the age distribution of the SMC cluster sample.

We introduce and test a new and highly efficient method for treating the thermal and radiative effects influencing the energy equation in SPH simulations of star formation. The method uses the density, temperature and gravitational potential of each particle to estimate a mean optical depth, which then regulates the particle's heating and cooling. The method captures -- at minimal computational cost -- the effects of (i) the rotational and vibrational degrees of freedom of H2, H2 dissociation, H0 ionisation, (ii) opacity changes due to ice mantle melting, sublimation of dust, molecular lines, H-, bound-free and free-free processes and electron scattering; (iv) external irradiation; and (v) thermal inertia. The new algorithm reproduces the results of previous authors and/or known analytic solutions. The computational cost is comparable to a standard SPH simulation with a simple barotropic equation of state. The method is easy to implement, can be applied to both particle- and grid-based codes, and handles optical depths 0<tau<10^{11}.

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.

Context: The principal mechanism by which brown dwarfs form, and its relation to the formation of higher-mass (i.e. hydrogen-burning) stars, is poorly understood. Aims: We advocate a new model for the formation of brown dwarfs. Methods: In this model, brown dwarfs are initially binary companions, formed by gravitational fragmentation of the outer parts (R > 100 AU) of protostellar discs around low-mass hydrogen-burning stars. Most of these binaries are then gently disrupted by passing stars to create a largely single population of brown dwarfs and low-mass hydrogen-burning stars. Results: This idea is consistent with the excess of binaries found in low-density pre-main sequence populations, like that in Taurus, where they should survive longer than in denser clusters. Conclusions: If brown dwarfs form in this way, as companions to more massive stars, the difficulty of forming very low-mass prestellar cores is avoided. Since the disrupted binaries will tend to be those involving low-mass components and wide orbits, and since disruption will be due to the gentle tides of passing stars (rather than violent N-body interactions in small-N sub-clusters), the liberated brown dwarfs will have velocity dispersions and spatial distributions very similar to higher-mass stars, and they will be able to retain discs, and thereby to sustain accretion and outflows. Thus the problems associated with the ejection and turbulence mechanisms can be avoided. This model implies that most, possibly all, stars and brown dwarfs form in binary or multiple systems.

Context: There is a need to improve the fidelity of SPH simulations of self-gravitating gas dynamics. Aims: We remind users of SPH that, if smoothing lengths are adjusted so as to keep the number of neighbours, ${\cal N}$, in the range ${\cal N}_{_{\rm NEIB}}\pm\Delta{\cal N}_{_{\rm NEIB}}$, the tolerance, $\Delta{\cal N}_{_{\rm NEIB}}$, should be set to zero, as first noted by Nelson & Papaloizou. We point out that this is a very straightforward and computationally inexpensive constraint to implement. Methods: We demonstrate this by simulating acoustic oscillations of a self-gravitating isentropic monatomic gas-sphere (cf. Lucy), using ${\cal N}_{_{\rm TOT}}\sim6,000$ particles and ${\cal N}_{_{\rm NEIB}}=50$. Results: We show that there is a marked reduction in the rates of numerical dissipation and diffusion as $\Delta{\cal N}_{_{\rm NEIB}}$ is reduced from 10 to zero. Moreover this reduction incurs a very small computational overhead. Conclusions: We propose that this should become a standard test for codes used in simulating star formation. It is a highly relevant test, because pressure waves generated by the switch from approximate isothermality to approximate adiabaticity play a critical role in the fragmentation of collapsing prestellar cores. Since many SPH simulations in the literature use ${\cal N}_{_{\rm NEIB}}=50$ and $\Delta{\cal N}_{_{\rm NEIB}}\geq10$, their results must be viewed with caution.

We study, analytically and numerically, the energy input into dark matter mini-haloes by interactions with stars. We find that the fractional energy input in simulations of Plummer spheres agrees well with the impulse approximation for small and large impact parameters, with a rapid transition between these two regimes. Using the impulse approximation the fractional energy input at large impact parameters is fairly independent of the mass and density profile of the mini-halo, however low-mass mini-haloes experience a greater fractional energy input in close encounters. We formulate a fitting function which encodes these results and use it to estimate the disruption timescales of mini-haloes, taking into account the stellar velocity dispersion and mass distribution. For mini-haloes with mass M< {\cal O}(10^{-7} M_{\odot} on typical orbits which pass through the disc, we find that the estimated disruption timescales are independent of mini-halo mass, and are of order the age of the Milky Way. For more massive mini-haloes the estimated disruption timescales increase rapidly with increasing mass.

We examine the luminosity and dynamical mass estimates for young massive stellar clusters. For many young (<50 Myr) clusters, the luminosity and dynamical mass estimates differ by a significant amount. We explain this as being due to many young clusters being out of virial equilibrium (which is assumed in dynamical mass estimates) because the clusters are undergoing violent relaxation after expelling gas not used in star formation. We show that, if we assume that luminous mass estimates are correct (for a standard IMF), at least 50 per cent of young clusters for which dynamical masses are known are likely to be destroyed within a few 10s Myr of their formation. Even clusters which will retain a bound core may lose a large fraction of their stellar mass. We also show that the core radius and other structural parameters change significantly during the violent relaxation that follows gas expulsion and that they should be considered instantaneous values only, not necessarily reflecting the final state of the cluster. In particular we note that the increasing core radii observed in young LMC/SMC clusters can be well explained as an effect of rapid gas loss.

Many recent works have attempted to constrain the stellar initial mass function (IMF) inside massive clusters by comparing their dynamical mass estimates to the measured light. These studies have come to different conclusions, with some claiming standard Kroupa-type IMFs, while others have claimed extreme non-standard IMFs. However, the results appear to be correlated with the age of the clusters, as older clusters (>80 Myr) all appear to be well fit by a Kroupa-type IMF whereas younger clusters display significant scatter in their best fitting IMF. Here we show that this is likely due to the fact that young clusters are out of virial equilibrium and therefore cannot be used for such studies. Hence only the older clusters are suitable for IMF studies. Using only these clusters we find that the IMF does not vary significantly. The youngest clusters can be used instead to constrain the star-formation efficiency (SFE) within clusters. We find that the SFE varies between 20 and 60% and we conclude that approximately 60% of young clusters are unbound and will not survive for more than a few 10's of Myr (i.e. "infant mortality").

We present detailed luminosity profiles of the young massive clusters M82-F, NGC 1569-A, and NGC 1705-1 which show significant departures from equilibrium (King and EFF) profiles. We compare these profiles with those from N-body simulations of clusters which have undergone the rapid removal of a significant fraction of their mass due to gas expulsion. We show that the observations and simulations agree very well with each other suggesting that these young clusters are undergoing violent relaxation and are also losing a significant fraction of their stellar mass. That these clusters are not in equilibrium can explain the discrepant mass-to-light ratios observed in many young clusters with respect to simple stellar population models without resorting to non-standard initial stellar mass functions as claimed for M82-F and NGC 1705-1. We also discuss the effect of rapid gas removal on the complete disruption of a large fraction of young massive clusters (``infant mortality''). Finally we note that even bound clusters may lose >50% of their initial stellar mass due to rapid gas loss (``infant weight-loss'').

Almost all young stars are found in multiple systems. This suggests that protostellar cores almost always fragment into multiple objects. The observed properties of multiple systems such as their separation distribution and mass ratios provide strong constraints on star formation theories. We review the observed properties of young and old multiple systems and find that the multiplicity of stars changes. Such an evolution is probably due to (a) the dynamical decay of small-N systems and/or (b) the destruction of multiple systems within dense clusters. We review simulations of the fragmentation of rotating and turbulent molecular cores. Such models almost always produce multiple systems, however the properties of those systems do not match observations at all well. Magnetic fields appear to supress fragmentation, prehaps suggesting that they are not dynamically important in the formation of multiple systems. We finish by discussing possible reasons why theory fails to match observation, and the future prospects for this field.

We investigate the effect of the turbulent power spectrum (P(k) \propto k^{-n}, with n=3, 4 or 5) on the fragmentation of low-mass cores, by means of SPH simulations. We adopt initial density profiles and low levels of turbulence based on observation, and for each n-value we conduct an ensemble of simulations with different initial seeds for the turbulent velocity field, so as to obtain reasonable statistics. We find that when power is concentrated at larger scales (i.e. for larger n), more protostellar objects form and there is a higher proportion of low-mass stars and brown dwarfs. This is in direct contrast with the recent results of Delgado Donate et al., presumably because they adopted much higher levels of turbulence.

Lada et al. recently presented a detailed near-infrared extinction map of Globule G2 in the Coalsack molecular cloud complex, showing that this starless core has a well-defined central extinction minimum. We propose a model for G2 in which a rapid increase in external pressure is driving an approximately symmetric compression wave into the core. The rapid increase in external pressure could arise because the core has recently been assimilated by the Coalsack cloud complex, or because the Coalsack has recently been created by two large-scale converging flows. The resulting compression wave has not yet converged on the centre of the core, so there is a central rarefaction. The compression wave has increased the density in the swept-up gas by about a factor of ten, and accelerated it inwards to speeds of order $0.4 {\rm km} {\rm s}^{-1}$. It is shown that even small levels of initial turbulence destroy the ring seen in projection almost completely. In the scenario of strong external compression that we are proposing this implies that the initial turbulent energy in this globule is such that $E_{{\rm turb}} / E_{{\rm grav}} \le 2 %$. Protostar formation should occur in about $40,000 {\rm years}$.

We examine the spatial distribution of brown dwarfs produced by the decay of small-N stellar systems as expected from the embryo ejection scenario. We model a cluster of several hundred stars grouped into 'cores' of a few stars/brown dwarfs. These cores decay, preferentially ejecting their lowest-mass members. Brown dwarfs are found to have a wider spatial distribution than stars, however once the effects of limited survey areas and unresolved binaries are taken into account it can be difficult to distinguish between clusters with many or no ejections. A large difference between the distributions probably indicates that ejections have occurred, however similar distributions sometimes arise even with ejections. Thus the spatial distribution of brown dwarfs is not necessarily a good discriminator between ejection and non-ejection scenarios.

We review four mechanisms for forming brown dwarfs: (i) turbulent fragmentation (producing very low-mass prestellar cores); (ii) gravitational instabilities in discs; (iii) dynamical ejection of stellar embryos from their placental cores; and (iv) photo-erosion of pre-existing cores in HII regions. We argue (a) that these are simply the mechanisms of low-mass star formation, and (b) that they are not mutually exclusive. If, as seems possible, all four mechanisms operate in nature, their relative importance may eventually be constrained by their ability to reproduce the binary statistics of brown dwarfs, but this will require fully 3-D radiative magneto-hydrodynamic simulations.

Jeans showed analytically that, in an infinite uniform-density isothermal gas, plane-wave perturbations collapse to dense sheets if their wavelength, $\lambda$, satisfies $\lambda > \lambda_{_{\rm JEANS}} = (\pi a^2 / G \rho_{_0})^{1/2}$ (where $a$ is the isothermal sound speed and $\rho_{_0}$ is the unperturbed density); in contrast, perturbations with smaller $\lambda$ oscillate about the uniform density state. Here we show that Smoothed Particle Hydrodynamics reproduces these results well, even when the diameters of the SPH particles are twice the wavelength of the perturbation. Our simulations are performed in 3-D with initially settled (i.e. non-crystalline) distributions of particles. Therefore there exists the seed noise for artificial fragmentation, but it does not occur. We conclude that, although there may be -- as with any numerical scheme -- `skeletons in the SPH cupboard', a propensity to fragment artificially is evidently not one of them.

Class 0 objects, which are thought to be the youngest protostars, are identified in terms of NIR or radio emission and/or the presence of molecular outflows. We present combined hydrodynamic and radiative transfer simulations of the collapse of a star-forming molecular core, which suggest two criteria for identifying dense cores with deeply embedded very young protostars that may not be observable in the NIR or radio with current telescopes. We find that cores with protostars are relatively warm (T>15K) with their SEDs peaking at wavelengths <170 micron, and they tend to appear circular.

Most stars - especially young stars - are observed to be in multiple systems. Dynamical evolution is unable to pair stars efficiently, which leads to the conclusion that star-forming cores must usually fragment into \geq 2 stars. However, the dynamical decay of systems with \geq 3 or 4 stars would result in a large single-star population that is not seen in the young stellar population. Additionally, ejections would produce a significant population of hard binaries that are not observed. This leads to a strong constraint on star formation theories that cores must typically produce only 2 or 3 stars. This conclusion is in sharp disagreement with the results of currently available numerical simulations that follow the fragmentation of molecular cores and typically predict the formation of 5--10 seeds per core. In addition, open cluster remnants may account for the majority of observed highly hierarchical higher-order multiple systems in the field.

We study the transition from a prestellar core to a Class 0 protostar, using SPH to simulate the dynamical evolution, and a Monte Carlo radiative transfer code to generate the SED and isophotal maps. For a prestellar core illuminated by the standard interstellar radiation field, the luminosity is low and the SED peaks at ~190 micron. Once a protostar has formed, the luminosity rises (due to a growing contribution from accretion onto the protostar) and the peak of the SED shifts to shorter wavelengths (~80-100 micron). However, by the end of the Class 0 phase, the accretion rate is falling, the luminosity has decreased, and the peak of the SED shifts back towards longer wavelengths (90-150 micron). In our simulations, the density of material around the protostar remains sufficiently high well into the Class 0 phase that the protostar only becomes visible in the NIR if it is displaced from the centre dynamically. Raw submm/mm maps of Class 0 protostars tend to be much more centrally condensed than those of prestellar cores. However, when convolved with a typical telescope beam, the difference in central concentration is less marked, although the Class 0 protostars appear more circular. Our results suggest that, if a core is deemed to be prestellar on the basis of having no associated IRAS source, no cm radio emission, and no outflow, but it has a circular appearance and an SED which peaks at wavelengths below ~170 micron, it may well contain a very young Class 0 protostar.

We present new spectroscopic and photometric observations of the young Galactic open cluster Westerlund 1 (Wd 1) that reveal a unique population of massive evolved stars. We identify ~200 cluster members and present spectroscopic classifications for ~25% of these. We find that all stars so classified are unambiguously post-Main Sequence objects, consistent with an apparent lack of an identifiable Main Sequence in our photometric data to V~20. We are able to identify rich populations of Wolf Rayet stars, OB supergiants and short lived transitional objects. Of these, the latter group consists of both hot (Luminous Blue Variable and extreme B supergiants) and cool (Yellow Hypergiant and Red Supergiant) objects - we find that half the known Galactic population of YHGs resides within Wd1. We obtain a mean V-M_V ~25 mag from the cluster Yellow Hypergiants, implying a Main Sequence turnoff at or below M_V =-5 (O7 V or later). Based solely on the masses inferred for the 53 spectroscopically classified stars, we determine an absolute minimum mass of \~1.5 x 10^3 Msun for Wd 1. However, considering the complete photometrically and spectroscopically selected cluster population and adopting a Kroupa IMF we infer a likely mass for Wd 1 of ~10^5 Msun, noting that inevitable source confusion and incompleteness are likely to render this an underestimate. As such, Wd 1 is the most massive compact young cluster yet identified in the Local Group, with a mass exceeding that of Galactic Centre clusters such as the Arches and Quintuplet. Indeed, the luminosity, inferred mass and compact nature of Wd 1 are comparable with those of Super Star Clusters - previously identified only in external galaxies - and is consistent with expectations for a Globular Cluster progenitor.

We review and evaluate four mechanisms for forming brown dwarves: (i) dynamical ejection of a stellar embryo from its placental prestellar core; (ii) opacity-limited fragmentation of a shock-compressed layer; (iii) gravitational instabilities in discs, triggered by impulsive interactions with other discs or naked stars; and (iv) photo-erosion of pre-existing cores. All these mechanisms can produce free-floating brown dwarves, but only (ii) and (iii) are likely to produce brown dwarves in multiple systems, and (i) has difficulty delivering brown dwarves with discs.

Changing the form of the star cluster mass function (CMF) can effectively change the upper end of the stellar initial mass function. The yields of from supernovae are very sensitive to the mass of the progenitor star. We show that by changing the parameters of the CMF, it is possible to change the yields of oxygen and magnesium by a factor of $\sim 1.5$ and of metals in general by a factor of $\sim 1.8$.

The non-dissipative relativistic force-free condition should be a good approximation to describe the electromagnetic field in much of the pulsar magnetosphere, but we may plausibly expect it to break down in singular domains. Self-consistent magnetospheric solutions are found with field lines closing both at and within the light-cylinder. In general, the detailed properties of the solutions may be affected critically by the physics determining the appropriate choice of equatorial boundary condition beyond the light-cylinder.

We examine the initial mass functions (IMFs) of stars produced by different molecular core mass functions. Simulations suggest that more massive cores produce more stars, so we propose a model in which the average number of stars formed in a core is equal to the initial number of Jeans masses in that core. Small-$N$ systems decay through dynamical interactions, ejecting low-mass stars and brown dwarfs which populate the low-mass tail of the IMF. Stars which remain in cores are able to competitively accrete more gas and become more massive. We deduce the forms of the core mass functions required to explain the IMFs of Taurus, Orion, IC 348 and NGC 2547. These core mass functions fall into two categories - one which peaks at a few $M_\odot$ to explain Taurus and NGC 2547, and one that peaks at around $0.2 M_\odot$ to explain Orion and IC 348.

(Abridged) We explore, by means of a large ensemble of SPH simulations, how the level of turbulence affects the collapse and fragmentation of a star-forming core. All our simulated cores have the same, except that we vary (a) the initial level of turbulence (as measured by the ratio of turbulent to gravitational energy, $\alpha_{\rm turb} \equiv U_{\rm turb}/|\Omega| = 0, 0.01, 0.025, 0.05, 0.10 {\rm and} 0.25$) and (b), for fixed $\alpha_{\rm turb}$, the details of the initial turbulent velocity field (so as to obtain good statistics). A low level of turbulence ($\alpha_{\rm turb} \sim 0.05$) suffices to produce multiple systems. As $\alpha_{\rm turb}$ is increased, the number of objects formed and the companion frequency both increase. The mass function is bimodal, with a flat low-mass segment representing single objects ejected from the core before they can accrete much, and a Gaussian high-mass segment representing objects which because they remain in the core grow by accretion and tend to pair up in multiple systems.

In comparison with other well studied star formation regions, Taurus is unusual in several respects. (i) Its stellar initial mass function (IMF) peaks at relatively high mass ($\sim 0.8 M_\odot$), but contains very few stars much more massive than $1 M_\odot$, and is relatively deficient in brown dwarfs. (ii) It has a higher binary fraction, particularly at large separations. (iii) Its core mass function is strongly peaked at a few $M_\odot$, and the cores have extended envelopes and relatively low levels of turbulence. We present here the results of an ensemble of hydrodynamic simulations which suggest that the unusual stellar IMF in Taurus is a direct consequence of the unusual properties of the cores there. By following the collapse and fragmentation of cores having properties typical of Taurus, we find that roughly 50% of the objects formed in a core, predominantly the low-mass ones, are ejected from the core to form a population of low-mass stars and brown dwarfs with a flat mass function. The remaining objects form multiple systems within the core, accreting until their masses approach $1 M_\odot$; this produces a population of intermediate-mass stars whose mass function peaks at $\sim 0.8 M_\odot$. Together these two populations reproduce the IMF in Taurus very well. This demonstrates, for the first time, a direct causal link between the core mass function and the stellar IMF in a star formation region.

We investigate numerically and semi-analytically the collapse of low-mass, rotating prestellar cores. Initially, the cores are in approximate equilibrium with low rotation (the initial ratio of thermal to gravitational energy is $\alpha_0 \simeq 0.5$, and the initial ratio of rotational to gravitational energy is $\beta_0 = 0.02 {\rm to} 0.05$). They are then subjected to a steady increase in external pressure. Fragmentation is promoted -- in the sense that more protostars are formed -- both by more rapid compression, and by higher rotation (larger $\beta_0$). In general, the large-scale collapse is non-homologous, and follows the pattern described in Paper I for non-rotating clouds, viz. a compression wave is driven into the cloud, thereby increasing the density and the inflow velocity. The effects of rotation become important at the centre, where the material with low angular momentum forms a central primary protostar (CPP), whilst the material with higher angular momentum forms an accretion disc around the CPP. More rapid compression drives a stronger compression wave and delivers material more rapidly into the outer parts of the disc.

We simulate the dynamics of fractal star clusters, in order to investigate the evolution of substructure in recently formed clusters. The velocity dispersion is found to be the key parameter determining the survival of substructure. In clusters with a low initial velocity dispersion, the ensuing collapse of the cluster tends to erase substructure, although some substructure may persist beyond the collapse phase. In clusters with virial ratios of 0.5 or higher, initial density substructure survives for several crossing times, in virtually all cases. Even an initially homogeneous cluster can develop substructure, if it is born with coherent velocity dispersion. These results suggest that the simple initial conditions used for many sophisticated $N$-body simulations could be missing a very important and dramatic phase of star cluster evolution.

We present the results of an ensemble of simulations of the collapse and fragmentation of dense star-forming cores. We show that even with very low levels of turbulence the outcome is usually a binary, or higher-order multiple, system. We take as the initial conditions for these simulations a typical low-mass core, based on the average properties of a large sample of observed cores. All the simulated cores start with a mass of $M = 5.4 M_{\odot}$, a flattened central density profile, a ratio of thermal to gravitational energy $\alpha_{\rm therm} = 0.45$ and a ratio of turbulent to gravitational energy $\alpha_{\rm turb} = 0.05 $. Even this low level of turbulence is sufficient to produce multiple star formation in 80% of the cores; the mean number of stars and brown dwarfs formed from a single core is 4.55, and the maximum is 10. At the outset, the cores have no large-scale rotation. The only difference between each individual simulation is the detailed structure of the turbulent velocity field. The multiple systems formed in the simulations have properties consistent with observed multiple systems. Dynamical evolution tends preferentially to eject lower mass stars and brown dwarves whilst hardening the remaining binaries so that the median semi-major axis of binaries formed is $\sim 30$ au. Ejected objects are usually single low-mass stars and brown dwarfs, yielding a strong correlation between mass and multiplicity. Our simulations suggest a natural mechanism for forming binary stars that does not require large-scale rotation, capture, or large amounts of turbulence.

We evaluate the intrinsic three dimensional shapes of molecular cores, by analysing their projected shapes. We use the recent catalogue of molecular line observations of Jijina et al. and model the data by the method originally devised for elliptical galaxies. Our analysis broadly supports the conclusion of Jones et al. that molecular cores are better represented by triaxial intrinsic shapes (ellipsoids) than biaxial intrinsic shapes (spheroids). However, we find that the best fit to all of the data is obtained with more extreme axial ratios ($1:0.8:0.4$) than those derived by Jones et al. More surprisingly, we find that starless cores have more extreme axial ratios than protostellar cores -- starless cores appear more `flattened'. This is the opposite of what would be expected from modeling the freefall collapse of triaxial ellipsoids. The collapse of starless cores would be expected to proceed most swiftly along the shortest axis - as has been predicted by every modeller since Zel'dovich - which should produce more flattened cores around protostars, the opposite of what is seen.

We present the results of a detailed non-LTE analysis of the UV and optical spectrum of the O6.5Iaf+ star HD153919 - the mass donor in the high-mass X-ray binary 4U1700-37. Given the eclipsing nature of the system these results allow us to determine the most likely masses of both components of the binary via Monte Carlo simulations. These suggest a mass for HD153919 of 58+/-11M_sun - implying the initial mass of the companion was rather high (>60 M_sun). The most likely mass for the compact companion is found to be 2.44+/-0.27M_sun, with only 3.5 per cent of the trials resulting in a mass less than 2.0M_sun and none less than 1.65M_sun. Our observational data is inconsistent with the canonical neutron star mass and the lowest black hole mass observed (>4.4M_sun; Nova Vel). Significantly changing observational parameters can force the compact object mass into either of these regimes but this results in the O-star mass changing by factors of greater than 2, well beyond the limits determined from its evolutionary state and surface gravity. The low mass of the compact object implies that it is difficult to form high mass black holes through both the Case A & B mass transfer channels and, if the compact object is a neutron star, would significantly constrain the high density nuclear EoS

The Hydra $N$-body hydrodynamics code has been modified to model, from the end of the Sedov phase, the effects of supernovae on the surrounding medium. The motivation is to investigate the feedback of energy into the interstellar/intergalactic medium. We compare our results for supernova remnants (SNRs) in a uniform medium to previous detailed work on the late evolution of SNRs. The code is found to reproduce the bulk characteristics of SNRs well. Results on the effects of a single central SNR on Plummer clouds are presented. The feedback of kinetic energy and the percentage mass loss can be parameterised in terms of the cloud mass and characteristic radius in a simple way. The kinetic energy fraction returned to the ISM from a SNR is $<3$ per cent. The removal of gas from the cold, dense phase and the addition of energy due to the lowering of the potential energy of a cloud is at least as significant, if not much more so, than the kinetic energy leaving a cloud.

The effects of differences between the local and global values of the Hubble parameter on the cosmologies consistent with studies of high-redshift Type Ia supernov\ae are discussed. It is found that with a local Hubble parameter around 10 per cent higher than the global value then open cosmological models (such as $\Omega_{M} = 0.3, \Omega_{\Lambda}=0$) are prefered and if the local value is around 20 per cent higher then standard cosmological models ($\Omega_{M} = 1, \Omega_{\Lambda}=0$) can be recovered. Even in the case where the Hubble parameter ratio is 1, low $\Omega_M$ open cosmologies with $\Omega_{\Lambda} = 0$ are not rejected at the 95 per cent confidence level.

Measurements of the linear diameters of 12 nearby spiral galaxies with distances determined from primary indicators suggest that both the Milky Way Galaxy and M31 are in the middle of the range of sizes for such galaxies. By comparing the measured linear diameters of these nearby systems with the inferred diameters of a sample of more than 3000 spirals with known redshifts, we conclude that the most likely value of the Hubble Parameter lies in the range 50 - 55 km/s Mpc^{-1}.

It has been thought for many years that the Milky Way is an overly large spiral galaxy. Using Cephied distances to 17 spiral galaxies we calculate the true linear diameters of those galaxies. These diameters are then compared to that of the Milky Way which is found to be, at most, an averagely sized spiral galaxy. When compared to galaxies of approximately the same Hubble type ($2<T<6$) the Milky Way is found to be slightly undersized. This suggests that the Hubble parameter is at the lower end of the currently accepted range of possibilities.

The globular clusters of the LMC are found to be significantly more elliptical than Galactic globular clusters, but very similar in virtually all other respects. The ellipticity of the LMC globular clusters is shown not be correlated with the age or mass of those clusters. It is proposed that the ellipticity differences are caused by the different strengths of the tidal fields in the LMC and the Galaxy. The strong Galactic tidal field erases initial velocity anisotropies and removes angular momentum from globular clusters making them more spherical. The tidal field of the LMC is not strong enough to perform these tasks and its globular clusters remain close to their initial states.

N-body simulations are used to model the early evolution of globular clusters. These simulations include residual gas which was not turned into stars which is expelled from the globular cluster by the actions of massive stars. The results of these simulations are compared to observations of 8 LMC globular clusters less than 100 Myr old. These observations are used to constrain the initial conditions that may have produced these clusters. It is found that the entire variety of young LMC globular clusters may be explained in a model where they form from a fairly uniform population of roughly spherical, relaxed proto-cluster clouds very similar to Giant Molecular Clouds in the Galaxy, with star formation efficiencies between 25% and 60%. This paper has been accepted for publication in MNRAS.

The results of N-body simulations of the expulsion of residual gas from young globular clusters are presented. Globular clusters with a variety of initial masses, Galactocentric radii, concentration and initial mass function slope with star formation efficiencies (SFEs) below 50% were investigated. The residual gas in the clusters was simply treated as an external potential acting upon the star particles which was reduced in a number of ways to simulate the loss of the gas. The states of the clusters after ~50 Myr was compared to the results of Chernoff & Shapiro (1987) in order to estimate if the clusters would be able to survive for a Hubble time or would be disrupted. The rapid expulsion of large amounts of a cluster's initial mass is found to considerably effect the structures of clusters, but they may be able to survive with SFEs far lower than 50%. It is suggested that the Galactic globular cluster population could have formed from a proto-cluster populations with SFEs ~40% and central densities approximately the same as those found in giant molecular clouds.