Source author record

M. Giersz

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

13 published item(s)

preprint2022arXiv

A Monte Carlo study of early gas expulsion and evolution of star clusters: new simulations with the MOCCA code in the AMUSE framework

We introduce a new prescription for the evolution of globular clusters (GCs) during the initial embedded gas phase into a Monte Carlo method. With a simplified version of the Monte Carlo MOCCA code embedded in the AMUSE framework, we study the survival of GCs after the removal of primordial gas. We first test our code and show that our results for the evolution of mass and Lagrangian radii are in good agreement with those obtained with N-body simulations. The Monte Carlo code enables a more rapid exploration of the evolution of systems with a larger number of stars than N-body simulations. We have carried out a new survey of simulations to explore the evolution of globular clusters with up to $N = 500000$ stars for a range of different star formation efficiencies and half-mass radii. Our study shows the range of initial conditions leading to the clusters' dissolution and those for which the clusters can survive this early evolutionary phase.

preprint2022arXiv

MOCCA-Survey Database: Extra Galactic Globular Clusters. II. Milky Way and Andromeda

A comprehensive study of the co-evolution of globular cluster systems (GCS) in galaxies requires the ability to model both the large scale dynamics (0.01 - 10 kpc) regulating their orbital evolution, and the small scale dynamics (sub-pc - AU) regulating the internal dynamics of each globular cluster (GC). In this work we present a novel method that combine semi-analytic models of GCS with fully self-consistent Monte Carlo models to simultaneously evolve large GCSs. We use the population synthesis code MASinGa and the MOCCA-Survey Database I to create synthetic GC populations aimed at representing the observed features of GCs in the Milky Way (MW) and Andromeda (M31). Our procedure enables us to recover the spatial and mass distribution of GCs in such galaxies, and to constrain the amount of mass that GCs left either in the halo as dispersed debris, or in the galactic centre, where they can contribute to the formation of a nuclear star cluster (NSC) and can bring stellar and possibly intermediate mass black holes there. The final masses reported by our simulations are of a few order of magnitudes smaller than the observed values. These differences show that mass build-up of a NSC and central BHs in galaxies like MW and M31 cannot be solely explained by the infalling GC scenario. This build-up is likely to depend on the interplay between interactions and mergers of infalling GCs and gas. The latter can contribute to both in-situ star formation in the NSC and growth of the central BH.

preprint2022arXiv

The impact of stellar evolution on rotating star clusters: the gravothermal-gravogyro catastrophe and the formation of a bar of black holes

We present results from a suite of eight direct N-body simulations, performed with \textsc{Nbody6++GPU}, representing realistic models of rotating star clusters with up to $1.1\times 10^5$ stars. Our models feature primordial (hard) binaries, a continuous mass spectrum, differential rotation, and tidal mass loss induced by the overall gravitational field of the host galaxy. We explore the impact of rotation and stellar evolution on the star cluster dynamics. In all runs for rotating star clusters we detect a previously predicted mechanism: an initial phase of violent relaxation followed by the so-called gravogyro catastrophe. We find that the gravogyro catastrophe reaches a finite amplitude, which depends in strength on the level of the bulk rotation, and then levels off. After this phase the angular momentum is transferred from high-mass to low-mass particles in the cluster (both stars and compact objects). Simultaneously, the system becomes gravothermally unstable and collapses, thus undergoing the so-called gravothermal-gravogyro catastrophe. Comparing models with and without stellar evolution, we find an interesting difference. When stellar evolution is not taken into account, the whole process proceeds at a faster pace. The population of heavy objects tend to form a triaxial structure that rotates in the cluster centre. When stellar evolution is taken into account, we find that such a {\it rotating bar} is populated by stellar black holes and their progenitors. The triaxial structure becomes axisymmetric over time, but we also find that the models without stellar evolution suffer repeated gravogyro catastrophes as sufficient angular momentum and mass are removed by the tidal field.

preprint2021arXiv

Dynamical Evolution of Multiple-Population Globular Clusters

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.

preprint2021arXiv

