Showing 1–49 of 49 results for author: Webb, J J

Stream Members Only: Data-Driven Characterization of Stellar Streams with Mixture Density Networks

Authors: Nathaniel Starkman, Jacob Nibauer, Jo Bovy, Jeremy J. Webb, Kiyan Tavangar, Adrian Price-Whelan, Ana Bonaca

Abstract: Stellar streams are sensitive probes of the Milky Way's gravitational potential. The mean track of a stream constrains global properties of the potential, while its fine-grained surface density constrains galactic substructure. A precise characterization of streams from potentially noisy data marks a crucial step in inferring galactic structure, including the dark matter, across orders of magnitud… ▽ More Stellar streams are sensitive probes of the Milky Way's gravitational potential. The mean track of a stream constrains global properties of the potential, while its fine-grained surface density constrains galactic substructure. A precise characterization of streams from potentially noisy data marks a crucial step in inferring galactic structure, including the dark matter, across orders of magnitude in mass scales. Here we present a new method for constructing a smooth probability density model of stellar streams using all of the available astrometric and photometric data. To characterize a stream's morphology and kinematics, we utilize mixture density networks to represent its on-sky track, width, stellar number density, and kinematic distribution. We model the photometry for each stream as a single-stellar population, with a distance track that is simultaneously estimated from the stream's inferred distance modulus (using photometry) and parallax distribution (using astrometry). We use normalizing flows to characterize the distribution of background stars. We apply the method to the stream GD-1, and the tidal tails of Palomar 5. For both streams we obtain a catalog of stellar membership probabilities that are made publicly available. Importantly, our model is capable of handling data with incomplete phase-space observations, making our method applicable to the growing census of Milky Way stellar streams. When applied to a population of streams, the resulting membership probabilities from our model form the required input to infer the Milky Way's dark matter distribution from the scale of the stellar halo down to subhalos. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: 35 pages, 13 figures, 4 tables, fully open-source and reproducible using ShowYourWork

Hierarchical Bayesian Inference of Globular Cluster Properties

Authors: Robin Y. Wen, Joshua S. Speagle, Jeremy J. Webb, Gwendolyn M. Eadie

Abstract: We present a hierarchical Bayesian inference approach to estimating the structural properties and the phase space center of a globular cluster (GC) given the spatial and kinematic information of its stars based on lowered isothermal cluster models. As a first step towards more realistic modelling of GCs, we built a differentiable, accurate emulator of the lowered isothermal distribution function u… ▽ More We present a hierarchical Bayesian inference approach to estimating the structural properties and the phase space center of a globular cluster (GC) given the spatial and kinematic information of its stars based on lowered isothermal cluster models. As a first step towards more realistic modelling of GCs, we built a differentiable, accurate emulator of the lowered isothermal distribution function using interpolation. The reliable gradient information provided by the emulator allows the use of Hamiltonian Monte Carlo methods to sample large Bayesian models with hundreds of parameters, thereby enabling inference on hierarchical models. We explore the use of hierarchical Bayesian modelling to address several issues encountered in observations of GC including an unknown GC center, incomplete data, and measurement errors. Our approach not only avoids the common technique of radial binning but also incorporates the aforementioned uncertainties in a robust and statistically consistent way. Through demonstrating the reliability of our hierarchical Bayesian model on simulations, our work lays out the foundation for more realistic and complex modelling of real GC data. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 16 pages, 12 figures, and 2 tables

A catalogue of Galactic GEMS: Globular cluster Extra-tidal Mock Stars

Authors: Steffani M. Grondin, Jeremy J. Webb, James M. M. Lane, Joshua S. Speagle, Nathan W. C. Leigh

Abstract: This work presents the Globular cluster Extra-tidal Mock Star (GEMS) catalogue of extra-tidal stars and binaries created via three-body dynamical encounters in globular cluster cores. Using the particle-spray code Corespray, we sample N=50,000 extra-tidal stars and escaped recoil binaries for 159 Galactic globular clusters. Sky positions, kinematics, stellar properties and escape information are p… ▽ More This work presents the Globular cluster Extra-tidal Mock Star (GEMS) catalogue of extra-tidal stars and binaries created via three-body dynamical encounters in globular cluster cores. Using the particle-spray code Corespray, we sample N=50,000 extra-tidal stars and escaped recoil binaries for 159 Galactic globular clusters. Sky positions, kinematics, stellar properties and escape information are provided for all simulated stars. Stellar orbits are integrated in seven different static and time-varying Milky Way gravitational potential models where the structure of the disc, perturbations from the Large Magellanic Cloud and the mass and sphericity of the Milky Way's dark matter halo are all investigated. We find that the action coordinates of the mock extra-tidal stars are largely Galactic model independent, where minor offsets and broadening of the distributions between models are likely due to interactions with substructure. Importantly, we also report the first evidence for stellar stream contamination by globular cluster core stars and binaries for clusters with pericentre radii larger than five kiloparsecs. Finally, we provide a quantitative tool that uses action coordinates to match field stars to host clusters with probabilities. Ultimately, combining data from the GEMS catalogue with information of observed stars will allow for association of extra-tidal field stars with any Galactic globular cluster; a requisite tool for understanding population-level dynamics and evolution of clusters in the Milky Way. △ Less

Submitted 16 January, 2024; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: 23 pages, 12 figures, 1 table. Accepted for publication in Monthly Notices of the Royal Astronomical Society. For the associated GEMS catalogue, please visit https://zenodo.org/records/8436703

The dominant mechanism(s) for populating the outskirts of star clusters with neutron star binaries

Authors: Nathan W. C. Leigh, Claire S. Ye, Steffani M. Grondin, Giacomo Fragione, Jeremy J. Webb, Craig O. Heinke

Abstract: It has been argued that heavy binaries composed of neutron stars (NSs) and millisecond pulsars (MSPs) can end up in the outskirts of star clusters via an interaction with a massive black hole (BH) binary expelling them from the core. We argue here, however, that this mechanism will rarely account for such observed objects. Only for primary masses $\lesssim$ 100 M$_{\odot}$ and a narrow range of or… ▽ More It has been argued that heavy binaries composed of neutron stars (NSs) and millisecond pulsars (MSPs) can end up in the outskirts of star clusters via an interaction with a massive black hole (BH) binary expelling them from the core. We argue here, however, that this mechanism will rarely account for such observed objects. Only for primary masses $\lesssim$ 100 M$_{\odot}$ and a narrow range of orbital separations should a BH-BH binary be both dynamically hard and produce a sufficiently low recoil velocity to retain the NS binary in the cluster. Hence, BH binaries are in general likely to eject NSs from clusters. We explore several alternative mechanisms that would cause NS/MSP binaries to be observed in the outskirts of their host clusters after a Hubble time. The most likely mechanism is a three-body interaction involving the NS/MSP binary and a normal star. We compare to Monte Carlo simulations of cluster evolution for the globular clusters NGC 6752 and 47 Tuc, and show that the models not only confirm that normal three-body interactions involving all stellar-mass objects are the dominant mechanism for putting NS/MSP binaries into the cluster outskirts, they also reproduce the observed NS/MSP binary radial distributions without needing to invoke the presence of a massive BH binary. Higher central densities and an episode of core-collapse can broaden the radial distributions of NSs/MSPs and NS/MSP binaries due to three-body interactions, making these clusters more likely to host NSs in the cluster outskirts. △ Less

Submitted 22 September, 2023; originally announced September 2023.

Comments: 13 pages, 7 figures, 2 tables, submitted to MNRAS

$clustertools$: A Python Package for Analyzing Star Cluster Simulations

Authors: Jeremy J. Webb

Abstract: $clustertools$ is a Python package for analyzing star cluster simulations. The package is built around the $StarCluster$ class, which stores all data read in from the snapshot of a given model star cluster. The package contains functions for loading data from commonly used $N$-body codes, generic snapshots, and software for generating initial conditions. All operations and functions within $cluste… ▽ More $clustertools$ is a Python package for analyzing star cluster simulations. The package is built around the $StarCluster$ class, which stores all data read in from the snapshot of a given model star cluster. The package contains functions for loading data from commonly used $N$-body codes, generic snapshots, and software for generating initial conditions. All operations and functions within $clustertools$ are then designed to act on a $StarCluster$. $clustertools$ can be used for unit and coordinate transformations, the calculation of key structural and kinematic parameters, analysis of the cluster's orbit and tidal tails, and measuring common cluster properties like its mass function, density profile, and velocity dispersion profile (among others). While originally designed with star clusters in mind, $clustertools$ can be used to study other types of $N$-body systems, including stellar streams and dark matter sub-halos. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Comments: Accepted for Publication in the Journal of Open Source Software, also see https://clustertools.readthedocs.io/en/latest/ for complete documentation

Journal ref: Journal of Open Source Software, 8(85), 4483

Made-to-Measure Modelling of Globular Clusters

Authors: Jeremy J. Webb, Jason A. S. Hunt, Jo Bovy

