Source author record

Milos Milosavljevic

Milos Milosavljevic 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

37works
10topics
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

37 published item(s)

preprint2022arXiv

Collapsible Linear Blocks for Super-Efficient Super Resolution

With the advent of smart devices that support 4K and 8K resolution, Single Image Super Resolution (SISR) has become an important computer vision problem. However, most super resolution deep networks are computationally very expensive. In this paper, we propose Super-Efficient Super Resolution (SESR) networks that establish a new state-of-the-art for efficient super resolution. Our approach is based on linear overparameterization of CNNs and creates an efficient model architecture for SISR. With theoretical analysis, we uncover the limitations of existing overparameterization methods and show how the proposed method alleviates them. Detailed experiments across six benchmark datasets demonstrate that SESR achieves similar or better image quality than state-of-the-art models while requiring 2x to 330x fewer Multiply-Accumulate (MAC) operations. As a result, SESR can be used on constrained hardware to perform x2 (1080p to 4K) and x4 (1080p to 8K) SISR. Towards this, we estimate hardware performance numbers for a commercial Arm mobile-Neural Processing Unit (NPU) for 1080p to 4K (x2) and 1080p to 8K (x4) SISR. Our results highlight the challenges faced by super resolution on AI accelerators and demonstrate that SESR is significantly faster (e.g., 6x-8x higher FPS) than existing models on mobile-NPU. Finally, SESR outperforms prior models by 1.5x-2x in latency on Arm CPU and GPU when deployed on a real mobile device. The code for this work is available at https://github.com/ARM-software/sesr.

preprint2019arXiv

Constraining Cosmology with Big Data Statistics of Cosmological Graphs

By utilizing large-scale graph analytic tools implemented in the modern Big Data platform, Apache Spark, we investigate the topological structure of gravitational clustering in five different universes produced by cosmological $N$-body simulations with varying parameters: (1) a WMAP 5-year compatible $Λ$CDM cosmology, (2) two different dark energy equation of state variants, and (3) two different cosmic matter density variants. For the Big Data calculations, we use a custom build of stand-alone Spark/Hadoop cluster at Korea Institute for Advanced Study (KIAS) and Dataproc Compute Engine in Google Cloud Platform (GCP) with the sample size ranging from 7 millions to 200 millions. We find that among the many possible graph-topological measures, three simple ones: (1) the average of number of neighbors (the so-called average vertex degree) $α$, (2) closed-to-connected triple fraction (the so-called transitivity) $τ_Δ$, and (3) the cumulative number density $n_{s\ge5}$ of subcomponents with connected component size $s \ge 5$, can effectively discriminate among the five model universes. Since these graph-topological measures are in direct relation with the usual $n$-points correlation functions of the cosmic density field, graph-topological statistics powered by Big Data computational infrastructure opens a new, intuitive, and computationally efficient window into the dark Universe.

preprint2016arXiv

Toward $ab\,initio$ extremely metal poor stars

Extremely metal poor stars have been the focus of much recent attention owing to the expectation that their chemical abundances can shed light on the metal and dust yields of the earliest supernovae. We present our most realistic simulation to date of the astrophysical pathway to the first metal enriched stars. We simulate the radiative and supernova hydrodynamic feedback of a $60\,M_\odot$ Population III star starting from cosmological initial conditions realizing Gaussian density fluctuations. We follow the gravitational hydrodynamics of the supernova remnant at high spatial resolution through its freely-expanding, adiabatic, and radiative phases, until gas, now metal-enriched, has resumed runaway gravitational collapse. Our findings are surprising: while the Population III progenitor exploded with a low energy of $10^{51}\,\text{erg}$ and injected an ample metal mass of $6\,M_\odot$, the first cloud to collapse after the supernova explosion is a dense surviving primordial cloud on which the supernova blastwave deposited metals only superficially, in a thin, unresolved layer. The first metal-enriched stars can form at a very low metallicity, of only $2-5\times10^{-4}\,Z_\odot$, and can inherit the parent cloud's highly elliptical, radially extended orbit in the dark matter gravitational potential.

preprint2015arXiv

Abundance anomalies in metal-poor stars from Population III supernova ejecta hydrodynamics

We present a simulation of the long-term evolution of a Population III supernova remnant in a cosmological minihalo. Employing passive Lagrangian tracer particles, we investigate how chemical stratification and anisotropy in the explosion can affect the abundances of the first low-mass, metal-enriched stars. We find that reverse shock heating can leave the inner mass shells at entropies too high to cool, leading to carbon-enhancement in the re-collapsing gas. This hydrodynamic selection effect could explain the observed incidence of carbon-enhanced metal-poor (CEMP) stars at low metallicity. We further explore how anisotropic ejecta distributions, recently seen in direct numerical simulations of core-collapse explosions, may translate to abundances in metal-poor stars. We find that some of the observed scatter in the Population II abundance ratios can be explained by an incomplete mixing of supernova ejecta, even in the case of only one contributing enrichment event. We demonstrate that the customary hypothesis of fully-mixed ejecta clearly fails if post-explosion hydrodynamics prefers the recycling of some nucleosynthetic products over others. Furthermore, to fully exploit the stellar-archaeological program of constraining the Pop III initial mass function from the observed Pop II abundances, considering these hydrodynamical transport effects is crucial. We discuss applications to the rich chemical structure of ultra-faint dwarf satellite galaxies, to be probed in unprecedented detail with upcoming spectroscopic surveys.

preprint2015arXiv

An Increasing Stellar Baryon Fraction in Bright Galaxies at High Redshift

Recent observations have shown that the characteristic luminosity of the rest-frame ultraviolet (UV) luminosity function does not significantly evolve at 4 < z < 7 and is approximately M*_UV ~ -21. We investigate this apparent non-evolution by examining a sample of 178 bright, M_UV < -21 galaxies at z=4 to 7, analyzing their stellar populations and host halo masses. Including deep Spitzer/IRAC imaging to constrain the rest-frame optical light, we find that M*_UV galaxies at z=4-7 have similar stellar masses of log(M/Msol)=9.6-9.9 and are thus relatively massive for these high redshifts. However, bright galaxies at z=4-7 are less massive and have younger inferred ages than similarly bright galaxies at z=2-3, even though the two populations have similar star formation rates and levels of dust attenuation. We match the abundances of these bright z=4-7 galaxies to halo mass functions from the Bolshoi Lambda-CDM simulation to estimate the halo masses. We find that the typical halo masses in ~M*_UV galaxies decrease from log(M_h/Msol)=11.9 at z=4 to log(M_h/Msol)=11.4 at z=7. Thus, although we are studying galaxies at a similar mass across multiple redshifts, these galaxies live in lower mass halos at higher redshift. The stellar baryon fraction in units of the cosmic mean Omega_b/Omega_m rises from 5.1% at z=4 to 11.7% at z=7; this evolution is significant at the ~3-sigma level. This rise does not agree with simple expectations of how galaxies grow, and implies that some effect, perhaps a diminishing efficiency of feedback, is allowing a higher fraction of available baryons to be converted into stars at high redshifts.

