The distribution of pairwise peculiar velocities in the nonlinear regime

The distribution of pairwise, relative peculiar velocities, $f(u;r)$, on small nonlinear scales, $r$, is derived from the Press--Schechter approach. This derivation assumes that Press--Schechter clumps are virialized and isothermal. The virialized assumption requires that the circular velocity, $V_c \propto M^{1/3}$, where $M$ denotes the mass of the clump. The isothermal assumption means that the circular velocity is independent of radius. Further, it is assumed that the velocity distribution within a clump is Maxwellian, that the pairwise relative velocity distribution is isotropic, and that on nonlinear scales clump-clump motions are unimportant when calculating the distribution of velocity differences. Comparison with $N$-body simulations shows that, on small nonlinear scales, all these assumptions are accurate. For most power spectra of interest, the resulting line of sight, pairwise, relative velocity distribution, $f(u_{\rm r})$, is well approximated by an exponential, rather than a Gaussian distribution. This simple Press--Schechter model is also able to provide a natural explanation for the observed, non-Maxwellian shape of $f(v)$, the distribution of peculiar velocities.Comment: (MNRAS, in press) 16 pages, uuencode

Deuterium Abundance in the Most Metal-Poor Damped Lyman alpha System: Converging on Omega_baryons

The most metal-poor DLA known to date, at z = 2.61843 in the spectrum of the QSO Q0913+072, with an oxygen abundance only about 1/250 of the solar value, shows six well resolved D I Lyman series transitions in high quality echelle spectra recently obtained with the ESO VLT. We deduce a value of the deuterium abundance log (D/H) = -4.56+/-0.04 which is in good agreement with four out of the six most reliable previous determinations of this ratio in QSO absorbers. We find plausible reasons why in the other two cases the 1 sigma errors may have been underestimated by about a factor of two. The addition of this latest data point does not change significantly the mean value of the primordial abundance of deuterium, suggesting that we are now converging to a reliable measure of this quantity. We conclude that = -4.55+/-0.03 and Omega_b h^2 (BBN) = 0.0213+/-0.0010 (68% confidence limits). Including the latter as a prior in the analysis of the five year data of WMAP leads to a revised best-fitting value of the power-law index of primordial fluctuations n_s = 0.956+/-0.013 (1 sigma) and n_s < 0.990 with 99% confidence. Considering together the constraints provided by WMAP 5, (D/H)_p, baryon oscillations in the galaxy distribution, and distances to Type Ia supernovae, we arrive at the current best estimates Omega_b h^2 = 0.0224+/-0.0005 and n_s = 0.959+/-0.013.Comment: 13 pages, 8 Figures. Revised version following referee's comments. Accepted for publication in Monthly Notices of the Royal Astronomical Society. A few typos correcte

The host galaxies of strong CaII QSO absorption systems at z<0.5

We present new imaging and spectroscopic observations of the fields of five QSOs with very strong intervening CaII absorption systems at redshifts z<0.5 selected from the Sloan Digital Sky Survey. Recent studies of these very rare absorbers indicate that they may be related to damped Lyman alpha systems (DLAs). In all five cases we identify a galaxy at the redshift of the CaII system with impact parameters up to ~24 kpc. In four out of five cases the galaxies are luminous (L ~L*), metal-rich (Z ~Zsun), massive (velocity dispersion, sigma ~100 km/s) spirals. Their star formation rates, deduced from Halpha emission, are high, in the range SFR = 0.3 - 30 Msun/yr. In our analysis, we paid particular attention to correcting the observed emission line fluxes for stellar absorption and dust extinction. We show that these effects are important for a correct SFR estimate; their neglect in previous low-z studies of DLA-selected galaxies has probably led to an underestimate of the star formation activity in at least some DLA hosts. We discuss possible links between CaII-selected galaxies and DLAs and outline future observations which will help clarify the relationship between these different classes of QSO absorbers.Comment: Accepted for publication in MNRAS, 14 pages, 9 figures. Version with full resolution images available at http://www.ast.cam.ac.uk/~bjz/papers/Zych_etal_2007a.pd

C, N, O Abundances in the Most Metal-Poor Damped Lyman alpha Systems

