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 {\Lambda}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.

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, $\sigma$, of their host galaxy. In the simplest regime it is $M_\bullet \sim \sigma^{3/2}\sqrt{M_\star t / G} \sim 10^{6}M_\odot (\sigma / 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.

We present $N$-body simulations of globular clusters, exploring the effect of different galactic potentials on cluster sizes, $r_h$. For various galactocentric distances, $R_G$, we assess how cluster sizes change when we vary the virial mass and concentration of the host galaxy's dark-matter halo. We show that sizes of GCs are determined by the local galactic mass density rather than the virial mass of the host galaxy. We find that clusters evolving in the inner haloes of less concentrated galaxies are significantly more extended than those evolving in more concentrated ones, while the sizes of those orbiting in the outer halo are almost independent of concentration. Adding a baryonic component to our galaxy models does not change these results much, since its effect is only significant in the very inner halo. Our simulations suggest that there is a relation between $r_h$ and $R_G$, which systematically depends on the physical parameters of the halo. Hence, observing such relations in individual galaxies can put a new observational constraint on dark-matter halo characteristics. However, by varying the halo mass in a wide range of $10^{9}\leq M_{vir}/\msun \leq10^{13}$, we find that the $r_h-R_G$ relationship will be nearly independent of halo mass, if one assumes $M_{vir}$ and $c_{vir}$ as two correlated parameters, as is suggested by cosmological simulations.

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.

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.

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.

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.

We quantify to what extent tidal erosion of globular clusters (GCs) has contributed to the observed u-shaped relation between GC specific frequencies S_N and host galaxy luminosity M_V. We used our MUESLI code to calculate GC survival rates for typical early-type galaxy potentials covering a wide range of observed galaxy properties. We do this for isotropic and radially anisotropic GC velocity distributions. We find that the calculated GC survival fraction, f_s, depends linearly on the logarithm of the 3D mass density, rho_3D, within the galaxy's half light radius, with f_s proportional to (rho_3D)^(-0.17). For a given galaxy, survival rates are lower for radially anisotropic configurations than for the isotropic GC cases. We apply these relations to a literature sample of 219 early-type galaxies from Harris et al. (2013) in the range M_V=[-24.5:-15.5] mag. The expected GC survival fraction ranges from ~50% for the most massive galaxies with the largest radii to ~10% for the most compact galaxies. We find that intermediate luminosity galaxies M_V=[-20.5:-17.5] mag have the strongest expected GC erosion. Within the considered literature sample, the predicted GC survival fraction therefore defines a u-shaped relation with M_V, similar to the relation between specific frequency S_N and M_V. As a consequence, the u-shape of S_N vs. M_V gets erased almost entirely when correcting the S_N values for the effect of GC erosion. We conclude that tidal erosion is an important contributor to the u-shaped relation between GC specific frequency and host galaxy luminosity. It must be taken into account when inferring primordial star cluster formation efficiencies from observations of GC systems in the nearby universe.

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.

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.

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.

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/(\Omega^2-\partial^2\Phi/\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.

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)

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.

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)

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.

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.

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.