Initial Populations of Black Holes in Star Clusters

Using an updated population synthesis code we study the formation and evolution of black holes (BHs) in young star clusters following a massive starburst. This study continues and improves on the initial work described by Belczynski, Sadowski & Rasio (2004). In our new calculations we account for the possible ejections of BHs and their progenitors from clusters because of natal kicks imparted by supernovae and recoil following binary disruptions. The results indicate that the properties of both retained BHs in clusters and ejected BHs (forming a field population) depend sensitively on the depth of the cluster potential. In particular, most BHs ejected from binaries are also ejected from clusters with central escape speeds Vesc < 100 km/s. Conversely, most BHs remaining in binaries are retained by clusters with Vesc > 50 km/s. BHs from single star evolution are also affected significantly: about half of the BHs originating from primordial single stars are ejected from clusters with Vesc < 50 km/s. Our results lay a foundation for theoretical studies of the formation of BH X-ray binaries in and around star clusters, including possible ultra-luminous sources, as well as merging BH--BH binaries detectable with future gravitational-wave observatories.Comment: 35 pages, 8 tables, 17 figures; resubmitted to ApJ (revised version

Intermediate-mass black holes in colliding clusters: Implications for lower-frequency gravitational-wave astronomy

Observations suggest that star clusters often form in binaries or larger bound groups. Therefore, mergers between two clusters are likely to occur. If these clusters both harbor an intermediate-mass black hole (IMBH; 10^{2-4} Msun) in their center, they can become a strong source of gravitational waves when the black holes merge with each other. In order to understand the dynamical processes that operate in such a scenario, one has to study the evolution of the merger of two such young massive star clusters, and more specifically, their respective IMBHs. We employ the direct-summation Nbody4 numerical tool on special-purpose GRAPE6 hardware to simulate a merger of two stellar clusters each containing 63,000 particles and a central IMBH. This allows us to study accurately the orbital evolution of the colliding clusters and the embedded massive black holes. Within ~7 Myr the clusters have merged and the IMBHs constitute a hard binary. The final coalescence happens in ~10^8 yrs. The implication of our analysis is that intermediate-mass black holes merging as the result of coalescence of young dense clusters could provide a source for the Laser Interferometer Space Antenna (LISA) space-based gravitational wave detector mission. We find that interactions with stars increase the eccentricity of the IMBH binary to about 0.8. Although the binary later circularizes by emission of gravitational waves, the residual eccentricity can be detectable through its influence on the phase of the waves if the last few years of inspiral are observed

Relativistic encounters in dense stellar systems

Two coalescing black holes (BHs) represent a conspicuous source of gravitational waves (GWs). The merger involves 17 parameters in the general case of Kerr BHs, so that a successful identification and parameter extraction of the information encoded in the waves will provide us with a detailed description of the physics of BHs. A search based on matched-filtering for characterization and parameter extraction requires the development of some $10^{15}$ waveforms. If a third additional BH perturbed the system, the waveforms would not be applicable, and we would need to increase the number of templates required for a valid detection. In this letter, we calculate the probability that more than two BHs interact in the regime of strong relativity in a dense stellar cluster. We determine the physical properties necessary in a stellar system for three black holes to have a close encounter in this regime and also for an existing binary of two BHs to have a strong interaction with a third hole. In both cases the event rate is negligible. While dense stellar systems such as galactic nuclei, globular clusters and nuclear stellar clusters are the breeding grounds for the sources of gravitational waves that ground-based detectors like Advanced LIGO and Advanced VIRGO will be exploring, the analysis of the waveforms in full general relativity needs only to evaluate the two-body problem. This reduces the number of templates of waveforms to create by orders of magnitude

Physical Processes in Star-Gas Systems

First we present a recently developed 3D chemodynamical code for galaxy evolution from the K**2 collaboration. It follows the evolution of all components of a galaxy such as dark matter, stars, molecular clouds and diffuse interstellar matter (ISM). Dark matter and stars are treated as collisionless N-body systems. The ISM is numerically described by a smoothed particle hydrodynamics (SPH) approach for the diffuse (hot) gas and a sticky particle scheme for the (cool) molecular clouds. Physical processs such as star formation, stellar death or condensation and evaporation processes of clouds interacting with the ISM are described locally. An example application of the model to a star forming dwarf galaxy will be shown for comparison with other codes. Secondly we will discuss new kinds of exotic chemodynamical processes, as they occur in dense gas-star systems in galactic nuclei, such as non-standard ``drag''-force interactions, destructive and gas producing stellar collisions. Their implementation in 1D dynamical models of galactic nuclei is presented. Future prospects to generalize these to 3D are work in progress and will be discussed.Comment: 4 pages, 4 figures, "The 5th Workshop on Galactic Chemodynamics" - Swinburne University (9-11 July 2003). To be published in the Publications of the Astronomical Society of Australia in 2004 (B.K. Gibson and D. Kawata, eds.). Accepted version, minor changes relative to origina

Binary Encounters With Supermassive Black Holes: Zero-Eccentricity LISA Events

Current simulations of the rate at which stellar-mass compact objects merge with supermassive black holes (called extreme mass ratio inspirals, or EMRIs) focus on two-body capture by emission of gravitational radiation. The gravitational wave signal of such events will likely involve a significant eccentricity in the sensitivity range of the Laser Interferometer Space Antenna (LISA). We show that tidal separation of stellar-mass compact object binaries by supermassive black holes will instead produce events whose eccentricity is nearly zero in the LISA band. Compared to two-body capture events, tidal separations have a high cross section and result in orbits that have a large pericenter and small apocenter. Therefore, the rate of interactions per binary is high and the resulting systems are very unlikely to be perturbed by other stars into nearly radial plunges. Depending on the fraction of compact objects that are in binaries within a few parsecs of the center, the rate of low-eccentricity LISA events could be comparable to or larger than the rate of high-eccentricity events.Comment: Final accepted version: ApJ Letters 2005, 631, L11

Anoxic nitrification in marine sediments

Nitrate peaks are found in pore-water profiles in marine sediments at depths considerably below the conventional zone of oxic nitrification. These have been interpreted to represent nonsteady- state effects produced by the activity of nitrifying bacteria, and suggest that nitrification occurs throughout the anoxic sediment region. In this study, ΣNO3 peaks and molecular analysis of DNA and RNA extracted from anoxic sediments of Loch Duich, an organic-rich marine fjord, are consistent with nitrification occurring in the anoxic zone. Analysis of ammonia oxidiser 16S rRNA gene fragments amplified from sediment DNA indicated the abundance of autotrophic ammonia-oxidising bacteria throughout the sediment depth sampled (40 cm), while RT-PCR analysis indicated their potential activity throughout this region. A large non-steady-state pore-water ΣNO3 peak at ~21 cm correlated with discontinuities in this ammonia-oxidiser community. In addition, a subsurface nitrate peak at ~8 cm below the oxygen penetration depth, correlated with the depth of a peak in nitrification rate, assessed by transformation of 15N-labelled ammonia. The source of the oxidant required to support nitrification within the anoxic region is uncertain. It is suggested that rapid recycling of N is occurring, based on a coupled reaction involving Mn oxides (or possibly highly labile Fe oxides) buried during small-scale slumping events. However, to fully investigate this coupling, advances in the capability of high-resolution pore-water techniques are required

On Convergence of the Inexact Rayleigh Quotient Iteration with the Lanczos Method Used for Solving Linear Systems

For the Hermitian inexact Rayleigh quotient iteration (RQI), the author has established new local general convergence results, independent of iterative solvers for inner linear systems. The theory shows that the method locally converges quadratically under a new condition, called the uniform positiveness condition. In this paper we first consider the local convergence of the inexact RQI with the unpreconditioned Lanczos method for the linear systems. Some attractive properties are derived for the residuals, whose norms are $\xi_{k+1}$'s, of the linear systems obtained by the Lanczos method. Based on them and the new general convergence results, we make a refined analysis and establish new local convergence results. It is proved that the inexact RQI with Lanczos converges quadratically provided that $\xi_{k+1}\leq\xi$ with a constant $\xi\geq 1$. The method is guaranteed to converge linearly provided that $\xi_{k+1}$ is bounded by a small multiple of the reciprocal of the residual norm $\|r_k\|$ of the current approximate eigenpair. The results are fundamentally different from the existing convergence results that always require $\xi_{k+1}<1$, and they have a strong impact on effective implementations of the method. We extend the new theory to the inexact RQI with a tuned preconditioned Lanczos for the linear systems. Based on the new theory, we can design practical criteria to control $\xi_{k+1}$ to achieve quadratic convergence and implement the method more effectively than ever before. Numerical experiments confirm our theory.Comment: 20 pages, 8 figures. arXiv admin note: text overlap with arXiv:0906.223

Modest-2: A Summary

This is a summary paper of MODEST-2, a workshop held at the Astronomical Institute ``Anton Pannekoek'' in Amsterdam, 16-17 December 2002. MODEST is a loose collaboration of people interested in MOdelling DEnse STellar systems, particularly those interested in modelling these systems using all the available physics (stellar dynamics, stellar evolution, hydrodynamics and the interplay between the three) by defining interfaces between different codes. In this paper, we summarize 1) the main advances in this endeavour since MODEST-1; 2) the main science goals which can be and should be addressed by these types of simulations; and 3) the most pressing theoretical and modelling advances that we identified.Comment: Accepted by New Astronom