Disruption Of Globular Clusters Research Articles

Properties of Globular Clusters in Galaxy Clusters: Sensitivity from the Formation and Evolution of Globular Clusters

We investigate the properties of globular clusters (GCs) in a galaxy cluster, using the particle tagging method with a semianalytical approach in a cosmological context. We assume GCs form from dark matter halo mergers and their metallicity is assigned based on the stellar mass of the host dark matter halos and the formation redshift of GCs. Dynamical evolution and disruption of GCs are considered using semianalytical approaches, controlled by several free parameters. In this paper, we investigate how our results are changed by the choice of free parameters. We compare our fiducial results with representative observations, including the mass ratio between the GC system and its host galaxy, the GC occupancy, the number fraction of blue GCs, and the metallicity gradient with the GC mass. Because we can know the positions of GCs with time, comparison with additional observations is possible, e.g., the median radii of the GC system in individual galaxies, the mean projected density profiles of intracluster GCs, and the metallicity and age gradients of GCs with a clustercentric radius. We also find that the specific mass of the GC system in each galaxy is different with a clustercentric radius.

APOGEE discovery of a chemically atypical star disrupted from NGC 6723 and captured by the Milky Way bulge

The central (‘bulge’) region of the Milky Way is teeming with a significant fraction of mildly metal-deficient stars with atmospheres that are strongly enriched in cyanogen (12C14N). Some of these objects, which are also known as nitrogen-enhanced stars, are hypothesised to be relics of the ancient assembly history of the Milky Way. Although the chemical similarity of nitrogen-enhanced stars to the unique chemical patterns observed in globular clusters has been observed, a direct connection between field stars and globular clusters has not yet been proven. In this work, we report on high-resolution, near-infrared spectroscopic observations of the bulge globular cluster NGC 6723, and the serendipitous discovery of a star, 2M18594405−3651518, located outside the cluster (near the tidal radius) but moving on a similar orbit, providing the first clear piece of evidence of a star that was very likely once a cluster member and has recently been ejected. Its nitrogen abundance ratio ([N/Fe] ≳ + 0.94) is well above the typical Galactic field-star levels, and it exhibits noticeable enrichment in the heavys-process elements (Ce, Nd, and Yb), along with moderate carbon enrichment; all characteristics are known examples in globular clusters. This result suggests that some of the nitrogen-enhanced stars in the bulge likely originated from the tidal disruption of globular clusters.

An excess of globular clusters in Ultra-Diffuse Galaxies formed through tidal heating

ABSTRACT To investigate the origin of elevated globular cluster (GC) abundances observed around Ultra-Diffuse Galaxies (UDGs), we simulate GC populations hosted by UDGs formed through tidal heating. Specifically, GC formation is modelled as occurring in regions of dense star formation. Because star formation-rate densities are higher at high redshift, dwarf galaxies in massive galaxy clusters, which formed most of their stars at high redshift, form a large fraction of their stars in GCs. Given that UDGs formed through environmental processes are more likely to be accreted at high redshift, these systems have more GCs than non-UDGs. In particular, our model predicts that massive UDGs have twice the GC mass of non-UDGs of similar stellar mass, in rough agreement with observations. Although this effect is somewhat diminished by GC disruption, we find that the relationship between GC mass fraction and cluster-centric distance, and the relationship between GC mass fraction and galaxy half-light radius are remarkably similar to observations. Among our model objects, both UDGs and non-UDGs present a correlation between halo mass and GC mass, although UDGs have lower dynamical masses at a given GC mass. Furthermore, because of the effectiveness of GC disruption, we predict that GCs around UDGs should have a more top heavy mass function than GCs around non-UDGs. This analysis suggests that dwarfs with older stellar populations, such as UDGs, should have higher GC mass fractions than objects with young stellar populations, such as isolated dwarfs.

Discovery of Extended Tidal Tails around the Globular Cluster Palomar 13