Abstract: We present the first application of the made-to-measure method for modelling dynamical systems to globular clusters. Through the made-to-measure algorithm, the masses of individual particles within a model cluster are adjusted while the system evolves forward in time via a gravitational $N$-body code until the model cluster is able to reproduce select properties of an observed cluster. The method… ▽ More We present the first application of the made-to-measure method for modelling dynamical systems to globular clusters. Through the made-to-measure algorithm, the masses of individual particles within a model cluster are adjusted while the system evolves forward in time via a gravitational $N$-body code until the model cluster is able to reproduce select properties of an observed cluster. The method is first applied to observations of mock isotropic and anisotropic clusters while fitting against the cluster's three dimensional or projected density profile, density weighted mean-squared velocity profile, or its density profile with individual mean-squared velocity profiles. We find that a cluster's three-dimensional density profile can easily be reproduced by the made-to-measure method, with minor discrepancies in the outer regions if fitting against a cluster's projected surface density or projected kinematic properties. If an observed cluster is anisotropic, only fitting against the cluster's density profile and individual mean-squared velocity profiles will fully recover the full degree of anisotropy. Partial anisotropy can be recovered as long as two kinematic properties are included in the fit. We further apply the method to observations of the Galactic globular cluster M4 and generate a complete six-dimensional representation of the cluster that reproduces observations of its surface density profile, mean-squared proper motion velocity profile, and mean-squared line of sight velocity profile. The M2M method predicts M4 is primarily isotropic with a mass of $9.2 \pm 0.4 \times 10^4\, M_{\odot}$ and a half-mass radius of $3.7 \pm 0.1$ pc. △ Less

Submitted 13 December, 2022; originally announced December 2022.

Comments: 11 pages, 10 figures, submitted to MNRAS for publication

On the Fast Track: Rapid construction of stellar stream paths

Authors: Nathaniel Starkman, Jo Bovy, Jeremy J. Webb, Daniela Calvetti, Erkki Somersalo

Abstract: Stellar streams are sensitive probes of the Galactic potential. The likelihood of a stream model given stream data is often assessed using simulations. However, comparing to simulations is challenging when even the stream paths can be hard to quantify. Here we present a novel application of Self-Organizing Maps and first-order Kalman Filters to reconstruct a stream's path, propagating measurement… ▽ More Stellar streams are sensitive probes of the Galactic potential. The likelihood of a stream model given stream data is often assessed using simulations. However, comparing to simulations is challenging when even the stream paths can be hard to quantify. Here we present a novel application of Self-Organizing Maps and first-order Kalman Filters to reconstruct a stream's path, propagating measurement errors and data sparsity into the stream path uncertainty. The technique is Galactic-model independent, non-parametric, and works on phase-wrapped streams. With this technique, we can uniformly analyze and compare data with simulations, enabling both comparison of simulation techniques and ensemble analysis with stream tracks of many stellar streams. Our method is implemented in the public Python package TrackStream, available at https://github.com/nstarman/trackstream. △ Less

Submitted 1 December, 2022; originally announced December 2022.

Comments: 16 pages, 11 figures, preprint

Searching for the extra-tidal stars of globular clusters using high-dimensional analysis and a core particle spray code

Authors: Steffani M. Grondin, Jeremy J. Webb, Nathan W. C. Leigh, Joshua S. Speagle, Reem J. Khalifeh

Abstract: Three-body interactions can eject stars from the core of a globular cluster, causing them to enter the Galactic halo as extra-tidal stars. While finding extra-tidal stars is imperative for understanding cluster evolution, connecting isolated extra-tidal field stars back to their birth cluster is extremely difficult. In this work, we present a new methodology consisting of high-dimensional data ana… ▽ More Three-body interactions can eject stars from the core of a globular cluster, causing them to enter the Galactic halo as extra-tidal stars. While finding extra-tidal stars is imperative for understanding cluster evolution, connecting isolated extra-tidal field stars back to their birth cluster is extremely difficult. In this work, we present a new methodology consisting of high-dimensional data analysis and a particle spray code to identify extra-tidal stars of any Galactic globular cluster using M3 as a case study. Using the t-Stochastic Neighbour Embedding (t-SNE) and Uniform Manifold Approximation and Projection (UMAP) machine learning dimensionality reduction algorithms, we first identify a set of 103 extra-tidal candidates in the APOGEE DR17 data catalogue with chemical abundances similar to M3 stars. To confirm each candidate's extra-tidal nature, we introduce corespray - a new Python-based three-body particle spray code that simulates extra-tidal stars for any Galactic globular cluster. Using Gaia EDR3 proper motions and APOGEE DR17 radial velocities, we apply multivariate Gaussian modelling and an extreme deconvolution to identify the extra-tidal candidates that are more likely to be associated with a distribution of corespray-simulated M3 extra-tidal stars than the field. Through these methods, we identify 10 new high-probability extra-tidal stars produced via three-body interactions in M3. We also explore whether any of our extra-tidal candidates are consistent with being ejected from M3 through different dynamical processes. Future applications of corespray will yield better understandings of core dynamics, star formation histories and binary fractions in globular clusters. △ Less

Submitted 16 November, 2022; v1 submitted 22 July, 2022; originally announced July 2022.

Comments: 16 pages, 8 figures, 3 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

The thermodynamics of stellar multiplicity: dynamical evolution of binary star populations in dense stellar environments

Authors: N. W. C. Leigh, N. C. Stone, J. J. Webb, W. Lyra

Abstract: We recently derived, using the density-of-states approximation, analytic distribution functions for the outcomes of direct single-binary scatterings (Stone & Leigh 2019). Using these outcome distribution functions, we present in this paper a self-consistent statistical mechanics-based analytic model obtained using the Fokker-Planck limit of the Boltzmann equation. Our model quantifies the dominant… ▽ More We recently derived, using the density-of-states approximation, analytic distribution functions for the outcomes of direct single-binary scatterings (Stone & Leigh 2019). Using these outcome distribution functions, we present in this paper a self-consistent statistical mechanics-based analytic model obtained using the Fokker-Planck limit of the Boltzmann equation. Our model quantifies the dominant gravitational physics, combining both strong and weak single-binary interactions, that drives the time evolution of binary orbital parameter distributions in dense stellar environments. We focus in particular the distributions of binary orbital energies and eccentricities. We find a novel steady state distribution of binary eccentricities, featuring strong depletions of both the highest and the lowest eccentricity binaries. In energy space, we compare the predictions of our analytic model to the results of numerical N-body simulations, and find that the agreement is good for the initial conditions considered here. This work is a first step toward the development of a fully self-consistent semi-analytic model for dynamically evolving binary star populations in dense stellar environments due to direct few-body interactions. △ Less

Submitted 30 May, 2022; originally announced May 2022.

Comments: 12 pages, 8 figures, 2 appendices

Bayesian Inference of Globular Cluster Properties Using Distribution Functions

Authors: Gwendolyn M. Eadie, Jeremy J. Webb, Jeffrey S. Rosenthal

Abstract: We present a Bayesian inference approach to estimating the cumulative mass profile and mean squared velocity profile of a globular cluster given the spatial and kinematic information of its stars. Mock globular clusters with a range of sizes and concentrations are generated from lowered isothermal dynamical models, from which we test the reliability of the Bayesian method to estimate model paramet… ▽ More We present a Bayesian inference approach to estimating the cumulative mass profile and mean squared velocity profile of a globular cluster given the spatial and kinematic information of its stars. Mock globular clusters with a range of sizes and concentrations are generated from lowered isothermal dynamical models, from which we test the reliability of the Bayesian method to estimate model parameters through repeated statistical simulation. We find that given unbiased star samples, we are able to reconstruct the cluster parameters used to generate the mock cluster and the cluster's cumulative mass and mean velocity squared profiles with good accuracy. We further explore how strongly biased sampling, which could be the result of observing constraints, may affect this approach. Our tests indicate that if we instead have biased samples, then our estimates can be off in certain ways that are dependent on cluster morphology. Overall, our findings motivate obtaining samples of stars that are as unbiased as possible. This may be achieved by combining information from multiple telescopes (e.g., Hubble and Gaia), but will require careful modeling of the measurement uncertainties through a hierarchical model, which we plan to pursue in future work. △ Less

Submitted 30 August, 2021; originally announced August 2021.

Comments: submitted to ApJ; 21 pages, 11 figures

Variation in the Stellar Mass Function Along Stellar Streams

Authors: Jeremy J. Webb, Jo Bovy

Abstract: Stellar streams are the inevitable end product of star cluster evolution, with the properties of a given stream being related to its progenitor. We consider how the dynamical history of a progenitor cluster, as traced by the evolution of its stellar mass function, is reflected in the resultant stream. We generate model streams by evolving star clusters with a range of initial half-mass relaxation… ▽ More Stellar streams are the inevitable end product of star cluster evolution, with the properties of a given stream being related to its progenitor. We consider how the dynamical history of a progenitor cluster, as traced by the evolution of its stellar mass function, is reflected in the resultant stream. We generate model streams by evolving star clusters with a range of initial half-mass relaxation times and dissolution times via direct N-body simulations. Stellar streams that dissolve quickly show no variation in the stellar mass function along the stream. Variation is, however, observed along streams with progenitor clusters that dissolve after several relaxation times. The mass function at the edges of a stream is approximately primordial as it is populated by the first stars to escape the cluster before segregation occurs. Moving inwards the mass function steepens as the intermediate parts of the stream consist of mostly low-mass stars that escaped the cluster after some segregation has occurred. The centre of the stream is then marked by a flatter mass function, as the region is dominated by high-mass stars that quickly segregated to the progenitor cluster's centre and were the last stars to become unbound. We further find that the maximum slope of the mass function along the stream and the rate at which it decreases with distance from the dissolved progenitor serve as proxies for the dynamical state reached by the progenitor cluster before dissolution; this may be able to be applied to observed streams with near-future observations. △ Less

