Showing 1–37 of 37 results for author: Hypki, A

Multiband gravitational wave observations of eccentric escaping binary black holes from globular clusters

Authors: Yuetong Zhao, Abbas Askar, Youjun Lu, Zhoujian Cao, Mirek Giersz, Grzegorz Wiktorowicz, Arkadiusz Hypki, Lucas Hellstrom, Sohaib Ali, Wei-Tou Ni

Abstract: Stellar-mass binary black holes (sBBHs) formed in globular clusters (GCs) are promising sources for multiband gravitational wave (GW) observations, particularly with low- and middle-frequency detectors. These sBBHs can retain detectable eccentricities when they enter the sensitivity bands of low-frequency GW observatories. We study multiband GW observations of eccentric sBBHs that escape from GC m… ▽ More Stellar-mass binary black holes (sBBHs) formed in globular clusters (GCs) are promising sources for multiband gravitational wave (GW) observations, particularly with low- and middle-frequency detectors. These sBBHs can retain detectable eccentricities when they enter the sensitivity bands of low-frequency GW observatories. We study multiband GW observations of eccentric sBBHs that escape from GC models simulated with the MOCCA code, focusing on how low- and middle-frequency detectors can constrain their eccentricities and other parameters. Using Monte Carlo simulations, we generate ten realizations of cosmic sBBHs by combining the MOCCA sample with a cosmological model for GC formation and evolution. We then assess their detectability and the precision of parameter estimation. Our results show that LISA, Taiji, the LISA-Taiji network (LT), and AMIGO could detect $0.8\pm0.7$, $11.6\pm2.0$, $15.4\pm2.7$, and $7.9\pm1.3$ escaping sBBHs, respectively, over four years, while LT-AMIGO could detect $20.6\pm3.0$ multiband sBBHs in the same period. LT and AMIGO can measure initial eccentricities with relative errors of approximately $10^{-6}-2\times10^{-4}$ and $10^{-3}-0.7$, respectively. Joint LT-AMIGO observations have a similar ability to estimate eccentricities as LT alone. △ Less

Submitted 15 December, 2025; originally announced December 2025.

Comments: 19 pages, 3 tables, and 10 figures, submitted to Astrophysical Journal (ApJ)

Simulations of Globular Cluster Evolution with Multiple Stellar Populations

Authors: Mirek Giersz, Abbas Askar, Arkadiusz Hypki, Jongsuk Hong, Grzegorz Wiktorowicz, Lucas Hellstrom

