Source author record

Simon P. Goodwin

Simon P. Goodwin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2022arXiv

Few Observation Binary Orbit Solver (FOBOS) from two (or more) astrometric observations

We have developed a new, fast method of estimating the orbital properties of a binary or triple system using as few as two epochs of astrometric data. FOBOS (Few Observation Binary Orbit Solver) uses a flat prior brute force Monte Carlo method to produce probability density functions of the likely orbital parameters. We test the code on fake observations and show that it can (fairly often) constrain the semi-major axis to within a factor of 2-3, and the inclination to within $\sim$20$^{\circ}$ from only two astrometric observations. We also show that the 68 and 95 per cent confidence intervals are statistically reliable. Applying this method to triple systems allows the relative inclination of the secondary and tertiary star orbits to be constrained. FOBOS can usually find a statistically significant number of possible matches in CPU minutes for binary systems, and CPU hours for triple systems.

preprint2020arXiv

Making top-heavy IMFs from canonical IMFs near the Galactic Centre

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.

preprint2020arXiv

The velocity structure of Cygnus OB2

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.

preprint2016arXiv

Discs in misaligned binary systems

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.

preprint2016arXiv

Mass segregation in star clusters is not energy equipartition

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.

preprint2015arXiv

Comparisons between different techniques for measuring mass segregation

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 $Λ_{\rm MSR}$ method; the local surface density $Σ_{\rm LDR}$ method; and the $Ω_{\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 $Ω_{\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 $Λ_{\rm MSR}$ and $Σ_{\rm LDR}$ methods can sometimes produce apparently contradictory results, because they use different definitions of mass segregation. We conclude that only $Λ_{\rm MSR}$ measures mass segregation in the classical sense (without the need for defining the centre of the region), although $Σ_{\rm LDR}$ does place limits on the amount of previous dynamical evolution in a star-forming region.

preprint2014arXiv

Multiplicity in Early Stellar Evolution

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.

preprint2013arXiv

Constraints on Massive Star Formation: Cygnus OB2 was always an Association

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.

preprint2013arXiv

Dynamical evolution of star forming regions

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.

preprint2013arXiv

No compelling evidence of significant early star cluster disruption in the Large Magellanic Cloud

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.

preprint2012arXiv

A natural formation scenario for misaligned and short-period eccentric extrasolar planets

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.

preprint2012arXiv

Binary mass ratios: system mass not primary mass

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.

preprint2012arXiv

On the Origin of Mass Segregation in NGC 3603

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 Γ=-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 $Λ$ minimum spanning tree technique confirms significant mass segregation down to 30 M$_{\odot}$. The dependence of $Λ$ 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.

preprint2012arXiv

Testing the universality of star formation - I. Multiplicity in nearby star-forming regions

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.

preprint2012arXiv

Testing the universality of star formation - II. Comparing separation distributions of nearby star-forming regions and the field

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.

preprint2012arXiv

The same, but different: Stochasticity in binary destruction

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.

preprint2011arXiv

Collisions in young triple systems

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.

preprint2011arXiv

The evolution of binary populations in cool, clumpy star clusters

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.

preprint2010arXiv

Life and Death of Very Massive Stars

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.

preprint2010arXiv

On the mass segregation of stars and brown dwarfs in Taurus

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.

preprint2010arXiv

The binary fraction of the young cluster NGC 1818 in the Large Magellanic Cloud

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

preprint2010arXiv

The dynamical evolution of very-low mass binaries in open clusters

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.

preprint2010arXiv

The early dynamical evolution of cool, clumpy star clusters

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.

preprint2010arXiv

Tidally induced brown dwarf and planet formation in circumstellar discs

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.

preprint2009arXiv

Binaries in star clusters and the origin of the field stellar population

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.

preprint2009arXiv

Do binaries in clusters form in the same way as in the field?

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.

preprint2009arXiv

The role of cluster evolution in disrupting planetary systems and disks: the Kozai mechanism

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.

preprint2005arXiv

Limits on the primordial stellar multiplicity

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.

preprint2004arXiv

Simulating star formation in molecular cores II. The effects of different levels of turbulence

(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, $α_{\rm turb} \equiv U_{\rm turb}/|Ω| = 0, 0.01, 0.025, 0.05, 0.10 {\rm and} 0.25$) and (b), for fixed $α_{\rm turb}$, the details of the initial turbulent velocity field (so as to obtain good statistics). A low level of turbulence ($α_{\rm turb} \sim 0.05$) suffices to produce multiple systems. As $α_{\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.