BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Mirek Giersz

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

ResearcherUnclaimed source record
astro-ph.GAastro-ph.SRastro-ph.HEastro-phastro-ph.IMphysics.comp-phastro-ph.EP

Catalog footprint

What is connected

33works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

33 published item(s)

preprint2025arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Forming short period sub-stellar companions in 47 Tucanae -- II. Analytic expressions for the orbital evolution of planets in dense environments

Short period, massive planets, known as hot Jupiters (HJs), have been discovered around $\sim 1$ percent of local field stars. The inward migration necessary to produce HJs may be `low eccentricity', due to torques in the primordial disc, or `high eccentricity' (HEM). The latter involves exciting high orbital eccentricity, allowing sufficiently close passages with the host star to raise circularising tides in the planet. We present an analytic framework for quantifying the role of dynamical encounters in high density environments during HEM. We show that encounters can enhance or suppress HEM, depending on the local stellar density and the initial semi-major axis $a_0$. For moderate densities, external perturbations can excite large eccentricities that allow a planet to circularise over the stellar lifetime. At extremely high densities, these perturbations can instead result in tidal disruption of the planet, thus yielding no HJ. This may explain the apparent excess of HJs in M67 compared with their local field star abundance versus their apparent deficit in 47 Tuc. Applying our analytic framework, we demonstrate that for an initial massive planet population similar to the field, the expected HJ occurrence rate in 47 Tuc is $f_\mathrm{HJ}=2.2\times 10^{-3}$, which remains consistent with present constraints. Future large (sample sizes $\gtrsim 10^5$) or sensitive transit surveys of stars in globular clusters are required to refute the hypothesis that the initial planet population is similar to the solar neighbourhood average. Non-detection in such surveys would have broad consequences for planet formation theory, implying planet formation rates in globular clusters must be suppressed across a wide range of $a_0$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

MOCCA-SURVEY Database -- I. Tidal disruption events of white dwarfs in globular clusters and young massive clusters

We exploit more than 1000 star cluster Monte Carlo models from the MOCCA-SURVEY Database I, to infer the local rate density of white dwarf (WD) tidal disruption events (TDEs) in globular clusters (GCs) and young massive clusters (YMCs). We suggest that the WD TDE rate for GCs and YMCs in the local universe is $\sim 90$-$500~{\rm yr}^{-1}~{\rm Gpc}^{-3}$, with 90 % of WD TDEs occurring in GCs. The total WD TDE rate density is $\sim 9$-$50$ times larger than estimated before. Our results show that thermonuclear explosions induced by WD TDEs can be observed at a rate of $\sim 100$-$550~{\rm yr}^{-1}$ by the next generation optical surveys, such as the Legacy Survey of Space & Time by the Vera C. Rubin Observatory. We also find that massive WDs are preferentially disrupted due to mass segregation, and that 20 % of exploding WDs have $\gtrsim 1.0 M_\odot$ despite of the small population of such WDs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Black hole mergers in compact star clusters and massive black hole formation beyond the mass-gap

We present direct N-body simulations, carried out with Nbody6++GPU, of young and compact low metallicity star clusters with $1.1\times 10^5$ stars, a velocity dispersion of $\sim$ 10 $\mathrm{km\,s^{-1}}$, a half mass radius $R_h=0.6$ pc, and a binary fraction of $10\%$ including updated evolution models for stellar winds and pair-instability supernovae (PISNe). Within the first tens of megayears of evolution, each cluster hosts several black hole (BH) merger events which nearly cover the complete mass range of primary and secondary BH masses for current LIGO/Virgo/Kagra gravitational wave detections. The importance of gravitational recoil is estimated statistically. We present several possible formation paths of massive BHs above the assumed lower PISNe mass-gap limit ($45 M_\odot$) into the intermediate-mass BH (IMBH) regime ($> 100 M_\odot$) which include collisions of stars and BHs as well as the direct collapse of stellar merger remnants with low mass cores. The stellar evolution updates result in the early formation of higher mass stellar BHs than for the previous model. The resulting higher collision rates with massive stars support the rapid formation of massive BHs. For models assuming a high accretion efficiency for star-BH mergers, we present a first-generation formation scenario for GW190521-like events, a merger of two BHs in the PISN mass-gap, which is dominated by star-BH mergers. This IMBH formation path is independent of gravitational recoil and therefore conceivable in dense stellar systems with low escape velocities. One simulated cluster even forms an IMBH binary (153$M_\odot$,173$M_\odot$) which is expected to merge within a Hubble time.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the absence of symbiotic stars in globular clusters

