BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Alison Sills

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

ResearcherUnclaimed source record
astro-ph.GAastro-ph.SRastro-phastro-ph.CO

Catalog footprint

What is connected

49works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

49 published item(s)

preprint2022arXiv

Constraining the shape of dark matter haloes with globular clusters

We explore how diffuse stellar light and globular clusters (GCs) can be used to trace the matter distribution of their host halo using an observational methodology. For this, we use 117 simulated dark matter (DM) haloes from the $(34.4~\rm cMpc)^3$ periodic volume of the E-MOSAICS project. For each halo, we compare the stellar surface brightness and GC projected number density maps to the surface densities of DM and total mass. We find that the dominant structures identified in the stellar light and in the GCs correspond closely with those from the DM and total mass. Our method is unaffected by the presence of satellites and its precision improves with fainter GC samples. We recover tight relations between the profiles of stellar surface brightness and GC number density to those of the DM, suggesting that the profile of DM can be accurately recovered from the stars and GCs ($σ\leq0.5~$dex). We quantify the projected morphology of DM, stars and GCs, and find that the stars and GCs are more flattened than the DM. Additionally, the semi-major axes of the distribution of stars and GCs are typically misaligned by $\sim 10~$degrees from that of DM. We demonstrate that deep imaging of diffuse stellar light and GCs can place constraints on the shape, profile and orientation of their host halo. These results extend down to haloes with central galaxies $M_{\star}\geq10^{10}~M_{\odot}$, and the analysis will be applicable to future data from the Euclid, Roman and the Rubin observatories.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Introducing EMP-Pathfinder: modelling the simultaneous formation and evolution of stellar clusters in their host galaxies

The formation and evolution of stellar clusters is intimately linked to that of their host galaxies. To study this connection, we present the EMP-Pathfinder suite of cosmological zoom-in Milky Way-mass simulations. These simulations contain a sub-grid description for stellar cluster formation and evolution, allowing us to study the simultaneous formation and evolution of stellar clusters alongside their host galaxies across cosmic time. As a key ingredient in these simulations, we include the physics of the multi-phase nature of the interstellar medium (ISM), which enables studies of how the presence of a cold, dense ISM affects cluster formation and evolution. We consider two different star formation prescriptions: a constant star formation efficiency per free-fall time, as well as an environmentally-dependent, turbulence-based prescription. We identify two key results drawn from these simulations. Firstly, we find that tidal shock-driven disruption caused by the graininess of the cold ISM produces old ($τ>10~$Gyr) stellar cluster populations with properties that are in excellent agreement with the observed populations in the Milky Way and M31. Importantly, the addition of the cold ISM addresses the areas of disagreement found in previous simulations that lacked the cold gas phase. Secondly, the formation of stellar clusters is extremely sensitive to the baryonic physics that govern the properties of the cold, dense gas reservoir in the galaxy. This implies that the demographics of stellar cluster populations represent an important diagnostic tool for constraining baryonic physics models in upcoming galaxy formation simulations that also include a description of the cold ISM.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Modelling Star Cluster Formation: Mergers

Star cluster formation in giant molecular clouds involves the local collapse of the cloud into small gas-rich subclusters, which can then subsequently collide and merge to build up the final star cluster(s). In this paper, we simulate collisions between these subclusters, using coupled smooth particle hydrodynamics for the gas and N-body dynamics for the stars. We are guided by previous radiation hydrodynamics simulations of molecular cloud collapse which provide the global properties of the colliding clusters, such as their stellar and gas masses, and their initial positions and velocities. The subclusters in the original simulation were treated as sink particles which immediately merged into a single entity after the collision. We show that the more detailed treatment provides a more complex picture. At collisional velocities above ~ 10 km/s, the stellar components of the cluster do not form a monolithic cluster within 3 Myr, although the gas may do so. At lower velocities, the clusters do eventually merge but over timescales that may be longer than the time for a subsequent collision. The structure of the resultant cluster is not well-fit by any standard density distribution, and the clusters are not in equilibrium but continue to expand over our simulation time. We conclude that the simple sink particle treatment of subcluster mergers in large-scale giant molecular cloud simulations provides an upper limit on the final cluster properties.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Implementing Primordial Binaries in Simulations of Star Cluster Formation with a Hybrid MHD and Direct N-Body Method

The fraction of stars in binary systems within star clusters is important for their evolution, but what proportion of binaries form by dynamical processes after initial stellar accretion remains unknown. In previous work, we showed that dynamical interactions alone produced too few low-mass binaries compared to observations. We therefore implement an initial population of binaries in the coupled MHD and direct N-body star cluster formation code Torch. We compare simulations with, and without, initial binary populations and follow the dynamical evolution of the binary population in both sets of simulations, finding that both dynamical formation and destruction of binaries take place. Even in the first few million years of star formation, we find that an initial population of binaries is needed at all masses to reproduce observed binary fractions for binaries with mass ratios above the $q \geq 0.1$ detection limit. Our simulations also indicate that dynamical interactions in the presence of gas during cluster formation modify the initial distributions towards binaries with smaller primary masses, larger mass ratios, smaller semi-major axes and larger eccentricities. Systems formed dynamically do not have the same properties as the initial systems, and systems formed dynamically in the presence of an initial population of binaries differ from those formed in simulations with single stars only. Dynamical interactions during the earliest stages of star cluster formation are important for determining the properties of binary star systems.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Constraining Mass Transfer Histories of Blue Straggler Stars with COS Spectroscopy of White Dwarf Companions

Recent studies show that the majority of blue straggler stars (BSSs) in old open clusters are formed through mass transfer from an evolved star onto a main-sequence companion, resulting in a BSS and white dwarf (WD) in a binary system. We present constraints on the mass transfer histories of two BSS-WD binaries in the open cluster NGC 188, using WD temperatures and surface gravities measured with HST COS far-ultraviolet spectroscopy. Adopting a Gaia-based cluster distance of $1847\pm107$ pc, we determine that one system, WOCS 4540, formed through Case C mass transfer resulting in a CO-core white dwarf with $T_{\text{eff}}=17000^{+140}_{-200}$ K and a $\log g=7.80^{+0.06}_{-0.06}$, corresponding to a mass of $0.53^{+0.03}_{-0.03}$ $M_{\odot}$ and a cooling age of $105^{+6}_{-5}$ Myr. The other system, WOCS 5379, formed through Case B mass transfer resulting in a He-core white dwarf with $T_{\text{eff}}=15500^{+170}_{-150}$ K and a $\log g=7.50^{+0.06}_{-0.05}$, corresponding to a mass of $0.42^{+0.02}_{-0.02}$ $M_{\odot}$ and an age of $250^{+20}_{-20}$ Myr. The WD parameters are consistent across four different cluster distance assumptions. We determine possible progenitor binary systems with a grid of accretion models using MESA, and investigate whether these systems would lead to stable or unstable mass transfer. WOCS 4540 likely resulted from stable mass transfer during periastron passage in an eccentric binary system, while WOCS 5379 challenges our current understanding of the expected outcomes for mass transfer from red giant branch stars. Both systems are examples of the value in using detailed analyses to fine-tune our physical understanding of binary evolutionary processes.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Emergence of Multiple Populations in Radiation Hydrodynamics Simulations of Cluster Formation