Abstract: The formation of stars with light-element abundance variations in globular clusters and the subsequent dynamical evolution of these multiple populations remains an open question. One of the most widely discussed is the AGB scenario, in which chemically processed material from the envelopes of AGB stars mixes with re-accreted primordial gas flowing into the center of the cluster. Based on this scen… ▽ More The formation of stars with light-element abundance variations in globular clusters and the subsequent dynamical evolution of these multiple populations remains an open question. One of the most widely discussed is the AGB scenario, in which chemically processed material from the envelopes of AGB stars mixes with re-accreted primordial gas flowing into the center of the cluster. Based on this scenario, more than two hundred MOCCA simulations of cluster evolution have been carried out, incorporating additional physical processes related to the external environment of globular clusters and the initial properties of multiple stellar populations. Analysis of the simulations shows that most observed properties of multiple stellar populations and the global parameters of Milky Way clusters are well reproduced, with the exception of the correlation between cluster mass and the fraction of second-population stars. We present a speculative scenario of globular cluster evolution that may account for the observed properties of Milky Way clusters, including the correlation between cluster mass and the fraction of enriched stars. The scenario further predicts that, under certain conditions, the pristine first population can be more centrally concentrated than the enriched second population, as observed in some clusters. \end{abstract △ Less

Submitted 8 October, 2025; originally announced October 2025.

Comments: Comments: 10 pages, 6 figures, 1 table. To appear in the Proceedings of IAU Symposium 398 and MODEST-25: "Compact Objects and Binaries in Dense Stellar Systems" (Seoul, South Korea, June 2025)

Formation Channels of Gravitationally Resolvable Double White Dwarf Binaries Inside Globular Clusters

Authors: Lucas Hellström, Mirosław Giersz, Abbas Askar, Arkadiusz Hypki, Yuetong Zhao, Youjun Lu, Siqi Zhang, Verónica Vázquez-Aceves, Grzegorz Wiktorowicz

Abstract: Current gravitational wave detectors are sensitive to coalescing black holes and neutron stars. However, double white dwarfs (DWDs) have long been recognized as promising sources of gravitational waves, and upcoming detectors like LISA will be able to observe these systems in abundance. DWDs are expected to be the dominant gravitational wave (GW) sources in parts of the LISA frequency range, makin… ▽ More Current gravitational wave detectors are sensitive to coalescing black holes and neutron stars. However, double white dwarfs (DWDs) have long been recognized as promising sources of gravitational waves, and upcoming detectors like LISA will be able to observe these systems in abundance. DWDs are expected to be the dominant gravitational wave (GW) sources in parts of the LISA frequency range, making it crucial to understand their formation for future detections. The Milky Way contains many white dwarfs (WDs) in both the field and star clusters, promising a rich population of DWDs for LISA. However, the large number of sources may make it difficult to resolve individual binaries, and DWDs in the field and clusters often have similar properties, complicating the identification of their origins from GW signals alone. In this work, we focus on eccentric and tight DWDs, which cannot form in the field, but require dynamical interactions in dense clusters to increase their eccentricity after circularization through mass transfer phases and common-envelope evolution during binary evolution. These binaries may also form in three- and four-body dynamical interactions where a DWD binary may directly form with high eccentricity and low separation. Our results show that we should expect eccentric and tight DWDs in clusters that can provide meaningful GW signal, however, the number is low; with an upper limit of 10-15 in the MW. These can be used to independently obtain distances of their host cluster. △ Less

Submitted 16 June, 2025; originally announced June 2025.

Comments: 10 pages, 5 figures, 2 tables

Report number: 702

Journal ref: A&A 702, A131 (2025)

Rapid formation of a very massive star >50000 $M_\odot$ and subsequently an IMBH from runaway collisions. Direct N-body and Monte Carlo simulations of dense star clusters

Authors: Marcelo C. Vergara, Abbas Askar, Albrecht W. H. Kamlah, Rainer Spurzem, Francesco Flammini Dotti, Dominik R. G. Schleicher, Manuel Arca Sedda, Arkadiusz Hypki, Mirek Giersz, Jarrod Hurley, Peter Berczik, Andres Escala, Nils Hoyer, Nadine Neumayer, Xiaoying Pang, Ataru Tanikawa, Renyue Cen, Thorsten Naab

Abstract: Context. We present simulations of a massive young star cluster using \textsc{Nbody6++GPU} and \textsc{MOCCA}. The cluster is initially more compact than previously published models, with one million stars, a total mass of $5.86 \times 10^5~\mathrm{M}_{\odot}$, and a half-mass radius of $0.1~\mathrm{pc}$. Aims. We analyse the formation and growth of a very massive star (VMS) through successive s… ▽ More Context. We present simulations of a massive young star cluster using \textsc{Nbody6++GPU} and \textsc{MOCCA}. The cluster is initially more compact than previously published models, with one million stars, a total mass of $5.86 \times 10^5~\mathrm{M}_{\odot}$, and a half-mass radius of $0.1~\mathrm{pc}$. Aims. We analyse the formation and growth of a very massive star (VMS) through successive stellar collisions and investigate the subsequent formation of an intermediate-mass black hole (IMBH) in the core of a dense star cluster. Methods. We use both direct \textit{N}-body and Monte Carlo simulations, incorporating updated stellar evolution prescriptions (SSE/BSE) tailored to massive stars and VMSs. These include revised treatments of stellar radii, rejuvenation, and mass loss during collisions. While the prescriptions represent reasonable extrapolations into the VMS regime, the internal structure and thermal state of VMSs formed through stellar collisions remain uncertain, and future work may require further refinement. Results. We find that runaway stellar collisions in the cluster core produce a VMS exceeding $5 \times 10^4~\mathrm{M}_{\odot}$ within 5 Myr, which subsequently collapses into an IMBH. Conclusions. Our model suggests that dense stellar environments may enable the formation of very massive stars and massive black hole seeds through runaway stellar collisions. These results provide a potential pathway for early black hole growth in star clusters and offer theoretical context for interpreting recent JWST observations of young, compact clusters at high redshift. △ Less

Submitted 17 December, 2025; v1 submitted 12 May, 2025; originally announced May 2025.

Comments: 21 pages, 11 figures, Accepted for publication in A&A

Multiple stellar populations in MOCCA globular cluster models: Transient spatial over-concentration of pristine red giant stars driven by strong dynamical encounters

Authors: M. Giersz, A. Askar, A Hypki, J. Hong, G. Wiktorowicz, L. Hellström

Abstract: Recent findings show that, in some Milky Way globular clusters (GCs), pristine red giant branch (RGB) stars are more centrally concentrated than enriched RGB stars. This contradicts most multiple-population formation scenarios, which predict that the enriched population 2p should initially be more concentrated than the pristine population 1P. We analyze a MOCCA GC model that exhibits a higher spat… ▽ More Recent findings show that, in some Milky Way globular clusters (GCs), pristine red giant branch (RGB) stars are more centrally concentrated than enriched RGB stars. This contradicts most multiple-population formation scenarios, which predict that the enriched population 2p should initially be more concentrated than the pristine population 1P. We analyze a MOCCA GC model that exhibits a higher spatial concentration of 1P RGB stars than 2P RGB stars at 13 Gyr. The MOCCA models assume the asymptotic giant branch (AGB) pollution scenario, where 2P stars are initially more concentrated than 1P stars. Our results indicate that the observed spatial distributions of multiple populations, and possibly their kinematics, are significantly shaped by dynamical interactions. These interactions preferentially eject 2P RGB progenitors from the central regions, leading to a transient over-concentration of 1P RGB stars at late times. This effect is particularly relevant for GCs with present-day of a few $10^5 M_{\odot}$, which have retained only about 10 - 20 percent of their initial mass. Such clusters may appear dynamically young due to heating from a black hole subsystem, even if they have undergone significant mass loss and dynamical evolution. Additionally, the relatively small number of RGB stars in these clusters suggests that interpreting the spatial distributions of multiple populations solely from RGB stars may lead to biased conclusions about the overall distribution of 2P and 1P. The apparent over-concentration of the 1P relative to the 2P is likely a transient effect driven by the preferential removal of 2P RGB progenitors via strong dynamical encounters. MOCCA models of multiple stellar populations based on the AGB scenario may explain anomalous features observed in some Galactic GCs, such as NGC 3201 and NGC 6101. △ Less

Submitted 26 May, 2025; v1 submitted 24 February, 2025; originally announced February 2025.

Comments: 7 pages, 6 figures, accepted to A&A letters

Journal ref: A&A 698, L11 (2025)

Ultraluminous X-ray sources in Globular Clusters

Authors: Grzegorz Wiktorowicz, Mirek Giersz, Abbas Askar, Arkadiusz Hypki, Lucas Helstrom

Abstract: This paper investigates the formation, populations, and evolutionary paths of UltraLuminous X-ray Sources (ULXs) within Globular Clusters (GCs). ULXs, characterised by their extreme X-ray luminosities, present a challenge to our understanding of accretion physics and compact object formation. While previous studies have largely focused on field populations, this research examines the unique enviro… ▽ More This paper investigates the formation, populations, and evolutionary paths of UltraLuminous X-ray Sources (ULXs) within Globular Clusters (GCs). ULXs, characterised by their extreme X-ray luminosities, present a challenge to our understanding of accretion physics and compact object formation. While previous studies have largely focused on field populations, this research examines the unique environment of GCs, where dynamical interactions play a significant role. Using the MOCCA Monte Carlo code, we explore how dynamics influences ULX populations within these dense stellar clusters. Our findings reveal that dynamical processes, such as binary hardening and exchanges, can both facilitate and impede ULX formation in GCs. The study explores the impact of parameters including the initial binary fraction, tidal filling, and multiple stellar populations on the evolution of ULXs. We find that non-tidally filling clusters exhibit significantly larger ULX populations compared to tidally filling ones. The results indicate that the apparent scarcity of ULXs in GCs may be related to the older stellar populations of GCs relative to the field. Furthermore, the study identifies a population of "escaper" ULXs, which originate in GCs but are ejected and emit X-rays outside the cluster. These escapers may significantly contribute to the observed field ULX population. △ Less

Submitted 10 April, 2025; v1 submitted 10 January, 2025; originally announced January 2025.

Comments: 12 pages, 3 figures, published in A&A

Journal ref: A&A 696, A90 (2025)

MOCCA: Effects of pristine gas accretion and cluster migration on globular cluster evolution, global parameters, and multiple stellar populations

Authors: Mirek Giersz, Abbas Askar, Arkadiusz Hypki, Jongsuk Hong, Grzegorz Wiktorowicz, Lucas Hellstrom

Abstract: Using the MOCCA code, we study the evolution of globular clusters (GCs) with multiple stellar populations. For this purpose, the MOCCA code has been significantly extended to take into account the formation of an enriched population of stars from re-accreted gas with a time delay after the formation of the pristine population of stars. The possibility of cluster migration in the host galaxy and th… ▽ More Using the MOCCA code, we study the evolution of globular clusters (GCs) with multiple stellar populations. For this purpose, the MOCCA code has been significantly extended to take into account the formation of an enriched population of stars from re-accreted gas with a time delay after the formation of the pristine population of stars. The possibility of cluster migration in the host galaxy and the fact that the pristine population can be described by a model not in virial equilibrium are also taken into account. Gas re-accretion and cluster migration have a decisive impact on the observational parameters of clusters and the ratio of the number of objects between the pristine and enriched populations. The obtained results, together with observational data, suggest a speculative refinement of the AGB scenario that makes it possible to explain some observational data, such as the ratio of the pristine to the enriched populations, the observational fact that for some GCs the pristine population is more concentrated than the enriched one, and possibly a correlation between the ratio of the number of enriched stars to the total number of stars and the mass of the cluster. In this scenario, it is important to take into account the environment in which the cluster lives, the conditions in the galaxy when it formed, and the fact that a significant part of the GCs associated with the Galaxy come from dwarf galaxies that merged with the Milky Way. The initial conditions of GCs in our simulations differ from the widely used typical models, as they require GCs to fill the Roche lobe rather than being highly concentrated within it, imposing strong constraints on their formation locations within the galaxy. △ Less

Submitted 26 May, 2025; v1 submitted 10 November, 2024; originally announced November 2024.

Comments: 17 pages, 13 figures, accepted A&A

Journal ref: A&A 699, A76 (2025)

Energy equipartition in multiple-population globular clusters

Authors: A. R. Livernois, F. I. Aros, E. Vesperini, A. Askar, A. Bellini, M. Giersz, J. Hong, A. Hypki, M. Libralato, T. Ziliotto

Abstract: We present the results of Monte Carlo simulations aimed at exploring the evolution towards energy equipartition of first- (1G) and second-generation (2G) stars in multiple-population globular clusters and how this evolution is affected by the initial differences between the spatial distributions of the two populations. Our results show that these initial differences have fundamental implications f… ▽ More We present the results of Monte Carlo simulations aimed at exploring the evolution towards energy equipartition of first- (1G) and second-generation (2G) stars in multiple-population globular clusters and how this evolution is affected by the initial differences between the spatial distributions of the two populations. Our results show that these initial differences have fundamental implications for the evolution towards energy equipartition of the two populations. We find that 2G stars, which are assumed to be initially more centrally concentrated than 1G stars, are generally characterized by a more rapid evolution towards energy equipartition. The evolution towards energy equipartition depends on the velocity dispersion component and is more rapid for the tangential velocity dispersion. The extent of the present-day differences between the degree of energy equipartition of 2G and 1G stars depends on the cluster's dynamical age and may be more significant in the tangential velocity dispersion and at intermediate distances from the cluster's center around the half-mass radius. △ Less

Submitted 16 October, 2024; originally announced October 2024.

Comments: 13 pages, 13 figures, Published in MNRAS

Gravitational Wave Astronomy With TianQin

Authors: En-Kun Li, Shuai Liu, Alejandro Torres-Orjuela, Xian Chen, Kohei Inayoshi, Long Wang, Yi-Ming Hu, Pau Amaro-Seoane, Abbas Askar, Cosimo Bambi, Pedro R. Capelo, Hong-Yu Chen, Alvin J. K. Chua, Enrique Condés-Breña, Lixin Dai, Debtroy Das, Andrea Derdzinski, Hui-Min Fan, Michiko Fujii, Jie Gao, Mudit Garg, Hongwei Ge, Mirek Giersz, Shun-Jia Huang, Arkadiusz Hypki , et al. (28 additional authors not shown)

Abstract: The opening of the gravitational wave window has significantly enhanced our capacity to explore the universe's most extreme and dynamic sector. In the mHz frequency range, a diverse range of compact objects, from the most massive black holes at the farthest reaches of the Universe to the lightest white dwarfs in our cosmic backyard, generate a complex and dynamic symphony of gravitational wave sig… ▽ More The opening of the gravitational wave window has significantly enhanced our capacity to explore the universe's most extreme and dynamic sector. In the mHz frequency range, a diverse range of compact objects, from the most massive black holes at the farthest reaches of the Universe to the lightest white dwarfs in our cosmic backyard, generate a complex and dynamic symphony of gravitational wave signals. Once recorded by gravitational wave detectors, these unique fingerprints have the potential to decipher the birth and growth of cosmic structures over a wide range of scales, from stellar binaries and stellar clusters to galaxies and large-scale structures. The TianQin space-borne gravitational wave mission is scheduled for launch in the 2030s, with an operational lifespan of five years. It will facilitate pivotal insights into the history of our universe. This document presents a concise overview of the detectable sources of TianQin, outlining their characteristics, the challenges they present, and the expected impact of the TianQin observatory on our understanding of them. △ Less

Submitted 2 December, 2024; v1 submitted 29 September, 2024; originally announced September 2024.

Comments: TianQin Gravitational Wave Whitepaper, 72 pages, 30 figures

MOCCA: Global properties of tidally filling and underfilling globular star clusters with multiple stellar populations

Authors: Arkadiusz Hypki, Enrico Vesperini, Mirek Giersz, Jongsuk Hong, Abbas Askar, Magdalena Otulakowska-Hypka, Lucas Hellstrom, Grzegorz Wiktorowicz

Abstract: We explore the evolution of various properties of multiple-population globular clusters (GCs) for a broad range of initial conditions. We simulated over 200 GC models using the MOCCA Monte Carlo code and find that present-day properties (core and half-light radii, ratio of the number of second-generation (SG) stars to the total number of stars, NSG/NTOT) of these models cover the observed values o… ▽ More We explore the evolution of various properties of multiple-population globular clusters (GCs) for a broad range of initial conditions. We simulated over 200 GC models using the MOCCA Monte Carlo code and find that present-day properties (core and half-light radii, ratio of the number of second-generation (SG) stars to the total number of stars, NSG/NTOT) of these models cover the observed values of these quantities for Milky Way GCs. Starting with a relatively small value of the SG fraction (NSG/NTOT ~ 0.25) and a SG system concentrated in the inner regions of the cluster, we find, in agreement with previous studies, that systems in which the first-generation (FG) is initially tidally filling or slightly tidally underfilling best reproduce the observed ratios of NSG/NTOT and have values of the core and half-light radii typical of those of many Galactic globular clusters. Models in which the FG is initially tidally underfilling retain values of NSG/NTOT close to their initial values. These simulations expand previous investigations and serve to further constrain the viable range of initial parameters and better understand their influence on present-day GC properties. The results of this investigation also provide the basis for our future survey aimed at building specific models to reproduce the observed trends (or lack thereof) between the properties of multiple stellar populations and other clusters properties. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: Accepted to Astronomy and Astrophysics (A&A)

Journal ref: A&A 693, A41 (2025)

Double white dwarf binary population in MOCCA star clusters -- Comparisons with observations of close and wide binaries

Authors: Lucas Hellström, Mirosław Giersz, Arkadiusz Hypki, Diogo Belloni, Abbas Askar, Grzegorz Wiktorowicz

Abstract: There could be a significant population of double white dwarf binaries (DWDs) inside globular clusters (GCs), however, these are often too faint to be individually observed. We have utilized a large number GC models evolved with the Monte Carlo Cluster Simulator (MOCCA) code, to create a large statistical dataset of DWDs. These models include multiple-stellar populations, resulting in two distinct… ▽ More There could be a significant population of double white dwarf binaries (DWDs) inside globular clusters (GCs), however, these are often too faint to be individually observed. We have utilized a large number GC models evolved with the Monte Carlo Cluster Simulator (MOCCA) code, to create a large statistical dataset of DWDs. These models include multiple-stellar populations, resulting in two distinct initial populations: one dense and another less dense. Due to the lower density of one population, a large number of objects escape during the early GC evolution, leading to a high mass-loss rate. In this dataset we have analysed three main groups of DWDs, namely in-cluster binaries, escaped binaries, and isolated evolution of primordial binaries. We compared the properties of these groups to observations of close and wide binaries. We find that the number of escaping DWDs is significantly larger than the number of in-cluster binaries and those that form via the isolated evolution of all promiridial binaries in our GC models. This suggests that dynamics play an important role in the formation of DWDs. For close binaries, we found a good agreement in the separations of escaped binaries and isolated binaries, but in-cluster binaries showed slight differences. We could not reproduce the observed extremely low mass WDs due to the limitations of our stellar and binary evolution prescriptions. For wide binaries, we also found a good agreement in the separations and masses, after accounting for observational selection effects. We conclude that, even though the current observational samples of DWDs are extremely biased and incomplete, our results compare reasonably well with observations. △ Less

Submitted 3 January, 2025; v1 submitted 7 May, 2024; originally announced May 2024.

Comments: 13 pages, 8 figures

Journal ref: A&A 690, A112 (2024)

New Parameters for Star Cluster Dynamics: the role of clusters initial conditions

Authors: Bhavana Bhat, Barbara Lanzoni, Enrico Vesperini, Francesco R. Ferraro, Francisco I. Aros, Abbas Askar, Arkadiusz Hypki

Abstract: We recently introduced three new parameters that describe the shape of the normalized cumulative radial distribution (nCRD) of the innermost stars in globular clusters and trace the clusters dynamical evolution. Here we extend our previous investigations to the case of a large set of Monte Carlo simulations of globular clusters, started from a broad range of initial conditions. All the models are… ▽ More We recently introduced three new parameters that describe the shape of the normalized cumulative radial distribution (nCRD) of the innermost stars in globular clusters and trace the clusters dynamical evolution. Here we extend our previous investigations to the case of a large set of Monte Carlo simulations of globular clusters, started from a broad range of initial conditions. All the models are analyzed at the same age of 13 Gyr, when they have reached different evolutionary phases. The sample of models is well representative of the structural properties of the observed population of Galactic globular clusters. We confirm that the three nCRD parameters are powerful tools to distinguish systems in early stages of dynamical evolution, from those that already experienced core collapse. They might also help disentangle clusters hosting a low-mass intermediate-mass black hole of a few hundred solar masses, from cases with large concentrations of dark remnants in their centers. With respect to other dynamical indicators, the nCRD parameters offer the advantage of being fully empirical and easier to measure from observational data. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJ, 15 pages, 10 figures, 1 table

MOCCA: Dynamics and evolution of binary stars of multiple stellar populations in tidally filling and underfilling globular star clusters

Authors: Arkadiusz Hypki, Mirek Giersz, Jongsuk Hong, Agostino Leveque, Abbas Askar, Diogo Belloni, Magdalena Otulakowska-Hypka

Abstract: We present an upgraded version of the \MOCCA code for the study of dynamical evolution of globular clusters (GCs) and its first application to the study of evolution of multiple stellar populations. We explore initial conditions spanning different structural parameters for the first (FG) and second generation of stars (SG) and we analyze their effect on the binary dynamics and survival. Here, we f… ▽ More We present an upgraded version of the \MOCCA code for the study of dynamical evolution of globular clusters (GCs) and its first application to the study of evolution of multiple stellar populations. We explore initial conditions spanning different structural parameters for the first (FG) and second generation of stars (SG) and we analyze their effect on the binary dynamics and survival. Here, we focus on the number ratio of FG and SG binaries, its spatial variation, and the way their abundances are affected by various cluster initial properties. We find that present-day SG stars are more abundant in clusters that were initially tidally filling. Conversely, FG stars stay more abundant in clusters that were initially tidally underfilling. We find that the ratio between binary fractions is not affected by the way we calculate these fractions (e.g. only main-sequence binaries (MS) or observational binaries, i.e. MS stars $> 0.4 M_{\odot}$ mass ratios $> 0.5$). This implies that the MS stars themselves are a very good proxy for probing entire populations of FG and SG. We also discuss how it relates to the observations of Milky Way GCs. We show that \MOCCA models are able to reproduce the observed range of SG fractions for Milky Way GCs for which we know these fractions. We show how the SG fractions depend on the initial conditions and provide some constraints for the initial conditions to have more numerous FG or SG stars at the Hubble time. △ Less

Submitted 3 October, 2022; v1 submitted 11 May, 2022; originally announced May 2022.

Comments: Acceped to MNRAS

MOCCA-SURVEY Database II -- Properties of Intermediate Mass Black Holes escaping from star clusters

Authors: Konrad Maliszewski, Mirek Giersz, Dorota Gondek-Rosińska, Abbas Askar, Arkadiusz Hypki

Abstract: In this work we investigate properties of intermediate-mass black holes (IMBHs) that escape from star clusters due to dynamical interactions. The studied models were simulated as part of the preliminary second survey carried out using the MOCCA code (MOCCA-SURVEY Database II), which is based on the Monte Carlo N-body method and does not include gravitational wave recoil kick prescriptions of the b… ▽ More In this work we investigate properties of intermediate-mass black holes (IMBHs) that escape from star clusters due to dynamical interactions. The studied models were simulated as part of the preliminary second survey carried out using the MOCCA code (MOCCA-SURVEY Database II), which is based on the Monte Carlo N-body method and does not include gravitational wave recoil kick prescriptions of the binary black hole merger product. We have found that IMBHs are more likely to be formed and ejected in models where both initial central density and central escape velocities have high values. Most of our studied objects escape in a binary with another black hole (BH) as their companion and have masses between $100$ and $140\: M_{\odot}$. Escaping IMBHs tend to build-up mass most effectively through repeated mergers in a binary with BHs due to gravitational wave emission. Binaries play a key role in their ejection from the system as they allow these massive objects to gather energy needed for escape. The binaries in which IMBHs escape tend to have very high binding energy at the time of escape and the last interaction is strong but does not involve a massive intruder. These IMBHs gain energy needed to escape the cluster gradually in successive dynamical interactions. We present specific examples of the history of IMBH formation and escape from star cluster models. We also discuss the observational implications of our findings as well as the potential influence of the gravitational wave recoil kicks on the process. △ Less

Submitted 18 June, 2022; v1 submitted 17 November, 2021; originally announced November 2021.

Comments: 11 pages, 11 figures, accepted for publish to MNRAS

Testing the presence of a dormant black hole inside HR 6819

Authors: A. Romagnolo, A. Olejak, A. Hypki, G. Wiktorowicz, K. Belczynski

Abstract: HR 6819 was recently reported to be a triple system with a dormant black hole (BH). The inner binary system was defined as a star estimated to be at the end of its main sequence and a dormant BH. As the inner binary is not resolved, the third component may actually just be spatially coinciding with the inner binary. In this study we test whether the system inner binary can be reconstructed using t… ▽ More HR 6819 was recently reported to be a triple system with a dormant black hole (BH). The inner binary system was defined as a star estimated to be at the end of its main sequence and a dormant BH. As the inner binary is not resolved, the third component may actually just be spatially coinciding with the inner binary. In this study we test whether the system inner binary can be reconstructed using the isolated binary evolution in the Galactic field or through the dynamical evolution within globular star clusters. Our goal is to understand the formation of the HR 6819 inner binary. To simulate the inner binary evolution we assumed that the influence of the third body on the inner binary is negligible. We created synthetic populations of BH main-sequence binaries for the Galactic disc and the Galactic globular clusters to compare to the reported parameters of the HR 6819 inner binary. We have adopted very optimistic input physics, in terms of common envelope evolution and BH formation, for the formation of binaries similar to the reported inner HR 6819 binary. Despite our optimistic assumptions we cannot form systems like the inner HR 6819 binary in globular clusters. Even with our extreme assumptions, the formation of an HR 6819-like binary in the Galactic field population is not expected. We argue that if a dormant BH actually exists in the reported configuration inside HR 6819, its presence cannot easily be explained by our models based on isolated and dynamical binary evolution. △ Less

Submitted 9 September, 2022; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 11 pages

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

Dynamical Evolution of Multiple-Population Globular Clusters

Authors: E. Vesperini, J. Hong, M. Giersz, A. Hypki

Abstract: We have carried out a set of Monte Carlo simulations to study a number of fundamental aspects of the dynamical evolution of multiple stellar populations in globular clusters with different initial masses, fractions of second generation (2G) stars, and structural properties. Our simulations explore and elucidate: 1) the role of early and long-term dynamical processes and stellar escape in the evolu… ▽ More We have carried out a set of Monte Carlo simulations to study a number of fundamental aspects of the dynamical evolution of multiple stellar populations in globular clusters with different initial masses, fractions of second generation (2G) stars, and structural properties. Our simulations explore and elucidate: 1) the role of early and long-term dynamical processes and stellar escape in the evolution of the fraction of 2G stars and the link between the evolution of the fraction of 2G stars and various dynamical parameters; 2) the link between the fraction of 2G stars inside the cluster and in the population of escaping stars during a cluster's dynamical evolution; 3) the dynamics of the spatial mixing of the first-generation (1G) and 2G stars and the details of the structural properties of the two populations as they evolve toward mixing; 4) the implications of the initial differences between the spatial distribution of 1G and 2G stars for the evolution of the anisotropy in the velocity distribution and the expected radial profile of the 1G and 2G anisotropy for clusters at different stages of their dynamical history; 5) the variation of the degree of energy equipartition of the 1G and the 2G populations as a function of the distance from the cluster's centre and the cluster's evolutionary phase. △ Less

