Globular Cluster Systems in Brightest Cluster Galaxies. III: Beyond Bimodality

We present new deep photometry of the rich globular cluster (GC) systems around the Brightest Cluster Galaxies UGC 9799 (Abell 2052) and UGC 10143 (Abell 2147), obtained with the HST ACS and WFC3 cameras. For comparison, we also present new reductions of similar HST/ACS data for the Coma supergiants NGC 4874 and 4889. All four of these galaxies have huge cluster populations (to the radial limits of our data, comprising from 12000 to 23000 clusters per galaxy). The metallicity distribution functions (MDFs) of the GCs can still be matched by a bimodal-Gaussian form where the metal-rich and metal-poor modes are separated by ~0.8 dex, but the internal dispersions of each mode are so large that the total MDF becomes very broad and nearly continuous from [Fe/H] = -2.4 to Solar. There are, however, significant differences between galaxies in the relative numbers of \emph{metal-rich} clusters, suggesting that they underwent significantly different histories of mergers with massive, gas-rich halos. Lastly, the proportion of metal-poor GCs rises especially rapidly outside projected radii R > 4 R_eff, suggesting the importance of accreted dwarf satellites in the outer halo. Comprehensive models for the formation of GCs as part of the hierarchical formation of their parent galaxies will be needed to trace the systematic change in structure of the MDF with galaxy mass, from the distinctly bimodal form in smaller galaxies up to the broad continuum that we see in the very largest systems.Comment: In press for Astrophysical Journa

Cosmological Reionization Around the First Stars: Monte Carlo Radiative Transfer

We study the evolution of ionization fronts around the first proto-galaxies by using high resolution numerical cosmological (Lambda+CDM model) simulations and Monte Carlo radiative transfer methods. We present the numerical scheme in detail and show the results of test runs from which we conclude that the scheme is both fast and accurate. As an example of interesting cosmological application, we study the reionization produced by a stellar source of total mass M=2 10^8 M_\odot turning on at z=12, located at a node of the cosmic web. The study includes a Spectral Energy Distribution of a zero-metallicity stellar population, and two Initial Mass Functions (Salpeter/Larson). The expansion of the I-front is followed as it breaks out from the galaxy and it is channeled by the filaments into the voids, assuming, in a 2D representation, a characteristic butterfly shape. The ionization evolution is very well tracked by our scheme, as realized by the correct treatment of the channeling and shadowing effects due to overdensities. We confirm previous claims that both the shape of the IMF and the ionizing power metallicity dependence are important to correctly determine the reionization of the universe.Comment: 8 pages, 8 figures. Revised version, accepted for publication by MNRA

On the Globular Cluster IMF below 1 Solar Mass

(Abridged) Accurate luminosity functions (LF) for a dozen globular clusters have now been measured at or just beyond their half-light radius using HST. They span almost the entire cluster main sequence below ~ 0.75 Msolar. All these clusters exhibit LF that rise continuously from an absolute I magnitude M_I ~ 6 to a peak at M_I ~ 8.5-9 and then drop with increasing M_I. Transformation of the LF into mass functions (MF) by means of the most recent mass luminosity relations that are consistent with all presently available data on the physical properties of low mass, low metallicity stars shows that all the LF observed so far can be obtained from MF having the shape of a log-normal distribution with characteristic mass m_c=0.33 +/- 0.03 Msolar and standard deviation sigma = 1.81 +/- 0.19. After correction for the effects of mass segregation, the variation of the ratio of the number of higher to lower mass stars with cluster mass or any simple orbital parameter or the expected time to disruption recently computed for these clusters shows no statistically significant trend over a range of this last parameter of more than a factor of 100. We conclude that the global MF of these clusters have not been measurably modified by evaporation and tidal interactions with the Galaxy and, thus, should reflect the initial distribution of stellar masses. Since the log-normal function that we find is also very similar to the one obtained independently for much younger clusters and to the form expected theoretically, the implication seems to be unavoidable that it represents the true stellar IMF for this type of stars in this mass range.Comment: Accepted for publication in The Astrophysical Journal. Contains 28 pages with 6 figure

Constraints on the mass and abundance of black holes in the Galactic halo: the high mass limit

We establish constraints on the mass and abundance of black holes in the Galactic halo by determining their impact on globular clusters which are conventionally considered to be little evolved. Using detailed Monte Carlo simulations and simple analytic estimates, we conclude that, at Galactocentric radius R~8 kpc, black holes with masses M_bh >~(1-3) x 10^6 M_sun can comprise no more than a fraction f_bh ~ 0.025-0.05 of the total halo density. This constraint significantly improves those based on disk heating and dynamical friction arguments as well as current lensing results. At smaller radius, the constraint on f_bh strengthens, while, at larger radius, an increased fraction of black holes is allowed.Comment: 13 pages, 10 figures, revised version, in press, Monthly Notice

Cosmological Magnetogenesis driven by Radiation Pressure

The origin of large scale cosmological magnetic fields remains a mystery, despite the continuous efforts devoted to that problem. We present a new model of magnetic field generation, based on local charge separation provided by an anisotropic and inhomogeneous radiation pressure. In the cosmological context, the processes we explore take place at the epoch of the reionisation of the Universe. Under simple assumptions, we obtain results (i) in terms of the order of magnitude of the field generated at large scales and (ii) in terms of its power spectrum. The amplitudes obtained (B ~ 8.10^(-6) micro-Gauss) are considerably higher than those obtained in usual magnetogenesis models and provide suitable seeds for amplification by adiabatic collapse and/or dynamo during structure formation.Comment: 9 pages, 2 figure