We present a new approach to understanding star-to-star helium abundance variations within globular clusters. We begin with detailed radiation hydrodynamics simulations of cluster formation within giant molecular clouds, and investigate the conditions under which multiple populations could be created. Chemical enrichment occurs dynamically as the cluster is assembled. We test two extreme mechanisms for injection of enriched gas within the clusters, and find that realistic multiple populations can be formed in both mechanisms. The stochastic cluster formation histories are dictated by the inherent randomness of the timing and location of the formation of small clusters, which rapidly merge to build up the larger cluster, in combination with continual accretion of gas from the cloud. These cluster formation histories naturally produce a diversity of abundance patterns across the massive cluster population. We conclude that multiple populations are a natural outcome of the typical mode of star cluster formation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Globular Cluster Scale Sizes in Giant Galaxies: Orbital Anisotropy and Tidally Under-filling Clusters in M87, NGC 1399, and NGC 5128

We investigate the shallow increase in globular cluster half-light radii with projected galactocentric distance $R_{gc}$ observed in the giant galaxies M87, NGC 1399, and NGC 5128. To model the trend in each galaxy, we explore the effects of orbital anisotropy and tidally under-filling clusters. While a strong degeneracy exists between the two parameters, we use kinematic studies to help constrain the distance $R_β$ beyond which cluster orbits become anisotropic, as well as the distance $R_{fα}$ beyond which clusters are tidally under-filling. For M87 we find $R_β> 27$ kpc and $20 < R_{fα} < 40$ kpc and for NGC 1399 $R_β> 13$ kpc and $10 < R_{fα} < 30$ kpc. The connection of $R_{fα}$ with each galaxy's mass profile indicates the relationship between size and $R_{gc}$ may be imposed at formation, with only inner clusters being tidally affected. The best fitted models suggest the dynamical histories of brightest cluster galaxies yield similar present-day distributions of cluster properties. For NGC 5128, the central giant in a small galaxy group, we find $R_β> 5$ kpc and $R_{fα} > 30$ kpc. While we cannot rule out a dependence on $R_{gc}$, NGC 5128 is well fitted by a tidally filling cluster population with an isotropic distribution of orbits, suggesting it may have formed via an initial fast accretion phase. Perturbations from the surrounding environment may also affect a galaxy's orbital anisotropy profile, as outer clusters in M87 and NGC 1399 have primarily radial orbits while outer NGC 5128 clusters remain isotropic.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Hubble Space Telescope Study of the Enigmatic Milky Way Halo Globular Cluster Crater

We analyze the resolved stellar populations of the faint stellar system, Crater, based on deep optical imaging taken with the Hubble Space Telescope. The HST/ACS-based color-magnitude diagram (CMD) of Crater extends $\sim$4 magnitudes below the oldest main sequence turnoff, providing excellent leverage on Crater's physical properties. Structurally, Crater has a half-light radius of $\sim$20 pc and shows no evidence for tidal distortions. Crater is well-described by a simple stellar population with an age of $\sim$7.5 Gyr, [M/H]$\sim-1.65$, a M$_{\star}\sim10^4$ M$_{\odot}$, M$_{\rm V}\sim -5.3$, located at a distance of d$\sim$ 145 kpc, with modest uncertainties in these properties due to differences in the underlying stellar evolution models. The sparse sampling of stars above the turnoff and sub-giant branch are likely to be 1.0-1.4 M$_{\odot}$ binary star systems (blue stragglers) and their evolved descendants, as opposed to intermediate age main sequence stars. Confusion of these populations highlights a substantial challenge in accurately characterizing sparsely populated stellar systems. Our analysis shows that Crater is not a dwarf galaxy, but instead is an unusually young cluster given its location in the Milky Way's very outer stellar halo. Crater is similar to SMC cluster Lindsay 38, and its position and velocity are in good agreement with observations and models of the Magellanic stream debris, suggesting it may have accreted from the Magellanic Clouds. However, its age and metallicity are also in agreement with the age-metallicity relationships of lower mass dwarf galaxies such as Leo I or Carina. Despite uncertainty over its progenitor system, Crater appears to have been incorporated into the Galaxy more recently than $z\sim1$ (8 Gyr ago), providing an important new constraint on the accretion history of the Milky Way. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Implications for the Formation of Blue Straggler Stars from HST Ultraviolet Observations of NGC 188

We present results of a Hubble Space Telescope far-ultraviolet (FUV) survey searching for white dwarf (WD) companions to blue straggler stars (BSSs) in open cluster NGC 188. The majority of NGC 188 BSSs (15 of 21) are single-lined binaries with properties suggestive of mass-transfer formation via Roche lobe overflow, specifically through an asymptotic giant branch star transferring mass to a main sequence secondary, yielding a BSS binary with a WD companion. In NGC 188, a BSS formed by this mechanism within the past 400 Myr will have a WD companion hot and luminous enough to be directly detected as a FUV photometric excess with HST. Comparing expected BSS FUV emission to observed photometry reveals four BSSs with WD companions above 12,000 K (younger than 250 Myr) and three WD companions with temperatures between 11,000-12,000 K. These BSS+WD binaries all formed through recent mass transfer. The location of the young BSSs in an optical color-magnitude diagram (CMD) indicates that distance from the zero-age main sequence does not necessarily correlate with BSS age. There is no clear CMD separation between mass transfer-formed BSSs and those likely formed through other mechanisms, such as collisions. The seven detected WD companions place a lower limit on the mass-transfer formation frequency of 33%. We consider other possible formation mechanisms by comparing properties of the BSS population to theoretical predictions. We conclude that 14 BSS binaries likely formed from mass transfer, resulting in an inferred mass-transfer formation frequency of approximately 67%.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Dynamics of Multiple Populations in the Globular Cluster NGC 6362

