Search

Circumstellar discs are essential for high mass star formation, while multiplicity, in particular binarity, appears to be an inevitable outcome since the vast majority of massive stars (> 8 Msun) are found in binaries (up to 100%). We constrain the sizes of the dust and ionised gas (Brgamma) emission of the innermost regions towards a sample of six MYSOs, and provide high-mass binary statistics of young stars at 2-300 au scales using VLTI (GRAVITY, AMBER) observations. We determine the inner radius of the dust emission and place MYSOs with K-band measurements in a size-luminosity diagram for the first time, and compare our findings to T Tauris and Herbig AeBes. We also compare the observed K-band sizes to the sublimation radius predicted by three different disc scenarios. Lastly, we apply binary geometries to trace close binarity among MYSOs. The inner sizes of MYSOs, Herbig AeBe and T Tauri stars appear to follow a universal trend at which the sizes scale with the square-root of the stellar luminosity. The Brgamma emission originates from somewhat smaller and co-planar area compared to the 2.2 {\mu}m continuum emission. We discuss this new finding with respect to disc-wind or jet origin. Finally, we report an MYSO binary fraction of 17-25% at milli-arcsecond separations (2-300 au). The size-luminosity diagram indicates that the inner regions of discs around young stars scale with luminosity independently of the stellar mass. At the targeted scales (2-300 au), the MYSO binary fraction is lower than what was previously reported for the more evolved main sequence massive stars, which, if further confirmed, could implicate the predictions from massive binary formation theories. Lastly, we spatially resolve the crucial star/disc interface in a sample of MYSOs, showing that au-scale discs are prominent in high-mass star formation and similar to their low-mass equivalents., Comment: Accepted for publication in A&A. 26 pages, 16 Figures, 7 Tables, full abstract to appear in A&A

We present results of a high angular resolution survey of massive OB stars in the Cygnus OB2 association that we conducted with the NIRI camera and ALTAIR adaptive optics system of the Gemini North telescope. We observed 74 O- and early B-type stars in Cyg OB2 in the $JHK$ infrared bands in order to detect binary and multiple companions. The observations are sensitive to equal-brightness pairs at separations as small as 0.08 \arcsec, and progressively fainter companions are detectable out to $\Delta$ K = 9 mag at a separation of 2 arcsec. This faint contrast limit due to readnoise continues out to 10 arcsec near the edge of the detector. We assigned a simple probability of chance alignment to each companion based upon its separation and magnitude difference from the central target star and upon areal star counts for the general star field of Cyg OB2. Companion stars with a field membership probability of less than 1% are assumed to be physical companions. This assessment indicates that 47% of the targets have at least one resolved companion that is probably gravitationally bound. Including known spectroscopic binaries, our sample includes 27 binary, 12 triple, and 9 systems with four or more components. These results confirm studies of high mass stars in other environments that find that massive stars are born with a high multiplicity fraction. The results are important for the placement of the stars in the H-R diagram, the interpretation of their spectroscopic analyses, and for future mass determinations through measurement of orbital motion., Comment: 53 pages, 5 figures, Accepted to AJ

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., Comment: In press in Astronomy and Astrophysics. arXiv admin note: text overlap with arXiv:1507.05926

