Statistical distribution of gravitational-lensing excursion angles: winding ways to us from the deep Universe

We investigate statistical distributions of differences in gravitational-lensing deflections between two light rays, the so-called lensing excursion angles. A probability distribution function of the lensing excursion angles, which plays a key role in estimates of lensing effects on angular clustering of objects (such as galaxies, quasi-stellar objects and also the cosmic microwave background temperature map), is known to consist of two components: a Gaussian core and an exponential tail. We use numerical gravitational-lensing experiments in a ΛCDM cosmology for quantifying these two components. We especially focus on the physical processes responsible for generating these two components. We develop a simple empirical model for the exponential tail which allows us to explore its origin. We find that the tail is generated by the coherent lensing scatter by massive haloes with   M > 10¹⁴  h ⁻¹ M_⊙  at   z < 1  and that its exponential shape arises due to the exponential cut-off of the halo mass function at that mass range. On scales larger than 1 arcmin, the tail does not have a practical influence on the lensing effects on the angular clustering. Our model predicts that the coherent scatter may have non-negligible effects on angular clustering at subarcminute scales.     

ZEN2: a narrow J-band search for z∼ 9 Lyα emitting galaxies directed towards three lensing clusters

We present the results of a continuing survey to detect Lyα emitting galaxies at redshifts   z ∼ 9  : the ' z equals nine' (ZEN) survey. We have obtained deep VLT Infrared Spectrometer and Array Camera observations in the narrow J -band filter NB119 directed towards three massive lensing clusters: Abell clusters 1689, 1835 and 114. The foreground clusters provide a magnified view of the distant Universe and permit a sensitive test for the presence of very high redshift galaxies. We search for   z ∼ 9 Lyα  emitting galaxies displaying a significant narrow-band excess relative to accompanying J -band observations that remain undetected in Hubble Space Telescope ( HST )/Advanced Camera for Surveys (ACS) optical images of each field. No sources consistent with this criterion are detected above the unlensed 90 per cent point-source flux limit of the narrow-band image,   F _NB= 3.7 × 10⁻¹⁸ erg s⁻¹ cm⁻²  . To date, the total coverage of the ZEN survey has sampled a volume at   z ∼ 9  of approximately 1700 comoving Mpc³ to a Lyα emission luminosity of  10⁴³ erg s⁻¹  . We conclude by considering the prospects for detecting   z ∼ 9 Lyα  emitting galaxies in light of both observed galaxy properties at   z < 7  and simulated populations at   z > 7  .     

The dynamical state of RX J1347.5−1145 from a combined strong lensing and X-ray analysis

We perform a combined X-ray and strong lensing analysis of RX J1347.5−1145, one of the most luminous galaxy clusters at X-ray wavelengths. We show that evidence from strong lensing alone, based on published Very Large Telescope (VLT) and new Hubble Space Telescope ( HST ) data, strongly argues in favour of a complex structure. The analysis takes into account arc positions, shapes and orientations, and is done thoroughly in the image plane. The cluster inner regions are well fitted by a bimodal mass distribution, with a total projected mass of   M _tot= (9.9 ± 0.3) × 10¹⁴ M_⊙  h ⁻¹  within a radius of 360 kpc  h ⁻¹ (1.5 arcmin). Such a complex structure could be a signature of a recent major merger as further supported by X-ray data. A temperature map of the cluster, based on deep Chandra observations, reveals a hot front located between the first main component and an X-ray emitting south-eastern subclump. The map also unveils a filament of cold gas in the innermost regions of the cluster, most probably a cooling wake caused by the motion of the cD inside the cool core region. A merger scenario in the plane of the sky between two dark matter subclumps is consistent with both our lensing and X-ray analyses, and can explain previous discrepancies with mass estimates based on the virial theorem.     

Microlensing by cosmic strings

