Showing 1–35 of 35 results for author: Zocchi, A

The GALAH survey: Tracing the Milky Way's formation and evolution through RR Lyrae stars

Authors: Valentina D'Orazi, Nicholas Storm, Andrew R. Casey, Vittorio F. Braga, Alice Zocchi, Giuseppe Bono, Michele Fabrizio, Christopher Sneden, Davide Massari, Riano E. Giribaldi, Maria Bergemann, Simon W. Campbell, Luca Casagrande, Richard de Grijs, Gayandhi De Silva, Maria Lugaro, Daniel B. Zucker, Angela Bragaglia, Diane Feuillet, Giuliana Fiorentino, Brian Chaboyer, Massimo Dall'Ora, Massimo Marengo, Clara E. Martínez-Vázquez, Noriyuki Matsunaga , et al. (17 additional authors not shown)

Abstract: Stellar mergers and accretion events have been crucial in shaping the evolution of the Milky Way (MW). These events have been dynamically identified and chemically characterised using red giants and main-sequence stars. RR Lyrae (RRL) variables can play a crucial role in tracing the early formation of the MW since they are ubiquitous, old (t$\ge$10 Gyr) low-mass stars and accurate distance indicat… ▽ More Stellar mergers and accretion events have been crucial in shaping the evolution of the Milky Way (MW). These events have been dynamically identified and chemically characterised using red giants and main-sequence stars. RR Lyrae (RRL) variables can play a crucial role in tracing the early formation of the MW since they are ubiquitous, old (t$\ge$10 Gyr) low-mass stars and accurate distance indicators. We exploited Data Release 3 of the GALAH survey to identify 78 field RRLs suitable for chemical analysis. Using synthetic spectra calculations, we determined atmospheric parameters and abundances of Fe, Mg, Ca, Y, and Ba. Most of our stars exhibit halo-like chemical compositions, with an iron peak around [Fe/H]$\approx -$1.40, and enhanced Ca and Mg content. Notably, we discovered a metal-rich tail, with [Fe/H] values ranging from $-$1 to approximately solar metallicity. This sub-group includes almost ~1/4 of the sample, it is characterised by thin disc kinematics and displays sub-solar $α$-element abundances, marginally consistent with the majority of the MW stars. Surprisingly, they differ distinctly from typical MW disc stars in terms of the s-process elements Y and Ba. We took advantage of similar data available in the literature and built a total sample of 535 field RRLs for which we estimated kinematical and dynamical properties. We found that metal-rich RRLs (1/3 of the sample) likely represent an old component of the MW thin disc. We also detected RRLs with retrograde orbits and provided preliminary associations with the Gaia-Sausage-Enceladus, Helmi, Sequoia, Sagittarius, and Thamnos stellar streams. △ Less

Submitted 7 May, 2024; originally announced May 2024.

Comments: Accepted for publication in MNRAS. 29 pages, 20 figures

Differential reddening in 48 globular clusters: An end to the quest for the intracluster medium

Authors: E. Pancino, A. Zocchi, M. Rainer, M. Monaci, D. Massari, M. Monelli, L. K. Hunt, L. Monaco, C. E. Martínez-Vázquez, N. Sanna, S. Bianchi, P. B. Stetson

Abstract: For decades, it has been theorized that a tenuous but detectable intracluster medium should be present in globular clusters, which is continuously replenished by the gas and dust ejected by bright giants and periodically cleared by interactions with the Galactic disk. However, dedicated searches, especially in infrared and radio wavelengths, have returned mostly upper limits, which are lower than… ▽ More For decades, it has been theorized that a tenuous but detectable intracluster medium should be present in globular clusters, which is continuously replenished by the gas and dust ejected by bright giants and periodically cleared by interactions with the Galactic disk. However, dedicated searches, especially in infrared and radio wavelengths, have returned mostly upper limits, which are lower than theoretical expectations by several orders of magnitude. We profited from recent wide-field photometry for 48 Galactic globular clusters to compute high-resolution maps of differential reddening, which can be used to correct any photometric catalog in these areas for reddening variations. Using 3D reddening maps from the literature, we evaluated the amount of foreground extinction. This allowed us to estimate the masses of the intracluster medium in our sample clusters, with an accuracy of one order of magnitude. Our estimates agree with the few available literature detections and with theoretical expectations. Because the discrepancy between observations and expectations only concerns literature upper limits, we explored possible reasons why they could be underestimated and we show that two recent discoveries can explain the discrepancy. The first is the recent discovery that the intracluster medium in 47 Tuc is not centrally concentrated. This is also supported by our maps, which in the majority of cases do not show a central reddening concentration. The second is the discovery that the dust in metal-poor ([Fe/H] less than about -1 dex) globular clusters is dominated by iron grains rather than silicates, which undermines previous dust mass estimates from observed upper limits. We conclude that current evidence, including our maps, does not contradict theoretical expectations and the problem of the missing intracluster medium is no longer an issue. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: Accepted by A&A, 20 pages, 14 figures

Investigating black hole accretion disks as potential polluter sources for the formation of enriched stars in globular clusters

Authors: Laurane Fréour, Alice Zocchi, Glenn van de Ven, Elena Pancino

Abstract: Accretion disks surrounding stellar mass black holes (BHs) have been suggested as potential locations for the nucleosynthesis of light elements, which are our primary observational discriminant of multiple stellar populations within globular clusters. The population of enriched stars in globular clusters are enhanced in N14, Na23, and sometimes in Al27 and/or in K39. In this study, our aim is to i… ▽ More Accretion disks surrounding stellar mass black holes (BHs) have been suggested as potential locations for the nucleosynthesis of light elements, which are our primary observational discriminant of multiple stellar populations within globular clusters. The population of enriched stars in globular clusters are enhanced in N14, Na23, and sometimes in Al27 and/or in K39. In this study, our aim is to investigate the feasibility of initiating nucleosynthesis for these four elements in BH accretion disks, considering various internal parameters such as the temperature of the gas and timescale of the accretion. To achieve this, we employed a 132-species reaction network. We used the slim disk model, suitable for the Super-Eddington mass accretion rate and for geometrically and optically thick disks. We explored the conditions related to the mass, mass accretion rate, viscosity, and radius of the BH-accretion disk system that would allow for the creation of N14, Na23, Al27, and K39 before the gas is accreted onto the central object. Our findings reveal that there is no region in the parameter space where the formation of Na23 can occur and only a very limited region where the formation of N14, Al27, and K39 is plausible. Specifically, this occurs for BHs with masses lower than 10 solar masses, with a preference toward even lower mass values and extremely low viscosity parameters ($α<10^{-3}$). Such values are highly unlikely based on current observations of stellar mass BHs. However, such low mass BHs could actually exist in the early universe, as so-called primordial BHs. In conclusion, our study suggests that the nucleosynthesis within BH accretion disks of four elements of interest for the multiple stellar populations is improbable, but not impossible, using the slim disk model. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: 11 pages, 7 figures, 2 tables, Accepted by A&A

HRMOS White Paper: Science Motivation

Authors: Laura Magrini, Thomas Bensby, Anna Brucalassi, Sofia Randich, Robin Jeffries, Gayandhi de Silva, Asa Skuladottir, Rodolfo Smiljanic, Oscar Gonzalez, Vanessa Hill, Nadege Lagarde, Eline Tolstoy, Jose' Maria Arroyo-Polonio, Martina Baratella, John R. Barnes, Giuseppina Battaglia, Holger Baumgardt, Michele Bellazzini, Katia Biazzo, Angela Bragaglia, Bradley Carter, Giada Casali, Gabriele Cescutti, Camilla Danielski, Elisa Delgado Mena , et al. (30 additional authors not shown)