Preparing the next gravitational million-body simulations: Evolution of single and binary stars in Nbody6++GPU, MOCCA and McLuster

We present the implementation of updated stellar evolution recipes in the codes \texttt{Nbody6++GPU, MOCCA} and \texttt{McLuster}. We test them through numerical simulations of star clusters containing $1.1\times 10^5$ stars (with $2.0\times 10^4$ in primordial hard binaries) performing high-resolution direct $N$-body (\texttt{Nbody6++GPU}) and Monte-Carlo (\texttt{MOCCA}) simulations to an age of 10~Gyr. We compare models implementing either delayed or core-collapse supernovae mechanisms, a different mass ratio distribution for binaries, and white dwarf natal kicks enabled/disabled. Compared to \texttt{Nbody6++GPU}, the \texttt{MOCCA} models appear to be denser, with a larger scatter in the remnant masses, and a lower binary fraction on average. The \texttt{MOCCA} models produce more black holes (BHs) and helium white dwarfs (WDs), whilst \texttt{Nbody6++GPU} models are characterised by a much larger amount of WD-WD binaries. The remnant kick velocity and escape speed distributions are similar for the BHs and neutron stars (NSs), and some NSs formed via electron-capture supernovae, accretion-induced collapse or merger-induced collapse escape the cluster in all simulations. The escape speed distributions for the WDs, on the other hand, are very dissimilar. We categorise the stellar evolution recipes available in \texttt{Nbody6++GPU}, \texttt{MOCCA} and \texttt{Mcluster} into four levels: the one implemented in previous \texttt{Nbody6++GPU} and \texttt{MOCCA} versions (\texttt{level A}), state-of-the-art prescriptions (\texttt{level B}), some in a testing phase (\texttt{level C}), and those that will be added in future versions of our codes.

preprint2020arXiv

The evolutionary roads leading to low effective spins, high black hole masses, and O1/O2 rates of LIGO/Virgo binary black holes