We report on the adsorption dynamics of phospholipid membranes on graphene-coated substrates using the quartz crystal microbalance with dissipation monitoring (QCM-D) technique. We compare the lipid vescle interaction and membranne formation on gold and silicon dioxide QCM crystal surfaces with their graphene oxide (GO) and reduced (r)GO coated counterparts, and report on the different lipid structures obtained. We establish graphene derivative coatings as support surfaces with tuneable hydrophobicity for the formation of controllable lipid structures. One structure of interest formed are lipid monolayer membrannes which were formed on rGO, which are otherwise challenging to produce. We also demonstrate and monitor biotin-avidin binding on such a membranne, which will then serve as a platform for a wide range of biosensing applications. The QCM-D technique could be extended to both fundamental studies and applications of other covalent and non-covalent interactions in 2-dimensional materials.

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.

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., Comment: 38 pages, 15 figures

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., Comment: 9 pages, 4 figures. Accepted for publication in MNRAS

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., Comment: Accepted for publication in the Monthly Notices of the Royal Astronomical Society, 2015 February 20. 34 pages and 21 figures, with 7 tables presented here in their entirety. All tables and the full set of figures in the Appendix to appear in the electronic version of the journal

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., Comment: 16 pages, 19 figures, 4 tables, submitted to MNRAS

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., Comment: 13 pages, 7 figures, accepted for publication in MNRAS

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., Comment: 15 pages, 4 figures, MNRAS in press

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., Comment: 10 pages, 4 figures

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., Comment: 19 pages, 9 figures, 1 table, accepted for publication in MNRAS

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., Comment: 20 pages, 8 figures, accepted for publication in MNRAS

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., Comment: 14 pages, 13 figures, 5 tables; Accepted for publication by A&A

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., Comment: Accepted in MNRAS

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., Comment: 5 pages, 3 figures, to be published in the proceedings of the 9th Pacific-Rim Conference on Stellar Astrophysics

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., Comment: 9 pages, 5 figures

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., Comment: 14 pages, 12 figures, accepted by MNRAS

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., Comment: 4 pages, 1 figure, to appear in the proceedings of IAU Symposium 266, eds. R. de Grijs & J.R.D. Lepine

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., Comment: To be published in the proceedings of IAU Symposium 262 "Stellar Populations - Planning for the Next Decade"

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., Comment: Accepted for publication in MNRAS, minor changes made in proof stage

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., Comment: 4 pages, to appear in the proceedings of "BINKEY: Binaries - a Key to Comprehension of the Universe", held 8-12 June 2009 in Brno, Czech Republic

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., Comment: ApJ Letter, accepted

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., Comment: 14 pages, 8 figures. Accepted for publication by A&A

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., Comment: 4 pages, 2 figures, in: Galactic and stellar dynamics in the era of high-resolution surveys, Boily C., Combes F., Hensler G., eds., Strasbourg (France), March 2008, in press (Astron. Nachr.)

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., Comment: 43 pages, 18 figures, accepted for publication in A&A

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., Comment: 9 pages, 5 figures, accepted MNRAS

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)., Comment: 9 pages, 3 figures, accepted for publication in A&A (replacement of v1 in which our results were based on a circular argument)

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., Comment: accepted in ApJ

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

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., Comment: 5 pages, 2 figures, accepted for publication in A&A

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")., Comment: 4 pages, contribution to "Globular Clusters: Guides to Galaxies", March 6th-10th, 2006

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., Comment: 16 pages, 6 figures. To appear in 'Protostars and Planets V'. Includes proto.cls style file. Figures degraded in quality, full version available by anonymous ftp (ftp:astro.shef.ac.uk, pub/spg/GoodwinPPV.pdf.gz)

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., Comment: 7 pages, 3 figures. To appear in A&A

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'')., Comment: 6 pages, 3 figures, MNRAS letters, accepted

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}$., Comment: Accepted for publication in A&A

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., Comment: 7 pages, 1 figure. To appear in Astron.Nachrichten. Includes an.cls file

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., Comment: 6 pages, 2 figures, accepted by A&A

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., Comment: Accepted in the 'Astronomical Notes' for the special issue with the proceedings of the "Ultra low-mass star formation and evolution" workshop

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., Comment: Accepted by A&A (avaliable with high-res images at http://carina.astro.cf.ac.uk/pub/Dimitrios.Stamatellos/publications)

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., Comment: A&A in press, 24 pages with 19 figures. For full resolution version see http://www.edpsciences.org/articles/aa/pdf/press-releases/PRAA200506.pdf (now correct url)

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., Comment: 17 pages, accepted for publication in MNRAS

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, Comment: 13 pages, 7 figures, accepted for publication in A&A

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}., Comment: 9 pages LaTeX file and 4 figures. Also available by e-mailing spg, or by ftp in the directory ftp://star-www.maps.susx.ac.uk/pub/papers/spg/H0/

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