Researcher profile

Jarrod R. Hurley

Jarrod R. Hurley contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

20 published item(s)

preprint2020arXiv

BSE versus StarTrack: implementations of new wind, remnant-formation, and natal-kick schemes in NBODY7 and their astrophysical consequences

The masses of stellar-remnant black holes (BH), as a result of their formation via massive single- and binary-stellar evolution, is of high interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star (BNS) mergers. Here we present new developments in the N-body evolution program NBODY7 in regards to its stellar-remnant formation and related schemes. We demonstrate that the newly-implemented stellar-wind and remnant-formation schemes in the NBODY7 code's BSE sector, such as the 'rapid' and the 'delayed' supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produces neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main sequence (ZAMS) mass and metallicity, with those from the StarTrack population-synthesis program. We also demonstrate the new implementations of various natal-kick mechanisms on NSs and BHs such as the 'convection-asymmetry-driven', 'collapse-asymmetry-driven', and 'neutrino-emission-driven' kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. We find that the SN material fallback causes the convection-asymmetry kick to effectively retain similar number and mass of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry kick would cause nearly all BHs to retain in clusters irrespective of remnant formation model and metallicity, whereas the inference of a large number of BHs in GCs would potentially rule out the neutrino-driven kick mechanism. Pre-SN mergers of massive primordial binaries would cause BH masses to deviate from the single-star ZAMS mass-remnant mass relation. Such mergers, at low metallicities, can produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster.

preprint2020arXiv

Modelling Double Neutron Stars: Radio and Gravitational Waves

We have implemented prescriptions for modelling pulsars in the rapid binary population synthesis code COMPAS. We perform a detailed analysis of the double neutron star (DNS) population, accounting for radio survey selection effects. The surface magnetic field decay timescale (${\sim}1000$\,Myr) and mass scale (${\sim}0.02$\,M$_\odot$) are the dominant uncertainties in our model. Mass accretion during common envelope evolution plays a non-trivial role in recycling pulsars. We find a best-fit model that is in broad agreement with the observed Galactic DNS population. Though the pulsar parameters (period and period derivative) are strongly biased by radio selection effects, the observed orbital parameters (orbital period and eccentricity) closely represent the intrinsic distributions. The number of radio observable DNSs in the Milky Way at present is about 2500 in our model, corresponding to approximately 10\% of the predicted total number of DNSs in the galaxy. Using our model calibrated to the Galactic DNS population, we make predictions for DNS mergers observed in gravitational waves. The DNS chirp mass distribution varies from 1.1M$_\odot$ to 2.1M$_\odot$ and the median is found to be 1.14\,M$_\mathrm{\odot}$. The expected effective spin $χ_\mathrm{eff}$ for isolated DNSs is $\lesssim$0.03 from our model. We predict that 34\% of the current Galactic isolated DNSs will merge within a Hubble time, and have a median total mass of 2.7\,M$_\mathrm{\odot}$. Finally, we discuss implications for fast radio bursts and post-merger remnant gravitational-waves.

preprint2016arXiv

A dynamical gravitational wave source in a dense cluster

Making use of a new N-body model to describe the evolution of a moderate-size globular cluster we investigate the characteristics of the population of black holes within such a cluster. This model reaches core-collapse and achieves a peak central density typical of the dense globular clusters of the Milky Way. Within this high-density environment we see direct confirmation of the merging of two stellar remnant black-holes in a dynamically-formed binary, a gravitational wave source. We describe how the formation, evolution and ultimate ejection/destruction of binary systems containing black holes impacts the evolution of the cluster core. Also, through comparison with previous models of lower density, we show that the period distribution of black hole binaries formed through dynamical interactions in this high-density model favours the production of gravitational wave sources. We confirm that the number of black holes remaining in a star cluster at late times and the characteristics of the binary black hole population depend on the nature of the star cluster, critically on the number density of stars and by extension the relaxation timescale.

preprint2015arXiv

Different Dynamical Ages for the Two Young and Coeval LMC Star Clusters NGC 1805 and NGC 1818 Imprinted on Their Binary Populations