We consider the signature and detectability of gravitational microlensing of distant quasars by cosmic strings. Because of the simple image configuration such events will have a characteristic lightcurve, in which a source would appear to brighten by exactly a factor of 2, before reverting to its original apparent brightness. We calculate the optical depth and event rate, and conclude that current predictions and limits on the total length of strings on the sky imply optical depths of ≲ 10⁻⁸ and event rates of fewer than one event per 10⁹ sources per year. Disregarding those predictions but replacing them with limits on the density of cosmic strings from the cosmic microwave background fluctuation spectrum, leaves only a small region of parameter space (in which the sky contains about 3 × 10⁵ strings with deficit angle of the order of 0.3 milli-seconds) for which a microlensing survey of exposure 10⁷ source years, spanning a 20–40-year period, might reveal the presence of cosmic strings.     

Sakurai's Object: characterizing the near-infrared CO ejecta between 2003 and 2007

We present observations of Sakurai's Object obtained at 1–5 μm between 2003 and 2007. By fitting a radiative transfer model to an echelle spectrum of CO fundamental absorption features around  4.7 μm  , we determine the excitation conditions in the line-forming region. We find  ¹²C/¹³C = 3.5^+2.0_−1.5  , consistent with CO originating in ejecta processed by the very late thermal pulse, rather than in the pre-existing planetary nebula. We demonstrate the existence of  2.2 × 10⁻⁶≤ M _CO≤ 2.7 × 10⁻⁶ M_⊙  of CO ejecta outside the dust, forming a high-velocity wind of  500 ± 80 km s⁻¹  . We find evidence for significant weakening of the CO band and cooling of the dust around the central star between 2003 and 2005. The gas and dust temperatures are implausibly high for stellar radiation to be the sole contributor.     

Relativistic ionized accretion disc models of MCG–6-30-15

We present BeppoSAX observations of Nova Velorum 1999 (V382 Vel), carried out in a broad X-ray band covering 0.1–300 keV only 15 d after the discovery and again after 6 months. The nova was detected at day 15 with the BeppoSAX instruments which measured a flux F _x≃1.8×10⁻¹¹ erg cm⁻² s⁻¹ in the 0.1–10 keV range and a 2 σ upper limit F _x<6.7×10⁻¹² erg cm⁻² s⁻¹ in the 15–60 keV range. We attribute the emission to shocked nebular ejecta at a plasma temperature kT ≃6 keV . At six months no bright component emerged in the 15–60 keV range, but a bright central supersoft X-ray source appeared. The hot nebular component previously detected had cooled to a plasma temperature kT <1 keV . There was strong intrinsic absorption of the ejecta in the first observation and not in the second, because the column density of neutral hydrogen decreased from N (H)≃1.7×10²³ to N (H)≃10²¹ cm⁻² (close to the interstellar value). The unabsorbed X-ray flux also decreased from F _x=4.3×10⁻¹¹ to F _x≃10⁻¹² erg cm⁻² s⁻¹ .     

Neutral sodium atoms release from the surface of the Moon induced by meteoroid impacts

In this work, we calculate the vapour and neutral Na production rates on the Moon, as due to the impacts of meteoroids in the radius range of 10⁻⁸–0.15 m. We limit our calculations to this size range, since meteoroids with radius larger than 0.15 m have not been found to be important for the production of the exosphere in a time interval comparable with that of the observations.
We have considered a new dynamical model of the flux of meteoroids at the heliocentric distance of the Moon, regarding objects in the radius range of 10⁻²–0.15 m. Instead, the flux of smaller meteoroids (radius range 10⁻⁸–10⁻² m) has been calculated using the two distributions adopted by Cintala and Love & Brownlee.
The results of our model are that (i) the neutral Na production rate is  ∼3–4.9 × 10⁴ atoms cm⁻² s⁻¹  , slightly larger than the previous estimates  (∼2–3 × 10⁴ atoms cm⁻² s⁻¹)  , and (ii) only about 6 per cent of neutral Na is produced by the impacts of meteoroids in the size range 10⁻³–0.15 m, whereas about 94 per cent of the Na comes from the  10⁻⁵–10⁻³ m  size range.     

Collapse and fragmentation of rotating magnetized clouds – I. Magnetic flux–spin relation

