Source author record

Mario Pasquato

Mario Pasquato 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

22works
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

22 published item(s)

preprint2022arXiv

3D Morphology of Open Clusters in the Solar Neighborhood with Gaia EDR3 II: Hierarchical Star Formation Revealed by Spatial and Kinematic Substructures

We identify members of 65 open clusters in the solar neighborhood using the machine-learning algorithm StarGO based on Gaia EDR3 data. After adding members of twenty clusters from previous studies (Pang et al. 2021a,b; Li et al. 2021) we obtain 85 clusters, and study their morphology and kinematics. We classify the substructures outside the tidal radius into four categories: filamentary (f1) and fractal (f2) for clusters $<100$ Myr, and halo (h) and tidal-tail (t) for clusters $>100$ Myr. The kinematical substructures of f1-type clusters are elongated; these resemble the disrupted cluster Group X. Kinematic tails are distinct in t-type clusters, especially Pleiades. We identify 29 hierarchical groups in four young regions (Alessi 20, IC 348, LP 2373, LP 2442); ten among these are new. The hierarchical groups form filament networks. Two regions (Alessi 20, LP 2373) exhibit global "orthogonal" expansion (stellar motion perpendicular to the filament), which might cause complete dispersal. Infalling-like flows (stellar motion along the filament) are found in UBC 31 and related hierarchical groups in the IC 348 region. Stellar groups in the LP 2442 region (LP 2442 gp 1-5) are spatially well-mixed but kinematically coherent. A merging process might be ongoing in the LP 2442 subgroups. For younger systems ($\lesssim30$ Myr), the mean axis ratio, cluster mass and half-mass radius tend to increase with age values. These correlations between structural parameters may imply two dynamical processes occurring in the hierarchical formation scenario in young stellar groups: (1) filament dissolution and (2) sub-group mergers.

preprint2022arXiv

Exploring X-ray variability with unsupervised machine learning I. Self-organizing maps applied to XMM-Newton data

XMM-Newton provides unprecedented insight into the X-ray Universe, recording variability information for hundreds of thousands of sources. Manually searching for interesting patterns in light curves is impractical, requiring an automated data-mining approach for the characterization of sources. Straightforward fitting of temporal models to light curves is not a sure way to identify them, especially with noisy data. We used unsupervised machine learning to distill a large data set of light-curve parameters, revealing its clustering structure in preparation for anomaly detection and subsequent searches for specific source behaviors (e.g., flares, eclipses). Self-organizing maps (SOMs) achieve dimensionality reduction and clustering within a single framework. They are a type of artificial neural network trained to approximate the data with a two-dimensional grid of discrete interconnected units, which can later be visualized on the plane. We trained our SOM on temporal-only parameters computed from more than 100,000 detections from the EXTraS catalog. The resulting map reveals that about 2500 most variable sources are clustered based on temporal characteristics. We find distinctive regions of the SOM map associated with flares, eclipses, dips, linear light curves, and others. Each group contains sources that appear similar by eye. We single out a handful of interesting sources for further study. The condensed view of our dataset provided by SOMs allowed us to identify groups of similar sources, speeding up manual characterization by orders of magnitude. Our method also highlights problems with fitting simple temporal models to light curves and can be used to mitigate them to an extent. This will be crucial for fully exploiting the high data volume expected from upcoming X-ray surveys, and may also help with interpreting supervised classification models.

preprint2022arXiv

NIHAO XXVIII: Collateral effects of AGN on dark matter concentration and stellar kinematics

Although active galactic nuclei (AGN) feedback is required in simulations of galaxies to regulate star formation, further downstream effects on the dark matter distribution of the halo and stellar kinematics of the central galaxy can be expected. We combine simulations of galaxies with and without AGN physics from the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) to investigate the effect of AGN on the dark matter profile and central stellar rotation of the host galaxies. Specifically, we study how the concentration-halo mass ($c-M$) relation and the stellar spin parameter ($λ_R$) are affected by AGN feedback. We find that AGN physics is crucial to reduce the central density of simulated massive ($\gtrsim 10^{12}$ M$_\odot$) galaxies and bring their concentration to agreement with results from the Spitzer Photometry & Accurate Rotation Curves (SPARC) sample. Similarly, AGN feedback has a key role in reproducing the dichotomy between slow and fast rotators as observed by the ATLAS$^{3\text{D}}$ survey. Without star formation suppression due to AGN feedback, the number of fast rotators strongly exceeds the observational constraints. Our study shows that there are several collateral effects that support the importance of AGN feedback in galaxy formation, and these effects can be used to constrain its implementation in numerical simulations.

