Intrinsic correlation of galaxy shapes: implications for weak lensing measurements

Weak gravitational lensing is now established as a powerful method to measure mass fluctuations in the universe. It relies on the measurement of small coherent distortions of the images of background galaxies. Even low-level correlations in the intrinsic shapes of galaxies could however produce a significant spurious lensing signal. These correlations are also interesting in their own right, since their detection would constrain models of galaxy formation. Using     haloes found in N -body simulations, we compute the correlation functions of the intrinsic ellipticity of spiral galaxies assuming that the disc is perpendicular to the angular momentum of the dark matter halo. We also consider a simple model for elliptical galaxies, in which the shape of the dark matter halo is assumed to be the same as that of the light. For deep lensing surveys with median redshifts ∼1, we find that intrinsic correlations of ∼10⁻⁴ on angular scales     are generally below the expected lensing signal, and contribute only a small fraction of the excess signals reported on these scales. On larger scales we find limits to the intrinsic correlation function at a level ∼10⁻⁵, which gives a (model-dependent) range of separations for which the intrinsic signal is about an order of magnitude below the ellipticity correlation function expected from weak lensing. Intrinsic correlations are thus negligible on these scales for dedicated weak lensing surveys. For wider but shallower surveys such as SuperCOSMOS, APM and SDSS, we cannot exclude the possibility that intrinsic correlations could dominate the lensing signal. We discuss how such surveys could be used to calibrate the importance of this effect, as well as study spin–spin correlations of spiral galaxies.     

Large-Scale Coherent Dipole Anisotropy?

We have reanalysed and compared the dipoles of the 1.2-Jy and 0.6-Jy (QDOT) IRAS galaxy samples. We find strong indications from both samples for (a) significant contributions to the gravitational field that shapes the Local Group motion from depths up to ∼ 170 h ⁻¹ Mpc and (b) a large-scale coherence of the dipole anisotropy, indications provided mainly by the fact that the differential dipoles of large equal-volume shells are aligned with the CMB dipole and exhibit significant dipole signals. The two IRAS dipoles are indistinguishable within 50 h ⁻¹ Mpc and beyond ∼ 130 h ⁻¹ Mpc while the QDOT dipole, having a lower flux limit, continues growing with respect to the 1.2-Jy sample up to ∼ 100 h ⁻¹ Mpc in agreement with Rowan-Robinson et al.     

Peculiar velocity reconstruction with the fast action method: tests on mock redshift surveys

We present extensive tests of the fast action method (FAM) for recovering the past orbits of mass tracers in an expanding universe from their redshift-space coordinates at the present epoch. The tests focus on the reconstruction of present-day peculiar velocities using mock catalogues extracted from high-resolution N -body simulations. The method allows for a self-consistent treatment of redshift-space distortions by direct minimization of a modified action for a cosmological gravitating system. When applied to ideal, volume-limited catalogues, FAM recovers unbiased peculiar velocities with a one-dimensional, 1σ error of ∼220 km  s⁻¹  , if velocities are smoothed on a scale of  5 h ⁻¹  Mpc. Alternatively, when no smoothing is applied, FAM predicts nearly unbiased velocities for objects residing outside the highest density regions. In this second case the 1σ error decreases to a level of ∼150 km  s⁻¹  . The correlation properties of the peculiar velocity fields are also correctly recovered on scales larger than  5 h ⁻¹  Mpc. Similar results are obtained when FAM is applied to flux-limited catalogues mimicking the IRAS PSC z survey. In this case FAM reconstructs peculiar velocities with similar intrinsic random errors, while velocity–velocity correlation properties are well reproduced beyond scales of  ∼8 h ⁻¹  Mpc. We also show that FAM provides better velocity predictions than other, competing methods based on linear theory or the Zel'dovich approximation. These results indicate that FAM can be successfully applied to presently available galaxy redshift surveys such as IRAS PSC z .     

Full-sky maps for gravitational lensing of the cosmic microwave background

