The U-band Galaxy Luminosity Function of Nearby Clusters

Despite the great potential of the U-band galaxy luminosity function (GLF) to constrain the history of star formation in clusters, to clarify the question of variations of the GLF across filter bands, to provide a baseline for comparisons to high-redshift studies of the cluster GLF, and to estimate the contribution of bound systems of galaxies to the extragalactic near-UV background, determinations have so far been hampered by the generally low efficiency of detectors in the U-band and by the difficulty of constructing both deep and wide surveys. In this paper, we present U-band GLFs of three nearby, rich clusters to a limit of M_U=-17.5 (M*_U+2). Our analysis is based on a combination of separate spectroscopic and R-band and U-band photometric surveys. For this purpose, we have developed a new maximum-likelihood algorithm for calculating the luminosity function that is particularly useful for reconstructing the galaxy distribution function in multi-dimensional spaces (e.g., the number of galaxies as a simultaneous function of luminosity in different filter bands, surface brightness, star formation rate, morphology, etc.), because it requires no prior assumptions as to the shape of the distribution function. The composite luminosity function can be described by a Schechter function with characteristic magnitude M*_U=-19.82+/-0.27 and faint end slope alpha_U=-1.09+/-0.18. The total U-band GLF is slightly steeper than the R-band GLF, indicating that cluster galaxies are bluer at fainter magnitudes. Quiescent galaxies dominate the cumulative U-band flux for M_U<-14. The contribution of galaxies in nearby clusters to the U-band extragalactic background is <1% Gyr^-1 for clusters of masses ~3*10^14 to 2*10^15 M_solar.Comment: 44 pages, 11 figures, accepted for publication in Ap

Inference on co-integration parameters in heteroskedastic vector autoregressions

We consider estimation and hypothesis testing on the coefficients of the co-integrating relations and the adjustment coefficients in vector autoregressions driven by shocks which display both conditional and unconditional heteroskedasticity of a quite general and unknown form. We show that the conventional results in Johansen (1996) for the maximum likelihood estimators and associated likelihood ratio tests derived under homoskedasticity do not in general hold under heteroskedasticity. As a result, standard confidence intervals and hypothesis tests on these coefficients are potentially unreliable. Solutions based on Wald tests (using a “sandwich” estimator of the variance matrix) and on the use of the wild bootstrap are discussed. These do not require the practitioner to specify a parametric model for volatility. We establish the conditions under which these methods are asymptotically valid. A Monte Carlo simulation study demonstrates that significant improvements in finite sample size can be obtained by the bootstrap over the corresponding asymptotic tests in both heteroskedastic and homoskedastic environments. An application to the term structure of interest rates in the US illustrates the difference between standard and bootstrap inferences regarding hypotheses on the co-integrating vectors and adjustment coefficients

Large Eddy Simulation/Conditional Moment Closure modeling of swirl-stabilized non-premixed flames with local extinction

This is the accepted manuscript. The final version is available from Elsevier at http://www.sciencedirect.com/science/article/pii/S1540748914000558.The Large Eddy Simulation (LES)/three-dimensional Conditional Moment Closure (3D-CMC) model with detailed chemistry and finite-volume formulation is employed to simulate a swirl-stabilized nonpremixed flame with local extinction. The results demonstrate generally good agreement with the measurements concerning velocity, flame shape, and statistics of flame lift-off, but the penetration of fuel jet into the recirculation zone is under-predicted possibly due to the over-predicted swirl velocities in the chamber. Localized extinctions are seen in the LES, in agreement with experiment. The local extinction event is shown by very low heat release rate and hydroxyl mass fraction and reduced temperature, and is accompanied by relatively high scalar dissipation. In mixture fraction space, CMC cells with strong turbulence-chemistry interaction and local extinction show relatively large fluctuations between fully burning and intermediate distributions. The probability density functions of conditional reactedness, which shows how far the conditionally-filtered scalars are from reference fully burning profiles, indicate that for CMC cells with local extinction, some reactive scalars demonstrate pronounced bimodality while for those cells with strong reactivity the PDFs are very narrow.HZ acknowledges the financial support from EPSRC through a Dorothy Hodgkin Postgraduate Award

