The clustering of barred galaxies in the local Universe

We study the clustering properties of barred galaxies using data from the Sloan Digital Sky Survey. We compute projected redshift-space two-point cross-correlation functions   w _p( r _p)  for a sample of nearly 1000 galaxies for which we have performed detailed structural decompositions using the methods described in Gadotti. The sample includes 286 barred galaxies. The clustering of barred and unbarred galaxies of similar stellar mass is indistinguishable over all the scales probed (∼20 kpc–30 Mpc). This result also holds even if the sample is restricted to bars with bluer   g − i   colours (and hence younger ages). Our result also does not change if we split our sample of barred galaxies according to bar-to-total luminosity ratio, bar boxyness, effective surface brightness, length or the shape of the surface density profile within the bar. There is a hint that red, elliptical bars are more strongly clustered than red and less elliptical bars, on scales  ≳1 Mpc  , although the statistical significance is not high. We conclude that there is no significant evidence that bars are a product of mergers or interactions. We tentatively interpret the stronger clustering of the more elliptical bars as evidence that they are located in older galaxies, which reside in more massive haloes.     

Interaction-induced star formation in a complete sample of 105 nearby star-forming galaxies

We investigate the clustering properties of a complete sample of 10⁵ star-forming galaxies drawn from the data release 4 (DR4) of the Sloan Digital Sky Survey. On scales less than 100 kpc, the amplitude of the correlation function exhibits a strong dependence on the specific star formation rate (SSFR) of the galaxy. We interpret this as the signature of enhanced star formation induced by tidal interactions. We then explore how the average star formation rate (SFR) in a galaxy is enhanced as the projected separation r _p between the galaxy and its companions decreases. We find that the enhancement strongly depends on r _p, but very weakly on the relative luminosity of the companions. The enhancement is also stronger in low-mass galaxies than in high-mass galaxies. In order to explore whether a tidal interaction is not only sufficient, but also necessary to trigger enhanced star formation in a galaxy, we compute background subtracted neighbour counts for the galaxies in our sample. The average number of close neighbours around galaxies with low to average values of SFR/ M _* is close to zero. At the highest SSFRs, however, more than 40 per cent of the galaxies in our sample have a companion within a projected radius of 100 kpc. Visual inspection of the highest SFR/ M _* galaxies without companions reveals that more than 50 per cent of these are clear interacting or merging systems. We conclude that tidal interactions are the dominant trigger of enhanced star formation in the most strongly star-forming systems. Finally, we find clear evidence that tidal interactions not only lead to enhanced star formation in galaxies, but also cause structural changes such as an increase in concentration.     

The one-point PDF of the initial conditions of our local Universe from the IRAS PSC redshift catalogue

The algorithm ztrace of Monaco & Efstathiou is applied to the IRAS PSCz catalogue to reconstruct the initial conditions of our local Universe with a resolution down to ~5  h ¹ Mpc. The one-point probability distribution function (PDF) of the reconstructed initial conditions is consistent with the assumptions that: (i) IRAS galaxies trace mass on scales of ~5  h ¹ Mpc and (ii) the statistics of the primordial density fluctuations are Gaussian. We use simulated PSCz catalogues, constructed from N -body simulations with Gaussian initial conditions, to show that local non-linear bias can cause the recovered initial PDF (assuming no bias) to be non-Gaussian. However, for plausible bias models, the distortions of the recovered PDF would be difficult to detect using the volume finely sampled by the PSCz catalogue. So, for Gaussian initial conditions, a range of bias models remain compatible with our PSCz reconstruction results.     

Non-linear redshift distortions: the two-point correlation function

We consider a situation where the density and peculiar velocities in real space are linear, and we calculate ξ _s , the two-point correlation function in redshift space, incorporating all non-linear effects which arise as a consequence of the map from real to redshift space. Our result is non-perturbative and it includes the effects of possible multi-streaming in redshift space. We find that the deviations from the predictions of the linear redshift distortion analysis increase for the higher spherical harmonics of ξ _s . While the deviations are insignificant for the monopole ξ ₀, the hexadecapole ξ ₄ exhibits large deviations from the linear predictions. For a COBE normalized     ,     cold dark matter (CDM) power spectrum, our results for ξ ₄ deviate from the linear predictions by a factor of two on the scale of ∼10  h ⁻¹ Mpc. The deviations from the linear predictions depend separately on f (Ω) and b . This holds the possibility of removing the degeneracy that exists between these two parameters in the linear analysis of redshift surveys which yields only     .
We also show that the commonly used phenomenological model, where the non-linear redshift two-point correlation function is calculated by convolving the linear redshift correlation function with an isotropic pair velocity distribution function, is a limiting case of our result.     

The limits of bound structures in the accelerating Universe

