Source author record

Mark Gieles

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

45works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

45 published item(s)

preprint2022arXiv

An X-ray quiet black hole born with a negligible kick in a massive binary within the Large Magellanic Cloud

Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5-sigma confidence level. The minimum companion mass implies that it is a black hole. No other X-ray quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the Cosmos.

preprint2022arXiv

On the response of a star cluster to a tidal perturbation

We study the response of star clusters to individual tidal perturbations using controlled $N$-body simulations. We consider perturbations by a moving point mass and by a disc, and vary the duration of the perturbation as well as the cluster density profile. For fast perturbations (i.e. `shocks'), the cluster gains energy in agreement with theoretical predictions in the impulsive limit. For slow disc perturbations, the energy gain is lower, and this has previously been attributed to adiabatic damping. However, the energy gain due to slow perturbations by a point-mass is similar to, or larger than that due to fast shocks, which is not expected because adiabatic damping should be almost independent of the nature of the tides. We show that the geometric distortion of the cluster during slow perturbations is of comparable importance for the energy gain as adiabatic damping, and that the combined effect can qualitatively explain the results. The half-mass radius of the bound stars after a shock increases up to $\sim$7\% for low-concentration clusters, and decreases $\sim$3\% for the most concentrated ones. The fractional mass loss is a non-linear function of the energy gain, and depends on the nature of the tides and most strongly on the cluster density profile, making semi-analytic model predictions for cluster lifetimes extremely sensitive to the adopted density profile.

preprint2021arXiv

MUSE narrow field mode observations of the central kinematics of M15

We present observations of the stellar kinematics of the centre of the core collapsed globular cluster M15 obtained with the MUSE integral field spectrograph on the VLT operating in narrow field mode. Thanks to the use of adaptive optics, we obtain a spatial resolution of 0.1arcsec and are able to reliably measure the radial velocities of 864 stars within 8 arcsec of the centre of M15 thus providing the largest sample of radial velocities ever obtained for the innermost regions of this system. Combined with previous observations of M15 using MUSE in wide field mode and literature data, we find that the central kinematics of M15 are complex with the rotation axis of the core of M15 offset from the rotation axis of the bulk of the cluster. While this complexity has been suggested by previous work, we confirm it at higher significance and in more detail.

preprint2021arXiv

The dynamics of the globular cluster NGC3201 out to the Jacobi radius

As part of a chemo-dynamical survey of five nearby globular clusters with 2dF/AAOmega on the AAT, we have obtained kinematic information for the globular cluster NGC3201. Our new observations confirm the presence of a significant velocity gradient across the cluster which can almost entirely be explained by the high proper motion of the cluster. After subtracting the contribution of this perspective rotation, we found a remaining rotation signal with an amplitude of $\sim1\ km/s$ around a different axis to what we expect from the tidal tails and the potential escapers, suggesting that this rotation is internal and can be a remnant of its formation process. At the outer part, we found a rotational signal that is likely a result from potential escapers. The proper motion dispersion at large radii reported by Bianchini et al. has previously been attributed to dark matter. Here we show that the LOS dispersion between 0.5-1 Jacobi radius is lower, yet above the predictions from an N-body model of NGC3201 that we ran for this study. Based on the simulation, we find that potential escapers cannot fully explain the observed velocity dispersion. We also estimate the effect on the velocity dispersion of different amounts of stellar-mass black holes and unbound stars from the tidal tails with varying escape rates and find that these effects cannot explain the difference between the LOS dispersion and the N-body model. Given the recent discovery of tidal tail stars at large distances from the cluster, a dark matter halo is an unlikely explanation. We show that the effect of binary stars, which is not included in the N-body model, is important and can explain part of the difference in dispersion. We speculate that the remaining difference must be the result of effects not included in the N-body model, such as initial cluster rotation, velocity anisotropy and Galactic substructure.

preprint2021arXiv

Weighing stars from birth to death: mass determination methods across the HRD

The mass of a star is the most fundamental parameter for its structure, evolution, and final fate. It is particularly important for any kind of stellar archaeology and characterization of exoplanets. There exists a variety of methods in astronomy to estimate or determine it. In this review we present a significant number of such methods, beginning with the most direct and model-independent approach using detached eclipsing binaries. We then move to more indirect and model-dependent methods, such as the quite commonly used isochrone or stellar track fitting. The arrival of quantitative asteroseismology has opened a completely new approach to determine stellar masses and to complement and improve the accuracy of other methods. We include methods for different evolutionary stages, from the pre-main sequence to evolved (super)giants and final remnants. For all methods uncertainties and restrictions will be discussed. We provide lists of altogether more than 200 benchmark stars with relative mass accuracies between $[0.3,2]\%$ for the covered mass range of $M\in [0.1,16]\,\msun$, $75\%$ of which are stars burning hydrogen in their core and the other $25\%$ covering all other evolved stages. We close with a recommendation how to combine various methods to arrive at a "mass-ladder" for stars.

preprint2020arXiv

Merger rate of black hole binaries from globular clusters: theoretical error bars and comparison to gravitational wave data from GWTC-2

Black hole binaries formed dynamically in globular clusters are believed to be one of the main sources of gravitational waves in the Universe. Here, we use our new population synthesis code, cBHBd, to determine the redshift evolution of the merger rate density and masses of black hole binaries formed in globular clusters. We simulate $\sim 2$ million models to explore the parameter space that is relevant to real clusters and over all mass scales. We show that when uncertainties on the initial cluster mass function and density are properly taken into account, they become the two dominant factors in setting the theoretical error bars on merger rates. Other model parameters (e.g., natal kicks, black hole masses, metallicity) have virtually no effect on the local merger rate density, although they affect the masses of the merging black holes. Modelling the merger rate density as a function of redshift as $R(z)=R_0(1+z)^κ$ at $z<2$, and marginalizing over uncertainties, we find: $R_0=7.2^{+21.5}_{-5.5}{\rm Gpc^{-3}yr^{-1}}$ and $κ=1.6^{+0.4}_{-0.6}$ ($90\%$ credibility). The rate parameters for binaries that merge inside the clusters are ${R}_{\rm 0,in}=1.6^{+1.9}_{-1.0}{\rm Gpc^{-3}yr^{-1}}$ and $κ_{\rm in}=2.3^{+1.3}_{-1.0}$; $\sim 20\%$ of these form as the result of a gravitational-wave capture, implying that eccentric mergers from globular clusters contribute $\lesssim 0.4 \rm Gpc^{-3}yr^{-1}$ to the local rate. A comparison to the merger rate reported by LIGO-Virgo shows that a scenario in which most of the detected black hole mergers are formed in globular clusters is consistent with current constraints, and requires initial cluster half-mass densities $\gtrsim 10^4 M_\odot \rm pc^{-3}$. Such models also reproduce the inferred primary black hole mass distribution for masses $13-30 M_\odot$, but under-predict the data outside this range.

preprint2020arXiv

Population synthesis of black hole binary mergers from star clusters

Black hole (BH) binary mergers formed through dynamical interactions in dense star clusters are believed to be one of the main sources of gravitational waves for Advanced LIGO and Virgo. Here we present a fast numerical method for simulating the evolution of star clusters with BHs, including a model for the dynamical formation and merger of BH binaries. Our method is based on Hénon's principle of balanced evolution, according to which the flow of energy within a cluster must be balanced by the energy production inside its core. Because the heat production in the core is powered by the BHs, one can then link the evolution of the cluster to the evolution of its BH population. This allows us to construct evolutionary tracks of the cluster properties including its BH population and its effect on the cluster and, at the same time, determine the merger rate of BH binaries as well as their eccentricity distributions. The model is publicly available and includes the effects of a BH mass spectrum, mass-loss due to stellar evolution, the ejection of BHs due to natal and dynamical kicks, and relativistic corrections during binary-single encounters. We validate our method using direct $N$-body simulations, and find it to be in excellent agreement with results from recent Monte Carlo models of globular clusters. This establishes our new method as a robust tool for the study of BH dynamics in star clusters and the modelling of gravitational wave sources produced in these systems. Finally, we compute the rate and eccentricity distributions of merging BH binaries for a wide range of cluster initial conditions, spanning more than two orders of magnitude in mass and radius.

preprint2020arXiv

Star clusters near and far; tracing star formation across cosmic time

Star clusters are fundamental units of stellar feedback and unique tracers of their host galactic properties. In this review, we will first focus on their constituents, i.e.\ detailed insight into their stellar populations and their surrounding ionised, warm, neutral, and molecular gas. We, then, move beyond the Local Group to review star cluster populations at various evolutionary stages, and in diverse galactic environmental conditions accessible in the local Universe. At high redshift, where conditions for cluster formation and evolution are more extreme, we are only able to observe the integrated light of a handful of objects that we believe will become globular clusters. We therefore discuss how numerical and analytical methods, informed by the observed properties of cluster populations in the local Universe, are used to develop sophisticated simulations potentially capable of disentangling the genetic map of galaxy formation and assembly that is carried by globular cluster populations.

preprint2020arXiv

The Physics of Star Cluster Formation and Evolution

Star clusters form in dense, hierarchically collapsing gas clouds. Bulk kinetic energy is transformed to turbulence with stars forming from cores fed by filaments. In the most compact regions, stellar feedback is least effective in removing the gas and stars may form very efficiently. These are also the regions where, in high-mass clusters, ejecta from some kind of high-mass stars are effectively captured during the formation phase of some of the low mass stars and effectively channeled into the latter to form multiple populations. Star formation epochs in star clusters are generally set by gas flows that determine the abundance of gas in the cluster. We argue that there is likely only one star formation epoch after which clusters remain essentially clear of gas by cluster winds. Collisional dynamics is important in this phase leading to core collapse, expansion and eventual dispersion of every cluster. We review recent developments in the field with a focus on theoretical work.

preprint2020arXiv

The star cluster age function in the Galactic disc with Gaia DR2: Fewer old clusters and a low cluster formation efficiency

We perform a systematic reanalysis of the age distribution of Galactic open star clusters. Using a catalogue of homogeneously determined ages for 834 open clusters contained in a 2 kpc cylinder around the Sun and characterised with astrometric and photometric data from the Gaia satellite, we find that it is necessary to revise earlier works that relied on data from the Milky Way Star Cluster survey. After establishing age-dependent completeness limits for our sample, we find that the cluster age function in the range $6.5 < \log t<10$ is compatible with a Schechter-type or broken power-law function, whose parameters we determine by MCMC fitting. Our best-fit values indicate an earlier drop of the age function (by a factor of $2-3$) with respect to the results obtained in the last five years, and are instead more compatible with results obtained in the early 2000s and radio observations of inner-disc clusters. Furthermore, we find a typical destruction time-scale of $\sim1.5$ Gyr for a $10^4\, {\rm M}_{\odot}$ cluster and a present-day cluster-formation rate of $0.55_{-0.15}^{+0.19}$ Myr$^{-1}$kpc$^{-2}$, suggesting that only $16_{-8}^{+11}$ \% of all stars born in the solar neighbourhood form in bound clusters. Accurate cluster-mass measurements are now needed to place more precise constraints on open-cluster formation and evolution models.

preprint2019arXiv

On the black hole content and initial mass function of 47 Tuc

The globular cluster (GC) 47 Tuc has recently been proposed to host an intermediate-mass black hole (IMBH) or a population of stellar-mass black holes (BHs). To shed light on its dark content, we present an application of self-consistent multimass models with a varying mass function and content of stellar remnants, which we fit to various observational constraints. Our best-fitting model successfully matches the observables and correctly predicts the radial distribution of millisecond pulsars and their gravitational accelerations inferred from long-term timing observations. The data favours a population of BHs with a total mass of $430^{+386}_{-301}$ $M_{\odot}$, but the most likely model has very few BHs. Since our models do not include a central IMBH and accurately reproduce the observations, we conclude that there is currently no need to invoke the presence of an IMBH in 47 Tuc. The global present-day mass function inferred is significantly depleted in low-mass stars (power-law slope $α=-0.52^{+0.17}_{-0.16}$). Given the orbit and predicted mass-loss history of this massive GC, the dearth of low-mass stars is difficult to explain with a standard initial mass function (IMF) followed by long-term preferential escape of low-mass stars driven by two-body relaxation, and instead suggests that 47 Tuc may have formed with a bottom-light IMF. We discuss alternative evolutionary origins for the flat mass function and ways to reconcile this with the low BH retention fraction. Finally, by capturing the effect of dark remnants, our method offers a new way to probe the IMF in a GC above the current main-sequence turn-off mass, for which we find a slope of $-2.49\pm0.08$.

preprint2019arXiv

Supermassive stars as the origin of the multiple populations in globular clusters

Globular clusters (GCs) display anomalous light element abundances (HeCNONaMgAl), resembling the yields of hot-hydrogen burning, but there is no consensus yet on the origin of these ubiquitous multiple populations. We present a model in which a super-massive star (SMS, >10^3 Msun) forms via stellar collisions during GC formation and pollutes the intra-cluster medium. The growth of the SMS finds a balance with the wind mass loss rate, such that the SMS can produce a significant fraction of the total GC mass in processed material, thereby overcoming the so-called mass-budget problem that plagues other models. Because of continuous rejuvenation, the SMS acts as a `conveyer-belt' of hot-hydrogen burning yields with (relatively) low He abundances, in agreement with empirical constraints. Additionally, the amount of processed material per unit of GC mass correlates with GC mass, addressing the specific mass budget problem. We discuss uncertainties and tests of this new self-enrichment scenario.

preprint2016arXiv

A family of lowered isothermal models

We present a family of self-consistent, spherical, lowered isothermal models, consisting of one or more mass components, with parameterised prescriptions for the energy truncation and for the amount of radially biased pressure anisotropy. The models are particularly suited to describe the phase-space density of stars in tidally limited, mass-segregated star clusters in all stages of their life-cycle. The models extend a family of isotropic, single-mass models by Gomez-Leyton and Velazquez, of which the well-known Woolley, King and Wilson (in the non-rotating and isotropic limit) models are members. We derive analytic expressions for the density and velocity dispersion components in terms of potential and radius, and introduce a fast model solver in PYTHON (LIMEPY), that can be used for data fitting or for generating discrete samples.

preprint2016arXiv

A stellar-mass black hole population in the globular cluster NGC 6101?

Dalessandro et al. observed a similar distribution for blue straggler stars and main-sequence turn-off stars in the Galactic globular cluster NGC 6101, and interpreted this feature as an indication that this cluster is not mass-segregated. Using direct N-body simulations, we find that a significant amount of mass segregation is expected for a cluster with the mass, radius and age of NGC 6101. Therefore, the absence of mass segregation cannot be explained by the argument that the cluster is not yet dynamically evolved. By varying the retention fraction of stellar-mass black holes, we show that segregation is not observable in clusters with a high black hole retention fraction (>50% after supernova kicks and >50% after dynamical evolution). Yet all model clusters have the same amount of mass segregation in terms of the decline of the mean mass of stars and remnants with distance to the centre. We also discuss how kinematics can be used to further constrain the presence of a stellar-mass black hole population and distinguish it from the effect of an intermediate-mass black hole. Our results imply that the kick velocities of black holes are lower than those of neutron stars. The large retention fraction during its dynamical evolution can be explained if NGC 6101 formed with a large initial radius in a Milky Way satellite.

preprint2016arXiv

If it does not kill them, it makes them stronger: collisional evolution of star clusters with tidal shocks

The radii of young (<100 Myr) star clusters correlate only weakly with their masses. This shallow relation has been used to argue that impulsive tidal perturbations, or `shocks', by passing giant molecular clouds (GMCs) preferentially disrupt low-mass clusters. We show that this mass-radius relation is in fact the result of the combined effect of two-body relaxation and repeated tidal shocks. Clusters in a broad range of environments including those like the solar neighbourhood evolve towards a typical radius of a few parsecs, as observed, independent of the initial radius. This equilibrium mass-radius relation is the result of a competition between expansion by relaxation and shrinking due to shocks. Interactions with GMCs are more disruptive for low-mass clusters, which helps to evolve the globular cluster mass function (GCMF). However, the properties of the interstellar medium in high-redshift galaxies required to establish a universal GCMF shape are more extreme than previously derived, challenging the idea that all GCs formed with the same power-law mass function.

preprint2016arXiv

Properties of the cluster population of NGC 1566 and their implications

We present results of a photometric study into the cluster population of NGC 1566, a nearby grand design spiral galaxy, sampled out to a Galactocentric radius of $\approx 5.5$ kpc. The shape of the mass-limited age distribution shows negligible variation with radial distance from the centre of the galaxy, and demonstrates three separate sections, with a steep beginning, flat middle and steep end. The luminosity function can be approximated by a power law at lower luminosities with evidence of a truncation at higher luminosity. The power law section of the luminosity function of the galaxy is best fitted by an index $\approx -2$, in agreement with other studies, and is found to agree with a model luminosity function, which uses an underlying Schechter mass function. The recovered power law slope of the mass distribution shows a slight steepening as a function of galactocentric distance, but this is within error estimates. It also displays a possible truncation at the high mass end. Additionally, the cluster formation efficiency ($Γ$) and the specific U-band luminosity of clusters ($T_L(U)$) are calculated for NGC 1566 and are consistent with values for similar galaxies. A difference in NGC 1566, however, is that the fairly high star formation rate is in contrast with a low $Σ_{SFR}$ and $Γ$, indicating that $Γ$ can only be said to depend strongly on $Σ_{SFR}$, not the star formation rate.

preprint2016arXiv

Testing lowered isothermal models with direct N-body simulations of globular clusters

Several self-consistent models have been proposed, aiming at describing the phase space distribution of stars in globular clusters. This study explores the ability of the recently proposed LIMEPY models (Gieles & Zocchi) to reproduce the dynamical properties of direct N-body models of a cluster in a tidal field, during its entire evolution. These dynamical models include prescriptions for the truncation and the degree of radially-biased anisotropy contained in the system, allowing us to explore the interplay between the role of anisotropy and tides in various stages of the life of star clusters. We show that the amount of anisotropy in an initially tidally underfilling cluster increases in the pre-collapse phase, and then decreases with time, due to the effect of the external tidal field on its spatial truncation. This is reflected in the correspondent model parameters, and the best-fit models reproduce the main properties of the cluster at all stages of its evolution, except for the phases immediately preceding and following core collapse. We also notice that the best-fit LIMEPY models are significantly different from isotropic King models, especially in the first part of the evolution of the cluster. Our results put limits on the amount of radial anisotropy that can be expected for clusters evolving in a tidal field, which is important to understand other factors that could give rise to similar observational signatures, such as the presence of an intermediate-mass black hole.

preprint2016arXiv

The origin of the Milky Way globular clusters

We present a cosmological zoom-in simulation of a Milky Way-like galaxy used to explore the formation and evolution of star clusters. We investigate in particular the origin of the bimodality observed in the colour and metallicity of globular clusters, and the environmental evolution through cosmic times in the form of tidal tensors. Our results self-consistently confirm previous findings that the blue, metal-poor clusters form in satellite galaxies which are accreted onto the Milky Way, while the red, metal-rich clusters form mostly in situ or, to a lower extent in massive, self-enriched galaxies merging with the Milky Way. By monitoring the tidal fields these populations experience, we find that clusters formed in situ (generally centrally concentrated) feel significantly stronger tides than the accreted ones, both in the present-day, and when averaged over their entire life. Furthermore, we note that the tidal field experienced by Milky Way clusters is significantly weaker in the past than at present-day, confirming that it is unlikely that a power-law cluster initial mass function like that of young massive clusters, is transformed into the observed peaked distribution in the Milky Way with relaxation-driven evaporation in a tidal field.

preprint2015arXiv

A flexible method to evolve collisional systems and their tidal debris in external potentials

We introduce a numerical method to integrate tidal effects on collisional systems, using any definition of the external potential as a function of space and time. Rather than using a linearisation of the tidal field, this new method follows a differential technique to numerically evaluate the tidal acceleration and its time derivative. Theses are then used to integrate the motions of the components of the collisional systems, like stars in star clusters, using a predictor-corrector scheme. The versatility of this approach allows the study of star clusters, including their tidal tails, in complex, multi-components, time-evolving external potentials. The method is implemented in the code nbody6 (Aarseth 2003).

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]

preprint2015arXiv

Biases in the inferred mass-to-light ratio of globular clusters: no need for variations in the stellar mass function

From a study of the integrated light properties of 200 globular clusters (GCs) in M31, Strader et al. found that the mass-to-light ratios are lower than what is expected from simple stellar population (SSP) models with a `canonical' stellar initial mass function (IMF), with the discrepancy being larger at high metallicities. We use dynamical multi-mass models, that include a prescription for equipartition, to quantify the bias in the inferred dynamical mass as the result of the assumption that light follows mass. For a universal IMF and a metallicity dependent present day mass function we find that the inferred mass from integrated light properties systematically under estimates the true mass, and that the bias is more important at high metallicities, as was found for the M31 GCs. We show that mass segregation and a flattening of the mass function have opposing effects of similar magnitude on the mass inferred from integrated properties. This makes the mass-to-light ratio as derived from integrated properties an inadequate probe of the low-mass end of the stellar mass function. There is, therefore, no need for variations in the IMF, nor the need to invoke depletion of low-mass stars, to explain the observations. Finally, we find that the retention fraction of stellar-mass black holes (BHs) is an equally important parameter in understanding the mass segregation bias. We speculatively put forward to idea that kinematical data of GCs can in fact be used to constrain the total mass in stellar-mass BHs in GCs.

preprint2015arXiv

Evolution of star clusters on eccentric orbits

We study the evolution of star clusters on circular and eccentric orbits using direct $N$-body simulations. We model clusters with initially $N=8{\rm k}$ and $N=16{\rm k}$ single stars of the same mass, orbiting around a point-mass galaxy. For each orbital eccentricity that we consider, we find the apogalactic radius at which the cluster has the same lifetime as the cluster with the same $N$ on a circular orbit. We show that then, the evolution of bound particle number and half-mass radius is approximately independent of eccentricity. Secondly, when we scale our results to orbits with the same semi-major axis, we find that the lifetimes are, to first order, independent of eccentricity. When the results of Baumgardt and Makino for a singular isothermal halo are scaled in the same way, the lifetime is again independent of eccentricity to first order, suggesting that this result is independent of the Galactic mass profile. From both sets of simulations we empirically derive the higher order dependence of the lifetime on eccentricity. Our results serve as benchmark for theoretical studies of the escape rate from clusters on eccentric orbits. Finally, our results can be useful for generative models for cold streams and cluster evolution models that are confined to spherical symmetry and/or time-independent tides, such as Fokker-Planck models, Monte Carlo models, and (fast) semi-analytic models.

preprint2015arXiv

Inverting the dynamical evolution of globular clusters: clues to their origin

Scaling relations for globular clusters (GC) differ from scaling relations for pressure supported (elliptical) galaxies. We show that two-body relaxation is the dominant mechanism in shaping the bivariate dependence of density on mass and Galactocentric distance for Milky Way GCs with masses <10^6 Msun, and it is possible, but not required, that GCs formed with similar scaling relations as ultra-compact dwarf galaxies. We use a fast cluster evolution model to fit a parameterised model for the initial properties of Milky Way GCs to the observed present-day properties. The best-fit cluster initial mass function is substantially flatter (power-law index alpha =- 0.6+/-0.2) than what is observed for young massive clusters (YMCs) forming in the nearby Universe (alpha =~-2). A slightly steeper CIMF is allowed when considering the metal-rich GCs separately (alpha =~-1.2+/-0.4$). If stellar mass loss and two-body relaxation in the Milky Way tidal field are the dominant disruption mechanisms, then GCs formed differently from YMCs.

preprint2015arXiv

Multiple populations in globular clusters: the distinct kinematic imprints of different formation scenarios

Several scenarios have been proposed to explain the presence of multiple stellar populations in globular clusters. Many of them invoke multiple generations of stars to explain the observed chemical abundance anomalies, but it has also been suggested that self-enrichment could occur via accretion of ejecta from massive stars onto the circumstellar disc of low-mass pre-main sequence stars. These scenarios imply different initial conditions for the kinematics of the various stellar populations. Given some net angular momentum initially, models for which a second generation forms from gas that collects in a cooling flow into the core of the cluster predict an initially larger rotational amplitude for the polluted stars compared to the pristine stars. This is opposite to what is expected from the accretion model, where the polluted stars are the ones crossing the core and are on preferentially radial (low-angular momentum) orbits, such that their rotational amplitude is lower. Here we present the results of a suite of $N$-body simulations with initial conditions chosen to capture the distinct kinematic properties of these pollution scenarios. We show that initial differences in the kinematics of polluted and pristine stars can survive to the present epoch in the outer parts of a large fraction of Galactic globular clusters. The differential rotation of pristine and polluted stars is identified as a unique kinematic signature that could allow us to distinguish between various scenarios, while other kinematic imprints are generally very similar from one scenario to the other.

preprint2015arXiv

The effect of secular galactic growth on the evolution of star clusters

The growth of galaxies through adiabatic accretion of dark matter is one of the main drivers of galaxy evolution. By isolating it from other processes like mergers, we analyse how it affects the evolution of star clusters. Our study comprises a fast and approximate exploration of the orbital and intrinsic cluster parameter space, and more detailed monitoring of their evolution, through N-body simulations for a handful of cases. We find that the properties of present-day star clusters and their tidal tails differ very little, whether the clusters are embedded in a growing galactic halo for 12 Gyr, or in a static one.

preprint2014arXiv

A prescription and fast code for the long-term evolution of star clusters - III. Unequal masses and stellar evolution

We present a new version of the fast star cluster evolution code Evolve Me A Cluster of StarS (EMACSS). While previous versions of EMACSS reproduced clusters of single-mass stars, this version models clusters with an evolving stellar content. Stellar evolution dominates early evolution, and leads to: (1) reduction of the mean mass of stars due to the mass loss of high-mass stars; (2) expansion of the half-mass radius; (3) for (nearly) Roche Volume filling clusters, the induced escape of stars. Once sufficient relaxation has occurred (~ 10 relaxation times-scales), clusters reach a second, 'balanced' state whereby the core releases energy as required by the cluster as a whole. In this state: (1) stars escape due to tidal effects faster than before balanced evolution; (2) the half-mass radius expands or contracts depending on the Roche volume filling factor; and (3) the mean mass of stars increases due to the preferential ejection of low-mass stars. We compare the EMACSS results of several cluster properties against N-body simulations of clusters spanning a range of initial number of stars, mass, half-mass radius, and tidal environments, and show that our prescription accurately predicts cluster evolution for this database. Finally, we consider applications for EMACSS, such as studies of galactic globular cluster populations in cosmological simulations.

preprint2014arXiv

The inefficiency of satellite accretion in forming extended star clusters

The distinction between globular clusters and dwarf galaxies has been progressively blurred by the recent discoveries of several extended star clusters, with size (20-30 pc) and luminosity (-6 < Mv < -2) comparable to the one of faint dwarf spheroidals. In order to explain their sparse structure, it has been suggested that they formed as star clusters in dwarf galaxy satellites that later accreted onto the Milky Way. If these clusters form in the centre of dwarf galaxies, they evolve in a tidally-compressive environment where the contribution of the tides to the virial balance can become significant, and lead to a super-virial state and subsequent expansion of the cluster, once removed. Using N-body simulations, we show that a cluster formed in such an extreme environment undergoes a sizable expansion, during the drastic variation of the external tidal field due to the accretion process. However, we show that the expansion due to the removal of the compressive tides is not enough to explain the observed extended structure, since the stellar systems resulting from this process are always more compact than the corresponding clusters that expand in isolation due to two-body relaxation. We conclude that an accreted origin of extended globular clusters is unlikely to explain their large spatial extent, and rather favor the hypothesis that such clusters are already extended at the stage of their formation.

preprint2013arXiv

A prescription and fast code for the long-term evolution of star clusters - II. Unbalanced and core evolution

We introduce version two of the fast star cluster evolution code Evolve Me A Cluster of StarS (EMACSS). The first version (Alexander & Gieles) assumed that cluster evolution is balanced for the majority of the life-cycle, meaning that the rate of energy generation in the core of the cluster equals the diffusion rate of energy by two-body relaxation, which makes the code suitable for modelling clusters in weak tidal fields. In this new version we extend the model to include an unbalanced phase of evolution to describe the pre-collapse evolution and the accompanying escape rate such that clusters in strong tidal fields can also be modelled. We also add a prescription for the evolution of the core radius and density and a related cluster concentration parameter. The model simultaneously solves a series of first-order ordinary differential equations for the rate of change of the core radius, half-mass radius and the number of member stars N. About two thousand integration steps in time are required to solve for the entire evolution of a star cluster and this number is approximately independent of N. We compare the model to the variation of these parameters following from a series of direct N-body calculations of single-mass clusters and find good agreement in the evolution of all parameters. Relevant time-scales, such as the total lifetimes and core collapse times, are reproduced with an accuracy of about 10% for clusters with various initial half-mass radii (relative to their Jacobi radii) and a range of different initial N up to N = 65536. We intend to extend this framework to include more realistic initial conditions, such as a stellar mass spectrum and mass loss from stars. The EMACSS code can be used in star cluster population studies and in models that consider the co-evolution of (globular) star clusters and large scale structures.

preprint2013arXiv

Constraining the initial conditions of globular clusters using their radius distribution

Studies of extra-galactic globular clusters have shown that the peak size of the globular cluster (GC) radius distribution (RD) depends only weakly on galactic environment, and can be used as a standard ruler. We model RDs of GC populations using a simple prescription for a Hubble time of relaxation driven evolution of cluster mass and radius, and explore the conditions under which the RD can be used as a standard ruler. We consider a power-law cluster initial mass function (CIMF) with and without an exponential truncation, and focus in particular on a flat and a steep CIMF (power-law indices of 0 and -2, respectively). For the initial half-mass radii at birth we adopt either Roche-lobe filling conditions ('filling',meaning that the ratio of half-mass to Jacobi radius is approximately rh/rJ ~ 0.15) or strongly Roche-lobe under-filling conditions ('under-filling', implying that initially rh/rJ << 0.15). Assuming a constant orbital velocity about the galaxy centre we find for a steep CIMF that the typical half-light radius scales with galactocentric radius RG as RG^1/3. This weak scaling is consistent with observations, but this scenario has the (well known) problem that too many low-mass clusters survive. A flat CIMF with 'filling' initial conditions results in the correct mass function at old ages, but with too many large (massive) clusters at large RG. An 'underfilling' GC population with a flat CIMF also results in the correct mass function, and can also successfully reproduce the shape of the RD, with a peak size that is (almost) independent of RG. In this case, the peak size depends (almost) only on the peak mass of the GC mass function. The (near) universality of the GC RD is therefore because of the (near) universality of the CIMF. There are some extended GCs in the outer halo of the Milky Way that cannot be explained by this model.

preprint2013arXiv

The evolution of the global stellar mass function of star clusters: an analytic description

The evolution of the global stellar mass function (MF) of star clusters is studied based on a large set of N-body simulations of clusters with a range of initial masses, initial concentrations, in circular or elliptical orbits in different tidal environments. Models with and without initial mass segregation are included. The depletion of low mass stars in initially Roche-volume (tidal) filling clusters starts typically on a time scale of the order of the core collapse time. In clusters that are initially underfilling their Roche-volume it takes longer because the clusters have to expand to their tidal radii before dynamical mass loss becomes important. We introduce the concept of the differential mass function (DMF), which describes the changes with respect to the initial mass function (IMF). We show that the evolution of the DMF can be described by a set of very simple analytic expressions that are valid for a wide range of initial cluster parameters and for different IMFs. The agreement between this description and the models is very good, except for initially Roche-volume underfilling clusters that are severely mass segregated.

preprint2013arXiv

The mass and radius evolution of globular clusters in tidal fields

We present a simple theory for the evolution of initially compact clusters in a tidal field. The fundamental ingredient of the model is that a cluster conducts a constant fraction of its own energy through the half-mass radius by two-body interactions every half-mass relaxation time. This energy is produced in a self-regulative way in the core by an (unspecified) energy source. We find that the half-mass radius increases during the first part (roughly half) of the evolution and decreases in the second half, while the escape rate is constant and set by the tidal field. We present evolutionary tracks and isochrones for clusters in terms of cluster half-mass density, cluster mass and galacto-centric radius. We find substantial agreement between model isochrones and Milky Way globular cluster parameters, which suggests that there is a balance between the flow of energy and the central energy production for almost all globular clusters. We also find that the majority of the globular clusters are still expanding towards their tidal radius. Finally, a fast code for cluster evolution is presented.

preprint2013arXiv

The role of galaxy mergers on the evolution of star clusters

Interacting galaxies favor the formation of star clusters but are also suspected to affect their evolution through an intense and rapidly varying tidal field. Treating this complex behaviour remains out-of-reach of (semi-)analytical studies. By computing the tidal field from galactic models and including it into star-by-star N-body simulations of star clusters, we monitor the structure and mass evolution of a population of clusters in a galaxy major merger, taking the Antennae galaxies (NGC 4038/39) as a prototype. On the long timescale (~ 10^9 yr), the merger only indirectly affects the evolution of clusters by modifying their orbits in or around the galaxies: the mass-loss of clusters in the merger remnant is faster, while clusters ejected in the tidal debris survive much longer, compared to in an isolated galaxy. The tidal perturbations of the galactic collisions themselves are too short lived and not strong enough to significantly influence the structure and dissolution of realistically dense/massive star clusters.

preprint2012arXiv

A prescription and fast code for the long-term evolution of star clusters

We introduce the star cluster evolution code Evolve Me A Cluster of StarS (EMACSS), a simple yet physically motivated computational model that describes the evolution of some fundamental properties of star clusters in static tidal fields. We base our prescription upon the flow of energy within the cluster, which is a constant fraction of the total energy per half-mass relaxation time. According to Henon's predictions, this flow is independent of the precise mechanisms for energy production within the core, and therefore does not require a complete description of the many-body interactions therein. For a cluster of equal-mass stars, we thence use dynamical theory and analytic descriptions of escape mechanisms to construct a series of coupled differential equations expressing the time evolution of cluster mass and radius. These equations are numerically solved using a fourth-order Runge-Kutta integration kernel, and the results benchmarked against a data base of direct N-body simulations. We use simulations containing a modest initial number of stars (1024 <= N <= 65 536) and point-mass tidal fields of various strengths. Our prescription is publicly available and reproduces the N-body results to within ~10 per cent accuracy for the entire post-collapse evolution of star clusters.

preprint2012arXiv

Do all stars in the solar neighbourhood form in clusters? A cautionary note on the use of the distribution of surface densities

Bressert et al. recently showed that the surface density distribution of low-mass, young stellar objects (YSOs) in the solar neighbourhood is approximately log-normal. The authors conclude that the star formation process is hierarchical and that only a small fraction of stars form in dense star clusters. Here, we show that the peak and the width of the density distribution are also what follow if all stars form in bound clusters which are not significantly affected by the presence of gas and expand by two-body relaxation. The peak of the surface density distribution is simply obtained from the typical ages (few Myr) and cluster membership number (few hundred) typifying nearby star-forming regions. This result depends weakly on initial cluster sizes, provided that they are sufficiently dense (initial half mass radius of <0.3 pc) for dynamical evolution to be important at an age of a few Myr. We conclude that the degeneracy of the YSO surface density distribution complicates its use as a diagnostic of the stellar formation environment.

preprint2012arXiv

Mass loss of stars in star clusters: an energy source for dynamical evolution

Dense star clusters expand until their sizes are limited by the tidal field of their host galaxy. During this expansion phase the member stars evolve and lose mass. We show that for clusters with short initial relaxation time scales (<~100 Myr) the dynamical expansion is largely powered by mass loss from stars in the core, but happens on a relaxation time scale. That is, the energy release following stellar mass loss is in balance with the amount of energy that is transported outward by two-body relaxation.

preprint2011arXiv

Dynamical evolution of stellar clusters

The evolution of star clusters is determined by several internal and external processes. Here we focus on two dominant internal effects, namely energy exchange between stars through close encounters (two-body relaxation) and mass-loss of the member stars through stellar winds and supernovae explosions. Despite the fact that the former operates on the relaxation timescale of the cluster and the latter on a stellar evolution timescale, these processes work together in driving a nearly self-similar expansion, without forming (hard) binaries. Low-mass clusters expand more, such that after some time the radii of clusters depend very little on their masses, even if all clusters have the same (surface) density initially. Throughout it is assumed that star clusters are in virial equilibrium and well within their tidal boundary shortly after formation, motivated by observations of young (few Myrs) clusters. We start with a discussion on how star clusters can be distinguished from (unbound) associations at these young ages.

preprint2011arXiv

Evolution of star clusters in arbitrary tidal fields

We present a novel and flexible tensor approach to computing the effect of a time-dependent tidal field acting on a stellar system. The tidal forces are recovered from the tensor by polynomial interpolation in time. The method has been implemented in a direct-summation stellar dynamics integrator (NBODY6) and test-proved through a set of reference calculations: heating, dissolution time and structural evolution of model star clusters are all recovered accurately. The tensor method is applicable to arbitrary configurations, including the important situation where the background potential is a strong function of time. This opens up new perspectives in stellar population studies reaching to the formation epoch of the host galaxy or galaxy cluster, as well as for star-burst events taking place during the merger of large galaxies. A pilot application to a star cluster in the merging galaxies NGC 4038/39 (the Antennae) is presented.

preprint2011arXiv

Outer density profiles of 19 Galactic globular clusters from deep and wide-field imaging

Using deep photometric data from WFC@INT and WFI@ESO2.2m we measure the outer number density profiles of 19 stellar clusters located in the inner region of the Milky Way halo (within a Galactocentric distance range of 10-30 kpc) in order to assess the impact of internal and external dynamical processes on the spatial distribution of stars. Adopting power-law fitting templates, with index $-γ$ in the outer region, we find that the clusters in our sample can be divided in two groups: a group of massive clusters ($ \ge 10^5 $ M_sun) that has relatively flat profiles with $2.5 < γ< 4$ and a group of low-mass clusters ($ \le 10^5 $ M_sun), with steep profiles ($γ> 4$) and clear signatures of interaction with the Galactic tidal field. We refer to these two groups as 'tidally unaffected' and 'tidally affected', respectively. Our results also show a clear trend between the slope of the outer parts and the half-mass density of these systems, which suggests that the outer density profiles may retain key information on the dominant processes driving the dynamical evolution of Globular Clusters.

preprint2011arXiv

The life cycle of star cluster in a tidal field

The evolution of globular clusters due to 2-body relaxation results in an outward flow of energy and at some stage all clusters need a central energy source to sustain their evolution. Henon provided the insight that we do not need to know the details of the energy production in order to understand the relaxation-driven evolution of the cluster, at least outside the core. He provided two self-similar solutions for the evolution of clusters based on the view that the cluster as a whole determines the amount of energy that is produced in the core: steady expansion for isolated clusters, and homologous contraction for clusters evaporating in a tidal field. We combine these models: the half-mass radius increases during the first half of the evolution, and decreases in the second half; while the escape rate approaches a constant value set by the tidal field. We refer to these phases as `expansion dominated' and `evaporation dominated'. These simple analytical solutions immediately allow us to construct evolutionary tracks and isochrones in terms of cluster half-mass density, cluster mass and galacto-centric radius. From a comparison to the Milky Way globular clusters we find that roughly 1/3 of them are in the second, evaporation-dominated phase and for these clusters the density inside the half-mass radius varies with the galactocentric distance R as rho_h ~ 1/R^2. The remaining 2/3 are still in the first, expansion-dominated phase and their isochrones follow the environment-independent scaling rho_h ~ M^2; that is, a constant relaxation time-scale. We find substantial agreement between Milky Way globular cluster parameters and the isochrones, which suggests that there is, as Henon suggested, a balance between the flow of energy and the central energy production for almost all globular clusters.

preprint2010arXiv

GLIMPSE-CO1: the most massive intermediate-age stellar cluster in the Galaxy

The stellar cluster GLIMPSE-C01 is a dense stellar system located in the Galactic Plane. Though often referred to in the literature as an old globular cluster traversing the Galactic disk, previous observations do not rule out that it is an intermediate age (less than a few Gyr) disk-borne cluster. Here, we present high-resolution near-infrared spectroscopy of over 50 stars in the cluster. We find an average radial velocity is consistent with being part of the disk, and determine the cluster's dynamical mass to be (8 \pm 3)x10^4 Msun. Analysis of the cluster's M/L ratio, the location of the Red Clump, and an extremely high stellar density, all suggest an age of 400-800Myr for GLIMPSE-C01, much lower than for a typical globular cluster. This evidence therefore leads us to conclude that GLIMPSE-C01 is part of the disk population, and is the most massive Galactic intermediate-age cluster discovered to date.

preprint2010arXiv

On the mass-radius relation of hot stellar systems

Most globular clusters have half-mass radii of a few pc with no apparent correlation with their masses. This is different from elliptical galaxies, for which the Faber-Jackson relation suggests a strong positive correlation between mass and radius. Objects that are somewhat in between globular clusters and low-mass galaxies, such as ultra-compact dwarf galaxies, have a mass-radius relation consistent with the extension of the relation for bright ellipticals. Here we show that at an age of 10 Gyr a break in the mass-radius relation at ~10^6 Msun is established because objects below this mass, i.e. globular clusters, have undergone expansion driven by stellar evolution and hard binaries. From numerical simulations we find that the combined energy production of these two effects in the core comes into balance with the flux of energy that is conducted across the half-mass radius by relaxation. An important property of this `balanced' evolution is that the cluster half-mass radius is independent of its initial value and is a function of the number of bound stars and the age only. It is therefore not possible to infer the initial mass-radius relation of globular clusters and we can only conclude that the present day properties are consistent with the hypothesis that all hot stellar systems formed with the same mass-radius relation and that globular clusters have moved away from this relation because of a Hubble time of stellar and dynamical evolution.

preprint2010arXiv

The distinction between star clusters and associations

In Galactic studies a distinction is made between (open) star clusters and associations. For barely resolved objects at a distance of several Mpc this distinction is not trivial to make. Here we provide an objective definition by comparing the age of the stars to the crossing time of nearby stellar agglomerates. We find that a satisfactory separation can be made where this ratio equals unity. Stellar agglomerates for which the age of the stars exceeds the crossing time are bound, and are referred to as star clusters. Alternatively, those for which the crossing time exceeds the stellar age are unbound and are referred to as associations. This definition is useful whenever reliable measurements for the mass, radius and age are available.

preprint2010arXiv

Young massive star clusters

Young massive clusters are dense aggregates of young stars that form the fundamental building blocks of galaxies. Several examples exist in the Milky Way Galaxy and the Local Group, but they are particularly abundant in starburst and interacting galaxies. The few young massive clusters that are close enough to resolve are of prime interest for studying the stellar mass function and the ecological interplay between stellar evolution and stellar dynamics. The distant unresolved clusters may be effectively used to study the star-cluster mass function, and they provide excellent constraints on the formation mechanisms of young cluster populations. Young massive clusters are expected to be the nurseries for many unusual objects, including a wide range of exotic stars and binaries. So far only a few such objects have been found in young massive clusters, although their older cousins, the globular clusters, are unusually rich in stellar exotica. In this review we focus on star clusters younger than $\sim100$ Myr, more than a few current crossing times old, and more massive than $\sim10^4$ \Msun, irrespective of cluster size or environment. We describe the global properties of the currently known young massive star clusters in the Local Group and beyond, and discuss the state of the art in observations and dynamical modeling of these systems. In order to make this review readable by observers, theorists, and computational astrophysicists, we also review the cross-disciplinary terminology.

preprint2009arXiv

Evidence for two populations of Galactic globular clusters from the ratio of their half-mass to Jacobi radii

We investigate the ratio between the half-mass radii r_h of Galactic globular clusters and their Jacobi radii r_J given by the potential of the Milky Way and show that clusters with galactocentric distances R_{GC}>8 kpc fall into two distinct groups: one group of compact, tidally-underfilling clusters with r_h/r_J<0.05 and another group of tidally filling clusters which have 0.1 < r_h/r_J<0.3. We find no correlation between the membership of a particular cluster to one of these groups and its membership in the old or younger halo population. Based on the relaxation times and orbits of the clusters, we argue that compact clusters and most clusters in the inner Milky Way were born compact with half-mass radii r_h < 1 pc. Some of the tidally-filling clusters might have formed compact as well, but the majority likely formed with large half-mass radii. Galactic globular clusters therefore show a similar dichotomy as was recently found for globular clusters in dwarf galaxies and for young star clusters in the Milky Way. It seems likely that some of the tidally-filling clusters are evolving along the main sequence line of clusters recently discovered by Kuepper et al. (2008) and are in the process of dissolution.

preprint2006arXiv

Clusters in the solar neighbourhood: how are they destroyed?

We predict the survival time of initially bound star clusters in the solar neighbourhood taking into account: (1) stellar evolution, (2) tidal stripping, (3) shocking by spiral arms and (4) encounters with giant molecular clouds. We find that the predicted dissolution time is t_dis= 1.7 (M_i/10^4 M_sun)^0.67 Gyr for clusters in the mass range of 10^2 < M_i < 10^5 M_sun, where M_i is the initial mass of the cluster.. The resulting predicted shape of the logarithmic age distribution agrees very well with the empirical one, derived from a complete sample of clusters in the solar neighbourhood within 600 pc. The required scaling factor implies a star formation rate of 400 M_sun/Myr within 600 pc from the Sun or a surface formation rate of 3.5 10^-10 M_sun/(yr pc^2) for stars in bound clusters with an initial mass in the range of 10^2 to 3 10^4 M_sun.