Even though plenty of symbiotic stars (SySts) have been found in the Galactic field and nearby galaxies, not a single one has ever been confirmed in a Galactic globular cluster (GC). We investigate the lack of such systems in GCs for the first time by analysing 144 GC models evolved with the MOCCA code, which have different initial properties and are roughly representative of the Galactic GC population. We focus here on SySts formed through the wind-accretion channel, which can be consistently modelled in binary population synthesis codes. We found that the orbital periods of the majority of such SySts are sufficiently long (${\gtrsim10^3}$ d) so that, for very dense GC models, dynamical interactions play an important role in destroying their progenitors before the present day (${\sim11-12}$ Gyr). In less dense GC models, some SySts are still predicted to exist. However, these systems tend to be located far from the central parts (${\gtrsim70}$ per cent are far beyond the half-light radius) and are sufficiently rare (${\lesssim1}$ per GC per Myr), which makes their identification rather difficult in observational campaigns. We propose that future searches for SySts in GCs should be performed in the outskirts of nearby low-density GCs with sufficiently long half-mass relaxation times and relatively large Galactocentric distances. Finally, we obtained spectra of the candidate proposed in $ω$ Cen (SOPS IV e-94) and showed that this object is most likely not a SySt.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Are most Cataclysmic Variables in Globular Clusters dynamically formed?

We have been investigating populations of cataclysmic variables (CVs) in a set of more than 300 globular cluster (GC) models evolved with the MOCCA code. One of the main questions we have intended to answer is whether most CVs in GCs are dynamically formed or not. Contrary to what has been argued for a long time, we found that dynamical destruction of primordial CV progenitors is much stronger in GCs than dynamical formation of CVs. In particular, we found that, on average, the detectable CV population is predominantly composed of CVs formed via a typical common envelope phase (> 70 per cent). However, core-collapsed models tend to have higher fractions of bright CVs than non-core-collapsed ones, which suggests then that the formation of CVs is indeed slightly favoured through strong dynamical interactions in core-collapsed GCs, due to the high stellar densities in their cores.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamical Formation of Cataclysmic Variables in Globular Clusters

The formation and evolution of X-ray sources in globular clusters is likely to be affected by the cluster internal dynamics and the stellar interactions in the cluster dense environment.Several observational studies have revealed a correlation between the number of X-ray sources and the stellar encounter rate and provided evidence of the role of dynamics in the formation of X-ray binaries. We have performed a survey of Monte-Carlo simulations aimed at exploring the connection between the dynamics and formation of cataclysmic variables (CVs) and the origin of the observed correlation between the number of these objects, $N_{\rm cv}$, and the stellar encounter rate, $Γ$.The results of our simulations show a correlation between $N_{\rm cv}$ and $Γ$ as found in observational data, illustrate the essential role played by dynamics, and shed light on the dynamical history behind this correlation. CVs in our simulations are more centrally concentrated than single stars with masses close to those of turn-off stars, although this trend is stronger for CVs formed from primordial binaries undergoing exchange encounters, which include a population of more massive CVs absent in the group of CVs formed from binaries not suffering any component exchange.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Formation Mechanisms of IMBH in Globular Clusters

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

MOCCA-SURVEY Database I: Coalescing Binary Black Holes Originating From Globular Clusters