The two Large Magellanic Cloud star clusters NGC 1805 and NGC 1818 are approximately the same chronological age (~30 Myr), but show different radial trends in binary frequency. The F-type stars (1.3 - 2.2 MSun) in NGC 1818 have a binary frequency that decreases towards the core, while the binary frequency for stars of similar mass in NGC 1805 is flat with radius, or perhaps bimodal (with a peak in the core). We show here, through detailed N-body modeling, that both clusters could have formed with the same primordial binary frequency and with binary orbital elements and masses drawn from the same distributions (defined from observations of open clusters and the field of our Galaxy). The observed radial trends in binary frequency for both clusters are best matched with models that have initial substructure. Furthermore, both clusters may be evolving along a very similar dynamical sequence, with the key difference that NGC 1805 is dynamically older than NGC 1818. The F-type binaries in NGC 1818 still show evidence of an initial period of rapid dynamical disruptions (which occur preferentially in the core), while NGC 1805 has already begun to recover a higher core binary frequency, owing to mass segregation (which will eventually produce a distribution in binary frequency that rises only towards the core, as is observed in old Milky Way star clusters). This recovery rate increases for higher-mass binaries, and therefore even at one age in one cluster, we predict a similar dynamical sequence in the radial distribution of the binary frequency as a function of binary primary mass.

preprint2015arXiv

Tracking Cluster Debris (TraCD)- I. Dissolution of clusters and searching the solar cradle

The capability to reconstruct dissolved stellar systems in dynamical and chemical space is a key factor in improving our understanding of the evolution of the Milky Way. Here we concentrate on the dynamical aspect and given that a significant portion of the stars in the Milky Way have been born in stellar associations or clusters that have lived a few Myr up to several Gyr, we further restrict our attention to the evolution of star clusters. We have carried out our simulations in two steps: (1) we create a simulation of dissolution and mixing processes which yields a close fit to the present-day Milky Way dynamics and (2) we have evolved three sets of stellar clusters with masses of 400, 1000 and 15000 M$_{\odot}$ to dissolution. The birth location of these sets was 4, 6, 8 and 10 kpc for the 400 and 1000 M$_{\odot}$ clusters and 4, 6, 8, 10 and 12 kpc for the 15000 M$_{\odot}$. We have focused our efforts on studying the state of the escapers from these clusters after 4.5 Gyr of evolution with particular attention to stars that reach the Solar annulus, i.e. 7.5 $\le$ R$_{\rm{gc}}$ $\le$ 8.5 kpc. We give results for Solar twins and siblings over a wide range of radii and cluster masses for two dissolution mechanisms. From kinematics alone, we conclude that the Sun was $\sim$50 per cent more likely to have been born near its current Galactocentric radius, rather than have migrated (radially) $\sim$2 kpc since birth. We conclude our analysis by calculating magnitudes and colours of our single stars for comparison with the samples that the $Gaia$, $Gaia$-ESO and GALAH-AAO surveys will obtain. In terms of reconstructing dissolved star clusters we find that on short time-scales we cannot rely on kinematic evolution alone and thus it will be necessary to extend our study to include information on chemical space.

preprint2014arXiv

Merging Binary Stars and the magnetic white dwarfs

A magnetic dynamo driven by differential rotation generated when stars merge can explain strong fields in certain classes of magnetic stars, including the high field magnetic white dwarfs (HFMWDs). In their case the site of the differential rotation has been variously proposed to be within a common envelope, the massive hot outer regions of a merged degenerate core or an accretion disc formed by a tidally disrupted companion that is subsequently incorporated into a degenerate core. We synthesize a population of binary systems to investigate the stellar merging hypothesis for observed single HFMWDs. Our calculations provide mass distribution and the fractions of white dwarfs that merge during a common envelope phase or as double degenerate systems in a post common envelope phase. We vary the common envelope efficiency parameter alpha and compare with observations. We find that this hypothesis can explain both the observed incidence of magnetism and the mass distribution of HFMWDs for a wide range of alpha. In this model, the majority of the HFMWDs are of the Carbon Oxygen type and merge within a common envelope. Less than about a quarter of a per cent of HFMWDs originate from double degenerate stars that merge after common envelope evolution and these populate the high-mass tail of the HFMWD mass distribution.

preprint2014arXiv