preprint2015arXiv

Confined Population III Enrichment and the Prospects for Prompt Second-Generation Star Formation

It is widely recognized that nucleosynthetic output of the first, Population III supernovae was a catalyst defining the character of subsequent stellar generations. Most of the work on the earliest enrichment was carried out assuming that the first stars were extremely massive and that the associated supernovae were unusually energetic, enough to completely unbind the baryons in the host cosmic minihalo and disperse the synthesized metals into the intergalactic medium. Recent work, however, suggests that the first stars may in fact have been somewhat less massive, with a characteristic mass scale of a few tens of solar masses. We present a cosmological simulation following the transport of the metals synthesized in a Population III supernova assuming that it had an energy of 1e51 ergs, compatible with standard Type II supernovae. A young supernova remnant is inserted in the first star's relic HII region in the free expansion phase and is followed for 40 Myr employing adaptive mesh refinement and Lagrangian tracer particle techniques. The supernova remnant remains partially trapped within the minihalo and the thin snowplow shell develops pronounced instability and fingering. Roughly half of the ejecta turn around and fall back toward the center of the halo, with 1% of the ejecta reaching the center in 30 kyr and 10% in 10 Myr. The average metallicity of the combined returning ejecta and the pristine filaments feeding into the halo center from the cosmic web is 0.001 - 0.01 Z_sun, but the two remain unmixed until accreting onto the central hydrostatic core that is unresolved at the end of the simulation. We conclude that if Population III stars had less extreme masses, they promptly enriched the host minihalos with metals and triggered Population II star formation.

preprint2015arXiv

Metal Transport and Chemical Heterogeneity in Early Star Forming Systems

To constrain the properties of the first stars with the chemical abundance patterns observed in metal-poor stars, one must identify any non-trivial effects that the hydrodynamics of metal dispersal can imprint on the abundances. We use realistic cosmological hydrodynamic simulations to quantify the distribution of metals resulting from one Population III supernova and from a small number of such supernovae exploding in close succession. Overall, supernova ejecta are highly inhomogeneously dispersed throughout the simulations. When the supernova bubbles collapse, quasi-virialized metal-enriched clouds, fed by fallback from the bubbles and by streaming of metal-free gas from the cosmic web, grow in the centers of the dark matter halos. Partial turbulent homogenization on scales resolved in the simulation is observed only in the densest clouds where the vortical time scales are short enough to ensure true homogenization on subgrid scales. However, the abundances in the clouds differ from the gross yields of the supernovae. Continuing the simulations until the cloud have gone into gravitational collapse, we predict that the abundances in second-generation stars will be deficient in the innermost mass shells of the supernova (if only one has exploded) or in the ejecta of the latest supernovae (when multiple have exploded). This indicates that hydrodynamics gives rise to biases complicating the identification of nucleosynthetic sources in the chemical abundance spaces of the surviving stars.

preprint2015arXiv

Radiation pressure driving of a dusty atmosphere

Radiation pressure can be dynamically important in star-forming environments such as ultra-luminous infrared and submillimeter galaxies. Whether and how radiation drives turbulence and bulk outflows in star formation sites is still unclear. The uncertainty in part reflects the limitations of direct numerical schemes that are currently used to simulate radiation transfer and radiation-gas coupling. An idealized setup in which radiation is introduced at the base of a dusty atmosphere in a gravitational field has recently become the standard test for radiation-hydrodynamics methods in the context of star formation. To a series of treatments featuring the flux-limited-diffusion approximation as well as a short-characteristics tracing and M1 closure for the variable Eddington tensor approximation, we here add another treatment that is based on the Implicit Monte Carlo radiation transfer scheme. Consistent with all previous treatments, the atmosphere undergoes Rayleigh-Taylor instability and readjusts to a near-Eddington-limited state. We detect late-time net acceleration in which the turbulent velocity dispersion matches that reported previously with the short-characteristics-based radiation transport closure, the most accurate of the three preceding treatments. Our technical result demonstrates the importance of accurate radiation transfer in simulations of radiative feedback.

preprint2015arXiv

Star Formation for Predictive Primordial Galaxy Formation

The elegance of inflationary cosmology and cosmological perturbation theory ends with the formation of the first stars and galaxies, the initial sources of light that launched the phenomenologically rich process of cosmic reionization. Here we review the current understanding of early star formation, emphasizing unsolved problems and technical challenges. We begin with the first generation of stars to form after the Big Bang and trace how they influenced subsequent star formation. The onset of chemical enrichment coincided with a sharp increase in the overall physical complexity of star forming systems. Ab-initio computational treatments are just now entering the domain of the predictive and are establishing contact with local observations of the relics of this ancient epoch.

preprint2015arXiv

Star formation in the first galaxies - III. Formation, evolution, and characteristics of the first stellar cluster

We simulate the formation of a low metallicity (0.01 Zsun) stellar cluster in a dwarf galaxy at redshift z~14. Beginning with cosmological initial conditions, the simulation utilizes adaptive mesh refinement and sink particles to follow the collapse and evolution of gas past the opacity limit for fragmentation, thus resolving the formation of individual protostellar cores. A time- and location-dependent protostellar radiation field, which heats the gas by absorption on dust, is computed by integration of protostellar evolutionary tracks with the MESA code. The simulation also includes a robust non-equilibrium chemical network that self-consistently treats gas thermodynamics and dust-gas coupling. The system is evolved for 18 kyr after the first protostellar source has formed. In this time span, 30 sink particles representing protostellar cores form with a total mass of 81 Msun. Their masses range from ~0.1 Msun to 14.4 Msun with a median mass ~0.5-1 Msun. Massive protostars grow by competitive accretion while lower-mass protostars are stunted in growth by close encounters and many-body ejections. In the regime explored here, the characteristic mass scale is determined by the temperature floor set by the cosmic microwave background and by the onset of efficient dust-gas coupling. It seems unlikely that host galaxies of the first bursts of metal-enriched star formation will be detectable with the James Webb Space Telescope or other next-generation infrared observatories. Instead, the most promising access route to the dawn of cosmic star formation may lie in the scrutiny of metal-poor, ancient stellar populations in the Galactic neighborhood. The observable targets that correspond to the system simulated here are ultra-faint dwarf satellite galaxies such as Bootes II, Segue I and II, and Willman I.