Submitted 24 November, 2021; v1 submitted 4 August, 2021; originally announced August 2021.

Comments: 11 pages, 11 figures, accepted by MNRAS for publication

The Effects of $Λ$CDM Dark Matter Substructure on the Orbital Evolution of Star Clusters

Authors: Nicholas Pavanel, Jeremy J. Webb

Abstract: We present a comprehensive study on how perturbations due to a distribution of $Λ$CDM dark matter subhalos can lead to star clusters deviating from their orbits. Through a large suite of massless test particle simulations, we find that (1) subhalos with masses less than $10^8 M_{\odot}$ negligibly affect test particle orbits, (2) perturbations lead to orbital deviations only in environments with s… ▽ More We present a comprehensive study on how perturbations due to a distribution of $Λ$CDM dark matter subhalos can lead to star clusters deviating from their orbits. Through a large suite of massless test particle simulations, we find that (1) subhalos with masses less than $10^8 M_{\odot}$ negligibly affect test particle orbits, (2) perturbations lead to orbital deviations only in environments with substructure fractions $f_{sub} \geq 1\%$, (3) perturbations from denser subhalos produce larger orbital deviations, and (4) subhalo perturbations that are strong relative to the background tidal field lead to larger orbital deviations. To predict how the variation in test particle orbital energy $σ_e(t)$ increases with time, we test the applicability of theory derived from single-mass subhalo populations to populations where subhalos have a mass spectrum. We find $σ_e(t)$ can be predicted for test particle evolution within a mass spectrum of subhalos by assuming subhalos all have masses equal to the mean subhalo mass and by using the local mean subhalo separation to estimate the change in test particle velocities due to subhalo interactions. Furthermore, the orbital distance variation at an orbital distance $r$ can be calculated via $σ_r=2.98 \times 10^{-5} \pm 8 \times 10^{-8} (\rm kpc^{-1} km^{-2} s^{2}) \times r \times σ_e$ with a dispersion about the line of best fit equalling 0.08 kpc. Finally, we conclude that clusters that orbit within 100 kpc of Milky Way-like galaxies experience a change no greater than $2\%$ in their dissolution times. △ Less

Submitted 22 February, 2021; v1 submitted 12 February, 2021; originally announced February 2021.

Comments: 9 pages, 7 figures, 1 table, accepted for publication in MNRAS

The Likelihood of Undiscovered Globular Clusters in the Outskirts of the Milky Way

Authors: Jeremy J. Webb, Raymond G. Carlberg

Abstract: The currently known Galactic globular cluster population extends out to a maximum galactocentric distance of $\sim$ 145 kpc, with the peculiarity that the outermost clusters predominantly have an inward velocity. Orbit averaging finds that this configuration occurs by chance about $6\%$ of the time, suggesting that several globular clusters with positive radial velocities remain undiscovered. We e… ▽ More The currently known Galactic globular cluster population extends out to a maximum galactocentric distance of $\sim$ 145 kpc, with the peculiarity that the outermost clusters predominantly have an inward velocity. Orbit averaging finds that this configuration occurs by chance about $6\%$ of the time, suggesting that several globular clusters with positive radial velocities remain undiscovered. We evaluate the expected number of undiscovered clusters at large distances under the assumption that the cluster population has a smooth radial distribution and is in equilibrium within the Milky Way's virial radius. By comparing the present day kinematic properties of outer clusters to random orbital configurations of the Galactic globular cluster system through orbit averaging, we estimate a likelihood of $73\%$ of there being at least one undiscovered globular cluster within the Milky Way. This estimate assumes the current population is complete out to 50 kpc, and increases to $91\%$ if the population is complete out to 150 kpc. The likelihood of there being two undiscovered clusters is between $60\%$ and $70\%$, with the likelihood of there being three undiscovered clusters being on the order of $50\%$. The most likely scenario is that the undiscovered clusters are moving outwards, which results in the outer cluster population being consistent with an equilibrium state. Searches for distant and possibly quite low concentration and very low metallicity globular clusters will be enabled with upcoming deep imaging surveys. △ Less

Submitted 15 December, 2020; originally announced December 2020.

Comments: 13 pages, 9 figures, 4 tables, resubmitted to MNRAS for publication

The Initial Properties of Young Star Clusters in M83

Authors: Jeremy J. Webb, Alison Sills

Abstract: The initial sizes and masses of massive star clusters provide information about the cluster formation process and also determine how cluster populations are modified and destroyed, which have implications for using clusters as tracers of galaxy assembly. Young massive cluster populations are often assumed to be unchanged since cluster formation, and therefore their distribution of masses and radii… ▽ More The initial sizes and masses of massive star clusters provide information about the cluster formation process and also determine how cluster populations are modified and destroyed, which have implications for using clusters as tracers of galaxy assembly. Young massive cluster populations are often assumed to be unchanged since cluster formation, and therefore their distribution of masses and radii are used as the initial values. However, the first few hundred million years of cluster evolution does change both cluster mass and cluster radius, through both internal and external processes. In this paper, we use a large suite of $N$-body cluster simulations in an appropriate tidal field to determine the best initial mass and initial size distributions of young clusters in the nearby galaxy M83. We find that the initial masses follow a power-law distribution with a slope of -2.7 $\pm$ 0.4, and the half-mass radii follow a log-normal distribution with a mean of 2.57 $\pm$ 0.04 pc and a dispersion of 1.59 $\pm$ 0.01 pc. The corresponding initial projected half-light radius function has a mean of 2.7 $\pm$ 0.3 pc and a dispersion of 1.7 $\pm$ 0.2 pc. The evolution of the initial mass and size distribution functions are consistent with mass loss and expansion due to stellar evolution, independent of the external tidal field and the cluster's initial density profile. Observed cluster sizes and masses should not be used as the initial values, even when clusters are only a few hundred million years old. △ Less

Submitted 14 December, 2020; originally announced December 2020.

Comments: 8 pages, 6 figures, accepted for publication in MNRAS

Radial variation of the stellar mass functions in the globular clusters M15 and M30: clues of a non-standard IMF?

Authors: Mario Cadelano, Emanuele Dalessandro, Jeremy J. Webb, Enrico Vesperini, Daniele Lattanzio, Giacomo Beccari, Matias Gomez, Lorenzo Monaco

Abstract: We exploit a combination of high-resolution Hubble Space Telescope and wide-field ESO-VLT observations to study the slope of the global mass function (alphaG) and its radial variation (alpha(r)) in the two dense, massive and post core-collapse globular clusters M15 and M30. The available data-set samples the clusters' Main Sequence down to 0.2 Msun and the photometric completeness allows the study… ▽ More We exploit a combination of high-resolution Hubble Space Telescope and wide-field ESO-VLT observations to study the slope of the global mass function (alphaG) and its radial variation (alpha(r)) in the two dense, massive and post core-collapse globular clusters M15 and M30. The available data-set samples the clusters' Main Sequence down to 0.2 Msun and the photometric completeness allows the study of the mass function between 0.40 Msun and 0.75 Msun from the central regions out to their tidal radii. We find that both clusters show a very similar variation in alpha(r) as a function of clustercentric distance. They both exhibit a very steep variation in alpha(r) in the central regions, which then attains almost constant values in the outskirts. Such a behavior can be interpreted as the result of long-term dynamical evolution of the systems driven by mass-segregation and mass-loss processes. We compare these results with a set of direct N-body simulations and find that they are only able to reproduce the observed values of alpha(r) and alphaG at dynamical ages (t/trh) significantly larger than those derived from the observed properties of both clusters. We investigate possible physical mechanisms responsible for such a discrepancy and argue that both clusters might be born with a non-standard (flatter/bottom-lighter) initial mass function. △ Less

Submitted 4 September, 2020; originally announced September 2020.

Comments: Accepted for publication by MNRAS

The Effects of Dwarf Galaxies on the Orbital Evolution of Galactic Globular Clusters

Authors: Turner Garrow, Jeremy J. Webb, Jo Bovy