The Effect of Orbital Eccentricity on the Dynamical Evolution of Star Clusters

We use N-body simulations to explore the influence of orbital eccentricity on the dynamical evolution of star clusters. Specifically we compare the mass loss rate, velocity dispersion, relaxation time, and the mass function of star clusters on circular and eccentric orbits. For a given perigalactic distance, increasing orbital eccentricity slows the dynamical evolution of a cluster due to a weaker mean tidal field. However, we find that perigalactic passes and tidal heating due to an eccentric orbit can partially compensate for the decreased mean tidal field by energizing stars to higher velocities and stripping additional stars from the cluster, accelerating the relaxation process. We find that the corresponding circular orbit which best describes the evolution of a cluster on an eccentric orbit is much less than its semi-major axis or time averaged galactocentric distance. Since clusters spend the majority of their lifetimes near apogalacticon, the properties of clusters which appear very dynamically evolved for a given galactocentric distance can be explained by an eccentric orbit. Additionally we find that the evolution of the slope of the mass function within the core radius is roughly orbit-independent, so it could place additional constraints on the initial mass and initial size of globular clusters with solved orbits. We use our results to demonstrate how the orbit of Milky Way globular clusters can be constrained given standard observable parameters like galactocentric distance and the slope of the mass function. We then place constraints on the unsolved orbits of NGC 1261,NGC 6352, NGC 6496, and NGC 6304 based on their positions and mass functions.

preprint2014arXiv

The Effects of Orbital Inclination on the Scale Size and Evolution of Tidally Filling Star Clusters

We have performed N-body simulations of tidally filling star clusters with a range of orbits in a Milky Way-like potential to study the effects of orbital inclination and eccentricity on their structure and evolution. At small galactocentric distances Rgc, a non-zero inclination results in increased mass loss rates. Tidal heating and disk shocking, the latter sometimes consisting of two shocking events as the cluster moves towards and away from the disk, help remove stars from the cluster. Clusters with inclined orbits at large Rgc have decreased mass loss rates than the non-inclined case, since the strength the disk potential decreases with Rgc. Clusters with inclined and eccentric orbits experience increased tidal heating due to a constantly changing potential, weaker disk shocks since passages occur at higher Rgc, and an additional tidal shock at perigalacticon. The effects of orbital inclination decrease with orbital eccentricity, as a highly eccentric cluster spends the majority of its lifetime at a large Rgc. The limiting radii of clusters with inclined orbits are best represented by the rt of the cluster when at its maximum height above the disk, where the cluster spends the majority of its lifetime and the rate of change in rt is a minimum. Conversely, the effective radius is independent of inclination in all cases.

preprint2014arXiv

The impact of galaxy geometry and mass evolution on the survival of star clusters

Direct N-body simulations of globular clusters in a realistic Milky Way-like potential are carried out using the code NBODY6 to determine the impact of the host galaxy disk mass and geometry on the survival of star clusters. A relationship between disk mass and star cluster dissolution timescale is derived. These N-body models show that doubling the mass of the disk from 5x10^10 solar masses to 10x10^10 solar masses halves the dissolution time of a satellite star cluster orbiting the host galaxy at 6 kpc from the galactic center. Different geometries in a disk of identical mass can determine either the survival or dissolution of a star cluster orbiting within the inner 6 kpc of the galactic center. Furthermore, disk geometry has measurable effects on the mass loss of star clusters up to 15 kpc from the galactic center. N-body simulations performed with a fine output time step show that at each disk crossing the outer layers of star clusters experience an increase in velocity dispersion of ~5% of the average velocity dispersion in the outer section of star clusters. This leads to an enhancement of mass-loss -- a clearly discernable effect of disk shocking. By running models with different inclinations we determine that star clusters with an orbit perpendicular to the Galactic plane have larger mass loss rates than both clusters evolving in the Galactic plane or in an inclined orbit.

preprint2013arXiv

Consequences of Dynamical Disruption and Mass Segregation for the Binary Frequencies of Star Clusters