preprint2015arXiv

The first galaxies: simulating their feedback-regulated assembly

We investigate the formation of a galaxy reaching a virial mass of $~ 10^8$ solar mass at $z=10$ by carrying out a zoomed radiation-hydrodynamical cosmological simulation. This simulation traces Population~III (Pop~III) star formation, characterized by a modestly top-heavy initial mass function (IMF), and considers stellar feedback such as photoionization heating from Pop III and Population~II (Pop~II) stars, mechanical and chemical feedback from supernovae (SNe), and X-ray feedback from accreting black holes (BHs) and high-mass X-ray binaries (HMXBs). We self-consistently impose a transition in star formation mode from top-heavy Pop III to low-mass Pop~II, and find that the star formation rate in the computational box is dominated by Pop~III until $z=13$, and by Pop~II thereafter. The simulated galaxy experiences bursty star formation, with a substantially reduced gas content due to photoionization heating from Pop~III and Pop~II stars, together with SN feedback. All the gas within the simulated galaxy is metal-enriched above $10^{-5}$ solar, such that there are no remaining pockets of primordial gas. The simulated galaxy has an estimated observed flux of $~10^{-3} nJy$, which is too low to be detected by the James Webb Space Telescope (JWST) without strong lensing amplification. We also show that our simulated galaxy is similar in terms of stellar mass to Segue 2, the least luminous dwarf known in the Local Group.

preprint2014arXiv

Dwarf Spheroidal Satellite Formation in a Reionized Local Group

Dwarf spheroidal satellite galaxies have emerged a powerful probe of small-scale dark matter clustering and of cosmic reionization. They exhibit structural and chemical continuity with dwarf irregular galaxies in the field and with spheroidal galaxies in high-density environments. By combining empirical constraints derived for star formation at low gas column densities and metallicities in the local universe with a model for dark matter and baryonic mass assembly, we provide an analytical description of how the dwarf spheroidals acquired their stellar content. Their progenitors formed stars until the gas content, initially reduced from the cosmic average by the thermal pressure of the reionized intergalactic medium, was finally ram pressure stripped during the progenitors' accretion on to the host galaxy. Dwarf spheroidal satellites of differing luminosities seem to share very similar most massive progenitor histories that reach thresholds for gas cooling by atomic line emission at epochs at which the Lagrangian volume of the Local Group should have been reionized. We hypothesize that dwarf spheroidals formed the bulk of their stars in partially rotationally supported HI disks in a reionized universe. This model provides an explanation for the "common mass scale" relation and reproduces the empirical luminosity-size and luminosity-metallicity relations. Explosive feedback phenomena, such as outflows driven by the concerted action of supernovae, need not have been significant in the dwarf spheroidals' formation. We further speculate that the true pre-reionization fossils should exhibit a structure distinct from that of the dwarf spheroidals, e.g., in the form of dense isolated or nuclear star clusters.

preprint2014arXiv

Formation of the First Low-Mass Stars from Cosmological Initial Conditions

We simulate the formation of a metal-poor (10^-2 Zsun) stellar cluster in one of the first galaxies to form in the early Universe, specifically a high-redshift atomic cooling halo (z~14). This is the first calculation that resolves the formation of individual metal-enriched stars in simulations starting from realistic cosmological initial conditions. We follow the evolution of a single dense clump among several in the parent halo. The clump forms a cluster of ~40 stars and sub-stellar objects within 7000 years and could continue forming stars ~5 times longer. Protostellar dust heating has a negligible effect on the star formation efficiency, at least during the early evolutionary stages, but it moderately suppresses gaseous fragmentation and brown dwarf formation. We observe fragmentation in thin gaseous filaments and sustained accretion in larger, rotating structures as well as ejections by binary interactions. The stellar initial mass function above 0.1 Msun, evaluated after ~10^4 years of fragmentation and accretion, seems in agreement with the recent measurement in ultra-faint dwarf spheroidal Galactic satellites of Geha et al. (2013).

preprint2014arXiv

Radiative Feedback from high mass X-ray binaries on the formation of the first galaxies and early reionization

Recent work suggests that the first generation of stars, the so-called Population III (Pop III), could have formed primarily in binaries or as members of small multiple systems. Here we investigate the impact of X-ray feedback from High-Mass X-ray Binaries (HMXBs) left behind in stellar binary systems after the primary forms a black hole (BH), accreting gas at a high rate from the companion, a process that is thought to be favored at the low metallicities characteristic of high-redshift gas. Thanks to their large mean free path, X-rays are capable of preionizing and preheating the gas in the intergalactic medium (IGM) and in haloes long before the reionization of the Universe is complete, and thus could have strongly affected the formation of subsequent generations of stars as well as reionization. We have carried out zoomed hydrodynamical cosmological simulations of minihaloes, accounting for the formation of Pop III stars and their collapse into BHs and HMXBs, and the associated radiation-hydrodynamic feedback from UV and X-ray photons. We find no strong net feedback from HMXBs on the simulated star formation history. On the other hand, the preheating of the IGM by HMXBs leads to a strong suppression of small-scale structures and significantly lowers the recombination rate in the IGM, thus yielding a net positive feedback on reionization. We further show that X-ray feedback from HMXBs can augment the ionizing feedback from the Pop III progenitor stars to suppress gas accretion onto the first BHs, limiting their growth into supermassive BHs. Finally, we show that X-ray ionization by HMXBs leaves distinct signatures in the properties of the high-redshift hydrogen that may be probed in upcoming observations of the redshifted 21cm spin-flip line.

preprint2014arXiv

Recovery from population III supernova explosions and the onset of second generation star formation