Abstract We use photometry from the DECam Legacy Survey to detect candidate tidal tails extending ∼5° on either side of the Palomar 13 globular cluster. The tails are aligned with the proper motion of Palomar 13 and are consistent with its old, metal-poor stellar population. We identify three RR Lyrae stars (RRLs) that are plausibly associated with the tails, in addition to four previously known in the cluster. From these RRLs, we find that the mean distance to the cluster and tails is 23.6 ± 0.2 kpc and estimate the total (initial) luminosity of the cluster to be , consistent with previous claims that its initial luminosity was higher than its current luminosity. Combined with previously determined proper motion and radial velocity measurements of the cluster, we find that Palomar 13 is on a highly eccentric orbit (e ∼ 0.8) with a pericenter of ∼9 kpc and an apocenter of ∼69 kpc, and a recent pericentric passage of the cluster ∼75 Myr ago. We note a prominent linear structure in the interstellar dust map that runs parallel to the candidate tidal features, but conclude that reddening due to dust is unlikely to account for the structure that we observe. If confirmed, the Palomar 13 stellar stream would be one of very few streams with a known progenitor system, making it uniquely powerful for studying the disruption of globular clusters, the formation of the stellar halo, and the distribution of matter within our Galaxy.

The mass fraction of halo stars contributed by the disruption of globular clusters in the E-MOSAICS simulations

ABSTRACT Globular clusters (GCs) have been posited, alongside dwarf galaxies, as significant contributors to the field stellar population of the Galactic halo. In order to quantify their contribution, we examine the fraction of halo stars formed in stellar clusters in the suite of 25 present-day Milky Way-mass cosmological zoom simulations from the E-MOSAICS project. We find that a median of 2.3 and 0.3 per cent of the mass in halo field stars formed in clusters and GCs, defined as clusters more massive than 5 × 103 and 105 M⊙, respectively, with the 25–75th percentiles spanning 1.9–3.0 and 0.2–0.5 per cent being caused by differences in the assembly histories of the host galaxies. Under the extreme assumption that no stellar cluster survives to the present day, the mass fractions increase to a median of 5.9 and 1.8 per cent. These small fractions indicate that the disruption of GCs plays a subdominant role in the build-up of the stellar halo. We also determine the contributed halo mass fraction that would present signatures of light-element abundance variations considered to be unique to GCs, and find that clusters and GCs would contribute a median of 1.1 and 0.2 per cent, respectively. We estimate the contributed fraction of GC stars to the Milky Way halo, based on recent surveys, and find upper limits of 2–5 per cent (significantly lower than previous estimates), suggesting that models other than those invoking strong mass loss are required to describe the formation of chemically enriched stellar populations in GCs.

Kinematics of the Palomar 5 Stellar Stream from RR Lyrae Stars

Abstract Thin stellar streams, formed from the tidal disruption of globular clusters, are important gravitational tools, sensitive to both global and small-scale properties of dark matter. The Palomar 5 stellar stream (Pal 5) is an exemplar stream within the Milky Way: its ∼20° tidal tails connect back to the progenitor cluster, and the stream has been used to study the shape, total mass, and substructure fraction of the dark matter distribution of the Galaxy. However, most details of the phase-space distribution of the stream are not fully explained, and dynamical models that use the stream for other inferences are therefore incomplete. Here we aim to measure distance and kinematic properties along the Pal 5 stream in order to motivate improved models of the system. We use a large catalog of RR Lyrae-type stars (RRLs) with astrometric data from the Gaia mission to probabilistically identify RRLs in the Pal 5 stream. RRLs are useful because they are intrinsically luminous standard candles and their distances can be inferred with small relative precision (∼3%). By building a probabilistic model of the Pal 5 cluster and stream in proper motion and distance, we find 27 RRLs consistent with being members of the cluster (10) and stream (17). Using these RRLs, we detect gradients in distance and proper motion along the stream, and provide an updated measurement of the distance to the Pal 5 cluster using the RRLs, . We provide a catalog of Pal 5 RRLs with inferred membership probabilities for future modeling work.

