Search

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., Comment: 20 pages, 16 figures

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., Comment: 7 pages, 3 figures. Accepted for publication in MNRAS

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., Comment: 18 pages, 7 figures, MNRAS, in press

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], Comment: 25 pages (including 5 pages appendix), 13 figures. Accepted for publication in ApJ

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., Comment: 13 pages, 8 figures

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., Comment: 21 pages, 19 figures, submitted to ApJ on 4 Nov 2008

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., Comment: ApJ, accepted (26 pages, 16 figures, 7 tables)

(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., Accepted by The Astrophysical Journal Supplement Series. This submission contains low resolution figures; we strongly recommend downloading the original version of the paper from the ACSVCS project website: http://www.cadc.hia.nrc.gc.ca/community/ACSVCS/publications.html#acsvcs6

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., Comment: 6 pages, 1 figure

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., Comment: 6 pages, 1 figure

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), 25 pages, 16 figures, 7 tables. Accepted for publication in the Astrophysical Journal. Version with Figures 1,5-9,11-16 in color is available at http://www.physics.rutgers.edu/~pcote/acs/acsvcs_9.pdf

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., Comment: 4 pages, 2 figures, Astrophys. J. Lett., in press

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., 14 pages, 5 figures; discussion no longer aperture-specific (TMT -> ELT), matches ApJ version

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., Comment: 4 pages, ApJ Letters, in press

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., Comment: 4 pages, 3 postscript figures, uses revtex4.sty

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., Comment: 35 pages, 19 postscript figures, uses emulateapj.sty, onecolfloat.sty