preprint2021arXiv

Introducing a new multi-particle collision method for the evolution of dense stellar systems II. Core collapse

In a previous paper we introduced a new method for simulating collisional gravitational $N$-body systems with linear time scaling on $N$, based on the Multi-Particle Collision (MPC) approach. This allows us to simulate globular clusters with a realistic number of stellar particles in a matter of hours on a typical workstation. We evolve star clusters containing up to $10^6$ stars to core collapse and beyond. We quantify several aspects of core collapse over multiple realizations and different parameters, while always resolving the cluster core with a realistic number of particles. We run a large set of N-body simulations with our new code. The cluster mass function is a power-law with no stellar evolution, allowing us to clearly measure the effects of the mass spectrum. Leading up to core collapse, we find a power-law relation between the size of the core and the time left to core collapse. Our simulations thus confirm the theoretical self-similar contraction picture but with a dependence on the slope of the mass function. The time of core collapse has a non-monotonic dependence on the slope, which is well fit by a parabola. This holds also for the depth of core collapse and for the dynamical friction timescale of heavy particles. Cluster density profiles at core collapse show a broken power law structure, suggesting that central cusps are a genuine feature of collapsed cores. The core bounces back after collapse, and the inner density slope evolves to an asymptotic value. The presence of an intermediate-mass black hole inhibits core collapse. We confirm and expand on several predictions of star cluster evolution before, during, and after core collapse. Such predictions were based on theoretical calculations or small-size direct $N$-body simulations. Here we put them to the test on MPC simulations with a much larger number of particles, allowing us to resolve the collapsing core.

preprint2020arXiv

An astrophysically motivated ranking criterion for low-latency electromagnetic follow-up of gravitational wave events

We investigate the properties of the host galaxies of compact binary mergers across cosmic time. To this end, we combine population synthesis simulations together with galaxy catalogues from the hydrodynamical cosmological simulation EAGLE to derive the properties of the host galaxies of binary neutron star (BNS), black hole-neutron star (BHNS) and binary black hole (BBH) mergers. Within this framework, we derive the host galaxy probability, i.e., the probability that a galaxy hosts a compact binary coalescence as a function of its stellar mass, star formation rate, $K_s$ magnitude and $B$ magnitude. This quantity is particularly important for low-latency searches of gravitational wave (GW) sources as it provides a way to rank galaxies lying inside the credible region in the sky of a given GW detection, hence reducing the number of viable host candidates. Furthermore, even if no electromagnetic counterpart is detected, the proposed ranking criterion can still be used to classify the galaxies contained in the error box. Our results show that massive galaxies (or equivalently galaxies with a high luminosity in $K_s$ band) have a higher probability of hosting BNS, BHNS, and BBH mergers. We provide the probabilities in a suitable format to be implemented in future low-latency searches.

preprint2020arXiv

Binary black holes in young star clusters: the impact of metallicity

