Engineering topics

Umbreit, Stefan

Publications and source records attributed to Umbreit, Stefan.

3 records

DOE OSTI2015

THE DYNAMICAL EVOLUTION OF STELLAR BLACK HOLES IN GLOBULAR CLUSTERS

Our current understanding of the stellar initial mass function and massive star evolution suggests that young globular clusters (GCs) may have formed hundreds to thousands of stellar-mass black holes (BHs), the remnants of stars with initial masses from ∼20-100 M {sub ☉}. Birth kicks from supernova explosions may eject some BHs from their birth clusters, but most should be retained. Using a Monte Carlo method we investigate the long-term dynamical evolution of GCs containing large numbers of stellar BHs. We describe numerical results for 42 models, covering a broad range of realistic initial conditions, including up to 1.6 × 10{sup 6} stars. In almost all models we find that significant numbers of BHs (up to ∼10{sup 3}) are retained all the way to the present. This is in contrast to previous theoretical expectations that most BHs should be ejected dynamically within a few gigayears The main reason for this difference is that core collapse driven by BHs (through the Spitzer {sup m}ass segregation instability{sup )} is easily reverted through three-body processes, and involves only a small number of the most massive BHs, while lower-mass BHs remain well-mixed with ordinary stars far from the central cusp. Thus the rapid segregation of stellar BHs does not lead to a long-term physical separation of most BHs into a dynamically decoupled inner core, as often assumed previously. Combined with the recent detections of several BH X-ray binary candidates in Galactic GCs, our results suggest that stellar BHs could still be present in large numbers in many GCs today, and that they may play a significant role in shaping the long-term dynamical evolution and the present-day dynamical structure of many clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

CONSTRAINING INTERMEDIATE-MASS BLACK HOLES IN GLOBULAR CLUSTERS

Decades after the first predictions of intermediate-mass black holes (IMBHs) in globular clusters (GCs) there is still no unambiguous observational evidence for their existence. The most promising signatures for IMBHs are found in the cores of GCs, where the evidence now comes from the stellar velocity distribution, the surface density profile, and, for very deep observations, the mass-segregation profile near the cluster center. However, interpretation of the data, and, in particular, constraints on central IMBH masses, require the use of detailed cluster dynamical models. Here we present results from Monte Carlo cluster simulations of GCs that harbor IMBHs. As an example of application, we compare velocity dispersion, surface brightness and mass-segregation profiles with observations of the GC M10, and constrain the mass of a possible central IMBH in this cluster. We find that, although M10 does not seem to possess a cuspy surface density profile, the presence of an IMBH with a mass up to 0.75% of the total cluster mass, corresponding to about 600 M{sub Sun }, cannot be excluded. This is also in agreement with the surface brightness profile, although we find it to be less constraining, as it is dominated by the light of giants, causing it to fluctuate significantly. We also find that the mass-segregation profile cannot be used to discriminate between models with and without IMBH. The reason is that M10 is not yet dynamically evolved enough for the quenching of mass segregation to take effect. Finally, detecting a velocity dispersion cusp in clusters with central densities as low as in M10 is extremely challenging, and has to rely on only 20-40 bright stars. It is only when stars with masses down to 0.3 M{sub Sun} are included that the velocity cusp is sampled close enough to the IMBH for a significant increase above the core velocity dispersion to become detectable.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

A PARALLEL MONTE CARLO CODE FOR SIMULATING COLLISIONAL N-BODY SYSTEMS

We present a new parallel code for computing the dynamical evolution of collisional N-body systems with up to N {approx} 10{sup 7} particles. Our code is based on the Henon Monte Carlo method for solving the Fokker-Planck equation, and makes assumptions of spherical symmetry and dynamical equilibrium. The principal algorithmic developments involve optimizing data structures and the introduction of a parallel random number generation scheme as well as a parallel sorting algorithm required to find nearest neighbors for interactions and to compute the gravitational potential. The new algorithms we introduce along with our choice of decomposition scheme minimize communication costs and ensure optimal distribution of data and workload among the processing units. Our implementation uses the Message Passing Interface library for communication, which makes it portable to many different supercomputing architectures. We validate the code by calculating the evolution of clusters with initial Plummer distribution functions up to core collapse with the number of stars, N, spanning three orders of magnitude from 10{sup 5} to 10{sup 7}. We find that our results are in good agreement with self-similar core-collapse solutions, and the core-collapse times generally agree with expectations from the literature. Also, we observe good total energy conservation, within {approx}< 0.04% throughout all simulations. We analyze the performance of the code, and demonstrate near-linear scaling of the runtime with the number of processors up to 64 processors for N = 10{sup 5}, 128 for N = 10{sup 6} and 256 for N = 10{sup 7}. The runtime reaches saturation with the addition of processors beyond these limits, which is a characteristic of the parallel sorting algorithm. The resulting maximum speedups we achieve are approximately 60 Multiplication-Sign , 100 Multiplication-Sign , and 220 Multiplication-Sign , respectively.

79 ASTRONOMY AND ASTROPHYSICS↗