We investigate how the Milky Way tidal field can affect the spatial mixing of multiple stellar populations in the globular cluster NGC 6362. We use $N$-body simulations of multiple population clusters on the orbit of this cluster around the Milky Way. Models of the formation of multiple populations in globular clusters predict that the second population should initially be more centrally concentrated than the first. However, NGC 6362 is comprised of two chemically distinct stellar populations having the same radial distribution. We show that the high mass loss rate experienced on this cluster's orbit significantly accelerates the spatial mixing of the two populations expected from two body relaxation. We also find that for a range of initial second population concentrations, cluster masses, tidal filling factors and fraction of first population stars, a cluster with two populations should be mixed when it has lost 70-80 per cent of its initial mass. These results fully account for the complete spatial mixing of NGC 6362, since, based on its shallow present day mass function, independent studies estimate that the cluster has lost 85 per cent of its initial mass.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Properties of Bound and Unbound Molecular Cloud Populations Formed in Galactic Disc Simulations

We explore the effect of galactic environment on properties of molecular clouds. Using clouds formed in a large-scale galactic disc simulation, we measure the observable properties from synthetic column density maps. We confirm that a significant fraction of unbound clouds forms naturally in a galactic disc environment and that a mixed population of bound and unbound clouds can match observed scaling relations and distributions for extragalactic molecular clouds. By dividing the clouds into inner and outer disc populations, we compare their distributions of properties and test whether there are statistically significant differences between them. We find that clouds in the outer disc have lower masses, sizes, and velocity dispersions as compared to those in the inner disc for reasonable choices of the inner/outer boundary. We attribute the differences to the strong impact of galactic shear on the disc stability at large galactocentric radii. In particular, our Toomre analysis of the disc shows a narrowing envelope of unstable masses as a function of radius, resulting in the formation of smaller, lower mass fragments in the outer disc. We also show that the star formation rate is affected by the environment of the parent cloud, and is particularly influenced by the underlying surface density profile of the gas throughout the disc. Our work highlights the strengths of using galaxy-scale simulations to understand the formation and evolution of cloud properties - and the star formation within them - in the context of their environment.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Spatial Structure of Young Stellar Clusters. III. Physical Properties and Evolutionary States

We analyze the physical properties of stellar clusters that are detected in massive star-forming regions in the MYStIX project--a comparative, multiwavelength study of young stellar clusters within 3.6 kpc that contain at least one O-type star. Tabulated properties of subclusters in these regions include physical sizes and shapes, intrinsic numbers of stars, absorptions by the molecular clouds, and median subcluster ages. Physical signs of dynamical evolution are present in the relations of these properties, including statistically significant correlations between subcluster size, central density, and age, which are likely the result of cluster expansion after gas removal. We argue that many of the subclusters identified in Paper I are gravitationally bound because their radii are significantly less than what would be expected from freely expanding clumps of stars with a typical initial stellar velocity dispersion of ~3 km/s for star-forming regions. We explore a model for cluster formation in which structurally simpler clusters are built up hierarchically through the mergers of subclusters--subcluster mergers are indicated by an inverse relation between the numbers of stars in a subcluster and their central densities (also seen as a density vs. radius relation that is less steep than would be expected from pure expansion). We discuss implications of these effects for the dynamical relaxation of young stellar clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Detection of white dwarf companions to blue stragglers in the open cluster NGC 188: direct evidence for recent mass transfer

Several possible formation pathways for blue straggler stars have been developed recently, but no one pathway has yet been observationally confirmed for a specific blue straggler. Here we report the first findings from a Hubble Space Telescope ACS/SBC far-UV photometric program to search for white dwarf companions to blue straggler stars. We find three hot and young white dwarf companions to blue straggler stars in the 7-Gyr open cluster NGC 188, indicating that mass transfer in these systems ended less than 300 Myr ago. These companions are direct and secure observational evidence that these blue straggler stars were formed through mass transfer in binary stars. Their existence in a well-studied cluster environment allows for observational constraints of both the current binary system and the progenitor binary system, mapping the entire mass transfer history.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Evolving Molecular Cloud Structure and the Column Density Probability Distribution Function

The structure of molecular clouds can be characterized with the probability distribution function (PDF) of the mass surface density. In particular, the properties of the distribution can reveal the nature of the turbulence and star formation present inside the molecular cloud. In this paper, we explore how these structural characteristics evolve with time and also how they relate to various cloud properties as measured from a sample of synthetic column density maps of molecular clouds. We find that, as a cloud evolves, the peak of its column density PDF will shift to surface densities below the observational threshold for detection, resulting in an underlying lognormal distribution which has been effectively lost at late times. Our results explain why certain observations of actively star-forming, dynamically older clouds, such as the Orion molecular cloud, do not appear to have any evidence of a lognormal distribution in their column density PDFs. We also study the evolution of the slope and deviation point of the power-law tails for our sample of simulated clouds and show that both properties trend towards constant values, thus linking the column density structure of the molecular cloud to the surface density threshold for star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Models of Individual Blue Stragglers

This chapter describes the current state of models of individual blue stragglers. Stellar collisions, binary mergers (or coalescence), and partial or ongoing mass transfer have all been studied in some detail. The products of stellar collisions retain memory of their parent stars and are not fully mixed. Very high initial rotation rates must be reduced by an unknown process to allow the stars to collapse to the main sequence. The more massive collision products have shorter lifetimes than normal stars of the same mass, while products between low mass stars are long-lived and look very much like normal stars of their mass. Mass transfer can result in a merger, or can produce another binary system with a blue straggler and the remnant of the original primary. The products of binary mass transfer cover a larger portion of the colour-magnitude diagram than collision products for two reasons: there are more possible configurations which produce blue stragglers, and there are differing contributions to the blended light of the system. The effects of rotation may be substantial in both collision and merger products, and could result in significant mixing unless angular momentum is lost shortly after the formation event. Surface abundances may provide ways to distinguish between the formation mechanisms, but care must be taking to model the various mixing mechanisms properly before drawing strong conclusions. Avenues for future work are outlined.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modifying two-body relaxation in N-body systems by gas accretion