We use cosmological simulations to assess how the explosion of the first stars in supernovae (SNe) influences early cosmic history. Specifically, we investigate the impact by SNe on the host systems for Population~III (Pop~III) star formation and explore its dependence on halo environment and Pop~III progenitor mass. We then trace the evolution of the enriched gas until conditions are met to trigger second-generation star formation. To this extent, we quantify the recovery timescale, which measures the time delay between a Pop~III SN explosion and the appearance of cold, dense gas, out of which second-generation stars can form. We find that this timescale is highly sensitive to the Pop~III progenitor mass, and less so to the halo environment. For more massive progenitors, including those exploding in pair instability SNe, second-generation star formation is delayed significantly, for up to a Hubble time. The dependence of the recovery time on the mass of the SN progenitor is mainly due to the ionizing impact of the progenitor star. Photoionization heating increases the gas pressure and initiates a hydrodynamical response that reduces the central gas density, an effect that is stronger in more massive. The gas around lower mass Pop~III stars remains denser and hence the SN remnants cool more rapidly, facilitating the subsequent re-condensation of the gas and formation of a second generation of stars. In most cases, the second-generation stars are already metal-enriched to ~2-5 X 10^{-4}\zsun, thus belonging to Population~II. The recovery timescale is a key quantity governing the nature of the first galaxies, able to host low-mass, long-lived stellar systems. These in turn are the target of future deep-field campaigns with the James Webb Space Telescope.

preprint2013arXiv

A Flux Rope Network and Particle Acceleration in Three Dimensional Relativistic Magnetic Reconnection

We investigate guide-field magnetic reconnection and particle acceleration in relativistic pair plasmas with three-dimensional particle-in-cell (PIC) simulations of a kinetic-scale current sheet in a periodic geometry at low magnetizations. The tearing instability is the dominant mode in the current sheet for all guide field strengths, while the linear kink mode is less important even without guide field. Oblique modes seem to be suppressed entirely. In its nonlinear evolution, the reconnection layer develops a network of interconnected and interacting magnetic flux ropes. As smaller flux ropes merge into larger ones, the reconnection layer evolves toward a three-dimensional, disordered state in which the resulting flux rope segments contain magnetic substructure on plasma skin depth scales. Embedded in the flux ropes, we detect spatially and temporally intermittent sites of dissipation reflected in peaks in the parallel electric field. Magnetic dissipation and particle acceleration persist until the end of the simulations, with simulations with higher magnetization and lower guide field strength exhibiting greater and faster energy conversion and particle energization. At the end of our largest simulation, the particle energy spectrum attains a tail extending to high Lorentz factors that is best modeled with a combination of two additional thermal components. We confirm that the primary energization mechanism is acceleration by the electric field in the X-line region. We discuss the implications of our results for macroscopic reconnection sites, and which of our results may be expected to hold in systems with higher magnetizations.

preprint2013arXiv

Star Formation in the First Galaxies - II: Clustered Star Formation and the Influence of Metal Line Cooling

Population III stars are believed to have been more massive than typical stars today and to have formed in relative isolation. The thermodynamic impact of metals is expected to induce a transition leading to clustered, low-mass Population II star formation. In this work, we present results from three cosmological simulations, only differing in gas metallicity, that focus on the impact of metal fine-structure line cooling on the formation of stellar clusters in a high-redshift atomic cooling halo. Introduction of sink particles allows us to follow the process of gas hydrodynamics and accretion onto cluster stars for 4 Myr corresponding to multiple local free-fall times. At metallicities at least $10^{-3}\, Z_{\odot}$, gas is able to reach the CMB temperature floor and fragment pervasively resulting in a stellar cluster of size $\sim1$ pc and total mass $\sim1000\, M_{\odot}$. The masses of individual sink particles vary, but are typically $\sim100\, M_{\odot}$, consistent with the Jeans mass when gas cools to the CMB temperature, though some solar mass fragments are also produced. At the low metallicity of $10^{-4}\, Z_{\odot}$, fragmentation is completely suppressed on scales greater than 0.01 pc and total stellar mass is lower by a factor of 3 than in the higher metallicity simulations. The sink particle accretion rates, and thus their masses, are determined by the mass of the gravitationally unstable gas cloud and the prolonged gas accretion over many Myr. The simulations thus exhibit features of both monolithic collapse and competitive accretion. Even considering possible dust induced fragmentation that would occur at higher densities, the formation of a bona fide stellar cluster seems to require metal line cooling and metallicities of at least $10^{-3}\, Z_{\odot}$.

preprint2013arXiv

The First Galaxies: Assembly under Radiative Feedback from the First Stars

We investigate how radiative feedback from the first stars affects the assembly of the first dwarf galaxies. We perform cosmological zoomed SPH simulations of a dwarf galaxy assembling inside a halo of virial mass 10^9 solar at z = 10. The simulations follow the non-equilibrium chemistry/cooling of primordial gas and the conversion of the gas into metal-free stars. To quantify the radiative feedback, we compare a simulation in which stars emit both molecular hydrogen dissociating and hydrogen/helium ionizing radiation with a simulation in which stars emit only dissociating radiation, and with a simulation in which stars remain dark. Photodissociation and -ionization exert a strong negative feedback on the assembly of the simulated galaxy. Gas condensation is strongly impeded, and star formation is strongly suppressed in comparison with the simulation in which stars remain dark. The feedback on the gas implies a suppression of the central dark matter densities in the minihalo progenitor by factors of up to a few, which is a significant deviation from the singular isothermal density profile characterizing the dark matter distribution in the absence of radiative feedback. The evolution of gas densities, star formation rates, and the distribution of dark matter becomes insensitive to the inclusion of dissociating radiation in the late stages of the minihalo assembly, and it becomes insensitive to the inclusion of ionizing radiation once the minihalo turns into an atomically cooling galaxy. The formation of an extended disk inside the dwarf galaxy is a robust outcome not affected by the inclusion of radiation. We estimate that dwarf galaxies such as simulated here will be among the faintest galaxies the upcoming James Webb Space Telescope will detect. Our conclusions are subject to our neglect of feedback from supernovae and chemical enrichment as well as to cosmic variance. [abridged]

preprint2012arXiv

Simulations of Accretion Powered Supernovae in the Progenitors of Gamma Ray Bursts