Young star clusters are the most common birth-place of massive stars and are dynamically active environments. Here, we study the formation of black holes (BHs) and binary black holes (BBHs) in young star clusters, by means of 6000 N-body simulations coupled with binary population synthesis. We probe three different stellar metallicities (Z=0.02, 0.002 and 0.0002) and two initial density regimes (density at the half-mass radius $ρ_{\rm h}\ge{}3.4\times10^4$ and $\ge{1.5\times10^2}$ M$_\odot$ pc$^{-3}$ in dense and loose star clusters, respectively). Metal-poor clusters tend to form more massive BHs than metal-rich ones. We find $\sim{}6$, $\sim{}2$, and $<1$% of BHs with mass $m_{\rm BH}>60$ M$_\odot$ at Z=0.0002, 0.002 and 0.02, respectively. In metal-poor clusters, we form intermediate-mass BHs with mass up to $\sim{}320$ M$_\odot$. BBH mergers born via dynamical exchanges (exchanged BBHs) can be more massive than BBH mergers formed from binary evolution: the former (latter) reach total mass up to $\sim{}140$ M$_\odot$ ($\sim{}80$ M$_\odot$). The most massive BBH merger in our simulations has primary mass $\sim{}88$ M$_\odot$, inside the pair-instability mass gap, and a mass ratio of $\sim{}0.5$. Only BBHs born in young star clusters from metal-poor progenitors can match the masses of GW170729, the most massive event in O1 and O2, and those of GW190412, the first unequal-mass merger. We estimate a local BBH merger rate density $\sim{}110$ and $\sim{}55$ Gpc$^{-3}$ yr$^{-1}$, if we assume that all stars form in loose and dense star clusters, respectively.

preprint2020arXiv

Different Fates of Young Star Clusters After Gas Expulsion

We identify structures of the young star cluster NGC 2232 in the solar neighborhood (323.0 pc), and a newly discovered star cluster LP 2439 (289.1 pc). Member candidates are identified using the Gaia DR2 sky position, parallax and proper motion data, by an unsupervised machine learning method, \textsc{StarGO}. Member contamination from the Galactic disk is further removed using the color magnitude diagram. The four identified groups (NGC 2232, LP 2439 and two filamentary structures) of stars are coeval with an age of 25 Myr and were likely formed in the same giant molecular cloud. We correct the distance asymmetry from the parallax error with a Bayesian method. The 3D morphology shows the two spherical distributions of clusters NGC 2232 and LP 2439. Two filamentary structures are spatially and kinematically connected to NGC 2232. Both NGC 2232 and LP 2439 are expanding. The expansion is more significant in LP 2439, generating a loose spatial distribution with shallow volume number and mass density profiles. The expansion is suggested to be mainly driven by gas expulsion. NGC 2232, with 73~percent of the cluster mass bound, is currently experiencing a process of re-virialization, However, LP 2439, with 52 percent cluster mass being unbound, may fully dissolve in the near future. The different survivability traces different dynamical states of NGC 2232 and LP 2439 prior to the onset of gas expulsion. NGC 2232 may have been substructured and subvirial, while LP 2439 may either have been virial/supervirial, or it has experienced a much faster rate of gas removal.

preprint2020arXiv

Taking apart the dynamical clock. Fat-tailed dynamical kicks shape the blue-straggler star bimodality

In globular clusters, blue straggler stars are heavier than the average star, so dynamical friction strongly affects them. The radial distribution of BSS, normalized to a reference population, appears bimodal in a fraction of Galactic GCs, with a density peak in the core, a prominent zone of avoidance at intermediate radii, and again higher density in the outskirts. The zone of avoidance appears to be located at larger radii the more relaxed the host cluster, acting as a sort of dynamical clock. We use a new method to compute the evolution of the BSS radial distribution under dynamical friction and diffusion. We evolve our BSS in the mean cluster potential under dynamical friction plus a random fluctuating force, solving the Langevin equation with the Mannella quasi symplectic scheme. This amounts to a new simulation method which is much faster and simpler than direct N-body codes but retains their main feature: diffusion powered by strong, if infrequent, kicks. We compute the radial distribution of initially unsegregated BSS normalized to a reference population as a function of time. We trace the evolution of its minimum, corresponding to the zone of avoidance. We compare the evolution under kicks extracted from a Gaussian distribution to that obtained using a Holtsmark distribution. The latter is a fat tailed distribution which correctly models the effects of close gravitational encounters. We find that the zone of avoidance moves outwards over time, as expected based on observations, only when using the Holtsmark distribution. Thus the correct representation of near encounters is crucial to reproduce the dynamics of the system. We confirm and extend earlier results that showed how the dynamical clock indicator depends both on dynamical friction and effective diffusion powered by dynamical encounters.

preprint2019arXiv

Analytical solutions for the dynamical clock A+ indicator in a toy model of pure dynamical friction