We discuss the evolution of the magnetic flux density and angular velocity in a molecular cloud core, on the basis of three-dimensional numerical simulations, in which a rotating magnetized cloud fragments and collapses to form a very dense optically thick core of  >5 × 10¹⁰ cm⁻³  . As the density increases towards the formation of the optically thick core, the magnetic flux density and angular velocity converge towards a single relationship between the two quantities. If the core is magnetically dominated its magnetic flux density approaches  1.5( n /5 × 10¹⁰ cm⁻³)^1/2 mG  , while if the core is rotationally dominated the angular velocity approaches  2.57 × 10⁻³ ( n /5 × 10¹⁰ cm⁻³)^1/2 yr⁻¹  , where n is the density of the gas. We also find that the ratio of the angular velocity to the magnetic flux density remains nearly constant until the density exceeds  5 × 10¹⁰ cm⁻³  . Fragmentation of the very dense core and emergence of outflows from fragments will be shown in the subsequent paper.     

Chandra measurements of the distribution of mass in the luminous lensing cluster Abell 2390

We present spatially resolved X-ray spectroscopy of the luminous lensing cluster Abell 2390, using observations made with the Chandra observatory. The temperature of the X-ray gas rises with increasing radius within the central ∼ 200 kpc of the cluster, and then remains approximately isothermal, with kT =11.5_−1.6^+1.5 keV , out to the limits of the observations at r ∼1.0 Mpc . The total mass profile determined from the Chandra data has a form in good agreement with the predictions from numerical simulations. Using the parametrization of Navarro, Frenk and White, we measure a scale radius r _s∼0.8 Mpc and a concentration parameter c ∼3 . The best-fitting X-ray mass model is in good agreement with independent gravitational lensing results and optical measurements of the galaxy velocity dispersion in the cluster. The X-ray gas to total mass ratio rises with increasing radius with f _gas∼21 per cent at r =0.9 Mpc . The azimuthally averaged 0.3–7.0 keV surface brightness profile exhibits a small core radius and a clear 'break' at r ∼500 kpc , where the slope changes from S _X ∼^∝  r ^−1.5 to S _X ∼^∝  r ^−3.6 . The data for the central region of the cluster indicate the presence of a cooling flow with a mass deposition rate of 200–300 M_⊙ yr⁻¹ and an effective age of 2–3 Gyr .     

An enlarged cosmic shear survey with the William Herschel Telescope

We report the results of a cosmic shear survey using the 4.2-m William Herschel Telescope on La Palma, to a depth of   R = 25.8 ( z ≈ 0.8)  , over 4 deg². The shear correlation functions are measured on scales from 1 to 15 arcmin, and are used to constrain cosmological parameters. We ensure that our measurements are free from instrumental systematic effects by performing a series of tests, including a decomposition of the signal into E - and B -modes. We also reanalyse the data independently, using the shear measurement pipeline developed for the COMBO-17 survey. This confirms our results and also highlights various effects introduced by different implementations of the basic 'Kaiser–Squires–Broadhurst' shear measurement method. We find that the normalization of the matter power spectrum on 8  h ⁻¹ Mpc scales is  σ₈= (1.02 ± 0.15)(0.3/Ω _m )^1/2  , where the 68 per cent confidence limit error includes noise, sample variance, covariance between angular scales, systematic effects, redshift uncertainty and marginalization over other parameters. We compare these results with other cosmic shear surveys and with recent constraints from the Wilkinson Microwave Anisotropy Probe experiment.     

Understanding the halo-mass and galaxy-mass cross-correlation functions