Observational evidence suggests a link between long duration gamma ray bursts (LGRBs) and Type Ic supernovae. Here, we propose a potential mechanism for Type Ic supernovae in LGRB progenitors powered solely by accretion energy. We present spherically-symmetric hydrodynamic simulations of the long-term accretion of a rotating gamma-ray burst progenitor star, a "collapsar," onto the central compact object, which we take to be a black hole. The simulations were carried out with the adaptive mesh refinement code FLASH in one spatial dimension and with rotation, an explicit shear viscosity, and convection in the mixing length theory approximation. Once the accretion flow becomes rotationally supported outside of the black hole, an accretion shock forms and traverses the stellar envelope. Energy is carried from the central geometrically thick accretion disk to the stellar envelope by convection. Energy losses through neutrino emission and nuclear photodisintegration are calculated but do not seem important following the rapid early drop of the accretion rate following circularization. We find that the shock velocity, energy, and unbound mass are sensitive to convective efficiency, effective viscosity, and initial stellar angular momentum. Our simulations show that given the appropriate combinations of stellar and physical parameters, explosions with energies ~5x10^50 ergs, velocities 3000 km/s, and unbound material masses >6 solar masses are possible in a rapidly rotating 16 solar mass main sequence progenitor star. Further work is needed to constrain the values of these parameters, to identify the likely outcomes in more plausible and massive LRGB progenitors, and to explore nucleosynthetic implications.

preprint2012arXiv

Star Formation in the First Galaxies I: Collapse Delayed by Lyman-Werner Radiation

We investigate the process of metal-free star formation in the first galaxies with a high-resolution cosmological simulation. We consider the cosmologically motivated scenario in which a strong molecule-destroying Lyman-Werner (LW) background inhibits effective cooling in low-mass haloes, delaying star formation until the collapse or more massive haloes. Only when molecular hydrogen (H2) can self-shield from LW radiation, which requires a halo capable of cooling by atomic line emission, will star formation be possible. To follow the formation of multiple gravitationally bound objects, at high gas densities we introduce sink particles which accrete gas directly from the computational grid. We find that in a 1 Mpc^3 (comoving) box, runaway collapse first occurs in a 3x10^7 M_sun dark matter halo at z~12 assuming a background intensity of J21=100. Due to a runaway increase in the H2 abundance and cooling rate, a self-shielding, supersonically turbulent core develops abruptly with ~10^4 M_sun in cold gas available for star formation. We analyze the formation of this self-shielding core, the character of turbulence, and the prospects for star formation. Due to a lack of fragmentation on scales we resolve, we argue that LW-delayed metal-free star formation in atomic cooling haloes is very similar to star formation in primordial minihaloes, although in making this conclusion we ignore internal stellar feedback. Finally, we briefly discuss the detectability of metal-free stellar clusters with the James Webb Space Telescope.

preprint2012arXiv

The First Galaxies: Assembly with Black Hole Feedback

We study how the first galaxies were assembled under feedback from the accretion onto a central black hole (BH) that is left behind by the first generation of metal-free stars through self-consistent, cosmological simulations. X-ray radiation from the accretion of gas onto BH remnants of Population III (Pop III) stars, or from high-mass X-ray binaries (HMXBs), again involving Pop III stars, influences the mode of second generation star formation. We track the evolution of the black hole accretion rate and the associated X-ray feedback starting with the death of the Pop III progenitor star inside a minihalo and following the subsequent evolution of the black hole as the minihalo grows to become an atomically cooling galaxy. We find that X-ray photoionization heating from a stellar-mass BH is able to quench further star formation in the host halo at all times before the halo enters the atomic cooling phase. X-ray radiation from a HMXB, assuming a luminosity close to the Eddington value, exerts an even stronger, and more diverse, feedback on star formation. It photoheats the gas inside the host halo, but also promotes the formation of molecular hydrogen and cooling of gas in the intergalactic medium and in nearby minihalos, leading to a net increase in the number of stars formed at early times. Our simulations further show that the radiative feedback from the first BHs may strongly suppress early BH growth, thus constraining models for the formation of supermassive BHs.

preprint2012arXiv

The Source Density and Observability of Pair-Instability Supernovae from the First Stars

Theoretical models predict that some of the first stars ended their lives as extremely energetic pair-instability supernovae (PISNe). With energies approaching 10^53 ergs, these supernovae are expected to be within the detection limits of the upcoming James Webb Space telescope (JWST), allowing observational constraints to be placed on the properties of the first stars. We estimate the source density of PISNe using a semi-analytic halo mass function based approach, accounting for the effects of feedback from star formation on the PISN rates using cosmological simulations. We estimate an upper limit of ~0.2 PISNe per JWST field of view at any given time. Feedback can reduce this rate significantly, e.g., lowering it to as little as one PISN per 4000 JWST fields of view for the most pessimistic explosion models. We also find that the main obstacle to observing PISNe from the first stars is their scarcity, not their faintness; exposures longer than a few times 10^4s will do little to increase the number of PISNe found. Given this we suggest a mosaic style search strategy for detecting PISNe from the first stars. Even rather high redshift PISNe are unlikely to be missed by moderate exposures, and a large number of pointings will be required to ensure a detection.

preprint2011arXiv

The First Galaxies: Assembly of Disks and Prospects for Direct Detection

The James Webb Space Telescope (JWST) will enable observations of galaxies at redshifts z > 10 and hence allow to test our current understanding of structure formation at very early times. Previous work has shown that the very first galaxies inside halos with virial temperatures T < 10^4 K and masses M < 10^8 M_sun at z > 10 are probably too faint, by at least one order of magnitude, to be detected even in deep exposures with JWST. The light collected with JWST may therefore be dominated by radiation from galaxies inside ten times more massive halos. We use cosmological zoomed smoothed particle hydrodynamics simulations to investigate the assembly of such galaxies and assess their observability with JWST. We compare two simulations that are identical except for the inclusion of non-equilibrium H/D chemistry and radiative cooling by molecular hydrogen. In both simulations a large fraction of the halo gas settles in two nested, extended gas disks which surround a compact massive gas core. The presence of molecular hydrogen allows the disk gas to reach low temperatures and to develop marked spiral structure but does not qualitatively change its stability against fragmentation. We post-process the simulated galaxies by combining idealized models for star formation with stellar population synthesis models to estimate the luminosities in nebular recombination lines as well as in the ultraviolet continuum. We demonstrate that JWST will be able to constrain the nature of the stellar populations in galaxies such as simulated here based on the detection of the He1640 recombination line. Extrapolation of our results to halos with masses both lower and higher than those simulated shows that JWST may find up to a thousand star-bursting galaxies in future deep exposures of the z > 10 universe.

preprint2011arXiv

Two-stream-like instability in dilute hot relativistic beams and astrophysical relativistic shocks

