Source author record

R. Spurzem

R. Spurzem 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

15works
9topics
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

15 published item(s)

preprint2022arXiv

Stability of planetary systems within the S-star cluster: the Solar system analogues

A dynamically relaxed dense cluster comprised of about 40 stars (the so-called S-stars) inhabits the central region of the Galaxy. Their stars revolve around the Sgr A$^*$ massive object. To understand the dynamical evolution of planetary systems in a particular environment like that around Sgr A$^*$, we carry out direct N--body simulations of planetary systems embedded in the S--star cluster. In this work, we investigate the short-term stability of the planets orbiting around S-stars after their close interactions with the central massive black hole of our galaxy. We find that planetary systems go through encounters with the SMBH and the nearby stars. We determine the frequency and the strength of planetary systems' encounters with the nearby stars as these encounters remarkably increase for systems assigned to S-stars closer to the SMBH. The SMBH severely destabilizes the planetary systems, though we noted that the small oscillations in the mutual eccentricity and inclination of the planetary system could be caused by the planet--planet coupling and the near-resonance effect between the two planets. We obtain estimates of the fraction of survivor planets ($\sim 51\%$), and find that planets stripped from their hosting star are generally captured on close orbits around Sgr A$^*$. We notify while gas giants are tidally disrupted, terrestrial planets do not. We estimate that Sgr A$^*$ flares can be due to the tidal disruption events of \textit{starless} giant planets.

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

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.

preprint2015arXiv

Dynamical evolution of massive black holes in galactic-scale N-body simulations - introducing the regularized tree code "rVINE"

We present a hybrid code combining the OpenMP-parallel tree code VINE with an algorithmic chain regularization scheme. The new code, called "rVINE", aims to significantly improve the accuracy of close encounters of massive bodies with supermassive black holes in galaxy-scale numerical simulations. We demonstrate the capabilities of the code by studying two test problems, the sinking of a single massive black hole to the centre of a gas-free galaxy due to dynamical friction and the hardening of a supermassive black hole binary due to close stellar encounters. We show that results obtained with rVINE compare well with NBODY7 for problems with particle numbers that can be simulated with NBODY7. In particular, in both NBODY7 and rVINE we find a clear N-dependence of the binary hardening rate, a low binary eccentricity and moderate eccentricity evolution, as well as the conversion of the galaxy's inner density profile from a cusp to a a core via the ejection of stars at high velocity. The much larger number of particles that can be handled by rVINE will open up exciting opportunities to model stellar dynamics close to SMBHs much more accurately in a realistic galactic context. This will help to remedy the inherent limitations of commonly used tree solvers to follow the correct dynamical evolution of black holes in galaxy scale simulations.

preprint2013arXiv

The dynamical evolution of multi-planet systems in open clusters

The majority of stars form in star clusters and many are thought to have planetary companions. We demonstrate that multi-planet systems are prone to instabilities as a result of frequent stellar encounters in these star clusters much more than single-planet systems. The cumulative effect of close and distant encounters on these planetary systems are investigated using Monte Carlo scattering experiments. We consider two types of planetary configurations orbiting Sun-like stars: (i) five Jupiter-mass planets in the semi-major axis range 1-42 AU orbiting a Solar mass star, with orbits that are initially co-planar, circular, and separated by 10 mutual Hill radii, and (ii) the four gas giants of our Solar system. Planets with short orbital periods are not directly affected by encountering stars. However, secular evolution of perturbed systems may result in the ejection of the innermost planets or in physical collisions of the innermost planets with the host star, up to many thousands of years after a stellar encounter. The simulations of the Solar system-like systems indicate that Saturn, Uranus and Neptune are affected by both direct interactions with encountering stars, as well as planet-planet scattering. Jupiter, on the other hand, is almost only affected by direct encounters with neighbouring stars, as its mass is too large to be substantially perturbed by the other three planets. Our results indicate that stellar encounters can account for the apparent scarcity of exoplanets in star clusters, not only for those on wide-orbit that are directly affected by stellar encounters, but also planets close to the star which can disappear long after a stellar encounter has perturbed the planetary system.

preprint2011arXiv

A highly efficient measure of mass segregation in star clusters

Investigations of mass segregation are of vital interest for the understanding of the formation and dynamical evolution of stellar systems on a wide range of spatial scales. Our method is based on the minimum spanning tree (MST) that serves as a geometry-independent measure of concentration. Compared to previous such approaches we obtain a significant refinement by using the geometrical mean as an intermediate-pass. It allows the detection of mass segregation with much higher confidence and for much lower degrees of mass segregation than other approaches. The method shows in particular very clear signatures even when applied to small subsets of the entire population. We confirm with high significance strong mass segregation of the five most massive stars in the Orion Nebula Cluster (ONC). Our method is the most sensitive general measure of mass segregation so far and provides robust results for both data from simulations and observations. As such it is ideally suited for tracking mass segregation in young star clusters and to investigate the long standing paradigm of primordial mass segregation by comparison of simulations and observations.

preprint2011arXiv

Dynamics in Young Star Clusters: From Planets to Massive Stars