Blue straggler stars are more massive than the average star in globular clusters, as they originate from the merger of two stars. Consequently, they experience dynamical friction, progressively sinking to the cluster center. Recently, several indicators of the degree of dynamical relaxation of a globular cluster have been proposed, based on the observed radial distribution of blue straggler stars. The most successful is the Alessandrini indicator, or A+ for short, which is the integral of the cumulative distribution of the blue straggler stars minus that of a lighter reference population. A+ correlates with the dynamical age of a cluster both in realistic simulations and in observations. Here I calculate the temporal dependence of the A+ indicator analytically in a simplified model of the evolution of the blue straggler star distribution under dynamical friction only.

preprint2019arXiv

Clustering clusters: unsupervised machine learning on globular cluster structural parameters

Globular Clusters (GCs) have historically been subdivided in either two (disk/halo) or three (disk/inner-halo/outer-halo) groups based on their orbital, chemical and internal physical properties. The qualitative nature of this subdivision makes it impossible to determine whether the natural number of groups is actually two, three, or more. In this paper we use cluster analysis on the $(\log M, \log σ_0, \log R_e, [Fe/H], \log | Z |)$ space to show that the intrinsic number of GC groups is actually either $k=2$ or $k=3$, with the latter being favored albeit non-significantly. In the $k=2$ case, the Partitioning Around Medoids (PAM) clustering algorithm recovers a metal-poor halo GC group and a metal-rich disk GC group. With $k=3$ the three groups can be interpreted as disk/inner-halo/outer-halo families. For each group we obtain a medoid, i.e. a representative element (NGC $6352$, NGC $5986$, and NGC $5466$ for the disk, inner halo, and outer halo respectively), and a measure of how strongly each GC is associated to its group, the so-called silhouette width. Using the latter, we find a correlation with age for both disk and outer halo GCs where the stronger the association of a GC with the disk (outer halo) group, the younger (older) it is.

preprint2019arXiv

Further properties of the dynamical clock A+ indicator in a toy model of pure dynamical friction

The Alessandrini A$+$ indicator is a quantitative measure of star cluster dynamical evolution based on the mass-segregation of blue straggler stars. A$+$ is defined as the integral of the cumulative distribution of blue straggler stars with the radius measured in logarithmic scale, minus a term related to the reference population used. In a companion paper I introduced a simplified model of dynamical friction and calculated the A$+$ indicator analytically. Here I show further properties of the time evolution of A$+$, focusing in particular on the physical interpretation of its time derivative $d{A^+}/dt$. I find that $d{A^+}/dt$ is the mean of the reciprocal dyamical friction timescale, weighted by the density of blue straggler stars. I show that it is non-negative (as expected based on monotonicity) due to the density of blue-straggler stars being non-negative and that, for a dynamical friction timescale that is non-decreasing with radius, $d{A^+}/dt$ is also non-decreasing with time, making the A$+$ indicator a convex function.

preprint2019arXiv

Mass and star formation rate of the host galaxies of compact binary mergers across cosmic time

We investigate the properties of the host galaxies of compact binary mergers across cosmic time, by means of population-synthesis simulations combined with galaxy catalogues from the EAGLE suite. We analyze the merger rate per galaxy of binary neutron stars (BNSs), black hole--neutron star binaries (BHNSs) and binary black holes (BBHs) from redshift zero up to redshift six. The binary merger rate per galaxy strongly correlates with the stellar mass of the host galaxy at any redshift considered here. This correlation is significantly steeper for BNSs than for both BHNSs and BBHs. Moreover, we find that the merger rate per galaxy depends also on host galaxy's star formation rate and metallicity. We derive a robust fitting formula that relates the merger rate per galaxy with galaxy's star formation rate, stellar mass and metallicity at different redshifts. The typical masses of the host galaxies increase significantly as redshift decreases, as a consequence of the interplay between delay time distribution of compact binaries and cosmic assembly of galaxies. Finally, we study the evolution of the merger rate density with redshift. At low redshift ($z\leq{}0.1$) early-type galaxies give a larger contribution to the merger rate density than late-type galaxies. This trend reverts at $z\ge{}1$.