Relativistic collisionless shocks are believed to be efficient particle accelerators. Nonlinear outcome of the interaction of accelerated particles that run ahead of the shock, the so-called "precursor", with the unperturbed plasma of the shock upstream, is thought to facilitate additional acceleration of these particles and to possibly modify the hydrodynamic structure of the shock. We explore here the linear growth of kinetic modes appearing in the precursor-upstream interaction in relativistic shocks propagating in non and weakly magnetized plasmas: electrostatic two-stream parallel mode and electrostatic oblique modes. These modes are of particular interest because they are the fastest growing modes known in this type of system. Using a simplified distribution function for a dilute ultra-relativistic beam that is relativistically hot in its own rest frame, yet has momenta that are narrowly collimated in the frame of the cold upstream plasma into which it propagates, we identify the fastest growing mode in the full $k$-space and calculate its growth rate. We consider all types of plasma (pairs and ions-electrons) and beam (charged and charge-neutral). We find that unstable electrostatic modes are present in any type of plasma and for any shock parameters. We further find that two modes, one parallel ($k_\perp=0$) and the other one oblique ($k_\perp \sim k_\|$), are competing for dominance and that either one may dominate the growth rate in different regions of the phase space. The dominant mode is determined mostly by the perpendicular spread of the accelerated particle momenta in the upstream frame, which reflects the shock Lorentz factor. The parallel mode becomes more dominant in shocks with lower Lorentz factors (i.e., with larger momentum spreads). We briefly discuss possible implications of our results for external shocks in gamma-ray burst sources.

preprint2010arXiv

Fragmentation in the First Galaxies

We study the gravitational fragmentation of cold accretion streams flowing into a typical first galaxy. We use a one-zone hydrodynamical model to examine the thermal evolution of the gas entering a 10^8 M_sun DM halo at z=10. The goal is to find the expected fragmentation mass scale and thus a characteristic mass of the first population of stars to form by shock fragmentation at high redshift. Our model accurately describes the chemical and thermal evolution of the gas as we are specifically concerned with how the cooling of the gas alters its fragmentation properties. We find there to be a sharp drop in the fragmentation mass at a metallicity of ~10^-4 Z_sun when a strong molecule destroying, LW background is present. However, If molecules can efficiently form, they dominate the cooling at T < 10^4 K, demonstrating no 'critical metallicity'. Dust grains are not included in our chemical model, but we argue their inclusion would not significantly the results. We also find that this physical scenario allows for the formation of a cluster of solar mass fragments, or a single 10^4 M_sun fragment, possibly the precursors to primeval clusters and SMBHs. Lastly, we conclude that the usual assumption of isobaricity for galactic shocks breaks down in gas of sufficiently high metallicity, suggesting that metal cooling may lead to thermal instabilities.

preprint2010arXiv

Nuclear Star Clusters from Clustered Star Formation

Photometrically distinct nuclear star clusters (NSCs) are common in late-type-disk and spheroidal galaxies. The formation of NSCs is inevitable in the context of normal star formation in which a majority of stars form in clusters. A young, mass-losing cluster embedded in an isolated star-forming galaxy remains gravitationally bound over a period determined by its initial mass and the galactic tidal field. The cluster migrates radially toward the center of the galaxy and becomes integrated in the NSC if it reaches the center. The rate at which the NSC grows by accreting young clusters can be estimated from empirical cluster formation rates and dissolution times. We model cluster migration and dissolution and find that the NSCs in late-type disks and in spheroidals could have assembled from migrating clusters. The resulting stellar nucleus contains a small fraction of the stellar mass of the galaxy; this fraction is sensitive to the high-mass truncation of the initial cluster mass function (ICMF). The resulting NSC masses are consistent with the observed values, but generically, the final NSCs are surrounded by a spatially more extended excess over the inward-extrapolated exponential (or Sersic) law of the outer galaxy. We suggest that the excess can be related to the pseudobulge phenomenon in disks, though not all of the pseudobulge mass assembles this way. Comparison with observed NSC masses can be used to constrain the truncation mass scale of the ICMF and the fraction of clusters suffering prompt dissolution. We infer truncation mass scales of <~ 10^6 M_sun (>~ 10^5 M_sun) without (with 90%) prompt dissolution.

preprint2009arXiv

Accretion onto Seed Black Holes in the First Galaxies

The validity of the hypothesis that the massive black holes in high redshift quasars grew from stellar-sized "seeds" is contingent on a seed's ability to double its mass every few ten million years. This requires that the seed accrete at approximately the Eddington-limited rate. In the specific case of radiatively efficient quasiradial accretion in a metal-poor protogalactic medium, for which the Bondi accretion rate is often prescribed in cosmological simulations of massive black hole formation, we examine the effects of the radiation emitted near the black hole's event horizon on the structure of the surrounding gas flow. We find that the radiation pressure from photoionization significantly reduces the steady-state accretion rate and renders the quasiradial accretion flow unsteady and inefficient. The time-averaged accretion rates are a small fraction of the Eddington-limited accretion rate for Thomson scattering. The pressure of Ly-alpha photons trapped near the HII region surrounding the black hole may further attenuate the inflow. These results suggest that an alternative to quasiradial, radiatively efficient Bondi-like accretion should be sought to explain the rapid growth of quasar-progenitor seed black holes.

preprint2009arXiv

Collapsar Accretion and the Gamma-Ray Burst X-Ray Light Curve

We present axisymmetric hydrodynamical simulations of the long-term accretion of a rotating GRB progenitor star, a "collapsar," onto the central compact object. The simulations were carried out with the adaptive mesh refinement code FLASH in two spatial dimensions and with an explicit shear viscosity. The evolution of the central accretion rate exhibits phases reminiscent of the long GRB gamma-ray and X-ray light curve, which lends support to the proposal that the luminosity is modulated by the central accretion rate. After a few tens of seconds, an accretion shock sweeps outward through the star. The formation and outward expansion of the accretion shock is accompanied with a sudden and rapid power-law decline in the central accretion rate Mdot ~ t^{-2.8}, which resembles the L_X ~ t^{-3} decline observed in the X-ray light curves. The collapsed, shock-heated stellar envelope settles into a thick, low-mass equatorial disk embedded within a massive, pressure-supported atmosphere. After a few hundred seconds, the inflow of low-angular-momentum material in the axial funnel reverses into an outflow from the surface of the thick disk. Meanwhile, the rapid decline of the accretion rate slows down, or even settles a in steady state with Mdot ~ 5x10^{-5} Msun/s, which resembles the "plateau" phase in the X-ray light curve. While the duration of the "prompt" phase depends on the resolution in our simulations, we provide an analytical model taking into account neutrino losses that estimates the duration to be ~20 s. The model suggests that the steep decline in GRB X-ray light curves is triggered by the circularization of the infalling stellar envelope at radii where the virial temperature is below ~10^{10} K, such that neutrino cooling shuts off and an outward expansion of the accretion shock becomes imminent.

