Source author record

Holger Baumgardt

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

51works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

51 published item(s)

preprint2021arXiv

First direct dynamical detection of a dual supermassive black hole system at sub-kiloparsec separation

We investigated whether the two recently discovered nuclei in NGC7727 both host a super-massive black hole (SMBH). We used the high spatial resolution mode of the integral-field spectrograph MUSE on the VLT in adaptive optics mode to resolve the stellar kinematics within the sphere of influence of both putative black holes. We combined the kinematic data with an HST-based mass model and used Jeans models to measure their SMBH mass. We report the discovery of a dual SMBH system in NGC7727. We detect a SMBH in the photometric center of the galaxy in Nucleus 1, with a mass of $M_{\rm SMBH}=1.54^{+0.18}_{-0.15} \times10^{8}M_{\odot}$. In the second nucleus, which is 500pc offset from the main nucleus, we also find a clear signal for a SMBH with a mass of $M_{\rm BH}=6.33^{+3.32}_{-1.40}\times10^{6}M_{\odot}$. Both SMBHs are detected at high significance. The off-axis nature of Nucleus 2 makes modeling the system challenging; but a number of robustness tests suggest that a black hole is required to explain the kinematics. The SMBH in the offset Nucleus 2 makes up 3.0% of its total mass, which means its SMBH is over-massive compared to the $M_{\rm BH}-M_{\rm Bulge}$ scaling relation. This confirms it as the surviving nuclear star cluster of a galaxy that has merged with NGC7727. This discovery is the first dynamically confirmed dual SMBH system with a projected separation of less than a kiloparsec and the nearest dynamically confirmed dual SMBH at a distance of 27.4Mpc. The second Nucleus is in an advanced state of inspiral, and it will eventually result in a 1:24 mass ratio SMBH merger. Optical emission lines suggest Nucleus 2 is a Seyfert galaxy, making it a low-luminosity Active Galactic Nuclei (AGN). There are likely many more quiescent SMBHs as well as dual SMBH pairs in the local Universe that have been missed by surveys that focus on bright accretion signatures.

preprint2021arXiv

The dynamics of the globular cluster NGC3201 out to the Jacobi radius

As part of a chemo-dynamical survey of five nearby globular clusters with 2dF/AAOmega on the AAT, we have obtained kinematic information for the globular cluster NGC3201. Our new observations confirm the presence of a significant velocity gradient across the cluster which can almost entirely be explained by the high proper motion of the cluster. After subtracting the contribution of this perspective rotation, we found a remaining rotation signal with an amplitude of $\sim1\ km/s$ around a different axis to what we expect from the tidal tails and the potential escapers, suggesting that this rotation is internal and can be a remnant of its formation process. At the outer part, we found a rotational signal that is likely a result from potential escapers. The proper motion dispersion at large radii reported by Bianchini et al. has previously been attributed to dark matter. Here we show that the LOS dispersion between 0.5-1 Jacobi radius is lower, yet above the predictions from an N-body model of NGC3201 that we ran for this study. Based on the simulation, we find that potential escapers cannot fully explain the observed velocity dispersion. We also estimate the effect on the velocity dispersion of different amounts of stellar-mass black holes and unbound stars from the tidal tails with varying escape rates and find that these effects cannot explain the difference between the LOS dispersion and the N-body model. Given the recent discovery of tidal tail stars at large distances from the cluster, a dark matter halo is an unlikely explanation. We show that the effect of binary stars, which is not included in the N-body model, is important and can explain part of the difference in dispersion. We speculate that the remaining difference must be the result of effects not included in the N-body model, such as initial cluster rotation, velocity anisotropy and Galactic substructure.

preprint2019arXiv

Galactic globular clusters: a new catalog of masses, structural parameters, velocity dispersion profiles, proper motions and space orbits

We collected radial velocities of more than 50.000 individual stars in 156 Galactic globular clusters (GGC) and matched them with HST photometry and Gaia DR2 proper motions. This allowed us to derive the GGC's mean proper motions and space velocities. By fitting a large set of N-body simulations to their velocity dispersion and surface density profiles, combined with new measurements of their internal radially dependent mass functions, we have determined their present-day masses and structural parameters, and for 144 GGCs their internal kinematics. We also derive the initial cluster masses by calculating the cluster orbits backwards in time applying suitable recipes to account for mass-loss and dynamical friction. The new fundamental parameters of GGCs are publicly available via an online database, which will regularly be updated.

preprint2016arXiv

Constraining ultra-compact dwarf galaxy formation with galaxy clusters in the local universe

We compare the predictions of a semi-analytic model for ultra-compact dwarf galaxy (UCD) formation by tidal stripping to the observed properties of globular clusters (GCs) and UCDs in the Fornax and Virgo clusters. For Fornax we find the predicted number of stripped nuclei agrees very well with the excess number of GCs$+$UCDs above the GC luminosity function. GCs$+$UCDs with masses $>10^{7.3}$ M$_\odot$ are consistent with being entirely formed by tidal stripping. Stripped nuclei can also account for Virgo UCDs with masses $>10^{7.3}$ M$_\odot$ where numbers are complete by mass. For both Fornax and Virgo, the predicted velocity dispersions and radial distributions of stripped nuclei are consistent with that of UCDs within $\sim$50-100 kpc but disagree at larger distances where dispersions are too high and radial distributions too extended. Stripped nuclei are predicted to have radially biased anisotropies at all radii, agreeing with Virgo UCDs at clustercentric distances larger than 50 kpc. However, ongoing disruption is not included in our model which would cause orbits to become tangentially biased at small radii. We find the predicted metallicities and central black hole masses of stripped nuclei agree well with the metallicities and implied black hole masses of UCDs for masses $>10^{6.5}$ M$_\odot$. The predicted black hole masses also agree well with that of M60-UCD1, the first UCD with a confirmed central black hole. These results suggest that observed GC$+$UCD populations are a combination of genuine GCs and stripped nuclei, with the contribution of stripped nuclei increasing toward the high-mass end.

preprint2016arXiv

Limits on the significant mass-loss scenario based on the globular clusters of the Fornax dwarf spheroidal galaxy

Many of the scenarios proposed to explain the origin of chemically peculiar stars in globular clusters (GCs) require significant mass-loss ($\ge95\%$) to explain the observed fraction of such stars. In the GCs of the Fornax dwarf galaxy significant mass-loss could be a problem. Larsen et al. (2012) showed that there is a large ratio of GCs to metal-poor field stars in Fornax and about $20-25\%$ of all the stars with ${\rm [Fe/H]}<-2$ belong to the four metal-poor GCs. This imposes an upper limit of $\sim80\%$ mass-loss that could have happened in Fornax GCs. In this paper, we propose a solution to this problem by suggesting that stars can leave the Fornax galaxy. We use a series of $N$-body simulations, to determine the limit of mass-loss from Fornax as a function of the initial orbital radii of GCs and the speed with which stars leave Fornax GCs. We consider a set of cored and cuspy density profiles for Fornax. Our results show that with a cuspy model for Fornax, the fraction of stars which leave the galaxy, can be as high as $\sim90\%$, when the initial orbital radii of GCs are $R=2-3\,\rm{kpc}$ and the initial speed of stars is $v>20\,\rm{km\,s^{-1}}$. We show that such large velocities can be achieved by gas expulsion induced mass-loss but not stellar evolution induced mass-loss. Our results imply that one cannot interpret the metallicity distribution of Fornax field stars as evidence against significant mass-loss in Fornax GCs, if mass-loss is due to gas expulsion.

preprint2016arXiv

N-body modeling of globular clusters: Masses, mass-to-light ratios and intermediate-mass black holes

