Optical and X-ray Observations of the Afterglow to XRF030723

The X-ray-flash XRF030723 was detected by the HETE satellite and rapidly disseminated, allowing for an optical transient to be detected ~1 day after the burst. We discuss observations in the optical with Magellan, which confirmed the fade of the optical transient. In a 2-epoch ToO observation with Chandra, we discovered a fading X-ray source spatially coincident with the optical transient. We present spectral fits to the X-ray data. We also discuss the possibility that the source underwent a rebrightening in the X-rays, as was observed in the optical. We find that the significance of a possible rebrightening is very low (~1 sigma).Comment: 4 pages, 2 figures, to appear in Santa Fe GRB Conference Proceedings, 200

Observation and implications of the Epeak - Eiso correlation in Gamma-Ray Bursts

The availability of a few dozen GRB redshifts now allows studies of the intrinsic properties of these high energy transients. Amati et al. recently discovered a correlation between Epeak, the intrinsic peak energy of the $\nu f \nu$ spectrum, and Eiso, the isotropic equivalent energy radiated by the source. Lamb et al. have shown that HETE-2 data confirm and extend this correlation. We discuss here one of the consequences of this correlation: the existence of a 'spectral standard candle', which can be used to construct a simple redshift indicator for GRBs.Comment: Proceedings of the GRB 2003 Conference in SantaFe, 5 pages, 4 figure

The Halo Mass Function: High-Redshift Evolution and Universality