We consider the effects that accretion from the interstellar medium onto the particles of an N-body system has on the rate of two-body relaxation. To this end, we derive an accretion-modified relaxation time by adapting Spitzer's two-component model to include the damping effects of accretion. We consider several different mass-dependencies and efficiency factors for the accretion rate, as well as different mass ratios for the two components of the model. The net effect of accretion is to accelerate mass segregation by increasing the average mass $\bar{m}$, since the relaxation time is inversely proportional to $\bar{m}$. Under the assumption that the accretion rate increases with the accretor mass, there are two additional effects that accelerate mass segregation. First, accretion acts to increase the range of any initial mass spectrum, quickly driving the heaviest members to even higher masses. Second, accretion acts to reduce the velocities of the accretors due to conservation of momentum, and it is the heaviest members that are affected the most. Using our two-component model, we quantify these effects as a function of the accretion rate, the total cluster mass, and the component masses. We conclude by discussing the implications of our results for the dynamical evolution of primordial globular clusters, primarily in the context of black holes formed from the most massive stellar progenitors.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Blue Straggler Star Population in NGC 1261: Evidence for a Post-Core-Collapse Bounce State

We present a multi-passband photometric study of the Blue Straggler Star (BSS) population in the Galactic globular cluster (GC) NGC\,1261, using available space- and ground-based survey data.~The inner BSS population is found to have two distinct sequences in the color-magnitude diagram, similar to double BSS sequences detected in other GCs. These well defined sequences are presumably linked to single short-lived events such as core collapse, which are expected to boost the formation of BSSs.~In agreement with this, we find a BSS sequence in NGC\,1261 which can be well reproduced individually by a theoretical model prediction of a 2 Gyr old population of stellar collision products, which are expected to form in the denser inner regions during short-lived core contraction phases.~Additionally, we report the occurrence of a group of BSSs with unusually blue colours in the CMD, which are consistent with a corresponding model of a 200 Myr old population of stellar collision products.~The properties of the NGC\,1261 BSS populations, including their spatial distributions, suggest an advanced dynamical evolutionary state of the cluster, but the core of this GC does not show the classical signatures of core-collapse.~We argue these apparent contradictions provide evidence for a post-core-collapse bounce state seen in dynamical simulations of old GCs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Celestial Buffet: multiple populations and globular cluster formation in dwarf galaxies

We present a framework that explains the commonly observed variation in light element abundances in globular clusters. If globular clusters form in the centres of dwarf galaxies, they will be pumped onto larger orbits as star formation progresses. The potential well will only retain the moderate velocity AGB ejecta, the expected source of enrichment, but not supernova ejecta. There is no need to increase the initial cluster mass, a requirement of self-enrichment scenarios, as all the stars within the dwarf can contribute. As the clusters move through the dwarf centre they sweep up a mix of AGB ejecta and in-falling pristine gas to form a second generation of stars. The specific mix will vary in time and is thus able to explain the spread in second generation abundances observed in different clusters. The globular clusters will survive to the present day or be stripped as part of the hierarchical merging process of larger galaxies. We illustrate how this process may operate using a high-resolution simulation of a dwarf galaxy at high redshift.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Size of Star Clusters Accreted by the Milky Way

We perform N-body simulations of a cluster that forms in a dwarf galaxy and is then accreted by the Milky Way to investigate how a cluster's structure is affected by a galaxy merger. We find that the cluster's half mass radius will respond quickly to this change in potential. When the cluster is placed on an orbit in the Milky Way with a stronger tidal field the cluster experiences a sharp decrease in size in response to increased tidal forces. Conversely, when placed on an orbit with a weaker tidal field the cluster expands since tidal forces decrease and stars moving outwards due to internal effects remain bound at further distances than before. In all cases, we find that the cluster's half mass radius will eventually be indistinguishable from a cluster that has always lived in the Milky Way on that orbit. These adjustments occur within 1-2 half mass relaxation times of the cluster in the dwarf galaxy. We also find this effect to be qualitatively independent of the time that the cluster is taken from the dwarf galaxy. In contrast to the half mass radius, we show the core radius of the cluster is not affected by the potential the cluster lives in. Our work suggests that structural properties of accreted clusters are not distinct from clusters born in the Milky Way. Other cluster properties, such as metallicity and horizontal branch morphology, may be the only way to identify accreted star clusters in the Milky Way.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Spatial Structure of Young Stellar Clusters. I. Subclusters

The clusters of young stars in massive star-forming regions show a wide range of sizes, morphologies, and numbers of stars. Their highly subclustered structures are revealed by the MYStIX project's sample of 31,754 young stars in nearby sites of star formation (regions at distances <3.6 kpc that contain at least one O-type star.) In 17 of the regions surveyed by MYStIX, we identify subclusters of young stars using finite mixture models -- collections of isothermal ellipsoids that model individual subclusters. Maximum likelihood estimation is used to estimate the model parameters, and the Akaike Information Criterion is used to determine the number of subclusters. This procedure often successfully finds famous subclusters, such as the BN/KL complex behind the Orion Nebula Cluster and the KW-object complex in M17. A catalog of 142 subclusters is presented, with 1 to 20 subclusters per region. The subcluster core radius distribution for this sample is peaked at 0.17 pc with a standard deviation of 0.43 dex, and subcluster core radius is negatively correlated with gas/dust absorption of the stars -- a possible age effect. Based on the morphological arrangements of subclusters, we identify four classes of spatial structure: long chains of subclusters, clumpy structures, isolated clusters with a core-halo structure, and isolated clusters well fit by a single isothermal ellipsoid.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The state of globular clusters at birth II: primordial binaries

(abridged) In this paper, we constrain the properties of primordial binary populations in Galactic globular clusters using the MOCCA Monte Carlo code for cluster evolution. Our results are compared to the observations of Milone et al. (2012) using the photometric binary populations as proxies for the true underlying distributions, in order to test the hypothesis that the data are consistent with an universal initial binary fraction near unity and the binary orbital parameter distributions of Kroupa (1995). With the exception of a few possible outliers, we find that the data are to first-order consistent with the universality hypothesis. Specifically, the present-day binary fractions inside the half-mass radius r$_{\rm h}$ can be reproduced assuming either high initial binary fractions near unity with a dominant soft binary component as in the Kroupa distribution combined with high initial densities (10$^4$-10$^6$ M$_{\odot}$ pc$^{-3}$), or low initial binary fractions ($\sim$ 5-10%) with a dominant hard binary component combined with moderate initial densities near their present-day values (10$^2$-10$^3$ M$_{\odot}$ pc$^{-3}$). This apparent degeneracy can be broken using the binary fractions outside r$_{\rm h}$- only high initial binary fractions with a significant soft component combined with high initial densities can contribute to reproducing the observed anti-correlation between the binary fractions outside r$_{\rm h}$ and the total cluster mass. We further illustrate using the simulated present-day binary orbital parameter distributions and the technique introduced in Leigh et al. (2012) that the relative fractions of hard and soft binaries can be used to further constrain the initial cluster density and mass-density relation. Our results favour an initial mass-density relation of the form r$_{\rm h} \propto$ M$_{\rm clus}^α$ with $α<$ 1/3.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The state of globular clusters at birth: emergence from the gas-embedded phase