We have determined the masses and mass-to-light ratios of 50 Galactic globular clusters by comparing their velocity dispersion and surface brightness profiles against a large grid of 900 N-body simulations of star clusters of varying initial concentration, size and central black hole mass fraction. Our models follow the evolution of the clusters under the combined effects of stellar evolution and two-body relaxation allowing us to take the effects of mass segregation and energy equipartition between stars self-consistently into account. For a subset of 16 well observed clusters we also derive their kinematic distances. We find an average mass-to-light ratio of Galactic globular clusters of $<M/L_V>=1.98 \pm 0.03$, which agrees very well with the expected M/L ratio if the initial mass function of the clusters was a standard Kroupa or Chabrier mass function. We do not find evidence for a decrease of the average mass-to-light ratio with metallicity. The surface brightness and velocity dispersion profiles of most globular clusters are incompatible with the presence of intermediate-mass black holes (IMBHs) with more than a few thousand $M_\odot$ in them. The only clear exception is $ω$ Cen, where the velocity dispersion profile provides strong evidence for the presence of a $\sim$40,000 $M_\odot$ IMBH in the centre of the cluster.

preprint2016arXiv

Probing the boundary between star clusters and dwarf galaxies: A MUSE view on the dynamics of Crater/Laevens I

We present MUSE observations of the debated ultra faint stellar system Crater. We spectroscopically confirm 26 member stars of this system via radial velocity measurements. We derive the systematic instrumental velocity uncertainty of MUSE spectra to be 2.27$\rm \,km\,s^{-1}$. This new dataset increases the confirmed member stars of Crater by a factor of 3. One out of three bright blue stars and a fainter blue star just above the main-sequence-turn-off are also found to be likely members of the system. The observations reveal that Crater has a systemic radial velocity of $v_{\rm sys}=148.18^{\rm +1.08}_{\rm -1.15}\rm \,km\,s^{-1}$, whereas the most likely velocity dispersion of this system is $σ_{\rm v}=2.04^{\rm +2.19}_{\rm -1.06} \rm \,km\,s^{-1}$. The total dynamical mass of the system, assuming dynamical equilibrium is then $M_{\rm tot}=1.50^{+4.9}_{-1.2}\cdot 10^{\rm 5}M_{\odot}$ implying a mass-to-light ratio of M/L$_{\rm V}$=8.52$^{+28.0}_{-6.5}\, M_{\odot}/L_{\odot}$, which is consistent with a purely baryonic stellar population within its errors and no significant evidence for the presence dark matter was found. We also find evidence for a velocity gradient in the radial velocity distribution. We conclude that our findings strongly support that Crater is a faint intermediate-age outer halo globular cluster and not a dwarf galaxy.

preprint2016arXiv

The Binary Fraction and Mass Segregation in Alpha Persei Open Cluster

We have obtained membership probabilities of stars within a field of radius $\sim3^\circ$ around the centre of the open cluster Alpha Persei using proper motions and photometry from the PPMXL and WISE catalogues. We have identified 810 possible stellar members of Alpha Persei. We derived the global and radial present-day mass function (MF) of the cluster and found that they are well matched by two-stage power-law relations with different slopes at different radii. The global MF of Alpha Persei shows a turnover at $m=0.62\,\mathrm{M}_{\odot}$ with low and high-mass slopes of $α_\mathrm{low}=0.50\pm0.09$ ($0.1<m/\mathrm{M}_{\odot}<0.62$) and $α_\mathrm{high}=2.32\pm0.14$ ($0.62\leq m/\mathrm{M}_{\odot}<4.68$) respectively. The high-mass slope of the cluster increases from $2.01$ inside $1\hbox{$.\!\!^\circ$}10$ to $2.63$ outside $2\hbox{$.\!\!^\circ$}2$, whereas the mean stellar mass decreases from $0.95$ to $0.57\,\mathrm{M}_{\odot}$ in the same regions, signifying clear evidence of mass segregation in the cluster. From an examination of the high-quality colour-magnitude data of the cluster and performing a series of Monte Carlo simulations we obtained a binary fraction of $f_{\rm bin}=34\pm12$ percent for stars with $0.70<m/\mathrm{M}_{\odot}<4.68$. This is significantly larger than the observed binary fraction, indicating that this open cluster contains a large population of unresolved binaries. Finally, we corrected the mass-function slopes for the effect of unresolved binaries and found low- and high-mass slopes of $α_\mathrm{low}=0.89\pm0.11$ and $α_\mathrm{high}=2.37\pm0.09$ and a total cluster mass of $352\,\mathrm{M}_{\odot}$.

preprint2016arXiv

The velocity dispersion and mass function of the outer halo globular cluster Palomar 4

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.

preprint2015arXiv

A re-evaluation of the central velocity-dispersion profile in NGC 6388

Recently, two independent groups found very different results when measuring the central velocity dispersion of the galactic globular cluster NGC 6388 with different methods. While Lützgendorf et al. (2011) found a rising profile and a high central velocity dispersion (23.3 km/s), measurements obtained by Lanzoni et al. (2013) showed a value 40% lower. The value of the central velocity dispersion has a serious impact on the mass and possible presence of an intermediate-mass black hole at the center of NGC 6388. We use a photometric catalog of NGC 6388 to create a simulated SINFONI and ARGUS dataset. The construction of the IFU data cube is done with different observing conditions reproducing the conditions reported for the original observations as closely as possible. In addition, we produce an N-body realization of a 10^6 M_SUN stellar cluster with the same photometric properties as NGC 6388 to account for unresolved stars. We find that the individual radial velocities, i.e. the measurements from the simulated SINFONI data, are systematically biased towards lower velocity dispersions. The reason is that due to the wings in the point spread function the velocities get biased towards the mean cluster velocity. This study shows that even with AO supported observations, individual radial velocities in crowded fields are likely to be biased. The ARGUS observations do not show this kind of bias but were found to have larger uncertainties than previously obtained. We find a bias towards higher velocity dispersions in the ARGUS pointing when fixing the extreme velocities of the three brightest stars but find those variations are within the determined uncertainties. We rerun Jeans models and fit the kinematic profile with the new uncertainties. This yields a BH mass of M_BH = (2.8 +- 0.4) x 10^4 M_SUN and M/L ratio M/L = (1.6 +- 0.1) M_SUN/L_SUN, consistent with our previous results.

preprint2015arXiv

Dynamical Constraints on the Origin of Multiple Stellar Populations in Globular Clusters

We have carried out a large grid of N-body simulations in order to investigate if mass-loss as a result of primordial gas expulsion can be responsible for the large fraction of second generation stars in globular clusters (GCs) with multiple stellar populations (MSPs). Our clusters start with two stellar populations in which $10\%$ of all stars are second generation stars. We simulate clusters with different initial masses, different ratios of the half-mass radius of first to second generation stars, different primordial gas fractions and Galactic tidal fields with varying strength. We then let our clusters undergo primordial gas-loss and obtain their final properties such as mass, half-mass radius and the fraction of second generation stars. Using our N-body grid we then perform a Monte Carlo analysis to constrain the initial masses, radii and required gas expulsion time-scales of GCs with MSPs. Our results can explain the present-day properties of GCs only if (1) a substantial amount of gas was present in the clusters after the formation of second generation stars and (2) gas expulsion time-scales were extremely short ($\lesssim 10^5$ yr). Such short gas expulsion time-scales are in agreement with recent predictions that dark remnants have ejected the primordial gas from globular clusters, and pose a potential problem for the AGB scenario. In addition, our results predict a strong anti-correlation between the number ratio of second-generation stars in GCs and the present-day mass of GCs. So far, the observational data show only a significantly weaker anti-correlation, if any at all.

preprint2015arXiv

Dynamical Evolution of Outer-Halo Globular Clusters

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.

preprint2015arXiv

Kinematics of dwarf galaxies in gas-rich groups, and the survival and detectability of tidal dwarf galaxies