preprint2019arXiv

Stellar-Mass Black Holes in Globular Clusters: Dynamical Consequences and Observational Signatures

Sizeable number of stellar-mass black holes (BHs) in globular clusters (GCs) can strongly influence the dynamical evolution and observational properties of their host cluster. Using results from a large set of numerical simulations, we identify the key ingredients needed to sustain a sizeable population of BHs in GCs up to a Hubble time. We find that while BH natal kick prescriptions are essential in determining the initial retention fraction of BHs in GCs, the long-term survival of BHs is determined by the size, initial central density and half-mass relaxation time of the GC. Simulated GC models that contain many BHs are characterized by relatively low central surface brightness, large half-light and core radii values. We also discuss novel ways to compare simulated results with available observational data to identify GCs that are most likely to contain many BHs.

preprint2016arXiv

Detecting intermediate mass black holes in globular clusters with machine learning

Mergers of stellar-mass black holes were recently observed in the gravitational wave window opened by LIGO. This puts the spotlight on dense stellar systems and their ability to create intermediate-mass black holes (IMBHs) through repeated merging. Unfortunately, attempts at direct and indirect IMBH detection in star clusters in the nearby universe have proven inconclusive as of now. Indirect detection methods attempt to constrain IMBHs through their effect on star cluster photometric and kinematic observables. They are usually based on looking for a specific, physically motivated signature. While this approach is justified, it may be suboptimal in its usage of the available data. Here I present a new indirect detection method, based on machine learning, that is unaffected by these restrictions. I reduce the scientific question whether a star cluster hosts an IMBH to a classification problem in the machine learning framework. I present preliminary results to illustrate how machine learning models are trained on simulated datasets and measure their performance on previously unseen, simulated data.

preprint2016arXiv

Globular clusters hosting intermediate-mass black-holes: no mass-segregation based candidates

Recently, both stellar mass-segregation and binary-fractions were uniformly measured on relatively large samples of Galactic Globular Clusters (GCs). Simulations show that both sizeable binary-star populations and Intermediate-Mass Black Holes (IMBHs) quench mass-segregation in relaxed GCs. Thus mass-segregation in GCs with a reliable binary-fraction measurement is a valuable probe to constrain IMBHs. In this paper we combine mass-segregation and binary-fraction measurements from the literature to build a sample of 33 GCs (with measured core-binary fractions), and a sample of 43 GCs (with a binary fraction measurement in the area between the core radius and the half-mass radius). Within both samples we try to identify IMBH-host candidates. These should have relatively low mass-segregation, a low binary fraction (< 5%), and short (< 1 Gyr) relaxation time. Considering the core binary fraction sample, no suitable candidates emerge. If the binary fraction between the core and the half-mass radius is considered, two candidates are found, but this is likely due to statistical fluctuations. We also consider a larger sample of 54 GCs where we obtained an estimate of the core binary fraction using a predictive relation based on metallicity and integrated absolute magnitude. Also in this case no suitable candidates are found. Finally, we consider the GC core- to half-mass radius ratio, that is expected to be larger for GCs containing either an IMBH or binaries. We find that GCs with large core- to half-mass radius ratios are less mass-segregated (and show a larger binary fraction), confirming the theoretical expectation that the energy sources responsible for the large core are also quenching mass-segregation

preprint2016arXiv

Merged or monolithic? Using machine-learning to reconstruct the dynamical history of simulated star clusters

Context. Machine-Learning (ML) solves problems by learning patterns from data, with limited or no human guidance. In Astronomy, it is mainly applied to large observational datasets, e.g. for morphological galaxy classification. Aims. We apply ML to gravitational N-body simulations of star clusters that are either formed by merging two progenitors or evolved in isolation, planning to later identify Globular Clusters (GCs) that may have a history of merging from observational data. Methods. We create mock-observations from simulated GCs, from which we measure a set of parameters (also called features in the machine-learning field). After dimensionality reduction on the feature space, the resulting datapoints are fed to various classification algorithms. Using repeated random subsampling validation we check whether the groups identified by the algorithms correspond to the underlying physical distinction between mergers and monolithically evolved simulations. Results. The three algorithms we considered (C5.0 trees, k-nearest neighbour, and support-vector machines) all achieve a test misclassification rate of about 10% without parameter tuning, with support-vector machines slightly outperforming the others. The first principal component of feature space correlates with cluster concentration. If we exclude it from the regression, the performance of the algorithms is only slightly reduced.