The massive (13,000-26,000 Msolar), young (15-30 Myr) Large Magellanic Cloud star cluster NGC 1818 reveals an unexpected increasing binary frequency with radius for F-type stars (1.3-2.2 Msolar). This is in contrast to many older star clusters that show a decreasing binary frequency with radius. We study this phenomenon with sophisticated N-body modeling, exploring a range of initial conditions, from smooth virialized density distributions to highly substructured and collapsing configurations. We find that many of these models can reproduce the cluster's observed properties, although with a modest preference for substructured initial conditions. Our models produce the observed radial trend in binary frequency through disruption of soft binaries (with semi-major axes, a > 3000 AU), on approximately a crossing time (~ 5.4 Myr), preferentially in the cluster core. Mass segregation subsequently causes the binaries to sink towards the core. After roughly one initial half-mass relaxation time (t_rh(0) ~ 340 Myr) the radial binary frequency distribution becomes bimodal, the innermost binaries having already segregated towards the core, leaving a minimum in the radial binary frequency distribution that marches outwards with time. After 4-6 t_rh(0), the rising distribution in the halo disappears, leaving a radial distribution that rises only towards the core. Thus, both a radial binary frequency distribution that falls towards the core (as observed for NGC 1818) and one that rises towards the core (as for older star clusters) can arise naturally from the same evolutionary sequence owing to binary disruption and mass segregation in rich star clusters.

preprint2013arXiv

The Influence of Orbital Eccentricity on Tidal Radii of Star Clusters

We have performed N-body simulations of star clusters orbiting in a spherically symmetric smooth galactic potential. The model clusters cover a range of initial half-mass radii and orbital eccentricities in order to test the historical assumption that the tidal radius of a cluster is imposed at perigalacticon. The traditional assumption for globular clusters is that since the internal relaxation time is larger than its orbital period, the cluster is tidally stripped at perigalacticon. Instead, our simulations show that a cluster with an eccentric orbit does not need to fully relax in order to expand. After a perigalactic pass, a cluster re-captures previously unbound stars, and the tidal shock at perigalacticon has the effect of energizing inner region stars to larger orbits. Therefore, instead of the limiting radius being imposed at perigalacticon, it more nearly traces the instantaneous tidal radius of the cluster at any point in the orbit. We present a numerical correction factor to theoretical tidal radii calculated at perigalacticon which takes into consideration both the orbital eccentricity and current orbital phase of the cluster.

preprint2012arXiv

A direct N-body model of core-collapse and core oscillations

We report on the results of a direct N-body simulation of a star cluster that started with N = 200 000, comprising 195 000 single stars and 5 000 primordial binaries. The code used for the simulation includes stellar evolution, binary evolution, an external tidal field and the effects of two-body relaxation. The model cluster is evolved to 12 Gyr, losing more than 80% of its stars in the process. It reaches the end of the main core-collapse phase at 10.5 Gyr and experiences core oscillations from that point onwards -- direct numerical confirmation of this phenomenon. However, we find that after a further 1 Gyr the core oscillations are halted by the ejection of a massive binary comprised of two black holes from the core, producing a core that shows no signature of the prior core-collapse. We also show that the results of previous studies with N ranging from 500 to 100 000 scale well to this new model with larger N. In particular, the timescale to core-collapse (in units of the relaxation timescale), mass segregation, velocity dispersion, and the energies of the binary population all show similar behaviour at different N.

preprint2012arXiv

Direct N-Body Modeling of the Old Open Cluster NGC 188: A Detailed Comparison of Theoretical and Observed Binary Star and Blue Straggler Populations

(Abridged) Following on from a recently completed radial-velocity survey of the old (7 Gyr) open cluster NGC 188 in which we study in detail the solar-type hard binaries and blue stragglers of the cluster, here we investigate the dynamical evolution of NGC 188 through a sophisticated N-body model. We employ the observed binary properties of the young (150 Myr) open cluster M35, where possible, to guide our choices for parameters of the initial binary population. At 7 Gyr the main-sequence solar-type hard-binary population in the model matches that of NGC 188 in both binary frequency and distributions of orbital parameters. This agreement between the model and observations is in a large part due to the similarities between the NGC 188 and M35 solar-type binaries. Indeed, among the 7 Gyr main-sequence binaries in the model, only those with P>1000 days show potentially observable evidence for modifications by dynamical encounters. This emphasizes the importance of defining accurate initial conditions for star cluster models, which we propose is best accomplished through comparisons with observations of young open clusters like M35. Furthermore, this suggests that observations of the present-day binaries in even old open clusters can provide valuable information on their primordial binary populations. However, despite the model's successes at matching the true cluster, the model underproduces blue stragglers and produces an overabundance of long-period circular main-sequence--white-dwarf binaries as compared to NGC 188. We conclude that improvements in the physics of mass transfer and common envelope may in fact solve both discrepancies with the observations. This project highlights the unique accessibility of open clusters to both comprehensive observational surveys and full-scale N-body simulations, and underscores the importance of open clusters to the study of star cluster dynamics.