Scaling Relations of Dwarf Galaxies without Supernova-Driven Winds

Nearby dwarf galaxies exhibit tight correlations between their global stellar and dynamical properties, such as circular velocity, mass-to-light ratio, stellar mass, surface brightness, and metallicity. Such correlations have often been attributed to gas or metal-rich outflows driven by supernova energy feedback to the interstellar medium. We use high-resolution cosmological simulations of high-redshift galaxies with and without energy feedback, as well as analytic modeling, to investigate whether the observed correlations can arise without supernova-driven outflows. We find that the simulated dwarf galaxies exhibit correlations similar to those observed as early as z~10, regardless of whether supernova feedback is included. We also show that the correlations can be well reproduced by our analytic model that accounts for realistic gas inflow but assumes no outflows, and star formation rate obeying the Kennicutt-Schmidt law with a critical density threshold. We argue that correlations in simulated galaxies arise due to the increasingly inefficient conversion of gas into stars in low-mass dwarf galaxies rather than supernova-driven outflows. We also show that the decrease of the observed effective yield in low-mass objects, often used as an indicator of gas and metal outflows, can be reasonably reproduced in our simulations without outflows. We show that this trend can arise if a significant fraction of metals in small galaxies is spread to the outer regions of the halo outside the stellar extent via mixing. In this case the effective yield can be significantly underestimated if only metals within the stellar radius are taken into account. Measurements of gas metallicity in the outskirts of gaseous disks of dwarfs would thus provide a key test of such explanation.Comment: accepted for publication in ApJ, 19 pages, 12 figures, uses emulateapj

Hydrogen Phases on the Surface of a Strongly Magnetized Neutron Star

The outermost layers of some neutron stars are likely to be dominated by hydrogen, as a result of fast gravitational settling of heavier elements. These layers directly mediate thermal radiation from the stars, and determine the characteristics of X-ray/EUV spectra. For a neutron star with surface temperature T\lo 10^6 K and magnetic field B\go 10^{12} G, various forms of hydrogen can be present in the envelope, including atom, poly-molecules, and condensed metal. We study the physical properties of different hydrogen phases on the surface of a strongly magnetized neutron star for a wide range of field strength $B$ and surface temperature $T$. Depending on the values of $B$ and $T$, the outer envelope can be either in a nondegenerate gaseous phase or in a degenerate metallic phase. For T\go 10^5 K and moderately strong magnetic field, B\lo 10^{13} G, the envelope is nondegenerate and the surface material gradually transforms into a degenerate Coulomb plasma as density increases. For higher field strength, $B>> 10^{13}$ G, there exists a first-order phase transition from the nondegenerate gaseous phase to the condensed metallic phase. The column density of saturated vapor above the metallic hydrogen decreases rapidly as the magnetic field increases or/and temperature decreases. Thus the thermal radiation can directly emerge from the degenerate metallic hydrogen surface. The characteristics of surface X-ray/EUV emission for different phases are discussed. A separate study concerning the possibility of magnetic field induced nuclear fusion of hydrogen on the neutron star surface is also presented.Comment: TeX, 35 pages including 6 postscript figures. To be published in Ap

Fast Large Volume Simulations of the 21 cm Signal from the Reionization and pre-Reionization Epochs

While limited to low spatial resolution, the next generation low-frequency radio interferometers that target 21 cm observations during the era of reionization and prior will have instantaneous fields-of-view that are many tens of square degrees on the sky. Predictions related to various statistical measurements of the 21 cm brightness temperature must then be pursued with numerical simulations of reionization with correspondingly large volume box sizes, of order 1000 Mpc on one side. We pursue a semi-numerical scheme to simulate the 21 cm signal during and prior to Reionization by extending a hybrid approach where simulations are performed by first laying down the linear dark matter density field, accounting for the non-linear evolution of the density field based on second-order linear perturbation theory as specified by the Zel'dovich approximation, and then specifying the location and mass of collapsed dark matter halos using the excursion-set formalism. The location of ionizing sources and the time evolving distribution of ionization field is also specified using an excursion-set algorithm. We account for the brightness temperature evolution through the coupling between spin and gas temperature due to collisions, radiative coupling in the presence of Lyman-alpha photons and heating of the intergalactic medium, such as due to a background of X-ray photons. The hybrid simulation method we present is capable of producing the required large volume simulations with adequate resolution in a reasonable time so a large number of realizations can be obtained with variations in assumptions related to astrophysics and background cosmology that govern the 21 cm signal.Comment: 14 pages and 15 figures. New version to match accepted version for MNRAS. Code available in: http://www.SimFast21.or

The binary fraction in the globular cluster M10 (NGC 6254): comparing core and outer regions

We study the binary fraction of the globular cluster M10 (NGC 6254) as a function of radius from the cluster core to the outskirts, by means of a quan- titative analysis of the color distribution of stars relative to the fiducial main sequence. By taking advantage of two data-sets, acquired with the Advanced Camera for Survey and the Wide Field Planetary Camera 2 on board the Hubble Space Telescope, we have studied both the core and the external regions of the cluster. The binary fraction is found to decrease from 14% within the core, to 1.5% in a region between 1 and 2 half-mass radii from the cluster centre. Such a trend and the derived values are in agreement with previous results ob- tained in clusters of comparable total magnitude. The estimated binary fraction is sufficient to account for the suppression of mass segregation observed in M10, without any need to invoke the presence of an intermediate-mass black hole in its centre.Comment: Accepted for publication in ApJ (22 pages, 7 figures, 3 tables