preprint2015arXiv

On the Use of the Number Count of Blue Horizontal--Branch Stars to Infer the Dominant Building Blocks of the Milky Way Halo

The formation of the Milky Way stellar halo is thought to be the result of merging and accretion of building blocks such as dwarf galaxies and massive globular clusters. Recently, Deason et al. (2015) suggested that the Milky Way outer halo formed mostly from big building blocks, such as dwarf spheroidal galaxies, based on the similar number ratio of blue straggler (BS) stars to blue horizontal-branch (BHB) stars. Here we demonstrate, however, that this result is seriously biased by not taking into detailed consideration on the formation mechanism of BHB stars from helium enhanced second-generation population. In particular, the high BS-to-BHB ratio observed in the outer halo fields is most likely due to a small number of BHB stars provided by GCs rather than to a large number of BS stars. This is supported by our dynamical evolution model of GCs which shows preferential removal of first generation stars in GCs. Moreover, there are a sufficient number of outer halo GCs which show very high BS-to-BHB ratio. Therefore, the BS-to-BHB number ratio is not a good indicator to use in arguing that more massive dwarf galaxies are the main building blocks of the Milky Way outer halo. Several lines of evidence still suggest that GCs can contribute a significant fraction of the outer halo stars.

preprint2013arXiv

Core collapse and horizontal-branch morphology in galactic globular clusters

Context. Stellar collision rates in globular clusters (GCs) do not appear to correlate with horizontal branch (HB) morphology, sug- gesting that dynamics does not play a role in the second-parameter problem. However, core densities and collision rates derived from surface-brightness may be significantly underestimated as the surface-brightness profile of GCs is not necessarily a good indicator of the dynamical state of GC cores. Core-collapse may go unnoticed if high central densities of dark remnants are present. Aims. We test whether GC HB morphology data supports a dynamical contribution to the so-called second-parameter effect. Methods. To remove first-parameter dependence we fitted the maximum effective temperature along the HB as a function of metal- licity in a sample of 54 Milky Way GCs. We plotted the residuals to the fit as a function of second-parameter candidates, namely dynamical age and total luminosity. We considered dynamical age (i.e. the ratio between age and half-light relaxation time) among possible second-parameters. We used a set of direct N-body simulations, including ones with dark remnants to illustrate how core density peaks, due to core collapse, in a dynamical-age range similar to that in which blue HBs are overabundant with respect to the metallicity expectation, especially for low-concentration initial conditions. Results. GC total luminosity shows nonlinear behavior compatible with the self-enrichment picture. However, the data are amenable to a different interpretation based on a dynamical origin of the second-parameter effect. Enhanced mass-stripping in the late red-giant- branch phase due to stellar interactions in collapsing cores is a viable candidate mechanism. In this picture, GCs with HBs bluer than expected based on metallicity are those undergoing core-collapse.

preprint2010arXiv

Croatian Black Hole School 2010 lecture notes on IMBHs in GCs

Black holes are fascinating objects. As a class of solutions to the Einstein equations they have been studied a great deal, yielding a wealth of theoretical results. But do they really exist? What do astronomers really mean when they claim to have observational evidence of their existence? To answer these questions, I will focus on a particular range of black-hole masses, approximately from 100 to 10000 solar masses. Black holes of this size are named Intermediate Mass Black Holes (IMBHs) and their existence is still heavily disputed, so they will be perfect for illustrating the observational challenges faced by a black hole hunter

preprint2010arXiv

On the Fundamental Line of Galactic and Extragalactic Globular Clusters