We use the large cosmological Millennium Simulation (MS) to construct the first all-sky maps of the lensing potential and the angle, aiming at gravitational lensing of the cosmic microwave background (CMB), with the goal of properly including small-scale non-linearities and non-Gaussianity. Exploiting the Born approximation, we implement a map-making procedure based on direct ray tracing through the gravitational potential of the MS. We stack the simulation box in redshift shells up to z ∼ 11, producing continuous all-sky maps with arcmin angular resolution. A randomization scheme avoids the repetition of structures along the line of sight, and structures larger than the MS box size are added to supply the missing contribution of large-scale (LS) structures to the lensing signal. The angular power spectra of the projected lensing potential and the deflection-angle modulus agree quite well with semi-analytic estimates on scales down to a few arcmin, while we find a slight excess of power on small scales, which we interpret as being due to non-linear clustering in the MS. Our map-making procedure, combined with the LS adding technique, is ideally suited for studying lensing of CMB anisotropies, for analysing cross-correlations with foreground structures, or other secondary CMB anisotropies such as the Rees–Sciama effect.     

The tidal disruption rate in dense galactic cusps containing a supermassive binary black hole

We consider the problem of tidal disruption of stars in the centre of a galaxy containing a supermassive binary black hole with unequal masses. We assume that over the separation distance between the black holes, the gravitational potential is dominated by the more massive primary black hole. Also, we assume that the number density of stars is concentric with the primary black hole and has a power-law cusp. We show that the bulk of stars with a small angular-momentum component normal to the black hole binary orbit can reach a small value of total angular momentum through secular evolution in the gravitational field of the binary, and hence they can be tidally disrupted by the larger black hole. This effect is analogous to the so-called Kozai effect well known in celestial mechanics. We develop an analytical theory for the secular evolution of the stellar orbits and calculate the rate of tidal disruption. We compare our analytical theory with a simple numerical model and find very good agreement.
Our results show that for a primary black hole mass of  ∼10⁶–10⁷ M_⊙  , the black hole mass-ratio   q > 10⁻²  , cusp size ∼1 pc, the tidal disruption rate can be as large as  ∼10⁻²–1 M_⊙ yr⁻¹  . This is at least 10²–10⁴ times larger than estimated for the case of a single supermassive black hole. The duration of the phase of enhanced tidal disruption is determined by the dynamical-friction time-scale, and it is rather short: ∼10⁵ yr. The dependence of the tidal disruption rate on the mass ratio, and on the size of the cusp, is also discussed.     

A measurement of the faint source correlation function in the GOODS and UDF surveys

We present a stable procedure for defining and measuring the two point angular autocorrelation function,   w (θ) =[θ/θ₀( V )]^−Γ  , of faint  (25 < V < 29)  , barely resolved and unresolved sources in the Hubble Space Telescope Great Observatories Origins Deep Survey and Ultra Deep Field data sets. We construct catalogues that include close pairs and faint detections. We show, for the first time, that, on subarcsec scales, the correlation function exceeds unity. This correlation function is well fit by a power law with index  Γ≈ 2.5  and a  θ₀= 10^{−0.1( V −25.8)} arcsec  . This is very different from the values of  Γ≈ 0.7  and  θ₀( r ) = 10^{−0.4( r −21.5)} arcsec  associated with the gravitational clustering of brighter galaxies. This observed clustering probably reflects the presence of giant star-forming regions within galactic-scale potential wells. Its measurement enables a new approach to measuring the redshift distribution of the faintest sources in the sky.     

Testing tidal-torque theory – I. Spin amplitude and direction