In this paper, we discuss the origin of the observed correlation between cluster concentration c and present-day mass function (PDMF) slope α reported by De Marchi, Paresce & Pulone. This relation can either be reproduced from universal initial conditions combined with some dynamical mechanism(s) that alter(s) the cluster structure and mass function over time, or it must arise early on in the cluster lifetime, such as during the gas-embedded phase of cluster formation. Using a combination of Monte Carlo and N-body models for globular cluster evolution performed with the MOCCA and NBODY6 codes, respectively, we explore a number of dynamical mechanisms that could affect the observed relation. For the range of initial conditions considered here, our results are consistent with an universal initial binary fraction \approx 10% (which does not, however, preclude 100%) and an universal initial stellar mass function resembling the standard Kroupa distribution. Most of the dispersion observed in the c-α relation can be attributed to two-body relaxation and Galactic tides. However, dynamical processes alone could not have reproduced the dispersion in concentration, and we require at least some correlation between the initial concentration and the total cluster mass. We argue that the origin of this trend could be connected to the gas-embedded phase of cluster evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Unbound Star-forming Molecular Clouds

We explore whether observed molecular clouds could include a substantial population of unbound clouds. Using simulations which include only turbulence and gravity, we are able to match observed relations and naturally reproduce the observed scatter in the cloud size-linewidth coefficient, at fixed surface density. We identify the source of this scatter as a spread in the intrinsic virial parameter. Thus these observational trends do not require that clouds exist in a state of dynamical equilibrium. We demonstrate that cloud virial parameters can be accurately determined observationally with an appropriate size estimator. All our simulated clouds eventually form collapsing cores, regardless of whether the cloud is bound overall. This supports the idea that molecular clouds do not have to be bound to form stars or to have observed properties like those of nearby low-mass clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Globular Cluster Scale Sizes in Giant Galaxies: The Case of M87 and the Role of Orbital Anisotropy and Tidal Filling

We present new Hubble Space Telescope imaging of the outer regions of M87 in order to study its globular cluster (GC) population out to large galactocentric distances. We discuss particularly the relationship between GC effective radii $r_h$ and projected galactocentric distance $R_{gc}$. The observations suggest a shallow trend $r_h \propto R_{gc}^{0.14}$ out to $R_{gc} \sim 100$ kpc, in agreement with studies of other giant elliptical galaxies. To theoretically reproduce this relationship we simulate GC populations with various distributions of orbits. For an isotropic distribution of cluster orbits we find a steeper trend of $r_h \propto R_{gc}^{0.4}$. Instead we suggest that (a) if the cluster system has an orbital anisotropy profile, where orbits become preferentially radial with increasing galactocentric distance, AND (b) if clusters become more tidally under-filling with galactocentric distance, the observed relationship can be recovered. We also apply this approach to the red and blue GC populations separately and predict that red clusters are preferentially under-filling at large $R_{gc}$ and have a more isotropic distribution of orbits than blue clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Mock Observations of Blue Stragglers in Globular Cluster Models

We created artificial color-magnitude diagrams of Monte Carlo dynamical models of globular clusters, and then used observational methods to determine the number of blue stragglers in those clusters. We compared these blue stragglers to various cluster properties, mimicking work that has been done for blue stragglers in Milky Way globular clusters to determine the dominant formation mechanism(s) of this unusual stellar population. We find that a mass-based prescription for selecting blue stragglers will choose approximately twice as many blue stragglers than a selection criterion that was developed for observations of real clusters. However, the two numbers of blue stragglers are well-correlated, so either selection criterion can be used to characterize the blue straggler population of a cluster. We confirm previous results that the simplified prescription for the evolution of a collision or merger product in the BSE code overestimates their lifetimes. We show that our model blue stragglers follow similar trends with cluster properties (core mass, binary fraction, total mass, collision rate) as the true Milky Way blue stragglers, as long as we restrict ourselves to model clusters with an initial binary fraction higher than 5%. We also show that, in contrast to earlier work, the number of blue stragglers in the cluster core does have a weak dependence on the collisional parameter Gamma in both our models and in Milky Way globular clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Overview of the Massive Young Star-Forming Complex Study in Infrared and X-ray (MYStIX) Project

MYStIX (Massive Young Star-Forming Complex Study in Infrared and Xray) seeks to characterize 20 OB-dominated young clusters and their environs at distances d < 4 kpc using imaging detectors on the Chandra X-ray Observatory, Spitzer Space Telescope, and the United Kingdom InfraRed Telescope. The observational goals are to construct catalogs of star-forming complex stellar members with well-defined criteria, and maps of nebular gas (particularly of hot X-ray emitting plasma) and dust. A catalog of MYStIX Probable Complex Members (MPCMs) with several hundred OB stars and > 30, 000 low mass premain sequence is assembled. This sample and related data products will be used to seek new empirical constraints on theoretical models of cluster formation and dynamics, mass segregation, OB star formation, star formation triggering on the periphery of HII regions, the survivability of protoplanetary disks in HII regions. This paper give an introduction and overview of the project, covering the data analysis methodology and application to two star forming regions, NGC 2264 and the Trifid Nebula.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stellar Collisions and Blue Straggler Stars in Dense Globular Clusters

Blue straggler stars (BSS) are abundantly observed in all Galactic globular clusters (GGC) where data exist. However, observations alone cannot reveal the relative importance of various formation channels or the typical formation times for this well studied population of anomalous stars. Using a state-of-the-art Hénon-type Monte Carlo code that includes all relevant physical processes, we create 128 models with properties typical of the observed GGCs. These models include realistic numbers of single and binary stars, use observationally motivated initial conditions, and span large ranges in central density, concentration, binary fraction, and mass. Their properties can be directly compared with those of observed GGCs. We can easily identify the BSSs in our models and determine their formation channels and birth times. We find that for central densities above ~10^3 Msun/pc^3 the dominant formation channel is stellar collisions while for lower density clusters, mass transfer in binaries provides a significant contribution (up to ~ 60% in our models). The majority of these collisions are binary-mediated, occurring during 3-body and 4-body interactions. As a result a strong correlation between the specific frequency of BSSs and the binary fraction in a cluster can be seen in our models. We find that the number of BSSs in the core shows only a weak correlation with the collision rate estimator Γtraditionally used by observers, in agreement with the latest Hubble Space Telescope (ACS) data. Using an idealized "full mixing" prescription for collision products, our models indicate that the BSSs observed today may have formed several Gyrs ago. However, denser clusters tend to have younger (~1 Gyr) BSSs.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Connecting the Dots: Analyzing Synthetic Observations of Star-Forming Clumps in Molecular Clouds