In this first of a series of papers, we utilize results for around two thousand star cluster models simulated using the MOCCA code for star cluster evolution (Survey Database I) to determine the astrophysical properties and local merger rate densities for coalescing binary black holes (BBHs) originating from globular clusters (GCs). We extracted information for all coalescing BBHs that escape the cluster models and subsequently merge within a Hubble time along with BBHs that are retained in our GC models and merge inside the cluster via gravitational wave (GW) emission. By obtaining results from a substantial number of realistic star cluster models that cover different initial parameters, we have an extremely large statistical sample of BBHs with stellar mass and massive stellar BH ($\lesssim 100M_{\odot}$) components that merge within a Hubble time. Using this data, we estimate local merger rate densities for these BBHs originating from GCs to be at least 5.4 ${\rm Gpc}^{-3}\,{\rm yr}^{-1}$

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

MOCCA-SURVEY database I. Accreting white dwarf binary systems in globular clusters -- II. Cataclysmic variables -- progenitors and population at birth

This is the second in a series of papers associated with cataclysmic variables (CVs) and related objects, formed in a suite of simulations for globular cluster evolution performed with the MOCCA Monte Carlo code. We study the properties of our simulated CV populations throughout the entire cluster evolution. We find that dynamics extends the range of binary CV progenitor properties, causing CV formation from binary progenitors that would otherwise not become CVs. The CV formation rate in our simulations can be separated into two regimes: an initial burst ($\lesssim$ 1 Gyr) connected with the formation of the most massive white dwarfs, followed by a nearly constant formation rate. This result holds for all models regardless of the adopted initial conditions, even when most CVs form dynamically. Given the cluster age-dependence of CV properties, we argue that direct comparisons to observed Galactic field CVs could be misleading, since cluster CVs can be up to 4 times older than their field counterparts. Our results also illustrate that, due mainly to unstable mass transfer, some CVs that form in our simulations are destroyed before the present-day. Finally, some field CVs might have originated from globular clusters, as found in our simulations, although the fraction of such escapers should be small relative to the entire Galactic field CV population.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The DRAGON simulations: globular cluster evolution with a million stars

Introducing the DRAGON simulation project, we present direct $N$-body simulations of four massive globular clusters (GCs) with $10^6$ stars and 5$\%$ primordial binaries at a high level of accuracy and realism. The GC evolution is computed with NBODY6++GPU and follows the dynamical and stellar evolution of individual stars and binaries, kicks of neutron stars and black holes, and the effect of a tidal field. We investigate the evolution of the luminous (stellar) and dark (faint stars and stellar remnants) GC components and create mock observations of the simulations (i.e. photometry, color-magnitude diagrams, surface brightness and velocity dispersion profiles). By connecting internal processes to observable features we highlight the formation of a long-lived 'dark' nuclear subsystem made of black holes (BHs), which results in a two-component structure. The inner core is dominated by the BH subsystem and experiences a core collapse phase within the first Gyr. It can be detected in the stellar (luminous) line-of-sight velocity dispersion profiles. The outer extended core - commonly observed in the (luminous) surface brightness profiles - shows no collapse features and is continuously expanding. We demonstrate how a King (1966) model fit to observed clusters might help identify the presence of post core-collapse BH subsystems. For global observables like core and half-mass radii the direct simulations agree well with Monte-Carlo models. Variations in the initial mass function can result in significantly different GC properties (e.g. density distributions) driven by varying amounts of early mass loss and the number of forming BHs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

MOCCA code for star cluster simulations - III. Stellar-mass black holes in the globular cluster M22

