Source author record

Jeremy J. Webb

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

19 published item(s)

preprint2021arXiv

Bayesian Inference of Globular Cluster Properties Using Distribution Functions

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.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2016arXiv

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

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

preprint2016arXiv

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

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.

preprint2016arXiv

Radial Variation in the Stellar Mass Functions of Star Clusters

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.

preprint2015arXiv

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

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.

preprint2015arXiv

The Dynamics of Multiple Populations in the Globular Cluster NGC 6362

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

preprint2014arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

The Size of Star Clusters Accreted by the Milky Way

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

preprint2014arXiv

The state of globular clusters at birth II: primordial binaries

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

preprint2014arXiv

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

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

preprint2013arXiv

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

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

preprint2013arXiv

The Influence of Orbital Eccentricity on Tidal Radii of Star Clusters

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

preprint2012arXiv

The Observational and Theoretical Tidal Radii of Globular Clusters in M87

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

preprint2012arXiv

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

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