preprint2012arXiv

Multiple stellar-mass black holes in globular clusters: theoretical confirmation

While tens or hundreds of stellar-remnant black holes are expected to form in globular star clusters, it is still unclear how many of those will be retained upon formation, and how many will be ejected through subsequent dynamical interactions. No such black holes have been found in any Milky Way globular cluster until the recent discovery of stellar-mass black holes in the globular cluster M22 (NGC 6656) with now an estimated population of 5-100 black holes. We present a direct N-body model of a star cluster of the same absolute and dynamical age as M22. Imposing an initial retention fraction of approx. 10% for black holes, 16 stellar-remnant black holes are retained at a cluster age of 12 Gyr, in agreement with the estimate for M22. Of those 16 BHs, two are in a binary system with a main sequence star each while also one pure black hole binary is present. We argue that multiple black holes can be present in any Milky Way cluster with an extended core radius, such as M22 or the model presented here.

preprint2012arXiv

N-body models of globular clusters: metallicity, half-light radii and mass-to-light ratios

Size differences of approx. 20% between red (metal-rich) and blue (metal-poor) sub-populations of globular clusters have been observed, generating an ongoing debate as to weather these originate from projection effects or the difference in metallicity. We present direct N-body simulations of metal-rich and metal-poor stellar populations evolved to study the effects of metallicity on cluster evolution. The models start with N = 100000 stars and include primordial binaries. We also take metallicity dependent stellar evolution and an external tidal field into account. We find no significant difference for the half-mass radii of those models, indicating that the clusters are structurally similar. However, utilizing observational tools to fit half-light (or effective) radii confirms that metallicity effects related to stellar evolution combined with dynamical effects such as mass segregation produce an apparent size difference of 17% on average. The metallicity effect on the overall cluster luminosity also leads to higher mass-to-light ratios for metal-rich clusters.

preprint2012arXiv

The Size Scale of Star Clusters

Direct N-body simulations of star clusters in a realistic Milky Way-like potential are carried out using the code NBODY6. Based on these simulations a new relationship between scale size and galactocentric distance is derived: the scale size of star clusters is proportional to the hyperbolic tangent of the galactocentric distance. The half-mass radius of star clusters increases systematically with galactocentric distance but levels off when star clusters orbit the galaxy beyond ~40 kpc. These simulations show that the half-mass radius of individual star clusters varies significantly as they evolve over a Hubble time, more so for clusters with shorter relaxation times, and remains constant through several relaxation times only in certain situations when expansion driven by the internal dynamics of the star cluster and the influence of the host galaxy tidal field balance each other. Indeed, the radius of a star cluster evolving within the inner 20 kpc of a realistic galactic gravitational potential is severely truncated by tidal interactions and does not remain constant over a Hubble time. Furthermore, the half-mass radius of star clusters measured with present day observations bears no memory of the original cluster size. Stellar evolution and tidal stripping are the two competing physical mechanisms that determine the present day size of globular clusters. These simulations also show that extended star clusters can form at large galactocentric distances while remaining fully bound to the host galaxy. There is thus no need to invoke accretion from an external galaxy to explain the presence of extended clusters at large galactocentric distances in a Milky Way-type galaxy.

preprint2010arXiv

N-body Models of Extended Clusters