We present DEIMOS multi-object spectroscopy (MOS) of 22 star-forming dwarf galaxies located in four gas-rich groups, including six newly-discovered dwarfs. Two of the galaxies are strong tidal dwarf galaxy (TDG) candidates based on our luminosity-metallicity relation definition. We model the rotation curves of these galaxies. Our sample shows low mass-to-light ratios (M/L=0.73$\pm0.39M_\odot/L_\odot$) as expected for young, star-forming dwarfs. One of the galaxies in our sample has an apparently strongly-falling rotation curve, reaching zero rotational velocity outside the turnover radius of $r_{turn}=1.2r_e$. This may be 1) a polar ring galaxy, with a tilted bar within a face-on disk; 2) a kinematic warp. These scenarios are indistinguishable with our current data due to limitations of slit alignment inherent to MOS-mode observations. We consider whether TDGs can be detected based on their tidal radius, beyond which tidal stripping removes kinematic tracers such as H$α$ emission. When the tidal radius is less than about twice the turnover radius, the expected falling rotation curve cannot be reliably measured. This is problematic for as much as half of our sample, and indeed more generally, galaxies in groups like these. Further to this, the H$α$ light that remains must be sufficiently bright to be detected; this is only the case for three (14%) galaxies in our sample. We conclude that the falling rotation curves expected of tidal dwarf galaxies are intrinsically difficult to detect.

preprint2015arXiv

Possible smoking-gun evidence for initial mass segregation in re-virialized post-gas expulsion globular clusters

We perform a series of direct $N$-body calculations to investigate the effect of residual gas expulsion from the gas-embedded progenitors of present-day globular clusters (GCs) on the stellar mass function (MF). Our models start either tidally filling or underfilling, and either with or without primordial mass segregation. We cover 100 Myr of the evolution of modeled clusters and show that the expulsion of residual gas from initially mass-segregated clusters leads to a significantly shallower slope of the stellar MF in the low- ($m\leq 0.50 M_\odot$) and intermediate-mass ($\simeq 0.50-0.85 M_\odot$) regime. Therefore, the imprint of residual gas expulsion and primordial mass segregation might be visible in the present-day MF. We find that the strength of the external tidal field, as an essential parameter, influences the degree of flattening, such that a primordially mass-segregated tidally-filling cluster with $r_h/r_t\geq 0.1$ shows a strongly depleted MF in the intermediate stellar mass range. Therefore, the shape of the present-day stellar MF in this mass range probes the birth place of clusters in the Galactic environment. We furthermore find that this flattening agrees with the observed correlation between the concentration of a cluster and its MF slope, as found by de Marchi et al.. We show that if the expansion through the residual gas expulsion in primordial mass segregated clusters is the reason for this correlation then GCs most probably formed in strongly fluctuating local tidal fields in the early proto-Milky Way potential, supporting the recent conclusion by Marks \& Kroupa.

preprint2015arXiv

Radial distributions of sub-populations in the globular cluster M15: a more centrally concentrated primordial population

We examine the radial distributions of stellar populations in the globular cluster (GC) M15, using HST/WFC3 photometry of red giants in the nitrogen-sensitive F343N-F555W color. Surprisingly, we find that giants with "primordial" composition (i.e., N abundances similar to those in field stars) are the most centrally concentrated within the WFC3 field. We then combine our WFC3 data with SDSS u, g photometry and find that the trend reverses for radii >1' (3 pc) where the ratio of primordial to N-enhanced giants increases outwards, as already found by Lardo et al. The ratio of primordial to enriched stars thus has a U-shaped dependency on radius with a minimum near the half-light radius. N-body simulations show that mass segregation might produce a trend resembling the observed one, but only if the N-enhanced giants are ~0.25 Mo less massive than the primordial giants, which requires extreme He enhancement (Y~0.40). However, such a large difference in Y is incompatible with the negligible optical color differences between primordial and enriched giants which suggest Delta Y < 0.03 and thus a difference in turn-off mass of Delta M < 0.04 Mo between the different populations. The radial trends in M15 are thus unlikely to be of dynamical origin and presumably reflect initial conditions, a result that challenges all current GC formation scenarios. We note that population gradients in the central regions of GCs remain poorly investigated and may show a more diverse behavior than hitherto thought.

preprint2015arXiv

Unveiling Gargantua: A new search strategy for the most massive central cluster black holes

We aim to unveil the most massive central cluster black holes in the universe. We present a new search strategy which is based on a black hole mass gain sensitive 'calorimeter' and which links the innermost stellar density profile of a galaxy to the adiabatic growth of its central SMBH. In a first step we convert observationally inferred feedback powers into SMBH growth rates by using reasonable energy conversion efficiency parameters, $ε$. In the main part of this paper we use these black hole growth rates, sorted in logarithmically increasing steps encompassing our whole parameter space, to conduct $N$-Body computations of brightest cluster galaxies with the newly developed MUESLI software. For the initial setup of galaxies we use core-Sersic models in order to account for SMBH scouring. We find that adiabatically driven core re-growth is significant at the highest accretion rates. As a result, the most massive black holes should be located in BCGs with less pronounced cores when compared to the predictions of empirical scaling relations which are usually calibrated in less extreme environments. For efficiency parameters $ε<0.1$, BCGs in the most massive, relaxed and X-ray luminous galaxy clusters might even develop steeply rising density cusps. Finally, we discuss several promising candidates for follow up investigations, among them the nuclear black hole in the Phoenix cluster. Based on our results, it might have a mass of the order of $10^{11} M_\odot$.

preprint2014arXiv

A prescription and fast code for the long-term evolution of star clusters - III. Unequal masses and stellar evolution

We present a new version of the fast star cluster evolution code Evolve Me A Cluster of StarS (EMACSS). While previous versions of EMACSS reproduced clusters of single-mass stars, this version models clusters with an evolving stellar content. Stellar evolution dominates early evolution, and leads to: (1) reduction of the mean mass of stars due to the mass loss of high-mass stars; (2) expansion of the half-mass radius; (3) for (nearly) Roche Volume filling clusters, the induced escape of stars. Once sufficient relaxation has occurred (~ 10 relaxation times-scales), clusters reach a second, 'balanced' state whereby the core releases energy as required by the cluster as a whole. In this state: (1) stars escape due to tidal effects faster than before balanced evolution; (2) the half-mass radius expands or contracts depending on the Roche volume filling factor; and (3) the mean mass of stars increases due to the preferential ejection of low-mass stars. We compare the EMACSS results of several cluster properties against N-body simulations of clusters spanning a range of initial number of stars, mass, half-mass radius, and tidal environments, and show that our prescription accurately predicts cluster evolution for this database. Finally, we consider applications for EMACSS, such as studies of galactic globular cluster populations in cosmological simulations.

preprint2014arXiv

A Supermassive Black Hole in an Ultracompact Dwarf Galaxy

Ultracompact dwarf galaxies (UCDs) are among the densest stellar systems in the universe. These systems have masses up to 200 million solar masses, but half light radii of just 3-50 parsecs. Dynamical mass estimates show that many UCDs are more massive than expected from their luminosity. It remains unclear whether these high dynamical mass estimates are due to the presence of supermassive black holes or result from a non-standard stellar initial mass function that causes the average stellar mass to be higher than expected. Here we present the detection of a supermassive black hole in a massive UCD. Adaptive optics kinematic data of M60-UCD1 show a central velocity dispersion peak above 100 km/s and modest rotation. Dynamical modeling of these data reveals the presence of a supermassive black hole with mass of 21 million solar masses. This is 15% of the object's total mass. The high black hole mass and mass fraction suggest that M60-UCD1 is the stripped nucleus of a galaxy. Our analysis also shows that M60-UCD1's stellar mass is consistent with its luminosity, implying many other UCDs may also host supermassive black holes. This suggests a substantial population of previously unnoticed supermassive black holes.

preprint2014arXiv

Contribution of stripped nuclear clusters to globular cluster and ultra-compact dwarf galaxy populations