Black Hole Mergers from an Evolving Population of Globular Clusters.

The high rate of black hole (BH) mergers detected by LIGO/Virgo opened questions on their astrophysical origin. One possibility is the dynamical channel, in which binary formation and hardening is catalyzed by dynamical encounters in globular clusters (GCs). Previous studies have shown that the BH merger rate from the present day GC density in the Universe is lower than the observed rate. In this Letter, we study the BH merger rate by accounting for the first time for the evolution of GCs within their host galaxies. The mass in GCs was initially ∼8×higher, which decreased to its present value due to evaporation and tidal disruption. Many BH binaries that were ejected long before their merger originated in GCs that no longer exist. We find that the comoving merger rate in the dynamical channel from GCs varies between 18 to 35 Gpc^{-3} yr^{-1} between redshift z=0.5 to 2, and the total rate is 1, 5, 24 events per day within z=0.5, 1, and 2, respectively. The cosmic evolution and disruption of GCs systematically increases the present-day merger rate by a factor ∼2 relative to isolated clusters. Gravitational wave detector networks offer an unique observational probe of the initial number of GC populations and their subsequent evolution across cosmic time.

The Little Engines That Could? Globular clusters contribute significantly to reionization-era star formation.

Metal-poor globular clusters (GCs) are both numerous and ancient, which indicates that they may be important contributors to ionizing radiation in the reionization era. Starting from the observed number density and stellar mass function of old GCs at z = 0, I compute the contribution of GCs to ultraviolet luminosity functions (UVLFs) in the high-redshift Universe (10 ≳ z ≳ 4). Even under absolutely minimal assumptions - no disruption of GCs and no reduction in GC stellar mass from early times to the present - GC star formation contributes non-negligibly to the UVLF at luminosities that are accessible to the Hubble Space Telescope (HST, M 1500 ≈ -17). If the stellar masses of GCs were significantly higher in the past, as is predicted by most models explaining GC chemical anomalies, then GCs dominate the UV emission from many galaxies in existing deep-field observations. On the other hand, it is difficult to reconcile observed UVLFs with models requiring stellar masses at birth that exceed present-day stellar masses by more than a factor of 5. The James Webb Space Telescope will be able to directly detect individual GCs at z ∼ 6 in essentially all bright galaxies, and many galaxies below the knee of the UVLF, for most of the scenarios considered here. The properties of a subset of high-redshift sources with -19 ≳ M 1500 ≲ -14 in HST lensing fields indicate that they may actually be GCs in formation.

Dating the Tidal Disruption of Globular Clusters with GAIA Data on Their Stellar Streams

Abstract The Gaia mission promises to deliver precision astrometry at an unprecedented level, heralding a new era for discerning the kinematic and spatial coordinates of stars in our Galaxy. Here, we present a new technique for estimating the age of tidally disrupted globular cluster streams using the proper motions and parallaxes of tracer stars. We evolve the collisional dynamics of globular clusters within the evolving potential of a Milky Way-like halo extracted from a cosmological ΛCDM simulation and analyze the resultant streams as they would be observed by Gaia. The simulations sample a variety of globular cluster orbits, and account for stellar evolution and the gravitational influence of the disk of the Milky Way. We show that a characteristic timescale, obtained from the dispersion of the proper motions and parallaxes of stars within the stream, is a good indicator for the time elapsed since the stream has been freely expanding away due to the tidal disruption of the globular cluster. This timescale, in turn, places a lower limit on the age of the cluster. The age can be deduced from astrometry using a modest number of stars, with the error on this estimate depending on the proximity of the stream and the number of tracer stars used.

Disrupted globular clusters and the gamma-ray excess in the Galactic Centre