We evaluate the success of linear tidal-torque theory (TTT) in predicting galactic-halo spin using a cosmological N -body simulation with thousands of well-resolved haloes. The protohaloes are identified by tracing today's haloes back to the initial conditions. The TTT predictions for the protohaloes match, on average, the spin amplitudes of the virialized haloes of today, if linear growth is assumed until ∼ t ₀/3, or  55–70  per cent of the halo effective turn-around time. This makes it a useful qualitative tool for understanding certain average properties of galaxies, such as total spin and angular momentum distribution within haloes, but with a random scatter of the order of the signal itself. Non-linear changes in spin direction cause a mean error of ∼50° in the TTT prediction at t ₀, such that the linear spatial correlations of spins on scales ≥1  h ⁻¹ Mpc are significantly weakened by non-linear effects. This questions the usefulness of TTT for predicting intrinsic alignments in the context of gravitational lensing. We find that the standard approximations made in TTT, including a second-order expansion of the Zel'dovich potential and a smoothing of the tidal field, provide close-to-optimal results.     

The supergalactic plane revisited with the Optical Redshift Survey

We re-examine the existence and extent of the planar structure in the local galaxy density field, the so-called supergalactic plane (SGP). This structure is studied here in three dimensions using both the new Optical Redshift Survey (ORS) and the IRAS 1.2-Jy redshift survey. The density contrast in a slab of thickness 20  h ⁻¹ Mpc and diameter 80 Mpc aligned with the standard de Vaucouleurs supergalactic coordinates is δ _sgp∼0.5 for both ORS and IRAS . The structure of the SGP is not well described by a homogeneous ellipsoid, although it does appear to be a flattened structure, which we quantify by calculating the moment of inertia tensor of the density field. The directions of the principal axes vary with radius, but the minor axis remains within θ _z ∼30° of the standard SGP Z -axis, out to a radius of 80  h ⁻¹ Mpc, for both ORS and IRAS . However, the structure changes shape with radius, varying between a flattened pancake and a dumbbell, the latter at a radius of ∼50  h ⁻¹ Mpc, where the Great Attractor and Perseus–Pisces superclusters dominate the distribution. This calls to question the connectivity of the 'plane' beyond ∼40  h ⁻¹ Mpc. The configuration found here can be viewed as part of a web of filaments and sheets, rather than as an isolated pancake-like structure. An optimal minimum variance reconstruction of the density field using Wiener filtering, which corrects for both redshift distortion and shot noise, yields a similar misalignment angle and behaviour of axes. The background-independent statistic of axes proposed here can be best used for testing cosmological models by comparison with N -body simulations.     

The effects of anticorrelation on gravitational clustering

We use non-linear scaling relations (NSRs) to investigate the effects arising from the existence of negative correlations on the evolution of gravitational clustering in an expanding universe. It turns out that such anticorrelated regions have important dynamical effects on all scales. In particular, the mere existence of negative values for the linear two-point correlation function ξ¯ _L over some range of scales starting from l = L ₀ implies that the non-linear correlation function is bounded from above at all scales x < L ₀ . This also results in the relation ξ¯   ∝  x ⁻³ , at these scales, at late times, independent of the original form of the correlation function. Current observations do not rule out the existence of negative ξ¯ for 200  h ⁻¹ Mpc≲ ξ¯ ≲1000  h ⁻¹ Mpc; the present work may thus have relevance for the real Universe. The only assumption made in the analysis is the existence of NSR; the results are independent of the form of the NSR as well as of the stable clustering hypothesis.     

Episodic accretion in magnetically layered protoplanetary discs

We study protoplanetary disc evolution assuming that angular momentum transport is driven by gravitational instability at large radii, and magnetohydrodynamic (MHD) turbulence in the hot inner regions. At radii of the order of 1 au such discs develop a magnetically layered structure, with accretion occurring in an ionized surface layer overlying quiescent gas that is too cool to sustain MHD turbulence. We show that layered discs are subject to a limit cycle instability, in which accretion on to the protostar occurs in ∼10⁴-yr bursts with ̇ ∼10⁻⁵ M_⊙ yr⁻¹, separated by quiescent intervals lasting ∼10⁵ yr where ̇ ≈10⁻⁸ M_⊙ yr⁻¹. Such bursts could lead to repeated episodes of strong mass outflow in young stellar objects. The transition to this episodic mode of accretion occurs at an early epoch ( t ≪1 Myr), and the model therefore predicts that many young pre-main-sequence stars should have low rates of accretion through the inner disc. At ages of a few Myr, the discs are up to an order of magnitude more massive than the minimum-mass solar nebula, with most of the mass locked up in the quiescent layer of the disc at r ∼1 au. The predicted rate of low-mass planetary migration is reduced at the outer edge of the layered disc, which could lead to an enhanced probability of giant planet formation at radii of 1–3 au.     