We study the formation of dark matter halos in the concordance LCDM model over a wide range of redshifts, from z=20 to the present. Our primary focus is the halo mass function, a key probe of cosmology. By performing a large suite of nested-box N-body simulations with careful convergence and error controls (60 simulations with box sizes from 4 to 256 Mpc/h, we determine the mass function and its evolution with excellent statistical and systematic errors, reaching a few percent over most of the considered redshift and mass range. Across the studied redshifts, the halo mass is probed over 6 orders of magnitude (10^7 - 10^13.5 M_sun/h). Historically, there has been considerable variation in the high redshift mass function as obtained by different groups. We have made a concerted effort to identify and correct possible systematic errors in computing the mass function at high redshift and to explain the discrepancies between some of the previous results. We discuss convergence criteria for the required force resolution, simulation box size, halo mass range, initial and final redshift, and time stepping. Because of conservative cuts on the mass range probed by individual boxes, our results are relatively insensitive to simulation volume, the remaining sensitivity being consistent with extended Press-Schechter theory. Previously obtained mass function fits near z=0, when scaled by linear theory, are in good agreement with our results at all redshifts, although a mild redshift dependence consistent with that found by Reed and collaborators exists at low redshifts.Comment: 20 pages, 15 figures. Minor changes to the text and figures; results and conclusions unchange

X-ray Observations of Optically Selected Giant Elliptical-Dominated Galaxy Groups

We present a combined optical and X-ray analysis of three optically selected X-ray bright groups with giant elliptical galaxies in the center. These massive ellipticals were targeted for XMM-Newton X-ray observations based on their large velocity dispersions and their proximity to a nearby ROSAT X-ray source. Additionally, these targets are significantly brighter in the optical than their nearest neighbors. We show that one of these systems meets the standard criteria for a fossil group. While the other two systems have a prominent magnitude gap in the E/S0 ridgeline, they do not appear to have reached the fossil-like final stage of group evolution.Comment: 8 pages, 6 figures, Accepted for publication in Ap

Robustness of Cosmological Simulations I: Large Scale Structure

The gravitationally-driven evolution of cold dark matter dominates the formation of structure in the Universe over a wide range of length scales. While the longest scales can be treated by perturbation theory, a fully quantitative understanding of nonlinear effects requires the application of large-scale particle simulation methods. Additionally, precision predictions for next-generation observations, such as weak gravitational lensing, can only be obtained from numerical simulations. In this paper, we compare results from several N-body codes using test problems and a diverse set of diagnostics, focusing on a medium resolution regime appropriate for studying many observationally relevant aspects of structure formation. Our conclusions are that -- despite the use of different algorithms and error-control methodologies -- overall, the codes yield consistent results. The agreement over a wide range of scales for the cosmological tests is test-dependent. In the best cases, it is at the 5% level or better, however, for other cases it can be significantly larger than 10%. These include the halo mass function at low masses and the mass power spectrum at small scales. While there exist explanations for most of the discrepancies, our results point to the need for significant improvement in N-body errors and their understanding to match the precision of near-future observations. The simulation results, including halo catalogs, and initial conditions used, are publicly available.Comment: 32 pages, 53 figures, data from the simulations is available at http://t8web.lanl.gov/people/heitmann/arxiv, accepted for publication in ApJS, several minor revisions, reference added, main conclusions unchange

Imposing a Lagrangian Particle Framework on an Eulerian Hydrodynamics Infrastructure in Flash

In many astrophysical simulations, both Eulerian and Lagrangian quantities are of interest. For example, in a galaxy cluster merger simulation, the intracluster gas can have Eulerian discretization, while dark matter can be modeled using particles. FLASH, a component-based scientific simulation code, superimposes a Lagrangian framework atop an adaptive mesh refinement Eulerian framework to enable such simulations. The discretization of the field variables is Eulerian, while the Lagrangian entities occur in many different forms including tracer particles, massive particles, charged particles in particle-in-cell mode, and Lagrangian markers to model fluid structure interactions. These widely varying roles for Lagrangian entities are possible because of the highly modular, flexible, and extensible architecture of the Lagrangian framework. In this paper, we describe the Lagrangian framework in FLASH in the context of two very different applications, Type Ia supernovae and galaxy cluster mergers, which use the Lagrangian entities in fundamentally different ways

Simulations of Astrophysical Fluid Instabilities

We present direct numerical simulations of mixing at Rayleigh-Taylor unstable interfaces performed with the FLASH code, developed at the ASCI/Alliances Center for Astrophysical Thermonuclear Flashes at the University of Chicago. We present initial results of single-mode studies in two and three dimensions. Our results indicate that three-dimensional instabilities grow significantly faster than two-dimensional instabilities and that grid resolution can have a significant effect on instability growth rates. We also find that unphysical diffusive mixing occurs at the fluid interface, particularly in poorly resolved simulations.Comment: 3 pages, 1 figure. To appear in the proceedings of the 20th Texas Symposium on Relativistic Astrophysic

Large-Scale Simulations of Clusters of Galaxies

We discuss some of the computational challenges encountered in simulating the evolution of clusters of galaxies. Eulerian adaptive mesh refinement (AMR) techniques can successfully address these challenges but are currently being used by only a few groups. We describe our publicly available AMR code, FLASH, which uses an object-oriented framework to manage its AMR library, physics modules, and automated verification. We outline the development of the FLASH framework to include collisionless particles, permitting it to be used for cluster simulation.Comment: 3 pages, 3 figures, to appear in Proceedings of the VII International Workshop on Advanced Computing and Analysis Techniques in Physics Research (ACAT 2000), Fermilab, Oct. 16-20, 200

Microscopic theory of glassy dynamics and glass transition for molecular crystals

We derive a microscopic equation of motion for the dynamical orientational correlators of molecular crystals. Our approach is based upon mode coupling theory. Compared to liquids we find four main differences: (i) the memory kernel contains Umklapp processes, (ii) besides the static two-molecule orientational correlators one also needs the static one-molecule orientational density as an input, where the latter is nontrivial, (iii) the static orientational current density correlator does contribute an anisotropic, inertia-independent part to the memory kernel, (iv) if the molecules are assumed to be fixed on a rigid lattice, the tensorial orientational correlators and the memory kernel have vanishing l,l'=0 components. The resulting mode coupling equations are solved for hard ellipsoids of revolution on a rigid sc-lattice. Using the static orientational correlators from Percus-Yevick theory we find an ideal glass transition generated due to precursors of orientational order which depend on X and p, the aspect ratio and packing fraction of the ellipsoids. The glass formation of oblate ellipsoids is enhanced compared to that for prolate ones. For oblate ellipsoids with X <~ 0.7 and prolate ellipsoids with X >~ 4, the critical diagonal nonergodicity parameters in reciprocal space exhibit more or less sharp maxima at the zone center with very small values elsewhere, while for prolate ellipsoids with 2 <~ X <~ 2.5 we have maxima at the zone edge. The off-diagonal nonergodicity parameters are not restricted to positive values and show similar behavior. For 0.7 <~ X <~ 2, no glass transition is found. In the glass phase, the nonergodicity parameters show a pronounced q-dependence.Comment: 17 pages, 12 figures, accepted at Phys. Rev. E. v4 is almost identical to the final paper version. It includes, compared to former versions v2/v3, no new physical content, but only some corrected formulas in the appendices and corrected typos in text. In comparison to version v1, in v2-v4 some new results have been included and text has been change