All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported from the O1/O2 runs have near zero effective spins. There are only three potential explanations of this fact. If the BH spin magnitudes are large then (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Or, (iii) the BH spin magnitudes are small. We test the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, include revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can match simultaneously the observed BH-BH merger rate density, BH masses, and effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be a key in better reproducing the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely overestimated if the merger GW170729 hosts a BH more massive than 50 Msun. We also estimate rate of BH-NS mergers from recent LIGO/Virgo observations. Our updated models of BH-BH, BH-NS and NS-NS mergers are now publicly available at www.syntheticuniverse.org.

preprint2016arXiv

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

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.

preprint2014arXiv

M4 Core Project with HST - III. Search for variable stars in the primary field

We present the results of a photometric search for variable stars in the core of the Galactic globular cluster M4. The input data are a large and unprecedented set of deep Hubble Space Telescope WFC3 images (large program GO-12911; 120 orbits allocated), primarily aimed at probing binaries with massive companions by detecting their astrometric wobbles. Though these data were not optimised to carry out a time-resolved photometric survey, their exquisite precision, spatial resolution and dynamic range enabled us to firmly detect 38 variable stars, of which 20 were previously unpublished. They include 19 cluster-member eclipsing binaries (confirming the large binary fraction of M4), RR Lyrae, and objects with known X-ray counterparts. We improved and revised the parameters of some among published variables.

preprint2014arXiv

The M4 Core Project with HST -- II. Multiple Stellar Populations at the Bottom of the Main Sequence

The M4 Core Project with HST is designed to exploit the Hubble Space Telescope to investigate the central regions of M4, the Globular Cluster closest to the Sun. In this paper we combine optical and near-infrared photometry to study multiple stellar populations in M4. We detected two sequences of M-dwarfs containing ~38% (MS_I) and ~62% (MS_II) of MS stars below the main-sequence (MS) knee. We compare our observations with those of NGC2808, which is the only other GCs where multiple MSs of very low-mass stars have been studied to date. We calculate synthetic spectra for M-dwarfs, assuming the chemical composition mixture inferred from spectroscopic studies of stellar populations along the red giant branch, and different Helium abundances, and we compare predicted and observed colors. Observations are consistent with two populations, one with primordial abundance and another with enhanced nitrogen and depleted oxygen.

preprint2013arXiv

The M4 Core Project with HST --- I. Overview and First-Epoch

We present an overview of the ongoing Hubble Space Telescope large program GO-12911. The program is focused on the core of M4, the nearest Galactic globular cluster, and the observations are designed to constrain the number of binaries with massive companions (black holes, neutron stars, or white dwarfs) by measuring the ``wobble'' of the luminous (main-sequence) companion around the center of mass of the pair, with an astrometric precision of ~50 micro-arcseconds. The high spatial resolution and stable medium-band PSFs of WFC3/UVIS will make these measurements possible. In this work we describe: (i) the motivation behind this study, (ii) our observing strategy, (iii) the many other investigations enabled by this unique data set, and which of those our team is conducting, and (iv) a preliminary reduction of the first-epoch data-set collected on October 10, 2012.

preprint2010arXiv

Compact Binaries in Star Clusters I - Black Hole Binaries Inside Globular Clusters

We study the compact binary population in star clusters, focusing on binaries containing black holes, using a self-consistent Monte Carlo treatment of dynamics and full stellar evolution. We find that the black holes experience strong mass segregation and become centrally concentrated. In the core the black holes interact strongly with each other and black hole-black hole binaries are formed very efficiently. The strong interactions, however, also destroy or eject the black hole-black hole binaries. We find no black hole-black hole mergers within our simulations but produce many hard escapers that will merge in the galactic field within a Hubble time. We also find several highly eccentric black hole-black hole binaries that are potential LISA sources, suggesting that star clusters are interesting targets for space-based detectors. We conclude that star clusters must be taken into account when predicting compact binary population statistics.

preprint2010arXiv

Compact Binaries in Star Clusters II - Escapers and Detection Rates

We use a self-consistent Monte Carlo treatment of stellar dynamics to investigate black hole binaries that are dynamically ejected from globular clusters to determine if they will be gravitational wave sources. We find that many of the ejected binaries have initially short periods and will merge within a Hubble time due to gravitational wave radiation. Thus they are potential sources for ground-based gravitational wave detectors. We estimate the yearly detection rate for current and advanced ground-based detectors and find a modest enhancement over the rate predicted for binaries produced by pure stellar evolution in galactic fields. We also find that many of the ejected binaries will pass through the longer wavelength Laser Interferometer Space Antenna (LISA) band and may be individually resolvable. We find a low probability that the Galaxy will contain a binary in the LISA band during its three-year mission. Some such binaries may, however, be detectable at Mpc distances implying that there may be resolvable stellar-mass LISA sources beyond our Galaxy. We conclude that globular clusters have a significant effect on the detection rate of ground-based detectors and may produce interesting LISA sources in local group galaxies.

preprint2003arXiv

A stochastic Monte Carlo approach to model real star cluster evolution, III. Direct integrations of three- and four-body interactions

Spherically symmetric equal mass star clusters containing a large amount of primordial binaries are studied using a hybrid method, consisting of a gas dynamical model for single stars and a Monte Carlo treatment for relaxation of binaries and the setup of close resonant and fly-by encounters of single stars with binaries and binaries with each other (three- and four-body encounters). What differs from our previous work is that each encounter is being integrated using a highly accurate direct few-body integrator which uses regularized variables. Hence we can study the systematic evolution of individual binary orbital parameters (eccentricity, semi-major axis) and differential and total cross sections for hardening, dissolution or merging of binaries (minimum distance) from a sampling of several ten thousands of scattering events as they occur in real cluster evolution including mass segregation of binaries, gravothermal collapse and reexpansion, binary burning phase and ultimately gravothermal oscillations. For the first time we are able to present empirical cross sections for eccentricity variation of binaries in close three- and four-body encounters. It is found that a large fraction of three-body and four-body encounters results in merging. Previous cross sections obtained by Spitzer and Gao for strong encounters can be reproduced, while for weak encounters non-standard processes like formation of hierarchical triples occur.