We use the Millennium II cosmological simulation combined with the semi-analytic galaxy formation model of Guo et al. (2011) to predict the contribution of galactic nuclei formed by the tidal stripping of nucleated dwarf galaxies to globular cluster (GC) and ultra-compact dwarf galaxy (UCD) populations of galaxies. We follow the merger trees of galaxies in clusters back in time and determine the absolute number and stellar masses of disrupted galaxies. We assume that at all times nuclei have a distribution in nucleus-to-galaxy mass and nucleation fraction of galaxies similar to that observed in the present day universe. Our results show stripped nuclei follow a mass function $N(M) \sim M^{-1.5}$ in the mass range $10^6 < M/M_\odot < 10^8$, significantly flatter than found for globular clusters. The contribution of stripped nuclei will therefore be most important among high-mass GCs and UCDs. For the Milky Way we predict between 1 and 3 star clusters more massive than $10^5 M_\odot$ come from tidally disrupted dwarf galaxies, with the most massive cluster formed having a typical mass of a few times $10^6 M_\odot$, like omega Centauri. For a galaxy cluster with a mass $7 \times 10^{13} M_\odot$, similar to Fornax, we predict $\sim$19 UCDs more massive than $2\times10^6 M_\odot$ and $\sim$9 UCDs more massive than $10^7 M_\odot$ within a projected distance of 300 kpc come from tidally stripped dwarf galaxies. The observed number of UCDs are $\sim$200 and 23, respectively. We conclude that most UCDs in galaxy clusters are probably simply the high mass end of the GC mass function.

preprint2014arXiv

Direct $N$-body simulations of globular clusters - II. Palomar 4

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.

preprint2014arXiv

Erosion of Globular Cluster Systems: The Influence of Radial Anisotropy, Central Black Holes and Dynamical Friction

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.

preprint2013arXiv

A comparison of the distribution of satellite galaxies around Andromeda and the results of $Λ$CDM simulations

Ibata et al. (2013) recently reported the existence of a vast thin plane of dwarf galaxies (VTPD) orbiting around Andromeda. We investigate whether such a configuration can be reproduced within the standard cosmological framework and search for similar planes of co-rotating satellite galaxies around Andromeda-like host haloes in data from the Millennium II simulation combined with a semi-analytic galaxy formation model. We apply a baryonic mass cut of $2.8\times 10^4 \text{M}_{\text{sun}}$ for the satellite haloes and restrict the data to a PAndAS like field. If we include the so-called orphan galaxies in our analysis, we find that planes with a rms lower than the VTPD are common in Millennium II. This is partially due to the strongly radially concentrated distribution of orphan galaxies. Excluding part of the orphan galaxies brings the radial distributions of Millennium II satellites into better agreement with the satellite distribution of Andromeda while still producing a significant fraction of planes with a lower rms than the VTPD. We also find haloes in Millennium II with an equal or higher number of co-rotating satellites than the VTPD. This demonstrates that the VTPD is not in conflict with the standard cosmological framework, although a definite answer of this question might require higher resolution cosmological simulations that do not have to consider orphan galaxies. Our results finally show that satellite planes in Millennium II are not stable structures, hence the VTPD might only be a statistical fluctuation of an underlying more spherical galaxy distribution.

preprint2013arXiv

A prescription and fast code for the long-term evolution of star clusters - II. Unbalanced and core evolution

We introduce version two of the fast star cluster evolution code Evolve Me A Cluster of StarS (EMACSS). The first version (Alexander & Gieles) assumed that cluster evolution is balanced for the majority of the life-cycle, meaning that the rate of energy generation in the core of the cluster equals the diffusion rate of energy by two-body relaxation, which makes the code suitable for modelling clusters in weak tidal fields. In this new version we extend the model to include an unbalanced phase of evolution to describe the pre-collapse evolution and the accompanying escape rate such that clusters in strong tidal fields can also be modelled. We also add a prescription for the evolution of the core radius and density and a related cluster concentration parameter. The model simultaneously solves a series of first-order ordinary differential equations for the rate of change of the core radius, half-mass radius and the number of member stars N. About two thousand integration steps in time are required to solve for the entire evolution of a star cluster and this number is approximately independent of N. We compare the model to the variation of these parameters following from a series of direct N-body calculations of single-mass clusters and find good agreement in the evolution of all parameters. Relevant time-scales, such as the total lifetimes and core collapse times, are reproduced with an accuracy of about 10% for clusters with various initial half-mass radii (relative to their Jacobi radii) and a range of different initial N up to N = 65536. We intend to extend this framework to include more realistic initial conditions, such as a stellar mass spectrum and mass loss from stars. The EMACSS code can be used in star cluster population studies and in models that consider the co-evolution of (globular) star clusters and large scale structures.

preprint2013arXiv

Intermediate-mass black holes in Globular Clusters

For a sample of nine Galactic globular clusters we measured the inner kinematic profiles with integral-field spectroscopy that we combined with existing outer kinematic measurements and HST luminosity profiles. With this information we are able to detect the crucial rise in the velocity-dispersion profile which indicates the presence of a central black hole. In addition, N-body simulations compared to our data will give us a deeper insight in the properties of clusters with black holes and stronger selection criteria for further studies. For the first time, we obtain a homogeneous sample of globular cluster integral- field spectroscopy which allows a direct comparison between clusters with and without an intermediate-mass black hole.

preprint2013arXiv

N-body simulations of globular clusters in tidal fields: Effects of intermediate-mass black holes

Intermediate-mass black holes (IMBHs) may provide the missing link to understanding the growth of supermassive black holes in the early Universe. Some formation scenarios predict that IMBHs could have formed by runaway collisions in globular clusters (GCs). Understanding the influence of an IMBH in the center of a GC on its environment might provide indirect detection methods. Our goal is to test the effects of different initial compositions of GCs on their evolution in a tidal field. We pin down the crucial observables that indicate the presence of an IMBH at the center of the cluster. In addition to central IMBHs, we also consider the effects of different stellar-mass black hole retention and stellar binary fractions. We performed a set of 22 N-body simulations and varied particle numbers, IMBH masses, stellar-mass black-hole retention fractions, and stellar binary fractions. These models are all run in an external tidal field to study the effect of black holes on the cluster mass loss, mass function, and life times. Finally, we compared our results with observational data. We found that a central massive black hole increases the escape rate of high-mass stars from a star cluster, implying that the relative depletion of the mass function at the low-mass end proceeds less rapidly. Furthermore, we found a similar behavior for a cluster hosting a high number of stellar-mass black holes instead of one massive central IMBH. The presence of an IMBH also weakly affects the fraction of the cluster mass that is constituted by stellar remnants, as does the presence of primordial binaries. We compared our simulations with observational data from the literature and found good agreement between our models and observed mass functions and structural parameters of GCs. We exploited this agreement to identify GCs that could potentially host IMBHs.

preprint2013arXiv

On central black holes in ultra-compact dwarf galaxies

CONTEXT: The dynamical mass-to-light (M/L) ratios of massive ultra-compact dwarf galaxies (UCDs) are about 50% higher than predicted by stellar population models. AIMS: Here we investigate the possibility that these elevated M/L ratios are caused by a central black hole (BH), heating up the internal motion of stars. We focus on a sample of ~50 extragalactic UCDs for which velocity dispersions and structural parameters have been measured. METHODS: Using up-to-date distance moduli and a consistent treatment of aperture and seeing effects, we calculate the ratio Psi=(M/L)_{dyn}/(M/L)_{pop} between the dynamical and the stellar population M/L of UCDs. For all UCDs with Psi>1 we estimate the mass of a hypothetical central BH needed to reproduce the observed integrated velocity dispersion. RESULTS: Massive UCDs (M>10^7 M_*) have an average Psi = 1.7 +-0.2, implying notable amounts of dark mass in them. We find that, on average, central BH masses of 10-15% of the UCD mass can explain these elevated dynamical M/L ratios. The implied BH masses in UCDs range from several 10^5 M_* to several 10^7 M_*. In the M_BH-Luminosity plane, UCDs are offset by about two orders of magnitude in luminosity from the relation derived for galaxies. Our findings can be interpreted such that massive UCDs originate from progenitor galaxies with masses around 10^9 M_*, and that those progenitors have SMBH occupation fractions of 60-100%. The suggested UCD progenitor masses agree with predictions from the tidal stripping scenario. Lower-mass UCDs (M<10^7 M_*) exhibit a bimodal distribution in Psi, suggestive of a coexistence of massive globular clusters and tidally stripped galaxies in this mass regime. CONCLUSIONS: Central BHs as relict tracers of tidally stripped progenitor galaxies are a plausible explanation for the elevated dynamical M/L ratios of UCDs.