Images, structural properties and metal abundances of galaxy clusters observed with Chandra ACIS-I at 0.1<z<1.3

We have assembled a sample of 115 galaxy clusters at 0.1<z<1.3 with archived Chandra ACIS-I observations. We present X-ray images of the clusters and make available region files containing contours of the smoothed X-ray emission. The structural properties of the clusters were investigated and we found a significant absence of relaxed clusters (as determined by centroid shift measurements) at z>0.5. The slope of the surface brightness profiles at large radii were steeper on average by 15% than the slope obtained by fitting a simple beta-model to the emission. This slope was also found to be correlated with cluster temperature, with some indication that the correlation is weaker for the clusters at z>0.5. We measured the mean metal abundance of the cluster gas as a function of redshift and found significant evolution, with the abundances dropping by 50% between z=0.1 and z~1. This evolution was still present (although less significant) when the cluster cores were excluded from the abundance measurements, indicating that the evolution is not solely due to the disappearance of relaxed, cool core clusters (which are known to have enhanced core metal abundances) from the population at z>0.5.Comment: 23 pages, 12 figures. Accepted for publication in ApJS. Updated to match published version. Redshifts of two clusters (RXJ1701 and CL0848) corrected and two observations of MACSJ0744.8 have been combined into one. Conclusions unchanged. A version with images of all of the clusters is available at http://hea-www.harvard.edu/~bmaughan/clusters.htm

Dark Energy and the mass of galaxy clusters

Up to now, Dark Energy evidences are based on the dynamics of the universe on very large scales, above 1 Gpc. Assuming it continues to behave like a cosmological constant $\Lambda$ on much smaller scales, I discuss its effects on the motion of non-relativistic test-particles in a weak gravitational field and I propose a way to detect evidences of $\Lambda \neq 0$ at the scale of about 1 Mpc: the main ingredient is the measurement of galaxy cluster masses.Comment: 5 pages, no figures, references adde

Self-Consistent Theory of Halo Mergers

The rate of merging of dark-matter halos is an absolutely essential ingredient for studies of both structure and galaxy formation. Remarkably, however, our quantitative understanding of the halo merger rate is still quite limited, and current analytic descriptions based upon the extended Press-Schechter formalism are fundamentally flawed. We show that a mathematically self-consistent merger rate must be consistent with the evolution of the halo abundance in the following sense: The merger rate must, when inserted into the Smoluchowski coagulation equation, yield the correct evolution of the halo abundance. We then describe a numerical technique to find merger rates that are consistent with this evolution. We present results from a preliminary study in which we find merger rates that reproduce the evolution of the halo abundance according to Press-Schechter for power-law power spectra. We discuss the limitations of the current approach and outline the questions that must still be answered before we have a fully consistent and correct theory of halo merger rates.Comment: 13 pages, 8 figures, submitted to MNRAS. Version with full resolution figures available at http://www-astro.physics.ox.ac.uk/~abenson/Papers/smoluchow.pd

Negative differential conductance in quantum dots in theory and experiment

Experimental results for sequential transport through a lateral quantum dot in the regime of spin blockade induced by spin dependent tunneling are compared with theoretical results obtained by solving a master equation for independent electrons. Orbital and spin effects in electron tunneling in the presence of a perpendicular magnetic field are identified and discussed in terms of the Fock-Darwin spectrum with spin. In the nonlinear regime, a regular pattern of negative differential conductances is observed. Electrical asymmetries in tunnel rates and capacitances must be introduced in order to account for the experimental findings. Fast relaxation of the excited states in the quantum dot have to be assumed, in order to explain the absence of certain structures in the transport spectra.Comment: 4 pages, 4 figure

X-ray Over-Luminous Elliptical Galaxies: A New Class of Mass Concentrations in the Universe?