Abstract: We investigate the effect that dwarf galaxies have on the orbits, tidal histories, and assumed formation environment of Milky Way globular clusters. We determine the orbits of the Milky Way's 150 globular clusters in a gravitational potential both with and without dwarf galaxies. We find that the presence of a small number of satellite galaxies can affect the orbits of many of the globular cluster… ▽ More We investigate the effect that dwarf galaxies have on the orbits, tidal histories, and assumed formation environment of Milky Way globular clusters. We determine the orbits of the Milky Way's 150 globular clusters in a gravitational potential both with and without dwarf galaxies. We find that the presence of a small number of satellite galaxies can affect the orbits of many of the globular clusters. Over 12 Gyr, we find that the semi-major axis and orbital eccentricity of individual clusters fluctuate with dispersions on the order of $\sim 10\%$ and $\sim 4\%$, respectively. Outer clusters are more strongly affected by dwarf galaxies than inner clusters, with their semi-major axis and orbital eccentricities fluctuating by more than $\sim 15\%$ and $\sim 5\%$, respectively. Using detailed $N$-body simulations of select clusters, we find that altering their orbital histories can lead to different mass loss rates and structural evolution. Furthermore, we caution against using kinematics alone to identify whether a Galactic cluster formed in-situ or was accreted during a past merger event as these values are no longer conserved. The presence of dwarf galaxies causes the orbital energies and actions of individual clusters to evolve over time, spanning a wider range than that coming from random uncertainties in a cluster's proper motions and radial velocity. △ Less

Submitted 20 December, 2020; v1 submitted 27 July, 2020; originally announced July 2020.

Comments: 10 pages, 6 figures, 2 tables, Accepted to MNRAS

And In The Darkness Unbind Them: High-Resolution Simulations of Dark Matter Subhalo Disruption in a Milky Way-like Tidal Field

Authors: Jeremy J. Webb, Jo Bovy

Abstract: We compare the results of high-resolution simulations of individual dark matter subhalos evolving in external tidal fields with and without baryonic bulge and disk components, where the average dark matter particle mass is three orders of magnitude smaller than cosmological zoom-in simulations of galaxy formation. The Via Lactea II simulation is used to setup our initial conditions and provides a… ▽ More We compare the results of high-resolution simulations of individual dark matter subhalos evolving in external tidal fields with and without baryonic bulge and disk components, where the average dark matter particle mass is three orders of magnitude smaller than cosmological zoom-in simulations of galaxy formation. The Via Lactea II simulation is used to setup our initial conditions and provides a basis for our simulations of subhalos in a dark matter-only tidal field, while an observationally motivated model for the Milky Way is used for the tidal field that is comprised of a dark matter halo, a stellar disk, and a stellar bulge. Our simulations indicate that including stellar components in the tidal field results in the number of subhalos in Milky Way-like galaxies being only $65\%$ of what is predicted by $Λ$ Cold Dark Matter ($Λ$CDM). For subhalos with small pericentres $(r_p \lesssim 25$ kpc), the subhalo abundance is reduced further to $40\%$, with the surviving subhalos being less dense and having a tangentially-anisotropic orbital distribution. Conversely, subhalos with larger pericentres are minimally affected by the inclusion of a stellar component in the tidal field, with the total number of outer subhalos $\approx 75\%$ of the $Λ$CDM prediction. The densities of outer subhalos are comparable to predictions from $Λ$CDM, with the subhalos having an isotropic distribution of orbits. These ratios are higher than those found in previous studies that include the effects baryonic matter, which are affected by spurious disruption caused by low resolution. △ Less

Submitted 16 December, 2020; v1 submitted 11 June, 2020; originally announced June 2020.

Comments: 14 pages, 14 figures, 1 table, Accepted for publication in MNRAS

Mass loss from massive globular clusters in tidal fields

Authors: Yohai Meiron, Jeremy J. Webb, Jongsuk Hong, Peter Berczik, Rainer Spurzem, Raymond G. Carlberg

Abstract: Massive globular clusters lose stars via internal and external processes. Internal processes include mainly two-body relaxation, while external processes include interactions with the Galactic tidal field. We perform a suite of N-body simulations of such massive clusters using three different direct-summation N-body codes, exploring different Galactic orbits and particle numbers. By inspecting the… ▽ More Massive globular clusters lose stars via internal and external processes. Internal processes include mainly two-body relaxation, while external processes include interactions with the Galactic tidal field. We perform a suite of N-body simulations of such massive clusters using three different direct-summation N-body codes, exploring different Galactic orbits and particle numbers. By inspecting the rate at which a star's energy changes as it becomes energetically unbound from the cluster, we can neatly identify two populations we call kicks and sweeps, that escape through two-body encounters internal to the cluster and the external tidal field, respectively. We find that for a typical halo globular cluster on a moderately eccentric orbit, sweeps are far more common than kicks but the total mass loss rate is so low that these clusters can survive for tens of Hubble times. The different N-body codes give largely consistent results, but we find that numerical artifacts may arise in relation to the time step parameter of the Hermite integration scheme, namely that the value required for convergent results is sensitive to the number of particles. △ Less

Submitted 24 March, 2021; v1 submitted 2 June, 2020; originally announced June 2020.

Strong chemical tagging with APOGEE: 21 candidate star clusters that have dissolved across the Milky Way disc

Authors: Natalie Price-Jones, Jo Bovy, Jeremy J. Webb, Carlos Allende Prieto, Rachael Beaton, Joel R. Brownstein, Roger E. Cohen, Katia Cunha, John Donor, Peter M. Frinchaboy, D. A. García-Hernández, Richard R. Lane, Steven R. Majewski, David L. Nidever, Alexandre Roman-Lopes

Abstract: Chemically tagging groups of stars born in the same birth cluster is a major goal of spectroscopic surveys. To investigate the feasibility of such strong chemical tagging, we perform a blind chemical tagging experiment on abundances measured from APOGEE survey spectra. We apply a density-based clustering algorithm to the eight dimensional chemical space defined by [Mg/Fe], [Al/Fe], [Si/Fe], [K/Fe]… ▽ More Chemically tagging groups of stars born in the same birth cluster is a major goal of spectroscopic surveys. To investigate the feasibility of such strong chemical tagging, we perform a blind chemical tagging experiment on abundances measured from APOGEE survey spectra. We apply a density-based clustering algorithm to the eight dimensional chemical space defined by [Mg/Fe], [Al/Fe], [Si/Fe], [K/Fe], [Ti/Fe], [Mn/Fe], [Fe/H], and [Ni/Fe], abundances ratios which together span multiple nucleosynthetic channels. In a high quality sample of 182,538 giant stars, we detect twenty-one candidate clusters with more than fifteen members. Our candidate clusters are more chemically homogeneous than a population of non-member stars with similar [Mg/Fe] and [Fe/H], even in abundances not used for tagging. Group members are consistent with having the same age and fall along a single stellar-population track in logg vs. Teff space. Each group's members are distributed over multiple kpc, and the spread in their radial and azimuthal actions increases with age. We qualitatively reproduce this increase using N-body simulations of cluster dissolution in Galactic potentials that include transient winding spiral arms. Observing our candidate birth clusters with high-resolution spectroscopy in other wavebands to investigate their chemical homogeneity in other nucleosynthetic groups will be essential to confirming the efficacy of strong chemical tagging. Our initially spatially-compact but now widely dispersed candidate clusters will provide novel limits on chemical evolution and orbital diffusion in the Galactic disc, and constraints on star formation in loosely-bound groups. △ Less

Submitted 26 June, 2020; v1 submitted 8 April, 2020; originally announced April 2020.

Comments: 15 pages, 9 figures, accepted by MNRAS

Searching for Solar Siblings in APOGEE and $Gaia$ DR2 with N-body Simulations

Authors: Jeremy J. Webb, Natalie Price-Jones, Jo Bovy, Simon Portegies Zwart, Jason A. S. Hunt, J. Ted Mackereth, Henry W. Leung

Abstract: We make use of APOGEE and $Gaia$ data to identify stars that are consistent with being born in the same association or star cluster as the Sun. We limit our analysis to stars that match solar abundances within their uncertainties, as they could have formed from the same Giant Molecular Cloud (GMC) as the Sun. We constrain the range of orbital actions that solar siblings can have with a suite of si… ▽ More We make use of APOGEE and $Gaia$ data to identify stars that are consistent with being born in the same association or star cluster as the Sun. We limit our analysis to stars that match solar abundances within their uncertainties, as they could have formed from the same Giant Molecular Cloud (GMC) as the Sun. We constrain the range of orbital actions that solar siblings can have with a suite of simulations of solar birth clusters evolved in static and time-dependent tidal fields. The static components of each galaxy model are the bulge, disk, and halo, while the various time-dependent components include a bar, spiral arms, and GMCs. In galaxy models without GMCs, simulated solar siblings all have $J_R < 122$ km $\rm s^{-1}$ kpc, $990 < L_z < 1986$ km $\rm s^{-1}$ kpc, and $0.15 < J_z < 0.58$ km $\rm s^{-1}$ kpc. Given the actions of stars in APOGEE and $Gaia$, we find 104 stars that fall within this range. One candidate in particular, Solar Sibling 1, has both chemistry and actions similar enough to the solar values that strong interactions with the bar or spiral arms are not required for it to be dynamically associated with the Sun. Adding GMCs to the potential can eject solar siblings out of the plane of the disk and increase their $J_z$, resulting in a final candidate list of 296 stars. The entire suite of simulations indicate that solar siblings should have $J_R < 122$ km $\rm s^{-1}$ kpc, $353 < L_z < 2110$ km $\rm s^{-1}$ kpc, and $J_z < 0.8$ km $\rm s^{-1}$ kpc. Given these criteria, it is most likely that the association or cluster that the Sun was born in has reached dissolution and is not the commonly cited open cluster M67. △ Less