preprint2013arXiv

The evolution of the global stellar mass function of star clusters: an analytic description

The evolution of the global stellar mass function (MF) of star clusters is studied based on a large set of N-body simulations of clusters with a range of initial masses, initial concentrations, in circular or elliptical orbits in different tidal environments. Models with and without initial mass segregation are included. The depletion of low mass stars in initially Roche-volume (tidal) filling clusters starts typically on a time scale of the order of the core collapse time. In clusters that are initially underfilling their Roche-volume it takes longer because the clusters have to expand to their tidal radii before dynamical mass loss becomes important. We introduce the concept of the differential mass function (DMF), which describes the changes with respect to the initial mass function (IMF). We show that the evolution of the DMF can be described by a set of very simple analytic expressions that are valid for a wide range of initial cluster parameters and for different IMFs. The agreement between this description and the models is very good, except for initially Roche-volume underfilling clusters that are severely mass segregated.

preprint2013arXiv

The M_BH - sigma relation for intermediate-mass black holes in globular clusters

For galaxies hosting supermassive black holes (SMBHs), it has been observed that the mass of the central black hole (M_BH) tightly correlates with the effective or central velocity dispersion (sigma) of the host galaxy. The origin of this M_BH - sigma scaling relation is assumed to lie in the merging history of the galaxies but many open questions about its origin and the behavior in different mass ranges still need to be addressed. The goal of this work is to study the black-hole scaling relations for low black-hole masses, where the regime of intermediate-mass black holes (IMBHs) in globular clusters (GCs) is entered. We collect all existing reports of dynamical black-hole measurements in globular clusters, providing black-hole masses or upper limits for 14 candidates. We plot the black-hole masses versus different cluster parameters including total mass, velocity dispersion, concentration and half-mass radius. We search for trends and test the correlations in order to quantify their significance using a set of different statistical approaches. For correlations showing a large significance we perform a linear fit, accounting for uncertainties and upper limits. We find a clear correlation between the mass of the IMBH and the velocity dispersion of the globular cluster. As expected, the total mass of the globular cluster then also correlates with the mass of the IMBH. While the slope of the M_BH - sigma correlation differs strongly from the one observed for SMBHs, the other scaling relations M_BH - M_TOT, and M_BH - L are similar to the correlations in galaxies. Significant correlations of black-hole mass with other cluster properties were not found in the present sample.

preprint2013arXiv

The stellar mass function, binary content and radial structure of the open cluster Praesepe derived from PPMXL and SDSS data