Submitted 2 February, 2021; originally announced February 2021.

Comments: 16 pages, 19 figures. Accepted for publication in MNRAS

MOCCA SURVEY Database I: Binary Black Hole Mergers from Globular Clusters with Intermediate Mass Black Holes

Authors: Jongsuk Hong, Abbas Askar, Mirek Giersz, Arkadiusz Hypki, Suk-Jin Yoon

Abstract: The dynamical formation of black hole binaries in globular clusters that merge due to gravitational waves occurs more frequently in higher stellar density. Meanwhile, the probability to form intermediate mass black holes (IMBHs) also increases with the density. To explore the impact of the formation and growth of IMBHs on the population of stellar mass black hole binaries from globular clusters, w… ▽ More The dynamical formation of black hole binaries in globular clusters that merge due to gravitational waves occurs more frequently in higher stellar density. Meanwhile, the probability to form intermediate mass black holes (IMBHs) also increases with the density. To explore the impact of the formation and growth of IMBHs on the population of stellar mass black hole binaries from globular clusters, we analyze the existing large survey of Monte-Carlo globular cluster simulation data (MOCCA SURVEY Database I). We show that the number of binary black hole mergers agrees with the prediction based on clusters' initial properties when the IMBH mass is not massive enough or the IMBH seed forms at a later time. However, binary black hole formation and subsequent merger events are significantly reduced compared to the prediction when the present-day IMBH mass is more massive than $\sim10^4 \rm M_{\odot}$ or the present-day IMBH mass exceeds about 1 per cent of cluster's initial total mass. By examining the maximum black hole mass in the system at the moment of black hole binary escaping, we find that $\sim$ 90 per cent of the merging binary black holes escape before the formation and growth of the IMBH. Furthermore, large fraction of stellar mass black holes are merged into the IMBH or escape as single black holes from globular clusters in cases of massive IMBHs, which can lead to the significant under-population of binary black holes merging with gravitational waves by a factor of 2 depending on the clusters' initial distributions. △ Less