The Ω dependence in equations of motion

The equations of motion governing the evolution of a collisionless gravitating system of particles in an expanding universe can be cast in a form which is almost independent of the cosmological density parameter, Ω, and the cosmological constant, Λ. The new equations are expressed in terms of a time variable τ=ln D , where D is the linear rate of growth of density fluctuations. The dependence on the density parameter is proportional to ε=Ω^−0.2−1 times the difference between the peculiar velocity (with respect to τ) of particles and the gravity field (minus the gradient of the potential); or, before shell-crossing, times the sum of the density contrast and the velocity divergence. In a one-dimensional collapse or expansion, the equations are fully independent of Ω and Λ before shell crossing. In the general case, the effect of this weak Ω dependence is to enhance the rate of evolution of density perturbations in dense regions. In a flat universe with Λ7ne;0, this enhancement is less pronounced than in an open universe with Λ=0 and the same Ω. Using the spherical collapse model, we find that the increase of the rms density fluctuations in a low-Ω universe relative to that in a flat universe with the same linear normalization is ∼0.01ε(Ω)〈δ³〉, where δ is the density field in the flat universe. The equations predict that the smooth average velocity field scales like Ω^0.6, while the local velocity dispersion (rms value) scales, approximately, like Ω^0.5. High-resolution N -body simulations confirm these results and show that density fields, when smoothed on scales slightly larger than clusters, are insensitive to the cosmological model. Haloes in an open model simulation are more concentrated than haloes of the same M /Ω in a flat model simulation.     

Seed magnetic Fields Generated by Primordial Supernova Explosions

The origin of the magnetic field in galaxies is an open question in astrophysics. Several mechanisms have been proposed related, in general, to the generation of small seed fields amplified by a dynamo mechanism. In general, these mechanisms have difficulty in satisfying both the requirements of a sufficiently high strength for the magnetic field and the necessary large coherent scales. We show that the formation of dense and turbulent shells of matter, in the multiple explosion scenario of Miranda &38; Opher for the formation of the large-scale structures of the Universe, can naturally act as a seed for the generation of a magnetic field. During the collapse and explosion of Population III objects, a temperature gradient not parallel to a density gradient can naturally be established, producing a seed magnetic field through the Biermann battery mechanism. We show that seed magnetic fields ∼ 10⁻¹²–10⁻¹⁴ G can be produced in this multiple explosion scenario on scales of the order of clusters of galaxies (with coherence length L  ∼ 1.8 Mpc) and up to ∼ 4.5 × 10⁻¹⁰ G on scales of galaxies ( L  ∼ 100 kpc).     

Reconstruction of cosmological initial conditions from galaxy redshift catalogues

We present and test a new method for the reconstruction of cosmological initial conditions from a full-sky galaxy catalogue. This method, called ZTRACE, is based on a self-consistent solution of the growing mode of gravitational instabilities according to the Zel'dovich approximation and higher order in Lagrangian perturbation theory. Given the evolved redshift-space density field, smoothed on some scale, ZTRACE finds, via an iterative procedure, an approximation to the initial density field for any given set of cosmological parameters; real-space densities and peculiar velocities are also reconstructed. The method is tested by applying it to N -body simulations of an Einstein–de Sitter and an open cold dark matter universe. It is shown that errors in the estimate of the density contrast dominate the noise of the reconstruction. As a consequence, the reconstruction of real-space density and peculiar velocity fields using non-linear algorithms is little improved over those based on linear theory. The use of a mass-preserving adaptive smoothing, equivalent to a smoothing in Lagrangian space, allows an unbiased (although noisy) reconstruction of initial conditions, as long as the (linearly extrapolated) density contrast does not exceed unity. The probability distribution function of the initial conditions is recovered to high precision, even for Gaussian smoothing scales of ∼5  h ⁻¹ Mpc, except for the tail at δ ≥1. This result is insensitive to the assumptions of the background cosmology.     

