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

P. Berczik

P. Berczik 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.COastro-phastro-ph.SR

Catalog footprint

What is connected

12works
4topics
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

12 published item(s)

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.

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Simulation of the loss-cone instability in spherical systems. I. Dominating harmonic potential

A new so-called `gravitational loss-cone instability' in stellar systems has recently been investigated theoretically in the framework of linear perturbation theory and proved to be potentially important in understanding the physical processes in centres of galaxies, star clusters, and the Oort comet cloud. Using N-body simulations of a toy model, we confirm previous findings for the harmonic dominating potential and go beyond the linear theory. Unlike the well-known instabilities, the new one shows no notable change in the spherical geometry of the cluster, but it significantly accelerates the speed of diffusion of particles in phase space leading to an early instability saturation.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Simulation of the loss-cone instability in spherical systems. II. Dominating Keplerian potential

A new so-called `gravitational loss-cone instability' in stellar systems has recently been investigated theoretically in the framework of linear perturbation theory and proved to be potentially important in understanding the physical processes in centres of galaxies, star clusters, and the Oort comet cloud. Using N-body simulations, we confirm previous findings and go beyond the linear theory. Unlike the well-known instabilities, the new one shows no notable change in spherical geometry of the cluster, but it significantly accelerates the speed of diffusion of particles in phase space leading to a repopulation of the loss cone and early instability saturation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Roche volume filling and the dissolution of open star clusters

From direct N-body simulations we find that the dynamical evolution of star clusters is strongly influenced by the Roche volume filling factor. We present a parameter study of the dissolution of open star clusters with different Roche volume filling factors and different particle numbers. We study both Roche volume underfilling and overfilling models and compare with the Roche volume filling case. We find that in the Roche volume overfilling limit of our simulations two-body relaxation is no longer the dominant dissolution mechanism but the changing cluster potential. We call this mechnism "mass-loss driven dissolution" in contrast to "two-body relaxation driven dissolution" which occurs in the Roche volume underfilling regime. We have measured scaling exponents of the dissolution time with the two-body relaxation time. In this experimental study we find a decreasing scaling exponent with increasing Roche volume filling factor. The evolution of the escaper number in the Roche volume overfilling limit can be described by a log-logistic differential equation. We report the finding of a resonance condition which may play a role for the evolution of star clusters and may be calibrated by the main periodic orbit in the large island of retrograde quasiperiodic orbits in the Poincaré surfaces of section. We also report on the existence of a stability curve which may be of relevance with respect to the structure of star clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simulations of the Hyades

Context: Using the recent observational data of Röser et al. we present $N$-body simulations of the Hyades open cluster. Aims: We make an attempt to determine initial conditions of the Hyades cluster at the time of its formation in order to reproduce the present-day cumulative mass profile, stellar mass and luminosity function (LF). Methods: We performed direct $N$-body simulations of the Hyades in an analytic Milky Way potential that account for stellar evolution and include primordial binaries in a few models. Furthermore, we applied a Kroupa (2001) IMF and used extensive ensemble-averaging. Results: We find that evolved single-star King initial models with King parameters $W_0 = 6-9$ and initial particle numbers $N_0 = 3000$ provide good fits to the observational present-day cumulative mass profile within the Jacobi radius. The best-fit King model has an initial mass of $1721\ M_\odot$ and an average mass loss rate of $-2.2 \ M_\odot/\mathrm{Myr}$. The K-band LFs of models and observations show a reasonable agreement. Mass segregation is detected in both observations and models. If 33% primordial binaries are included the initial particle number is reduced by 5% as compared to the model without primordial binaries. Conclusions: The present-day properties of the Hyades can be well reproduced by a standard King or Plummer initial model when choosing appropriate initial conditions. The degeneracy of good-fitting models can be quite high due to the large dimension of the parameter space. More simulations with different Roche-lobe filling factors and primordial binary fractions are required to explore this degeneracy in more detail.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Calibration of radii and masses of open clusters with a simulation

Context: A recent new approach to apply a simple dynamical mass estimate of tidally limited star clusters is based on the identification of the tidal radius in a King profile with the dynamical Jacobi radius. The application to an unbiased open cluster catalogue yields significantly higher cluster masses compared to the classical methods. Aims: We quantify the bias in the mass determination as function of projection direction and cluster age by analysing a simulated star cluster. Methods: We use direct $N$-body simulations of a star cluster including stellar evolution in an analytic Milky Way potential and apply a best fit to the projected number density of cluster stars. Results: We obtain significantly overestimated star cluster masses which depend strongly on the viewing direction. The overestimation is typically in the range of 10-50 percent and reaches a factor of 3.5 for young clusters. Mass segregation reduces the derived limiting radii systematically.

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

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Quantitative analysis of clumps in the tidal tails of star clusters

Tidal tails of star clusters are not homogeneous but show well defined clumps in observations as well as in numerical simulations. Recently an epicyclic theory for the formation of these clumps was presented. A quantitative analysis was still missing. We present a quantitative derivation of the angular momentum and energy distribution of escaping stars from a star cluster in the tidal field of the Milky Way and derive the connection to the position and width of the clumps. For the numerical realization we use star-by-star $N$-body simulations. We find a very good agreement of theory and models. We show that the radial offset of the tidal arms scales with the tidal radius, which is a function of cluster mass and the rotation curve at the cluster orbit. The mean radial offset is 2.77 times the tidal radius in the outer disc. Near the Galactic centre the circumstances are more complicated, but to lowest order the theory still applies. We have also measured the Jacobi energy distribution of bound stars and showed that there is a large fraction of stars (about 35%) above the critical Jacobi energy at all times, which can potentially leave the cluster. This is a hint that the mass loss is dominated by a self-regulating process of increasing Jacobi energy due to the weakening of the potential well of the star cluster, which is induced by the mass loss itself.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

The study of colliding molecular clumps evolution

The results of study of the gravitational fragmentation in the interstellar medium (ISM) by clump-clump collisions are presented. We suggest, that collision of clumps, that are subparts of Giant Molecular Clouds (GMC) may be on of the basic mechanism, which result to ISM fragmentation and define the dynamical as well as statistical characteristics (e.g. the mass spectra) of protostellar condensation. In the present paper, we describe our 3D SPH-modeling, in isothermal approximation, of supersonic collisions of two identical clumps with a few variants of initial impact parameters ($β$), that cover the wide range. Our results shown, that at all $β$ in system began intensive fragmentation. The resulting fragments mass function depend from initial impact parameter. The obtained mass spectra have the slopes in a good enough agreement with observational data for our Galaxy -- especially for large impact parameters, which are more realistic as for large clumps ensembles.

0 citations0 reviewsNo code linkedOpen access