Submitted 25 August, 2020; originally announced August 2020.

Comments: 9 pages, 8 figures, Accepted for publication in MNRAS

MOCCA-SURVEY Database I: Dissolution of tidally filling star clusters harboring black hole subsystem

Authors: Mirek Giersz, Abbas Askar, Long Wang, Arkadiusz Hypki, Agostino Leveque, Rainer Spurzem

Abstract: We investigate the dissolution process of star clusters embedded in an external tidal field and harboring a subsystem of stellar-mass black hole. For this purpose we analyzed the MOCCA models of real star clusters contained in the Mocca Survey Database I. We showed that the presence of a stellar-mass black hole subsystem in tidally filling star cluster can lead to abrupt cluster dissolution connec… ▽ More We investigate the dissolution process of star clusters embedded in an external tidal field and harboring a subsystem of stellar-mass black hole. For this purpose we analyzed the MOCCA models of real star clusters contained in the Mocca Survey Database I. We showed that the presence of a stellar-mass black hole subsystem in tidally filling star cluster can lead to abrupt cluster dissolution connected with the loss of cluster dynamical equilibrium. Such cluster dissolution can be regarded as a third type of cluster dissolution mechanism. We additionally argue that such a mechanism should also work for tidally under-filling clusters with a top-heavy initial mass function. △ Less