The Fermi Large Area Telescope has provided the most detailed view toward the Galactic Centre (GC) in high-energy gamma rays. Besides the interstellar emission and point-source contributions, the data suggest a residual diffuse gamma-ray excess. The similarity of its spatial distribution with the expected profile of dark matter has led to claims that this may be evidence for dark matter particle annihilation. Here, we investigate an alternative explanation that the signal originates from millisecond pulsars (MSPs) formed in dense globular clusters and deposited at the GC as a consequence of cluster inspiral and tidal disruption. We use a semi-analytical model to calculate the formation, migration, and disruption of globular clusters in the Galaxy. Our model reproduces the mass of the nuclear star cluster and the present-day radial and mass distribution of globular clusters. For the first time, we calculate the evolution of MSPs from disrupted globular clusters throughout the age of the Galaxy and consistently include the effect of the MSP spin-down due to magnetic-dipole breaking. The final gamma-ray amplitude and spatial distribution are in good agreement with the Fermi observations and provide a natural astrophysical explanation for the GC excess.

A synoptic map of halo substructures from the Pan-STARRS1 3π survey

We present a panoramic map of the entire Milky Way halo north of dec~-30 degrees (~30,000 deg^2), constructed by applying the matched-filter technique to the Pan-STARRS1 3Pi Survey dataset. Using single-epoch photometry reaching to g~22, we are sensitive to stellar substructures with heliocentric distances between 3.5 and ~35 kpc. We recover almost all previously-reported streams in this volume and demonstrate that several of these are significantly more extended than earlier datasets have indicated. In addition, we also report five new candidate stellar streams. One of these features appears significantly broader and more luminous than the others and is likely the remnant of a dwarf galaxy. The other four streams are consistent with a globular cluster origin, and three of these are rather short in projection (<10 degrees), suggesting that streams like Ophiuchus may not be that rare. Finally, a significant number of more marginal substructures are also revealed by our analysis; many of these features can also be discerned in matched-filter maps produced by other authors from SDSS data, and hence they are very likely to be genuine. However, the extant 3Pi data is currently too shallow to determine their properties or produce convincing CMDs. The global view of the Milky Way provided by Pan-STARRS1 provides further evidence for the important role of both globular cluster disruption and dwarf galaxy accretion in building the Milky Way's stellar halo.

The gamma-ray pulsar population of globular clusters: implications for the GeV excess

It has been suggested that the GeV excess, observed from the region surrounding the Galactic Center, might originate from a population of millisecond pulsars that formed in globular clusters. With this in mind, we employ the publicly available Fermi data to study the gamma-ray emission from 157 globular clusters, identifying a statistically significant signal from 25 of these sources (ten of which are not found in existing gamma-ray catalogs). We combine these observations with the predicted pulsar formation rate based on the stellar encounter rate of each globular cluster to constrain the gamma-ray luminosity function of millisecond pulsars in the Milky Way's globular cluster system. We find that this pulsar population exhibits a luminosity function that is quite similar to those millisecond pulsars observed in the field of the Milky Way (i.e. the thick disk). After pulsars are expelled from a globular cluster, however, they continue to lose rotational kinetic energy and become less luminous, causing their luminosity function to depart from the steady-state distribution. Using this luminosity function and a model for the globular cluster disruption rate, we show that millisecond pulsars born in globular clusters can account for only a few percent or less of the observed GeV excess. Among other challenges, scenarios in which the entire GeV excess is generated from such pulsars are in conflict with the observed mass of the Milky Way's Central Stellar Cluster.

Chaotic dispersal of tidal debris

