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

Rainer Spurzem

Rainer Spurzem contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph.SRastro-ph.EPastro-ph.HEastro-ph.IMgr-qc

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
6topics
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

13 published item(s)

preprint2025arXiv

From Nbody1 to Nbody7: the Growth of Sverres Industry

From NBODY1 to NBODY6 : The Growth of an Industry is the title of a 1999 invited review by Sverre Aarseth, for Publications of the Astronomical Society of the Pacific (PASP). I took this as an inspiration for the title of this paper; it describes how Sverres NBODY Industry has further grown since 90s of the previous century, and how it is further flourishing and hopefully developing, in his spirit, even after the sad news of his passing away reached us. My contact and friendship with Sverre started a few decades ago being sent to Cambridge to learn NBODY5, counting input parameters, and learning about the fact that even a sophisticated code (which had already at that time quite a history) requires permanent maintenance and bug fixes. Managed by Sverre, who relentlessly ran his code and responded to the widely spread crowd of customer colleagues. There has been a phase of massive and fast development and improvements due to vectorization, parallelization, GRAPE and GPU acceleration, and Sverre has been always on top of it if not ahead, but also fully adopting ideas of collaborators, once they tested well. NBODY6++GPU and NBODY7 entered the scene, and also recent new competitors, such as PETAR or BIFROST . We all have learnt a lot from Sverre, and strive to continue in his open-minded spirit, for open source and exchange. A striking evidence for the further growth of the industry is the number of papers here (and two of them follow in this session, but also in other sessions) using and further developing the aforementioned codes, as well as the occurrence of new and competing codes, which keep the field alive.

0 citations0 reviewsNo code linkedOpen access
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

Birth cluster simulations of planetary systems with multiple super-Earths: initial conditions for white dwarf pollution drivers

Previous investigations have revealed that eccentric super-Earths represent a class of planets which are particularly effective at transporting minor bodies towards white dwarfs and subsequently polluting their atmospheres with observable chemical signatures. However, the lack of discoveries of these planets beyond a few astronomical units from their host stars prompts a better understanding of their orbital architectures from their nascent birth cluster. Here, we perform stellar cluster simulations of 3-planet and 7-planet systems containing super-Earths on initially circular, coplanar orbits. We adopt the typical stellar masses of main-sequence progenitors of white dwarfs ($1.5\,\mathrm{M}_{\odot}$-$2.5\,\mathrm{M}_{\odot}$) as host stars and include 8,000 main-sequence stars following a Kroupa initial mass function in our clusters. Our results reveal that about 30 per cent of the simulated planets generate eccentricities of at least 0.1 by the time of cluster dissolution, which would aid white dwarf pollution. We provide our output parameters to the community for potential use as initial conditions for subsequent evolution 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 'final-parsec problem', 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 'leftover star' for it; two novel effects occur in the simulation: (1) variation of the leftover star's mass and radius, (2) variation of the leftover star'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
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

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's BSE sector, such as the 'rapid' and the 'delayed' 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 'convection-asymmetry-driven', 'collapse-asymmetry-driven', and 'neutrino-emission-driven' 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
preprint2020arXiv

Linking the formation and fate of exo-Kuiper belts within solar system analogues

Escalating observations of exo-minor planets and their destroyed remnants both passing through the solar system and within white dwarf planetary systems motivate an understanding of the orbital history and fate of exo-Kuiper belts and planetesimal discs. Here we explore how the structure of a 40-1000 au annulus of planetesimals orbiting inside of a solar system analogue that is itself initially embedded within a stellar cluster environment varies as the star evolves through all of its stellar phases. We attempt this computationally challenging link in four parts: (1) by performing stellar cluster simulations lasting 100 Myr, (2) by making assumptions about the subsequent quiescent 11 Gyr main-sequence evolution, (3) by performing simulations throughout the giant branch phases of evolution, and (4) by making assumptions about the belt's evolution during the white dwarf phase. Throughout these stages, we estimate the planetesimals' gravitational responses to analogues of the four solar system giant planets, as well as to collisional grinding, Galactic tides, stellar flybys, and stellar radiation. We find that the imprint of stellar cluster dynamics on the architecture of $\gtrsim 100$ km-sized exo-Kuiper belt planetesimals is retained throughout all phases of stellar evolution unless violent gravitational instabilities are triggered either (1) amongst the giant planets, or (2) due to a close ($\ll 10^3$ au) stellar flyby. In the absence of these instabilities, these minor planets simply double their semimajor axis while retaining their primordial post-cluster eccentricity and inclination distributions, with implications for the free-floating planetesimal population and metal-polluted white dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the survival of resonant and non-resonant planetary systems in star clusters