Submitted 1 August, 2019; originally announced August 2019.

Comments: 4 pages, 1 figure. Contribution to the proceedings of the IAU Symposium No. 351: "Star Clusters: from the Milky Way to the Early Universe"

MOCCA Survey Database I: Dissolution of tidally filling star clusters harbouring BH subsystems

Authors: Mirek Giersz, Abbas Askar, Long Wang, Arkadiusz Hypki, Agostino Leveque, Rainer Spurzem

Abstract: We investigate the dissolution process for dynamically evolving star clusters embedded in an external tidal field by exploring the MOCCA Survey Database I, with focus on the presence and evolution of a stellar-mass black hole subsystem. We argue that the presence of a black hole subsystem can lead to the dissolution of tidally filling star clusters and this can be regarded as a third type of clust… ▽ More We investigate the dissolution process for dynamically evolving star clusters embedded in an external tidal field by exploring the MOCCA Survey Database I, with focus on the presence and evolution of a stellar-mass black hole subsystem. We argue that the presence of a black hole subsystem can lead to the dissolution of tidally filling star clusters and this can be regarded as a third type of cluster dissolution mechanism (in addition to well-known mechanisms connected with strong mass loss due to stellar evolution and mass loss connected with the relaxation process). This third process is characterized by abrupt cluster dissolution connected with the loss of dynamical equilibrium. The abrupt dissolution is powered by strong energy generation from a stellar-mass black hole subsystem accompanied by tidal stripping. Additionally, we argue that such a mechanism should also work for even tidally under-filling clusters with top-heavy initial mass function. Observationally, star clusters which undergo dissolution powered by the third mechanism would look as a 'dark cluster' i.e. composed of stellar mass black holes surrounded by an expanding halo of luminous stars (Banerjee & Kroupa 2011), and they should be different from 'dark clusters' harbouring intermediate mass black holes as discussed by Askar et al. (2017a). An additional observational consequence of an operation of the third dissolution mechanism should be a larger than expected abundance of free floating black holes in the Galactic halo. △ Less

Submitted 2 August, 2019; v1 submitted 2 April, 2019; originally announced April 2019.

Comments: 14 pages, 14 figures, accepted to MNRAS

Gaia Data Release 2: Observational Hertzsprung-Russell diagrams

Authors: Gaia Collaboration, C. Babusiaux, F. van Leeuwen, M. A. Barstow, C. Jordi, A. Vallenari, D. Bossini, A. Bressan, T. Cantat-Gaudin, M. van Leeuwen, A. G. A. Brown, T. Prusti, J. H. J. de Bruijne, C. A. L. Bailer-Jones, M. Biermann, D. W. Evans, L. Eyer, F. Jansen, S. A. Klioner, U. Lammers, L. Lindegren, X. Luri, F. Mignard, C. Panem, D. Pourbaix , et al. (428 additional authors not shown)

Abstract: We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the select… ▽ More We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies. △ Less

Submitted 13 August, 2018; v1 submitted 25 April, 2018; originally announced April 2018.

Comments: Published in the A&A Gaia Data Release 2 special issue. Tables 2 and A.4 corrected. Tables available at http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/616/A10

Journal ref: A&A 616, A10 (2018)

Gaia Data Release 1. Testing the parallaxes with local Cepheids and RR Lyrae stars

Authors: Gaia Collaboration, G. Clementini, L. Eyer, V. Ripepi, M. Marconi, T. Muraveva, A. Garofalo, L. M. Sarro, M. Palmer, X. Luri, R. Molinaro, L. Rimoldini, L. Szabados, I. Musella, R. I. Anderson, T. Prusti, J. H. J. de Bruijne, A. G. A. Brown, A. Vallenari, C. Babusiaux, C. A. L. Bailer-Jones, U. Bastian, M. Biermann, D. W. Evans, F. Jansen , et al. (566 additional authors not shown)

Abstract: Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by… ▽ More Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018. △ Less

Submitted 1 May, 2017; originally announced May 2017.

Comments: 29 pages, 25 figures. Accepted for publication by A&A

Journal ref: A&A 605, A79 (2017)

Gaia Data Release 1. Open cluster astrometry: performance, limitations, and future prospects

Authors: Gaia Collaboration, F. van Leeuwen, A. Vallenari, C. Jordi, L. Lindegren, U. Bastian, T. Prusti, J. H. J. de Bruijne, A. G. A. Brown, C. Babusiaux, C. A. L. Bailer-Jones, M. Biermann, D. W. Evans, L. Eyer, F. Jansen, S. A. Klioner, U. Lammers, X. Luri, F. Mignard, C. Panem, D. Pourbaix, S. Randich, P. Sartoretti, H. I. Siddiqui, C. Soubiran , et al. (567 additional authors not shown)

Abstract: Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the ast… ▽ More Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs. △ Less

Submitted 3 March, 2017; originally announced March 2017.

Comments: Accepted for publication by A&A. 21 pages main text plus 46 pages appendices. 34 figures main text, 38 figures appendices. 8 table in main text, 19 tables in appendices

Journal ref: A&A 601, A19 (2017)

Formation Mechanisms of IMBH in Globular Clusters

Authors: Mirek Giersz, Nathan Leigh, Arkadiusz Hypki, Abbas Askar, Nora Lützgendorf

Abstract: We very briefly discuss proposed in the literature possible scenarios for intermediate mass black holes (IMBH) formation in globular clusters. We also discuss the results of the MOCCA simulations of about 2000 models (BigSurvey) regarding the distribution of events connected with electromagnetic and gravitational radiations, namely: mass transfer on IMBH, collisions and mergers with IMBH and merge… ▽ More We very briefly discuss proposed in the literature possible scenarios for intermediate mass black holes (IMBH) formation in globular clusters. We also discuss the results of the MOCCA simulations of about 2000 models (BigSurvey) regarding the distribution of events connected with electromagnetic and gravitational radiations, namely: mass transfer on IMBH, collisions and mergers with IMBH and mergers with IMBH due to gravitational radiation. The rates of these events are very small, so their observation is very improbable. △ Less