preprint2008arXiv

Aspherical Core-Collapse Supernovae in Red Supergiants Powered by Nonrelativistic Jets

We explore the observational characteristics of jet-driven supernovae by simulating bipolar-jet-driven explosions in a red supergiant progenitor. We present results of four models in which we hold the injected kinetic energy at a constant $10^{51}$ ergs across all jet models but vary the specific characteristics of the jets to explore the influence of the nature of jets on the structure of the supernova ejecta. We evolve the explosions past shock-breakout and into quasi-homologous expansion of the supernova envelope into a red supergiant wind. The oppositely-directed, nickel-rich jets give a large-scale asymmetry that may account for the non-spherical excitation and substructure of spectral lines such as H$α$ and He I 10830Å. Jets with a large fraction of kinetic to thermal energy punch through the progenitor envelope and give rise to explosions that would be observed to be asymmetric from the earliest epochs, inconsistent with spectropolarimetric measurements of Type II supernovae. Jets with higher thermal energy fractions result in explosions that are roughly spherical at large radii but are significantly elongated at smaller radii, deep inside the ejecta, in agreement with the polarimetric observations. We present shock breakout light curves that indicate that strongly aspherical shock breakouts are incompatible with recent {\it GALEX} observations of shock breakout from red supergiant stars. Comparison with observations indicates that jets must deposit their kinetic energy efficiently throughout the ejecta while in the hydrogen envelope. Thermal energy-dominated jets satisfy this criterion and yield many of the observational characteristics of Type II supernovae.

preprint2007arXiv

An Eccentric Circumbinary Accretion Disk and the Detection of Binary Massive Black Holes

We present a two-dimensional grid-based hydrodynamic simulation of a thin, viscous, locally-isothermal corotating disk orbiting an equal-mass Newtonian binary point mass on a fixed circular orbit. We study the structure of the disk after multiple viscous times. The binary maintains a central hole in the viscously-relaxed disk with radius equal to about twice the binary semimajor axis. Disk surface density within the hole is reduced by orders of magnitude relative to the density in the disk bulk. The inner truncation of the disk resembles the clearing of a gap in a protoplanetary disk. An initially circular disk becomes elliptical and then eccentric. Disturbances in the disk contain a component that is stationary in the rotating frame in which the binary is at rest; this component is a two-armed spiral density wave. We measure the distribution of the binary torque in the disk and find that the strongest positive torque is exerted inside the central low-density hole. We make connection with the linear theory of disk forcing at outer Lindblad resonances (OLRs) and find that the measured torque density distribution is consistent with forcing at the 3:2 (m=2) OLR, well within the central hole. We also measure the time dependence of the rate at which gas accretes across the hole and find quasi-periodic structure. We discuss implications for variability and detection of active galactic nuclei containing a binary massive black hole.

preprint2007arXiv

Cluster Merger Shock Constraints on Particle Acceleration and Nonthermal Pressure in the Intracluster Medium

X-ray observations of galaxy cluster merger shocks can be used to constrain nonthermal processes in the intracluster medium (ICM). The presence of nonthermal pressure components in the ICM, as well as the shock acceleration of particles and their escape, all affect shock jump conditions in distinct ways. Therefore, these processes can be constrained using X-ray surface brightness and temperature maps of merger shock fronts. Here we use these observations to place constraints on particle acceleration efficiency in intermediate Mach number (M ~ 2-3) shocks and explore the potential to constrain the contribution of nonthermal components (e.g., cosmic rays, magnetic field, and turbulence) to ICM pressure in cluster outskirts. We model the hydrodynamic jump conditions in merger shocks discovered in the galaxy clusters A520 (M ~ 2) and 1E 0657-56 (M ~ 3) using a multifluid model comprised of a thermal plasma, a nonthermal plasma, and a magnetic field. Based on the published X-ray spectroscopic data alone, we find that the fractional contribution of cosmic rays accelerated in these shocks is lower than about 10% of the shock downstream pressure. Current observations do not constrain the fractional contribution of nonthermal components to the pressure of the undisturbed shock upstream. Future X-ray observations, however, have the potential to either detect particle acceleration in these shocks through its effect on the shock dynamics, or to place a lower limit on the nonthermal pressure contributions in the undisturbed ICM. We briefly discuss implications for models of particle acceleration in collisionless shocks and the estimates of galaxy cluster masses derived from X-ray and Sunyaev-Zel'dovich effect observations.

preprint2005arXiv

Cluster Merger Variance and the Luminosity Gap Statistic

The presence of multiple luminous galaxies in clusters can be explained by the finite time over which a galaxy sinks to the center of the cluster and merges with the the central galaxy. The simplest measurable statistic to quantify the dynamical age of a system of galaxies is the luminosity (magnitude) gap, which is the difference in photometric magnitude between the two most luminous galaxies. We present a simple analytical estimate of the luminosity gap distribution in groups and clusters as a function of dark matter halo mass. The luminosity gap is used to define "fossil" groups; we expect the fraction of fossil systems to exhibit a strong and model-independent trend with mass: ~1-3% of massive clusters and ~5-40% of groups should be fossil systems. We also show that, on cluster scales, the observed intrinsic scatter in the central galaxy luminosity-halo mass relation can be ascribed to dispersion in the merger histories of satellites within the cluster. We compare our predictions to the luminosity gap distribution in a sample of 730 clusters in the Sloan Digital Sky Survey C4 Catalog and find good agreement. This suggests that theoretical excursion set merger probabilities and the standard theory of dynamical segregation are valid on cluster scales.

preprint2005arXiv

Dissipationless Merging and the Assembly of Central Galaxies

