BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Peter Berczik

Peter Berczik contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph.HEastro-ph.SRgr-qc

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

8 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
preprint2023arXiv

Statistical Analysis of the Probability of Interaction of Globular Clusters with Each Other and with the Galactic Center on the Cosmological Time Scale According to Gaia DR2 Data

This study is aimed at investigating the dynamic evolution of the orbits of stellar globular clusters (GCs). To integrate the orbits backward in time, the authors use models of the time-varying potentials derived from cosmological simulations, which are closest to the potential of the Galaxy. This allows for estimating the probability of close passages (collisions) of GCs with respect to each other and the Galactic center (GalC) in the Galaxy undergoing dynamic changes in the past. To reproduce the dynamics of the Galaxy in time, five potentials selected from the IllustrisTNG-100 large-scale cosmological database, which are similar in their characteristics to the current physical parameters of the Milky Way, are used. With these time-varying potentials, we have reproduced the orbital trajectories of 143 GCs 10 Gyr back in time using our original phi-GPU N-body code. Each GC was treated as a single physical particle with the assigned position and velocity of the cluster center from the Gaia DR2 observations. For each of the potentials, 1000 initial conditions were generated with randomized initial velocities of GCs within the errors of the observational data. In this study, we consider close passages to be passages with a relative distance of less than 100 pc and a relative speed of less than 250 km/s. To select clusters that pass at close distances from the GalC, the following criterion is applied based only on the relative distance: it must be less than 100 pc. Applying the above criteria, the authors obtained statistically significant rates of close passages of GCs with respect to each other and to the GalC. It has been determined that GCs during their evolution have approximately 10 intersecting trajectories with each other on the average and approximately 3 to 4 close passages near the GalC in 1 Gyr at a distance of 50 pc for each of the chosen potentials.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Eccentricity evolution of massive black hole binaries from formation to coalescence

Coalescing supermassive black hole binaries (BHBs) are expected to be the loudest sources of gravitational waves (GWs) in the Universe. Detection rates for ground or space-based detectors based on cosmological simulations and semi-analytic models are highly uncertain. A major difficulty stems from the necessity to model the BHB from the scale of the merger to that of inspiral. Of particular relevance to the GW merger timescale is the binary eccentricity. Here we present a self-consistent numerical study of the eccentricity of BHBs formed in massive gas-free mergers from the early stages of the merger to the hardening phase, followed by a semi-analytical model down to coalescence. We find that the early eccentricity of the unbound black hole pair is largely determined by the initial orbit. It systematically decreases during the dynamical friction phase. The eccentricity at binary formation is affected by stochasticity and noise owing to encounters with stars, but preserves a strong correlation with the initial orbital eccentricity. Binding of the black holes is a phase characterised by strong perturbations, and we present a quantitative definition of the time of binary formation. During hardening the eccentricity increases in minor mergers, unless the binary is approximately circular, but remains largely unchanged in major mergers, in agreement with predictions from semi-analytical models based on isotropic scattering experiments. Coalescence times due to hardening and GW emission in gas-poor non-rotating ellipticals are <~0.5 Gyr for the large initial eccentricities (0.5 < e < 0.9) typical of galaxy mergers in cosmological simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Merging of unequal mass binary black holes in non-axisymmetric galactic nuclei

In this work, we study the stellar-dynamical hardening of unequal mass supermassive black hole (SMBH) binaries in the central regions of merging galactic nuclei. We present a comprehensive set of direct $N$-body simulations of the problem, varying both the total mass and the mass ratio of the SMBH binary (SMBHB). Simulations were carried out with the $φ-$GPU $N$-body code, which enabled us to fully exploit supercomputers equipped with graphic processing units (GPUs). As a model for the galactic nuclei, we adopted initial axisymmetric, rotating models, aimed at reproducing the properties of a galactic nucleus emerging from a galaxy merger event, containing two SMBHs which were unbound initially. We found no &#39;final-parsec problem&#39;, as our SMBHs tend to pair and shrink without showing significant signs of stalling. This confirms earlier results and extends them to large particle numbers and rotating systems. We find that the SMBHB hardening depends on the binary-reduced mass ratio via a single parameter function. Our results suggest that, at a fixed value for the SMBHB primary mass, the merger time of highly asymmetric binaries is up to four order of magnitudes smaller than the equal-mass binaries. This can significantly affect the population of SMBHs potentially detectable as gravitational wave sources.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Revisit the rate of tidal disruption events: the role of the partial tidal disruption event