Submitted 28 July, 2016; originally announced July 2016.

Comments: 4 pages, 1 figure

MOCCA-SURVEY Database I: Is NGC 6535 a Dark Star Cluster Harboring an IMBH?

Authors: Abbas Askar, Paolo Bianchini, Ruggero de Vita, Mirek Giersz, Arkadiusz Hypki, Sebastian Kamann

Abstract: We describe the dynamical evolution of a unique type of dark star cluster model in which the majority of the cluster mass at Hubble time is dominated by an intermediate-mass black hole (IMBH). We analyzed results from about 2000 star cluster models (Survey Database I) simulated using the Monte-Carlo code MOCCA and identified these dark star cluster models. Taking one of these models, we apply the… ▽ More We describe the dynamical evolution of a unique type of dark star cluster model in which the majority of the cluster mass at Hubble time is dominated by an intermediate-mass black hole (IMBH). We analyzed results from about 2000 star cluster models (Survey Database I) simulated using the Monte-Carlo code MOCCA and identified these dark star cluster models. Taking one of these models, we apply the method of simulating realistic "mock observations" by utilizing the COCOA and SISCO codes to obtain the photometric and kinematic observational properties of the dark star cluster model at 12 Gyr. We find that the perplexing Galactic globular cluster NGC 6535 closely matches the observational photometric and kinematic properties of the dark star cluster model presented in this paper. Based on our analysis and currently observed properties of NGC 6535, we suggest that this globular cluster could potentially harbour an IMBH. If it exists, the presence of this IMBH can be detected robustly with proposed kinematic observations of NGC 6535. △ Less

Submitted 4 October, 2016; v1 submitted 27 July, 2016; originally announced July 2016.

Comments: 11 pages, 10 figures, 5 tables, Accepted for publication in MNRAS

Cataclysmic variables in Globular clusters: First results on the analysis of the MOCCA simulations database

Authors: Diogo Belloni, Mirek Giersz, Abbas Askar, Arkadiusz Hypki

Abstract: In this first investigation of the MOCCA database with respect to cataclysmic variables, we found that for models with Kroupa initial distributions, considering the standard value of the efficiency of the common-envelope phase adopted in BSE, no single cataclysmic variable was formed only via binary stellar evolution, i. e., in order to form them, strong dynamical interactions have to take place.… ▽ More In this first investigation of the MOCCA database with respect to cataclysmic variables, we found that for models with Kroupa initial distributions, considering the standard value of the efficiency of the common-envelope phase adopted in BSE, no single cataclysmic variable was formed only via binary stellar evolution, i. e., in order to form them, strong dynamical interactions have to take place. Our results also indicate that the population of cataclysmic variables in globular clusters are, mainly, in the last stage of their evolution and observational selection effects can change drastically the expected number and properties of observed cataclysmic variables. △ Less

Submitted 26 July, 2016; originally announced July 2016.

Comments: 4 pages, 3 figures. Presented at the MODEST 16/Cosmic Lab conference in Bologna, Italy, April 18-22 2016. To be pusblished in Mem. S. A. It. Conference Series

MOCCA-SURVEY database I. Accreting white dwarf binary systems in globular clusters -- I. cataclysmic variables -- present-day population

Authors: Diogo Belloni, Mirek Giersz, Abbas Askar, Nathan Leigh, Arkadiusz Hypki

Abstract: In this paper, which is the first in a series of papers associated with cataclysmic variables and related objects, we introduce the CATUABA code, a numerical machinery written for analysis of the MOCCA simulations, and show some first results by investigating the present-day population of cataclysmic variables in globular clusters. Emphasis was given on their properties and the observational selec… ▽ More In this paper, which is the first in a series of papers associated with cataclysmic variables and related objects, we introduce the CATUABA code, a numerical machinery written for analysis of the MOCCA simulations, and show some first results by investigating the present-day population of cataclysmic variables in globular clusters. Emphasis was given on their properties and the observational selection effects when observing and detecting them. In this work we analysed in this work six models, including three with Kroupa distributions of the initial binaries. We found that for models with Kroupa initial distributions, considering the standard value of the efficiency of the common envelope phase adopted in BSE, no single cataclysmic variable was formed only via binary stellar evolution, i. e., in order to form them, strong dynamical interactions have to take place. We show and explain why this is inconsistent with observational and theoretical results. Our results indicate that the population of cataclysmic variables in globular clusters is, mainly, in the last stage of their evolution and observational selection effects can drastically change the expected number of observed cataclysmic variables. We show that the probability of observing them during the outbursts is extremely small and conclude that the best way of looking for cataclysmic variables in globular clusters is by searching for variabilities during quiescence, instead of during outbursts. For that, one would need a very deep observation which could reach magnitudes $\gtrsim$ 27 mag. Finally, we argue that cataclysmic variables in globular clusters are not necessarily magnetic. △ Less

Submitted 11 August, 2016; v1 submitted 26 July, 2016; originally announced July 2016.

Comments: 21 pages, 5 figures, 2 tables; minor changes made to the text, accepted for publication in MNRAS

MOCCA code for star cluster simulations - VI. Bimodal spatial distribution of blue stragglers

Authors: Arkadiusz Hypki, Mirek Giersz

Abstract: The paper presents an analysis of formation mechanism and properties of spatial distributions of blue stragglers in evolving globular clusters, based on numerical simulations done with the MOCCA code. First, there are presented N-body and MOCCA simulations which try to reproduce the simulations presented by Ferraro (2012). Then, the agreement between N-body and the MOCCA code is shown. Finally, we… ▽ More The paper presents an analysis of formation mechanism and properties of spatial distributions of blue stragglers in evolving globular clusters, based on numerical simulations done with the MOCCA code. First, there are presented N-body and MOCCA simulations which try to reproduce the simulations presented by Ferraro (2012). Then, the agreement between N-body and the MOCCA code is shown. Finally, we discuss the formation process of the bimodal distribution. We report that so-called bimodal spatial distribution of blue stragglers is a very transient feature. It is formed for one snapshot in time and it can easily vanish in the next one. Moreover, we show that the radius of avoidance proposed by Ferraro (2012) goes out of sync with the apparent minimum of the bimodal distribution after about two half-mass relaxation times. This finding creates a real challenge for the dynamical clock, which uses this radius to determine the dynamical age of globular clusters. Additionally, the paper discusses a few important problems concerning the apparent visibilities of the bimodal distributions which have to be taken into account while studying the spatial distributions of blue stragglers. △ Less

Submitted 24 April, 2016; originally announced April 2016.

Comments: submitted to MNRAS, comments are welcome

MOCCA code for star cluster simulations - V. Initial globular cluster conditions influence on blue stragglers

Authors: Arkadiusz Hypki, Mirek Giersz

Abstract: The paper presents an analysis of properties of populations of blue stragglers (BSs) in evolving globular clusters, based on numerical simulations done with the MOCCA code for various initial globular clusters conditions. We find that various populations of BSs strongly depend on the initial semi-major axes distributions. With a significant number of compact binaries, the number of evolutionary… ▽ More The paper presents an analysis of properties of populations of blue stragglers (BSs) in evolving globular clusters, based on numerical simulations done with the MOCCA code for various initial globular clusters conditions. We find that various populations of BSs strongly depend on the initial semi-major axes distributions. With a significant number of compact binaries, the number of evolutionary BSs can be also significant. In turn, for semi-major axes distributions preferring binaries with wider orbits, dynamical BSs are the dominant ones. Their formation scenario is very distinct: for wide binaries the number of dynamical interactions is significantly larger. Most interactions are weak and increase only slightly the eccentricities. However, due to a large number of such interactions, the eccentricities of a number of binaries finally get so large that the stars collide. We study how larger initial clusters' concentrations influence the BSs. Besides the expected increase of the number of dynamically created BSs (for denser GCs the probabilities of strong dynamical interactions and collisions are higher), we find that the number of the evolutionary BSs is not affected even by very high initial concentrations. This has a very important implication on observations - it supports the theory that the evolutionary BSs are the result of the unperturbed evolution of the primordial binaries. In addition, the paper presents the evolution of the ratio between the number of BSs in binaries and as single stars R_B/S. For a vast diversity of models, the ratio R_B/S approaches the value ~0.4. Additionally, we identified two subgroups which differ in the initial semi-major axes distributions. The first group starts with a high ratio R_B/S, it decreases with time and settles around 0.4. The second group starts with lower values of the ratio R_B/S ... . (etc., abstract continues) △ Less

