Alice Zocchi

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.

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 per cent after supernova kicks and >50 per cent 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.

The properties of energetically unbound stars in stellar clusters

Several Milky Way star clusters show a roughly flat velocity dispersion profile at large radii, which is not expected from models with a tidal cut-off energy. Possible explanations for this excess velocity include: the effects of a dark matter halo, modified gravity theories and energetically unbound stars inside of clusters. These stars are known as potential escapers (PEs) and can exist indefinitely within clusters which are on circular orbits. Through a series of N-body simulations of star cluster systems, where we vary the galactic potential, orbital eccentricity and stellar mass function, we investigate the properties of the PEs and their effects on the kinematics. We derive a prediction for the scaling of the velocity dispersion at the Jacobi surface due to PEs, as a function of cluster mass, angular velocity of the cluster orbit, and slope of the mass profile of the host galaxy. We see a tentative signal of the mass and orbital velocity dependence in kinematic data of globular clusters from literature. We also find that the fraction of PEs depends sensitively on the galactic mass profile, reaching as high as 40% in the cusp of a Navarro-Frenk-White profile and as the velocity anisotropy also depends on the slope of the galactic mass profile, we conclude that PEs provide an independent way of inferring the properties of the dark matter mass profile at the galactic radius of (globular) clusters in the Gaia-era

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

Radial anisotropy in ω Cen limiting the room for an intermediate-mass black hole

Finding an intermediate-mass black hole (IMBH) in a globular cluster (or proving its absence) would provide valuable insights into our understanding of galaxy formation and evolution. However, it is challenging to identify a unique signature of an IMBH that cannot be accounted for by other processes. Observational claims of IMBH detection are indeed often based on analyses of the kinematics of stars in the cluster core, the most common signature being a rise in the velocity dispersion profile towards the centre of the system. Unfortunately, this IMBH signal is degenerate with the presence of radially-biased pressure anisotropy in the globular cluster. To explore the role of anisotropy in shaping the observational kinematics of clusters, we analyse the case of ω Cen by comparing the observed profiles to those calculated from the family of LIMEPY models, that account for the presence of anisotropy in the system in a physicallymotivated way. The best-fit radially anisotropicmodels reproduce the observational profiles well, and describe the central kinematics as derived from Hubble Space Telescope proper motions without the need for an IMBH.

Concurrent formation of supermassive stars and globular clusters : implications for early self-enrichment

We present a model for the concurrent formation of globular clusters (GCs) and supermassive stars (SMSs, ≳103 M⊙) to address the origin of the HeCNONaMgAl abundance anomalies in GCs. GCs form in converging gas flows and accumulate low-angular momentum gas, which accretes on to protostars. This leads to an adiabatic contraction of the cluster and an increase of the stellar collision rate. A SMS can form via runaway collisions if the cluster reaches sufficiently high density before two-body relaxation halts the contraction. This condition is met if the number of stars ≳106 and the gas accretion rate ≳105 M⊙ Myr−1, reminiscent of GC formation in high gas-density environments, such as – but not restricted to – the early Universe. The strong SMS wind mixes with the inflowing pristine gas, such that the protostars accrete diluted hot-hydrogen burning yields of the SMS. Because of continuous rejuvenation, the amount of processed material liberated by the SMS can be an order of magnitude higher than its maximum mass. This ‘conveyor-belt’ production of hot-hydrogen burning products provides a solution to the mass budget problem that plagues other scenarios. Additionally, the liberated material is mildly enriched in helium and relatively rich in other hot-hydrogen burning products, in agreement with abundances of GCs today. Finally, we find a super-linear scaling between the amount of processed material and cluster mass, providing an explanation for the observed increase of the fraction of processed material with GC mass. We discuss open questions of this new GC enrichment scenario and propose observational tests.

Mass models of NGC 6624 without an intermediate-mass black hole

An intermediate-mass black hole (IMBH) was recently reported to reside in the centre of the Galactic globular cluster (GC) NGC 6624, based on timing observations of a millisecond pulsar (MSP) located near the cluster centre in projection. We present dynamical models with multiple mass components of NGC 6624 – without an IMBH – which successfully describe the surface brightness profile and proper motion kinematics from the Hubble Space Telescope (HST) and the stellar mass function at different distances from the cluster centre. The maximum line-of-sight acceleration at the position of the MSP accommodates the inferred acceleration of the MSP, as derived from its first period derivative. With discrete realisations of the models we show that the higher-order period derivatives – which were previously used to derive the IMBH mass – are due to passing stars and stellar remnants, as previously shown analytically in literature. We conclude that there is no need for an IMBH to explain the timing observations of this MSP.

Testing lowered isothermal models with direct N-body simulations of globular clusters – II. Multimass models

Lowered isothermal models, such as the multimass Michie–King models, have been successful in describing observational data of globular clusters. In this study, we assess whether such models are able to describe the phase space properties of evolutionary N-body models. We compare the multimass models as implemented in LIMEPY (Gieles & Zocchi) to N-body models of star clusters with different retention fractions for the black holes and neutron stars evolving in a tidal field. We find that multimass models successfully reproduce the density and velocity dispersion profiles of the different mass components in all evolutionary phases and for different remnants retention. We further use these results to study the evolution of global model parameters. We find that over the lifetime of clusters, radial anisotropy gradually evolves from the low- to the high-mass components and we identify features in the properties of observable stars that are indicative of the presence of stellar-mass black holes. We find that the model velocity scale depends on mass as m−δ, with δ ≃ 0.5 for almost all models, but the dependence of central velocity dispersion on m can be shallower, depending on the dark remnant content, and agrees well with that of the N-body models. The reported model parameters, and correlations amongst them, can be used as theoretical priors when fitting these types of mass models to observational data.

Asking gender questions: Results from a survey of gender and question asking among UK Astronomers at NAM2014

Jonathan Pritchard, Karen Masters, James Allen, Filippo Contenta, Leo Huckvale, Stephen Wilkins and Alice Zocchi report on a survey of the gender of astronomers attending and asking questions at this years UK National Astronomy Meeting.

Asking gender questions

Jonathan Pritchard, Karen Masters, James Allen, Filippo Contenta, Leo Huckvale, Stephen Wilkins and Alice Zocchi report on a survey of the gender of astronomers attending and asking questions at this years UK National Astronomy Meeting.