The young star clusters we observe today are the building blocks of a new generation of stars and planets in our Galaxy and beyond. Despite their fundamental role we still lack knowledge about the conditions under which star clusters form and the impact of these often harsh environments on the evolution of their stellar and substellar members. We demonstrate the vital role numerical simulations play to uncover both key issues. Using dynamical models of different star cluster environments we show the variety of effects stellar interactions potentially have. Moreover, our significantly improved measure of mass segregation reveals that it can occur rapidly even for star clusters without substructure. This finding is a critical step to resolve the controversial debate on mass segregation in young star clusters and provides strong constraints on their initial conditions.

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.

preprint2010arXiv

Dynamical evolution of rotating dense stellar systems with embedded black holes

Evolution of self-gravitating rotating dense stellar systems (e.g. globular clusters, galactic nuclei) with embedded black holes is investigated. The interaction between the black hole and stellar component in differential rotating flattened systems is followed. The interplay between velocity diffusion due to relaxation and black hole star accretion is investigated together with cluster rotation using 2D+1 Fokker-Planck numerical methods. The models can reproduce the Bahcall-Wolf solution $f \propto E^{1/4}$ ($n \propto r^{-7/4}$) inside the zone of influence of the black hole. Gravo-gyro and gravothermal instabilities conduce the system to a faster evolution leading to shorter collapse times with respect to the non-rotating systems. Angular momentum transport and star accretion support the development of central rotation in relaxation time scales. We explore system dissolution due to mass-loss in the presence of an external tidal field (e.g. globular clusters in galaxies).

preprint2010arXiv

Dynamical friction of massive objects in galactic centres

Dynamical friction leads to an orbital decay of massive objects like young compact star clusters or Massive Black Holes in central regions of galaxies. The dynamical friction force can be well approximated by Chandrasekhar's standard formula, but recent investigations show, that corrections to the Coulomb logarithm are necessary. With a large set of N-body simulations we show that the improved formula for the Coulomb logarithm fits the orbital decay very well for circular and eccentric orbits. The local scale-length of the background density distribution serves as the maximum impact parameter for a wide range of power-law indices of -1 ... -5. For each type of code the numerical resolution must be compared to the effective minimum impact parameter in order to determine the Coulomb logarithm. We also quantify the correction factors by using self-consistent velocity distribution functions instead of the standard Maxwellian often used. These factors enter directly the decay timescale and cover a range of 0.5 ... 3 for typical orbits. The new Coulomb logarithm combined with self-consistent velocity distribution functions in the Chandrasekhar formula provides a significant improvement of orbital decay times with correction up to one order of magnitude compared to the standard case. We suggest the general use of the improved formula in parameter studies as well as in special applications.

preprint2010arXiv

On the isolated dwarf galaxies: from cuspy to flat dark matter density profiles and metallicity gradients

The chemodynamical evolution of spherical multi-component self-gravitating models for isolated dwarf galaxies is studied. We compare their evolution with and without feedback effects from star formation processes. We find that initially cuspy dark matter profiles flatten with time as a result of star formation, without any special tuning conditions. Thus the seemingly flattened profiles found in many dwarfs do not contradict the cuspy profiles predicted by cosmological models. We also calculate the chemical evolution of stars and gas, to permit comparisons with observational data.

preprint2010arXiv

Orbital evolution of the Carina dwarf galaxy and self-consistent star formation history determination

We present a new study of the evolution of the Carina dwarf galaxy that includes a simultaneous derivation of its orbit and star formation history. The structure of the galaxy is constrained through orbital parameters derived from the observed distance, proper motions, radial velocity and star formation history. The different orbits admitted by the large proper motion errors are investigated in relation to the tidal force exerted by an external potential representing the Milky Way (MW). Our analysis is performed with the aid of fully consistent N-body simulations that are able to follow the dynamics and the stellar evolution of the dwarf system in order to determine self-consistently the star formation history of Carina. We find a star formation history characterized by several bursts, partially matching the observational expectation. We find also compatible results between dynamical projected quantities and the observational constraints. The possibility of a past interaction between Carina and the Magellanic Clouds is also separately considered and deemed unlikely.

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.

preprint1999arXiv

Towards the million body problem on the computer - no news since the three-body problem?

Astrophysical Challenges which demand the solution of the one million (or more) gravitating body problem are briefly discussed for the fields of cosmology, galactic nuclei and globular star clusters. Results from the classical three-body problem are to be combined with methods of statistical mechanics or thermodynamics in order to provide a reasonable model for astrophysically relevant particle numbers of the order 10**4-10**8. For high-precision modelling of star clusters and galactic nuclei over the required hundreds and thousands of crossing times specially tailored ``brute-force'' N-body integration methods are used together with special purpose hardware. A recent implementation of such code on a general purpose parallel computer is presented. After a comparison of the relevance of different N-body algorithms a new concept for a more flexible customized special purpose computer based on a combination of GRAPE and FPGA is proposed. It is an ideal machine for all kinds of N-body simulations using neighbour schemes, as the Ahmad-Cohen direct N-body codes and smoothed particle hydrodynamics (SPH) for systems including interstellar gas dynamics.