Submitted 23 March, 2020; v1 submitted 3 October, 2019; originally announced October 2019.

Comments: 12 pages, 9 figures, 3 tables, Accepted for publication in MNRAS

Modelling the Effects of Dark Matter Substructure on Globular Cluster Evolution with the Tidal Approximation

Authors: Jeremy J. Webb, Jo Bovy, Raymond G. Carlberg, Mark Gieles

Abstract: We present direct $N$-body simulations of tidally filling 30,000 ${\rm M}_\odot$ star clusters orbiting between 10 kpc and 100 kpc in galaxies with a range of dark matter substructure properties. The time-dependent tidal force is determined based on the combined tidal tensor of the galaxy's smooth and clumpy dark matter components, the latter of which causes fluctuations in the tidal field that ca… ▽ More We present direct $N$-body simulations of tidally filling 30,000 ${\rm M}_\odot$ star clusters orbiting between 10 kpc and 100 kpc in galaxies with a range of dark matter substructure properties. The time-dependent tidal force is determined based on the combined tidal tensor of the galaxy's smooth and clumpy dark matter components, the latter of which causes fluctuations in the tidal field that can heat clusters. The strength and duration of these fluctuations are sensitive to the local dark matter density, substructure fraction, sub-halo mass function, and the sub-halo mass-size relation. Based on the cold dark matter framework, we initially assume sub-halos are Hernquist spheres following a power-law mass function between $10^5$ and $10^{11} {\rm M}_\odot$ and find that tidal fluctuations are too weak and too short to affect star cluster evolution. Treating sub-halos as point masses, to explore how denser sub-halos affect clusters, we find that only sub-halos with masses greater than $10^{6} {\rm M}_\odot$ will cause cluster dissolution times to decrease. These interactions can also decrease the size of a cluster while increasing the velocity dispersion and tangential anisotropy in the outer regions via tidal heating. Hence increased fluctuations in the tidal tensor, especially fluctuations that are due to low-mass halos, do not necessarily translate into mass loss. We further conclude that the tidal approximation can be used to model cluster evolution in the tidal fields of cosmological simulations with a minimum cold dark matter sub-halo mass of $10^{6} {\rm M}_\odot$, as the effect of lower-mass sub-halos on star clusters is negligible. △ Less

Submitted 30 July, 2019; originally announced July 2019.

Comments: 17 pages, 10 figures, Accepted for publication in MNRAS

The evolution of kicked stellar-mass black holes in star cluster environments II. Rotating star clusters

Authors: Jeremy J. Webb, Nathan W. C. Leigh, Roberto Serrano, Jillian Bellovary, K. E. Saavik Ford, Barry McKernan, Mario Spera, Alessandro A. Trani

Abstract: In this paper, we continue our study on the evolution of black holes (BHs) that receive velocity kicks at the origin of their host star cluster potential. We now focus on BHs in rotating clusters that receive a range of kick velocities in different directions with respect to the rotation axis. We perform N-body simulations to calculate the trajectories of the kicked BHs and develop an analytic fra… ▽ More In this paper, we continue our study on the evolution of black holes (BHs) that receive velocity kicks at the origin of their host star cluster potential. We now focus on BHs in rotating clusters that receive a range of kick velocities in different directions with respect to the rotation axis. We perform N-body simulations to calculate the trajectories of the kicked BHs and develop an analytic framework to study their motion as a function of the host cluster and the kick itself. Our simulations indicate that for a BH that is kicked outside of the cluster's core, as its orbit decays in a rotating cluster the BH will quickly gain angular momentum as it interacts with stars with high rotational frequencies. Once the BH decays to the point where its orbital frequency equals that of local stars, its orbit will be circular and dynamical friction becomes ineffective since local stars will have low relative velocities. After circularization, the BH's orbit decays on a longer timescale than if the host cluster was not rotating. Hence BHs in rotating clusters will have longer orbital decay times. The timescale for orbit circularization depends strongly on the cluster's rotation rate and the initial kick velocity, with kicked BHs in slowly rotating clusters being able to decay into the core before circularization occurs. The implication of the circularization phase is that the probability of a BH undergoing a tidal capture event increases, possibly aiding in the formation of binaries and high-mass BHs. △ Less

Submitted 9 July, 2019; originally announced July 2019.

Comments: 13 pages, 11 figures, Accepted for publication in MNRAS

The orbital anisotropy profiles of nearby globular clusters from Gaia Data Release 2

Authors: Abhinav Jindal, Jeremy J. Webb, Jo Bovy

Abstract: Gaia Data Release 2 provides a wealth of data to study the internal structure of nearby globular clusters. We use this data to investigate the internal kinematics of 11 nearby globular clusters, with a particular focus on their poorly-studied outer regions. We apply a strict set of selection criteria to remove contaminating sources and create pure cluster-member samples over a significant fraction… ▽ More Gaia Data Release 2 provides a wealth of data to study the internal structure of nearby globular clusters. We use this data to investigate the internal kinematics of 11 nearby globular clusters, with a particular focus on their poorly-studied outer regions. We apply a strict set of selection criteria to remove contaminating sources and create pure cluster-member samples over a significant fraction of the radial range of each cluster. We confirm previous measurements of rotation (or a lack thereof) in the inner regions of several clusters, while extending the detection of rotation well beyond where it was previously measured and finding a steady decrease in rotation with radius. We also determine the orbital anisotropy profile and determine that clusters have isotropic cores, are radially anisotropic out to $\approx$ 4 half-light radii or $35\%$ of their limiting radii, and are then isotropic out to the limits of our datasets. We detect for the first time the presence of radial anisotropy in M 22, while confirming previous detections of radial anisotropy in 47 Tuc, M 3, M 13, M 15, and $ω$ Cen's innermost regions. The implications of these measurements are that clusters can be separated into two categories: 1) clusters with observed radial anisotropy that likely formed tidally under-filling or are dynamically young, and 2) clusters that are primarily isotropic that likely formed tidally filling or are dynamically old. △ Less

Submitted 26 March, 2019; originally announced March 2019.

Comments: 10 pages, 5 figures, Submitted to MNRAS

Spatial mixing of binary stars in multiple-population globular clusters

Authors: Jongsuk Hong, Saahil Patel, Enrico Vesperini, Jeremy J. Webb, Emanuele Dalessandro

Abstract: We present the results of a study aimed at investigating the effects of dynamical evolution on the spatial distribution and mixing of primordial binary stars in multiple-population globular clusters. Multiple stellar population formation models predict that second-generation (SG) stars form segregated in the inner regions of a more extended first-generation (FG) cluster. Our study, based on the… ▽ More We present the results of a study aimed at investigating the effects of dynamical evolution on the spatial distribution and mixing of primordial binary stars in multiple-population globular clusters. Multiple stellar population formation models predict that second-generation (SG) stars form segregated in the inner regions of a more extended first-generation (FG) cluster. Our study, based on the results of a survey of N-body simulations, shows that the spatial mixing process for binary stars is more complex than that of single stars since additional processes such as binary ionization, recoil and ejection following binary-single and binary-binary interactions play a key role in determining the spatial distribution of the population of surviving binaries. The efficiency and relative importance of these additional effects depends on the binary binding energy and determines the timescale of the spatial mixing of FG and SG binaries. Our simulations illustrate the role of ionization, recoil and ejection combined with the effects of mass segregation driven by two-body relaxation and show that the complex interplay of all these processes results in a significant extension of the time needed for the complete spatial mixing of FG and SG binaries compared to that of single stars. Clusters in which FG and SG single stars have already reached complete spatial mixing might be characterized by a significant radial gradient in the ratio of the FG-to-SG binary fraction. The implications of the delayed mixing of FG and SG binaries for the differences between the kinematics of the two populations are discussed. △ Less

Submitted 4 December, 2018; originally announced December 2018.

Comments: 8 pages, 7 figures, 1 table; accepted for publication in MNRAS

A systematic analysis of star cluster disruption by tidal shocks -- I. Controlled N-body simulations and a new theoretical model

Authors: Jeremy J. Webb, Marta Reina-Campos, J. M. Diederik Kruijssen

Abstract: Understanding the evolution of stellar clusters in an evolving tidal field is critical for studying the disruption of stellar clusters in a cosmological context. We systematically characterise the response of stellar clusters to tidal shocks using controlled N-body simulations of clusters with various properties that are subjected to different types of shocks. We find that the strength of the shoc… ▽ More Understanding the evolution of stellar clusters in an evolving tidal field is critical for studying the disruption of stellar clusters in a cosmological context. We systematically characterise the response of stellar clusters to tidal shocks using controlled N-body simulations of clusters with various properties that are subjected to different types of shocks. We find that the strength of the shock and the density of the cluster within the half-mass radius are the dominant properties that drive the amount of mass lost by the cluster, with the shape of the cluster profile being of minor influence. When the shock is applied as two separate sub-shocks, the amount of mass loss during the second sub-shock is sensitive to the gap time between them. Clusters that experience successive sub-shocks separated by less than their crossing time attain the same masses and sizes at the end of the simulation. However, clusters subjected to sub-shocks separated by more than a crossing time experience different evolutionary histories. The amount of mass lost in the N-body models and its scaling with shock and cluster properties differs from that predicted by classical tidal disruption theory. We demonstrate that the discrepancy is alleviated by including a dependence on the escape time-scale of unbound stars, analogously to mass loss driven by two-body relaxation. With our new theoretical model for shock-driven mass loss, the predicted relative amounts of mass loss agree with the results of the N-body simulations to ~0.3 dex across the full suite of simulations. △ Less