This study focuses on some of the most metal-poor damped Lyman alpha absorbers known in the spectra of high redshift QSOs, using new and archival observations obtained with UV-sensitive echelle spectrographs on the Keck and VLT telescopes. The weakness and simple velocity structure of the absorption lines in these systems allows us to measure the abundances of several elements, and in particular those of C, N, and O, a group that is difficult to study in DLAs of more typical metallicities. We find that when the oxygen abundance is less than about 1/100 of solar, the C/O ratio in high redshift DLAs and sub-DLAs matches that of halo stars of similar metallicity and shows higher values than expected from galactic chemical evolution models based on conventional stellar yields. Furthermore, there are indications that at these low metallicities the N/O ratio may also be above simple expectations and may exhibit a minimum value, as proposed by Centurion and her collaborators in 2003. Both results can be interpreted as evidence for enhanced production of C and N by massive stars in the first few episodes of star formation, in our Galaxy and in the distant proto-galaxies seen as QSO absorbers. The higher stellar yields implied may have an origin in stellar rotation which promotes mixing in the stars' interiors, as considered in some recent model calculations. We briefly discuss the relevance of these results to current ideas on the origin of metals in the intergalactic medium and the universality of the stellar initial mass function.Comment: 17 pages, 9 Figures, Accepted for publication in Monthly Notices of the Royal Astronomical Societ

Relativistic MHD with Adaptive Mesh Refinement

This paper presents a new computer code to solve the general relativistic magnetohydrodynamics (GRMHD) equations using distributed parallel adaptive mesh refinement (AMR). The fluid equations are solved using a finite difference Convex ENO method (CENO) in 3+1 dimensions, and the AMR is Berger-Oliger. Hyperbolic divergence cleaning is used to control the $\nabla\cdot {\bf B}=0$ constraint. We present results from three flat space tests, and examine the accretion of a fluid onto a Schwarzschild black hole, reproducing the Michel solution. The AMR simulations substantially improve performance while reproducing the resolution equivalent unigrid simulation results. Finally, we discuss strong scaling results for parallel unigrid and AMR runs.Comment: 24 pages, 14 figures, 3 table

Discovery of GeV Emission from the Circinus galaxy with the Fermi-LAT

We report the discovery of gamma-ray emission from the Circinus galaxy using the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope. Circinus is a nearby (~4 Mpc) starburst with a heavily obscured Seyfert-type active nucleus, bipolar radio lobes perpendicular to the spiral disk, and kpc-scale jet-like structures. Our analysis of 0.1-100 GeV events collected during 4 years of LAT observations reveals a significant (~ 7.3 sigma) excess above the background. We find no indications of variability or spatial extension beyond the LAT point-spread function. A power-law model used to describe the 0.1-100 GeV gamma-ray spectrum yields a flux of (18.8+/-5.8)x10^{-9} ph cm^{-2} s^{-1} and photon index 2.19+/-0.12, corresponding to an isotropic gamma-ray luminosity of 3 x 10^{40} erg s^{-1}. This observed gamma-ray luminosity exceeds the luminosity expected from cosmic-ray interactions in the interstellar medium and inverse Compton radiation from the radio lobes. Thus the origin of the GeV excess requires further investigation.Comment: 7 pages, 7 figures, accepted for publication in the Astrophysical Journa

The relativistic Sagnac Effect: two derivations

The phase shift due to the Sagnac Effect, for relativistic matter and electromagnetic beams, counter-propagating in a rotating interferometer, is deduced using two different approaches. From one hand, we show that the relativistic law of velocity addition leads to the well known Sagnac time difference, which is the same independently of the physical nature of the interfering beams, evidencing in this way the universality of the effect. Another derivation is based on a formal analogy with the phase shift induced by the magnetic potential for charged particles travelling in a region where a constant vector potential is present: this is the so called Aharonov-Bohm effect. Both derivations are carried out in a fully relativistic context, using a suitable 1+3 splitting that allows us to recognize and define the space where electromagnetic and matter waves propagate: this is an extended 3-space, which we call "relative space". It is recognized as the only space having an actual physical meaning from an operational point of view, and it is identified as the 'physical space of the rotating platform': the geometry of this space turns out to be non Euclidean, according to Einstein's early intuition.Comment: 49 pages, LaTeX, 3 EPS figures. Revised (final) version, minor corrections; to appear in "Relativity in Rotating Frames", ed. G. Rizzi and M.L. Ruggiero, Kluwer Academic Publishers, Dordrecht, (2003). See also http://digilander.libero.it/solciclo