We reanalyze the galaxy-mass correlation function measured by the Sloan Digital Sky Survey to obtain host dark matter halo masses at galaxy and galaxy group scales. We extend the data to galaxy clusters in the 2MASS catalog and study the relation between central galaxy luminosity and halo mass. While the central galaxy luminosity scales as ~M^{0.7-0.8} at low masses, the relation flattens to ~M^{<0.3} above ~4x10^{13} M_sun. The total luminosity of galaxies in the halo, however, continues to grow as a power-law ~M^{0.8-0.9}. Starting from the hypothesis that the central galaxies grow by merging ("galactic cannibalism"), we develop a simple model for the evolution of their luminosities as a consequence of the accretion of satellite galaxies. The luminosity-mass relation flattens when the time scale on which dynamical friction induces orbital decay in the satellite galaxies exceeds the age of the dark matter halo. Then, the growth of the central galaxy is suppressed as it can cannibalize only the rare, massive satellite galaxies. The model takes the dependence of the total luminosity of galaxies in a halo on its mass and the global galaxy luminosity function as input, and reproduces the observed central galaxy luminosity-mass relation over three decades in halo mass, (10^{12}-10^{15}) M_sun. The success of the model suggests that gas cooling and subsequent star formation did not play an important role in the final assembly of central galaxies from sub-L_star precursors.

preprint2005arXiv

Massive Black Hole Binary Evolution

Coalescence of binary supermassive black holes (SBHs) would constitute the strongest sources of gravitational waves to be observed by LISA. While the formation of binary SBHs during galaxy mergers is almost inevitable, coalescence requires that the separation between binary components first drop by a few orders of magnitude, due presumably to interaction of the binary with stars and gas in a galactic nucleus. This article reviews the observational evidence for binary SBHs and discusses how they would evolve. No completely convincing case of a bound, binary SBH has yet been found, although a handful of systems (e.g. interacting galaxies; remnants of galaxy mergers) are now believed to contain two SBHs at projected separations of roughly one kpc. N-body studies of binary evolution in gas-free galaxies have reached large enough particle numbers to reproduce the slow, "diffusive" refilling of the binary's loss cone that is believed to characterize binary evolution in real galactic nuclei. While some of the results of these simulations -- e.g. the binary hardening rate and eccentricity evolution -- are strongly $N$-dependent, others -- e.g. the "damage" inflicted by the binary on the nucleus -- are not. Luminous early-type galaxies often exhibit depleted cores with masses of 1-2 times the mass of their nuclear SBHs, consistent with the predictions of the binary model. Studies of the interaction of massive binaries with gas are still in their infancy, although much progress is expected in the near future. Binary coalescence has a large influence on the spins of SBHs, even for mass ratios as extreme as 10:1, and evidence of spin-flips may have been observed.

preprint2005arXiv

The ACS Virgo Cluster Survey IX: The Color Distributions of Globular Cluster Systems in Early-Type Galaxies

We present the color distributions of globular cluster (GC) systems for 100 Virgo cluster early-type galaxies observed in the ACS Virgo Cluster Survey. The color distributions of individual GC systems are consistent with continuous trends across galaxy luminosity, color, and stellar mass. On average, almost all galaxies possess a component of metal-poor GCs, with the average fraction of metal-rich GCs ranging from 15 to 60%. The colors of both subpopulations correlate with host galaxy luminosity and color, with the red GCs having a steeper slope. To convert color to metallicity, we also introduce a preliminary (g-z)-[Fe/H] relation calibrated to Galactic, M49 and M87 GCs. This relation is nonlinear with a steeper slope for [Fe/H] < -0.8. As a result, the metallicities of the metal-poor and metal-rich GCs vary similarly with respect to galaxy luminosity and stellar mass, with relations of [Fe/H]_MP ~ L^0.16 ~ M_star^0.17 and [Fe/H]_MR ~ L^0.26 ~ M_star^0.22, respectively. Although these relations are shallower than the mass-metallicity relation predicted by wind models and observed for dwarf galaxies, they are very similar to the mass-metallicity relation for star forming galaxies in the same mass range. The offset between the two GC populations varies slowly (~ M_star^0.05) and is approximately 1 dex across three orders of magnitude in mass, suggesting a nearly universal amount of enrichment between the formation of the two populations of GCs. We also find that although the metal-rich GCs show a larger dispersion in color, it is the *metal-poor GCs* that have an equal or larger dispersion in metallicity. Like the color-magnitude relation, these relations derived from globular clusters present stringent constraints on the formation and evolution of early-type galaxies. (Abridged)

preprint2005arXiv

The ACS Virgo Cluster Survey. VII. Resolving the Connection Between Globular Clusters and Ultra-Compact Dwarf Galaxies

We investigate the connection between globular clusters and ultra-compact dwarf galaxies (UCDs) by examining the properties of several compact objects associated with M87, all of which were previously classified as globular clusters. Combining imaging from the Hubble Space Telescope with ground-based Keck spectroscopy, we find two objects to have half-light radii, velocity dispersions and mass-to-light ratios that are consistent with the predictions of population synthesis models for old, metal-rich, luminous globular clusters. Three other objects are much larger, with half-light radii of approximately 20pc, and have V-band mass-to-light ratios in the range 6-9. These objects, which we consider to be UCDs, resemble the nuclei of nucleated dwarf elliptical galaxies in Virgo, having similar mass-to-light ratios, luminosities and colors. These UCDs are found to obey the extrapolated scaling relations of galaxies more closely than those of Galactic globular clusters. There appears to be a transition between the two types of stellar systems at a mass of about two million solar masses. If the UCDs are gravitationally bound, then we suggest that the presence of dark matter is the fundamental property distinguishing globular clusters from UCDs. More than half of the UCD candidates uncovered in the ACS Virgo Cluster Survey are associated with a single galaxy -- M87 -- which suggests that proximity to the Virgo center may be of critical importance for the formation of these objects. These results show that distinguishing bonafide UCDs from bright globular clusters requires a careful analysis of their detailed structural and dynamical properties, particularly their mass-to-light ratios. (ABRIDGED)

preprint2004arXiv

Consequences of gravitational radiation recoil

Coalescing binary black holes experience an impulsive kick due to anisotropic emission of gravitational waves. We discuss the dynamical consequences of the recoil accompanying massive black hole mergers. Recoil velocities are sufficient to eject most coalescing black holes from dwarf galaxies and globular clusters, which may explain the apparent absence of massive black holes in these systems. Ejection from giant elliptical galaxies would be rare, but coalescing black holes are displaced from the center and fall back on a time scale of order the half-mass crossing time. Displacement of the black holes transfers energy to the stars in the nucleus and can convert a steep density cusp into a core. Radiation recoil calls into question models that grow supermassive black holes from hierarchical mergers of stellar-mass precursors.