-
The Fates of Merging Supermassive Black Holes and a Proposal for a New Class of X-Ray Sources
Authors:
Charles Zivancev,
Jeremiah Ostriker,
Andreas H. W. Kupper
Abstract:
We perform N-body simulations on some of the most massive galaxies extracted from a cosmological simulation of hierarchical structure formation with total masses in the range $10^{12} M_{\odot} < M_{tot} < 3\times 10^{13} M_{\odot}$ from $4\geq z \geq 0$. After galactic mergers, we track the dynamical evolution of the infalling black holes (BHs) around their host's central BHs. From 11 different s…
▽ More
We perform N-body simulations on some of the most massive galaxies extracted from a cosmological simulation of hierarchical structure formation with total masses in the range $10^{12} M_{\odot} < M_{tot} < 3\times 10^{13} M_{\odot}$ from $4\geq z \geq 0$. After galactic mergers, we track the dynamical evolution of the infalling black holes (BHs) around their host's central BHs. From 11 different simulations, we find that, of the 86 infalling BHs with masses > $10^4 M_{\odot}$, 36 merge with their host's central BH, 13 are ejected from their host galaxy, and 37 are still orbiting at $z=0$. Across all galaxies, 33 BHs are kicked to a higher orbit after close interactions with the central BH binary or multiple, after which only one of them merged with their hosts. These orbiting BHs should be detectable by their anomalous (not Low Mass X-ray Binary) spectra. The X-ray luminosities of the orbiting massive BHs at z=0 are in the range $10^{28}-10^{43}$ $\mathrm{erg}~\mathrm{s}^{-1}$, with a currently undetectable median value of $10^{33}$ $\mathrm{erg}~\mathrm{s}^{-1}$. However, the most luminous $\sim$5\% should be detectable by existing X-ray facilities.
△ Less
Submitted 19 August, 2020; v1 submitted 13 April, 2020;
originally announced April 2020.
-
Quantifying tidal stream disruption in a simulated Milky Way
Authors:
Emily Sandford,
Andreas H. W. Kuepper,
Kathryn V. Johnston,
Juerg Diemand
Abstract:
Simulations of tidal streams show that close encounters with dark matter subhalos induce density gaps and distortions in on-sky path along the streams. Accordingly, observing disrupted streams in the Galactic halo would substantiate the hypothesis that dark matter substructure exists there, while in contrast, observing collimated streams with smoothly varying density profiles would place strong up…
▽ More
Simulations of tidal streams show that close encounters with dark matter subhalos induce density gaps and distortions in on-sky path along the streams. Accordingly, observing disrupted streams in the Galactic halo would substantiate the hypothesis that dark matter substructure exists there, while in contrast, observing collimated streams with smoothly varying density profiles would place strong upper limits on the number density and mass spectrum of subhalos. Here, we examine several measures of stream "disruption" and their power to distinguish between halo potentials with and without substructure and with different global shapes. We create and evolve a population of 1280 streams on a range of orbits in the Via Lactea II simulation of a Milky Way-like halo, replete with a full mass range of ΛCDM subhalos, and compare it to two control stream populations evolved in smooth spherical and smooth triaxial potentials, respectively. We find that the number of gaps observed in a stellar stream is a poor indicator of the halo potential, but that (i) the thinness of the stream on-sky, (ii) the symmetry of the leading and trailing tails, and (iii) the deviation of the tails from a low-order polynomial path on-sky ("path regularity") distinguish between the three potentials more effectively. We find that globular cluster streams on low-eccentricity orbits far from the galactic center (apocentric radius ~ 30-80 kpc) are most powerful in distinguishing between the three potentials. If they exist, such streams will shortly be discoverable and mapped in high dimensions with near-future photometric and spectroscopic surveys.
△ Less
Submitted 19 May, 2017;
originally announced May 2017.
-
Direct N-body simulations of globular clusters -- III. Palomar\,4 on an eccentric orbit
Authors:
Akram Hasani Zonoozi,
Hosein Haghi,
Pavel Kroupa,
Andreas H. W. Küpper,
Holger Baumgardt
Abstract:
Palomar 4 is a low-density globular cluster with a current mass $\approx30000 M_{\odot}$ in the outer halo of the Milky Way with a two-body relaxation time of the order of a Hubble time. Yet, it is strongly mass segregated and contains a stellar mass function depleted of low-mass stars. Pal 4 was either born this way or it is a result of extraordinary dynamical evolution. Since two-body relaxation…
▽ More
Palomar 4 is a low-density globular cluster with a current mass $\approx30000 M_{\odot}$ in the outer halo of the Milky Way with a two-body relaxation time of the order of a Hubble time. Yet, it is strongly mass segregated and contains a stellar mass function depleted of low-mass stars. Pal 4 was either born this way or it is a result of extraordinary dynamical evolution. Since two-body relaxation cannot explain these signatures alone, enhanced mass loss through tidal shocking may have had a strong influence on Pal 4. Here, we compute a grid of direct N-body simulations to model Pal 4 on various eccentric orbits within the Milky Way potential to find likely initial conditions that reproduce its observed mass, half-light radius, stellar MF-slope and line-of-sight velocity dispersion. We find that Pal 4 is most likely orbiting on an eccentric orbit with an eccentricity of $e\approx 0.9$ and pericentric distance of $R_p\approx5$ kpc. In this scenario, the required 3D half-mass radius at birth is similar to the average sizes of typical GCs ($R_h\approx4-5$ pc), while its birth mass is about $M_0\approx10^5 M_{\odot}$. We also find a high degree of primordial mass segregation among the cluster stars, which seems to be necessary in every scenario we considered. Thus, using the tidal effect to constrain the perigalactic distance of the orbit of Pal 4, we predict that the proper motion of Pal 4 should be in the range $-0.52\leqμ_δ\leq-0.38$ mas\,yr$^{-1}$ and $-0.30\leqμ_{α\cosδ}\leq-0.15$ mas\,yr$^{-1}$.
△ Less
Submitted 22 January, 2017;
originally announced January 2017.
-
Exploding Satellites -- The Tidal Debris of the Ultra-Faint Dwarf Galaxy Hercules
Authors:
Andreas H. W. Küpper,
Kathryn V. Johnston,
Steffen Mieske,
Michelle L. M. Collins,
Erik J. Tollerud
Abstract:
The ultra-faint satellite galaxy Hercules has a strongly elongated and irregular morphology with detections of tidal features up to 1.3 deg (3 kpc) from its center. This suggests that Hercules may be dissolving under the Milky Way's gravitational influence, and hence could be a tidal stream in formation rather than a bound, dark-matter dominated satellite. Using Bayesian inference in combination w…
▽ More
The ultra-faint satellite galaxy Hercules has a strongly elongated and irregular morphology with detections of tidal features up to 1.3 deg (3 kpc) from its center. This suggests that Hercules may be dissolving under the Milky Way's gravitational influence, and hence could be a tidal stream in formation rather than a bound, dark-matter dominated satellite. Using Bayesian inference in combination with N-body simulations, we show that Hercules has to be on a very eccentric orbit (epsilon~0.95) within the Milky Way in this scenario. On such an orbit, Hercules "explodes" as a consequence of the last tidal shock at pericenter 0.5 Gyr ago. It is currently decelerating towards apocenter of its orbit with a velocity of V=157 km/s -- of which 99% is directed radially outwards. Due to differential orbital precession caused by the non-spherical nature of the Galactic potential, its debris fans out nearly perpendicular to its orbit. This explains why Hercules has an elongated shape without showing a distance gradient along its main body: it is in fact a stream that is significantly broader than it is long. In other words, it is moving perpendicular to its apparent major axis. In this scenario, there is a spike in the radial velocity profile created by the dominant debris component that formed through the last pericenter passage. This is similar to kinematic substructure that is observed in the real Hercules. Modeling a satellite on such a highly eccentric orbit is strongly dependent on the form of the Galactic potential. We therefore propose that detailed kinematic investigation of Hercules and other exploding satellite candidates can yield strong constraints on the potential of the Milky Way.
△ Less
Submitted 17 August, 2016;
originally announced August 2016.
-
Formation of Massive Black Holes in Galactic Nuclei: Runaway Tidal Encounters
Authors:
Nicholas C. Stone,
Andreas H. W. Kuepper,
Jeremiah P. Ostriker
Abstract:
Nuclear star clusters (NSCs) and supermassive black holes (SMBHs) both inhabit galactic nuclei, coexisting in a range of bulge masses, but excluding each other in the largest or smallest galaxies. We propose that the transformation of NSCs into SMBHs occurs via runaway tidal captures, once NSCs exceed a certain critical central density and velocity dispersion. The bottleneck in this process, as wi…
▽ More
Nuclear star clusters (NSCs) and supermassive black holes (SMBHs) both inhabit galactic nuclei, coexisting in a range of bulge masses, but excluding each other in the largest or smallest galaxies. We propose that the transformation of NSCs into SMBHs occurs via runaway tidal captures, once NSCs exceed a certain critical central density and velocity dispersion. The bottleneck in this process, as with all collisional runaways, is growing the first e-fold in black hole mass. The growth of a stellar mass black hole past this bottleneck occurs as tidally captured stars are consumed in repeated episodes of mass transfer at pericenter. Tidal captures may turn off as a growth channel once the black hole reaches a mass ~100-1000 solar masses, but tidal disruption events will continue and appear capable of growing the seed SMBH to larger sizes. The runaway slows (becomes sub-exponential) once the seed SMBH consumes the core of its host NSC. While the bulk of the cosmic mass density in SMBHs is ultimately produced (via the Soltan-Paczynski argument) by episodic gaseous accretion in very massive galaxies, the smallest SMBHs have probably grown from strong tidal encounters with NSC stars. SMBH seeds that grow for a time $t$ entirely through this channel will follow simple power law relations with the velocity dispersion, $σ$, of their host galaxy. In the simplest regime it is $M_\bullet \sim σ^{3/2}\sqrt{M_\star t / G} \sim 10^{6}M_\odot (σ/ 50~{\rm km~s}^{-1})^{3/2}(t/10^{10}~{\rm yr})^{1/2}$, but the exponents and prefactor can differ slightly depending on the details of loss cone refilling. Current tidal disruption event rates predicted from this mechanism are consistent with observations.
△ Less
Submitted 6 June, 2016;
originally announced June 2016.
-
Dynamical Evolution of Outer-Halo Globular Clusters
Authors:
Andreas H. W. Küpper,
Akram H. Zonoozi,
Hosein Haghi,
Nora Luetzgendorf,
Steffen Mieske,
Matthias Frank,
Holger Baumgardt,
Pavel Kroupa
Abstract:
Outer-halo globular clusters show large half-light radii and flat stellar mass functions, depleted in low-mass stars. Using N-body simulations of globular clusters on eccentric orbits within a Milky Way-like potential, we show how a cluster's half-mass radius and its mass function develop over time. The slope of the central mass function flattens proportionally to the amount of mass a cluster has…
▽ More
Outer-halo globular clusters show large half-light radii and flat stellar mass functions, depleted in low-mass stars. Using N-body simulations of globular clusters on eccentric orbits within a Milky Way-like potential, we show how a cluster's half-mass radius and its mass function develop over time. The slope of the central mass function flattens proportionally to the amount of mass a cluster has lost, and the half-mass radius grows to a size proportional to the average strength of the tidal field. The main driver of these processes is mass segregation of dark remnants. We conclude that the extended, depleted clusters observed in the Milky Way must have had small half-mass radii in the past, and that they expanded due to the weak tidal field they spend most of their lifetime in. Moreover, their mass functions must have been steeper in the past but flattened significantly as a cause of mass segregation and tidal mass loss.
△ Less
Submitted 18 December, 2015;
originally announced December 2015.
-
Globular Cluster Streams as Galactic High-Precision Scales
Authors:
A. H. W. Küpper,
E. Balbinot,
A. Bonaca,
K. V. Johnston,
D. W. Hogg,
P. Kroupa,
B. X. Santiago
Abstract:
Tidal streams of globular clusters are ideal tracers of the Galactic gravitational potential. Compared to the few known, complex and diffuse dwarf-galaxy streams, they are kinematically cold, have thin morphologies and are abundant in the halo of the Milky Way. Their coldness and thinness in combination with potential epicyclic substructure in the vicinity of the stream progenitor turns them into…
▽ More
Tidal streams of globular clusters are ideal tracers of the Galactic gravitational potential. Compared to the few known, complex and diffuse dwarf-galaxy streams, they are kinematically cold, have thin morphologies and are abundant in the halo of the Milky Way. Their coldness and thinness in combination with potential epicyclic substructure in the vicinity of the stream progenitor turns them into high-precision scales. With the example of Palomar 5, we demonstrate how modeling of a globular cluster stream allows us to simultaneously measure the properties of the disrupting globular cluster, its orbital motion, and the gravitational potential of the Milky Way.
△ Less
Submitted 30 September, 2015; v1 submitted 30 September, 2015;
originally announced September 2015.
-
Chaotic Dispersal of Tidal Debris
Authors:
Adrian M. Price-Whelan,
Kathryn V. Johnston,
Monica Valluri,
Sarah Pearson,
Andreas H. W. Kupper,
David W. Hogg
Abstract:
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 or…
▽ More
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.
△ Less
Submitted 6 April, 2016; v1 submitted 30 July, 2015;
originally announced July 2015.
-
Globular Cluster Streams as Galactic High-Precision Scales - The Poster Child Palomar 5
Authors:
Andreas H. W. Küpper,
Eduardo Balbinot,
Ana Bonaca,
Kathryn V. Johnston,
David W. Hogg,
Pavel Kroupa,
Basilio X. Santiago
Abstract:
Using the example of the tidal stream of the Milky Way globular cluster Palomar 5 (Pal 5), we demonstrate how observational data on streams can be efficiently reduced in dimensionality and modeled in a Bayesian framework. Our approach combines detection of stream overdensities by a Difference-of-Gaussians process with fast streakline models, a continuous likelihood function built from these models…
▽ More
Using the example of the tidal stream of the Milky Way globular cluster Palomar 5 (Pal 5), we demonstrate how observational data on streams can be efficiently reduced in dimensionality and modeled in a Bayesian framework. Our approach combines detection of stream overdensities by a Difference-of-Gaussians process with fast streakline models, a continuous likelihood function built from these models, and inference with MCMC. By generating $\approx10^7$ model streams, we show that the geometry of the Pal 5 debris yields powerful constraints on the solar position and motion, the Milky Way and Pal 5 itself. All 10 model parameters were allowed to vary over large ranges without additional prior information. Using only SDSS data and a few radial velocities from the literature, we find that the distance of the Sun from the Galactic Center is $8.30\pm0.25$ kpc, and the transverse velocity is $253\pm16$ km/s. Both estimates are in excellent agreement with independent measurements of these quantities. Assuming a standard disk and bulge model, we determine the Galactic mass within Pal 5's apogalactic radius of 19 kpc to be $(2.1\pm0.4)\times10^{11}$ M$_\odot$. Moreover, we find the potential of the dark halo with a flattening of $q_z = 0.95^{+0.16}_{-0.12}$ to be essentially spherical within the radial range that is effectively probed by Pal 5. We also determine Pal 5's mass, distance and proper motion independently from other methods, which enables us to perform vital cross-checks. We conclude that with more observational data and by using additional prior information, the precision of this method can be significantly increased.
△ Less
Submitted 9 February, 2015;
originally announced February 2015.
-
Tidal Stream Morphology as an Indicator of Dark Matter Halo Geometry: the Case of Palomar 5
Authors:
Sarah Pearson,
Andreas H. W. Küpper,
Kathryn V. Johnston,
Adrian M. Price-Whelan
Abstract:
This paper presents an example where the morphology of a single stellar stream can be used to rule out a specific galactic potential form without the need for velocity information. We investigate the globular cluster Palomar5 (Pal 5), which is tidally disrupting into a cold, thin stream mapped over 22 degrees on the sky with a typical width of 0.7 degrees. We generate models of this stream by fixi…
▽ More
This paper presents an example where the morphology of a single stellar stream can be used to rule out a specific galactic potential form without the need for velocity information. We investigate the globular cluster Palomar5 (Pal 5), which is tidally disrupting into a cold, thin stream mapped over 22 degrees on the sky with a typical width of 0.7 degrees. We generate models of this stream by fixing Pal 5's present-day position, distance and radial velocity via observations, while allowing its proper motion to vary. In a spherical dark matter halo we easily find models that fit the observed morphology. However, no plausible Pal 5 model could be found in the triaxial potential of Law & Majewski (2010), which has been proposed to explain the properties of the Sagittarius stream. In this case, the long, thin and curved morphology of the Pal5 stream alone can be used to rule out such a potential configuration. Pal 5 like streams in this potential are either too straight, missing the curvature of the observations, or show an unusual morphology which we dub stream-fanning: a signature sensitive to the triaxiality of a potential. We conclude that the mere existence of other thin tidal streams must provide broad constraints on the orientation and shape of the dark matter halo they inhabit.
△ Less
Submitted 31 March, 2016; v1 submitted 13 October, 2014;
originally announced October 2014.
-
The Puzzling Negative Orbit-Period Derivative of the Low-Mass X-Ray Binary 4U 1820-30 in NGC 6624
Authors:
Miklos Peuten,
Michael Brockamp,
Andreas H. W. Kuepper,
Pavel Kroupa
Abstract:
4U 1820-30 is a low-mass X-ray binary near the center of the globular cluster NGC 6624 consisting of, at least, one neutron star and one helium white dwarf. Analyzing 16 years of data from the Rossi X-ray Timing Explorer (RXTE) allows us to measure its orbital period and its time derivative with unprecedented accuracy to be P = 685.01197 +- 0.00003 s and dP/dt /P = -5.3 +- 0.3x10^-8 yr^-1. Hence,…
▽ More
4U 1820-30 is a low-mass X-ray binary near the center of the globular cluster NGC 6624 consisting of, at least, one neutron star and one helium white dwarf. Analyzing 16 years of data from the Rossi X-ray Timing Explorer (RXTE) allows us to measure its orbital period and its time derivative with unprecedented accuracy to be P = 685.01197 +- 0.00003 s and dP/dt /P = -5.3 +- 0.3x10^-8 yr^-1. Hence, we confirm that the period derivative is significantly negative at the >17 sigma level, contrary to theoretical expectations for an isolated X-ray binary. We discuss possible scenarios that could explain this discrepancy, and conclude that the center of NGC 6624 most likely contains large amounts of non-luminous matter such as dark remnants. We also discuss the possibility of an IMBH inside NGC 6624, or that a dark remnant close to 4U 1820-30 causes the observed shift.
△ Less
Submitted 21 December, 2014; v1 submitted 16 September, 2014;
originally announced September 2014.
-
The effect of primordial mass segregation on the size scale of globular clusters
Authors:
Hosein Haghi,
Seyed Mohammad Hoseini-Rad,
Akram Hasani Zonoozi,
Andreas H. W. Kuepper
Abstract:
We use direct $N$-body calculations to investigate the impact of primordial mass segregation on the size scale and mass-loss rate of star clusters in a galactic tidal field. We run a set of simulations of clusters with varying degrees of primordial mass segregation at various galactocentric radii and show that, in primordially segregated clusters, the early, impulsive mass-loss from stellar evolut…
▽ More
We use direct $N$-body calculations to investigate the impact of primordial mass segregation on the size scale and mass-loss rate of star clusters in a galactic tidal field. We run a set of simulations of clusters with varying degrees of primordial mass segregation at various galactocentric radii and show that, in primordially segregated clusters, the early, impulsive mass-loss from stellar evolution of the most massive stars in the innermost regions of the cluster leads to a stronger expansion than for initially non-segregated clusters. Therefore, models in stronger tidal fields dissolve faster due to an enhanced flux of stars over the tidal boundary. Throughout their lifetimes, the segregated clusters are more extended by a factor of about 2, suggesting that (at least) some of the very extended globular clusters in the outer halo of the Milky Way may have been born with primordial mass segregation. We finally derive a relation between star-cluster dissolution time, $T_{diss}$, and galactocentric radius, $R_G$, and show how it depends on the degree of primordial mass segregation.
△ Less
Submitted 4 September, 2014;
originally announced September 2014.
-
Mass segregation in the outer halo globular cluster Palomar 14
Authors:
Matthias J. Frank,
Eva K. Grebel,
Andreas H. W. Kuepper
Abstract:
We present evidence for mass segregation in the outer-halo globular cluster Palomar 14, which is intuitively unexpected since its present-day two-body relaxation time significantly exceeds the Hubble time. Based on archival Hubble Space Telescope imaging, we analyze the radial dependence of the stellar mass function in the cluster's inner 39.2 pc in the mass range of 0.53-0.80 M_sun, ranging from…
▽ More
We present evidence for mass segregation in the outer-halo globular cluster Palomar 14, which is intuitively unexpected since its present-day two-body relaxation time significantly exceeds the Hubble time. Based on archival Hubble Space Telescope imaging, we analyze the radial dependence of the stellar mass function in the cluster's inner 39.2 pc in the mass range of 0.53-0.80 M_sun, ranging from the main-sequence turn-off down to a V-band magnitude of 27.1 mag. The mass function at different radii is well approximated by a power law and rises from a shallow slope of 0.6+/-0.2 in the cluster's core to a slope of 1.6+/-0.3 beyond 18.6 pc. This is seemingly in conflict with the finding by Beccari et al. (2011), who interpret the cluster's non-segregated population of (more massive) blue straggler stars, compared to (less massive) red giants and horizontal branch stars, as evidence that the cluster has not experienced dynamical segregation yet. We discuss how both results can be reconciled. Our findings indicate that the cluster was either primordially mass-segregated and/or used to be significantly more compact in the past. For the latter case, we propose tidal shocks as the mechanism driving the cluster's expansion, which would imply that Palomar 14 is on a highly eccentric orbit. Conversely, if the cluster formed already extended and with primordial mass segregation, this could support an accretion origin of the cluster.
△ Less
Submitted 10 July, 2014;
originally announced July 2014.
-
Milky Way Mass and Potential Recovery Using Tidal Streams in a Realistic Halo
Authors:
Ana Bonaca,
Marla Geha,
Andreas H. W. Kuepper,
Juerg Diemand,
Kathryn V. Johnston,
David W. Hogg
Abstract:
We present a new method for determining the Galactic gravitational potential based on forward modeling of tidal stellar streams. We use this method to test the performance of smooth and static analytic potentials in representing realistic dark matter halos, which have substructure and are continually evolving by accretion. Our FAST-FORWARD method uses a Markov Chain Monte Carlo algorithm to compar…
▽ More
We present a new method for determining the Galactic gravitational potential based on forward modeling of tidal stellar streams. We use this method to test the performance of smooth and static analytic potentials in representing realistic dark matter halos, which have substructure and are continually evolving by accretion. Our FAST-FORWARD method uses a Markov Chain Monte Carlo algorithm to compare, in 6D phase space, an "observed" stream to models created in trial analytic potentials. We analyze a large sample of streams evolved in the Via Lactea II (VL2) simulation, which represents a realistic Galactic halo potential. The recovered potential parameters are in agreement with the best fit to the global, present-day VL2 potential. However, merely assuming an analytic potential limits the dark matter halo mass measurement to an accuracy of 5 to 20%, depending on the choice of analytic parametrization. Collectively, mass estimates using streams from our sample reach this fundamental limit, but individually they can be highly biased. Individual streams can both under- and overestimate the mass, and the bias is progressively worse for those with smaller perigalacticons, motivating the search for tidal streams at galactocentric distances larger than 70 kpc. We estimate that the assumption of a static and smooth dark matter potential in modeling of the GD-1 and Pal5-like streams introduces an error of up to 50% in the Milky Way mass estimates.
△ Less
Submitted 23 June, 2014;
originally announced June 2014.
-
Direct $N$-body simulations of globular clusters - II. Palomar 4
Authors:
Akram Hasani Zonoozi,
Hosein Haghi,
Andreas H. W. Kuepper,
Holger Baumgardt,
Matthias J. Frank,
Pavel Kroupa
Abstract:
We use direct $N$-body calculations to study the evolution of the unusually extended outer halo globular cluster Palomar 4 (Pal~4) over its entire lifetime in order to reproduce its observed mass, half-light radius, velocity dispersion and mass function slope at different radii.
We find that models evolving on circular orbits, and starting from a non-mass segregated, canonical initial mass funct…
▽ More
We use direct $N$-body calculations to study the evolution of the unusually extended outer halo globular cluster Palomar 4 (Pal~4) over its entire lifetime in order to reproduce its observed mass, half-light radius, velocity dispersion and mass function slope at different radii.
We find that models evolving on circular orbits, and starting from a non-mass segregated, canonical initial mass function (IMF) can reproduce neither Pal 4's overall mass function slope nor the observed amount of mass segregation. Including either primordial mass segregation or initially flattened IMFs does not reproduce the observed amount of mass segregation and mass function flattening simultaneously. Unresolved binaries cannot reconcile this discrepancy either. We find that only models with both a flattened IMF and primordial segregation are able to fit the observations. The initial (i.e. after gas expulsion) mass and half-mass radius of Pal~4 in this case are about 57000 M${\odot}$ and 10 pc, respectively. This configuration is more extended than most globular clusters we observe, showing that the conditions under which Pal~4 formed must have been significantly different from that of the majority of globular clusters. We discuss possible scenarios for such an unusual configuration of Pal~4 in its early years.
△ Less
Submitted 23 April, 2014;
originally announced April 2014.
-
Erosion of Globular Cluster Systems: The Influence of Radial Anisotropy, Central Black Holes and Dynamical Friction
Authors:
Michael Brockamp,
Andreas H. W. Kuepper,
Ingo Thies,
Holger Baumgardt,
Pavel Kroupa
Abstract:
We present the adaptable MUESLI code for investigating dynamics and erosion processes of globular clusters (GCs) in galaxies. MUESLI follows the orbits of individual clusters and applies internal and external dissolution processes to them. Orbit integration is based on the self-consistent field method in combination with a time-transformed leapfrog scheme, allowing us to handle velocity-dependent…
▽ More
We present the adaptable MUESLI code for investigating dynamics and erosion processes of globular clusters (GCs) in galaxies. MUESLI follows the orbits of individual clusters and applies internal and external dissolution processes to them. Orbit integration is based on the self-consistent field method in combination with a time-transformed leapfrog scheme, allowing us to handle velocity-dependent forces like triaxial dynamical friction. In a first application, the erosion of globular cluster systems (GCSs) in elliptical galaxies is investigated. Observations show that massive ellipticals have rich, radially extended GCSs, while some compact dwarf ellipticals contain no GCs at all. For several representative examples, spanning the full mass scale of observed elliptical galaxies, we quantify the influence of radial anisotropy, galactic density profiles, SMBHs, and dynamical friction on the GC erosion rate. We find that GC number density profiles are centrally flattened in less than a Hubble time, naturally explaining observed cored GC distributions. The erosion rate depends primarily on a galaxy's mass, half-mass radius and radial anisotropy. The fraction of eroded GCs is nearly 100% in compact, M 32 like galaxies and lowest in extended and massive galaxies. Finally, we uncover the existence of a violent tidal disruption dominated phase which is important for the rapid build-up of halo stars.
△ Less
Submitted 19 March, 2014;
originally announced March 2014.
-
Mass segregation in the diffuse outer-halo globular cluster Palomar 14
Authors:
Matthias J. Frank,
Eva K. Grebel,
Andreas H. W. Küpper
Abstract:
We present an analysis of the radial dependence of the stellar mass function in the diffuse outer-halo globular cluster Palomar 14. Using archival HST/WFPC2 data of the cluster's central 39 pc (corresponding to ~0.85*r_h) we find that the mass function in the mass range of 0.55 to 0.85 solar masses is well approximated by a power-law at all radii. The mass function steepens with increasing radius,…
▽ More
We present an analysis of the radial dependence of the stellar mass function in the diffuse outer-halo globular cluster Palomar 14. Using archival HST/WFPC2 data of the cluster's central 39 pc (corresponding to ~0.85*r_h) we find that the mass function in the mass range of 0.55 to 0.85 solar masses is well approximated by a power-law at all radii. The mass function steepens with increasing radius, from a shallow power-law slope of 0.66+/-0.32 in the cluster's centre to a slope of 1.61+/-0.33 beyond the core radius, showing that the cluster is mass-segregated. This is seemingly in conflict with its long present-day half-mass relaxation time of ~20 Gyr, and with the recent finding by Beccari et al. (2011), who interpret the cluster's non-concentrated population of blue straggler stars as evidence that dynamical segregation has not affected the cluster yet. We discuss this apparent conflict and argue that the cluster must have either formed with primordial mass segregation, or that its relaxation time scale must have been much smaller in the past, i.e. that the cluster must have undergone a significant expansion.
△ Less
Submitted 21 January, 2013;
originally announced January 2013.
-
The velocity dispersion and mass function of the outer halo globular cluster Palomar 4
Authors:
Matthias J. Frank,
Michael Hilker,
Holger Baumgardt,
Patrick Côté,
Eva K. Grebel,
Hosein Haghi,
Andreas H. W. Küpper,
S. G. Djorgovski
Abstract:
We obtained precise line-of-sight radial velocities of 23 member stars of the remote halo globular cluster Palomar 4 (Pal 4) using the High Resolution Echelle Spectrograph (HIRES) at the Keck I telescope. We also measured the mass function of the cluster down to a limiting magnitude of V~28 mag using archival HST/WFPC2 imaging. We derived the cluster's surface brightness profile based on the WFPC2…
▽ More
We obtained precise line-of-sight radial velocities of 23 member stars of the remote halo globular cluster Palomar 4 (Pal 4) using the High Resolution Echelle Spectrograph (HIRES) at the Keck I telescope. We also measured the mass function of the cluster down to a limiting magnitude of V~28 mag using archival HST/WFPC2 imaging. We derived the cluster's surface brightness profile based on the WFPC2 data and on broad-band imaging with the Low-Resolution Imaging Spectrometer (LRIS) at the Keck II telescope. We find a mean cluster velocity of 72.55+/-0.22 km/s and a velocity dispersion of 0.87+/-0.18 km/s. The global mass function of the cluster, in the mass range 0.55<=M<=0.85 M_solar, is shallower than a Kroupa mass function and the cluster is significantly depleted in low-mass stars in its center compared to its outskirts. Since the relaxation time of Pal 4 is of the order of a Hubble time, this points to primordial mass segregation in this cluster. Extrapolating the measured mass function towards lower-mass stars and including the contribution of compact remnants, we derive a total cluster mass of 29800 M_solar. For this mass, the measured velocity dispersion is consistent with the expectations of Newtonian dynamics and below the prediction of Modified Newtonian Dynamics (MOND). Pal 4 adds to the growing body of evidence that the dynamics of star clusters in the outer Galactic halo can hardly be explained by MOND.
△ Less
Submitted 29 November, 2016; v1 submitted 11 May, 2012;
originally announced May 2012.
-
The Tidal Tails of 47 Tucanae
Authors:
Richard R. Lane,
Andreas H. W. Küpper,
Douglas C. Heggie
Abstract:
The Galactic globular cluster 47 Tucanae (47 Tuc) shows a rare increase in its velocity dispersion profile at large radii, indicative of energetic, yet bound, stars at large radii dominating the velocity dispersion and, potentially, of ongoing evaporation. Escaping stars will form tidal tails, as seen with several Galactic globular clusters, however, the tidal tails of 47 Tuc are yet to be uncover…
▽ More
The Galactic globular cluster 47 Tucanae (47 Tuc) shows a rare increase in its velocity dispersion profile at large radii, indicative of energetic, yet bound, stars at large radii dominating the velocity dispersion and, potentially, of ongoing evaporation. Escaping stars will form tidal tails, as seen with several Galactic globular clusters, however, the tidal tails of 47 Tuc are yet to be uncovered. We model these tails of 47 Tuc using the most accurate input data available, with the specific aim of determining their locations, as well as the densities of the epicyclic overdensities within the tails. The overdensities from our models show an increase of 3-4% above the Galactic background and, therefore, should be easily detectable using matched filtering techniques. We find that the most influential parameter with regard to both the locations and densities of the epicyclic overdensities is the Heliocentric distance to the cluster. Hence, uncovering these tidal features observationally will contribute greatly to the ongoing problem of determining the distance to 47 Tuc, tightly constraining the distance of the cluster independent of other methods. Using our streakline method for determining the locations of the tidal tails and their overdensities, we show how, in principle, the shape and extent of the tidal tails of any Galactic globular cluster can be determined without resorting to computationally expensive N-body simulations.
△ Less
Submitted 3 July, 2012; v1 submitted 11 April, 2012;
originally announced April 2012.
-
More on the structure of tidal tails
Authors:
Andreas H. W. Kuepper,
Richard R. Lane,
Douglas C. Heggie
Abstract:
We investigate the epicyclic motion of stars escaping from star clusters. Using streaklines, we visualise the path of escaping stars and show how epicyclic motion leads to over- and underdensities in tidal tails of star clusters moving on circular and eccentric orbits about a galaxy. Additionally, we investigate the effect of the cluster mass on the tidal tails, by showing that their structure is…
▽ More
We investigate the epicyclic motion of stars escaping from star clusters. Using streaklines, we visualise the path of escaping stars and show how epicyclic motion leads to over- and underdensities in tidal tails of star clusters moving on circular and eccentric orbits about a galaxy. Additionally, we investigate the effect of the cluster mass on the tidal tails, by showing that their structure is better matched when the perturbing effect of the cluster mass is included. By adjusting streaklines to results of N-body computations we can accurately and quickly reproduce all observed substructure, especially the streaky features often found in simulations which may be interpreted in observations as multiple tidal tails. Hence, we can rule out tidal shocks as the origin of such substructures. Finally, from the adjusted streakline parameters we can verify that for the star clusters we studied escape mainly happens from the tidal radius of the cluster, given by x_L = (GM/(Ω^2-\partial^2Φ/\partial R^2))^{1/3}. We find, however, that there is another limiting radius, the "edge" radius, which gives the smallest radius from which a star can escape during one cluster orbit about the galaxy. For eccentric cluster orbits the edge radius shrinks with increasing orbital eccentricity (for fixed apocentric distance) but is always significantly larger than the respective perigalactic tidal radius. In fact, the edge radii of the clusters we investigated, which are extended and tidally filling, agree well with their (fitted) King radii, which may indicate a fundamental connection between these two quantities.
△ Less
Submitted 15 December, 2011; v1 submitted 21 November, 2011;
originally announced November 2011.
-
Mass segregation and fractal substructure in young massive clusters: (I) the McLuster code and method calibration
Authors:
A. H. W. Kuepper,
Th. Maschberger,
P. Kroupa,
H. Baumgardt
Abstract:
By analysing models of the young massive cluster R136 in 30 Doradus, set-up using the herewith introduced and publicly made available code McLuster, we investigate and compare different methods for detecting and quantifying mass segregation and substructure in non-seeing limited N-body data. For this purpose we generate star cluster models with different degrees of mass segregation and fractal sub…
▽ More
By analysing models of the young massive cluster R136 in 30 Doradus, set-up using the herewith introduced and publicly made available code McLuster, we investigate and compare different methods for detecting and quantifying mass segregation and substructure in non-seeing limited N-body data. For this purpose we generate star cluster models with different degrees of mass segregation and fractal substructure and analyse them. We quantify mass segregation by measuring, from the projected 2d model data, the mass function slope in radial annuli, by looking for colour gradients in radial colour profiles, by measuring Allison's Lambda parameter, and by determining the local stellar surface density around each star. We find that these methods for quantifying mass segregation often produce ambiguous results. Most reliable for detecting mass segregation is the mass function slope method, whereas the colour gradient method is the least practical in an R136-like configuration. The other two methods are more sensitive to low degrees of mass segregation but are computationally much more demanding. We also discuss the effect of binaries on these measures. Moreover, we quantify substructure by looking at the projected radial stellar density profile, by comparing projected azimuthal stellar density profiles, and by determining Cartwright & Whitworth's Q parameter. We find that only high degrees of substructure affect the projected radial density profile, whereas the projected azimuthal density profile is very sensitive to substructure. The Q parameter is also sensitive to substructure but its absolute value shows a dependence on the radial density gradient of the cluster and is strongly influenced by binaries. (abridged)
△ Less
Submitted 16 August, 2011; v1 submitted 12 July, 2011;
originally announced July 2011.
-
The Curious Case of Palomar 13: The Influence of the Orbital Phase on the Appearance of Galactic Satellites
Authors:
A. H. W. Kuepper,
S. Mieske,
P. Kroupa
Abstract:
We investigate the dynamical status of the low-mass globular cluster Palomar 13 by means of N-body computations to test whether its unusually high mass-to-light ratio of about 40 and its peculiarly shallow surface density profile can be caused by tidal shocking. Alternatively, we test - by varying the assumed proper motion - if the orbital phase of Palomar 13 within its orbit about the Milky Way c…
▽ More
We investigate the dynamical status of the low-mass globular cluster Palomar 13 by means of N-body computations to test whether its unusually high mass-to-light ratio of about 40 and its peculiarly shallow surface density profile can be caused by tidal shocking. Alternatively, we test - by varying the assumed proper motion - if the orbital phase of Palomar 13 within its orbit about the Milky Way can influence its appearance and thus may be the origin of these peculiarities, as has been suggested by Kuepper et al. (2010). We find that, of these two scenarios, only the latter can explain the observed mass-to-light ratio and surface density profile. We note, however, that the particular orbit that best reproduces those observed parameters has a proper motion inconsistent with the available literature value. We discuss this discrepancy and suggest that it may be caused by an underestimation of the observational uncertainties in the proper motion determination. We demonstrate that Palomar 13 is most likely near apogalacticon, which makes the cluster appear supervirial and blown-up due to orbital compression of its tidal debris. Since the satellites of the Milky Way are on average closer to apo- than perigalacticon, their internal dynamics may be influenced by the same effect, and we advocate that this needs to be taken into account when interpreting their kinematical data. Moreover, we briefly discuss the influence of a possible binary population on such measurements.
△ Less
Submitted 10 June, 2013; v1 submitted 14 December, 2010;
originally announced December 2010.
-
Direct N-body simulations of globular clusters: (I) Palomar 14
Authors:
Akram Hasani Zonoozi,
Andreas H. W. Kuepper,
Holger Baumgardt,
Hosein Haghi,
Pavel Kroupa,
Michael Hilker
Abstract:
We present the first ever direct $N$-body computations of an old Milky Way globular cluster over its entire life time on a star-by-star basis. Using recent GPU hardware at Bonn University, we have performed a comprehensive set of $N$-body calculations to model the distant outer halo globular cluster Palomar 14 (Pal 14). By varying the initial conditions we aim at finding an initial $N$-body model…
▽ More
We present the first ever direct $N$-body computations of an old Milky Way globular cluster over its entire life time on a star-by-star basis. Using recent GPU hardware at Bonn University, we have performed a comprehensive set of $N$-body calculations to model the distant outer halo globular cluster Palomar 14 (Pal 14). By varying the initial conditions we aim at finding an initial $N$-body model which reproduces the observational data best in terms of its basic parameters, i.e. half-light radius, mass and velocity dispersion. We furthermore focus on reproducing the stellar mass function slope of Pal 14 which was found to be significantly shallower than in most globular clusters. While some of our models can reproduce Pal 14's basic parameters reasonably well, we find that dynamical mass segregation alone cannot explain the mass function slope of Pal 14 when starting from the canonical Kroupa initial mass function (IMF). In order to seek for an explanation for this discrepancy, we compute additional initial models with varying degrees of primordial mass segregation as well as with a flattened IMF. The necessary degree of primordial mass segregation turns out to be very high. This modelling has shown that the initial conditions of Pal 14 after gas expulsion must have been a half-mass radius of about 20 pc, a mass of about 50000 M$_{\odot}$, and possibly some mass segregation or an already established non-canonical IMF depleted in low-mass stars. Such conditions might be obtained by a violent early gas-expulsion phase from an embedded cluster born with mass segregation. Only at large Galactocentric radii are clusters likely to survive as bound entities the destructive gas-expulsion process we seem to have uncovered for Pal 14. In addition we compute a model with a 5% primordial binary fraction to test if such a population has an effect on the cluster's evolution.
△ Less
Submitted 11 October, 2010;
originally announced October 2010.
-
Peculiarities in Velocity Dispersion and Surface Density Profiles of Star Clusters
Authors:
A. H. W. Kuepper,
P. Kroupa,
H. Baumgardt,
D. C. Heggie
Abstract:
Based on our recent work on tidal tails of star clusters (Kuepper et al. 2009) we investigate star clusters of a few 10^4 Msun by means of velocity dispersion profiles and surface density profiles. We use a comprehensive set of $N$-body computations of star clusters on various orbits within a realistic tidal field to study the evolution of these profiles with time, and ongoing cluster dissolution…
▽ More
Based on our recent work on tidal tails of star clusters (Kuepper et al. 2009) we investigate star clusters of a few 10^4 Msun by means of velocity dispersion profiles and surface density profiles. We use a comprehensive set of $N$-body computations of star clusters on various orbits within a realistic tidal field to study the evolution of these profiles with time, and ongoing cluster dissolution From the velocity dispersion profiles we find that the population of potential escapers, i.e. energetically unbound stars inside the Jacobi radius, dominates clusters at radii above about 50% of the Jacobi radius. Beyond 70% of the Jacobi radius nearly all stars are energetically unbound. The velocity dispersion therefore significantly deviates from the predictions of simple equilibrium models in this regime. We furthermore argue that for this reason this part of a cluster cannot be used to detect a dark matter halo or deviations from Newtonian gravity. By fitting templates to the about 10^4 computed surface density profiles we estimate the accuracy which can be achieved in reconstructing the Jacobi radius of a cluster in this way. We find that the template of King (1962) works well for extended clusters on nearly circular orbits, but shows significant flaws in the case of eccentric cluster orbits. This we fix by extending this template with 3 more free parameters. Our template can reconstruct the tidal radius over all fitted ranges with an accuracy of about 10%, and is especially useful in the case of cluster data with a wide radial coverage and for clusters showing significant extra-tidal stellar populations. No other template that we have tried can yield comparable results over this range of cluster conditions. All templates fail to reconstruct tidal parameters of concentrated clusters, however. (abridged)
△ Less
Submitted 25 May, 2010;
originally announced May 2010.
-
Is the distant globular cluster Pal 14 in a deep-freeze?
Authors:
Andreas H. W. Kuepper,
Pavel Kroupa
Abstract:
We investigate the velocity dispersion of Pal 14, an outer Milky-Way globular cluster at Galactocentric distance of 71 kpc with a very low stellar density (central density 0.1-0.2 Msun/pc^3). Due to this low stellar density the binary population of Pal 14 is likely to be close to the primordial binary population. Artificial clusters are generated with the observed properties of Pal 14 and the velo…
▽ More
We investigate the velocity dispersion of Pal 14, an outer Milky-Way globular cluster at Galactocentric distance of 71 kpc with a very low stellar density (central density 0.1-0.2 Msun/pc^3). Due to this low stellar density the binary population of Pal 14 is likely to be close to the primordial binary population. Artificial clusters are generated with the observed properties of Pal 14 and the velocity dispersion within these clusters is measured as Jordi et al. (2009) have done with 17 observed stars of Pal 14. We discuss the effect of the binary population on these measurements and find that the small velocity dispersion of 0.38 km/s which has been found by Jordi et al. (2009) would imply a binary fraction of less than 0.1, even though from the stellar density of Pal 14 we would expect a binary fraction of more than 0.5. We also discuss the effect of mass segregation on the velocity dispersion as possible explanation for this discrepancy, but find that it would increase the velocity dispersion further. Thus, either Pal 14 has a very unusual stellar population and its birth process was significantly different than we see in today's star forming regions, or the binary population is regular and we would have to correct the observed 0.38 km/s for binarity. In this case the true velocity dispersion of Pal 14 would be much smaller than this value and the cluster would have to be considered as "kinematically frigid", thereby possibly posing a challenge for Newtonian dynamics but in the opposite sense to MOND.
△ Less
Submitted 3 May, 2010;
originally announced May 2010.
-
Tidal Tails of Star Clusters
Authors:
A. H. W. Kuepper,
P. Kroupa,
H. Baumgardt,
D. C. Heggie
Abstract:
Based on recent findings of a formation mechanism of substructure in tidal tails by Kuepper, Macleod & Heggie (2008) we investigate a more comprehensive set of N-body models of star clusters on orbits about a Milky-Way-like potential. We find that the predicted epicyclic overdensities arise in any tidal tail no matter which orbit the cluster follows as long as the cluster lives long enough for t…
▽ More
Based on recent findings of a formation mechanism of substructure in tidal tails by Kuepper, Macleod & Heggie (2008) we investigate a more comprehensive set of N-body models of star clusters on orbits about a Milky-Way-like potential. We find that the predicted epicyclic overdensities arise in any tidal tail no matter which orbit the cluster follows as long as the cluster lives long enough for the overdensities to build up.
The distance of the overdensities along the tidal tail from the cluster centre depends for circular orbits only on the mass of the cluster and the strength of the tidal field, and therefore decreases monotonically with time, while for eccentric orbits the orbital motion influences the distance, causing a periodic compression and stretching of the tails and making the distance oscillate with time. We provide an approximation for estimating the distance of the overdensities in this case.
We describe an additional type of overdensity which arises in extended tidal tails of clusters on eccentric orbits, when the acceleration of the tidal field on the stellar stream is no longer homogeneous. Moreover, we conclude that a pericentre passage or a disk shock is not the direct origin of an overdensity within a tidal tail. Escape due to such tidal perturbations does not take place immediately after the perturbation but is rather delayed and spread over the orbit of the cluster. All observable overdensities are therefore of the mentioned two types. In particular, we note that substructured tidal tails do not imply the existence of dark-matter sub-structures in the haloes of galaxies.
△ Less
Submitted 14 October, 2009; v1 submitted 14 September, 2009;
originally announced September 2009.
-
The Main Sequence of Star Clusters
Authors:
Andreas H. W. Kuepper,
Pavel Kroupa,
Holger Baumgardt
Abstract:
A novel way of looking at the evolution of star clusters is presented. With a dynamical temperature, given by the mean kinetic energy of the cluster stars, and a dynamical luminosity, which is defined as the kinetic energy of the stars leaving the cluster in analogy to the energy of photons emitted by a star, the dissolution of star clusters is studied using a new dynamical temperature-luminosit…
▽ More
A novel way of looking at the evolution of star clusters is presented. With a dynamical temperature, given by the mean kinetic energy of the cluster stars, and a dynamical luminosity, which is defined as the kinetic energy of the stars leaving the cluster in analogy to the energy of photons emitted by a star, the dissolution of star clusters is studied using a new dynamical temperature-luminosity diagram for star clusters. The investigation contains a parameter-space study of open clusters of up to N = 32768 single-mass stars with different initial density distributions, half-mass radii, tidal conditions and binary fractions. The clusters show a strong correlation between dynamical temperature and dynamical luminosity and most of the investigated cluster families share a common sequence in such a dynamical temperature-luminosity diagram. Deviations from this sequence are analyzed and discussed. After core collapse, the position of a cluster within this diagram can be defined by three parameters: the mass, the tidal conditions and the binary fraction. Due to core collapse all initial conditions are lost and the remaining stars adjust to the given tidal conditions. Binaries as internal energy sources influence this adjustment. A further finding concerns the Lagrange radii of star clusters: Throughout the investigated parameter space nearly all clusters show a constant half-mass radius for the time after core collapse until dissolution. Furthermore, the ratio of half-mass radius to tidal radius evolves onto a common sequence which only depends on the mass left in the cluster.
△ Less
Submitted 24 June, 2008;
originally announced June 2008.
-
On the structure of tidal tails
Authors:
Andreas H. W. Küpper,
Andrew Macleod,
Douglas C. Heggie
Abstract:
We examine the longitudinal distribution of the stars escaping from a cluster along tidal tails. Using both theory and simulations, we show that, even in the case of a star cluster in a circular galactic orbit, when the tide is steady, the distribution exhibits maxima at a distance of many tidal radii from the cluster.
We examine the longitudinal distribution of the stars escaping from a cluster along tidal tails. Using both theory and simulations, we show that, even in the case of a star cluster in a circular galactic orbit, when the tide is steady, the distribution exhibits maxima at a distance of many tidal radii from the cluster.
△ Less
Submitted 15 April, 2008;
originally announced April 2008.