In a previous paper we found that the Globular Clusters of our Galaxy lie around a line in the log(Re), SBe, log(sigma) parameter space, with a moderate degree of scatter and remarkable axi-symmetry. This implies the existence of a purely photometric scaling law obtained by projecting such a line onto the log(Re), SBe plane. Such photometric quantities are readily available for large samples of clusters, as opposed to stellar velocity dispersion data. We study a sample of 129 Galactic and extragalactic clusters on such photometric plane in the V-band. We look for a linear relation between SBe and log(Re) and study how the scatter around it is influenced by age and dynamical environment. We interpret our results as a test on the evolutionary versus primordial origin of the Fundamental Line. We perform a detailed analysis of surface brightness profiles, which allows us to present a catalogue of structural properties, without relying on a given dynamical model. We find a linear relation between SBe and log(Re), in the form SBe = (5.25 +- 0.44) log(Re) + (15.58 +- 0.28), where SBe is measured in mag/arcsec^2 and Re in parsec. Both young and old clusters lie on the scaling law, with a scatter of approximately 1 mag in SBe. However, young clusters display more scatter and a clear trend of such scatter with age, which old clusters do not. Such trend becomes tighter if cluster age is measured in units of the cluster half-light relaxation time. Two-body relaxation therefore plays a major role, together with passive stellar population evolution, in shaping the relation between SBe, log(Re), and cluster age. We argue that the log(Re)-SBe relation and hence the Fundamental Line scaling law is not primordially set at cluster formation, but rather is the result of combined stellar evolution and collisional dynamical evolution.

preprint2009arXiv

Mass Segregation in NGC 2298: limits on the presence of an Intermediate Mass Black Hole

[abridged] Theoretical investigations have suggested the presence of Intermediate Mass Black Holes (IMBHs, with masses in the 100-10000 Msun range) in the cores of some Globular Clusters (GCs). In this paper we present the first application of a new technique to determine the presence or absence of a central IMBH in globular clusters that have reached energy equipartition via two-body relaxation. The method is based on the measurement of the radial profile for the average mass of stars in the system, using the fact that a quenching of mass segregation is expected when an IMBH is present. Here we measure the radial profile of mass segregation using main-sequence stars for the globular cluster NGC 2298 from resolved source photometry based on HST-ACS data. The observations are compared to expectations from direct N-body simulations of the dynamics of star clusters with and without an IMBH. The mass segregation profile for NGC 2298 is quantitatively matched to that inferred from simulations without a central massive object over all the radial range probed by the observations, that is from the center to about two half-mass radii. Profiles from simulations containing an IMBH more massive than ~ 300-500 Msun (depending on the assumed total mass of NGC 2298) are instead inconsistent with the data at about 3 sigma confidence, irrespective of the IMF and binary fraction chosen for these runs. While providing a null result in the quest of detecting a central black hole in globular clusters, the data-model comparison carried out here demonstrates the feasibility of the method which can also be applied to other globular clusters with resolved photometry in their cores.

preprint2009arXiv

Tidal disruption, global mass function and structural parameters evolution in star clusters

[abridged] We present a unified picture for the evolution of star clusters on the two-body relaxation timescale. We use direct N-body simulations of star clusters in a galactic tidal field starting from different multi-mass King models, up to 10% of primordial binaries and up to Ntot=65536 particles. An additional run also includes a central Intermediate Mass Black Hole. We find that for the broad range of initial conditions we have studied the stellar mass function of these systems presents a universal evolution which depends only on the fractional mass loss. The structure of the system, as measured by the core to half mass radius ratio, also evolves toward a universal state, which is set by the efficiency of heating on the visible population of stars induced by dynamical interactions in the core of the system. Interactions with dark remnants are dominant over the heating induced by a moderate population of primordial binaries (3-5%), especially under the assumption that most of the neutron stars and black holes are retained in the system. All our models without primordial binaries undergo a deep gravothermal collapse in the radial mass profile. However their projected light distribution can be well fitted by medium concentration King models (with parameter W0 ~ 8), even though there tends to be an excess over the best fit for the innermost points of the surface brightness. This excess is consistent with a shallow cusp in the surface brightness (mu(R) ~ R^{-v} with v ~ 0.4-0.7), like it has been observed for many globular clusters from high-resolution HST imaging. Classification of core-collapsed globular clusters based on their surface brightness profile is likely to fail in systems that have already bounced back to lower concentrations.