We detect four isolated, X-ray over-luminous (Lx>2e43 [h/0.5]**-2 erg/s) elliptical galaxies (OLEGs) in our 160 square degree ROSAT PSPC survey. The extent of their X-ray emission, total X-ray luminosity, total mass, and mass of the hot gas in these systems correspond to poor clusters, and the optical luminosity of the central galaxies (M_R<-22.5 + 5 lg h) is comparable to that of cluster cDs. However, there are no detectable fainter galaxy concentrations around the central elliptical. The mass-to-light ratio within the radius of detectable X-ray emission is in the range 250-450 Msun/Lsun, which is 2-3 times higher than typically found in clusters or groups. These objects can be the result of galaxy merging within a group. However, their high M/L values are difficult to explain in this scenario. OLEGs must have been undisturbed for a very long time, which makes them the ultimate examples of systmes in hydrostatic equilibrium. The number density of OLEGs is n=2.4(+3.1-1.2}x10**-7 (h/0.5)**-3 Mpc**-3 at the 90% confidence. They comprise 20% of all clusters and groups of comparable X-ray luminosity, and nearly all galaxies brighter than M_R=-22.5. The estimated contirubution of OLEGs to the total mass density in the Universe is close to that of T>7 keV clusters.Comment: 4 pages, 2 figures, uses emulateapj.sty, submitted to ApJ Letter

Massive Coronae of Galaxies

There is reason to suspect that about half of the baryons are in pressure-supported plasma in the halos of normal galaxies, drawn in by gravity along with about half of the dark matter. To be consistent with the observations this baryonic component, the corona, would have to be hotter than the kinetic temperature of the dark matter in the halo so as to produce acceptable central electron densities. We ascribe this hotter plasma temperature to the addition of entropy prior to and during assembly of the system, in an analogy to cluster formation. The plasma cooling time would be longer than the gravitational collapse time but, in the inner parts, shorter than the Hubble time, making the corona thermally unstable to the formation of a cloudy structure that may be in line with what is indicated by quasar absorption line systems. The corona of an isolated spiral galaxy would be a source of soft X-ray and recombination radiation, adding to the more commonly discussed effects of stars and supernovae. In this picture the mass in the corona is much larger than the mass in condensed baryons in a spiral galaxy. The corona thus would be a substantial reservoir of diffuse baryons that are settling and adding to the mass in interstellar matter and stars, so that star formation in isolated spirals will continue well beyond the present epoch.Comment: 14 pages, 5 figure

CIRS: Cluster Infall Regions in the Sloan Digital Sky Survey I. Infall Patterns and Mass Profiles

We use the Fourth Data Release of the Sloan Digital Sky Survey to test the ubiquity of infall patterns around galaxy clusters and measure cluster mass profiles to large radii. We match X-ray cluster catalogs with SDSS, search for infall patterns, and compute mass profiles for a complete sample of X-ray selected clusters. Very clean infall patterns are apparent in most of the clusters, with the fraction decreasing with increasing redshift due to shallower sampling. All 72 clusters in a well-defined sample limited by redshift (ensuring good sampling) and X-ray flux (excluding superpositions) show infall patterns sufficient to apply the caustic technique. This sample is by far the largest sample of cluster mass profiles extending to large radii to date. Similar to CAIRNS, cluster infall patterns are better defined in observations than in simulations. Further work is needed to determine the source of this difference. We use the infall patterns to compute mass profiles for 72 clusters and compare them to model profiles. Cluster scaling relations using caustic masses agree well with those using X-ray or virial mass estimates, confirming the reliability of the caustic technique. We confirm the conclusion of CAIRNS that cluster infall regions are well fit by NFW and Hernquist profiles and poorly fit by singular isothermal spheres. This much larger sample enables new comparisons of cluster properties with those in simulations. The shapes (specifically, NFW concentrations) of the mass profiles agree well with the predictions of simulations. The mass inside the turnaround radius is on average 2.19$\pm$0.18 times that within the virial radius. This ratio agrees well with recent predictions from simulations of the final masses of dark matter haloes.Comment: 34 pages, 24 figures, accepted for publication in AJ, full resolution version available at http://www.astro.yale.edu/krines