Sunyaev–Zel'dovich distortion from early galactic winds

We consider the distortion in the cosmic microwave background (CMB) resulting from galactic winds at high redshift. Winds outflowing from galaxies have been hypothesized to be possible sources of metals in the intergalactic medium, which is known to have been enriched to 10^−2.5 Z_⊙ at z ∼3. We model these winds as functions of mass of the parent galaxy and redshift, assuming that they activate at a common initial redshift, z _in, and calculate the mean y -distortion and the angular power spectrum of the distortion in the CMB. We find that the thermal Sunyaev–Zel'dovich (SZ) effect resulting from the winds is consistent with previous estimates. The distortion arising from the kinetic SZ (kSZ) effect is, however, found to be more important than the thermal SZ (tSZ) effect. We find that the distortion resulting from galactic winds is an important contribution to the power spectrum of distortion at very small angular scales ( l ∼10⁴). We also find that the power spectrum resulting from clustering dominates the Poisson power spectrum for l ≤(4–5)×10⁵. We show explicitly how the combined power spectrum from wind dominates over that of clusters at 217 GHz, relevant for PLANCK . We also show how these constraints change when the efficiency of the winds is varied.     

Gravitational wave signals from the collapse of the first stars

We study the gravitational wave emission from the first stars, which are assumed to be very massive objects (VMOs). We take into account various feedback (both radiative and stellar) effects regulating the collapse of objects in the early Universe and thus derive the VMO initial mass function and formation rate. If the final fate of VMOs is to collapse, leaving very massive black hole remnants, then the gravitational waves emitted during each collapse would be seen as a stochastic background. The predicted spectral strain amplitude in a critical density cold dark matter (CDM) universe peaks in the frequency range ν ≈5×10⁻⁴–5×10⁻³ Hz, where it has a value in the range ≈10⁻²⁰–10⁻¹⁹ Hz^−1/2, and might be detected by the Laser Interferometer Space Antenna ( LISA ). The expected emission rate is roughly 4000 event yr⁻¹, resulting in a stationary discrete sequence of bursts, i.e. a shot-noise signal.     

Neutralino annihilations and the gas temperature in the dark ages

Assuming that the dark matter is entirely made up of neutralinos, we re-visit the role of their annihilation on the temperature of diffuse gas in the high-redshift universe  ( z > 10)  , before the formation of luminous structures. We consider neutralinos of particle mass 36 and 100 GeV. The former is able to produce  ∼7  e ⁻ e ⁺  particles per annihilation through the fremionic channel, and the latter ∼53 particles assuming a purely bosonic channel. High-energy   e ⁻ e ⁺  particles up-scatter the cosmic microwave background (CMB) photons into higher energies via the inverse-Compton scattering. The process produces a power-law   e ⁻ e ⁺  energy spectrum of index −1 in the energy range of interest, independent of the initial energy distribution. The corresponding energy spectrum of the up-scattered photons is a power law of index −1/2, if absorption by the gas is not included. The scattered photons photoheat the gas by releasing electrons which deposit a fraction (14 per cent) of their energy as heat into the ambient medium. For uniformly distributed neutralinos, the heating is insignificant. The effect is greatly enhanced by the clumping of neutralinos into dense haloes. We use a time-dependent clumping model which takes into account the damping of density fluctuations on mass-scales smaller than  ∼10⁻⁶ M_⊙  . With this clumping model, the heating mechanism boosts the gas temperature above that of the CMB after a redshift of   z ∼ 30  . By   z ≈ 10  , the gas temperature is nearly 100 times its temperature when no heating is invoked. Similar increase is obtained for the two neutralino masses considered.     