Tidal disruption of stars in dense nuclear star clusters containing supermassive central black holes (SMBH) is modeled by high-accuracy direct N-body simulation. Stars getting too close to the SMBH are tidally disrupted and a tidal disruption event (TDE) happens. TDEs probe properties of SMBH, their accretion disks, and the surrounding nuclear stellar cluster. In this paper we compare rates of full tidal disruption events (FTDE) with partial tidal disruption events (PTDE). Since a PTDE does not destroy the star, a leftover object emerges; we use the term &#39;leftover star&#39; for it; two novel effects occur in the simulation: (1) variation of the leftover star&#39;s mass and radius, (2) variation of the leftover star&#39;s orbital energy. After switching on these two effects in our simulation, the number of FTDEs is reduced by roughly 28%, and the reduction is mostly due to the ejection of the leftover stars from PTDEs coming originally from relatively large distance. The number of PTDEs is about 75% higher than the simple estimation given by Stone et al. (2020), and the enhancement is mainly due to the multiple PTDEs produced by the leftover stars residing in the diffusive regime. We compute the peak mass fallback rate for the PTDEs and FTDEs recorded in the simulation, and find 58% of the PTDEs have peak mass fallback rate exceeding the Eddington limit, and the number of super-Eddington PTDEs is 2.3 times the number of super-Eddington FTDEs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bimodality of [α/Fe]-[Fe/H] distributions is a natural outcome of dissipative collapse and disc growth in Milky Way-type galaxies

We present a set of self-consistent chemo-dynamical simulations of MW-type galaxies formation to study the origin of the bimodality of $α$-elements in stellar populations. We explore how the bimodality is related to the geometrically and kinematically defined stellar discs, gas accretion and radial migration. We find that the two $α$-sequences are formed in quite different physical environments. The high-$α$ sequence is formed early from a burst of star formation (SF) in a turbulent, compact gaseous disc which forms a thick disc. The low-$α$ stellar populations is the result of quiescent SF supported by the slow accretion of enriched gas onto a radially extended thin disc. Stellar feedback-driven outflows during the formation of the thick disc are responsible for the enrichment of the surrounding gaseous halo, which subsequently feeds the disc on a longer time-scale. During the thin disc phase, chemical evolution reaches an equilibrium metallicity and abundance, where the stars pile-up. This equilibrium metallicity decreases towards the outer disc, generating the ridgeline that forms the low-$α$ sequence. We identify a second mechanism capable of creating a low-$α$ sequence in one of our simulations. Rapid shutdown of the SF, provoked by the feedback at the end of the thick disc phase, suppresses the chemical enrichment of the halo gas, which, once accreted onto the star-forming disc, dilutes the ISM at the beginning of the thin disc formation. Both mechanisms can operate in a galaxy, but the former is expected to occur when SF efficiency ceases to be dominated by the formation of the thick disc, while the latter can occur in the inner regions. Being the result of the presence of low and high gas density environments, the bimodality is independent of any particular merger history, suggesting that it could be much more widespread than has been claimed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

BSE versus StarTrack: implementations of new wind, remnant-formation, and natal-kick schemes in NBODY7 and their astrophysical consequences

The masses of stellar-remnant black holes (BH), as a result of their formation via massive single- and binary-stellar evolution, is of high interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star (BNS) mergers. Here we present new developments in the N-body evolution program NBODY7 in regards to its stellar-remnant formation and related schemes. We demonstrate that the newly-implemented stellar-wind and remnant-formation schemes in the NBODY7 code&#39;s BSE sector, such as the &#39;rapid&#39; and the &#39;delayed&#39; supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produces neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main sequence (ZAMS) mass and metallicity, with those from the StarTrack population-synthesis program. We also demonstrate the new implementations of various natal-kick mechanisms on NSs and BHs such as the &#39;convection-asymmetry-driven&#39;, &#39;collapse-asymmetry-driven&#39;, and &#39;neutrino-emission-driven&#39; kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. We find that the SN material fallback causes the convection-asymmetry kick to effectively retain similar number and mass of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry kick would cause nearly all BHs to retain in clusters irrespective of remnant formation model and metallicity, whereas the inference of a large number of BHs in GCs would potentially rule out the neutrino-driven kick mechanism. Pre-SN mergers of massive primordial binaries would cause BH masses to deviate from the single-star ZAMS mass-remnant mass relation. Such mergers, at low metallicities, can produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Star Clusters in the Galactic tidal field, from birth to dissolution

We study the evolution of star clusters in the Galactic tidal field starting from their birth in molecular clumps. Our model clusters form according to the local-density-driven cluster formation model in which the stellar density profile is steeper than that of gas. As a result, clusters resist the gas expulsion better than predicted by earlier models. We vary the impact of the Galactic tidal field λ, considering different Galactocentric distances (3-18 kpc), as well as different cluster sizes. Our model clusters survive the gas expulsion independent of λ. We investigated the relation between the cluster mass at the onset of secular evolution and their dissolution time. The model clusters formed with a high star-formation efficiency (SFE) follow a tight mass-dependent dissolution relation, in agreement with previous theoretical studies. However, the low-SFE models present a shallower mass-dependent relation than high-SFE clusters, and most dissolve before reaching 1 Gyr (cluster teenage mortality).

0 citations0 reviewsNo code linkedOpen access