Submitted 24 April, 2016; originally announced April 2016.

Comments: submitted to MNRAS, comments are welcome

Gaia archive

Authors: Arkadiusz Hypki, Anthony G. A. Brown

Abstract: The Gaia archive is being designed and implemented by the DPAC Consortium. The purpose of the archive is to maximize the scientific exploitation of the Gaia data by the astronomical community. Thus, it is crucial to gather and discuss with the community the features of the Gaia archive as much as possible. It is especially important from the point of view of the GENIUS project to gather the feedba… ▽ More The Gaia archive is being designed and implemented by the DPAC Consortium. The purpose of the archive is to maximize the scientific exploitation of the Gaia data by the astronomical community. Thus, it is crucial to gather and discuss with the community the features of the Gaia archive as much as possible. It is especially important from the point of view of the GENIUS project to gather the feedback and potential use cases for the archive. This paper presents very briefly the general ideas behind the Gaia archive and presents which tools are already provided to the community. △ Less

Submitted 23 March, 2016; originally announced March 2016.

Comments: Proceedings of the XXXVII Meeting of the Polish Astronomical Society

BEANS - a software package for distributed Big Data analysis

Authors: Arkadiusz Hypki

Abstract: BEANS software is a web based, easy to install and maintain, new tool to store and analyse data in a distributed way for a massive amount of data. It provides a clear interface for querying, filtering, aggregating, and plotting data from an arbitrary number of datasets. Its main purpose is to simplify the process of storing, examining and finding new relations in the so-called Big Data. Creation… ▽ More BEANS software is a web based, easy to install and maintain, new tool to store and analyse data in a distributed way for a massive amount of data. It provides a clear interface for querying, filtering, aggregating, and plotting data from an arbitrary number of datasets. Its main purpose is to simplify the process of storing, examining and finding new relations in the so-called Big Data. Creation of BEANS software is an answer to the growing needs of the astronomical community to have a versatile tool to store, analyse and compare the complex astrophysical numerical simulations with observations (e.g. simulations of the Galaxy or star clusters with the Gaia archive). However, this software was built in a general form and it is ready to use in any other research field or open source software. △ Less

Submitted 23 March, 2016; originally announced March 2016.

Comments: 14 pages, 6 figures, submitted to MNRAS, comments are welcome

MOCCA Code for Star Cluster Simulations - IV. A new Scenario for Intermediate Mass Black Hole Formation in Globular Clusters

Authors: M. Giersz, N. Leigh, A. Hypki, N. Lützgendorf, A. Askar

Abstract: We discuss a new scenario for the formation of intermediate mass black holes in dense star clusters. In this scenario, intermediate mass black holes are formed as a result of dynamical interactions of hard binaries containing a stellar mass black hole, with other stars and binaries. We discuss the necessary conditions to initiate the process of intermediate mass black hole formation and the influe… ▽ More We discuss a new scenario for the formation of intermediate mass black holes in dense star clusters. In this scenario, intermediate mass black holes are formed as a result of dynamical interactions of hard binaries containing a stellar mass black hole, with other stars and binaries. We discuss the necessary conditions to initiate the process of intermediate mass black hole formation and the influence of an intermediate mass black hole on the host global globular cluster properties. We discuss two scenarios for intermediate mass black hole formation. The SLOW and FAST scenarios. They occur later or earlier in the cluster evolution and require smaller or extremely large central densities, respectively. In our simulations, the formation of intermediate mass black holes is highly stochastic. In general, higher formation probabilities follow from larger cluster concentrations (i.e. central densities). We further discuss possible observational signatures of the presence of intermediate mass black holes in globular clusters that follow from our simulations. These include the spatial and kinematic structure of the host cluster, possible radio, X-ray and gravitational wave emissions due to dynamical collisions or mass-transfer and the creation of hypervelocity main sequence escapers during strong dynamical interactions between binaries and an intermediate mass black hole. All simulations discussed in this paper were performed with the MOCCA Monte Carlo code. MOCCA accurately follows most of the important physical processes that occur during the dynamical evolution of star clusters but, as with other dynamical codes, it approximates the dissipative processes connected with stellar collisions and binary mergers. △ Less

Submitted 4 August, 2016; v1 submitted 17 June, 2015; originally announced June 2015.

Comments: 16 pages, 1 table, 24 figures. printed in MNRAS

MOCCA Code for Star Cluster Simulation: Comparison with Optical Observations using COCOA

Authors: Abbas Askar, Mirek Giersz, Wojciech Pych, Arkadiusz Olech, Arkadiusz Hypki

Abstract: We introduce and present preliminary results from COCOA (Cluster simulatiOn Comparison with ObservAtions) code for a star cluster after 12 Gyrs of evolution simulated using the MOCCA code. The COCOA code is being developed to quickly compare results of numerical simulations of star clusters with observational data. We use COCOA to obtain parameters of the projected cluster model. For comparison, a… ▽ More We introduce and present preliminary results from COCOA (Cluster simulatiOn Comparison with ObservAtions) code for a star cluster after 12 Gyrs of evolution simulated using the MOCCA code. The COCOA code is being developed to quickly compare results of numerical simulations of star clusters with observational data. We use COCOA to obtain parameters of the projected cluster model. For comparison, a FITS file of the projected cluster was provided to observers so that they could use their observational methods and techniques to obtain cluster parameters. The results show that the similarity of cluster parameters obtained through numerical simulations and observations depends significantly on the quality of observational data and photometric accuracy. △ Less

Submitted 8 August, 2016; v1 submitted 29 December, 2014; originally announced January 2015.

Comments: 2 pages, 1 figure, Conference proceedings for poster presented at IAUS 312

Journal ref: Star Clusters and Black Holes in Galaxies across Cosmic Time, Proceedings of the International Astronomical Union, IAU Symposium, Volume 312, pp. 262-263 02/2016

Monte Carlo modelling of globular star clusters - many primordial binaries, IMBH formation

Authors: Mirek Giersz, Nathan Leigh, Michael Marks, Arkadiusz Hypki, Abbas Askar

Abstract: We will discuss the evolution of star clusters with an large initial binary fraction, up to 95%. The initial binary population is chosen to follow the invariant orbital-parameter distributions suggested by Kroupa (1995). The Monte Carlo MOCCA simulations of star cluster evolution are compared to the observations of Milone et al. (2012) for photometric binaries. It is demonstrated that the observed… ▽ More We will discuss the evolution of star clusters with an large initial binary fraction, up to 95%. The initial binary population is chosen to follow the invariant orbital-parameter distributions suggested by Kroupa (1995). The Monte Carlo MOCCA simulations of star cluster evolution are compared to the observations of Milone et al. (2012) for photometric binaries. It is demonstrated that the observed dependence on cluster mass of both the binary fraction and the ratio of the binary fractions inside and outside of the half mass radius are well recovered by the MOCCA simulations. This is due to a rapid decrease in the initial binary fraction due to the strong density-dependent destruction of wide binaries described by Marks, Kroupa & Oh (2011). We also discuss a new scenario for the formation of intermediate mass black holes in dense star clusters. In this scenario, intermediate mass black holes are formed as a result of dynamical interactions of hard binaries containing a stellar mass black hole, with other stars and binaries. We will discuss the necessary conditions to initiate the process of intermediate mass black hole formation and the dependence of its mass accretion rate on the global cluster properties. △ Less

Submitted 27 November, 2014; originally announced November 2014.