Using a Monte Carlo code, we construct a dynamic evolutionary model of the Galactic globular cluster M22 (NGC6656). The initial conditions are chosen so that, after about 12Gyr of stellar and dynamical evolution, the model is an approximate fit to the surface brightness and velocity dispersion profiles of the cluster, to its mass function, and to the current binary fraction. Depending on the distribution of black hole natal kicks, we predict that the present-day population of stellar-mass black holes ranges from about 40 (no kicks) down to essentially zero (kicks distributed like those of neutron stars). Provided that natal kicks do not eject all new black holes, it is suggested that clusters with a present-day half-mass relaxation time above about 1Gyr are the ones that may still retain an appreciable population of black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Revisiting the universality of (multiple) star formation in present-day star formation regions

Populations of multiple stars inside clustered regions are known to change through dynamical interactions. The efficiency of binary disruption is thought to be determined by stellar density. King and collaborators recently investigated the multiplicity properties in young star forming regions and in the Galactic field. They concluded that stellar density-dependent modification of a universal initial binary population (the standard or null hypothesis model) cannot explain the observations. We re-visit their results, analyzing the data within the framework of different model assumptions, namely non-universality without dynamical modification and universality with dynamics. We illustrate that the standard model does account for all known populations if regions were significantly denser in the past. Some of the effects of using present-day cluster properties as proxies for their past values are emphasized and that the degeneracy between age and density of a star forming region can not be omitted when interpreting multiplicity data. A new analysis of the Corona Australis region is performed within the standard model. It is found that this region is likely as unevolved as Taurus and an initial density of $\approx190M_{\odot}\;pc^{-3}$ is required to produce the presently observed binary population, which is close to its present-day density.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The state of globular clusters at birth II: primordial binaries

