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 $\Lambda$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) $\alpha$, (2) closed-to-connected triple fraction (the so-called transitivity) $\tau_\Delta$, 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.

We demonstrate physics-blind and purely data-based nucleosynthetic archetype discovery from heterogeneous stellar elemental abundance measurements. Archetype extraction is performed via non-negative matrix factorization (NMF). Because stellar elemental abundance measurements can be incomplete and sparsely sampled, we perform combined low-rank NMF and matrix completion on the Hypatia Catalog and JINAbase to extract archetypal elemental abundance patterns, the so-called endmembers, and impute the missing measurements. We qualitatively compare the discovered nucleosynthetic archetypes to patterns anticipated from direct astrophysical nucleosynthetic calculations. Our tentative results indicate that the physics-blind analysis relying entirely on the incomplete chemical mixing in star formation can discern nucleosynthetic archetypes associated with asymptotic giant branch (AGB) stellar sources of the s-process, thermonuclear supernovae, and non-thermonuclear supernovae with an r-process.

We model dust formation in the core collapse supernova explosion SN 1987A by treating the gas-phase formation of dust grain nuclei as a chemical process. To compute the synthesis of fourteen species of grains we integrate a non-equilibrium network of nucleating and related chemical reactions and follow the growth of the nuclei into grains via accretion and coagulation. The effects of the radioactive cobalt, titanium, and sodium on the thermodynamics and chemistry of the ejecta are taken into account. The grain temperature, which we allow to differ from the gas temperature, affects the surface-tension-corrected evaporation rate. We also account for He$^+$, Ne$^+$, Ar$^+$, and O weathering. We combine our dust synthesis model with a crude prescription for anisotropic radioactive nickel dredge-up into the core ejecta, the so-called `nickel bubbles', to compute the total dust mass and molecular-species-specific grain size distribution. The total mass varies between $0.41\,M_\odot$ and $0.73\,M_\odot$, depending on the bubble shell density contrast. In the decreasing order of abundance, the grain species produced are: magnesia, silicon, forsterite, iron sulfide, carbon, silicon dioxide, alumina, and iron. The combined grain size distribution is a power law $dN/da\propto a^{-4.39}$. Early ejecta compaction by expanding radioactive nickel bubbles strongly enhances dust synthesis. This underscores the need for improved understanding of hydrodynamic transport and mixing over the entire pre-homologous expansion.

Due to its accuracy and generality, Monte Carlo radiative transfer (MCRT) has emerged as the prevalent method for Ly$\alpha$ radiative transfer in arbitrary geometries. The standard MCRT encounters a significant efficiency barrier in the high optical depth, diffusion regime. Multiple acceleration schemes have been developed to improve the efficiency of MCRT but the noise from photon packet discretization remains a challenge. The discrete diffusion Monte Carlo (DDMC) scheme has been successfully applied in state-of-the-art radiation hydrodynamics (RHD) simulations. Still, the established framework is not optimal for resonant line transfer. Inspired by the DDMC paradigm, we present a novel extension to resonant DDMC (rDDMC) in which diffusion in space and frequency are treated on equal footing. We explore the robustness of our new method and demonstrate a level of performance that justifies incorporating the method into existing Ly$\alpha$ codes. We present computational speedups of $\sim 10^2$-$10^6$ relative to contemporary MCRT implementations with schemes that skip scattering in the core of the line profile. This is because the rDDMC runtime scales with the spatial and frequency resolution rather than the number of scatterings - the latter is typically $\propto \tau_0$ for static media, or $\propto (a \tau_0)^{2/3}$ with core-skipping. We anticipate new frontiers in which on-the-fly Ly$\alpha$ radiative transfer calculations are feasible in 3D RHD. More generally, rDDMC is transferable to any computationally demanding problem amenable to a Fokker-Planck approximation of frequency redistribution.

The physics of star formation at its extreme, in the nuclei of the densest and the most massive star clusters in the universe - potential massive black hole nurseries - has for decades eluded scrutiny. Spectroscopy of these systems has been scarce, whereas theoretical arguments suggest that radiation pressure on dust grains somehow inhibits star formation. Here, we harness an accelerated Monte Carlo radiation transport scheme to report a radiation hydrodynamical simulation of super star cluster formation in turbulent clouds. We find that radiation pressure reduces the global star formation efficiency by 30-35%, and the star formation rate by 15-50%, both relative to a radiation-free control run. Overall, radiation pressure does not terminate the gas supply for star formation and the final stellar mass of the most massive cluster is $\sim 1.3\times10^6\,M_\odot$. The limited impact as compared to in idealized theoretical models is attributed to a radiation-matter anti-correlation in the supersonically turbulent, gravitationally collapsing medium. In isolated regions outside massive clusters, where the gas distribution is less disturbed, radiation pressure is more effective in limiting star formation. The resulting stellar density at the cluster core is $\ge 10^{8}\,M_{\odot}\,\textrm{pc}^{-3}$, with stellar velocity dispersion $\gtrsim 70\,\text{km}\,\text{s}^{-1}$. We conclude that the super star cluster nucleus is propitious to the formation of very massive stars via dynamical core collapse and stellar merging. We speculate that the very massive star may avoid the claimed catastrophic mass loss by continuing to accrete dense gas condensing from a gravitationally-confined ionized phase.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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).

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}$.

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]

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.

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.

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.

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.

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.

The association of long-duration gamma-ray bursts (LGRBs) with Type Ic supernovae presents a challenge to supernova explosion models. In the collapsar model for LGRBs, gamma rays are produced in an ultrarelativistic jet launching from the magnetosphere of the black hole that forms in the aftermath of the collapse of a rotating progenitor star. The jet is collimated along the star's rotation axis, but the concomitant luminous supernova should be relatively--though certainly not entirely--spherical, and should synthesize a substantial mass of 56Ni. Our goal is to provide a qualitative assessment of the possibility that accretion of the progenitor envelope onto the black hole, which powers the LGRB, could also deposit sufficient energy and nickel mass in the envelope to produce a luminous supernova. For this, the energy dissipated near the black hole during accretion must be transported outward, where it can drive a supernova-like shockwave. Here we suggest that the energy is transported by convection and develop an analytical toy model, relying on global mass and energy conservation, for the dynamics of stellar collapse. The model suggests that a ~10,000 km/s shock can be driven into the envelope and that ~10^51 erg explosions are possible. The efficiency with which the accretion energy is being transferred to the envelope is governed by the competition of advection and convection at distances ~100-1,000 km from the black hole and is sensitive to the values of the convective mixing length, the magnitude of the effective viscous stress, and the specific angular momentum of the infalling envelope. Substantial masses of 56Ni may be synthesized in the convective accretion flow over the course of tens of seconds from the initial circularization of the infalling envelope around the black hole. The synthesized nickel is convectively mixed with a much larger mass of unburned ejecta.

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.

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.

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.

Shock breakout is the earliest, readily-observable emission from a core-collapse supernova explosion. Observing supernova shock breakout may yield information about the nature of the supernova shock prior to exiting the progenitor and, in turn, about the core-collapse supernova mechanism itself. X-ray Outburst 080109, later associated with SN 2008D, is a very well-observed example of shock breakout from a core-collapse supernova. Despite excellent observational coverage and detailed modeling, fundamental information about the shock breakout, such as the radius of breakout and driver of the light curve time scale, is still uncertain. The models constructed for explaining the shock breakout emission from SN 2008D all assume spherical symmetry. We present a study of the observational characteristics of {\it aspherical} shock breakout from stripped-envelope core-collapse supernovae. We conduct two-dimensional, jet-driven supernova simulations from stripped-envelope progenitors and calculate the resulting shock breakout X-ray spectra and light curves. The X-ray spectra evolve significantly in time as the shocks expand outward and are not well-fit by single-temperature and radius black bodies. The time scale of the X-ray burst light curve of the shock breakout is related to the shock crossing time of the progenitor, not the much shorter light crossing time that sets the light curve time scale in spherical breakouts. This could explain the long shock breakout light curve time scale observed for XRO 080109/SN 2008D.

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.

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.

As the closest example of a galactic nucleus, the Galactic center presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environs. Detailed studies of stellar dynamics deep in the potential well of a galaxy, with exisiting and future large ground-based telescopes, offer several exciting directions in the coming decade. First, it will be possible to obtain precision measurements of the Galaxy's central potential, providing both a unique test of General Relativity (GR) and a detection of the extended dark matter distribution that is predicted to exist around the SMBH. Tests of gravity have not previously been possible on the mass scale of a SMBH. Similarly, only upper limits on the extended matter distribution on small scales currently exist; detection of dark matter on these scales is an important test of Lambda-CDM and the detection of stellar remnants would reveal a population that may dominate the stellar dynamics on the smallest scales. Second, our detailed view of the SMBH and its local gas and stellar environment provides insight into how SMBHs at the centers of galaxies form, grow and interact with their environs as well as on the exotic processes at work in the densest stellar clusters in the Universe. The key questions, still unanswered, of when and how SMBHs formed in the early universe, and the myriad ways in which feedback from SMBHs can affect structure formation, can be informed by directly observing the physical processes operating at the SMBH.

We study the formation of disk-dominated galaxies in a LCDM universe. Their existence is considered to be a challenge for the LCDM cosmology, because galaxy mergers isotropize stellar disks and trigger angular momentum transport in gas disks, thus fostering the formation of central stellar spheroids. Here, we postulate that the formation of stellar spheroids from gasrich disks is controlled by two parameters that characterize galaxy mergers, the mass ratio of merging dark matter halos, and the virial velocity of the larger merging halo. We utilize merger histories generated from realizations of the cosmological density field to calculate the fraction of dark matter halos that have avoided spheroid formation, and compare the derived statistics with the spheroid occupation fractions in surveys of nearby galaxies. We find, for example, that the survival rate of disk-dominated galaxies in LCDM is just high enough to explain the observed fractional representation of disk-dominated galaxies in the universe if the only mergers which lead to central spheroid formation are those with mass ratios M2/M1 > 0.3 and virial velocities Vvir,1 > 55 km/s. We discuss the physical origin of this criterion, and show that the dependence of the disk-dominated fraction on galaxy mass provides a further test of the merger hypothesis. [For additional details, see, Koda et al. (2009).]

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.

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$\alpha$ and He I 10830\AA. 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.

We present the first results from two-dimensional simulations of radiatively-efficient accretion of metal-free gas onto intermediate-mass black holes. We fix the shape of the spectral energy distribution of the radiation produced near the event horizon and study the structure of the irradiated low-angular-momentum accretion flow over three orders of magnitude in radius from the black hole, 10^{14}-10^{17} cm for a 100 M_sun black hole. The luminosity of the central source is made to be proportional to the rate at which gas accretes across the inner boundary, which we set just inside the sonic radius. We find that photoionization heating and radiation pressure modify the structure of the flow. When the ambient gas density is 10^7 cm^{-3}, accretion is intermittent and on average reduced to 32% of the Eddington-limited rate, two orders of magnitude below the "Bondi" rate evaluated ignoring radiation, in agreement with simplified theoretical models. Even if the vicinity of the black hole is supplied with high density gas, accretion is rendered inefficient through heating and radiation pressure.

One of the fundamental astrobiology questions is how life has formed in our Solar System. In this context the formation and stability of abiotic organic molecules such as CH4, formic acid and amino acids, is important for understanding how organic material has formed and survived shocks and energetic particle impact from winds in the early Solar System. Shock waves have been suggested as a plausible scenario to create chondrules, small meteoritic components that have been completely molten by energetic events such as shocks and high velocity particle impacts. We study here the formation and destruction of certain gas-phase molecules such as methane and water during such shock events and compare the chemical timescales with the timescales for shocks arising from gravitational instabilities in a protosolar nebula.

We report new precision measurements of the properties of our Galaxy's supermassive black hole. Based on astrometric (1995-2007) and radial velocity (2000-2007) measurements from the W. M. Keck 10-meter telescopes, a fully unconstrained Keplerian orbit for the short period star S0-2 provides values for Ro of 8.0+-0.6 kpc, M_bh of 4.1+-0.6x10^6 Mo, and the black hole's radial velocity, which is consistent with zero with 30 km/s uncertainty. If the black hole is assumed to be at rest with respect to the Galaxy, we can further constrain the fit and obtain Ro = 8.4+-0.4 kpc and M_bh = 4.5+-0.4x10^6 Mo. More complex models constrain the extended dark mass distribution to be less than 3-4x10^5 Mo within 0.01 pc, ~100x higher than predictions from stellar and stellar remnant models. For all models, we identify transient astrometric shifts from source confusion and the assumptions regarding the black hole's radial motion as previously unrecognized limitations on orbital accuracy and the usefulness of fainter stars. Future astrometric and RV observations will remedy these effects. Our estimates of Ro and the Galaxy's local rotation speed, which it is derived from combining Ro with the apparent proper motion of Sgr A*, (theta0 = 229+-18 km/s), are compatible with measurements made using other methods. The increased black hole mass found in this study, compared to that determined using projected mass estimators, implies a longer period for the innermost stable orbit, longer resonant relaxation timescales for stars in the vicinity of the black hole and a better agreement with the M_bh-sigma relation.

We investigate the interaction of nonthermal ions (protons and nuclei) accelerated in an ultrarelativistic blastwave with the pre-existing magnetic field of the medium into which the blastwave propagates. While particle acceleration processes such as diffusive shock acceleration can accelerate ions and electrons, the accelerated electrons suffer larger radiative losses. Under certain conditions, the ions can attain higher energies and reach farther ahead of the shock than the electrons, and so the nonthermal particles can be partially charge-separated. To compensate for the charge separation, the upstream plasma develops a return current, which, as it flows across the magnetic field, drives transverse acceleration of the upstream plasma and a growth of density contrast in the shock upstream. If the density contrast is strong by the time the fluid is shocked, vorticity is generated at the shock transition. The resulting turbulence can amplify the post-shock magnetic field to the levels inferred from gamma-ray burst afterglow spectra and light curves. Therefore, since the upstream inhomogeneities are induced by the ions accelerated in the shock, they are generic even if the blastwave propagates into a medium of uniform density. We speculate about the global structure of the shock precursor, and delineate several distinct physical regimes that are classified by an increasing distance from the shock and, correspondingly, a decreasing density of nonthermal particles that reach that distance.

The merging galaxy cluster 1E 0657-56, known as the "bullet cluster," is one of the hottest clusters known. The X-ray emitting plasma exhibits bow-shock-like temperature and density jumps. The segregation of this plasma from the peaks of the mass distribution determined by gravitational lensing has been interpreted as a direct proof of collisionless dark matter. If the high shock speed inferred from the shock jump conditions equals the relative speed of the merging CDM halos, however, this merger is predicted to be such a rare event in a LCDM universe that observing it presents a possible conflict with the LCDM model. We examined this question using high resolution, 2D simulations of gas dynamics in cluster collisions to analyze the relative motion of the clusters, the bow shock, and the contact discontinuity, and relate these to the X-ray data for the bullet cluster. We find that the velocity of the fluid shock need not equal the relative velocity of the CDM components. An illustrative simulation finds that the present relative velocity of the CDM halos is 16% lower than that of the shock. While this conclusion is sensitive to the detailed initial mass and gas density profiles of the colliding clusters, such a decrease of the inferred halo relative velocity would significantly increase the likelihood of finding 1E 0657-56 in a LCDM universe. (Conference proceedings based on a poster at Bash Symposium 2007)

We study the formation of disk-dominated galaxies in a Lambda CDM universe. Their existence is considered to be a challenge for the Lambda CDM cosmology, because galaxy mergers isotropize stellar disks and trigger angular momentum transport in gas disks, thus fostering the formation of central stellar spheroids. Here, we postulate that the formation of stellar spheroids from gas-rich disks is controlled by two parameters that characterize galaxy mergers, the mass ratio of merging dark matter halos, and the virial velocity of the larger merging halo. We utilize merger histories generated from realizations of the cosmological density field to calculate the fraction of dark matter halos that have avoided spheroid formation, and compare the derived statistics with the spheroid occupation fractions in surveys of nearby galaxies. We find, for example, that the survival rate of disk-dominated galaxies in Lambda CDM is just high enough to explain the observed fractional representation of disk-dominated galaxies in the universe if the only mergers which lead to central spheroid formation are those with mass ratios M2/M1 > 0.3 and virial velocities Vvir,1 > 55 km/s. We discuss the physical origin of this criterion.

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.

Because conical segments of quasispherical ultrarelativistic blastwaves are causally disconnected on angular scales larger than the blastwave inverse Lorentz factor, astrophysical blastwaves can sustain initial anisotropy, imprinted by the process that drives the explosion, while they remain relativistic. We show that initial angular energy fluctuations in ultrarelativistic blastwaves imply a production of vorticity in the blastwave, and calculate the vortical energy production rate. In gamma-ray burst (GRB) afterglows, the number of vortical eddy turnovers as the shocked fluid crosses the blastwave shell is about unity for marginally nonlinear anisotropy. Thus the anisotropy must be nonlinear to explain the magnetic energy density inferred in measured GRB spectra.

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.

Shock waves driven in the intergalactic medium during the merging of galaxy clusters have been observed in X-ray imaging and spectroscopy. Fluid motions inferred from the shock strength and morphology can be compared to the cold dark matter (CDM) distribution inferred from gravitational lensing. A detailed reconstruction of the CDM kinematics, however, must take into account the nontrivial response of the fluid intracluster medium to the collisionless CDM motions. We have carried out two-dimensional simulations of gas dynamics in cluster collisions. We analyze the relative motion of the clusters, the bow shock wave, and the contact discontinuity and relate these to X-ray data. We focus on the "bullet cluster," 1E 0657-56, a near head-on collision of unequal-mass clusters, for which the gas density and temperature jumps across the prominent bow shock imply a high shock velocity 4,700 km/s. The velocity of the fluid shock has been widely interpreted as the relative velocity of the CDM components. This need not be the case, however. An illustrative simulation finds that the present relative velocity of the CDM halos is 16% lower than that of the shock. While this conclusion is sensitive to the detailed initial mass and gas density profile of the colliding clusters, such a decrease of the inferred halo relative velocity would increase the likelihood of finding 1E 0657-56 in a LambdaCDM universe.

Collisionless shocks are commonly argued to be the sites of cosmic ray (CR) acceleration. We study the influence of CRs on weakly magnetized relativistic collisionless shocks and apply our results to external shocks in gamma-ray burst (GRB) afterglows. The common view is that the transverse Weibel instability (TWI) generates a small-scale magnetic field that facilitates collisional coupling and thermalization in the shock transition. The TWI field is expected to decay rapidly, over a finite number of proton plasma skin depths from the transition. However, the synchrotron emission in GRB afterglows suggests that a strong and persistent magnetic field is present in the plasma that crosses the shock; the origin of this field is a key open question. Here we suggest that the common picture involving TWI demands revision. Namely, the CRs drive turbulence in the shock upstream on scales much larger than the skin depth. This turbulence generates a large-scale magnetic field that quenches TWI and produces a magnetized shock. The new field efficiently confines CRs and enhances the acceleration efficiency. The CRs modify the shocks in GRB afterglows at least while they remain relativistic. The origin of the magnetic field that gives rise to the synchrotron emission is plausibly in the CR-driven turbulence. We do not expect ultrahigh energy cosmic ray production in external GRB shocks.

(Abridged) The ACS Virgo Cluster Survey is an HST program to obtain high-resolution, g and z-band images for 100 early-type members of the Virgo Cluster, spanning a range of ~460 in blue luminosity. Based on this large, homogeneous dataset, we present a sharp upward revision in the frequency of nucleation in early-type galaxies brighter than M_B ~ -15 (66 < f_n < 82%), and find no evidence for nucleated dwarfs to be more concentrated to the center of Virgo than their non-nucleated counterparts. Resolved stellar nuclei are not present in galaxies brighter than M_B ~ -20.5, however, there is no clear evidence from the properties of the nuclei, or from the overall incidence of nucleation, for a change at M_B ~ -17.6, the traditional dividing point between dwarf and giant galaxies. On average, nuclei are ~3.5 mag brighter than a typical globular cluster and have a median half-light radius ~4.2 pc. Nuclear luminosities correlate with nuclear sizes and, in galaxies fainter than M_B ~ -17.6, nuclear colors. Comparing the nuclei to the "nuclear clusters" found in late-type spiral galaxies reveals a close match in terms of size, luminosity and overall frequency, pointing to a formation mechanism that is rather insensitive to the detailed properties of the host galaxy. The mean nuclear-to-galaxy luminosity ratio is indistinguishable from the mean SBH-to-bulge mass ratio, calculated in early-type galaxies with detected supermassive black holes (SBHs). We argue that compact stellar nuclei might be the low-mass counterparts of the SBHs detected in the bright galaxies, and that one should think in terms of "Central Massive Objects" -- either SBHs or compact stellar nuclei -- that accompany the formation of almost all early-type galaxies and contain a mean fraction ~0.3% of the total bulge mass.

The luminosity function of active galactic nuclei has been measured down to luminosities ~10^{42} erg/s in the soft and hard X-rays. Some fraction of this activity is associated with the accretion of the material liberated by the tidal disruption of stars by massive black holes. We estimate the contribution to the X-ray luminosity function from the tidal disruption process. While the contribution depends on a number poorly known parameters, it appears that it can account for the majority of X-ray selected AGN with soft or hard X-ray luminosities <~10^{43}-10^{44} erg/s. If this is correct, a picture emerges in which a significant portion of the X-ray luminosity function of AGN is comprised of sources powered by tidal-disruption at the faint end, while the sources at the bright end are powered by non-stellar accretion. Black holes with masses <~2x10^6 Msun could have acquired most of their present mass by an accretion of tidal debris. In view of the considerable theoretical uncertainty concerning the detailed shape of the light curves of tidal disruption events, we focus on power-law luminosity decay (as identified in candidate tidal disruption events), but we also discuss constant accretion rate models.

(Abridged) We present a detailed analysis of the morphology, isophotal parameters and surface brightness profiles for 100 early-type members of the Virgo Cluster, from dwarfs (M_B = -15.1 mag) to giants (M_B = -21.8 mag). Each galaxy has been imaged in two filters, closely resembling the Sloan g and z passbands, using the Advanced Camera for Surveys on board the Hubble Space Telescope. Dust and complex morphological structures are common, with kiloparsec-scale stellar disks, bars, and nuclear stellar disks seen in 60% of galaxies with intermediate luminosity (-20 < M_B < -17), and dust seen in 42% of galaxies brighter than M_B = -18.9 mag. Dust morphologies range from faint wisps and patches on tens of parsec scales, to regular, highly organized kpc-scale dust disks, often showing evidence of recent star formation. Surface brightness profiles and isophotal parameters are derived typically within 8 kpc from the center for the brightest galaxies, and 1.5 kpc for the faintest systems, with a resolution (FWHM) of 7 pc. Based on a parametrization of the surface brightness profiles in terms of a Sersic or core-Sersic model, we find that 1) there is no evidence of a bimodal behavior of the slope, gamma, of the profile in the innermost regions; 2) although the brightest galaxies have shallow inner profiles, the shallowest profiles (lowest gamma values) are found in faint dwarf systems; 3) the widely adopted separation of early-type galaxies between "core" and "power-law" types, which had originally been prompted by the claim of a clear bimodal distribution of gamma values, is untenable; and 4) there is no evidence of a structural dichothomy between dwarf and regular ellipticals.

We show that with a Next Generation Large Telescope one can detect the accelerated motions of ~100 stars orbiting the massive black hole at the Galactic center. The positions and velocities of these stars will be measured to astrometric and spectroscopic precision several times better than currently attainable enabling detailed measurements of the gravitational potential in the neighborhood of the massive black hole. We show that the monitoring of stellar motions with such a telescopes enables: (1) a measurement of the Galactic center distance R_0 to better than 0.1% accuracy, (2) a measurement of the extended matter distribution near the black hole, including that of the exotic dark matter, (3) a detection of general relativistic effects due to the black hole including the prograde precession of stars and possibly the black hole spin, and (4) a detection of gravitational encounters between monitored stars and stellar remnants that accumulate near the Galactic center. Such encounters probe the mass function of the remnants.

We discuss experiments achievable via monitoring of stellar dynamics near the massive black hole at the Galactic center with a next generation, extremely large telescope (ELT). Given the likely observational capabilities of an ELT and current knowledge of the stellar environment at the Galactic center, we synthesize plausible samples of stellar orbits around the black hole. We use the Markov Chain Monte Carlo method to evaluate the constraints that orbital monitoring places on the matter content near the black hole. Results are expressed as functions of the number N of stars with detectable orbital motions and the astrometric precision dtheta and spectroscopic precision dv at which stellar proper motions and radial velocities are monitored. For N = 100, dtheta = 0.5 mas, and dv = 10 km/s -- a conservative estimate of the capabilities of a 30 meter telescope -- the extended matter distribution enclosed by the orbits will produce measurable deviations from Keplerian motion if >1000 Msun is enclosed within 0.01 pc. The black hole mass and distance to the Galactic center will be measured to better than ~0.1%. Lowest-order relativistic effects, such as the prograde precession, will be detectable if dtheta < 0.5 mas. Higher-order effects, including frame dragging due to black hole spin, requires dtheta < 0.05 mas, or the favorable discovery of a compact, highly eccentric orbit. Finally, we calculate the rate at which monitored stars undergo detectable nearby encounters with background stars. Such encounters probe the mass function of stellar remnants that accumulate near the black hole. We find that ~30 encounters will be detected over a 10 yr baseline for dtheta = 0.5 mas.

We use HST/ACS imaging of 100 early-type galaxies in the ACS Virgo Cluster Survey to investigate the nature of diffuse star clusters (DSCs). Compared to globular clusters (GCs), these star clusters have moderately low luminosities (M_V > -8) and a broad distribution of sizes (3 < r_h < 30 pc), but they are principally characterized by their low mean surface brightnesses which can be more than three magnitudes fainter than a typical GC (mu_g > 20 mag arcsec^-2). The median colors of diffuse star cluster systems are red, 1.1 < g-z < 1.6, which is redder than metal-rich GCs and often as red as the galaxy itself. Most DSC systems thus have mean ages older than 5 Gyr or else have super-solar metallicities. We find that 12 galaxies in our sample contain a significant excess of diffuse star cluster candidates -- nine are lenticulars (S0s), and five visibly contain dust. We also find a substantial population of DSCs in the halo of the giant elliptical M49, associated with the companion galaxy VCC 1199. Most DSC systems appear to be both aligned with the galaxy light and associated with galactic disks, but at the same time many lenticulars do not host substantial DSC populations. Diffuse star clusters in our sample share similar characteristics to those identified in other nearby lenticular, spiral, and dwarf galaxies, and we suggest that DSCs preferentially form, survive, and coevolve with galactic disks. Their properties are broadly consistent with those of merged star cluster complexes, and we note that despite being 3 - 5 magnitudes brighter than DSCs, ultra-compact dwarfs have similar surface brightnesses. The closest Galactic analogs to the DSCs are the old open clusters. We suggest that if a diffuse star cluster population did exist in the disk of the Milky Way, it would be very difficult to find. (Abridged)

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.

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)

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.

Models for the synchrotron emission of gamma-ray burst afterglows suggest that the magnetic field is generated in the shock wave that forms as relativistic ejecta plow through the circum-burst medium. Transverse Weibel instability efficiently generates magnetic fields near equipartition with the post-shock energy density. The detailed saturated state of the instability, as seen in particle-in-cell simulations, consists of magnetically self-pinched current filaments. The filaments are parallel to the direction of propagation of the shock and are about a plasma skin depth in radius, forming a quasi--two-dimensional structure. We use a rudimentary analytical model to argue that the Weibel filaments are unstable to a kink-like mode, which destroys their quasi--two-dimensional structure. For wavelengths longer than than the skin depth, the instability grows at the rate equal to the speed of light divided by the wavelength. We calculate the transport of collisionless test particles in the filaments experiencing the instability and show that the particles diffuse in energy. This diffusion marks the beginning of thermalization in the shock transition layer, and causes initial magnetic field decay as particles escape from the filaments. We discuss the implications of these results for the structure of the shock and the polarization of the afterglow.

We have measured half-light radii, r_h, for globular clusters (GCs) belonging to the 100 early-type galaxies observed in the ACS Virgo Cluster Survey and the elliptical galaxy NGC 4697. An analysis of the dependencies of the measured r_h on both the properties of the GCs themselves and their host galaxies reveals that the average r_h increases with increasing galactocentric distance or, alternatively, with decreasing galaxy surface brightness. For the first time, we find that the average r_h decreases with the host galaxy color. We also show that there is no evidence for a variation of r_h with the luminosity of the GCs. Finally, we find in agreement with previous observations that the average r_h depends on the color of GCs, with red GCs being ~17% smaller than their blue counterparts. We show that this difference is probably a consequence of an intrinsic mechanism, rather than projection effects, and that it is in good agreement with the mechanism proposed in Jordan (2004). We discuss these findings in light of two simple pictures for the origin of the r_h of GCs and show that both lead to a behavior in rough agreement with the observations. After accounting for the dependencies found we show that the average GC half-light radii <r_h> can be successfully used as a standard ruler for distance estimation. We outline the methodology, and provide a calibration for its use. We find <r_h> = 2.7 +- 0.35 pc for GCs with (g-z)=1.2 mag in a galaxy with color (g-z)_{gal}=1.5 mag and at an underlying surface z-band brightness of mu_z = 21 mag arcsec^{-2}. Using this technique, we place an upper limit of 3.4 Mpc on the 1-sigma line-of-sight depth of the Virgo Cluster. Finally, we examine the form of the r_h distribution for our sample galaxies and provide an analytic expression which successfully describes this distribution. (Abridged)

It is generally accepted that magnetic fields generated in the nonlinear development of the transverse Weibel instability provide effective collisionality in unmagnetized collisionless shocks. Recently, extensive two and three dimensional simulations improved our understanding of the growth and saturation of the instability in colliding plasma shells. However, the steady-state structure of the shock wave transition layers remains poorly understood. We use basic physical considerations and order-of-magnitude arguments to study the steady state structure in relativistic unmagnetized collisionless shocks in pair plasmas. The shock contains an electrostatic layer resulting from the formation of stationary, magnetically-focused current filaments. The filaments form where the cold upstream plasma and the counterstreaming thermal plasma interpenetrate. The filaments are not entirely neutral and strong electrostatic fields are present. Most of the downstream particles cannot cross this layer into the upstream because they are trapped by the electrostatic field. We identify the critical location in the shock transition layer where the electromagnetic field ceases to be static. At this location, the degree of charge separation in the filaments reaches a maximum value, the current inside the filaments comes close to the Alfv\'en limit, and the phase space distribution function starts to isotropize. We argue that the radius of the current filaments upstream of the critical location is about twice the upstream plasma skin depth. Finally, we show that some downstream particles cross the electrostatic layer and run ahead of the shock into the preshock medium without causing instability. These particles may play an important role in particle acceleration.

As part of the Advanced Camera for Survey (ACS) Virgo Cluster Survey, we have measured Surface Brightness Fluctuations (SBF) in a sample of 100 early-type Virgo galaxies. Distances derived from these measurements are needed to explore the three-dimensional structure of the Virgo Cluster, study the intrinsic parameters of globular clusters associated with the program galaxies, and compare with the galaxy distances derived from globular cluster luminosity functions. Our SBF measurements have been performed in the F850LP bandpass of the Wide Field Channel of the ACS on the Hubble Space Telescope. These are the first measurements of this kind, and we present the first SBF calibration for this bandpass. The measured fluctuations depend on galaxy stellar population variations, which we quantify by galaxy color (g_475-z_850)_0, where g_475 and z_850 are the galaxy magnitudes, respectively, in the F475W and F850LP ACS filters. We derive the following calibration for the absolute SBF magnitude M_850: M_850 = -2.06 +/- 0.04 + (2.0 +/ 0.2) [ (g_475-z_850)_0-1.3 ] in the range 1.3 <(g_475-z_850)_0 <= 1.6, and M_850 = -2.06 +/- 0.04 + (0.9 +/- 0.2) [ (g_475-z_850)_0-1.3 ] in the range 1.0 <=(g_475-z_850)_0 <= 1.3. Physically, the two different color regimes correspond to different galaxy types: giant ellipticals and S0s at the red end, and early-type dwarfs at the blue end. For the first time in SBF studies, we are able to provide a firm empirical calibration of SBF in early-type dwarf galaxies. Our results agree with stellar population model predictions from Bruzual & Charlot (2003).

Using the empirical relations between the central galaxy luminosity and the halo mass, and between the total galaxy luminosity in a halo and the halo mass, we construct the galaxy luminosity function (LF). To the luminosity of the central galaxy in a halo of a given mass we assign log-normal scatter with a mean calibrated against the observations. In halos where the total galaxy luminosity exceeds that of the central galaxy, satellite galaxies are distributed as a power-law in luminosity. Combined with the halo mass function, this description reproduces the observed characteristics of the galaxy LF, including a shape consistent with the Schechter function. When all galaxies are included, regardless of the environment or the Hubble type, the Schechter L* is the luminosity scale above which the central galaxy luminosity-halo mass relation flattens; L* corresponds to ~10^{13}Msun on the halo mass scale. In surveys where central galaxies in massive clusters are neglected, either by design or because of the cosmic variance, L* is simply the mean luminosity of central galaxies in halos at the upper end of the selected mass range. The smooth, exponential decay of the Schechter function toward high luminosities reflects the intrinsic scatter in the central galaxy luminosity-halo mass relation. In addition to the LF, the model successfully reproduces the observed dependence of galaxy clustering bias on luminosity.

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.

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)

The final merger of a pair of massive black holes in a galactic nucleus is compelled by gravitational radiation. Gravitational waves from the mergers of black holes of masses (10^5-10^7)(1+z)^{-1} Msun at redshifts of 1-20 will be readily detectable by the Laser Interferometer Space Antenna (LISA), but an electromagnetic afterglow would be helpful in pinpointing the source and its redshift. Long before the merger, the binary "hollows out" any surrounding gas and shrinks slowly compared to the viscous timescale of a circumbinary disk. The inner gas disk is truncated at the radius where gravitational torque from the binary balances the viscous torque, and accretion onto the black holes is diminished. Initially, the inner truncation radius is able to follow the shrinking binary inward. But eventually the gravitational radiation timescale becomes shorter than the viscous timescale in the disk, leading to a merged black hole surrounded by a hollow disk of gas. We show that the subsequent viscous evolution of the hollow, radiation-pressure dominated disk will create a ~10^{43.5}(M/10^6Msun) ergs s^{-1} X-ray source on a timescale \~7(1+z)(M/10^6Msun)^{1.32} yr. This justifies follow-up monitoring of gravitational wave events with next-generation X-ray observatories. Analysis of the detailed light curve of these afterglows will yield new insights into the subtle physics of accretion onto massive black holes.

The Advanced Camera for Surveys (ACS) Virgo Cluster Survey is a large program to image 100 early-type Virgo galaxies using the F475W and F850LP bandpasses of the Wide Field Channel of the ACS instrument on the Hubble Space Telescope (HST). The scientific goals of this survey include an exploration of the three-dimensional structure of the Virgo Cluster and a critical examination of the usefulness of the globular cluster luminosity function as a distance indicator. Both of these issues require accurate distances for the full sample of 100 program galaxies. In this paper, we describe our data reduction procedures and examine the feasibility of accurate distance measurements using the method of surface brightness fluctuations (SBF) applied to the ACS Virgo Cluster Survey F850LP imaging. The ACS exhibits significant geometrical distortions due to its off-axis location in the HST focal plane; correcting for these distortions by resampling the pixel values onto an undistorted frame results in pixel correlations that depend on the nature of the interpolation kernel used for the resampling. This poses a major challenge for the SBF technique, which normally assumes a flat power spectrum for the noise. We investigate a number of different interpolation kernels and show through an analysis of simulated galaxy images having realistic noise properties that it is possible, depending on the kernel, to measure SBF distances using distortion-corrected ACS images without introducing significant additional error from the resampling. We conclude by showing examples of real image power spectra from our survey.

The ACS Virgo Cluster Survey is a large program to carry out multi-color imaging of 100 early-type members of the Virgo Cluster using the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. Deep F475W and F850LP images (~ SDSS g and z) are being used to study the central regions of the program galaxies, their globular cluster systems, and the three-dimensional structure of Virgo itself. In this paper, we describe in detail the data reduction procedures used for the survey, including image registration, drizzling strategies, the computation of weight images, object detection, the identification of globular cluster candidates, and the measurement of their photometric and structural parameters.

The ACIS instrument on board the Chandra X-ray Observatory has been used to carry out the first systematic study of low-mass X-ray binaries (LMXBs) in M87. We identify 174 X-ray point-sources, of which ~150 are likely LMXBs. This LMXB catalog is combined with deep F475W and F850LP images taken with ACS on HST to examine the connection between LMXBs and globular clusters in M87. Of the 1688 globular clusters in our catalog, f_X = 3.6 +- 0.5% contain a LMXB and we find that the metal-rich clusters are 3 +- 1 times more likely to harbor a LMXB than their metal-poor counterparts. In agreement with previous findings for other galaxies, we find that brighter, more metal-rich clusters are more likely to contain a LMXB. For the first time, however, we are able to demonstrate that the probability, p_X, that a given cluster will contain a LMXB depends sensitively on the dynamical properties of the host cluster. Specifically, we use the HST images to measure the half-light radius, concentration index and central density, \rho_0, for each globular, and define a parameter, \Gamma, which is related to the tidal capture and binary-neutron star exchange rate. Our preferred form for p_X is then p_X \propto \Gamma \rho_0^{-0.42\pm0.11} (Z/Z_{\odot})^{0.33\pm0.1}. We argue that if the form of p_X is determined by dynamical processes, then the observed metallicity dependence is a consequence of an increased number of neutron stars per unit mass in metal-rich globular clusters. Finally, we find no compelling evidence for a break in the luminosity distribution of resolved X-ray point sources. Instead, the LMXB luminosity function is well described by a power law with an upper cutoff at L_X ~ 10^39 erg/s. (abridged)

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.

The Virgo Cluster is the dominant mass concentration in the Local Supercluster and the largest collection of elliptical and lenticular galaxies in the nearby universe. In this paper, we present an introduction to the ACS Virgo Cluster Survey: a program to image, in the F475W and F850LP bandpasses (Sloan g and z), 100 early-type galaxies in the Virgo Cluster using the Advanced Camera for Surveys on the Hubble Space Telescope. We describe the selection of the program galaxies and their ensemble properties, the choice of filters, the field placement and orientation, the limiting magnitudes of the survey, coordinated parallel observations of 100 "intergalactic" fields with WFPC2, and supporting ground-based spectroscopic observations of the program galaxies. In terms of depth, spatial resolution, sample size and homogeneity, this represents the most comprehensive imaging survey to date of early-type galaxies in a cluster environment. We briefly describe the main scientific goals of the survey which include the measurement of luminosities, metallicities, ages, and structural parameters for the many thousands of globular clusters associated with these galaxies, a high-resolution isophotal analysis of galaxies spanning a factor of ~ 450 in luminosity and sharing a common environment, the measurement of accurate distances for the full sample of galaxies using the method of surface brightness fluctuations, and a determination of the three-dimensional structure of Virgo itself.

A large fraction of bulgeless disk galaxies contain young compact stellar systems at their centers, in spite of the local gravitational stability of these disks. We evaluate two contrasting hypotheses for the origin of the nuclear star clusters in late-type disk galaxies. The clusters could not have migrated from distant eccentric locations in the disk. Instead they must have formed in situ, requiring radial transport of gas toward the center of the disk. This transport could be a consequence of the development of the magnetorotational instability in the differentially rotating warm neutral medium. We evaluate the rate of gas transport into the disk center and find that it is sufficient to support continuous star formation in that location. Enhanced stellar surface brightness in the inner few hundred parsecs and the formation of a compact stellar system in the central few parsecs are unavoidable in dark matter halos with divergent density profiles. We illustrate our conclusions on a model of the nearest late-type disk galaxy M33.

The discovery of hundreds of young, bright stars within a parsec from the massive black hole at the center of the Galaxy presents a challenge to star formation theories. The requisite Roche densities for the gravitational collapse of gas clouds are most naturally achieved in accretion disks. The water maser sources in Keplerian rotation in the nuclei of NGC4258, NGC1068, and the Circinus galaxy indicate the presence of warm, extended, molecular accretion disks around black holes similar in mass to the one at the Galactic center. Here we argue that the current conditions in the maser nuclei, and those near the Galactic center, represent two consecutive, recurrent phases in the life cycle of the nucleus of a typical gas-rich spiral bulge. The warm molecular disks that give rise to the observed maser emission fragment into stellar-size objects. The stellar masses, their orbital geometry, and the total number of stars thus formed are consistent with the values identified at the Galactic center. The stars tend to form in compact groups resembling the IRS 13 complex that dominates the stellar light in the neighborhood of the black hole.

The long-term evolution of massive black hole binaries at the centers of galaxies is studied in a variety of physical regimes, with the aim of resolving the ``final parsec problem,'' i.e., how black hole binaries manage to shrink to separations at which emission of gravity waves becomes efficient. A binary ejects stars by the gravitational slingshot and carves out a loss cone in the host galaxy. Continued decay of the binary requires a refilling of the loss cone. We show that the standard treatment of loss cone refilling, derived for collisionally relaxed systems like globular clusters, can substantially underestimate the refilling rates in galactic nuclei. We derive expressions for non-equilibrium loss-cone dynamics and calculate time scales for the decay of massive black hole binaries following galaxy mergers, obtaining significantly higher decay rates than heretofore. Even in the absence of two-body relaxation, decay of binaries can persist due to repeated ejection of stars returning to the nucleus on eccentric orbits. We show that this recycling of stars leads to a gradual, approximately logarithmic dependence of the binary binding energy on time. We derive an expression for the loss cone refilling induced by the Brownian motion of a black hole binary. We also show that numerical N-body experiments are not well suited to probe these mechanisms over long times due to spurious relaxation.

The coalescence of massive black hole binaries is one of the main sources of low-frequency gravitational radiation that can be detected by LISA. When two galaxies containing massive black holes merge, a binary forms at the center of the new galaxy. We discuss the evolution of the binary after its separation decreases below one parsec. Whether or not stellar dynamical processes can drive the black holes to coalesce depends on the supply of stars that scatter against the binary. We discuss various mechanisms by which this supply can be replenished after the loss cone has been depleted.

Deep infra-red observations and long-term monitoring programs have provided dynamical evidence for a supermassive black hole of mass 3.e6 solar masses associated with the radio source Sagitarrius A* at the center of our Galaxy. The brightest stars orbiting within 0.1 parsecs of the black hole appear to be young, massive main sequence stars, n spite of an environment near the black hole that is hostile to star formation. We discuss mechanisms by which stars born outside the central parsec can sink towards the black hole and conclude that the drag coming from plausible stellar populations does not operate on the short timescales required by the stellar ages. We propose that these stars were dragged in by a second black hole of mass of 1.e3-1.e4 solar masses, which would be classified as an intermediate-mass black hole. We discuss the implications for the stellar populations and the kinematics in the Galactic center. Finally we note that continued astrometric monitoring of the central radio source offers the prospect for a direct detection of such objects.

We revisit the issue of cosmological parameter estimation in light of current and upcoming high-precision measurements of the cosmic microwave background power spectrum. Physical quantities which determine the power spectrum are reviewed, and their connection to familiar cosmological parameters is explicated. We present a set of physical parameters, analytic functions of the usual cosmological parameters, upon which the microwave background power spectrum depends linearly (or with some other simple dependence) over a wide range of parameter values. With such a set of parameters, microwave background power spectra can be estimated with high accuracy and negligible computational effort, vastly increasing the efficiency of cosmological parameter error determination. The techniques presented here allow calculation of microwave background power spectra $10^5$ times faster than comparably accurate direct codes (after precomputing a handful of power spectra). We discuss various issues of parameter estimation, including parameter degeneracies, numerical precision, mapping between physical and cosmological parameters, and systematic errors, and illustrate these considerations with an idealized model of the MAP experiment.

Formation and disruption of dark-matter cusps are reviewed. Accumulation of baryons at the center of a halo can displace the dark matter, converting singular density cusps into low-density cores. The displaced mass can be of order 10 M_b with M_b the mass of the infalling population. If M_b is identified with the masses of the black holes currently observed at the centers of bright galaxies, predicted core radii are ~a few hundred pc. Other mechanisms, such as early mass outflow, may explain the large dark-matter cores in dwarf and low-surface-brightness galaxies. Predictions of dark matter annihilation radiation from the center of the Milky Way galaxy are shown to be strongly dependent on the galaxy's merger history.

The annihilation rate of weakly interacting cold dark matter particles at the galactic center could be greatly enhanced by the growth of a density spike around the central supermassive black hole (SBH). Here we discuss the effects of hierarchical mergers on the central spike. Mergers between halos containing SBHs lead to the formation of SBH binaries which transfer energy to the dark matter particles, lowering their density. The predicted flux of annihiliation photons from the galactic center is several orders of magnitude smaller than in models that ignore the effects of SBHs and mergers. Measurement of the annihilation radiation could in principle be used to constrain the merger history of the galaxy.

We investigate the hypothesis that the cores of elliptical galaxies and bulges are created from the binding energy liberated by the coalescence of supermassive binary black holes during galaxy mergers. Assuming that the central density profiles of galaxies were initially steep power laws, we define the ``mass deficit'' as the mass in stars that had to be removed from the nucleus in order to produce the observed core. We use nonparametric deprojection to compute the mass deficit in a sample of 35 early-type galaxies with high-resolution imaging data. We find that the mass deficit correlates well with the mass of the nuclear black hole, consistent with the predictions of merger models. We argue that core sizes in halos of non-interacting dark matter particles should be comparable to those observed in the stars.

We investigate a model in which galactic nuclei form via the coalescence of pre-existing stellar systems containing supermassive black holes. Merger simulations are carried out using N-body algorithms that can follow the formation and decay of a black-hole binary and its effect on the surrounding stars down to sub-parsec scales. Our initial stellar systems have steep central density cusps similar to those in low-luminosity elliptical galaxies. Formation of a black-hole binary transfers energy to the stars and lowers the central density; continued decay of the binary creates a ~1/r density cusp similar to those observed in bright elliptical galaxies, with a break radius that extends well beyond the sphere of gravitational influence of the black holes. The decay of the black hole binary is followed over a factor of ~20 in separation after formation of a hard binary, considerably farther than in previous simulations. We see almost no dependence of the binary's decay rate on number of particles in the simulation, contrary to earlier studies in which a lower initial density of stars led to a more rapid depletion of the binary's loss cone. We nevertheless argue that the decay of a black hole binary in a real galaxy would be expected to stall at separations of 0.01-1 pc unless some additional mechanism is able to extract energy from the binary. Our results support a picture in which the observed dependence of nuclear cusp slope on galaxy luminosity is a consequence of galaxy interactions. We also discuss the implications of our results for the survivability of dark-matter cusps.

Space-based observations are beginning to yield detailed information about the stellar kinematics at the very centers of galaxies, within the sphere of gravitational influence of the black hole. The structure and dynamics of these regions is probably determined in part by the infall and coalescence of black holes during galaxy mergers. A goal of N-body simulations is to reproduce the kinematics near the black holes as well as the relations that exist between the nuclear and global properties of galaxies. However, the problem is computationally difficult due to the wide range of length and time scales, and no single N-body code can efficiently follow the evolution from kiloparsec to sub-parsec scales. We review existing N-body work on this problem and present the first, fully self-consistent merger simulations of galaxies containing dense stellar cusps and black holes.