We use direct N-body simulations to investigate the evolution of star clusters with large size-scales with the particular goal of understanding the so-called extended clusters observed in various Local Group galaxies, including M31 and NGC6822. The N-body models incorporate a stellar mass function, stellar evolution and the tidal field of a host galaxy. We find that extended clusters can arise naturally within a weak tidal field provided that the tidal radius is filled at the start of the evolution. Differences in the initial tidal filling-factor can produce marked differences in the subsequent evolution of clusters and the size-scales that would be observed. These differences are more marked than any produced by internal evolution processes linked to the properties of cluster binary stars or the action of an intermediate-mass black hole, based on models performed in this work and previous work to date. Models evolved in a stronger tidal field show that extended clusters cannot form and evolve within the inner regions of a galaxy such as M31. Instead our results support the suggestion many extended clusters found in large galaxies were accreted as members of dwarf galaxies that were subsequently disrupted. Our results also enhance the recent suggestion that star clusters evolve to a common sequence in terms of their size and mass.

preprint2010arXiv

On The Maximum Mass of Stellar Black Holes

We present the spectrum of compact object masses: neutron stars and black holes that originate from single stars in different environments. In particular, we calculate the dependence of maximum black hole mass on metallicity and on some specific wind mass loss rates (e.g., Hurley et al. and Vink et al.). Our calculations show that the highest mass black holes observed in the Galaxy M_bh = 15 Msun in the high metallicity environment (Z=Zsun=0.02) can be explained with stellar models and the wind mass loss rates adopted here. To reach this result we had to set Luminous Blue Variable mass loss rates at the level of about 0.0001 Msun/yr and to employ metallicity dependent Wolf-Rayet winds. With such winds, calibrated on Galactic black hole mass measurements, the maximum black hole mass obtained for moderate metallicity (Z=0.3 Zsun=0.006) is M_bh,max = 30 Msun. This is a rather striking finding as the mass of the most massive known stellar black hole is M_bh = 23-34 Msun and, in fact, it is located in a small star forming galaxy with moderate metallicity. We find that in the very low (globular cluster-like) metallicity environment the maximum black hole mass can be as high as M_bh,max = 80 Msun (Z=0.01 Zsun=0.0002). It is interesting to note that X-ray luminosity from Eddington limited accretion onto an 80 Msun black hole is of the order of about 10^40 erg/s and is comparable to luminosities of some known ULXs. We emphasize that our results were obtained for single stars only and that binary interactions may alter these maximum black hole masses (e.g., accretion from a close companion). This is strictly a proof-of-principle study which demonstrates that stellar models can naturally explain even the most massive known stellar black holes.

preprint2001arXiv

Free-Floating Planets: Not So Surprising

We find that free-floating planets can remain bound to a star cluster for much longer than was previously assumed: of the order of the cluster half-mass relaxation timescale as opposed to the crossing-time. This result is based on N-body simulations performed with the new GRAPE-6 special purpose hardware and is important in the context of the preliminary detection of a population of free-floating sub-stellar objects in the globular cluster M22. The planets in our N-body study are of Jupiter mass and are initially placed in circular orbits of between 0.05 and 50 AU about a parent star whose mass is chosen from a realistic initial mass function. The presence of the free-floating planets is the result of dynamical encounters between planetary systems and the cluster stars. Most planets are liberated from their parent star in, or near, the cluster core, and then drift outwards on a timescale of ~ 10^{8-9} year. This still implies the existence of many (~100) planets per star if the M22 result is confirmed.

preprint2000arXiv

Direct N-body Modelling of Stellar Populations: Blue Stragglers in M67

We present a state-of-the-art N-body code which includes a detailed treatment of stellar and binary evolution as well as the cluster dynamics. This code is ideal for investigating all aspects relating to the evolution of star clusters and their stellar populations. It is applicable to open and globular clusters of any age. We use the N-body code to model the blue straggler population of the old open cluster M67. Preliminary calculations with our binary population synthesis code show that binary evolution alone cannot explain the observed numbers or properties of the blue stragglers. On the other hand, our N-body model of M67 generates the required number of blue stragglers and provides formation paths for all the various types found in M67. This demonstrates the effectiveness of the cluster environment in modifying the nature of the stars it contains and highlights the importance of combining dynamics with stellar evolution. We also perform a series of N = 10000 simulations in order to quantify the rate of escape of stars from a cluster subject to the Galactic tidal field.