(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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Monte Carlo Simulations of Star Clusters - VII. The globular cluster 47 Tuc

We describe Monte Carlo models for the dynamical evolution of the massive globular cluster 47 Tuc (NGC 104). The code includes treatments of two-body relaxation, most kinds of three- and four-body interactions involving primordial binaries and those formed dynamically, the Galactic tide, and the internal evolution of both single and binary stars. We arrive at a set of initial parameters for the cluster which, after 12Gyr of evolution, gives a model with a fairly satisfactory match to surface brightness and density profiles, the velocity dispersion profile, the luminosity function in two fields, and the acceleration of pulsars. Our models appear to require a relatively steep initial mass function for stars above about turnoff, with an index of about 2.8 (where the Salpeter mass function has an index of 2.35), and a relatively flat initial mass function (index about 0.4) for the lower main sequence. According to the model, the current mass is estimated at 0.9 million solar masses, of which about 34% consists of remnants. We find that primordial binaries are gradually taking over from mass loss by stellar evolution as the main dynamical driver of the core. Despite the high concentration of the cluster, core collapse will take at least another 20Gyr.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

1Gyr in the Life of the Globular Cluster NGC 6397

M4 and NGC 6397 are two very similar galactic globular clusters, which differ mainly in their surface brightness profile. M4 has a classic King-like profile, whereas NGC 6397 has a more concentrated profile, which is often interpreted as that of a post-core collapse cluster. Heggie & Giersz (2008), however, found that M4 is also a post-core collapse cluster, and Giersz & Heggie (2009) concluded that the main reason for the difference between the two surface brightness profiles is fluctuations. This conclusion was reached on the basis of Monte Carlo models, however, and in the present Letter we verify that similar fluctuations occur in N-body models. The models were initialised by generating initial conditions from the Monte Carlo model of NGC 6397 at the simulated age of 12Gyr, and one was followed for 1Gyr. The new models help to clarify the nature of the fluctuations, which have the nature of semi-regular oscillations with a time scale of order 10^8yr. They are influenced by the dynamical role which is played by primordial binaries in the evolution of the core.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Monte Carlo simulations of star clusters - II. Tidally limited, multi-mass systems with stellar evolution

A revision of Stod\{'o}łkiewicz's Monte Carlo code is used to simulate evolution of large star clusters. A survey of the evolution of N-body systems influenced by the tidal field of a parent galaxy and by stellar evolution is presented. The results presented are in good agreement with theoretical expectations and the results of other methods (Fokker-Planck, Monte Carlo and N-body). The initial rapid mass loss, due to stellar evolution of the most massive stars, causes expansion of the whole cluster and eventually leads to the disruption of less bound systems ($W_0=3$). Models with larger $W_0$ survive this phase of evolution and then undergo core collapse and subsequent post-collapse expansion, like isolated models. The expansion phase is eventually reversed when tidal limitation becomes important. The results presented are the first major step in the direction of simulating evolution of real globular clusters by means of the Monte Carlo method.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Statistics of N-Body Simulations. III. Unequal Masses

We describe results from large numbers of $N$-body simulations containing from $250$ to $1000$ stars each. The distribution of stellar masses is a power law, and the systems are isolated. While the collapse of the core exhibits the expected segregation of different masses, we find that the post-collapse evolution is, at a first approximation, homologous. This is quite surprising because there is no reason for supposing that mass segregation should not continue to have a substantial effect on the evolution of the cluster. In fact the spatial distribution of the mean stellar mass is nearly static throughout the post-collapse regime, except for the overall expansion of the systems, and this helps to explain why the post-collapse evolution is nearly self-similar. Self-similarity is also exhibited by the distribution of anisotropy and the profile of departures from equipartition, which show little change during the post-collapse phase. The departures from energy equipartition and isotropy are small in the core and increase with radius. During post-collapse evolution massive stars (mainly) are removed from the system by binary activity. This effect dominates the preferential escape of low-mass stars due to standard two-body relaxation processes.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

Statistics of N-Body Simulations. II. Equal Masses after Core Collapse

This paper presents and analyses statistical results from a large number of $N$-body simulations of isolated systems with equal masses, in which $250\le N\le 2000$. It concentrates on the phase starting around the end of core collapse, when binaries start to play a crucial role. Interactions of hard binaries are consistent with the Spitzer (1987) cross section. The evolution of the half-mass radius after core collapse nearly follows classical theory, The evolution of the inner parts of the system around the time of core bounce is consistent with these simplified models provided that the continuous production of energy, as usually assumed, is replaced by a model of {\sl stochastic} energy production. Similarly, post-collapse evolution of the core requires a modest recalibration of the coefficient of energy generation, especially for small $N$. The distributions of velocity dispersion and anisotropy become remarkably homologous soon after core bounce. The bound mass of the systems very nearly follows a power law with time. A small number of escapers, presumed to be those associated with binary activity, dominate the energy which is carried off. The ``internal" energy of escaping binaries is consistent with theoretical expectations.

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

Comparing Direct N-Body Integration with Anisotropic Gaseous Models of Star Clusters

We compare the results for the dynamical evolution of star clusters derived from anisotropic gaseous models with the data from N-body simulations of isolated and one-component systems, each having modest number of stars. The statistical quality of N-body data was improved by averaging results from many N-body runs, each with the same initial parameters but with different sequences of random numbers used to initialize positions and velocities of the particles. We study the development of anisotropy, the spatial evolution and energy generation by three-body binaries and its N-dependence. We estimate the following free parameters of anisotropic gaseous models: the time scale for collisional anisotropy decay and the coefficient in the formulae for energy generation by three-body binaries. To achieve a fair agreement between N-body and gaseous models for the core in pre- as well as in post-collapse only the energy generation by binaries had to be varied by N. We find that anisotropy has considerable influence on the spatial structure of the cluster particularly for the intermediate and outer regions.

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

Statistics of N-body Simulations I. Equal Masess Before Core Collapse

We study the dynamical evolution of idealised stellar systems by averaging results from many $N$-body simulations, each having modest numbers of stars. For isolated systems with stars of uniform mass, we discuss aspects of evolution up to the point of core collapse: relaxation and its $N$-dependence, the evolution of the density profile, the development of the velocity dispersion and anisotropy, and the rate of stellar escape.

0 citations0 reviewsNo code linkedOpen access