Several long, dynamically cold stellar streams have been observed around the Milky Way Galaxy, presumably formed from the tidal disruption of globular clusters. In integrable potentials---where all orbits are regular---tidal debris phase-mixes close to the orbit of the progenitor system. However, the Milky Way's dark matter halo is expected not to be fully integrable; an appreciable fraction of orbits will be chaotic. This paper examines the influence of chaos on the phase-space morphology of cold tidal streams. Streams even in weakly chaotic regions look very different from those in regular regions. We find that streams can be sensitive to chaos on a much shorter time-scale than any standard prediction (from the Lyapunov or frequency-diffusion times). For example, on a weakly chaotic orbit with a chaotic timescale predicted to be $>$1000 orbital periods ($>$1000 Gyr), the resulting stellar stream is, after just a few 10's of orbits, substantially more diffuse than any formed on a nearby but regular orbit. We find that the enhanced diffusion of the stream stars can be understood by looking at the variance in orbital frequencies of orbit ensembles centered around the parent (progenitor) orbit. Our results suggest that long, cold streams around our Galaxy must exist only on regular (or very nearly regular) orbits; they potentially provide a map of the regular regions of the Milky Way potential. This suggests a promising new direction for the use of tidal streams to constrain the distribution of dark matter around our Galaxy.

Tidal disruption of globular clusters in dwarf galaxies with triaxial dark matter haloes

We use N-body simulations to study the tidal evolution of globular clusters (GCs) in dwarf spheroidal (dSph) galaxies. Our models adopt a cosmologically motivated scenario in which the dSph is approximated by a static NFW halo with a triaxial shape. We apply our models to five GCs spanning three orders of magnitude in stellar density and two in mass, chosen to represent the properties exhibited by the five GCs of the Fornax dSph. We show that only the object representing Fornax's least dense GC (F1) can be fully disrupted by Fornax's internal tidal field--the four denser clusters survive even if their orbits decay to the centre of Fornax. For a large set of orbits and projection angles we examine the spatial and velocity distribution of stellar debris deposited during the complete disruption of an F1-like GC. Our simulations show that such debris appears as shells, isolated clumps and elongated over-densities at low surface brightness (>26 mag/arcsec^2), reminiscent of substructure observed in several MW dSphs. Such features arise from the triaxiality of the galaxy potential and do not dissolve in time. The kinematics of the debris depends strongly on the progenitor's orbit. Debris associated with box and resonant orbits does not display stream motions and may appear "colder"/"hotter" than the dSph's field population if the viewing angle is perpendicular/parallel to progenitor's orbital plane. In contrast, debris associated with loop orbits shows a rotational velocity that may be detectable out to few kpc from the galaxy centre. Chemical tagging that can distinguish GC debris from field stars may reveal whether the merger of GCs contributed to the formation of multiple stellar components observed in dSphs.

The Birthplace of Low‐Mass X‐Ray Binaries: Field Versus Globular Cluster Populations

Recent Chandra studies of low-mass X-ray binaries (LMXBs) within early-type galaxies have found that LMXBs are commonly located within globular clusters of the galaxies. However, whether all LMXBs are formed within globular clusters has remained an open question. If all LMXBs formed within globular clusters, the summed X-ray luminosity of the LMXBs in a galaxy should be directly proportional to the number of globular clusters in the galaxy regardless of where the LMXBs currently reside. We have compared these two quantities over the same angular area for a sample of 12 elliptical and S0 galaxies observed with Chandra and found that the correlation between the two quantities is weaker than expected if all LMXBs formed within globular clusters. This indicates that a significant number of the LMXBs were formed in the field and naturally accounts for the spread in field-to-cluster fractions of LMXBs from galaxy to galaxy. We also find that the pollution of globular cluster LMXBs into the field has been minimal within elliptical galaxies, but there is evidence that roughly half of the LMXBs originally in the globular clusters of S0 galaxies in our sample have escaped into the field. This is likely due to higher globular cluster disruption rates in S0s, resulting from stronger gravitational shocks caused by the passage of globular clusters through the disks of S0 galaxies that are absent in elliptical galaxies.

A bound on galactic mass loss rates obtained from globular cluster dynamics