In this paper, we investigate the extent to which observations of molecular clouds can correctly identify and measure star-forming clumps. We produced a synthetic column density map and a synthetic spectral-line data cube from the simulated collapse of a 5000 M$_{\odot}$ molecular cloud. By correlating the clumps found in the simulation to those found in the synthetic observations, clump masses derived from spectral-line data cubes were found to be quite close to the true physical properties of the clumps. We also find that the `observed' clump mass function derived from the column density map is shifted by a factor of ~ 3 higher than the true clump mass function, due to projection of low-density material along the line of sight. Alves et al. (2007) first proposed that a shift of a clump mass function to higher masses by a factor of 3 can be attributed to a star formation efficiency of 30 %. Our results indicate that this finding may instead be due to an overestimate of clump masses determined from column density observations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quantifying the Universality of the Stellar Initial Mass Function in Old Star Clusters

We present a new technique to quantify cluster-to-cluster variations in the observed present-day stellar mass functions of a large sample of star clusters. Our method quantifies these differences as a function of both the stellar mass and the total cluster mass, and offers the advantage that it is insensitive to the precise functional form of the mass function. We applied our technique to data taken from the ACS Survey for Globular Clusters, from which we obtained completeness-corrected stellar mass functions in the mass range 0.25-0.75 M$_{\odot}$ for a sample of 27 clusters. The results of our observational analysis were then compared to Monte Carlo simulations for globular cluster evolution spanning a range of initial mass functions, total numbers of stars, concentrations, and virial radii. We show that the present-day mass functions of the clusters in our sample can be reproduced by assuming an universal initial mass function for all clusters, and that the cluster-to-cluster differences are consistent with what is expected from two-body relaxation. A more complete exploration of the initial cluster conditions will be needed in future studies to better constrain the precise functional form of the initial mass function. This study is a first step toward using our technique to constrain the dynamical histories of a large sample of old Galactic star clusters and, by extension, star formation in the early Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Observational and Theoretical Tidal Radii of Globular Clusters in M87

Globular clusters have linear sizes (tidal radii) which theory tells us are determined by their masses and by the gravitational potential of their host galaxy. To explore the relationship between observed and expected radii, we utilize the globular cluster population of the Virgo giant M87. Unusually deep, high signal-to-noise images of M87 are used to measure the effective and limiting radii of approximately 2000 globular clusters. To compare with these observations, we simulate a globular cluster population that has the same characteristics as the observed M87 cluster population. Placing these simulated clusters in the well-studied tidal field of M87, the orbit of each cluster is solved and the theoretical tidal radius of each cluster is determined. We compare the predicted relationship between cluster size and projected galactocentric distance to observations. We find that for an isotropic distribution of cluster velocities, theoretical tidal radii are approximately equal to observed limiting radii for Rgc < 10 kpc. However, the isotropic simulation predicts a steep increase in cluster size at larger radii, which is not observed in large galaxies beyond the Milky Way. To minimize the discrepancy between theory and observations, we explore the effects of orbital anisotropy on cluster sizes, and suggest a possible orbital anisotropy profile for M87 which yields a better match between theory and observations. Finally, we suggest future studies which will establish a stronger link between theoretical tidal radii and observed radii.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Origins of Blue Stragglers and Binarity in Globular Clusters

(abridged) We use newly available empirical binary fractions for globular clusters to carry out a direct test of the binary evolution hypothesis, and of collisional channels that involve binary stars. More specifically, using the previously reported correlation between blue straggler numbers and core mass as a benchmark, we test for correlations with the number of binary stars, as well as with the rates of single-single, single-binary, and binary-binary encounters. Surprisingly, we find that the simple correlation with core mass remains by far the strongest predictor of blue straggler population size, even in our joint models. This is despite the fact that the binary fractions themselves strongly anti-correlate with core mass, just as expected in the binary evolution model. At first sight, these results do not fit neatly with either binary evolution or collisional models in their simplest forms. Arguably the simplest and most intriguing possibility to explain this unexpected result is that observational errors on the core binary fractions are larger than the true intrinsic dispersion associated with their dependence on core mass. In the context of the binary evolution model, this would explain why the combination of binary fraction and core mass is a poorer predictor of blue straggler numbers than core mass alone. It would also imply that core mass is a remarkably clean predictor of core binary fractions. This would be of considerable importance for the dynamical evolution of globular clusters, and provides an important benchmark for models attempting to understand their present-day properties.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Size Difference Between Red And Blue Globular Clusters Is NOT Due To Projection Effects

Metal-rich (red) globular clusters in massive galaxies are, on average, smaller than metal-poor (blue) globular clusters. One of the possible explanations for this phenomenon is that the two populations of clusters have different spatial distributions. We test this idea by comparing clusters observed in unusually deep, high signal-to-noise images of M87 with a simulated globular cluster population in which the red and blue clusters have different spatial distributions, matching the observations. We compare the overall distribution of cluster effective radii as well as the relationship between effective radius and galactocentric distance for both the observed and simulated red and blue subpopulations. We find that the different spatial distributions does not produce a significant size difference between the red and blue subpopulations as a whole, or at a given galactocentric distance. These results suggest that the size difference between red and blue globular clusters is likely due to differences during formation or later evolution

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

An Analytic Model for Blue Straggler Formation in Globular Clusters

We present an analytic model for blue straggler formation in globular clusters. We assume that blue stragglers are formed only through stellar collisions and binary star evolution, and compare our predictions to observed blue straggler numbers taken from the catalogue of Leigh, Sills & Knigge (2011). We can summarize our key results as follows: (1) Binary star evolution consistently dominates blue straggler production in all our best-fitting models. (2) In order to account for the observed sub-linear dependence of blue straggler numbers on the core masses (Knigge, Leigh & Sills 2009), the core binary fraction must be inversely proportional to the total cluster luminosity and should always exceed at least a few percent. (3) In at least some clusters, blue straggler formation must be enhanced by dynamical encounters (either via direct collisions or by stimulating mass-transfer to occur by altering the distribution of binary orbital parameters) relative to what is expected by assuming a simple population of binaries evolving in isolation. (4) The agreement between the predictions of our model and the observations can be improved by including blue stragglers that form outside the core but later migrate in due to dynamical friction. (5) Longer blue straggler lifetimes are preferred in models that include blue stragglers formed outside the core since this increases the fraction that will have sufficient time to migrate in via dynamical friction.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dissecting the Colour-Magnitude Diagram: A Homogeneous Catalogue of Stellar Populations in Globular Clusters

We present a homogeneous catalogue for blue straggler, red giant branch, horizontal branch and main-sequence turn-off stars in a sample of 35 clusters taken from the ACS Survey for Globular Clusters. As a result of the superior photometry and relatively large field of view offered by the ACS data, this new catalogue is a significant improvement upon the one presented in Leigh, Sills & knigge (2007). Using our catalogue, we study and compare the radial distributions of the different stellar populations. We have confirmed our previous result (Knigge, Leigh & Sills 2009) that there is a clear, but sub-linear, correlation between the number of blue stragglers found in the cluster core and the total stellar mass contained within it. By considering a larger spatial extent than just the core, our results suggest that mass segregation is not the dominant effect contributing to the observed sub-linearity. We also investigate the radial distributions of the different stellar populations in our sample of clusters. Our results are consistent with a linear relationship between the number of stars in these populations and the total mass enclosed within the same radius. Therefore, we conclude that the cluster dynamics does not significantly affect the relative distributions of these populations in our sample.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Effect of Stellar Metallicity on the Sizes of Star Clusters

Observations indicate blue globular clusters have half-light radii systematically larger than those of red globular clusters. In this paper, we test whether the different metallicity-dependent stellar evolution timescales and mass-loss rates within the clusters can impact their early dynamical evolution. By means of N-body simulations including stellar evolution recipes we simulate the early evolution of small centrally concentrated clusters with and without primordial mass segregation. Our simulations include accurate metallicity-dependent mass loss from massive stars. We find blue clusters to be larger than red clusters regardless of whether the clusters have been primordially mass segregated. In addition, the size difference is found to be larger and consistent with observations for initial models with a low central concentration. These results indicate that the systematic size difference found between red and blue clusters can, at least in part, be attributed to the dynamical effects of differing stellar evolution histories, driven by metallicity.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

An Analytic Technique for Constraining the Dynamical Origins of Multiple Star Systems Containing Merger Products

We present a technique to identify the most probable dynamical formation scenario for an observed binary or triple system containing one or more merger products or, alternatively, to rule out the possibility of a dynamical origin. Our method relies on an analytic prescription for energy conservation during stellar encounters. With this, observations of the multiple star system containing the merger product(s) can be used to work backwards in order to constrain the initial orbital energies of any single, binary or triple systems that went into the encounter. The initial semi-major axes of the orbits provide an estimate for the collisional cross section and therefore the time-scale for the encounter to occur in its host cluster. We have applied our analytic prescription to observed binary and triple systems containing blue stragglers, in particular the triple system S1082 in M67 and the period distribution of the blue straggler binaries in NGC 188. We have shown that both S1082 and most of the blue straggler binaries in NGC 188 could have a dynamical origin, and that encounters involving triples could be a significant contributor to BS populations in old open clusters. In general, our results suggest that encounters involving triples could make up a significant fraction of those dynamical interactions that result in stellar mergers, in particular encounters that produce multiple star systems containing one or more blue stragglers.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Blue Straggler Formation in Clusters

Blue stragglers are thought to be formed from the merger or coalescence of two stars, but the details of their formation in clusters has been difficult to disentangle. We discuss the two main formation mechanisms for blue stragglers (stellar collisions or mass transfer in a binary system). We then look at the additional complications caused by the stars living in the dynamically active environment of a star cluster. We review the recent observational and theoretical work which addresses the question "which mechanism dominates?" and conclude that the most likely answer is that both mechanisms are at work, although with different importance in different environments and at different times in the cluster lifetime. We finish with a short discussion of some avenues for future work

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Mass Transfer in Binary Stars using SPH. I. Numerical Method

Close interactions and mass transfer in binary stars can lead to the formation of many different exotic stellar populations, but detailed modeling of mass transfer is a computationally challenging problem. Here, we present an alternate Smoothed Particle Hydrodynamics approach to the modeling of mass transfer in binary systems that allows a better resolution of the flow of matter between main-sequence stars. Our approach consists of modeling only the outermost layers of the stars using appropriate boundary conditions and ghost particles. We arbitrarily set the radius of the boundary and find that our boundary treatment behaves physically and conserves energy well. In particular, when used with our binary relaxation procedure, our treatment of boundary conditions is also shown to evolve circular binaries properly for many orbits. The results of our first simulation of mass transfer are also discussed and used to assess the strengths and limitations of our method. We conclude that it is well suited for the modeling of interacting binary stars. The method presented here represents a convenient alternative to previous hydrodynamical techniques aimed at modeling mass transfer in binary systems since it can be used to model both the donor and the accretor while maintaining the density profiles taken from realistic stellar models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Mass Transfer in Binary Stars using SPH. II. Eccentric Binaries

Despite numerous efforts to better understand binary star evolution, some aspects of it remain poorly constrained. In particular, the evolution of eccentric binaries has remained elusive mainly because the Roche lobe formalism derived for circular binaries does not apply. Here, we report the results of our Smoothed Particle Hydrodynamics simulations of mass transfer in eccentric binaries using an alternate method in which we model only the outermost layers of the stars with appropriate boundary conditions. Using this technique, along with properly relaxed model stars, we characterize the mass transfer episodes of binaries with various orbital parameters. In particular, we show that these episodes can be described by Gaussians with a FWMH of ~0.12 P_orb and that the peak rates occur after periastron, at an orbital phase of ~0.58, independently of the eccentricity and mass of the stars. The accreted material is observed to form a rather sparse envelope around either or both stars. Although the fate of this envelope is not modeled in our simulations, we show that a constant fraction (~5%) of the material transferred is ejected from the systems. We discuss this result in terms of the non-conservative mass transfer scenario. We suggest our results could be incorporated in analytical and binary population synthesis studies to help better understand the evolution of eccentric binaries and the formation of exotic stellar populations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multiple Populations in Globular Clusters: The Possible Contributions of Stellar Collisions

Globular clusters were thought to be simple stellar populations, but recent photometric and spectroscopic evidence suggests that the clusters' early formation history was more complicated. In particular, clusters show star-to-star abundance variations, and multiple sequences in their colour-magnitude diagrams. These effects seem to be restricted to globular clusters, and are not found in open clusters or the field. In this paper, we combine the two competing models for these multiple populations and include a consideration of the effects of stellar collisions. Collisions are one of the few phenomena which occur solely in dense stellar environments like (proto-)globular clusters. We find that runaway collisions between massive stars can produce material which has abundances comparable to the observed second generations, but that very little total mass is produced by this channel. We then add the contributions of rapidly-rotating massive stars (under the assumption that massive stars are spun up by collisions and interactions), and the contribution of asymptotic giant branch stars. We find that collisions can help produce the extreme abundances which are seen in some clusters. However, the total amount of material produced in these generations is still too small (by at least a factor of 10) to match the observations. We conclude with a discussion of the additional effects which probably need to be considered to solve this particular problem.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multiple stellar populations and their influence on blue stragglers

It has become clear in recent years that globular clusters are not simple stellar populations, but may host chemically distinct sub-populations, typically with an enhanced helium abundance. These helium-rich populations can make up a substantial fraction of all cluster stars. One of the proposed formation channels for blue straggler stars is the physical collision and merger of two stars. In the context of multiple populations, collisions between stars with different helium abundances should occur and contribute to the observed blue straggler population. This will affect the predicted blue straggler colour and luminosity function. We quantify this effect by calculating models of mergers resulting from collisions between stars with different helium abundances and using these models to model a merger population. We then compare these results to four observed clusters, NGC 1851, NGC 2808, NGC 5634 and NGC 6093. As in previous studies our models deviate from the observations, particularly in the colour distributions. However, our results are consistent with observations of multiple populations in these clusters. In NGC 2808, our best fitting models include normal and helium enhanced populations, in agreement with helium enhancement inferred in this cluster. The other three clusters show better agreement with models that do not include helium enhancement. We discuss future prospects to improve the modelling of blue straggler populations and the role that the models we present here can play in such a study.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stellar Encounters with Multiple Star Systems and the Blue Straggler Problem

We present a technique to identify the most probable dynamical formation scenario for an observed binary or triple system containing one or more merger products or, alternatively, to rule out the possibility of a dynamical origin. Our method relies on an analytic prescription for energy conservation during stellar encounters. With this, observations of the multiple star system containing the merger product(s) can be used to work backwards in order to constrain the initial orbital energies of any single, binary or triple systems that went into the encounter. The initial semi-major axes of the orbits provide an estimate for the collisional cross section and therefore the time-scale for the encounter to occur in its host cluster. We have applied our analytic prescription to observed binary and triple systems containing blue stragglers, in particular the triple system S1082 in M67 and the period distribution of the blue straggler binaries in NGC 188. We have shown that both S1082 and most of the blue straggler binaries in NGC 188 could have a dynamical origin, and that encounters involving triples are a significant contributor to blue straggler populations in old open clusters. In general, our results suggest that encounters involving triples could make up a significant fraction of those dynamical interactions that result in stellar mergers, in particular encounters that produce multiple star systems containing one or more blue stragglers.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Blue Stragglers After the Main Sequence

We study the post-main sequence evolution of products of collisions between main sequence stars (blue stragglers), with particular interest paid to the horizontal branch and asymptotic giant branch phases. We found that the blue straggler progeny populate the colour-magnitude diagram slightly blueward of the red giant branch and between 0.2 and 1 magnitudes brighter than the horizontal branch. We also found that the lifetimes of collision products on the horizontal branch is consistent with the numbers of so-called "evolved blue straggler stars" (E-BSS) identified by various authors in a number of globular clusters, and is almost independent of mass or initial composition profile. The observed ratio of the number of E-BSS to blue stragglers points to a main sequence lifetime for blue stragglers of approximately 1-2 Gyr on average.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Where the Blue Stragglers Roam: Searching for a Link Between Formation and Environment

The formation of blue stragglers is still not completely understood, particularly the relationship between formation environment and mechanism. We use a large, homogeneous sample of blue stragglers in the cores of 57 globular clusters to investigate the relationships between blue straggler populations and their environments. We use a consistent definition of "blue straggler" based on position in the color-magnitude diagram, and normalize the population relative to the number of red giant branch stars in the core. We find that the previously determined anti-correlation between blue straggler frequency and total cluster mass is present in the purely core population. We find some weak correlations with central velocity dispersion and with half-mass relaxation time. The blue straggler frequency does not show any trend with any other cluster parameter. Even though collisions may be expected to be a dominant blue straggler formation process in globular cluster cores, we find no correlation between the frequency of blue stragglers and the collision rate in the core. We also investigated the blue straggler luminosity function shape, and found no relationship between any cluster parameter and the distribution of blue stragglers in the color-magnitude diagram. Our results are inconsistent with some recent models of blue straggler formation that include collisional formation mechanisms, and may suggest that almost all observed blue stragglers are formed in binary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

The Dynamical Implications of Multiple Stellar Formation Events in Galactic Globular Clusters

Various galactic globular clusters display abundance anomalies that affect the morphology of their colour-magnitude diagrams. In this paper we consider the possibility of helium enhancement in the anomalous horizontal branch of NGC 2808. We examine the dynamics of a self-enrichment scenario in which an initial generation of stars with a top-heavy initial mass function enriches the interstellar medium with helium via the low-velocity ejecta of its asymptotic giant branch stars. This enriched medium then produces a second generation of stars which are themselves helium-enriched. We use a direct N-body approach to perform five simulations and conclude that such two-generation clusters are both possible and would not differ significantly from their single-generation counterparts on the basis of dynamics. We find, however, that the stellar populations of such clusters would differ from single-generation clusters with a standard initial mass function and in particular would be enhanced in white dwarf stars. We conclude, at least from the standpoint of dynamics, that two-generation globular clusters are feasible.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Evolution of Stellar Collision Products in Globular Clusters - II. Off-axis Collision

We continue our exploration of collisionally merged stars in the blue straggler region of the color-magnitude diagram. We report the results of new SPH calculations of parabolic collisions between two main-sequence stars, with the initial structure and composition profiles of the parent stars having been determined from stellar evolution calculations. Parallelization of the SPH code has permitted much higher numerical resolution of the hydrodynamics. We also present evolutionary tracks for the resulting collision products, which emerge as rapidly rotating blue stragglers. The rotating collision products are brighter, bluer and remain on the main sequence longer than their non-rotating counterparts. In addition, they retain their rapid rotation rates throughout their main sequence lifetime. Rotationally-induced mixing strongly affects the evolution of the collision products, although it is not sufficient to mix the entire star. We discuss the implications of these results for studies of blue straggler populations in clusters. This work shows that off-axis collision products cannot become blue stragglers unless they lose a large fraction of their initial angular momentum. The mechanism for this loss is not apparent, although some possibilities are discussed.

0 citations0 reviewsNo code linkedOpen access