We use the Millennium Simulation (MS) to measure the cross-correlation between halo centres and mass (or equivalently the average density profiles of dark haloes) in a Lambda cold dark matter (ΛCDM) cosmology. We present results for radii in the range  10  h ⁻¹ kpc < r < 30  h ⁻¹ Mpc  and for halo masses in the range  4 × 10¹⁰ < M ₂₀₀ < 4 × 10¹⁴  h ⁻¹ M_⊙  . Both at   z = 0  and at   z = 0.76  these cross-correlations are surprisingly well fitted if the inner region is approximated by a density profile of NFW or Einasto form, the outer region by a biased version of the linear mass autocorrelation function, and the maximum of the two is adopted where they are comparable. We use a simulation of galaxy formation within the MS to explore how these results are reflected in cross-correlations between galaxies and mass. These are directly observable through galaxy–galaxy lensing. Here also we find that simple models can represent the simulation results remarkably well, typically to ≲10 per cent. Such models can be used to extend our results to other redshifts, to cosmologies with other parameters, and to other assumptions about how galaxies populate dark haloes. Our galaxy formation simulation already reproduces current galaxy–galaxy lensing data quite well. The characteristic features predicted in the galaxy–galaxy lensing signal should provide a strong test of the ΛCDM cosmology as well as a route to understanding how galaxies form within it.     

Measurement of the Sunyaev–Zel'dovich increment in massive galaxy clusters

We have detected the Sunyaev–Zel'dovich (SZ) increment at 850 μm in two galaxy clusters (Cl 0016+16 and MS 1054.4−0321) using the Submillimetre Common User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope. Fits to the isothermal β model yield a central Compton y parameter of  (2.2 ± 0.7) × 10⁻⁴  and a central 850-μm flux of  Δ I ₀= 2.2 ± 0.7 mJy beam⁻¹  in Cl 0016. This can be combined with decrement measurements to infer   y = (2.38 ±^0.36_0.34) × 10⁻⁴  and   v _pec= 400±¹⁹⁰⁰₁₄₀₀ km s⁻¹  . In MS 1054 we find a peak 850-μm flux of  Δ I ₀= 2.0 ± 1.0 mJy beam⁻¹  and   y = (2.0 ± 1.0) × 10⁻⁴  . To be successful such measurements require large chop throws and non-standard data analysis techniques. In particular, the 450-μm data are used to remove atmospheric variations in the 850-μm data. An explicit annular model is fit to the SCUBA difference data in order to extract the radial profile, and separately fit to the model differences to minimize the effect of correlations induced by our scanning strategy. We have demonstrated that with sufficient care, SCUBA can be used to measure the SZ increment in massive, compact galaxy clusters.     

The continuing saga of Sakurai's object (V4334 Sgr): dust production and helium line emission

We report further UKIRT spectroscopic observations of Sakurai's object (V4334 Sgr) made in 1999 April/May in the 1–4.75 μm range, and find that the emission is dominated by amorphous carbon at T _d~600 K. The estimated maximum grain size is 0.6 μm, and the mass lower limit is 1.7±0.2×10⁻⁸ M_⊙ to 8.9±0.6×10⁻⁷ M_⊙ for distances of 1.1–8 kpc. For 3.8 kpc the mass is 2.0±0.1×10⁻⁷ M_⊙.
We also report strong He  i emission at 1.083 μm, in contrast to the strong absorption in this line in 1998. We conclude that the excitation is collisional, and is probably caused by a wind, consistent with the P Cygni profile observed by Eyres et al. in 1998.     

Hyperluminous infrared galaxies

39 galaxies are now known, from follow-up of faint IRAS sources and from submillimetre observations of high-redshift AGN, with far-infrared luminosities >10¹³ L_⊙. 13 of these, which have been found in 60- or 850-μm surveys, form an important unbiased subsample. 12 have been found by comparison of 60-μm surveys with quasar or radio galaxy catalogues, or from infrared surveys with colour selection biased towards AGN, while a further 14 have been found through submillimetre observations of known high-redshift AGN. In this paper I argue, on the basis of detailed modelling of the spectral energy distributions of hyperluminous galaxies with accurate radiative transfer models, and from evidence of high gas mass in several cases, that the bulk of the emission from these galaxies at rest frame wavelengths ≥50 μm is caused by star formation. Even after correction for the effects of lensing, hyperluminous galaxies with emission peaking at rest frame wavelengths ≥50 μm are therefore undergoing star formation at rates >10³ M_⊙ yr⁻¹ and are strong candidates for being primeval galaxies, in the process of a major episode of star formation.     