Comments: 12 pages and 7 figures. Presented at IAU312, to be published as Conference Proceedings

The state of globular clusters at birth II: primordial binaries

Authors: Nathan W. C. Leigh, Mirek Giersz, Michael Marks, Jeremy J. Webb, Arkadiusz Hypki, Craig O. Heinke, Pavel Kroupa, Alison Sills

Abstract: (abridged) In this paper, we constrain the properties of primordial binary populations in Galactic globular clusters using the MOCCA Monte Carlo code for cluster evolution. Our results are compared to the observations of Milone et al. (2012) using the photometric binary populations as proxies for the true underlying distributions, in order to test the hypothesis that the data are consistent with a… ▽ More (abridged) In this paper, we constrain the properties of primordial binary populations in Galactic globular clusters using the MOCCA Monte Carlo code for cluster evolution. Our results are compared to the observations of Milone et al. (2012) using the photometric binary populations as proxies for the true underlying distributions, in order to test the hypothesis that the data are consistent with an universal initial binary fraction near unity and the binary orbital parameter distributions of Kroupa (1995). With the exception of a few possible outliers, we find that the data are to first-order consistent with the universality hypothesis. Specifically, the present-day binary fractions inside the half-mass radius r$_{\rm h}$ can be reproduced assuming either high initial binary fractions near unity with a dominant soft binary component as in the Kroupa distribution combined with high initial densities (10$^4$-10$^6$ M$_{\odot}$ pc$^{-3}$), or low initial binary fractions ($\sim$ 5-10%) with a dominant hard binary component combined with moderate initial densities near their present-day values (10$^2$-10$^3$ M$_{\odot}$ pc$^{-3}$). This apparent degeneracy can be broken using the binary fractions outside r$_{\rm h}$- only high initial binary fractions with a significant soft component combined with high initial densities can contribute to reproducing the observed anti-correlation between the binary fractions outside r$_{\rm h}$ and the total cluster mass. We further illustrate using the simulated present-day binary orbital parameter distributions and the technique introduced in Leigh et al. (2012) that the relative fractions of hard and soft binaries can be used to further constrain the initial cluster density and mass-density relation. Our results favour an initial mass-density relation of the form r$_{\rm h} \propto$ M$_{\rm clus}^α$ with $α<$ 1/3. △ Less

Submitted 8 October, 2014; originally announced October 2014.

Comments: 17 pages; 5 figures; 1 table; accepted for publication in MNRAS

The state of globular clusters at birth: emergence from the gas-embedded phase

Authors: Nathan W. C. Leigh, Mirek Giersz, Jeremy J. Webb, Arkadiusz Hypki, Guido de Marchi, Pavel Kroupa, Alison Sills

Abstract: In this paper, we discuss the origin of the observed correlation between cluster concentration c and present-day mass function (PDMF) slope α reported by De Marchi, Paresce & Pulone. This relation can either be reproduced from universal initial conditions combined with some dynamical mechanism(s) that alter(s) the cluster structure and mass function over time, or it must arise early on in the clus… ▽ More In this paper, we discuss the origin of the observed correlation between cluster concentration c and present-day mass function (PDMF) slope α reported by De Marchi, Paresce & Pulone. This relation can either be reproduced from universal initial conditions combined with some dynamical mechanism(s) that alter(s) the cluster structure and mass function over time, or it must arise early on in the cluster lifetime, such as during the gas-embedded phase of cluster formation. Using a combination of Monte Carlo and N-body models for globular cluster evolution performed with the MOCCA and NBODY6 codes, respectively, we explore a number of dynamical mechanisms that could affect the observed relation. For the range of initial conditions considered here, our results are consistent with an universal initial binary fraction \approx 10% (which does not, however, preclude 100%) and an universal initial stellar mass function resembling the standard Kroupa distribution. Most of the dispersion observed in the c-α relation can be attributed to two-body relaxation and Galactic tides. However, dynamical processes alone could not have reproduced the dispersion in concentration, and we require at least some correlation between the initial concentration and the total cluster mass. We argue that the origin of this trend could be connected to the gas-embedded phase of cluster evolution. △ Less

Submitted 18 November, 2014; v1 submitted 26 September, 2013; originally announced September 2013.

Comments: 16 pages, 5 figures, 2 tables; accepted for publication in MNRAS; edited to match proof corrections

MOCCA code for star cluster simulations - I. Blue Stragglers, first results

Authors: Arkadiusz Hypki, Mirek Giersz

Abstract: We introduce an improved code for simulations of star clusters, called MOCCA. It combines the Monte Carlo method for star cluster evolution and the Fewbody code to perform scattering experiments. The Fewbody was added in order to track more precisely dynamical interactions between objects which can lead to creations of various exotic objects observed in the star clusters, like Blue Stragglers Star… ▽ More We introduce an improved code for simulations of star clusters, called MOCCA. It combines the Monte Carlo method for star cluster evolution and the Fewbody code to perform scattering experiments. The Fewbody was added in order to track more precisely dynamical interactions between objects which can lead to creations of various exotic objects observed in the star clusters, like Blue Stragglers Stars (BSS). The MOCCA code is currently one of the most advanced codes for simulating real size star clusters. It follows the star cluster evolution closely to N-body codes but is much faster. We show that the MOCCA code is able to follow the evolution of BSS with details. It is a suitable tool to perform full scale evolution of real star clusters and detail comparison with observations of exotic star cluster objects like BSS. This paper is the first one of the series of papers about properties of BSS in star clusters. This type of stars is particularly interesting today, because by studying them one can get important constrains on a link between the stellar and dynamical evolution of star clusters. We discuss here first results concerning BSS for an arbitrary chosen test model. We investigate properties of BSS which characterize different channels of formation like masses, semi-major axes, eccentricities, and orbital periods. We show how BSS from different channels change their types, and discuss initial and final positions of BSS, their bimodal distribution in the star cluster, lifetimes and more. △ Less

Submitted 28 July, 2012; originally announced July 2012.

Comments: 25 pages, 25 figures

MOCCA Code for Star Cluster Simulations - II. Comparison with N-body Simulations

Authors: Mirek Giersz, Douglas C. Heggie, Jarrod Hurley, Arkadiusz Hypki

Abstract: We describe a major upgrade of a Monte Carlo code which has previously been used for many studies of dense star clusters. We outline the steps needed in order to calibrate the results of the new Monte Carlo code against $N$-body simulations for large $N$ systems, up to $N=200000$. The new version of the Monte Carlo code (called MOCCA), in addition to the features of the old version, incorporates t… ▽ More We describe a major upgrade of a Monte Carlo code which has previously been used for many studies of dense star clusters. We outline the steps needed in order to calibrate the results of the new Monte Carlo code against $N$-body simulations for large $N$ systems, up to $N=200000$. The new version of the Monte Carlo code (called MOCCA), in addition to the features of the old version, incorporates the direct Fewbody integrator (Fregeau et al. 2004) for three- and four-body interactions, and a new treatment of the escape process based on Fukushige & Heggie (2000). Now stars which fulfil the escape criterion are not removed immediately, but can stay in the system for a certain time which depends on the excess of the energy of a star above the escape energy. They are called potential escapers. With the addition of the Fewbody integrator the code can follow all interaction channels which are important for the rate of creation of various types of objects observed in star clusters, and ensures that the energy generation by binaries is treated in a manner similar to the $N$-body model. There are at most three new parameters which have to be adjusted against $N$-body simulations for large $N$: two (or one, depending on the chosen approach) connected with the escape process, and one responsible for the determination of the interaction probabilities. The values adopted for the free parameters have at most a weak dependence on $N$. They allow MOCCA to reproduce $N$-body results with reasonable precision, not only for the rate of cluster evolution and the cluster mass distribution, but also for the detailed distributions of mass and binding energy of binaries. Additionally, the code can follow the rate of formation of blue stragglers and black hole - black hole binaries. △ Less

Submitted 17 December, 2013; v1 submitted 29 December, 2011; originally announced December 2011.

Comments: 15 pages, 24 figures. Matched to published version in MNRAS