Source author record

S. P. Goodwin

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

15works
5topics
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

15 published item(s)

preprint2020arXiv

A High Angular Resolution Survey of Massive Stars in Cygnus OB2: $JHK$ Adaptive Optics Results from the Gemini Near-InfraRed Imager

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

preprint2015arXiv

The dynamical evolution of low-mass hydrogen-burning stars, brown dwarfs and planetary-mass objects formed through disc fragmentation

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.

preprint2015arXiv

The M-dwarfs in Multiples (MinMs) survey - I. Stellar multiplicity among low-mass stars within 15 pc

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.

preprint2014arXiv

How does a low-mass cut-off in the stellar IMF affect the evolution of young star clusters?

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.

preprint2013arXiv

A hybrid SPH/N-body method for star cluster simulations

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.

preprint2013arXiv

Convergence of AMR and SPH simulations - I. Hydrodynamical resolution and convergence tests

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.

preprint2013arXiv

Early Disc Accretion as the Origin of Abundance Anomalies in Globular Clusters

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.

preprint2013arXiv

Mapping the core mass function onto the stellar IMF: multiplicity matters

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.

preprint2012arXiv

The VLT-FLAMES Tarantula Survey IV: Candidates for isolated high-mass star formation in 30 Doradus

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

preprint2011arXiv

The origin of very wide binary stars

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.

preprint2010arXiv

The formation of very wide binaries during the star cluster dissolution phase

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.

preprint2009arXiv

The mass ratio and formation mechanisms of Herbig Ae/Be star binary systems

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.

preprint2009arXiv

The origin of very wide binary systems

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.

preprint2009arXiv

What does the IMF really tell us about star formation?

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.

preprint2007arXiv

Adaptive smoothing lengths in SPH

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Δ{\cal N}_{_{\rm NEIB}}$, the tolerance, $Δ{\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 $Δ{\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 $Δ{\cal N}_{_{\rm NEIB}}\geq10$, their results must be viewed with caution.