We have determined possible cluster members of the nearby open cluster Praesepe (M44) based on J and K photometry and proper motions from the PPMXL catalogue and z photometry from the Sloan Digital Sky Survey (SDSS). In total we identified 893 possible cluster members down to a magnitude of J = 15.5 mag, corresponding to a mass of about 0.15 Msun for an assumed cluster distance modulus of (m-M)_0 = 6.30 mag (d ~ 182 pc), within a radius of 3.5° around the cluster centre. We derive a new cluster centre for Praesepe (α_centre = 8h 39m 37s, δ_centre = 19° 35' 02"). We also derive a total cluster mass of about 630 Msun and a 2D half-number and half-mass radius of 4.25 pc and 3.90 pc respectively. The global mass function (MF) of the cluster members shows evidence for a turnover around m = 0.65 Msun. While more massive stars can be fit by a power-law ξ(m) ~ m^(-α) with slope α = 2.88 +/- 0.22, stars less massive than m = 0.65 Msun are best fitted with α = 0.85 +/- 0.10. In agreement with its large dynamical age, we find that Praesepe is strongly mass segregated and that the mass function slope for high mass stars steepens from a value of α = 2.32 +/- 0.24 inside the half-mass radius to α = 4.90 +/- 0.51 outside the half-mass radius. We finally identify a significant population of binaries and triples in the colour-magnitude diagram of Praesepe. Assuming non-random pairing of the binary components, a binary fraction of about 35% for primaries in the mass range 0.6 < m/Msun < 2.20 is required to explain the observed number of binaries in the colour-magnitude diagram (CMD).

preprint2013arXiv

Ultra compact dwarf galaxy formation by tidal stripping of nucleated dwarf galaxies

Ultra Compact Dwarf Galaxies (UCDs) and dwarf galaxy nuclei have many common properties, such as internal velocity dispersions and colour-magnitude trends, suggesting tidally stripped dwarf galaxies as a possible UCD origin. However, UCDs typically have sizes more than twice as large as nuclei at the same luminosity. We use a GPU-enabled version of the particle-mesh code \textsc{superbox} to study the possibility of turning nucleated dwarf galaxies into UCDs by tidally stripping them in a Virgo-like galaxy cluster. We find that motion in spherical potentials, where close passages happen many times, leads to the formation of compact ($r_h \lesssim 20$ pc) star clusters/UCDs. In contrast, orbital motion where close passages happen only once or twice leads to the formation of extended objects which are large enough to account for the full range of observed UCD sizes. For such motion, we find that dwarf galaxies need close pericentre passages with distances less than 10 kpc to undergo strong enough stripping so that UCD formation is possible. As tidal stripping produces objects with similar properties to UCDs, and our estimates suggest dwarf galaxies have been destroyed in sufficient numbers to explain the observed number of UCDs in M87, we consider tidal stripping to be a likely origin of UCDs. However, comparison with cosmological simulations is needed to determine if the number and spatial distribution of UCDs formed by tidal stripping matches the observations of UCDs in galaxy clusters.

preprint2012arXiv

Catch me if you can: is there a runaway-mass black hole in the Orion Nebula Cluster?

We investigate the dynamical evolution of the Orion Nebula Cluster (ONC) by means of direct N-body integrations. A large fraction of residual gas was probably expelled when the ONC formed, so we assume that the ONC was much more compact when it formed compared to its current size, in agreement with the embedded cluster radius-mass relation from Marks & Kroupa (2012). Hence, we assume that few-body relaxation played an important role during the initial phase of evolution of the ONC. In particular, three body interactions among OB stars likely led to their ejection from the cluster and, at the same time, to the formation of a massive object via runaway physical stellar collisions. The resulting depletion of the high mass end of the stellar mass function in the cluster is one of the important points where our models fit the observational data. We speculate that the runaway-mass star may have collapsed directly into a massive black hole (Mbh > 100Msun). Such a dark object could explain the large velocity dispersion of the four Trapezium stars observed in the ONC core. We further show that the putative massive black hole is likely to be a member of a binary system with appr. 70 per cent probability. In such a case, it could be detected either due to short periods of enhanced accretion of stellar winds from the secondary star during pericentre passages, or through a measurement of the motion of the secondary whose velocity would exceed 10 km/s along the whole orbit.

preprint2012arXiv

Central kinematics of the globular cluster NGC 2808: Upper limit on the mass of an intermediate-mass black hole

Globular clusters are an excellent laboratory for stellar population and dynamical research. Recent studies have shown that these stellar systems are not as simple as previously assumed. With multiple stellar populations as well as outer rotation and mass segregation they turn out to exhibit high complexity. This includes intermediate-mass black holes which are proposed to sit at the centers of some massive globular clusters. Today's high angular resolution ground based spectrographs allow velocity-dispersion measurements at a spatial resolution comparable to the radius of influence for plausible IMBH masses, and to detect changes in the inner velocity-dispersion profile. Together with high quality photometric data from HST, it is possible to constrain black-hole masses by their kinematic signatures. We determine the central velocity-dispersion profile of the globular cluster NGC 2808 using VLT/FLAMES spectroscopy. In combination with HST/ACS data our goal is to probe whether this massive cluster hosts an intermediate-mass black hole at its center and constrain the cluster mass to light ratio as well as its total mass. We derive a velocity-dispersion profile from integral field spectroscopy in the center and Fabry Perot data for larger radii. High resolution HST data are used to obtain the surface brightness profile. Together, these data sets are compared to dynamical models with varying parameters such as mass to light ratio profiles and black-hole masses. Using analytical Jeans models in combination with variable M/L profiles from N-body simulations we find that the best fit model is a no black hole solution. After applying various Monte Carlo simulations to estimate the uncertainties, we derive an upper limit of the back hole mass of M_BH < 1 x 10^4 M_SUN (with 95 % confidence limits) and a global mass-to-light ratio of M/L_V = (2.1 +- 0.2) M_SUN/L_SUN.

preprint2012arXiv

High-velocity stars in the cores of globular clusters: The illustrative case of NGC 2808

We report the detection of five high-velocity stars in the core of the globular cluster NGC 2808. The stars lie on the the red giant branch and show total velocities between 40 and 45 km/s. For a core velocity dispersion sigma_c = 13.4 km/s, this corresponds to up to 3.4 sigma_c. These velocities are close to the estimated escape velocity (~ 50 km/s) and suggest an ejection from the core. Two of these stars have been confirmed in our recent integral field spectroscopy data and we will discuss them in more detail here. These two red giants are located at a projected distance of ~ 0.3 pc from the center. According to their positions on the color magnitude diagram, both stars are cluster members. We investigate several possible origins for the high velocities of the stars and conceivable ejection mechanisms. Since the velocities are close to the escape velocity, it is not obvious whether the stars are bound or unbound to the cluster. We therefore consider both cases in our analysis. We perform numerical simulations of three-body dynamical encounters between binaries and single stars and compare the resulting velocity distributions of escapers with the velocities of our stars. We compare the predictions for a single dynamical encounter with a compact object with those of a sequence of two-body encounters due to relaxation. If the stars are unbound, the encounter must have taken place recently, when the stars were already in the giant phase. After including binary fractions and black-hole retention fractions, projection effects, and detection probabilities from Monte-Carlo simulations, we estimate the expected numbers of detections for all the different scenarios. Based on these numbers, we conclude that the most likely scenario is that the stars are bound and were accelerated by a single encounter between a binary of main-sequence stars and a ~ 10 M_sun black hole.

preprint2012arXiv

Limits on intermediate-mass black holes in six Galactic globular clusters with integral-field spectroscopy

The formation of supermassive black holes at high redshift still remains a puzzle to astronomers. Their growth becomes reasonable only when starting from a massive seed black hole with mass of the order of 10^2 - 10^5 M_SUN. Intermediate-mass black holes (IMBHs) are therefore an important field of research. Especially the possibility of finding them in the centers of globular clusters has recently drawn attention. The search for IMBHs in the centers of globular clusters could therefore shed light on the process of black-hole formation and cluster evolution. We are investigating six galactic globular clusters for the presence of an IMBH at their centers. Based on their kinematic and photometric properties, we selected the globular clusters NGC 1851, NGC 1904 (M79), NGC 5694, NGC 5824, NGC 6093 (M80) and NGC 6266 (M62). We use integral field spectroscopy in order to obtain the central velocity-dispersion profile of each cluster. We compute the cluster photometric center and the surface brightness profile using HST data. After combining these datasets we compare them to analytic Jeans models. We use varying M/L_V profiles for clusters with enough data points in order to reproduce their kinematic profiles in an optimal way. Finally, we vary the mass of the central black hole and test whether the cluster is better fitted with or without an IMBH. We present the statistical significance, including upper limits, of the black-hole mass for each cluster. NGC 1904 and NGC 6266 provide the highest significance for a black hole. Jeans models in combination with a M/L_V profile obtained from N-body simulations (in the case of NGC 6266) predict a central black hole of M_BH = (3 +- 1) x 10^3 M_SUN for NGC 1904 and M_BH = (2 +- 1) x 10^3 M_SUN for NGC 6266. Furthermore, we discuss the possible influence of dark remnants and mass segregation at the center of the cluster on the detection of an IMBH.

preprint2012arXiv

The mass function and dynamical mass of young star clusters: Why their initial crossing-time matters crucially

We highlight the impact of cluster-mass-dependent evolutionary rates upon the evolution of the cluster mass function during violent relaxation, that is, while clusters dynamically respond to the expulsion of their residual star-forming gas. Mass-dependent evolutionary rates arise when the mean volume density of cluster-forming regions is mass-dependent. In that case, even if the initial conditions are such that the cluster mass function at the end of violent relaxation has the same shape as the embedded-cluster mass function (i.e. infant weight-loss is mass-independent), the shape of the cluster mass function does change transiently {\it during} violent relaxation. In contrast, for cluster-forming regions of constant mean volume density, the cluster mass function shape is preserved all through violent relaxation since all clusters then evolve at the same mass-independent rate. On the scale of individual clusters, we model the evolution of the ratio between the dynamical mass and luminous mass of a cluster after gas expulsion. Specifically, we map the radial dependence of the time-scale for a star cluster to return to equilibrium. We stress that fields-of-view a few pc in size only, typical of compact clusters with rapid evolutionary rates, are likely to reveal cluster regions which have returned to equilibrium even if the cluster experienced a major gas expulsion episode a few Myr earlier. We provide models with the aperture and time expressed in units of the initial half-mass radius and initial crossing-time, respectively, so that our results can be applied to clusters with initial densities, sizes, and apertures different from ours.

preprint2012arXiv

The star cluster formation history of the LMC

The Large Magellanic Cloud is one of the nearest galaxies to us and is one of only few galaxies where the star formation history can be determined from studying resolved stellar populations. We have compiled a new catalogue of ages, luminosities and masses of LMC star clusters and used it to determine the age distribution and dissolution rate of LMC star clusters. We find that the frequency of massive clusters with masses M>5000 Msun is almost constant between 10 and 200 Myr, showing that the influence of residual gas expulsion is limited to the first 10 Myr of cluster evolution or clusters less massive than 5000 Msun. Comparing the cluster frequency in that interval with the absolute star formation rate, we find that about 15% of all stars in the LMC were formed in long-lived star clusters that survive for more than 10 Myr. We also find that the mass function of LMC clusters younger than 1 Gyr can be fitted by a power-law mass function with slope α=-2.3, while older clusters follow a significantly shallower slope and interpret this is a sign of the ongoing dissolution of low-mass clusters. Our data shows that for ages older than 200 Myr, about 90% of all clusters are lost per dex of lifetime. The implied cluster dissolution rate is significantly faster than that based on analytic estimates and N-body simulations. Our cluster age data finally shows evidence for a burst in cluster formation about 1 Gyr ago, but little evidence for bursts at other ages.

preprint2011arXiv

Distant star clusters of the Milky Way in MOND

We determine the mean velocity dispersion of six Galactic outer halo globular clusters, AM 1, Eridanus, Pal 3, Pal 4, Pal 15, and Arp 2 in the weak acceleration regime to test classical vs. modified Newtonian dynamics (MOND). Owing to the non-linearity of MOND's Poisson equation, beyond tidal effects, the internal dynamics of clusters is affected by the external field in which they are immersed. For the studied clusters, particle accelerations are much lower than the critical acceleration a_0 of MOND, but the motion of stars is neither dominated by internal accelerations (a_i >> a_e) nor external accelerations (a_e >> a_i). We use the N-body code N-MODY in our analysis, which is a particle-mesh-based code with a numerical MOND potential solver developed by Ciotti, Londrillo, and Nipoti (2006) to derive the line-of-sight velocity dispersion by adding the external field effect. We show that Newtonian dynamics predicts a low-velocity dispersion for each cluster, while in modified Newtonian dynamics the velocity dispersion is much higher. We calculate the minimum number of measured stars necessary to distinguish between Newtonian gravity and MOND with the Kolmogorov-Smirnov test. We also show that for most clusters it is necessary to measure the velocities of between 30 to 80 stars to distinguish between both cases. Therefore the observational measurement of the line-of-sight velocity dispersion of these clusters will provide a test for MOND.

preprint2011arXiv

Testing Photometric Diagnostics for the Dynamical State and Possible IMBH presence in Globular Clusters

Surface photometry is a necessary tool to establish the dynamical state of stars clusters. We produce realistic HST-like images from N-body models of star clusters with and without central intermediate-mass black holes (IMBHs) in order to measure their surface brightness profiles. The models contain ~600,000 individual stars, black holes of various masses between 0% to 2% of the total mass, and are evolved for a Hubble time. We measure surface brightness and star count profiles for every constructed image in order to test the effect of intermediate mass black holes on the central logarithmic slope, the core radius, and the half-light radius. We use these quantities to test diagnostic tools for the presence of central black holes using photometry. We find that the the only models that show central shallow cusps with logarithmic slopes between -0.1 and -0.4 are those containing central black holes. Thus, the central logarithmic slope seems to be a good way to choose clusters suspect of containing intermediate-mass black holes. Clusters with steep central cusps can definitely be ruled out to host an IMBH. The measured r_c/r_h ratio has similar values for clusters that have not undergone core-collapse, and those containing a central black hole. We notice that observed Galactic globular clusters have a larger span of values for central slope and r_c/r_h than our modeled clusters, and suggest possible reasons that could account for this and contribute to improve future models.

preprint2011arXiv

Tidal mass loss in star clusters and treatment of escapers in Fokker-Planck models

This paper presents a new scheme to treat escaping stars in the orbit-averaged Fokker-Planck models of globular star clusters in a galactic tidal field. The existence of a large number of potential escapers, which have energies above the escape energy but are still within the tidal radius, is taken into account in the models. The models allow potential escapers to experience gravitational scatterings before they leave clusters and thus some of them may lose enough energy to be bound again. It is shown that the mass evolution of the Fokker-Planck models are in good agreement with that of N-body models including the full tidal-force field. The mass-loss time does not simply scale with the relaxation time due to the existence of potential escapers; it increases with the number of stars more slowly than the relaxation time, though it tends to be proportional to the relaxation time in the limit of a weak tidal field. The Fokker-Planck models include two parameters, the coefficient gamma in the Coulomb logarithm ln(gamma N) and the coefficient nu_e controlling the efficiency of the mass loss. The values of these parameters are determined by comparing the Fokker-Planck models with the N-body models. It is found that the parameter set (gamma, nu_e)=(0.11, 7) works well for both single-mass and multi-mass clusters, but that the parameter set (gamma, nu_e)=(0.02, 40) is another possible choice for multi-mass clusters.

preprint2010arXiv

Direct N-body simulations of globular clusters: (I) Palomar 14

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.

preprint2010arXiv

On the mass-radius relation of hot stellar systems

Most globular clusters have half-mass radii of a few pc with no apparent correlation with their masses. This is different from elliptical galaxies, for which the Faber-Jackson relation suggests a strong positive correlation between mass and radius. Objects that are somewhat in between globular clusters and low-mass galaxies, such as ultra-compact dwarf galaxies, have a mass-radius relation consistent with the extension of the relation for bright ellipticals. Here we show that at an age of 10 Gyr a break in the mass-radius relation at ~10^6 Msun is established because objects below this mass, i.e. globular clusters, have undergone expansion driven by stellar evolution and hard binaries. From numerical simulations we find that the combined energy production of these two effects in the core comes into balance with the flux of energy that is conducted across the half-mass radius by relaxation. An important property of this `balanced' evolution is that the cluster half-mass radius is independent of its initial value and is a function of the number of bound stars and the age only. It is therefore not possible to infer the initial mass-radius relation of globular clusters and we can only conclude that the present day properties are consistent with the hypothesis that all hot stellar systems formed with the same mass-radius relation and that globular clusters have moved away from this relation because of a Hubble time of stellar and dynamical evolution.

preprint2010arXiv

VLT Kinematics for omega Centauri: Further Support for a Central Black Hole

The Galactic globular cluster omega Centauri is a prime candidate for hosting an intermediate mass black hole. Recent measurements lead to contradictory conclusions on this issue. We use VLT-FLAMES to obtain new integrated spectra for the central region of omega Centauri. We combine these data with existing measurements of the radial velocity dispersion profile taking into account a new derived center from kinematics and two different centers from the literature. The data support previous measurements performed for a smaller field of view and show a discrepancy with the results from a large proper motion data set. We see a rise in the radial velocity dispersion in the central region to 22.8+-1.2 km/s, which provides a strong sign for a central black hole. Isotropic dynamical models for omega Centauri imply black hole masses ranging from 3.0 to 5.2x10^4 solar masses depending on the center. The best-fitted mass is 4.7+-1.0x10^4 solar masses.

preprint2009arXiv

Constraining the initial mass function of stars in the Galactic Centre

(abridged) Here we discuss the question whether the extreme circumstances in the centre of the Milky Way may be the reason for a significant variation of the IMF. By means of stellar evolution models using different codes we show that the observed luminosity in the central parsec is too high to be explained by a long-standing top-heavy IMF, considering the limited amount of mass inferred from stellar kinematics in this region. In contrast, continuous star formation over the Galaxy's lifetime following a canonical IMF results in a mass-to-light ratio and a total mass of stellar black holes (SBHs) consistent with the observations. Furthermore, these SBHs migrate towards the centre due to dynamical friction, turning the cusp of visible stars into a core as observed in the Galactic Centre. For the first time here we explain the luminosity and dynamical mass of the central cluster and both the presence and extent of the observed core, since the number of SBHs expected from a canonical IMF is just enough to make up for the missing luminous mass. We conclude that the Galactic Centre is consistent with the canonical IMF and do not suggest a systematic variation as a result of the region's properties such as high density, metallicity, strong tidal field etc.

preprint2009arXiv

Evidence for two populations of Galactic globular clusters from the ratio of their half-mass to Jacobi radii

We investigate the ratio between the half-mass radii r_h of Galactic globular clusters and their Jacobi radii r_J given by the potential of the Milky Way and show that clusters with galactocentric distances R_{GC}>8 kpc fall into two distinct groups: one group of compact, tidally-underfilling clusters with r_h/r_J<0.05 and another group of tidally filling clusters which have 0.1 < r_h/r_J<0.3. We find no correlation between the membership of a particular cluster to one of these groups and its membership in the old or younger halo population. Based on the relaxation times and orbits of the clusters, we argue that compact clusters and most clusters in the inner Milky Way were born compact with half-mass radii r_h < 1 pc. Some of the tidally-filling clusters might have formed compact as well, but the majority likely formed with large half-mass radii. Galactic globular clusters therefore show a similar dichotomy as was recently found for globular clusters in dwarf galaxies and for young star clusters in the Milky Way. It seems likely that some of the tidally-filling clusters are evolving along the main sequence line of clusters recently discovered by Kuepper et al. (2008) and are in the process of dissolution.

preprint2009arXiv

Globular cluster systems in nearby dwarf galaxies - II. Nuclear star clusters and their relation to massive Galactic globular clusters

(Abridged) Using luminosities and structural parameters of globular clusters (GCs) in the nuclear regions (nGCs) of low-mass dwarf galaxies from HST/ACS imaging we derive the present-day escape velocities (v_esc) of stellar ejecta to reach the cluster tidal radius and compare them with those of Galactic GCs with extended (hot) horizontal branches (EHBs-GCs). For EHB-GCs, we find a correlation between the present-day v_esc and their metallicity as well as (V-I)_0 colour. The similar v_esc, (V-I)_0 distribution of nGCs and EHB-GCs implies that nGCs could also have complex stellar populations. The v_esc-[Fe/H] relation could reflect the known relation of increasing stellar wind velocity with metallicity, which in turn could explain why more metal-poor clusters typically show more peculiarities in their stellar population than more metal-rich clusters of the same mass do. Thus the cluster v_esc can be used as parameter to describe the degree of self-enrichment. The nGCs populate the same Mv vs. rh region as EHB-GCs, although they do not reach the sizes of the largest EHB-GCs like wCen and NGC 2419. We argue that during accretion the rh of an nGC could increase due to significant mass loss in the cluster vicinity and the resulting drop in the external potential in the core once the dwarf galaxy dissolves. Our results support the scenario in which Galactic EHB-GCs have originated in the centres of pre-Galactic building blocks or dwarf galaxies that were later accreted by the Milky Way.

preprint2009arXiv

Influence of a stellar cusp on the dynamics of young stellar discs and the origin of the S-stars in the Galactic Centre

Observations of the Galactic Centre show evidence of one or two disc-like structures of very young stars orbiting the central super-massive black hole within a distance of a few 0.1 pc. A number of analyses have been carried out to investigate the dynamical behaviour and consequences of these discs, including disc thickness and eccentricity growth as well as mutual interaction and warping. However, most of these studies have neglected the influence of the stellar cusp surrounding the black hole, which is believed to be 1-2 orders of magnitude more massive than the disc(s). By means of N-body integrations using our bhint code, we study the impact of stellar cusps of different compositions. We find that although the presence of a cusp does have an important effect on the evolution of an otherwise isolated flat disc, its influence on the evolution of disc thickness and warping is rather mild in a two-disc configuration. However, we show that the creation of highly eccentric orbits strongly depends on the graininess of the cusp (i.e. the mean and maximum stellar masses): While Chang (2009) recently found that full cycles of Kozai resonance are prevented by the presence of an analytic cusp, we show that relaxation processes play an important role in such highly dense regions and support short-term resonances. We thus find that young disc stars on initially circular orbits can achieve high eccentricities by resonant effects also in the presence of a cusp of stellar remnants, yielding a mechanism to create S-stars and hyper-velocity stars. Furthermore, we discuss the underlying initial mass function (IMF) of the young stellar discs and find no definite evidence for a non-canonical IMF.

preprint2009arXiv

On the number of young stellar discs in the Galactic Centre

Observations of the Galactic Centre show evidence of disc-like structures of very young stars orbiting the central super-massive black hole within a distance of a few 0.1 pc. While it is widely accepted that about half of the stars form a relatively flat disc rotating clockwise on the sky, there is a substantial ongoing debate on whether there is a second, counter-clockwise disc of stars. By means of N-body simulations using our bhint code, we show that two highly inclined stellar discs with the observed properties cannot be recognised as two flat circular discs after 5 Myr of mutual interaction. Instead, our calculations predict a significant warping of the two discs, which we show to be apparent among the structures observed in the Galactic Centre. While the high eccentricities of the observed counter-clockwise orbits suggest an eccentric origin of this system, we show the eccentricity distribution in the inner part of the more massive clockwise disc to be perfectly consistent with an initially circular disc in which stellar eccentricities increase due to both non-resonant and resonant relaxation. We conclude that the relevant question to ask is therefore not whether there are two discs of young stars, but whether there were two such discs to begin with.

preprint2009arXiv

Stellar-mass black holes in star clusters: implications for gravitational wave radiation

We study the dynamics of stellar-mass black holes (BH) in star clusters with particular attention to the formation of BH-BH binaries, which are interesting as sources of gravitational waves (GW). We examine the properties of these BH-BH binaries through direct N-body simulations of star clusters using the GPU-enabled NBODY6 code. We perform simulations of N <= 10^5 Plummer clusters of low-mass stars with an initial population of BHs. Additionally, we do several calculations of star clusters confined within a reflective boundary mimicking only the core of a massive cluster. We find that stellar-mass BHs with masses ~ 10 solar mass segregate rapidly into the cluster core and form a sub-cluster of BHs within typically 0.2 - 0.5 pc radius, which is dense enough to form BH-BH binaries through 3-body encounters. While most BH binaries are ejected from the cluster by recoils received during super-elastic encounters with the single BHs, few of them harden sufficiently so that they can merge via GW emission within the cluster. We find that for clusters with $N \ga 5\times 10^4$, typically 1 - 2 BH-BH mergers occur within them during the first ~ 4 Gyr of evolution. Also for each of these clusters, there are a few escaping BH binaries that can merge within a Hubble time, most of the merger times being within a few Gyr. These results indicate that intermediate-age massive clusters constitute the most important class of candidates for producing dynamical BH-BH mergers. Old globular clusters cannot contribute significantly to the present-day BH-BH merger rate since most of the mergers from them would have occurred earlier. In contrast, young massive clusters are too young to produce significant number of BH-BH mergers. Our results imply significant BH-BH merger detection rates for the proposed "Advanced LIGO" GW detector. (Abridged)

preprint2009arXiv

Testing fundamental physics with distant star clusters: theoretical models for pressure-supported stellar systems

We investigate the mean velocity dispersion and the velocity dispersion profile of stellar systems in MOND, using the N-body code N-MODY, which is a particle-mesh based code with a numerical MOND potential solver developed by Ciotti, Londrillo and Nipoti (2006). We have calculated mean velocity dispersions for stellar systems following Plummer density distributions with masses in the range of $10^4 M_\odot$ to $10^9 M_\odot$ and which are either isolated or immersed in an external field. Our integrations reproduce previous analytic estimates for stellar velocities in systems in the deep MOND regime ($a_i, a_e \ll a_0$), where the motion of stars is either dominated by internal accelerations ($a_i \gg a_e$) or constant external accelerations ($a_e \gg a_i$). In addition, we derive for the first time analytic formulae for the line-of-sight velocity dispersion in the intermediate regime ($a_i \sim a_e \sim a_0$). This allows for a much improved comparison of MOND with observed velocity dispersions of stellar systems. We finally derive the velocity dispersion of the globular cluster Pal 14 as one of the outer Milky Way halo globular clusters that have recently been proposed as a differentiator between Newtonian and MONDian dynamics.

preprint2007arXiv

The core radius of a star cluster containing a massive black hole

(Simplified) We present a theoretical framework which establishes how the core radius of a star cluster varies with the mass of an assumed central black hole. Our result is that the ratio of core to half-mass radius varies as the 3/4 power of the ratio of the black hole to cluster mass, when the system is well relaxed. The theory compares favourably with a number of simulations of this problem, which extend to black hole masses of order 10% of the cluster mass. Though strictly limited as yet to clusters with stars of equal mass, our conclusion strengthens the view that clusters with large core radii are the most promising candidates in which to find a massive black hole.

preprint2004arXiv

The formation of massive black holes through collision runaway in dense young star clusters

A luminous X-ray source is associated with a cluster (MGG-11) of young stars \~200pc from the center of the starburst galaxy M82. The properties of the X-ray source are best explained by a black hole with a mass of at least 350Msun, which is intermediate between stellar-mass and supermassive black holes. A nearby but somewhat more massive star cluster (MGG-9) shows no evidence of such an intermediate mass black hole, raising the issue of just what physical characteristics of the clusters can account for this difference. Here we report numerical simulations of the evolution and the motions of stars within the clusters, where stars are allowed to mergers with each other. We find that for MGG-11 dynamical friction leads to the massive stars sinking rapidly to the center of the cluster to participate in a runaway collision, thereby producing a star of 800-3000Msun, which ultimately collapses to an black hole of intermediate mass. No such runaway occurs in the cluster MGG-9 because the larger cluster radius leads to a mass-segregation timescale a factor of five longer than for MGG-11.