Submitted 7 May, 2019; v1 submitted 30 November, 2018; originally announced December 2018.

Comments: 14 pages, 12 figures, 1 table, Accepted for publication in MNRAS

Searching for the GD-1 Stream Progenitor in Gaia DR2 with Direct N-body Simulations

Authors: Jeremy J. Webb, Jo Bovy

Abstract: We perform a large suite of direct N-body simulations aimed at revealing the location of the progenitor, or its remnant, of the GD-1 stream. Data from \gaia\ DR2 reveals the GD-1 stream extends over $\approx 100^\circ$, allowing us to determine the stream's leading and trailing ends. Our models suggest the length of the stream is consistent with a dynamical age of between 2-3 Gyr and the exact len… ▽ More We perform a large suite of direct N-body simulations aimed at revealing the location of the progenitor, or its remnant, of the GD-1 stream. Data from \gaia\ DR2 reveals the GD-1 stream extends over $\approx 100^\circ$, allowing us to determine the stream's leading and trailing ends. Our models suggest the length of the stream is consistent with a dynamical age of between 2-3 Gyr and the exact length, width and location of the GD-1 stream correspond to the stream's progenitor being located between $-30^\circ < φ_{1,\mathrm{pro}} < -45^{\circ}$ in the standard GD-1 coordinate system. The model stream density profiles reveal that intact progenitors leave a strong over-density, recently-dissolved progenitors appear as gaps in the stream as escaped stars continue to move away from the remnant progenitor's location, and long-dissolved progenitors leave no observational signature on the remaining stream. Comparing our models to the GD-1 stream yields two possible scenarios for its progenitor's history: a) the progenitor reached dissolution approximately 500 Myr ago during the cluster's previous perigalactic pass and is both located at and responsible for the observed gap at $φ_1=-40^{\circ}$ or b) the progenitor reached dissolution over 2.5 Gyr ago, the fully-dissolved remnant is at $-30^\circ < φ_1 < -45^{\circ}$, and an observational signature of its location no longer exists. That the dissolved progenitor is in the range $-30^\circ < φ_1 < -45^{\circ}$ implies that density fluctuations outside of this range, e.g., a deep gap at $φ_1 \approx -20^\circ$, are likely produced by compact baryonic or dark-matter perturbers. △ Less

Submitted 21 March, 2019; v1 submitted 16 November, 2018; originally announced November 2018.

Comments: 9 pages, 5 figures, Accepted for publication in MNRAS

Rediscovering the Tidal Tails of NGC 288 with Gaia DR2

Authors: Shaziana Kaderali, Jason A. S. Hunt, Jeremy J. Webb, Natalie Price-Jones, Raymond Carlberg

Abstract: NGC 288 is a Galactic globular cluster having observed extra tidal structure, without confirmed tidal tails. $Gaia$ DR2 provides photometric and astrometric data for many of the stars in NGC 288 and its extra tidal structure. To compare with $Gaia$ data, we simulate an $N$-body model of a star cluster with the same orbit as NGC 288 in a Milky Way potential. The simulation shows that the cluster fo… ▽ More NGC 288 is a Galactic globular cluster having observed extra tidal structure, without confirmed tidal tails. $Gaia$ DR2 provides photometric and astrometric data for many of the stars in NGC 288 and its extra tidal structure. To compare with $Gaia$ data, we simulate an $N$-body model of a star cluster with the same orbit as NGC 288 in a Milky Way potential. The simulation shows that the cluster forms tidal tails that are compressed along the cluster's orbit when it is at apocentre and are expected to be a diffuse bipolar structure. In this letter, we present a comparison between the simulation and observations from $Gaia$ DR2. We find that both the simulation and the observations share comparable trends in the position on the sky and proper motions of the extra-tidal stars, supporting the presence of tidal tails around NGC 288. △ Less

Submitted 1 February, 2019; v1 submitted 11 September, 2018; originally announced September 2018.

Comments: 5 pages, 5 figures

The Structural and Kinematic Evolution of Central Star Clusters in Dwarf Galaxies and Their Dependence on Dark Matter Halo Profiles

Authors: Jeremy J. Webb, Enrico Vesperini

Abstract: Through a suite of direct N-body simulations, we explore how the structural and kinematic evolution of a star cluster located at the center of a dwarf galaxy is affected by the shape of its host's dark matter density profile. The stronger central tidal fields of cuspier halos minimize the cluster's ability to expand in response to mass loss due to stellar evolution during its early evolutionary st… ▽ More Through a suite of direct N-body simulations, we explore how the structural and kinematic evolution of a star cluster located at the center of a dwarf galaxy is affected by the shape of its host's dark matter density profile. The stronger central tidal fields of cuspier halos minimize the cluster's ability to expand in response to mass loss due to stellar evolution during its early evolutionary stages and during its subsequent long-term evolution driven by two-body relaxation. Hence clusters evolving in cuspier dark matter halos are characterized by more compact sizes, higher velocity dispersions and remain approximately isotropic at all clustercentric distances. Conversely, clusters in cored halos can expand more and develop a velocity distribution profile that becomes increasingly radially anisotropic at larger clustercentric distances. Finally, the larger velocity dispersion of clusters evolving in cuspier dark matter profiles results in them having longer relaxation times. Hence clusters in cuspy galaxies relax at a slower rate and, consequently, they are both less mass segregated and farther from complete energy equipartition than cluster's in cored galaxies. Application of this work to observations allows for star clusters to be used as tools to measure the distribution of dark matter in dwarf galaxies and to distinguish isolated star clusters from ultra-faint dwarf galaxies. △ Less

Submitted 6 July, 2018; originally announced July 2018.

Comments: 8 pages, 7 figures, Accepted for publication in MNRAS

The Peculiar Radial Distribution of Multiple Populations in the massive globular cluster M80

Authors: E. Dalessandro, M. Cadelano, E. Vesperini, M. Salaris, F. R. Ferraro, B. Lanzoni, S. Raso, J. Hong, J. J. Webb, A. Zocchi

Abstract: We present a detailed analysis of the radial distribution of light-element multiple populations (LE-MPs) in the massive and dense globular cluster M80 based on the combination of UV and optical Hubble Space Telescope data. Surprisingly, we find that first generation stars (FG) are significantly more centrally concentrated than extreme second generation ones (SG) out to $\sim 2.5 r_h$ from the clus… ▽ More We present a detailed analysis of the radial distribution of light-element multiple populations (LE-MPs) in the massive and dense globular cluster M80 based on the combination of UV and optical Hubble Space Telescope data. Surprisingly, we find that first generation stars (FG) are significantly more centrally concentrated than extreme second generation ones (SG) out to $\sim 2.5 r_h$ from the cluster center. To understand the origin of such a peculiar behavior, we used a set of $N$-body simulations following the long-term dynamical evolution of LE-MPs. We find that, given the advanced dynamical state of the cluster, the observed difference does not depend on the primordial relative distributions of FG and SG stars. On the contrary, a difference of $\sim 0.05-0.10 M_{\odot}$ between the average masses of the two sub-populations is needed to account for the observed radial distributions. We argue that such a mass difference might be the result of the higher He abundance of SG stars (of the order of $ΔY\sim 0.05-0.06$) with respect to FG. Interestingly, we find that a similar He variation is necessary to reproduce the horizontal branch morphology of M80. These results demonstrate that differences in mass among LE-MPs, due to different He content, should be properly taken into account for a correct interpretation of their radial distribution, at least in dynamically evolved systems. △ Less

Submitted 9 April, 2018; originally announced April 2018.

Comments: 19 pages, 9 figures, 1 table; Accepted for publication in ApJ

Evolution of the Stellar Mass Function in Multiple-Population Globular Clusters

Authors: Enrico Vesperini, Jongsuk Hong, Jeremy J. Webb, Franca D'Antona, Annibale D'Ercole

Abstract: We present the results of a survey of N-body simulations aimed at studying the effects of the long-term dynamical evolution on the stellar mass function (MF) of multiple stellar populations in globular clusters. Our simulations show that if first-(1G) and second-generation (2G) stars have the same initial MF (IMF), the global MFs of the two populations are affected similarly by dynamical evolution… ▽ More We present the results of a survey of N-body simulations aimed at studying the effects of the long-term dynamical evolution on the stellar mass function (MF) of multiple stellar populations in globular clusters. Our simulations show that if first-(1G) and second-generation (2G) stars have the same initial MF (IMF), the global MFs of the two populations are affected similarly by dynamical evolution and no significant differences between the 1G and the 2G MFs arise during the cluster's evolution. If the two populations have different IMFs, dynamical effects do not completely erase memory of the initial differences. Should observations find differences between the global 1G and 2G MF, these would reveal the fingerprints of differences in their IMFs. Irrespective of whether the 1G and 2G populations have the same global IMF or not, dynamical effects can produce differences between the local (measured at various distances from the cluster centre) 1G and 2G MFs; these differences are a manifestation of the process of mass segregation in populations with different initial structural properties. In dynamically old and spatially mixed clusters, however, differences between the local 1G and 2G MFs can reveal differences between the 1G and 2G global MFs. In general, for clusters with any dynamical age, large differences between the local 1G and 2G MFs are more likely to be associated with differences in the global MF. Our study also reveals a dependence of the spatial mixing rate on the stellar mass, another dynamical consequence of the multiscale nature of multiple-population clusters. △ Less