According to the latest evidence, the Universe is entering an era of exponential expansion, where gravitationally bound structures will get disconnected from each other, forming isolated 'island universes'. In this scenario, we present a theoretical criterion to determine the boundaries of gravitationally bound structures and a physically motivated definition of superclusters as the largest bound structures in the Universe. We use the spherical collapse model self-consistently to obtain an analytical condition for the mean density enclosed by the last bound shell of the structure (2.36 times the critical density in the present Universe, assumed to be flat, with 30 per cent matter and 70 per cent cosmological constant, in agreement with the previous, numerical result of Chiueh & He). N -body simulations extended to the future show that this criterion, applied at the present cosmological epoch, defines a sphere that encloses ≈99.7 per cent of the particles that will remain bound to the structure at least until the scale parameter of the Universe is 100 times its present value. On the other hand, (28 ± 13) per cent of the enclosed particles are in fact not bound, so the enclosed mass overestimates the bound mass, in contrast with the previous, less rigorous criterion of, e.g. Busha and collaborators, which gave a more precise mass estimate. We also verify that the spherical collapse model estimate for the radial infall velocity of a shell enclosing a given mean density gives an accurate prediction for the velocity profile of infalling particles, down to very near the centre of the virialized core.     

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.     

Recovering the initial conditions of our local Universe from NOG and PSCz catalogues

We apply the ztrace algorithm to the optical NOG and infrared PSC z galaxy catalogues to reconstruct the pattern of primordial fluctuations that have generated our local Universe. We check that the density fields traced by the two catalogues are well correlated, and consistent with a linear relation [either in δ or in  log (1 +δ)  ] with relative bias (of NOG with respect to PSC z )   b _rel= 1.1 ± 0.1  . The relative bias relation is used to fill the optical zone of avoidance at  | b | < 20°  using the PSC z galaxy density field.
We perform extensive testing on simulated galaxy catalogues to optimize the reconstruction. The quality of the reconstruction is predicted to be good at large scales, up to a limiting wavenumber   k _lim≃ 0.4 h Mpc⁻¹  beyond which all information is lost. We find that the improvement arising from the denser sampling of the optical catalogue is compensated by the uncertainties connected to the larger zone of avoidance.
The initial conditions reconstructed from the NOG catalogue are found (analogously to those from the PSC z ) to be consistent with a Gaussian paradigm. We use the reconstructions to produce sets of initial conditions ready to be used for constrained simulations of our local Universe.     

Interactions, star formation and AGN activity

It has long been known that galaxy interactions are associated with enhanced star formation. In a companion paper, we explored this connection by applying a variety of statistics to Sloan Digital Sky Survey (SDSS) data. In particular, we showed that specific star formation rates of galaxies are higher if they have close neighbours. Here, we apply exactly the same techniques to active galactic nuclei (AGN) in the survey, showing that close neighbours are not associated with any similar enhancement of nuclear activity. Star formation is enhanced in AGN with close neighbours in exactly the same way as in inactive galaxies, but the accretion rate on to the black hole, as estimated from the extinction-corrected [O  iii ] luminosity, is not influenced by the presence or absence of companions. Previous work has shown that galaxies with more strongly accreting black holes contain more young stars in their inner regions. This leads us to conclude that star formation induced by a close companion and star formation associated with black hole accretion are distinct events. These events may be part of the same physical process, for example a merger, provided they are separated in time. In this case, accretion on to the black hole and its associated star formation would occur only after the two interacting galaxies have merged. The major caveat in this work is our assumption that the extinction-corrected [O  iii ] luminosity is a robust indicator of the bolometric luminosity of the central black hole. It is thus important to check our results using indicators of AGN activity at other wavelengths.     

The cluster distribution as a test of dark matter models — IV. Topology and geometry

We study the geometry and topology of the large-scale structure traced by galaxy clusters in numerical simulations of a box of side 320 h ⁻¹ Mpc, and compare them with available data on real clusters. The simulations we use are generated by the Zel'dovich approximation, using the same methods as we have used in the first three papers in this series. We consider the following models to see if there are measurable differences in the topology and geometry of the superclustering they produce: (i) the standard cold dark matter model (SCDM); (ii) a CDM model with Ω₀ = 0.2 (OCDM); (iii) a CDM model with a 'tilted' power spectrum having n  = 0.7 (TCDM); (iv) a CDM model with a very low Hubble constant, h  = 0.3 (LOWH); (v) a model with mixed CDM and HDM (CHDM); (vi) a flat low-density CDM model with Ω₀ = 0.2 and a non-zero cosmological Λ term (ΛCDM). We analyse these models using a variety of statistical tests based on the analysis of: (i) the Euler–Poincaré characteristic; (ii) percolation properties; (iii) the minimal spanning tree construction. Taking all these tests together we find that the best-fitting model is ΛCDM and, indeed, the others do not appear to be consistent with the data. Our results demonstrate that despite their biased and extremely sparse sampling of the cosmological density field, it is possible to use clusters to probe subtle statistical diagnostics of models, which go far beyond the low-order correlation functions usually applied to study superclustering.