Lensing, reddening and extinction effects of Mg ii absorbers from z= 0.4 to 2

Using a sample of almost 7000 strong Mg  ii absorbers with   W ₀ > 1 Å  and  0.4 < z < 2.2  detected in the SDSS DR4 data set, we investigate the gravitational lensing and dust extinction effects they induce on background quasars. After carefully quantifying several selection biases, we isolate the reddening effects as a function of redshift and absorber rest equivalent width, W ₀. We find the amount of dust to increase with cosmic time as  τ( z ) ∝ (1 + z )^−1.1±0.4  , following the evolution of cosmic star density or integrated star formation rate. We measure the reddening effects over a factor of 30 in E ( B − V ) and we find that  τ∝ ( W ₀)^1.9±0.1  , providing us with an important scaling for theoretical modelling of metal absorbers. We also measure the dust-to-metal ratio and find it similar to that of the Milky Way. In contrast to previous studies, we do not detect any gravitational magnification by Mg  ii systems. We measure the upper limit  μ < 1.10  and discuss the origin of the discrepancy. Finally, we estimate the fraction of absorbers missed due to extinction effects and show that it rises from 1 to 50 per cent in the range  1 < W ₀ < 6 Å  . We parametrize this effect and provide a correction for recovering the intrinsic  ∂ N /∂ W ₀  distribution.     

RR Pic (1925): a Chandra X-ray view

We present the Chandra ACIS-S3 data of the old classical nova RR Pic (1925). The source has a count rate of 0.067 ± 0.002 count s⁻¹ in the 0.3–5.0 keV energy range. We detect the orbital period of the underlying binary system in the X-ray wavelengths. We also find that the neutral hydrogen column density differs for orbital minimum and orbital maximum spectra with values  0.25^+0.23_−0.18× 10²²  and  0.64^+0.13_−0.14× 10²² cm⁻²  at 3σ confidence level. The X-ray spectrum of RR Pic can be represented by a composite model of bremsstrahlung with a photoelectric absorption, two absorption lines centered around 1.1–1.4 keV and five Gaussian lines centered at emission lines around 0.3–1.1 keV corresponding to various transitions of S, N, O, C, Ne and Fe. The bremsstrahlung temperature derived from the fits ranges from 0.99 to 1.60 keV and the unabsorbed X-ray flux is found to be  2.5^+0.4_−1.2× 10⁻¹³ erg  cm⁻² s⁻¹  in the 0.3–5.0 keV range with a luminosity of 1.1 ± 0.2  10³¹ erg s⁻¹  at 600 pc. We also detect excess emission in the spectrum possibly originating from the reverse shock in the ejecta. A fit with a cooling flow plasma emission model shows enhanced abundances of He, C, N, O and Ne in the X-ray emitting region indicating existence of diffusive mixing.     

Probing the submillimetre number counts at f850 μm < 2 mJy

We have conducted a submillimetre mapping survey of faint, gravitationally lensed sources, where we have targeted 12 galaxy clusters and additionally the New Technology Telescope (NTT) Deep Field. The total area surveyed is 71.5 arcmin² in the image plane; correcting for gravitational lensing, the total area surveyed is 40 arcmin² in the source plane for a typical source redshift z ≈ 2.5. In the deepest maps, an image plane depth of 1σ rms ∼0.8 mJy is reached. This survey is the largest survey to date to reach such depths. In total 59 sources were detected, including three multiply imaged sources. The gravitational lensing makes it possible to detect sources with flux density below the blank field confusion limit. The lensing-corrected fluxes range from 0.11 to 19 mJy. After correcting for multiplicity, there are 10 sources with fluxes <2 mJy of which seven have submJy fluxes, doubling the number of such sources known. Number counts are determined below the confusion limit. At 1 mJy, the integrated number count is  ∼10⁴ deg⁻²  , and at 0.5 mJy it is  ∼2 × 10⁴ deg⁻²  . Based on the number counts, at a source plan flux limit of 0.1 mJy, essentially all of the 850-μm background emission has been resolved. The dominant contribution (>50 per cent) to the integrated background arises from sources with fluxes S ₈₅₀ between 0.4 and 2.5 mJy, while the bright sources S ₈₅₀ > 6 mJy contribute only 10 per cent.     