Abstract: The High-Resolution Multi-Object Spectrograph (HRMOS) is a facility instrument that we plan to propose for the Very Large Telescope (VLT) of the European Southern Observatory (ESO), following the initial presentation at the VLT 2030 workshop held at ESO in June 2019. HRMOS provides a combination of capabilities that are essential to carry out breakthrough science across a broad range of active res… ▽ More The High-Resolution Multi-Object Spectrograph (HRMOS) is a facility instrument that we plan to propose for the Very Large Telescope (VLT) of the European Southern Observatory (ESO), following the initial presentation at the VLT 2030 workshop held at ESO in June 2019. HRMOS provides a combination of capabilities that are essential to carry out breakthrough science across a broad range of active research areas from stellar astrophysics and exoplanet studies to Galactic and Local Group archaeology. HRMOS fills a gap in capabilities amongst the landscape of future instrumentation planned for the next decade. The key characteristics of HRMOS will be high spectral resolution (R = 60000 - 80000) combined with multi-object (20-100) capabilities and long term stability that will provide excellent radial velocity precision and accuracy (10m/s). Initial designs predict that a SNR~100 will be achievable in about one hour for a star with mag(AB) = 15, while with the same exposure time a SNR~ 30 will be reached for a star with mag(AB) = 17. The combination of high resolution and multiplexing with wavelength coverage extending to relatively blue wavelengths (down to 380\,nm), makes HRMOS a spectrograph that will push the boundaries of our knowledge and that is envisioned as a workhorse instrument in the future. The science cases presented in this White Paper include topics and ideas developed by the Core Science Team with the contributions from the astronomical community, also through the wide participation in the first HRMOS Workshop (https://indico.ict.inaf.it/event/1547/) that took place in Firenze (Italy) in October 2021. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: 88 pages, 39 figures. Comments and expressions of interest are welcome by contacting members of the Core Science Team

Total mass slopes and enclosed mass constrained by globular cluster system dynamics

Authors: Tadeja Veršič, Sabine Thater, Glenn van de Ven, Laura L. Watkins, Prashin Jethwa, Ryan Leaman, Alice Zocchi

Abstract: The goal of this work is to probe the total mass distribution of early-type galaxies with globular clusters (GCs) as kinematic tracers, by constraining the parameters of the profile with a flexible modelling approach. To that end, we leverage the extended spatial distribution of GCs from the SLUGGS survey ($\langle R_{\rm GC,\ max} \rangle \sim 8R_{\rm e}$) in combination with discrete dynamical m… ▽ More The goal of this work is to probe the total mass distribution of early-type galaxies with globular clusters (GCs) as kinematic tracers, by constraining the parameters of the profile with a flexible modelling approach. To that end, we leverage the extended spatial distribution of GCs from the SLUGGS survey ($\langle R_{\rm GC,\ max} \rangle \sim 8R_{\rm e}$) in combination with discrete dynamical modelling. We use discrete Jeans anisotropic modelling in cylindrical coordinates to determine the velocity moments at the location of the GCs in our sample. We use a Bayesian framework to determine the best-fit parameters of the total mass density profile and orbital properties of the GC systems. We find that the orbital properties (anisotropy and rotation of the dispersion-dominated GC systems) minimally impact the measurements of the inner slope and enclosed mass, while a strong presence of dynamically-distinct subpopulations or low numbers of kinematic tracers can bias the results. Owing to the large spatial extent of the tracers our method is sensitive to the intrinsic inner slope of the total mass profile and we find $\overlineα = -1.88\pm 0.01$ for 12 galaxies with robust measurements. To compare our results with literature values we fit a single power-law profile to the resulting total mass density. In the radial range 0.1-4~$R_{\rm e}$ our measured slope has a value of $\langle γ_{\rm tot}\rangle = -2.22\pm0.14$ and is in good agreement with the literature. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 17 pages, 13 Figures, 8 Tables, Accepted for publication in A&A

Journal ref: A&A 681, A46 (2024)

Schwarzschild Modeling of Barred S0 Galaxy NGC 4371

Authors: Behzad Tahmasebzadeh, Ling Zhu, Juntai Shen, Dimitri A. Gadotti, Monica Valluri, Sabine Thater, Glenn van de Ven, Yunpeng Jin, Ortwin Gerhard, Peter Erwin, Prashin Jethwa, Alice Zocchi, Edward J. Lilley, Francesca Fragkoudi, Adriana de Lorenzo-Cáceres, Jairo Méndez-Abreu, Justus Neumann, Rui Guo

Abstract: We apply the barred Schwarzschild method developed by Tahmasebzadeh et al. (2022) to a barred S0 galaxy, NGC 4371, observed by IFU instruments from the TIMER and ATLAS3D projects. We construct the gravitational potential by combining a fixed black hole mass, a spherical dark matter halo, and stellar mass distribution deprojected from $3.6$ $μ$m S$^4$G image considering an axisymmetric disk and a t… ▽ More We apply the barred Schwarzschild method developed by Tahmasebzadeh et al. (2022) to a barred S0 galaxy, NGC 4371, observed by IFU instruments from the TIMER and ATLAS3D projects. We construct the gravitational potential by combining a fixed black hole mass, a spherical dark matter halo, and stellar mass distribution deprojected from $3.6$ $μ$m S$^4$G image considering an axisymmetric disk and a triaxial bar. We independently modeled kinematic data from TIMER and ATLAS3D. Both models fit the data remarkably well. We find a consistent bar pattern speed from the two sets of models, with $Ω_{\rm p} = 23.6 \pm 2.8 \hspace{.08cm} \mathrm{km \hspace{.04cm} s^{-1} \hspace{.04cm} kpc^{-1} }$. The dimensionless bar rotation parameter is determined to be $ R_{\rm cor}/R_{\rm bar}=1.9 \pm 0.3$, indicating a slow bar in NGC 4371. Additionally, our model predicts a high amount of dark matter within the bar region ($M_{\rm DM}/ M_{\rm total}$ $\sim 0.51 \pm 0.06$), which, aligned with the predictions of cosmological simulations, indicates that fast bars are generally found in baryon-dominated disks. Based on the best-fitting model, we further decompose the galaxy into multiple 3D orbital structures, including a BP/X bar, a classical bulge, a nuclear disk, and a main disk. The BP/X bar is not perfectly included in the input 3D density model, but BP/X-supporting orbits are picked through the fitting to the kinematic data. This is the first time a real barred galaxy has been modeled utilizing the Schwarzschild method including a 3D bar. △ Less

Submitted 20 May, 2024; v1 submitted 30 September, 2023; originally announced October 2023.

Comments: Resubmitted after addressing the referee's reports

Oxygen, sulfur, and iron radial abundance gradients of classical Cepheids across the Galactic thin disk

Authors: R. da Silva, V. D'Orazi, M. Palla, G. Bono, V. F. Braga, M. Fabrizio, B. Lemasle, E. Spitoni, F. Matteucci, H. Jonsson, V. Kovtyukh, L. Magrini, M. Bergemann, M. Dall'Ora, I. Ferraro, G. Fiorentino, P. Francois, G. Iannicola, L. Inno, R. -P. Kudritzki, N. Matsunaga, M. Monelli, M. Nonino, C. Sneden, J. Storm , et al. (3 additional authors not shown)

Abstract: Classical Cepheids (CCs) are solid distance indicators and tracers of young stellar populations. Our aim is to provide iron, oxygen, and sulfur abundances for the largest and most homogeneous sample of Galactic CCs ever analyzed. The current sample covers a wide range in Galactocentric distances (RG), pulsation modes and periods. High-resolution and high S/N spectra collected with different spectr… ▽ More Classical Cepheids (CCs) are solid distance indicators and tracers of young stellar populations. Our aim is to provide iron, oxygen, and sulfur abundances for the largest and most homogeneous sample of Galactic CCs ever analyzed. The current sample covers a wide range in Galactocentric distances (RG), pulsation modes and periods. High-resolution and high S/N spectra collected with different spectrographs were adopted to estimate the atmospheric parameters. Individual distances are based on Gaia trigonometric parallaxes or on near-infrared Period-Luminosity relations. We found that Fe and alpha-element radial gradients based on CCs display a well-defined change in the slope for RG larger than 12 kpc. Radial gradients based on open clusters, covering a wide range in age, display similar trends, meaning that the flattening in the outer disk is an intrinsic feature of the radial gradients since it is independent of age. Empirical evidence indicates that the radial gradient for S is steeper than for Fe. The difference in the slope is a factor of two in the linear fit. We also found that S is, on average, under-abundant compared with O. We performed a detailed comparison with Galactic chemical evolution models and we found that a constant Star Formation Efficiency for RG larger than 12 kpc takes account for the flattening in both Fe and alpha-elements. To further constrain the impact that predicted S yields for massive stars have on radial gradients, we adopted a "toy model" and we found that the flattening in the outermost regions requires a decrease of a factor of four in the current S predictions. Sulfur photospheric abundances, compared with other alpha-elements, have the key advantage of being a volatile element. Therefore, stellar S abundances can be directly compared with nebular S abundances in external galaxies. △ Less

Submitted 2 August, 2023; originally announced August 2023.

Comments: 21 pages, 13 figures, 3 tables; to be published in the A&A journal

Journal ref: A&A 678, A195 (2023)

Dynamical modelling of ATLAS$^{\rm 3D}$ galaxies

Authors: Sabine Thater, Prashin Jethwa, Edward J. Lilley, Alice Zocchi, Giulia Santucci, Glenn van de Ven

Abstract: Triaxial dynamical models of massive galaxies observed in the ATLAS3D project can provide new insights into the complex evolutionary processes that shape galaxies. The ATLAS3D survey is ideal as the sample comprises a good mix of fast and slow rotators with vastly different mass assembly histories. We present a detailed dynamical study with our triaxial modelling code DYNAMITE, which models galaxi… ▽ More Triaxial dynamical models of massive galaxies observed in the ATLAS3D project can provide new insights into the complex evolutionary processes that shape galaxies. The ATLAS3D survey is ideal as the sample comprises a good mix of fast and slow rotators with vastly different mass assembly histories. We present a detailed dynamical study with our triaxial modelling code DYNAMITE, which models galaxies as a superposition of their stellar orbits. The models allow us to constrain the intrinsic shape of the stellar component, the distributions of the visible and invisible matter and the orbit distribution in these nearby early-type galaxies and to relate it with different evolutionary scenarios. Triaxial modelling is essential for these galaxies to understand their complex kinematical features. △ Less

Submitted 16 May, 2023; originally announced May 2023.

Comments: 4 pages, 3 figures, Proceeding of IAU Symposium 379: Dynamical Masses of Local Group galaxies, ed. P. Bonifacio, M.-R. Cioni, F. Hammer, M. Pawlowski, and S. Taibi

A Wide View of the Galactic Globular Cluster NGC 2808: Red Giant and Horizontal Branch Star Spatial Distributions

Authors: Christian I. Johnson, Annalisa Calamida, Justin A. Kader, Ivan Ferraro, Catherine A. Pilachowski, Giuseppe Bono, Alessandra Mastrobuono-Battisti, Armin Rest, Alfredo Zenteno, Alice Zocchi

Abstract: Wide-field and deep DECam multi-band photometry, combined with HST data for the core of the Galactic globular cluster NGC 2808, allowed us to study the distribution of various stellar sub-populations and stars in different evolutionary phases out to the cluster tidal radius. We used the C_ugi = (u-g)-(g-i) index to identify three chemically distinct sub-populations along the red giant branch and c… ▽ More Wide-field and deep DECam multi-band photometry, combined with HST data for the core of the Galactic globular cluster NGC 2808, allowed us to study the distribution of various stellar sub-populations and stars in different evolutionary phases out to the cluster tidal radius. We used the C_ugi = (u-g)-(g-i) index to identify three chemically distinct sub-populations along the red giant branch and compared their spatial distributions. The most light-element enriched sub-population (P3) is more centrally concentrated; however, it shows a more extended distribution in the external regions of the cluster compared to the primordial (P1) and intermediate (P2) composition populations. Furthermore, the P3 sub-population centroid is off-center relative to those of the P1 and P2 groups. We also analyzed the spatial distribution of horizontal branch stars and found that the relative fraction of red horizontal branch stars increases for radial distances larger than ~ 1.5' while that of the blue and hotter stars decreases. These new observations, combined with literature spectroscopic measurements, suggest that the red horizontal branch stars are the progeny of all the stellar sub-populations in NGC 2808, i.e. primordial and light-element enhanced, while the blue stars are possibly the result of a combination of the "hot-flasher" and the "helium-enhanced" scenarios. A similar distribution of different red giant branch sub-populations and horizontal branch stars was also found for the most massive Galactic globular cluster, omega Cen, based on combined DECam and HST data, which suggests the two may share a similar origin. △ Less

Submitted 10 May, 2023; originally announced May 2023.

Comments: Accepted for publication in AJ; 22 pages; 16 figures; the full version of table 2 will be available in electronic form with the published version

On the robustness of triaxial Schwarzschild modelling: The effects of correcting the orbit mirroring

Authors: Sabine Thater, Prashin Jethwa, Behzad Tahmasebzadeh, Ling Zhu, Mark den Brok, Giulia Santucci, Yuchen Ding, Adriano Poci, Edward Lilley, P. Tim de Zeeuw, Alice Zocchi, Thomas I. Maindl, Fabio Rigamonti, Glenn van de Ven, Meng Yang, Katja Fahrion

Abstract: In the past 15 years, the triaxial Schwarzschild orbit-superposition code by van den Bosch et al. (2008) has been widely applied to study the dynamics of galaxies. Recently, Quenneville et al. (2022) reported a bug in the orbit calculation of this code, specifically in the mirroring procedure that is used to speed up the computation. We have fixed the incorrect mirroring in DYNAMITE, which is the… ▽ More In the past 15 years, the triaxial Schwarzschild orbit-superposition code by van den Bosch et al. (2008) has been widely applied to study the dynamics of galaxies. Recently, Quenneville et al. (2022) reported a bug in the orbit calculation of this code, specifically in the mirroring procedure that is used to speed up the computation. We have fixed the incorrect mirroring in DYNAMITE, which is the publicly-released successor of the triaxial Schwarzschild code by van den Bosch et al. (2008). In this study, we provide a thorough quantification of how this bug has affected the results of dynamical analyses performed with this code. We compare results obtained with the original and corrected versions of DYNAMITE, and discuss the differences in the phase-space distribution of a single orbit and in the global stellar orbit distribution, in the mass estimate of the central black hole in the highly triaxial galaxy PGC 46832, and in the measurement of intrinsic shape and enclosed mass for more than 50 galaxies. Focusing on the typical scientific applications of a Schwarzschild triaxial code, in all our tests we find that differences are negligible with respect to the statistical and systematic uncertainties. We conclude that previous results with the van den Bosch et al. (2008) triaxial Schwarzschild code are not significantly affected by the incorrect mirroring. △ Less

Submitted 9 May, 2022; originally announced May 2022.

Comments: 14 pages, 11 figures. Submitted to A&A

Journal ref: A&A 667, A51 (2022)

The proper motion of sub-populations in Omega Centauri

Authors: N. Sanna, E. Pancino, A. Zocchi, F. R. Ferraro, P. B. Stetson

Abstract: The galactic globular cluster Omega Centauri is the most massive of its kind, with a complex mix of multiple stellar populations and several kinematic and dynamical peculiarities. Different mean proper motions have been detected among the three main sub-populations, implying that the most metal-rich one is of accreted origin. This particular piece of evidence has been a matter of debate because th… ▽ More The galactic globular cluster Omega Centauri is the most massive of its kind, with a complex mix of multiple stellar populations and several kinematic and dynamical peculiarities. Different mean proper motions have been detected among the three main sub-populations, implying that the most metal-rich one is of accreted origin. This particular piece of evidence has been a matter of debate because the available data have either not been sufficiently precise or limited to a small region of the cluster to ultimately confirm or refute the result. Using astrometry from the second Gaia data release and recent high-quality, multi-band photometry, we are now in a position to resolve the controversy. We reproduced the original analysis using the Gaia data and found that the three populations have the same mean proper motion. Thus, there is no need to invoke an accreted origin for the most metal-rich sub-population. △ Less

Submitted 30 March, 2020; originally announced March 2020.

Comments: 6 pages, 2 tables, 5 figures. Accepted for publication on A&A

Journal ref: A&A 637, A46 (2020)

The not so simple stellar system Omega Cen. II. Evidence in support of a merging scenario

Authors: Annalisa Calamida, Alice Zocchi, Giuseppe Bono, Ivan Ferraro, Alessandra Mastrobuono-Battisti, Abhijit Saha, Giacinto Iannicola, Armin Rest, Giovanni Strampelli, Alfredo Zenteno

Abstract: We present multi-band photometry covering $\sim$ 5°$\times$ 5°across $ω$ Cen collected with the Dark Energy Camera, combined to Hubble Space Telescope and Wide Field Imager data for the central regions. The unprecedented photometric accuracy and field coverage allowed us to confirm the different spatial distribution of blue and red main-sequence stars, and of red-giant branch (RGB) stars with diff… ▽ More We present multi-band photometry covering $\sim$ 5°$\times$ 5°across $ω$ Cen collected with the Dark Energy Camera, combined to Hubble Space Telescope and Wide Field Imager data for the central regions. The unprecedented photometric accuracy and field coverage allowed us to confirm the different spatial distribution of blue and red main-sequence stars, and of red-giant branch (RGB) stars with different metallicities. The ratio of the number of blue to red main-sequence stars shows that the blue main-sequence sub-population has a more extended spatial distribution compared to the red main-sequence one, and the frequency of blue main-sequence stars increases at a distance of $\sim$ 20 arcmin from $ω$ Cen center. Similarly, the more metal-rich RGB stars show a more extended spatial distribution compared to the more metal-poor ones in the outskirts of the cluster. Moreover, the centers of the distributions of metal-rich and metal-poor RGB stars are shifted in different directions with respect to the geometrical center of $ω$ Cen. We constructed stellar density profiles for the blue and red main-sequence stars; they confirm that the blue main-sequence sub-population has a more extended spatial distribution compared to the red main-sequence one in the outskirts of $ω$ Cen, as found based on the star number ratio. We also computed the ellipticity profile of $ω$ Cen, which has a maximum value of 0.16 at a distance of $\sim$ 8 arcmin from the center, and a minimum of 0.05 at $\sim$ 30 arcmin; the average ellipticity is $\sim0.10$. The circumstantial evidence presented in this work suggests a merging scenario for the formation of the peculiar stellar system $ω$ Cen. △ Less

Submitted 17 February, 2020; originally announced February 2020.

Comments: 24 pages, 17 figures, accepted for publication in The Astrophysical Journal

Spherical models of star clusters with potential escapers

Authors: Ian Claydon, Mark Gieles, Anna Lisa Varri, Douglas C. Heggie, Alice Zocchi

Abstract: An increasing number of observations of the outer regions of globular clusters (GCs) have shown a flattening of the velocity dispersion profile and an extended surface density profile. Formation scenarios of GCs can lead to different explanations of these peculiarities, therefore the dynamics of stars in the outskirts of GCs are an important tool in tracing back the evolutionary history and format… ▽ More An increasing number of observations of the outer regions of globular clusters (GCs) have shown a flattening of the velocity dispersion profile and an extended surface density profile. Formation scenarios of GCs can lead to different explanations of these peculiarities, therefore the dynamics of stars in the outskirts of GCs are an important tool in tracing back the evolutionary history and formation of star clusters. One possible explanation for these features is that GCs are embedded in dark matter halos. Alternatively, these features are the result of a population of energetically unbound stars that can be spatially trapped within the cluster, known as potential escapers (PEs). We present a prescription for the contribution of these energetically unbound members to a family of self-consistent, distribution function-based models, which, for brevity, we call the Spherical Potential Escapers Stitched (SPES) models. We show that, when fitting to mock data of bound and unbound stars from an N-body model of a tidally-limited star cluster, the SPES models correctly reproduce the density and velocity dispersion profiles up to the Jacobi radius, and they are able to recover the value of the Jacobi radius itself to within 20%. We also provide a comparison to the number density and velocity dispersion profiles of the Galactic cluster 47 Tucanae. Such a case offers a proof of concept that an appropriate modelling of PEs is essential to accurately interpret Gaia data in the outskirts of GCs, and, in turn, to formulate meaningful present-day constraints for GC formation scenarios in the early universe. △ Less

Submitted 15 March, 2019; v1 submitted 14 March, 2019; originally announced March 2019.

Comments: 16 pages, 11 figures, submitted to MNRAS

Homogeneous photometry VII. Globular clusters in the Gaia era

Authors: P. B. Stetson, E. Pancino, A. Zocchi, N. Sanna, M. Monelli

Abstract: We present wide-field, ground-based Johnson-Cousins UBVRI photometry for 48 Galactic globular clusters based on almost 90000 public and proprietary images. The photometry is calibrated with the latest transformations obtained in the framework of our secondary standard project, with typical internal and external uncertainties of order a few millimagnitudes. These data provide a bridge between exist… ▽ More We present wide-field, ground-based Johnson-Cousins UBVRI photometry for 48 Galactic globular clusters based on almost 90000 public and proprietary images. The photometry is calibrated with the latest transformations obtained in the framework of our secondary standard project, with typical internal and external uncertainties of order a few millimagnitudes. These data provide a bridge between existing small-area, high-precision HST photometry and all sky-catalogues from large surveys like Gaia, SDSS, or LSST. For many clusters, we present the first publicly available photometry in some of the five bands (typically U and R). We illustrate the scientific potential of the photometry with examples of surface density and brightness profiles and of colour-magnitude diagrams, with the following highlights: (i) we study the morphology of NGC 5904, finding a varying ellipticity and position angle as a function of radial distance; (ii) we show U-based colour-magnitude diagrams and demonstrate that no cluster in our sample is free from multiple stellar populations, with the possible exception of a few clusters with high and differential reddening or field contamination, for which more sophisticated investigations are required. This is true even for NGC 5694 and Terzan 8, that were previously considered as (mostly) single-population candidates. △ Less

Submitted 26 February, 2019; originally announced February 2019.

Comments: 24 pages, 12 figures, accepted for publication by MNRAS

A MODEST review

Authors: Anna Lisa Varri, Maxwell Xu Cai, Francisca Concha-Ramírez, František Dinnbier, Nora Lützgendorf, Václav Pavlík, Sara Rastello, Antonio Sollima, Long Wang, Alice Zocchi

Abstract: We present an account of the state of the art in the fields explored by the research community invested in 'Modeling and Observing DEnse STellar systems'. For this purpose, we take as a basis the activities of the MODEST-17 conference, which was held at Charles University, Prague, in September 2017. Reviewed topics include recent advances in fundamental stellar dynamics, numerical methods for the… ▽ More We present an account of the state of the art in the fields explored by the research community invested in 'Modeling and Observing DEnse STellar systems'. For this purpose, we take as a basis the activities of the MODEST-17 conference, which was held at Charles University, Prague, in September 2017. Reviewed topics include recent advances in fundamental stellar dynamics, numerical methods for the solution of the gravitational N-body problem, formation and evolution of young and old star clusters and galactic nuclei, their elusive stellar populations, planetary systems, and exotic compact objects, with timely attention to black holes of different classes of mass and their role as sources of gravitational waves. Such a breadth of topics reflects the growing role played by collisional stellar dynamics in numerous areas of modern astrophysics. Indeed, in the next decade, many revolutionary instruments will enable the derivation of positions and velocities of individual stars in the Milky Way and its satellites and will detect signals from a range of astrophysical sources in different portions of the electromagnetic and gravitational spectrum, with an unprecedented sensitivity. On the one hand, this wealth of data will allow us to address a number of long-standing open questions in star cluster studies; on the other hand, many unexpected properties of these systems will come to light, stimulating further progress of our understanding of their formation and evolution. △ Less

Submitted 17 October, 2018; originally announced October 2018.

Comments: 42 pages; accepted for publication in 'Computational Astrophysics and Cosmology'. We are much grateful to the organisers of the MODEST-17 conference (Charles University, Prague, September 2017). We acknowledge the input provided by all MODEST-17 participants, and, more generally, by the members of the MODEST community

Mass modelling globular clusters in the Gaia era: a method comparison using mock data from an $N$-body simulation of M4

Authors: Vincent Hénault-Brunet, Mark Gieles, Antonio Sollima, Laura L. Watkins, Alice Zocchi, Ian Claydon, Elena Pancino, Holger Baumgardt

Abstract: As we enter a golden age for studies of internal kinematics and dynamics of Galactic globular clusters (GCs), it is timely to assess the performance of modelling techniques in recovering the mass, mass profile, and other dynamical properties of GCs. Here, we compare different mass-modelling techniques (distribution-function (DF)-based models, Jeans models, and a grid of N-body models) by applying… ▽ More As we enter a golden age for studies of internal kinematics and dynamics of Galactic globular clusters (GCs), it is timely to assess the performance of modelling techniques in recovering the mass, mass profile, and other dynamical properties of GCs. Here, we compare different mass-modelling techniques (distribution-function (DF)-based models, Jeans models, and a grid of N-body models) by applying them to mock observations from a star-by-star N-body simulation of the GC M 4 by Heggie. The mocks mimic existing and anticipated data for GCs: surface brightness or number density profiles, local stellar mass functions, line-of-sight velocities, and Hubble Space Telescope- and Gaia-like proper motions. We discuss the successes and limitations of the methods. We find that multimass DF-based models, Jeans, and N-body models provide more accurate mass profiles compared to single-mass DF-based models. We highlight complications in fitting the kinematics in the outskirts due to energetically unbound stars associated with the cluster ("potential escapers", not captured by truncated DF models nor by N-body models of clusters in isolation), which can be avoided with DF-based models including potential escapers, or with Jeans models. We discuss ways to account for mass segregation. For example, three-component DF-based models with freedom in their mass function are a simple alternative to avoid the biases of single-mass models (which systematically underestimate the total mass, half-mass radius, and central density), while more realistic multimass DF-based models with freedom in the remnant content represent a promising avenue to infer the total mass and the mass function of remnants. △ Less

Submitted 22 November, 2018; v1 submitted 11 October, 2018; originally announced October 2018.

Comments: 27 pages, 12 figures, MNRAS accepted

The effect of stellar-mass black holes on the central kinematics of omega Cen: a cautionary tale for IMBH interpretations

Authors: Alice Zocchi, Mark Gieles, Vincent Hénault-Brunet

Abstract: The search for intermediate-mass black holes (IMBHs) in the centre of globular clusters is often based on the observation of a central cusp in the surface brightness profile and a rise towards the centre in the velocity dispersion profiles. Similar signatures, however, could result from other effects, that need to be taken into account in order to determine the presence (or the absence) of an IMBH… ▽ More The search for intermediate-mass black holes (IMBHs) in the centre of globular clusters is often based on the observation of a central cusp in the surface brightness profile and a rise towards the centre in the velocity dispersion profiles. Similar signatures, however, could result from other effects, that need to be taken into account in order to determine the presence (or the absence) of an IMBH in these stellar systems. Following our previous exploration of the role of radial anisotropy in shaping these observational signatures, we analyse here the effects produced by the presence of a population of centrally concentrated stellar-mass black holes. We fit dynamical models to omega Cen data, and we show that models with ~5% of their mass in black holes (consistent with ~100% retention fraction after natal kicks) can reproduce the data. When simultaneously considering both radial anisotropy and mass segregation, the best-fit model includes a smaller population of remnants, and a less extreme degree of anisotropy with respect to the models that include only one of these features. These results underline that before conclusions about putative IMBHs can be made, the effects of stellar-mass black holes and radial anisotropy need to be properly accounted for. △ Less

Submitted 6 June, 2018; originally announced June 2018.

Comments: 13 pages, 5 figures. Accepted for publication in MNRAS

Concurrent formation of supermassive stars and globular clusters: implications for early self-enrichment

Authors: Mark Gieles, Corinne Charbonnel, Martin Krause, Vincent Henault-Brunet, Oscar Agertz, Henny Lamers, Nathan Bastian, Alessia Gualandris, Alice Zocchi, James Petts

Abstract: We present a model for the concurrent formation of globular clusters (GCs) and supermassive stars (SMSs, $>10^3\,{\rm M}_\odot$) to address the origin of the HeCNONaMgAl abundance anomalies in GCs. GCs form in converging gas flows and accumulate low-angular momentum gas, which accretes onto protostars. This leads to an adiabatic contraction of the cluster and an increase of the stellar collision r… ▽ More We present a model for the concurrent formation of globular clusters (GCs) and supermassive stars (SMSs, $>10^3\,{\rm M}_\odot$) to address the origin of the HeCNONaMgAl abundance anomalies in GCs. GCs form in converging gas flows and accumulate low-angular momentum gas, which accretes onto protostars. This leads to an adiabatic contraction of the cluster and an increase of the stellar collision rate. A SMS can form via runaway collisions if the cluster reaches sufficiently high density before two-body relaxation halts the contraction. This condition is met if the number of stars $\gtrsim10^6$ and the gas accretion rate $\gtrsim10^5\,{\rm M}_\odot$/Myr, reminiscent of GC formation in high gas-density environments, such as -- but not restricted to -- the early Universe. The strong SMS wind mixes with the inflowing pristine gas, such that the protostars accrete diluted hot-hydrogen burning yields of the SMS. Because of continuous rejuvenation, the amount of processed material liberated by the SMS can be an order of magnitude higher than its maximum mass. This `conveyor-belt' production of hot-hydrogen burning products provides a solution to the mass budget problem that plagues other scenarios. Additionally, the liberated material is mildly enriched in helium and relatively rich in other hot-hydrogen burning products, in agreement with abundances of GCs today. Finally, we find a super-linear scaling between the amount of processed material and cluster mass, providing an explanation for the observed increase of the fraction of processed material with GC mass. We discuss open questions of this new GC enrichment scenario and propose observational tests. △ Less

Submitted 24 April, 2018; v1 submitted 12 April, 2018; originally announced April 2018.

Comments: 19 pages, 5 figures, accepted to MNRAS (shortened abstract and included feedback from the community)

The Peculiar Radial Distribution of Multiple Populations in the massive globular cluster M80

Authors: E. Dalessandro, M. Cadelano, E. Vesperini, M. Salaris, F. R. Ferraro, B. Lanzoni, S. Raso, J. Hong, J. J. Webb, A. Zocchi

Abstract: We present a detailed analysis of the radial distribution of light-element multiple populations (LE-MPs) in the massive and dense globular cluster M80 based on the combination of UV and optical Hubble Space Telescope data. Surprisingly, we find that first generation stars (FG) are significantly more centrally concentrated than extreme second generation ones (SG) out to $\sim 2.5 r_h$ from the clus… ▽ More We present a detailed analysis of the radial distribution of light-element multiple populations (LE-MPs) in the massive and dense globular cluster M80 based on the combination of UV and optical Hubble Space Telescope data. Surprisingly, we find that first generation stars (FG) are significantly more centrally concentrated than extreme second generation ones (SG) out to $\sim 2.5 r_h$ from the cluster center. To understand the origin of such a peculiar behavior, we used a set of $N$-body simulations following the long-term dynamical evolution of LE-MPs. We find that, given the advanced dynamical state of the cluster, the observed difference does not depend on the primordial relative distributions of FG and SG stars. On the contrary, a difference of $\sim 0.05-0.10 M_{\odot}$ between the average masses of the two sub-populations is needed to account for the observed radial distributions. We argue that such a mass difference might be the result of the higher He abundance of SG stars (of the order of $ΔY\sim 0.05-0.06$) with respect to FG. Interestingly, we find that a similar He variation is necessary to reproduce the horizontal branch morphology of M80. These results demonstrate that differences in mass among LE-MPs, due to different He content, should be properly taken into account for a correct interpretation of their radial distribution, at least in dynamically evolved systems. △ Less

Submitted 9 April, 2018; originally announced April 2018.

Comments: 19 pages, 9 figures, 1 table; Accepted for publication in ApJ

Three candidate double clusters in the LMC: truth or dare?

Authors: E. Dalessandro, A. Zocchi, A. L. Varri, A. Mucciarelli, M. Bellazzini, F. R. Ferraro, B. Lanzoni, E. Lapenna, L. Origlia

Abstract: The Large Magellanic Cloud (LMC) hosts a large number of candidate stellar cluster pairs. Binary stellar clusters provide important clues about cluster formation processes and the evolutionary history of the host galaxy. However, to properly extract and interpret this information, it is crucial to fully constrain the fraction of real binary systems and their physical properties. Here we present a… ▽ More The Large Magellanic Cloud (LMC) hosts a large number of candidate stellar cluster pairs. Binary stellar clusters provide important clues about cluster formation processes and the evolutionary history of the host galaxy. However, to properly extract and interpret this information, it is crucial to fully constrain the fraction of real binary systems and their physical properties. Here we present a detailed photometric analysis based on ESO-FORS2 images of three candidate cluster multiplets in the LMC, namely SL349-SL353, SL385-SL387-NGC1922 and NGC1836-BRHT4b-NGC1839. For each cluster we derived ages, structural parameters and morphological properties. We have also estimated the degree of filling of their Roche lobe, as an approximate tool to measure the strength of the tidal perturbations induced by the LMC. We find that the members of the possible pairs SL349-SL353 and BRHT4b-NGC1839 have a similar age ($t = 1.00 \pm 0.12$ Gyr and $t = 140 \pm 15$ Myr, respectively), thus possibly hinting to a common origin of their member systems We also find that all candidate pairs in our sample show evidence of intra-cluster overdensities that can be a possible indication of real binarity. Particularly interesting is the case of SL349-SL353. In fact, SL353 is relatively close to the condition of critical filling, thus suggesting that these systems might actually constitute an energetically bound pair. It is therefore key to pursue a detailed kinematic screening of such clusters, without which, at present, we do not dare making a conclusive statement about the true nature of this putative pair. △ Less

Submitted 7 November, 2017; originally announced November 2017.

Comments: 14 pages, 12 figures, 3 tables; Accepted for publication in MNRAS

Mass models of NGC 6624 without an intermediate-mass black hole

Authors: Mark Gieles, Eduardo Balbinot, Rashid Yaaqib, Vincent Henault-Brunet, Alice Zocchi, Miklos Peuten, Peter G. Jonker

Abstract: An intermediate-mass black hole (IMBH) was recently reported to reside in the centre of the Galactic globular cluster (GC) NGC 6624, based on timing observations of a millisecond pulsar (MSP) located near the cluster centre in projection. We present dynamical models with multiple mass components of NGC 6624 - without an IMBH - which successfully describe the surface brightness profile and proper m… ▽ More An intermediate-mass black hole (IMBH) was recently reported to reside in the centre of the Galactic globular cluster (GC) NGC 6624, based on timing observations of a millisecond pulsar (MSP) located near the cluster centre in projection. We present dynamical models with multiple mass components of NGC 6624 - without an IMBH - which successfully describe the surface brightness profile and proper motion kinematics from the Hubble Space Telescope (HST) and the stellar mass function at different distances from the cluster centre. The maximum line-of-sight acceleration at the position of the MSP accommodates the inferred acceleration of the MSP, as derived from its first period derivative. With discrete realizations of the models we show that the higher-order period derivatives - which were previously used to derive the IMBH mass - are due to passing stars and stellar remnants, as previously shown analytically in literature. We conclude that there is no need for an IMBH to explain the timing observations of this MSP. △ Less

Submitted 29 November, 2017; v1 submitted 20 September, 2017; originally announced September 2017.

Comments: 8 pages, 7 figures, MNRAS. Updated to match final journal style

Journal ref: Mon. Not. R. Astron. Soc. 473, 4832-4839 (2018)

Testing lowered isothermal models with direct N-body simulations of globular clusters - II: Multimass models

Authors: Miklos Peuten, Alice Zocchi, Mark Gieles, Vincent Hénault-Brunet

Abstract: Lowered isothermal models, such as the multimass Michie-King models, have been successful in describing observational data of globular clusters. In this study we assess whether such models are able to describe the phase space properties of evolutionary $N$-body models. We compare the multimass models as implemented in (Gieles \& Zocchi) to $N$-body models of star clusters with different retention… ▽ More Lowered isothermal models, such as the multimass Michie-King models, have been successful in describing observational data of globular clusters. In this study we assess whether such models are able to describe the phase space properties of evolutionary $N$-body models. We compare the multimass models as implemented in (Gieles \& Zocchi) to $N$-body models of star clusters with different retention fractions for the black holes and neutron stars evolving in a tidal field. We find that multimass models successfully reproduce the density and velocity dispersion profiles of the different mass components in all evolutionary phases and for different remnants retention. We further use these results to study the evolution of global model parameters. We find that over the lifetime of clusters, radial anisotropy gradually evolves from the low-mass to the high-mass components and we identify features in the properties of observable stars that are indicative of the presence of stellar-mass black holes. We find that the model velocity scale depends on mass as $m^{-δ}$, with $δ\simeq0.5$ for almost all models, but the dependence of central velocity dispersion on $m$ can be shallower, depending on the dark remnant content, and agrees well with that of the $N$-body models. The reported model parameters, and correlations amongst them, can be used as theoretical priors when fitting these types of mass models to observational data. △ Less

Submitted 18 July, 2017; v1 submitted 6 February, 2017; originally announced February 2017.

Comments: 28 pages, 22 figures, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, 2017, Volume 470, Issue 3, p.2736-2761

Radial anisotropy in omega Cen limiting the room for an intermediate-mass black hole

Authors: Alice Zocchi, Mark Gieles, Vincent Hénault-Brunet

Abstract: Finding an intermediate-mass black hole (IMBH) in a globular cluster (or proving its absence) would provide valuable insights into our understanding of galaxy formation and evolution. However, it is challenging to identify a unique signature of an IMBH that cannot be accounted for by other processes. Observational claims of IMBH detection are indeed often based on analyses of the kinematics of sta… ▽ More Finding an intermediate-mass black hole (IMBH) in a globular cluster (or proving its absence) would provide valuable insights into our understanding of galaxy formation and evolution. However, it is challenging to identify a unique signature of an IMBH that cannot be accounted for by other processes. Observational claims of IMBH detection are indeed often based on analyses of the kinematics of stars in the cluster core, the most common signature being a rise in the velocity dispersion profile towards the centre of the system. Unfortunately, this IMBH signal is degenerate with the presence of radially-biased pressure anisotropy in the globular cluster. To explore the role of anisotropy in shaping the observational kinematics of clusters, we analyse the case of omega Cen by comparing the observed profiles to those calculated from the family of LIMEPY models, that account for the presence of anisotropy in the system in a physically motivated way. The best-fit radially anisotropic models reproduce the observational profiles well, and describe the central kinematics as derived from Hubble Space Telescope proper motions without the need for an IMBH. △ Less

Submitted 2 February, 2017; originally announced February 2017.

Comments: 14 pages, 7 figures. Accepted for publication in MNRAS

The properties of energetically unbound stars in stellar clusters

Authors: Ian Claydon, Mark Gieles, Alice Zocchi

Abstract: Several Milky Way star clusters show a roughly flat velocity dispersion profile at large radii, which is not expected from models with a tidal cut-off energy. Possible explanations for this excess velocity include: the effects of a dark matter halo, modified gravity theories and energetically unbound stars inside of clusters. These stars are known as potential escapers (PEs) and can exist indefini… ▽ More Several Milky Way star clusters show a roughly flat velocity dispersion profile at large radii, which is not expected from models with a tidal cut-off energy. Possible explanations for this excess velocity include: the effects of a dark matter halo, modified gravity theories and energetically unbound stars inside of clusters. These stars are known as potential escapers (PEs) and can exist indefinitely within clusters which are on circular orbits. Through a series of N-body simulations of star cluster systems, where we vary the galactic potential, orbital eccentricity and stellar mass function, we investigate the properties of the PEs and their effects on the kinematics. We derive a prediction for the scaling of the velocity dispersion at the Jacobi surface due to PEs, as a function of cluster mass, angular velocity of the cluster orbit, and slope of the mass profile of the host galaxy. We see a tentative signal of the mass and orbital velocity dependence in kinematic data of globular clusters from literature. We also find that the fraction of PEs depends sensitively on the galactic mass profile, reaching as high as 40% in the cusp of a Navarro-Frenk-White profile and as the velocity anisotropy also depends on the slope of the galactic mass profile, we conclude that PEs provide an independent way of inferring the properties of the dark matter mass profile at the galactic radius of (globular) clusters in the Gaia era. △ Less

Submitted 28 April, 2017; v1 submitted 7 December, 2016; originally announced December 2016.

Comments: 16 pages, 17 figures, MNRAS, updated to match final journal style

Journal ref: Mon. Not. R. Astron. Soc. 466, 3937-3950 (2017)

A stellar-mass black hole population in the globular cluster NGC 6101?

Authors: Miklos Peuten, Alice Zocchi, Mark Gieles, Alessia Gualandris, Vincent Henault-Brunet

Abstract: Dalessandro et al. observed a similar distribution for blue straggler stars and main-sequence turn-off stars in the Galactic globular cluster NGC 6101, and interpreted this feature as an indication that this cluster is not mass-segregated. Using direct N-body simulations, we find that a significant amount of mass segregation is expected for a cluster with the mass, radius and age of NGC 6101. Ther… ▽ More Dalessandro et al. observed a similar distribution for blue straggler stars and main-sequence turn-off stars in the Galactic globular cluster NGC 6101, and interpreted this feature as an indication that this cluster is not mass-segregated. Using direct N-body simulations, we find that a significant amount of mass segregation is expected for a cluster with the mass, radius and age of NGC 6101. Therefore, the absence of mass segregation cannot be explained by the argument that the cluster is not yet dynamically evolved. By varying the retention fraction of stellar-mass black holes, we show that segregation is not observable in clusters with a high black hole retention fraction (>50% after supernova kicks and >50% after dynamical evolution). Yet all model clusters have the same amount of mass segregation in terms of the decline of the mean mass of stars and remnants with distance to the centre. We also discuss how kinematics can be used to further constrain the presence of a stellar-mass black hole population and distinguish it from the effect of an intermediate-mass black hole. Our results imply that the kick velocities of black holes are lower than those of neutron stars. The large retention fraction during its dynamical evolution can be explained if NGC 6101 formed with a large initial radius in a Milky Way satellite. △ Less

Submitted 6 September, 2016; originally announced September 2016.

Comments: 10 pages, 7 figures

Journal ref: MNRAS (November 01, 2016) 462 (3): 2333

Testing lowered isothermal models with direct N-body simulations of globular clusters

Authors: Alice Zocchi, Mark Gieles, Vincent Hénault-Brunet, Anna Lisa Varri

Abstract: Several self-consistent models have been proposed, aiming at describing the phase space distribution of stars in globular clusters. This study explores the ability of the recently proposed LIMEPY models (Gieles & Zocchi) to reproduce the dynamical properties of direct N-body models of a cluster in a tidal field, during its entire evolution. These dynamical models include prescriptions for the trun… ▽ More Several self-consistent models have been proposed, aiming at describing the phase space distribution of stars in globular clusters. This study explores the ability of the recently proposed LIMEPY models (Gieles & Zocchi) to reproduce the dynamical properties of direct N-body models of a cluster in a tidal field, during its entire evolution. These dynamical models include prescriptions for the truncation and the degree of radially-biased anisotropy contained in the system, allowing us to explore the interplay between the role of anisotropy and tides in various stages of the life of star clusters. We show that the amount of anisotropy in an initially tidally underfilling cluster increases in the pre-collapse phase, and then decreases with time, due to the effect of the external tidal field on its spatial truncation. This is reflected in the correspondent model parameters, and the best-fit models reproduce the main properties of the cluster at all stages of its evolution, except for the phases immediately preceding and following core collapse. We also notice that the best-fit LIMEPY models are significantly different from isotropic King models, especially in the first part of the evolution of the cluster. Our results put limits on the amount of radial anisotropy that can be expected for clusters evolving in a tidal field, which is important to understand other factors that could give rise to similar observational signatures, such as the presence of an intermediate-mass black hole. △ Less

Submitted 6 May, 2016; originally announced May 2016.

Comments: 20 pages, 11 figures. Accepted for publication in MNRAS

A class of spherical, truncated, anisotropic models for application to globular clusters

Authors: Ruggero de Vita, Giuseppe Bertin, Alice Zocchi

Abstract: Recently, a class of non-truncated radially-anisotropic models (the so-called $f^{(ν)}$-models), originally constructed in the context of violent relaxation and modeling of elliptical galaxies, has been found to possess interesting qualities in relation to observed and simulated globular clusters. In view of new applications to globular clusters, we improve this class of models along two direction… ▽ More Recently, a class of non-truncated radially-anisotropic models (the so-called $f^{(ν)}$-models), originally constructed in the context of violent relaxation and modeling of elliptical galaxies, has been found to possess interesting qualities in relation to observed and simulated globular clusters. In view of new applications to globular clusters, we improve this class of models along two directions. To make them more suitable for the description of small stellar systems hosted by galaxies, we introduce a 'tidal' truncation (by means of a procedure that guarantees full continuity of the distribution function). The new $f_T^{(ν)}$-models are shown to provide a better fit to the observed photometric and spectroscopic profiles for a sample of 13 globular clusters studied earlier by means of non-truncated models; interestingly, the best-fit models also perform better with respect to the radial-orbit instability. Then we design a flexible but simple two-component family of truncated models, to study the separate issues of mass segregation and of multiple populations. We do not aim at a fully realistic description of globular clusters, to compete with the description currently obtained by means of dedicated simulations. The goal here is to try to identify the simplest models, that is, those with the smallest number of free parameters, but still able to provide a reasonable description for clusters that are evidently beyond the reach of one-component models: with this tool we aim at identifying the key factors that characterize mass segregation or the presence of multiple populations. To reduce the relevant parameter space, we formulate a few physical arguments (based on recent observations and simulations). A first application to two well-studied globular clusters is briefly described and discussed. △ Less

Submitted 18 March, 2016; originally announced March 2016.

Comments: 15 pages, 13 figures, Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 590, A16 (2016)

A family of lowered isothermal models

Authors: Mark Gieles, Alice Zocchi

Abstract: We present a family of self-consistent, spherical, lowered isothermal models, consisting of one or more mass components, with parameterised prescriptions for the energy truncation and for the amount of radially biased pressure anisotropy. The models are particularly suited to describe the phase-space density of stars in tidally limited, mass-segregated star clusters in all stages of their life-cyc… ▽ More We present a family of self-consistent, spherical, lowered isothermal models, consisting of one or more mass components, with parameterised prescriptions for the energy truncation and for the amount of radially biased pressure anisotropy. The models are particularly suited to describe the phase-space density of stars in tidally limited, mass-segregated star clusters in all stages of their life-cycle. The models extend a family of isotropic, single-mass models by Gomez-Leyton and Velazquez, of which the well-known Woolley, King and Wilson (in the non-rotating and isotropic limit) models are members. We derive analytic expressions for the density and velocity dispersion components in terms of potential and radius, and introduce a fast model solver in PYTHON (LIMEPY), that can be used for data fitting or for generating discrete samples. △ Less

Submitted 17 September, 2016; v1 submitted 9 August, 2015; originally announced August 2015.

Comments: 17 pages, 10 figures, 4 appendices, MNRAS, updated to match final journal style

Journal ref: Mon. Not. R. Astron. Soc. 454, 576-592 (2015)

Biases in the determination of dynamical parameters of star clusters: today and in the Gaia era

Authors: A. Sollima, H. Baumgardt, A. Zocchi, E. Balbinot, M. Gieles, V. Henault-Brunet, A. L. Varri

Abstract: The structural and dynamical properties of star clusters are generally derived by means of the comparison between steady-state analytic models and the available observables. With the aim of studying the biases of this approach, we fitted different analytic models to simulated observations obtained from a suite of direct N-body simulations of star clusters in different stages of their evolution and… ▽ More The structural and dynamical properties of star clusters are generally derived by means of the comparison between steady-state analytic models and the available observables. With the aim of studying the biases of this approach, we fitted different analytic models to simulated observations obtained from a suite of direct N-body simulations of star clusters in different stages of their evolution and under different levels of tidal stress to derive mass, mass function and degree of anisotropy. We find that masses can be under/over-estimated up to 50% depending on the degree of relaxation reached by the cluster, the available range of observed masses and distances of radial velocity measures from the cluster center and the strength of the tidal field. The mass function slope appears to be better constrainable and less sensitive to model inadequacies unless strongly dynamically evolved clusters and a non-optimal location of the measured luminosity function are considered. The degree and the characteristics of the anisotropy developed in the N-body simulations are not adequately reproduced by popular analytic models and can be detected only if accurate proper motions are available. We show how to reduce the uncertainties in the mass, mass-function and anisotropy estimation and provide predictions for the improvements expected when Gaia proper motions will be available in the near future. △ Less

Submitted 11 May, 2015; originally announced May 2015.

Comments: 14 pages, 8 figures, accepted for publication by MNRAS

On the uniqueness of kinematical signatures of intermediate-mass black holes in globular clusters

Authors: Alice Zocchi, Mark Gieles, Vincent Hénault-Brunet

Abstract: Finding an intermediate-mass black hole (IMBH) in a globular cluster (GC), or proving its absence, is a crucial ingredient in our understanding of galaxy formation and evolution. The challenge is to identify a unique signature of an IMBH that cannot be accounted for by other processes. Observational claims of IMBH detection are often based on analyses of the kinematics of stars, such as a rise in… ▽ More Finding an intermediate-mass black hole (IMBH) in a globular cluster (GC), or proving its absence, is a crucial ingredient in our understanding of galaxy formation and evolution. The challenge is to identify a unique signature of an IMBH that cannot be accounted for by other processes. Observational claims of IMBH detection are often based on analyses of the kinematics of stars, such as a rise in the velocity dispersion profile towards the centre. In this contribution we discuss the degeneracy between this IMBH signal and pressure anisotropy in the GC. We show that that by considering anisotropic models it is possible to partially explain the innermost shape of the projected velocity dispersion profile, even though models that do not account for an IMBH do not exhibit a cusp in the centre. △ Less

Submitted 21 January, 2015; originally announced January 2015.

Comments: 4 pages, 2 figures. To be published in the Proceedings IAU Symposium No. 312, Star Clusters and Black Holes in Galaxies Across Cosmic Time

Asking gender questions: Results from a survey of gender and question asking among UK Astronomers at NAM2014

Authors: Jonathan Pritchard, Karen Masters, James Allen, Filippo Contenta, Leo Huckvale, Stephen Wilkins, Alice Zocchi

Abstract: We report on a survey of astronomers asking questions at the most recent National Astronomy Meeting (NAM2014). The gender balance of both speakers and session chairs at NAM (31% and 29% women respectively) closely matched that of attendees (28% female). However, we find that women were under-represented among question askers (just 18% female). Women were especially under-represented in asking the… ▽ More We report on a survey of astronomers asking questions at the most recent National Astronomy Meeting (NAM2014). The gender balance of both speakers and session chairs at NAM (31% and 29% women respectively) closely matched that of attendees (28% female). However, we find that women were under-represented among question askers (just 18% female). Women were especially under-represented in asking the first question (only 14% of first questions asked by women), but when the Q&A session reached four or more questions, women and men were observed to ask roughly equal numbers of questions. We found a small, but statistically insignificant, increase in the fraction of questions from women in sessions where the chair was also female. On average $2.2\pm0.1$ questions were asked per talk, with no detectable difference in the number of questions asked of female and male speakers, but on average female chairs solicited slightly fewer questions than male chairs. We compare these results to a similar study by Davenport et al. (2014) for the AAS, who also found under-representation of women amongst question askers, but saw more pronounced gender effects when a session chair or speaker were female. We conclude with suggestions for improving the balance of questions at future astronomy meetings. △ Less

Submitted 15 December, 2014; originally announced December 2014.

Comments: 15 pages, 7 figures, comments welcome

Journal ref: Astronomy and Geophysics, 55(6),8-12 (2014)

Rotating Globular Clusters

Authors: P. Bianchini, A. L Varri, G. Bertin, A. Zocchi

Abstract: Internal rotation is considered to play a major role in the dynamics of some globular clusters. However, in only few cases it has been studied by quantitative application of realistic and physically justified global models. Here we present a dynamical analysis of the photometry and three-dimensional kinematics of omega Cen, 47 Tuc, and M15, by means of a recently introduced family of self-consiste… ▽ More Internal rotation is considered to play a major role in the dynamics of some globular clusters. However, in only few cases it has been studied by quantitative application of realistic and physically justified global models. Here we present a dynamical analysis of the photometry and three-dimensional kinematics of omega Cen, 47 Tuc, and M15, by means of a recently introduced family of self-consistent axisymmetric rotating models. The three clusters, characterized by different relaxation conditions, show evidence of differential rotation and deviations from sphericity. The combination of line-of-sight velocities and proper motions allows us to determine their internal dynamics, predict their morphology, and estimate their dynamical distance. The well-relaxed cluster 47 Tuc is very well interpreted by our model; internal rotation is found to explain the observed morphology. For M15, we provide a global model in good agreement with the data, including the central behavior of the rotation profile and the shape of the ellipticity profile. For the partially relaxed cluster omega Cen, the selected model reproduces the complex three-dimensional kinematics; in particular the observed anisotropy profile, characterized by a transition from isotropy, to weakly-radial anisotropy, and then to tangential anisotropy in the outer parts. The discrepancy found for the steep central gradient in the observed line-of-sight velocity dispersion profile and for the ellipticity profile is ascribed to the condition of only partial relaxation of this cluster and the interplay between rotation and radial anisotropy. △ Less

Submitted 26 May, 2013; originally announced May 2013.

Comments: 19 pages, 14 figures, accepted for publication in the Astrophysical Journal

On the density profile of the globular cluster M92

Authors: A. Di Cecco, A. Zocchi, A. L. Varri, M. Monelli, G. Bertin, G. Bono, P. B. Stetson, M. Nonino, R. Buonanno, I. Ferraro, G. Iannicola, A. Kunder, A. R. Walker

Abstract: We present new number density and surface brightness profiles for the globular cluster M92 (NGC 6341). These profiles are calculated from optical images collected with the CCD mosaic camera MegaCam at the Canada-France-Hawaii-Telescope and with the Advanced Camera for Surveys on the Hubble Space Telescope. The ground-based data were supplemented with the Sloan Digital Sky Survey photometric catalo… ▽ More We present new number density and surface brightness profiles for the globular cluster M92 (NGC 6341). These profiles are calculated from optical images collected with the CCD mosaic camera MegaCam at the Canada-France-Hawaii-Telescope and with the Advanced Camera for Surveys on the Hubble Space Telescope. The ground-based data were supplemented with the Sloan Digital Sky Survey photometric catalog. Special care was taken to discriminate candidate cluster stars from field stars and to subtract the background contamination from both profiles. By examining the contour levels of the number density, we found that the stellar distribution becomes clumpy at radial distances larger than about 13 arcminutes, and there is no preferred orientation of contours in space. We performed detailed fits of King and Wilson models to the observed profiles. The best-fit models underestimate the number density inside the core radius. Wilson models better represent the observations, in particular in the outermost cluster regions: the good global agreement of these models with the observations suggests that there is no need to introduce an extra-tidal halo to explain the radial distribution of stars at large radial distances. The best-fit models for the number density and the surface brightness profiles are different, even though they are based on the same observations. Additional tests support the evidence that this fact reflects the difference in the radial distribution of the stellar tracers that determine the observed profiles (main sequence stars for the number density, bright evolved stars for the surface brightness). △ Less

Submitted 7 February, 2013; originally announced February 2013.

Comments: 18 pages, 10 figures, Accepted by AJ

The flattening of globular clusters: internal rotation or velocity anisotropy?

Authors: Paolo Bianchini, Anna Lisa Varri, Giuseppe Bertin, Alice Zocchi

Abstract: Internal rotation is considered to play a major role in determining the structure and dynamics of some globular clusters. We present a dynamical analysis of the photometry and three-dimensional kinematics of 47 Tuc and omega Cen, by means of a new family of self-consistent axisymmetric rotating models. The combined use of line-of-sight velocities and proper motions allows us to obtain a global des… ▽ More Internal rotation is considered to play a major role in determining the structure and dynamics of some globular clusters. We present a dynamical analysis of the photometry and three-dimensional kinematics of 47 Tuc and omega Cen, by means of a new family of self-consistent axisymmetric rotating models. The combined use of line-of-sight velocities and proper motions allows us to obtain a global description of the internal dynamical structure of the objects together with an estimate of their dynamical distances. The well-relaxed cluster 47 Tuc is very well interpreted by our dynamical models; in particular, internal rotation is found to explain the observed morphology. For the partially relaxed cluster omega Cen, the selected model provides a good representation of its complex three-dimensional kinematics, in general qualitative agreement with the observed anisotropy profile, which is characterized by tangential anisotropy in the outer parts; discrepancies are found between the observed and the expected ellipticity profile and are ascribed to the presence of a high degree of radial anisotropy in the intermediate region and to its interplay with rotation. △ Less

Submitted 4 February, 2013; originally announced February 2013.

Comments: 4 pages, 3 figures. To appear in the conference proceedings "Reading the book of globular clusters with the lens of stellar evolution", Mem. S. A. It. Eds. P. Ventura, C. Charbonnel, M. Castellani and M. Di Criscienzo

A dynamical study of Galactic globular clusters under different relaxation conditions

Authors: A. Zocchi, G. Bertin, A. L. Varri

Abstract: We perform a systematic combined photometric and kinematic analysis of a sample of globular clusters under different relaxation conditions, based on their core relaxation time (as listed in available catalogs), by means of two well-known families of spherical stellar dynamical models. Systems characterized by shorter relaxation time scales are expected to be better described by isotropic King mode… ▽ More We perform a systematic combined photometric and kinematic analysis of a sample of globular clusters under different relaxation conditions, based on their core relaxation time (as listed in available catalogs), by means of two well-known families of spherical stellar dynamical models. Systems characterized by shorter relaxation time scales are expected to be better described by isotropic King models, while less relaxed systems might be interpreted by means of non-truncated, radially-biased anisotropic f^(ν) models, originally designed to represent stellar systems produced by a violent relaxation formation process and applied here for the first time to the study of globular clusters. The comparison between dynamical models and observations is performed by fitting simultaneously surface brightness and velocity dispersion profiles. For each globular cluster, the best-fit model in each family is identified, along with a full error analysis on the relevant parameters. Detailed structural properties and mass-to-light ratios are also explicitly derived. We find that King models usually offer a good representation of the observed photometric profiles, but often lead to less satisfactory fits to the kinematic profiles, independently of the relaxation condition of the systems. For some less relaxed clusters, f^(ν) models provide a good description of both observed profiles. Some derived structural characteristics, such as the total mass or the half-mass radius, turn out to be significantly model-dependent. The analysis confirms that, to answer some important dynamical questions that bear on the formation and evolution of globular clusters, it would be highly desirable to acquire larger numbers of accurate kinematic data-points, well distributed over the cluster field. △ Less

Submitted 6 January, 2012; originally announced January 2012.

Comments: 18 pages, 7 figures. Accepted for publication in Astronomy & Astrophysics