Bianchi model CMB polarization and its implications for CMB anomalies

We derive the cosmic microwave background (CMB) radiative transfer equation in the form of a multipole hierarchy in the nearly Friedmann–Robertson–Walker limit of homogeneous, but anisotropic, universes classified via their Bianchi type. Compared with previous calculations, this allows a more sophisticated treatment of recombination, produces predictions for the polarization of the radiation and allows for reionization. Our derivation is independent of any assumptions about the dynamical behaviour of the field equations, except that it requires anisotropies to be small back to recombination; this is already demanded by observations.
We calculate the polarization signal in the Bianchi VII _h case, with the parameters recently advocated to mimic the several large-angle anomalous features observed in the CMB. We find that the peak polarization signal is  ∼1.2 μK  for the best-fitting model to the temperature anisotropies, and is mostly confined to multipoles   l < 10  . Remarkably, the predicted large-angle EE and TE power spectra in the Bianchi model are consistent with Wilkinson Microwave Anisotropy Probe ( WMAP ) observations that are usually interpreted as evidence of early reionization. However, the power in B-mode polarization is predicted to be similar to the E-mode power and parity-violating correlations are also predicted by the model; the WMAP non-detection of either of these signals casts further strong doubts on the veracity of attempts to explain the large-angle anomalies with global anisotropy. On the other hand, given that there exist further dynamical degrees of freedom in the VII _h universes that are yet to be compared with CMB observations, we cannot at this time definitively reject the anisotropy explanation.     

The evolution of galactic discs

We use recent observations of high-redshift galaxies to study the evolution of galactic discs over the redshift range 0 <  z ≲1. The data are inconsistent with models in which discs were already assembled at z  = 1 and have evolved only in luminosity since that time. Assuming that disc properties change with redshift as powers of 1 +   z and analysing the observations assuming an Einstein–de Sitter universe, we find that for given rotation speed, disc scalelength decreases with z as ∼ (1 +  z )⁻¹, total B -band mass-to-light ratio decreases with z as ∼ (1 +  z )⁻¹, and disc luminosity (again in B ) depends only weakly on z . These scalings are consistent with current data on the evolution of disc galaxy abundance as a function of size and luminosity. Both the scalings and the abundance evolution are close to the predictions of hierarchical models for galaxy formation. If different cosmogonies are compared, the observed evolution in disc size and disc abundance favours a flat low-Ω₀ universe over an Einstein–de Sitter universe.     

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.     

The three-point function in large-scale structure – I. The weakly non-linear regime in N-body simulations

This paper presents a comparison of the predictions for the two- and three-point correlation functions of density fluctuations, ξ and ζ , in gravitational perturbation theory (PT) against large cold dark matter (CDM) simulations. This comparison is made possible for the first time on large weakly non-linear scales (>10  h ⁻¹ Mpc) thanks to the development of a new algorithm for estimating correlation functions for millions of points in only a few minutes. Previous studies in the literature comparing the PT predictions of the three-point statistics with simulations have focused mostly on Fourier space, angular space or smoothed fields. Results in configuration space, such as those presented here, were limited to small scales where leading-order PT gives a poor approximation. Here we also propose and apply a method for separating the first-order and subsequent contributions to PT by combining different output times from the evolved simulations. We find that in all cases there is a regime where simulations do reproduce the leading-order (tree-level) predictions of PT for the reduced three-point function   Q ₃∼ ζ / ξ ²  . For steeply decreasing correlations (such as the standard CDM model) deviations from the tree-level results are important even at relatively large scales, ≃20 Mpc  h ⁻¹. On larger scales ξ goes to zero and the results are dominated by sampling errors. In more realistic models (such as the ΛCDM cosmology) deviations from the leading-order PT become important at smaller scales   r ≃10 Mpc  h ^-1  , although this depends on the particular three-point configuration. We characterize the range of validity of this agreement and show the behaviour of the next-order (one-loop) corrections.