The coarse geometry of merger trees in Λ cold dark matter

We introduce the contour process to describe the geometrical properties of merger trees. The contour process translates a tree structure into a one-dimensional object: the contour walk. We characterize the contour walk measuring its length and action. The length is proportional to the number of progenitors in the tree, and the action is a proxy for the mean length of a branch in the tree.
We construct the contour walk for merger trees extracted from the public data base of the Millennium Run and also for merger trees constructed with a public Monte Carlo code which implements a Markovian algorithm. The trees correspond to haloes of final masses between 10¹¹ and  10¹⁴  h ⁻¹ M_⊙  . We study how the length and action scale with the mass of the final halo.
In all cases, except for the action measured from Markovian trees, we find a transitional scale around  1–6 × 10¹²  h ⁻¹ M_⊙  , which can be interpreted as the point where the halo growth starts to be dominated by mergers and not accretion. The transitional scale depends weakly on the mass resolution and time discretization of the merging trees.     

XMM–Newton and Chandra observations of SHEEP sources

We present Chandra and XMM–Newton observations of 12 bright  [ f (2–10 keV) > 10⁻¹³ erg cm⁻² s⁻¹]  sources from the ASCA search for the High Energy Extragalactic Population (SHEEP) survey. Most of these have been either not observed or not detected previously with the ROSAT mission, and therefore they constitute a sample biased towards hard sources. The Chandra observations are important in locating the optical counterpart of the X-ray sources with accuracy. Optical spectroscopic observations show that our sample is associated with both narrow-line (NL) (six objects) and broad-line (BL) active galactic nuclei (AGN) (five objects), with one source remaining unidentified. Our sources cover the redshift range 0.04–1.29, spanning luminosities from 10⁴² to  10⁴⁵ erg s⁻¹  (2–10 keV). The NL sources have preferentially lower redshift (and luminosity) compared to the BL ones. This can be most easily explained in a model where the NL AGN are intrinsically less luminous than the BL ones in line with the results of Steffen et al. The X-ray spectral fittings show a roughly equal number of obscured  ( N _H > 10²² cm⁻²)  and unobscured  ( N _H < 10²² cm⁻²)  sources. There is a clear tendency for obscured sources to be associated with NL AGN and unobscured sources with BL ones. However, there is a marked exception with the highest obscuring column observed at a BL AGN at a redshift of z = 0.5.     

Gravitational lens magnification by Abell 1689: distortion of the background galaxy luminosity function

Gravitational lensing magnifies the observed flux of galaxies behind the lens. We use this effect to constrain the total mass in the cluster Abell 1689 by comparing the lensed luminosities of background galaxies with the luminosity function of an undistorted field. Under the assumption that these galaxies are a random sample of luminosity space, this method is not limited by clustering noise. We use photometric redshift information to estimate galaxy distance and intrinsic luminosity. Knowing the redshift distribution of the background population allows us to lift the mass/background degeneracy common to lensing analysis. In this paper we use nine filters observed over 12 h with the Calar Alto 3.5-m telescope to determine the redshifts of 1000 galaxies in the field of Abell 1689. Using a complete sample of 146 background galaxies we measure the cluster mass profile. We find that the total projected mass interior to 0.25  h ⁻¹ Mpc is M _2D(<0.25  h ⁻¹ Mpc)=(0.48±0.16)×10¹⁵  h ⁻¹ M_⊙, where our error budget includes uncertainties from the photometric redshift determination, the uncertainty in the offset calibration and finite sampling. This result is in good agreement with that found by number-count and shear-based methods and provides a new and independent method to determine cluster masses.