Despite the discovery of thousands of exoplanets in recent years, the number of known exoplanets in star clusters remains tiny. This may be a consequence of close stellar encounters perturbing the dynamical evolution of planetary systems in these clusters. Here, we present the results from direct $N$-body simulations of multiplanetary systems embedded in star clusters containing $N = 8k, 16k, 32k$, and $64k$ stars. The planetary systems, which consist of the four Solar system giant planets Jupiter, Saturn, Uranus, and Neptune, are initialized in different orbital configurations, to study the effect of the system architecture on the dynamical evolution of the entire planetary system, and on the escape rate of the individual planets. We find that the current orbital parameters of the Solar system giants (with initially circular orbits, as well as with present-day eccentricities) and a slightly more compact configuration, have a high resilience against stellar perturbations. A configuration with initial mean-motion resonances of 3:2, 3:2, and 5:4 between the planets, which is inspired by the Nice model, and for which the two outermost planets are usually ejected within the first $10^5$ yr, is in many cases stabilized due to the removal of the resonances by external stellar perturbation and by the rapid ejection of at least one planet. Assigning all planets the same mass of 1 Jovian mass almost equalizes the survival fractions. Our simulations reproduce the broad diversity amongst observed exoplanet systems. We find not only many very wide and/or eccentric orbits, but also a significant number of (stable) retrograde orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Planetary Systems in a Star Cluster II: intermediate-mass black holes and planetary systems

Most stars form in dense stellar environments. It is speculated that some dense star clusters may host intermediate-mass black holes (IMBHs), which may have formed from runaway collisions between high-mass stars, or from the mergers of less massive black holes. Here, we numerically explore the evolution of populations of planets in star clusters with an IMBH. We study the dynamical evolution of single-planet systems and free-floating planets, over a period of 100~Myr, in star clusters without an IMBH, and in clusters with a central IMBH of mass $100~M_\odot$ or $200~M_\odot$. In the central region ($r\lesssim 0.2$~pc), the IMBH&#39;s tidal influence on planetary systems is typically 10~times stronger than the average neighbour star. For a star cluster with a $200M_\odot$ IMBH, the region in which the IMBH&#39;s influence is stronger within the virial radius ($\sim 1$~pc). The IMBH quenches mass segregation, and the stars in the core tend to move towards intermediate regions. The ejection rate of both stars and planets is higher when an IMBH is present. The rate at which planets are expelled from their host star rate is higher for clusters with higher IMBH masses, for $t<0.5 t_{rh}$, while remains mostly constant while the star cluster fills its Roche lobe, similar to a star cluster without an IMBH. The disruption rate of planetary systems is higher in initially denser clusters, and for wider planetary orbits, but this rate is substantially enhanced by the presence of a central IMBH.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The long-term evolution of main-sequence binaries in DRAGON simulations

We present a comprehensive investigation of main-sequence (MS) binaries in the DRAGON simulations, which are the first one-million particles direct $N$-body simulations of globular clusters. We analyse the orbital parameters of the binary samples in two of the DRAGON simulations, D1-R7- IMF93 and D2-R7-IMF01, focusing on their secular evolution and correlations up to 12 Gyr. These two models have different initial stellar mass functions: Kroupa 1993 (D1-R7-IMF93) and Kroupa 2001 (D2-R7-IMF01); and different initial mass ratio distributions: random paring (D1-R7-IMF93) and a power-law (D1-R7-IMF93). In general, the mass ratio of a population of binaries increases over time due to stellar evolution, which is less significant in D2-R7-IMF01. In D1-R7-IMF93, primordial binaries with mass ratio $q \approx$ 0.2 are most common, and the frequency linearly declines with increasing $q$ at all times. Dynamical binaries of both models have higher eccentricities and larger semi-major axes than primordial binaries. They are preferentially located in the inner part of the star cluster. Secular evolution of binary orbital parameters does not depend on the initial mass-ratio distribution, but is sensitive to the initial binary distribution of the system. At t = 12 Gyr, the binary fraction decreases radially outwards, and mass segregation is present. A color difference of 0.1 mag in $F330W-F814W$ and 0.2 mag in $NUV-y$ between the core and the outskirts of both clusters is seen, which is a reflection of the binary radial distribution and the mass segregation in the cluster. The complete set of data for primordial and dynamical binary systems at all snapshot intervals is made publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Interaction of Stars Hosting Planets with Sgr A* Black hole

We present some preliminary results of our ongoing project about planetary systems around S-stars in the vicinity of Sgr A* black hole. Since S-stars might have migrated in the Galactic Centre (GC) from elsewhere, they probably still keep their planetary systems throughout their voyage. In this work, we study the destiny of their putative planetary systems after close interaction with the central black hole of our galaxy.

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