Submitted 6 March, 2018; originally announced March 2018.

Comments: 12 pages, 11 figures. Accepted for publication in MNRAS

Kinematic fingerprint of core-collapsed globular clusters

Authors: P. Bianchini, J. J. Webb, A. Sills, E. Vesperini

Abstract: Dynamical evolution drives globular clusters toward core collapse, which strongly shapes their internal properties. Diagnostics of core collapse have so far been based on photometry only, namely on the study of the concentration of the density profiles. Here we present a new method to robustly identify core-collapsed clusters based on the study of their stellar kinematics. We introduce the \textit… ▽ More Dynamical evolution drives globular clusters toward core collapse, which strongly shapes their internal properties. Diagnostics of core collapse have so far been based on photometry only, namely on the study of the concentration of the density profiles. Here we present a new method to robustly identify core-collapsed clusters based on the study of their stellar kinematics. We introduce the \textit{kinematic concentration} parameter, $c_k$, the ratio between the global and local degree of energy equipartition reached by a cluster, and show through extensive direct $N$-body simulations that clusters approaching core collapse and in the post-core collapse phase are strictly characterized by $c_k>1$. The kinematic concentration provides a suitable diagnostic to identify core-collapsed clusters, independent from any other previous methods based on photometry. We also explore the effects of incomplete radial and stellar mass coverage on the calculation of $c_k$ and find that our method can be applied to state-of-art kinematic datasets. △ Less

Submitted 23 January, 2018; originally announced January 2018.

Comments: Accepted for publication in MNRAS Letter

The Evolution of Kicked Stellar-Mass Black Holes in Star Cluster Environments

Authors: Jeremy J. Webb, Nathan W. C. Leigh, Abhishek Singh, K. E. Saavik Ford, Barry McKernan, Jillian Bellovary

Abstract: We consider how dynamical friction acts on black holes that receive a velocity kick while located at the center of a gravitational potential, analogous to a star cluster, due to either a natal kick or the anisotropic emission of gravitational waves during a black hole-black hole merger. Our investigation specifically focuses on how well various Chandrasekhar-based dynamical friction models can pre… ▽ More We consider how dynamical friction acts on black holes that receive a velocity kick while located at the center of a gravitational potential, analogous to a star cluster, due to either a natal kick or the anisotropic emission of gravitational waves during a black hole-black hole merger. Our investigation specifically focuses on how well various Chandrasekhar-based dynamical friction models can predict the orbital decay of kicked black holes with $m_{bh} \lesssim 100 M_\odot$ due to an inhomogeneous background stellar field. In general, the orbital evolution of a kicked black hole follows that of a damped oscillator where two-body encounters and dynamical friction serve as sources of damping. However, we find models for approximating the effects of dynamical friction do not accurately predict the amount of energy lost by the black hole if the initial kick velocity $v_{k}$ is greater than the stellar velocity dispersion $σ$. For all kick velocities, we also find that two-body encounters with nearby stars can cause the energy evolution of a kicked BH to stray significantly from standard dynamical friction theory as encounters can sometimes lead to an energy gain. For larger kick velocities, we find the orbital decay of a black hole departs from classical theory completely as the black hole's orbital amplitude decays linearly with time as opposed to exponentially. Therefore, we have developed a linear decay formalism which scales linearly with black hole mass and $\frac{v_{k}}σ$ in order to account for the variations in the local gravitational potential. △ Less

Submitted 24 November, 2017; originally announced November 2017.

Comments: 13 pages, 9 figures, Accepted for publication in MNRAS

Modelling the Observed Stellar Mass Function and its Radial Variation in Galactic Globular Clusters

Authors: Jeremy J. Webb, Enrico Vesperini, Emanuele Dalessandro, Giacomo Beccari, Francesco R. Ferraro, Barbara Lanzoni

Abstract: We measure how the slope $α$ of the stellar mass function (MF) changes as a function of clustercentric distance $r$ in five Galactic globular clusters and compare $α(r)$ to predictions from direct $N$-body star cluster simulations. Theoretical studies predict that $α(r)$ (which traces the degree of mass segregation in a cluster) should steepen with time as a cluster undergoes two-body relaxation a… ▽ More We measure how the slope $α$ of the stellar mass function (MF) changes as a function of clustercentric distance $r$ in five Galactic globular clusters and compare $α(r)$ to predictions from direct $N$-body star cluster simulations. Theoretical studies predict that $α(r)$ (which traces the degree of mass segregation in a cluster) should steepen with time as a cluster undergoes two-body relaxation and that the amount by which the global MF can evolve from its initial state due to stellar escape is directly linked to $α(r)$. We find that the amount of mass segregation in M10, NGC 6218, and NGC 6981 is consistent with their dynamical ages, but only the global MF of M10 is consistent with its degree of mass segregation as well. NGC 5466 and NGC 6101 on the other hand appear to be less segregated than their dynamical ages would indicate. Furthermore, despite the fact that the escape rate of stars in non-segregated clusters is independent of stellar mass, both NGC 5466 and NGC 6101 have near-flat MFs. We discuss various mechanisms which could produce non-segregated clusters with near-flat MFs, including higher mass-loss rates and black hole retention, but argue that for some clusters (NGC 5466 and NGC 6101) explaining the present-day properties might require either a non-universal IMF or a much more complex dynamical history. △ Less

Submitted 20 July, 2017; originally announced July 2017.

Comments: 12 pages, 9 figures, Accepted for publication in MNRAS

The Early Evolution of Star Clusters in Compressive and Extensive Tidal Fields

Authors: Jeremy J. Webb, Saahil S. Patel, Enrico Vesperini

Abstract: We present N-body simulations of star clusters that initially evolve within a strong compressive tidal field and then transition into extensive tidal fields of varying strengths. While subject to compressive tides, clusters can undergo significant heating due to two-body interactions and mass loss due to stellar evolution. When the cluster transitions into an extensive tidal field it is super-viri… ▽ More We present N-body simulations of star clusters that initially evolve within a strong compressive tidal field and then transition into extensive tidal fields of varying strengths. While subject to compressive tides, clusters can undergo significant heating due to two-body interactions and mass loss due to stellar evolution. When the cluster transitions into an extensive tidal field it is super-virialized, which leads to a rapid expansion and significant mass loss before the cluster reaches virial equilibrium. After the transition, clusters are significantly less massive, more extended and therefore more tidally filling than clusters which have spent their entire lifetime in a similar extensive tidal field. △ Less

Submitted 10 February, 2017; originally announced February 2017.

Comments: 5 pages, 4 figures, Accepted for publication in MNRAS Letters

On the Link Between Energy Equipartition and Radial Variation in the Stellar Mass Function of Star Clusters

Authors: Jeremy J. Webb, Enrico Vesperini

Abstract: We make use of $N$-body simulations to determine the relationship between two observable parameters that are used to quantify mass segregation and energy equipartition in star clusters. Mass segregation can be quantified by measuring how the slope of a cluster's stellar mass function $α$ changes with clustercentric distance r, and then calculating $δ_α= \frac{d α(r)}{d ln(r/r_m)}$ where $r_m$ is t… ▽ More We make use of $N$-body simulations to determine the relationship between two observable parameters that are used to quantify mass segregation and energy equipartition in star clusters. Mass segregation can be quantified by measuring how the slope of a cluster's stellar mass function $α$ changes with clustercentric distance r, and then calculating $δ_α= \frac{d α(r)}{d ln(r/r_m)}$ where $r_m$ is the cluster's half-mass radius. The degree of energy equipartition in a cluster is quantified by $η$, which is a measure of how stellar velocity dispersion $σ$ depends on stellar mass m via $σ(m) \propto m^{-η}$. Through a suite of $N$-body star cluster simulations with a range of initial sizes, binary fractions, orbits, black hole retention fractions, and initial mass functions, we present the co-evolution of $δ_α$ and $η$. We find that measurements of the global $η$ are strongly affected by the radial dependence of $σ$ and mean stellar mass and the relationship between $η$ and $δ_α$ depends mainly on the cluster's initial conditions and the tidal field. Within $r_m$, where these effects are minimized, we find that $η$ and $δ_α$ initially share a linear relationship. However, once the degree of mass segregation increases such that the radial dependence of $σ$ and mean stellar mass become a factor within $r_m$, or the cluster undergoes core collapse, the relationship breaks down. We propose a method for determining $η$ within $r_m$ from an observational measurement of $δ_α$. In cases where $η$ and $δ_α$ can be measured independently, this new method offers a way of measuring the cluster's dynamical state. △ Less