Using tidal disruption of globular clusters by the galactic centre, we put limits on the total mass enclosed into orbits of observed globular clusters. Under the assumption that the rate of mass loss from the galaxy is steady, we then deduce a bound on this rate. In particular, this bound can be used to constrain the galactic gravitational wave luminosity.

Dynamical Constraints on Dark Matter in Compact Objects

Many of the baryons in the universe are dark and at least some of the dark baryons could be in the form of compact objects. Such objects could be in various locations—galactic disks, galactic halos, clusters of galaxies, or intergalactic space—and each of these is associated with a dark matter problem. For each site we consider the various dynamical constraints that can be placed on the fraction of the dark matter in compact objects of different mass. Small compact objects in the Galaxy are constrained byupper limits on their encounter rate with the Earth and solar system since they would resemble meteors or comets. Larger objects are constrained by the disruptive or disturbing effects they would have on various astronomical systems. For disk objects, the most interesting constraints come from the disruption of binary stars or open star clusters. For halo objects, they come from the disruption of globular clusters, the heating of the Galactic disk, and their accumulation in the Galactic nucleus as a result of dynamical friction. For cluster objects, they come from the tidal distortion and disruption of galaxies. For intergalactic objects, they come from the upper limit on the peculiar motions induced in galaxies. We also apply these limits to the situation in which the compact objects are clusters of smaller objects.

On the evolution of the Galactic globular cluster system

In this paper we address the issue of the origin of some observational properties of the galactic globular cluster system. After a preliminary study of some general properties of the main evolutionary processes, we investigate the evolution of systems of globular clusters located in a model of the Milky Way starting from different initial conditions. We study the role of the evolutionary processes in changing the spatial distribution and mass function of the cluster system, in establishing and/or preserving some of the observed correlations and trends between internal properties of globular clusters and between internal properties and location inside the host galaxy and we provide an estimate for the rates of core collapse and disruption of globular clusters. The initial mass function and spatial distribution of the cluster system evolve quite significantly in one Hubble time and the evolution is toward a final state similar to the observed one. If the mass function is initially taken to be a log-normal distribution similar to the one currently observed in our galaxy, its shape is not significantly altered during the entire simulation even though a significant number of clusters are disrupted before one Hubble time, which suggests that the present mass function might represent a sort of 'quasi-equilibrium' distribution.

X-Ray Binaries in Globular Clusters

X-ray binaries in globular clusters provide a powerful tool for the exploration of the evolution of compact binaries and their host globular clusters. Recent x-ray and optical studies of these systems have yielded long-sought binary periods and fundamental properties for two sources (in NGC 6624 and M 15). It appears that tidal capture formation of compact binaries in globular clusters can proceed by several different routes and lead to exotic systems such as the white dwarf-neutron star binary with an 11-minute period recently discovered in NGC 6624. Combined with previously reported long-term periods for several globular cluster (and field) x-ray sources, this suggests again that many of these systems may in fact be hierarchical triple systems. The prospects for forming these in the dense cores of clusters undergoing core collapse is discussed, and searches for color gradients in the cores of globular clusters showing cusps in their central surface brightness distribution are presented. A program to test for the high central density of binaries (and triples) expected in cusp clusters by searching for diffuse line emission from their constituent cataclysmic variables is briefly described. Finally, the case for globular cluster disruption and the formation of galactic x-ray burst source is reviewed in light of recent developments.

Tidal Stripping and Disruption of Globular Clusters

Tidal stripping and disruption of globular clusters may be responsible for the absence of low density clusters in the inner region of the Galaxy. We have studied these processes by integrating orbits of stars while the cluster is moving through a spherically symmetric galactic potential with constant circular velocity (Vcir =220 km/s). The response of the cluster to the tidal field of the galaxy is calculated in a selfconsistent manner with a collisionless N-body code with N=5000 (van Albada 1982, van Albada and Bontekoe in preparation).