Submitted 28 September, 2016; originally announced September 2016.

Comments: 9 pages, 5 figures, 2 tables, Accepted for publication in MNRAS

Radial Variation in the Stellar Mass Functions of Star Clusters

Authors: Jeremy J. Webb, Enrico Vesperini

Abstract: A number of recent observational studies of Galactic globular clusters have measured the variation in the slope of a cluster's stellar mass function $α$ with clustercentric distance $r$. In order to gather a deeper understanding of the information contained in such observations, we have explored the evolution of $α(r)$ for star clusters with a variety of initial conditions using a large suite of… ▽ More A number of recent observational studies of Galactic globular clusters have measured the variation in the slope of a cluster's stellar mass function $α$ with clustercentric distance $r$. In order to gather a deeper understanding of the information contained in such observations, we have explored the evolution of $α(r)$ for star clusters with a variety of initial conditions using a large suite of $N$-body simulations. We have specifically studied how the time evolution of $α(r)$ is affected by initial size, mass, binary fraction, primordial mass segregation, black hole retention, an external tidal field, and the initial mass function itself. Previous studies have shown that the evolution of $α_G$ is closely related to the amount of mass loss suffered by a cluster. Hence for each simulation we have also followed the evolution of the slope of the cluster's global stellar mass function, $α_G$, and have shown that clusters follow a well-defined track in the $α_G$-$dα(r)/d(ln(r/r_m))$ plane. The location of a cluster on the $α_G-dα(r)/d(ln(r/r_m))$ plane can therefore constrain its dynamical history and, in particular, constrain possible variations in the stellar initial mass function. The $α_G$-$dα(r)/d(ln(r/r_m))$ plane thus serves as a key tool for fully exploiting the information contained in wide field studies of cluster stellar mass functions. △ Less

Submitted 29 August, 2016; v1 submitted 25 August, 2016; originally announced August 2016.

Comments: 13 pages, 13 figures, 1 table, Accepted for publication in MNRAS

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

Authors: Jeremy J. Webb, Alison Sills, William E. Harris, Matías Gómez, Maurizio Paolillo, Kristin A. Woodley, Thomas H. Puzia

Abstract: 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… ▽ More 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. △ Less

Submitted 10 May, 2016; originally announced May 2016.

Comments: 16 pages, 7 figures, 4 tables, Accepted for publication in MNRAS

The Dynamical Evolution of Accreted Star Clusters in the Milky Way

Authors: Meghan Miholics, Jeremy J. Webb, Alison Sills

Abstract: We perform $N$-body simulations of star clusters in time-dependant galactic potentials. Since the Milky Way was built-up through mergers with dwarf galaxies, its globular cluster population is made up of clusters formed both during the initial collapse of the Galaxy and in dwarf galaxies that were later accreted. Throughout a dwarf-Milky Way merger, dwarf galaxy clusters are subject to a changing… ▽ More We perform $N$-body simulations of star clusters in time-dependant galactic potentials. Since the Milky Way was built-up through mergers with dwarf galaxies, its globular cluster population is made up of clusters formed both during the initial collapse of the Galaxy and in dwarf galaxies that were later accreted. Throughout a dwarf-Milky Way merger, dwarf galaxy clusters are subject to a changing galactic potential. Building on our previous work, we investigate how this changing galactic potential affects the evolution of a cluster's half mass radius. In particular, we simulate clusters on circular orbits around a dwarf galaxy that either falls into the Milky Way or evaporates as it orbits the Milky Way. We find that the dynamical evolution of a star cluster is determined by whichever galaxy has the strongest tidal field at the position of the cluster. Thus, clusters entering the Milky Way undergo changes in size as the Milky Way tidal field becomes stronger and that of the dwarf diminishes. We find that ultimately accreted clusters quickly become the same size as a cluster born in the Milky Way on the same orbit. Assuming their initial sizes are similar, clusters born in the Galaxy and those that are accreted cannot be separated based on their current size alone. △ Less

Submitted 15 November, 2015; v1 submitted 12 November, 2015; originally announced November 2015.

Comments: 9 pages, 4 figures, accepted for publication in MNRAS

The Dynamics of Multiple Populations in the Globular Cluster NGC 6362

Authors: Meghan Miholics, Jeremy J. Webb, Alison Sills

Abstract: 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 concentra… ▽ More 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. △ Less

Submitted 25 September, 2015; v1 submitted 8 September, 2015; originally announced September 2015.

Comments: 8 pages, 5 figures, Accepted for publication in MNRAS

Back to the Future: Estimating Initial Globular Cluster Masses from their Present Day Stellar Mass Functions

Authors: Jeremy J. Webb, Nathan W. C. Leigh

Abstract: We use N-body simulations to model the 12 Gyr evolution of a suite of star clusters with identical initial stellar mass functions over a range of initial cluster masses, sizes, and orbits. Our models reproduce the distribution of present-day global stellar mass functions that is observed in the Milky Way globular cluster population. We find that the slope of a star cluster's stellar mass function… ▽ More We use N-body simulations to model the 12 Gyr evolution of a suite of star clusters with identical initial stellar mass functions over a range of initial cluster masses, sizes, and orbits. Our models reproduce the distribution of present-day global stellar mass functions that is observed in the Milky Way globular cluster population. We find that the slope of a star cluster's stellar mass function is strongly correlated with the fraction of mass that the cluster has lost, independent of the cluster's initial mass, and nearly independent of its orbit and initial size. Thus, the mass function - initial mass relation can be used to determine a Galactic cluster's initial total stellar mass, if the initial stellar mass function is known. We apply the mass function - initial mass relation presented here to determine the initial stellar masses of 33 Galactic globular clusters, assuming an universal Kroupa initial mass function. Our study suggests that globular clusters had initial masses that were on average a factor of 4.5 times larger than their present day mass, with three clusters showing evidence for being 10 times more massive at birth. △ Less

Submitted 11 August, 2015; v1 submitted 3 August, 2015; originally announced August 2015.

Comments: Revised Version: 11 pages, 4 figures, Accepted for publication in MNRAS

The state of globular clusters at birth II: primordial binaries

Authors: Nathan W. C. Leigh, Mirek Giersz, Michael Marks, Jeremy J. Webb, Arkadiusz Hypki, Craig O. Heinke, Pavel Kroupa, Alison Sills

Abstract: (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 a… ▽ More (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. △ Less

Submitted 8 October, 2014; originally announced October 2014.

Comments: 17 pages; 5 figures; 1 table; accepted for publication in MNRAS

The Size of Star Clusters Accreted by the Milky Way

Authors: Meghan Miholics, Jeremy J. Webb, Alison Sills

Abstract: 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 decreas… ▽ More 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. △ Less

Submitted 17 September, 2014; originally announced September 2014.

Comments: 7 pages, 6 figures, Accepted for publication MNRAS

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

Authors: Jeremy J. Webb, Alison Sills, William E. Harris, Jarrod R. Hurley

Abstract: 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 e… ▽ More 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. △ Less

Submitted 2 September, 2014; originally announced September 2014.

Comments: 9 pages, 9 figures, Accepted for publication in MNRAS

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

Authors: Jeremy J. Webb, Nathan Leigh, Alison Sills, William E. Harris, Jarrod R. Hurley

Abstract: 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… ▽ More 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. △ Less

Submitted 4 July, 2014; v1 submitted 13 May, 2014; originally announced May 2014.

Comments: 10 pages, 7 figures, Accepted for publication in MNRAS

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

Authors: Jeremy J. Webb, Alison Sills, William E. Harris

Abstract: 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 agreemen… ▽ More 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. △ Less

Submitted 17 October, 2013; originally announced October 2013.

Comments: 10 pages, 11 figures, Accepted for publication in ApJ

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

Authors: Nathan W. C. Leigh, Mirek Giersz, Jeremy J. Webb, Arkadiusz Hypki, Guido de Marchi, Pavel Kroupa, Alison Sills

Abstract: 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 clus… ▽ More 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. △ Less

Submitted 18 November, 2014; v1 submitted 26 September, 2013; originally announced September 2013.

Comments: 16 pages, 5 figures, 2 tables; accepted for publication in MNRAS; edited to match proof corrections

The Influence of Orbital Eccentricity on Tidal Radii of Star Clusters

Authors: Jeremy J. Webb, William E. Harris, Alison Sills, Jarrod R. Hurley

Abstract: 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 t… ▽ More 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. △ Less

Submitted 3 January, 2013; originally announced January 2013.

Comments: 10 pages, 12 figures, Accepted for publication in ApJ

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

Authors: Jeremy J. Webb, William E. Harris, Alison Sills

Abstract: 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… ▽ More 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 △ Less

Submitted 9 October, 2012; originally announced October 2012.

Comments: 14 Pages, 5 figures, Accepted for publication in ApJL

The Observational and Theoretical Tidal Radii of Globular Clusters in M87

Authors: Jeremy J. Webb, Alison Sills, William E. Harris

Abstract: 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 a… ▽ More 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. △ Less

Submitted 4 January, 2012; originally announced January 2012.

Comments: 8 pages, 8 figures, Accepted for publication in ApJ