Large multi-object spectroscopic surveys require automated algorithms to optimise their observing strategy. One of the most ambitious upcoming spectroscopic surveys is the 4MOST survey. The 4MOST survey facility is a fibre-fed spectroscopic instrument on the VISTA telescope with a large enough field of view to survey a large fraction of the southern sky within a few years. Several Galactic and extragalactic surveys will be carried out simultaneously, so the combined target density will strongly vary. In this paper, we describe a new tiling algorithm that can naturally deal with the large target density variations on the sky and which automatically handles the different exposure times of targets. The tiling pattern is modelled as a marked point process, which is characterised by a probability density that integrates the requirements imposed by the 4MOST survey. The optimal tilling pattern with respect to the defined model is estimated by the tiles configuration that maximises the proposed probability density. In order to achieve this maximisation a simulated annealing algorithm is implemented. The algorithm automatically finds an optimal tiling pattern and assigns a tentative sky brightness condition and exposure time for each tile, while minimising the total execution time that is needed to observe the list of targets in the combined input catalogue of all surveys. Hence, the algorithm maximises the long-term observing efficiency and provides an optimal tiling solution for the survey. While designed for the 4MOST survey, the algorithm is flexible and can with simple modifications be applied to any other multi-object spectroscopic survey.

Torsion is a non-Riemannian geometrical extension of general relativity that allows including the spin of matter and the twisting of spacetime. Cosmological models with torsion have been considered in the literature to solve problems of either the very early (high redshift $z$) or the present-day Universe. This paper focuses on distinguishable observational signatures of torsion that could not be otherwise explained with a scalar field in pseudo-Riemannian geometry. We show that when torsion is present, the cosmic duality relation between the angular diameter distance, $D_{\mathrm A}$, and the luminosity distance, $D_{\mathrm L}$, is broken. We show how the deviation described by the parameter $\eta = D_{\mathrm L}/[ D_{\mathrm A}(1+z)^2] -1 $ is linked to torsion and how different forms of torsion lead to special-case parametrisations of $\eta$, including $\eta_0 z$, $\eta_0 z/(1+z)$, and $\eta_0 \ln (1+z)$. We also show that the effects of torsion could be visible in low-redshift data, inducing biases in supernovae-based $H_0$ measurements. We also show that torsion can impact the Clarkson-Bassett-Lu (CBL) function ${\cal C}(z) = 1 + H^2 (D D'' - D'^2) + H H' D D'$, where $D$ is the transverse comoving distance. If $D$ is inferred from the luminosity distance, then, in general non-zero torsion models, ${\cal C}(z) \ne 0$. For pseudo-Riemannian geometry, the Friedmann-Lemaitre-Robertson-Walker (FLRW) metric has ${\cal C}(z) \equiv 0$; thus, measurement of the CBL function could provide another diagnostic of torsion.

Context. Several new multi-object spectrographs are currently planned or under construction that are capable of observing thousands of Galactic and extragalactic objects simultaneously. Aims. In this paper we present a probabilistic fibre-to-target assignment algorithm that takes spectrograph targeting constraints into account and is capable of dealing with multiple concurrent surveys. We present this algorithm using the 4-metre Multi-Object Spectroscopic Telescope (4MOST) as an example. Methods. The key idea of the proposed algorithm is to assign probabilities to fibre-target pairs. The assignment of probabilities takes the fibre positioner's capabilities and constraints into account. Additionally, these probabilities include requirements from surveys and take the required exposure time, number density variation, and angular clustering of targets across each survey into account. The main advantage of a probabilistic approach is that it allows for accurate and easy computation of the target selection function for the different surveys, which involves determining the probability of observing a target, given an input catalogue. Results. The probabilistic fibre-to-target assignment allows us to achieve maximally uniform completeness within a single field of view. The proposed algorithm maximises the fraction of successfully observed targets whilst minimising the selection bias as a function of exposure time. In the case of several concurrent surveys, the algorithm maximally satisfies the scientific requirements of each survey and no specific survey is penalised or prioritised. Conclusions. The algorithm presented is a proposed solution for the 4MOST project that allows for an unbiased targeting of many simultaneous surveys. With some modifications, the algorithm may also be applied to other multi-object spectroscopic surveys.

Cosmological observational analysis frequently assumes that the Universe is spatially flat. We aim to non-perturbatively check the conditions under which a flat or nearly flat expanding dust universe, including the LCDM model if interpreted as strictly flat, forbids the gravitational collapse of structure. We quantify spatial curvature at turnaround. We use the Hamiltonian constraint to determine the pointwise conditions required for an overdensity to reach its turnaround epoch in an exactly flat spatial domain. We illustrate this with a plane-symmetric, exact, cosmological solution of the Einstein equation, extending earlier work. More generally, for a standard initial power spectrum, we use the relativistic Zel'dovich approximation implemented in 'inhomog' to numerically estimate how much positive spatial curvature is required for turnaround to be allowed at typical epochs/length scales in almost-EdS and almost-LCDM models with inhomogeneous curvature. We find that gravitational collapse in a spatially exactly flat, irrotational, expanding, dust universe is relativistically forbidden pointwise. In the spatially flat plane-symmetric model considered here, pancake collapse is excluded both pointwise and in averaged domains. In an almost-EdS or LCDM model, the per-domain average curvature in collapsing domains almost always becomes strongly positive prior to turnaround, with the expansion-normalised curvature functional reaching $\Omega_{\cal R}^{\cal D} \sim -5$. We show analytically that a special case gives $\Omega_{\cal R}^{\cal D} = -5$ exactly (if normalised using the EdS expansion rate) at turnaround. An interpretation of LCDM as literally 3-Ricci flat would forbid structure formation. The difference between relativistic cosmology and a strictly flat LCDM model is fundamental in principle, but we find that the geometrical effect is weak.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

The 4MOST Cosmology Redshift Survey (CRS) will perform stringent cosmological tests via spectroscopic clustering measurements that will complement the best lensing, cosmic microwave background and other surveys in the southern hemisphere. The combination of carefully selected samples of bright galaxies, luminous red galaxies, emission-line galaxies and quasars, totalling about 8 million objects over the redshift range $z = 0.15$ to $3.5$, will allow definitive tests of gravitational physics. Many key science questions will be addressed by combining CRS spectra of these targets with data from current or future facilities such as the Large Synoptic Survey Telescope, the Square Kilometre Array and the Euclid mission.

Standard cosmological $N$-body simulations have background scale factor evolution that is decoupled from non-linear structure formation. Prior to gravitational collapse, kinematical backreaction ($Q_D$) justifies this approach in a Newtonian context. However, the final stages of a gravitational collapse event are sudden; a globally imposed expansion rate thus forces at least one expanding region to suddenly decelerate. This is relativistically unrealistic. Instead, we allow non-collapsed domains to evolve in volume according to the $Q_D$ Zel'dovich Approximation (QZA). We study the inferred average expansion under this "silent" virialisation hypothesis. We set standard (mpgrafic) EdS cosmological $N$-body initial conditions. Using RAMSES, we call DTFE to estimate the initial values of the three invariants of the extrinsic curvature tensor in Lagrangian domains $D$. We integrate the Raychaudhuri equation in each domain using inhomog, adopt the stable clustering hypothesis (VQZA), and average spatially. We adopt an early-epoch--normalised EdS reference-model Hubble constant $H_1^{bg} = 37.7$ km/s/Mpc and an effective Hubble constant $H_0^{eff} = 67.7$ km/s/Mpc. From 2000 simulations at resolution $256^3$, a unity effective scale factor is reached at 13.8~Gyr (16% above EdS) for an averaging scale of $L_{13.8}=2.5^{+0.1}_{-0.4}$ Mpc/$h^{eff}$. Relativistically interpreted, this corresponds to strong average negative curvature evolution. The virialisation fraction and super-EdS expansion correlate strongly at fixed cosmological time. Thus, starting from EdS initial conditions and averaging on a typical non-linear structure formation scale, the VQZA dark-energy--free average expansion matches $\Lambda$CDM expansion to first order. The software packages used here are free-licensed.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

We discuss the potential of a next generation space-borne Cosmic Microwave Background (CMB) experiment for studies of extragalactic sources. Our analysis has particular bearing on the definition of the future space project, CORE, that has been submitted in response to ESA's call for a Medium-size mission opportunity as the successor of the Planck satellite. Even though the effective telescope size will be somewhat smaller than that of Planck, CORE will have a considerably better angular resolution at its highest frequencies, since, in contrast with Planck, it will be diffraction limited at all frequencies. The improved resolution implies a considerable decrease of the source confusion, i.e. substantially fainter detection limits. In particular, CORE will detect thousands of strongly lensed high-z galaxies distributed over the full sky. The extreme brightness of these galaxies will make it possible to study them, via follow-up observations, in extraordinary detail. Also, the CORE resolution matches the typical sizes of high-z galaxy proto-clusters much better than the Planck resolution, resulting in a much higher detection efficiency; these objects will be caught in an evolutionary phase beyond the reach of surveys in other wavebands. Furthermore, CORE will provide unique information on the evolution of the star formation in virialized groups and clusters of galaxies up to the highest possible redshifts. Finally, thanks to its very high sensitivity, CORE will detect the polarized emission of thousands of radio sources and, for the first time, of dusty galaxies, at mm and sub-mm wavelengths, respectively.

Self-consistent treatment of cosmological structure formation and expansion within the context of classical general relativity may lead to "extra" expansion above that expected in a structureless universe. We argue that in comparison to an early-epoch, extrapolated Einstein-de Sitter model, about 10-15% "extra" expansion is sufficient at the present to render superfluous the "dark energy" 68% contribution to the energy density budget, and that this is observationally realistic.

The current cosmological model ($\Lambda$CDM) with the underlying FLRW metric relies on the assumption of local isotropy, hence homogeneity of the Universe. Difficulties arise when one attempts to justify this model as an average description of the Universe from first principles of general relativity, since in general, the Einstein tensor built from the averaged metric is not equal to the averaged stress--energy tensor. In this context, the discrepancy between these quantities is called "cosmological backreaction" and has been the subject of scientific debate among cosmologists and relativists for more than $20$ years. Here we present one of the methods to tackle this problem, i.e. averaging the scalar parts of the Einstein equations, together with its application, the cosmological mass function of galaxy clusters.

In relativistic inhomogeneous cosmology, structure formation couples to average cosmological expansion. A conservative approach to modelling this assumes an Einstein--de Sitter model (EdS) at early times and extrapolates this forward in cosmological time as a "background model" against which average properties of today's Universe can be measured. This requires adopting an early-epoch--normalised background Hubble constant $H_1^{bg}$. Here, we show that the $\Lambda$CDM model can be used as an observational proxy to estimate $H_1^{bg}$ rather than choose it arbitrarily. We assume (i) an EdS model at early times; (ii) a zero dark energy parameter; (iii) bi-domain scalar averaging---division of the spatial sections into over- and underdense regions; and (iv) virialisation (stable clustering) of collapsed regions. We find $H_1^{bg}= 37.7 \pm 0.4$ km/s/Mpc (random error only) based on a Planck $\Lambda$CDM observational proxy. Moreover, since the scalar-averaged expansion rate is expected to exceed the (extrapolated) background expansion rate, the expected age of the Universe should be much less than $2/(3 H_1^{bg}) = 17.3$ Gyr. The maximum stellar age of Galactic Bulge microlensed low-mass stars (most likely: 14.7 Gyr; 68\% confidence: 14.0--15.0 Gyr) suggests an age about a Gyr older than the (no-backreaction) $\Lambda$CDM estimate.

Systematic effects in dual-beam, differential, radio observations of extended objects are discussed in the context of the One Centimeter Receiver Array (OCRA). We use simulated samples of Sunyaev-Zel'dovich (SZ) galaxy clusters at low ($z<0.4$) and intermediate ($0.4<z<1.0$) redshifts to study the implications of operating at a single frequency (30 GHz) on the accuracy of extracting SZ flux densities and of reconstructing comptonization parameters with OCRA. We analyze dependences on cluster mass, redshift, observation strategy, and telescope pointing accuracy. Using $Planck$ data to make primary cosmic microwave background (CMB) templates, we test the feasibility of mitigating CMB confusion effects in observations of SZ profiles at angular scales larger than the separation of the receiver beams.

We summarise some of the main observational challenges for the standard Friedmann-Lemaitre-Robertson-Walker cosmological model and describe how results recently presented in the parallel session `Large--scale Structure and Statistics' (DE3) at the `Fourteenth Marcel Grossman Meeting on General Relativity' are related to these challenges.

The general relativistic description of cosmological structure formation is an important challenge from both the theoretical and the numerical point of views. In this paper we present a brief prescription for a general relativistic treatment of structure formation and a resulting mass function on galaxy cluster scales in a highly generic scenario. To obtain this we use an exact scalar averaging scheme together with the relativistic generalization of Zel'dovich's approximation (RZA) that serves as a closure condition for the averaged equations.

For about a decade, the baryon acoustic oscillation (BAO) peak at about 105 Mpc/h has provided a standard ruler test of the LCDM cosmological model, a member of the Friedmann--Lemaitre--Robertson--Walker (FLRW) family of cosmological models---according to which comoving space is rigid. However, general relativity does not require comoving space to be rigid. During the virialisation epoch, when the most massive structures form by gravitational collapse, it should be expected that comoving space evolves inhomogeneous curvature as structure grows. The BAO peak standard ruler should also follow this inhomogeneous evolution if the comoving rigidity assumption is false. This "standard" ruler has now been detected to be flexible, as expected under general relativity.

Backreaction in the cosmological context is a longstanding problem that is especially important in the present era of precise cosmology. The standard model of a homogeneous background plus density perturbations is most probably oversimplified and is expected to fail to fully account for the near-future observations of sub-percent precision. From a theoretical point of view, the problem of backreaction is very complicated and deserves careful examination. Recently, Green and Wald claimed in a series of papers to have developed a formalism to properly describe the influence of density inhomogeneities on average properties of the Universe, i.e., the backreaction effect. A brief discussion of this framework is presented, focussing on its drawbacks and on misconceptions that have arisen during the "backreaction debate".

No. In a number of papers Green and Wald argue that the standard FLRW model approximates our Universe extremely well on all scales, except close to strong field astrophysical objects. In particular, they argue that the effect of inhomogeneities on average properties of the Universe (backreaction) is irrelevant. We show that this latter claim is not valid. Specifically, we demonstrate, referring to their recent review paper, that (i) their two-dimensional example used to illustrate the fitting problem differs from the actual problem in important respects, and it assumes what is to be proven; (ii) the proof of the trace-free property of backreaction is unphysical and the theorem about it fails to be a mathematically general statement; (iii) the scheme that underlies the trace-free theorem does not involve averaging and therefore does not capture crucial non-local effects; (iv) their arguments are to a large extent coordinate-dependent, and (v) many of their criticisms of backreaction frameworks do not apply to the published definitions of these frameworks. It is therefore incorrect to infer that Green and Wald have proven a general result that addresses the essential physical questions of backreaction in cosmology.

In the standard model of cosmology, the Universe is static in comoving coordinates; expansion occurs homogeneously and is represented by a global scale factor. The baryon acoustic oscillation (BAO) peak location is a statistical tracer that represents, in the standard model, a fixed comoving-length standard ruler. Recent gravitational collapse should modify the metric, rendering the effective scale factor, and thus the BAO standard ruler, spatially inhomogeneous. Using the Sloan Digital Sky Survey, we show to high significance (P < 0.001) that the spatial compression of the BAO peak location increases as the spatial paths' overlap with superclusters increases. Detailed observational and theoretical calibration of this BAO peak location environment dependence will be needed when interpreting the next decade's cosmological surveys.

The Friedmann-Lemaitre-Robertson-Walker (FLRW) metric assumes comoving spatial rigidity of metrical properties. The curvature term in comoving coordinates is environment-independent and cannot evolve. In the standard model, structure formation is interpreted accordingly: structures average out on the chosen metrical background, which remains rigid in comoving coordinates despite nonlinear structure growth. The latter claim needs to be tested, since it is a hypothesis that is not derived using general relativity. We introduce a test of the comoving rigidity assumption by measuring the two-point auto-correlation function on comoving scales---assuming FLRW comoving spatial rigidity---in order to detect shifts in the baryon acoustic oscillation (BAO) peak location for Large Red Galaxy (LRG) pairs of the Sloan Digital Sky Survey Data Release 7. In tangential directions, subsets of pairs overlapping with superclusters or voids show the BAO peak. The tangential BAO peak location for overlap with Nadathur & Hotchkiss superclusters is $4.3\pm1.6$ Mpc/h less than that for LRG pairs unselected for supercluster overlap, and $6.6\pm2.8$ Mpc/h less than that of the complementary pairs. Liivamagi et al. superclusters give corresponding differences of $3.7\pm2.9$ Mpc/h and $6.3\pm2.6$ Mpc/h, respectively. We have found moderately significant evidence (Kolmogorov--Smirnov tests suggest very significant evidence) that the BAO peak location for supercluster-overlapping pairs is compressed by about 6% compared to that of the complementary sample, providing a potential challenge to FLRW models and a benchmark for predictions from models based on an averaging approach that leaves the spatial metric a priori unspecified.

We discuss the potential of a next generation space-borne CMB experiment for studies of extragalactic sources with reference to COrE+, a project submitted to ESA in response to the M4 call. We consider three possible options for the telescope size: 1m, 1.5m and 2m (although the last option is probably impractical, given the M4 boundary conditions). The proposed instrument will be far more sensitive than Planck and will have a diffraction-limited angular resolution. These properties imply that even the 1m telescope option will perform substantially better than Planck for studies of extragalactic sources. The source detection limits as a function of frequency have been estimated by means of realistic simulations. The most significant improvements over Planck results are presented for each option. COrE+ will provide much larger samples of truly local star-forming galaxies, making possible analyses of the properties of galaxies (luminosity functions, dust mass functions, star formation rate functions, dust temperature distributions, etc.) across the Hubble sequence. Even more interestingly, COrE+ will detect, at |b|> 30 deg, thousands of strongly gravitationally lensed galaxies. Such large samples are of extraordinary astrophysical and cosmological value in many fields. Moreover, COrE+ high frequency maps will be optimally suited to pick up proto-clusters of dusty galaxies, i.e. to investigate the evolution of large scale structure at larger redshifts than can be reached by other means. Thanks to its high sensitivity COrE+ will also yield a spectacular advance in the blind detection of extragalactic sources in polarization. This will open a new window for studies of radio source polarization and of the global properties of magnetic fields in star forming galaxies and of their relationships with SFRs.

The concordance model is successful in explaining numerous observable phenomena at the price of introducing an exotic source of unknown origin: dark energy. Dark energy dominance occurs at recent epochs, when we expect most cosmological structures to have already formed, and thus, when the error induced by forcing the homogeneous FLRW metric onto the data is expected to be the most significant. We propose a way to quantify the impact of deviations from homogeneity on the evolution of cosmological parameters. Using a multi-scale partitioning approach and the virialization fraction estimated from numerical simulations in an Einstein-de Sitter model, we obtain an observationally realistic distance modulus over redshifts 0 < z < 3 by a relativistic correction of the FLRW metric.

The 3-torus (T^3) FLRW model better fits the nearly zero large-scale auto-correlation of the WMAP CMB sky maps than the infinite flat model. The T^3 model's WMAP parameters imply approximately equal-redshift topological lensing at z \sim 6. We investigate observational strategies for rejecting the T^3 solution or providing candidate topologically lensed galaxy pairs. T^3 holonomies are applied to (i) existing z \sim 6 observations and (ii) simulated observations, creating multiply connected catalogues. Corresponding simply connected catalogues are generated. Each catalogue is analysed using a successive filter method and collecting matched quadruples. Quadruple statistics between the multiply and simply connected catalogues are compared. The expected rejection of the hypothesis, or detection of candidate topologically lensed galaxies, is possible at a significance of 5% for a pair of T^3 axis-centred northern and southern surveys if photometric redshift accuracy is \sigma(\zphot) < 0.01 for a pair of nearly complete 100 deg^2 surveys with a total of > 500 galaxies over 4.3 < z < 6.6, or for a pair of 196 deg^2 surveys with > 400 galaxies and \sigma(\zphot) < 0.02 over 4<z<7. Dropping the maximum time interval in a pair from \Delta t =1 Gyr/h to \Delta t =0.1 Gyr/h yields a requirement of \sigma(\zphot) < 0.005 or \sigma(\zphot) < 0.01, respectively. Millions of z \sim 6 galaxies will be observed over fields of these sizes during the coming decades, implying much stronger constraints. The question is not if the hypothesis will be rejected or confirmed, it is when.

The geometry of the dark energy and cold dark matter dominated cosmological model (LambdaCDM) is commonly assumed to be given by a Friedmann-Lemaitre-Robertson-Walker (FLRW) metric, i.e. it assumes homogeneity in the comoving spatial section. The homogeneity assumption fails most strongly at (i) small distance scales and (ii) recent epochs, implying that the FLRW approximation is most likely to fail at these scales. We use the virialisation fraction to quantify (i) and (ii), which approximately coincide with each other on the observational past light cone. For increasing time, the virialisation fraction increases above 10% at about the same redshift (sim 1-3) at which Omega_Lambda grows above 10% (approx 1.8). Thus, instead of non-zero Omega_Lambda, we propose an approximate, general-relativistic correction to the matter-dominated (Omega_m =1, Omega_Lambda=0), homogeneous metric (Einstein de Sitter, EdS). A low-redshift effective matter-density parameter of Omega_m^eff(0) = 0.26 pm 0.05 is inferred. Over redshifts 0 < z < 3, the distance modulus of the virialisation-corrected EdS model approximately matches the LambdaCDM distance modulus. This rough approximation assumes "old physics" (general relativity), not "new physics". Thus, pending more detailed calculations, we strengthen the claim that "dark energy" should be considered as an artefact of emerging average curvature in the void-dominated Universe, via a novel approach that quantifies the relation between virialisation and average curvature evolution.

The standard model of cosmology is dominated - at the present epoch - by dark energy. Its voids are rigid and Newtonian within a relativistic background. The model prevents them from becoming hyperbolic. Observations of rapid velocity flows out of voids are normally interpreted within the standard model that is rigid in comoving coordinates, instead of allowing the voids' density parameter to drop below critical and their curvature to become negative. Isn't it time to advance beyond nineteenth century physics and relegate dark energy back to the "no significant evidence" box?

The approximate homogeneity of spatial sections of the Universe is well supported observationally, but the inhomogeneity of the spatial sections is even better supported. Here, we consider the implications of inhomogeneity in dust models for the connectedness of spatial sections at early times. We consider a non-global Lemaitre-Tolman-Bondi (LTB) model designed to match observations, a more general, heuristic model motivated by the former, and two specific, global LTB models. We propose that the generic class of solutions of the Einstein equations projected back in time from the spatial section at the present epoch includes subclasses in which the spatial section evolves (with increasing time) smoothly (i) from being disconnected to being connected, or (ii) from being simply connected to being multiply connected, where the coordinate system is comoving and synchronous. We show that (i) and (ii) each contain at least one exact solution. These subclasses exist because the Einstein equations allow non-simultaneous big bang times. The two types of topology evolution occur over time slices that include significantly post-quantum epochs if the bang time varies by much more than a Planck time. In this sense, it is possible for cosmic topology evolution to be "mostly" classical.

V383Sco was discovered to be an eclipsing binary at the beginning of the XX century. This system has one of the longest orbital periods known (13.5yr) and was initially classified as a zet_Aur-type variable. It was then forgotten for decades. This study provides a detailed look at the V383Sco, using new data obtained around the last eclipse in 2007/8. There was a suspicion that this system could be similar to eclipsing systems with extensive dusty disks like EECep and eps_Aur. This and other, alternative hypotheses are considered. The ASAS-3 VI light curves have been used to examine photometric changes. Low-(LRS) and high-res.(HRS) spectra have been used for spectral classification, to analyse line profiles, as well as to determine the reddening, radial velocities (RVs) and distance. The SED was analysed. Using original numerical code, we performed a simplified model of the eclipse, taking into account the pulsations of one of the components. The LRS shows traces of molecular bands, characteristic of an M-type supergiant. The presence of this star in the system is confirmed by SED, by a strong dependence of the eclipse depth on the photometric bands, and by pulsational changes. The presence of a low excitation nebula around the system has been inferred from [OI] 6300A emission. Analysis of the RVs, reddening, and P-L relation for Mira-type stars imply a distance to the V383Sco of 8.4+-0.6 kpc. The distance to the nearby V381Sco is 6.4+-0.8 kpc. The very different and oppositely directed RVs of these systems (89.8 vs -178.8 km/s) seem to be in agreement with a bulge/bar kinematic model of the Galactic centre and inconsistent with purely circular motion. We have found evidence for the presence of a pulsating M-type supergiant in the V383Sco which periodically obscures the much more luminous F0I-type star, causing the deep (possibly total) eclipses which vary in duration and shape.

Heuristic approaches in cosmology bypass more difficult calculations that would more strictly agree with the standard Einstein equation. These give us the well-known Friedmann-Lemaitre-Robertson-Walker (FLRW) models, and, more recently, the feedback effect of the global topology of spatial sections on the acceleration of test particles. Forcing the FLRW heuristic model on observations leads to dark energy, which, pending fully relativistic calculations, is best interpreted as an artefact. Could topological acceleration also be an artefact of using a heuristic approach? A multiply connected exact solution of the Einstein equation shows that topological acceleration is present in at least one fully relativistic case---it is not an artefact of Newtonian-like thinking.

Two observational campaigns were carried out during the eclipses of EE Cep in 2003 and 2008/9 to verify whether the eclipsing body in the system is indeed a dark disk and to understand the observed changes in the depth and durations of the eclipses. Multicolour photometric data and spectroscopic observations at both low and high resolution were collected. We numerically modelled the variations in brightness and colour during the eclipses. We tested models with different disk structure. We considered the possibility of disk precession. The complete set of observational data collected during the last three eclipses are made available to the astronomical community. Two blue maxima in the colour indices were detected during these two eclipses, one before and one after the photometric minimum. The first (stronger) blue maximum is simultaneous with a "bump" that is very clear in all the UBVRI light curves. Variations in the spectral line profiles seem to be recurrent during each cycle. NaI lines always show at least three absorption components during the eclipse minimum and strong absorption is superimposed on the H_alpha emission. These observations confirm that the eclipsing object in EE Cep system is indeed a dark, dusty disk around a low luminosity object. The primary appears to be a rapidly rotating Be star that is strongly darkened at the equator and brightened at the poles. Some of the conclusions of this work require verification in future studies: (i) a complex, possibly multi-ring structure of the disk in EE Cep; (ii) our explanation of the "bump" observed during the last two eclipses in terms of the different times of obscuration of the hot polar regions of the Be star by the disk; and (iii) our suggested period of the disk precession (~11-12 P_orb) and predicted depth of about 2 mag the forthcoming eclipse in 2014.

It has previously been shown heuristically that the topology of the Universe affects gravity, in the sense that a test particle near a massive object in a multiply connected universe is subject to a topologically induced acceleration that opposes the local attraction to the massive object. This effect distinguishes different comoving 3-manifolds, potentially providing a theoretical justification for the Poincar\'e dodecahedral space observational hypothesis and a dynamical test for cosmic topology. It is necessary to check if this effect occurs in a fully relativistic solution of the Einstein equations that has a multiply connected spatial section. A Schwarzschild-like exact solution that is multiply connected in one spatial direction is checked for analytical and numerical consistency with the heuristic result. The T$^1$ (slab space) heuristic result is found to be relativistically correct. For a fundamental domain size of $L$, a slow-moving, negligible-mass test particle lying at distance $x$ along the axis from the object of mass $M$ to its nearest multiple image, where $GM/c^2 \ll x \ll L/2$, has a residual acceleration away from the massive object of $4\zeta(3) G(M/L^3)\,x$, where $\zeta(3)$ is Ap\'ery's constant. For $M \sim 10^14 M_\odot$ and $L \sim 10$ to $20\hGpc$, this linear expression is accurate to $\pm10%$ over $3\hMpc \ltapprox x \ltapprox 2\hGpc$. Thus, at least in a simple example of a multiply connected universe, the topological acceleration effect is not an artefact of Newtonian-like reasoning, and its linear derivation is accurate over about three orders of magnitude in $x$.

One of the FLRW models that best fits the WMAP sky maps of the CMB is the Poincare dodecahedral space. The optimal fit of this model to WMAP data was recently found using an optimal cross-correlation method, but the systematic error in the estimate of the matched-circle angular radius \alpha, or equivalently, the (comoving) size of the Universe 2\rinj (twice the injectivity radius), might be much higher than the random error. In order to increase the falsifiability of the model, it would be useful to reduce the uncertainty in this estimate and to estimate the fraction of the sky where multiply imaged gravitationally bound objects should potentially be detectable. "Matched discs" are defined in order to describe a useful subset of multiply imaged objects. The cross-correlation method at \ltapprox 1 \hGpc is applied to WMAP 7-year data in order to improve the estimate of \alpha. The improved matched-circle radius estimate is \alpha = 23 \pm 1.4 deg, where the uncertainty represents systematic error dependent on the choices of galactic mask and all-sky map. This is equivalent to 2\rinj = 18.2\pm 0.5 \hGpc for matter density parameter Omega_m=0.28\pm 0.02. The lowest redshift of multiply imaged objects is z=106\pm18. Multiply imaged high overdensity (rare) peaks visible during 200>z>106 should be present in matched discs of radius 14.8\pm2.3 deg. The accuracy in the matched circle radius estimate is considerably improved by using the higher resolution signal. The predicted matched discs (over 200>z>106) project to about 20% of the full sky. Since any object located exactly in the discs would be multiply imaged at equal redshifts, evolutionary effects would be small for objects that are nearly located in the discs.

We present 30 GHz Sunyaev Zel'dovich observations of a statistically complete sample of galaxy clusters with OCRA-p. The clusters are the 18 most X-ray luminous clusters at z > 0.2 in the ROSAT Brightest Cluster Sample. We correct for contaminant radio sources via supplementary observations with the Green Bank Telescope, also at 30 GHz, and remove a cluster that is contaminated by an unresolved X-ray source. All 17 remaining clusters have central SZ effects with Comptonisation parameter y_0 exceeding 1.9x10^-4, and 13 are detected at significance > 3 sigma. We use our data to examine scalings between y_0 and X-ray temperature, X-ray luminosity, and the X-ray mass proxy Y_X, and find good agreement with predictions from self-similar models of cluster formation, with an intrinsic scatter in y_0 of about 25%. We also comment on the success of the observations in the face of the contaminant source population, and the implications for upcoming cm-wave surveys.

Knowledge of the population of radio sources in the range ~2-200 GHz is important for understanding their effects on measurements of the Cosmic Microwave Background power spectrum. We report measurements of the 30 GHz flux densities of 605 radio sources from the Combined Radio All-sky Targeted Eight-GHz Survey (CRATES), which have been made with the One Centimetre Receiver Array prototype (OCRA-p) on the Torun 32-m telescope. The flux densities of sources that were also observed by WMAP and previous OCRA surveys are in broad agreement with those reported here, however a number of sources display intrinsic variability. We find a good correlation between the 30 GHz and Fermi gamma-ray flux densities for common sources. We examine the radio spectra of all observed sources and report a number of Gigahertz-peaked and inverted spectrum sources. These measurements will be useful for comparison to those from the Low Frequency Instrument of the Planck satellite, which will make some of its most sensitive observations in the region covered here.

In the time-ordered data (TOD) files of the WMAP CMB observations, there is an undocumented timing offset of -25.6 ms between the spacecraft attitude and radio flux density timestamps. If the offset induced an error during calibration of the raw TOD, then estimates of the WMAP CMB quadrupole might be substantially in error. A timing error during calibration would not only induce an artificial quadrupole-like signal in the mean sky map, it would also add variance per pixel. This variance would be present in the calibrated TOD. Low-resolution map-making as a function of timing offset should show a minimum variance for the correct timing offset. Three years of the calibrated, filtered WMAP 3-year TOD are compiled into sky maps at HEALPix resolution N_side=8, individually for each of the K, Ka, Q, V and W band differencing assemblies (DA's), as a function of timing offset. The median per map of the temperature fluctuation variance per pixel is calculated and minimised against timing offset. Minima are clearly present. The timing offsets that minimise the median variance are -38 \pm 8 ms (K, Ka), -27 \pm 3 ms (Q), -43 \pm 8 ms (V), and -47 \pm 194 ms (W), i.e. an average of -30 \pm 3 ms, where the WMAP collaboration's preferred offset is 0 \pm 1.7 ms. A non-parametric bootstrap analysis rejects the latter at a significance of 99.999%. The hypothesis of a -25.6 ms offset, suggested by Liu, Xiong & Li from the TOD file timing offset, is consistent with these minima. It is difficult to avoid the conclusion that the WMAP calibrated TOD and inferred maps are wrongly calibrated. CMB quadrupole estimates (3/pi)C_2 based on the incorrectly calibrated TOD are overestimated by roughly 64 \pm 6% (KQ85 mask) to 94 \pm 10% (KQ75 mask). Ideally, the WMAP map-making pipelines should be redone starting from the uncalibrated TOD and using the -25.6 ms timing offset correction.

Compact, flat Friedmann-Lemaitre-Robertson-Walker (FLRW) models have recently regained interest as a good fit to the observed cosmic microwave background temperature fluctuations. However, it is generally thought that a globally, exactly-flat FLRW model is theoretically improbable. Here, in order to obtain a probability space on the set F of compact, comoving, 3-spatial sections of FLRW models, a physically motivated hypothesis is proposed, using the density parameter Omega as a derived rather than fundamental parameter. We assume that the processes that select the 3-manifold also select a global mass-energy and a Hubble parameter. The inferred range in Omega consists of a single real value for any 3-manifold. Thus, the obvious measure over F is the discrete measure. Hence, if the global mass-energy and Hubble parameter are a function of 3-manifold choice among compact FLRW models, then probability spaces parametrised by Omega do not, in general, give a zero probability of a flat model. Alternatively, parametrisation by the injectivity radius r_inj ("size") suggests the Lebesgue measure. In this case, the probability space over the injectivity radius implies that flat models occur almost surely (a.s.), in the sense of probability theory, and non-flat models a.s. do not occur.

The largest uncollapsed inhomogeneity in the observable Universe is statistically represented in the quadrupole signal of the cosmic microwave background (CMB) sky maps as observed by the Wilkinson Microwave Anisotropy Probe (WMAP). The constant temporal offset of -25.6 ms between the timestamps of the spacecraft attitude and observational data records in the time-ordered data (TOD) of the WMAP observations was suspected to imply that previously derived all-sky CMB maps are erroneous, and that the quadrupole is in large part an artefact. The optimal focussing of bright objects in the Galactic Plane plays a key role in showing that no error occurred at the step of mapmaking from the calibrated TOD. Instead, the error had an effect when the uncalibrated TOD were calibrated. Estimates of the high-latitude quadrupole based on the wrongly calibrated WMAP maps are overestimated by about 15-60%.

A large fraction of the previously estimated WMAP CMB quadrupole signal would be an artefact of incorrect Doppler dipole subtraction if the hypothesis of a small timing interpolation error were correct. Observations of bright foreground objects constitute part of the time-ordered-data (TOD). Scans of an object in different directions should be shifted by the would-be timing error, causing a blurring effect. Three half-years of the calibrated, filtered WMAP TOD are compiled individually for the four W band differencing assemblies (DA's), with no masking of bright objects, giving 12 maps for each timing offset. Percentiles of the temperature-fluctuation distribution in each map at HEALPix resolution N_side=2048 are used to determine the dependence of all-sky image sharpness on the timing offset. In the W band, the 99.999% percentile, i.e. the temperature fluctuation in the approx 503-rd brightest pixel, is the least noisy percentile. Using this statistic, the hypothesis that a -25.6 ms offset relative to the timing adopted by the WMAP collaboration gives a focus at least as sharp as the uncorrected timing is rejected at 4.6\sigma significance. The Q and V band maps also reject the -25.6 ms offset hypothesis at high statistical significance. The requirement that the correct choice of timing offset must maximise image sharpness implies that the hypothesis of a timing error in the WMAP collaboration's compilation of the WMAP calibrated, filtered TOD is rejected at high statistical significance in each of the Q, V and W wavebands. However, the hypothesis that a timing error was applied during calibration of the raw TOD, leading to a dipole-induced difference signal, is not excluded by this method.

The lack of structure greater than 10 Gpc/h in Wilkinson Microwave Anisotropy Probe (WMAP) observations of the cosmic microwave background (CMB) favours compact Friedmann-Lemaitre-Robertson-Walker (FLRW) models of the Universe. The present best candidates based on observations are the Poincare dodecahedral space S^3/I^* and the 3-torus T^3. The residual gravity effect favours the Poincare space, while a measure space argument where the density parameter is a derived parameter favours flat spaces almost surely.

Is the Doppler interpretation of galaxy redshifts in a Friedmann-Lemaitre-Robertson-Walker (FLRW) model valid in the context of the approach to comoving spatial sections pioneered by de Sitter, Friedmann, Lemaitre and Robertson, i.e. according to which the 3-manifold of comoving space is characterised by both its curvature and topology? Holonomy transformations for flat, spherical and hyperbolic FLRW spatial sections are proposed. By quotienting a simply-connected FLRW spatial section by an appropriate group of holonomy transformations, the Doppler interpretation in a non-expanding Minkowski space-time, obtained via four-velocity parallel transport along a photon path, is found to imply that an inertial observer is receding from herself at a speed greater than zero, implying contradictory world-lines. The contradiction in the multiply-connected case occurs for arbitrary redshifts in the flat and spherical cases, and for certain large redshifts in the hyperbolic case. The link between the Doppler interpretation of redshifts and cosmic topology can be understood physically as the link between parallel transport along a photon path and the fact that the comoving spatial geodesic corresponding to a photon's path can be a closed loop in an FLRW model of any curvature. Closed comoving spatial loops are fundamental to cosmic topology.

Observational evidence of dark energy that makes the Universe nearly flat at the present epoch is very strong. We study the link between spatial continuity and dark energy. We assume that comoving space is a compact 3-manifold of constant curvature, described by a homogeneous Friedman-Lemaitre-Robertson-Walker metric. We assume that spatial continuity cannot be violated, i.e. that the global topology of the comoving section of the Universe cannot change during post-quantum epochs. We find that if the Universe was flat and compact during early epochs, then the presently low values of the radiation and matter densities imply that dark energy was created as a spatial continuity effect. Moreover, if the Universe is compact, then Omega_tot=1 is dynamically stable, where Omega_tot is the total density parameter in units of the critical density. Dark energy was observationally detected as a geometrical phenomenon. It is difficult to imagine a simpler explanation for dark energy than spatial continuity, finiteness and homogeneity.

The perturbed Friedmann-Lemaitre-Robertson-Walker models allow many different possibilities for the 3-manifold of the comoving spatial section of the Universe. It used to be thought that global properties of the spatial section, including the topology of the space, have no feedback effect on dynamics. However, an elementary, weak-limit calculation shows that in the presence of a density perturbation, a gravitational feedback effect that is algebraically similar to dark energy does exist. Moreover, the effect differs between different 3-spaces. Among the well-proportioned spaces, the effect disappears down to third order in several cases, and down to fifth order for the Poincare dodecahedral space S^3/I^*. The Poincare space, that which also is preferred in many observational analyses, is better-balanced than the other spaces.

It has generally been thought that in perturbed Friedmann-Lemaitre-Robertson-Walker models of the Universe, global topology should not have any feedback effects on dynamics. However, a weak-field limit heuristical argument, assuming a finite particle horizon for the transmission of gravitational signals, shows that a residual acceleration effect can occur. The nature of this effect differs algebraically between different constant curvature 3-manifolds. This potentially provides a selection mechanism for the 3-manifold of comoving space.

Small angular scale (high l) studies of cosmic microwave background anisotropies require accurate knowledge of the statistical properties of extragalactic sources at cm-mm wavelengths. We have used a 30 GHz dual-beam receiver (OCRA-p) on the Torun 32-m telescope to measure the flux densities of 121 sources in VSA fields selected at 15 GHz with the Ryle Telescope. We have detected 57 sources above a limiting flux density of 5mJy, of which 31 sources have a flux density greater than 10mJy, which is our effective completeness limit. From these measurements we derive a surface density of sources above 10mJy at 30 GHz of 2.0+/-0.4 per square degree. This is consistent with the surface density obtained by Mason et al. (2009) who observed a large sample of sources selected at a much lower frequency (1.4 GHz). We have also investigated the dependence of the spectral index distribution on flux density by comparing our results with those for sources above 1 Jy selected from the WMAP 22 GHz catalogue. We conclude that the proportion of steep spectrum sources increases with decreasing flux density, qualitatively consistent with the predictions of deZotti et al. (2005). We find no evidence for an unexpected population of sources whose spectra rise towards high frequencies, which would affect our ability to interpret current high resolution CMB observations at 30 GHz and above.

In a flat space, the global topology of comoving space can induce a weak acceleration effect similar to dark energy. Does a similar effect occur in the case of the Poincare dodecahedral space S^3/I^*? Does the effect distinguish the Poincare space from other well-proportioned spaces? The residual acceleration effect in the Poincare space is studied here using a massive particle and a nearby test particle of negligible mass, in S^3 embedded in R^4. The weak limit gravitational attraction on a test particle at distance r is set \propto [r_C \sin(r/r_C)]^{-2}, where r_C = curvature radius, in order to satisfy Stokes' theorem. A finite particle horizon large enough to include the adjacent topological images of the massive particle is assumed. The regular, flat, 3-torus T^3 is re-examined, and two other well-proportioned spaces, S^3/T^* and S^3/O^*, are also studied. The residual gravity effect occurs in all four cases. In a perfectly regular 3-torus of side length L_a, and in S^3/T^* and S^3/O^*, the highest order term in the residual acceleration is the third order term in the Taylor expansion of r/L_a (3-torus), or r/r_C, respectively. However, the Poincare dodecahedral space is unique among the four spaces. The third order cancels, leaving the fifth order term \sim \pm 300 (r/r_C)^5 as the most significant. Not only are three of the four perfectly regular well-proportioned spaces better balanced than most other multiply connected spaces in terms of the residual acceleration effect by a factor of about a million (setting r/L_a = r/r_C \sim 10^{-3}), but the fourth of these spaces is about 10^4 times better balanced than the other three. This is the Poincare dodecahedral space. Is this unique dynamical property of the Poincare space a clue towards a theory of cosmic topology?

The missing fluctuations problem in cosmic microwave background observations is naturally explained by well-proportioned small universe models. Among the well-proportioned models, the Poincare dodecahedral space is empirically favoured. Does gravity favour this space? The residual gravity effect is the residual acceleration induced by weak limit gravity from multiple topological images of a massive object on a nearby negligible mass test object. At the present epoch, the residual gravity effect is about a million times weaker in three of the well-proportioned spaces than in ill-proportioned spaces. However, in the Poincare space, the effect is 10,000 times weaker still, i.e. the Poincare space is about 10^{10} times "better balanced" than ill-proportioned spaces. Both observations and weak limit dynamics select the Poincare space to be special.

We aim to improve the surface of last scattering (SLS) optimal cross-correlation method in order to refine estimates of the Poincar\'e dodecahedral space (PDS) cosmological parameters. We analytically derive the formulae required to exclude points on the sky that cannot be members of close SLS-SLS cross-pairs. These enable more efficient pair selection without sacrificing uniformity of the underlying selection process. In certain cases this decreases the calculation time and increases the number of pairs per separation bin. (i) We recalculate Monte Carlo Markov Chains (MCMC) on the five-year WMAP data; and (ii) we seek PDS solutions in a small number of Gaussian random fluctuation (GRF) simulations. For 5 < alpha/deg < 60, a calculation speed-up of 3-10 is obtained. (i) The best estimates of the PDS parameters for the five-year WMAP data are similar to those for the three-year data. (ii) Comparison of the optimal solutions found by the MCMC chains in the observational map to those found in the simulated maps yields a slightly stronger rejection of the simply connected model using $\alpha$ than using the twist angle $\phi$. The best estimate of $\alpha$ implies that_given a large scale auto-correlation as weak as that observed,_ the PDS-like cross-correlation signal in the WMAP data is expected with a probability of less than about 10%. The expected distribution of $\phi$ from the GRF simulations is approximately Gaussian around zero, it is not uniform on $[-\pi,\pi]$. We infer that for an infinite, flat, cosmic concordance model with Gaussian random fluctuations, the chance of finding_both_ (a) a large scale auto-correlation as weak as that observed,_and_ (b) a PDS-like signal similar to that observed is less than about 0.015% to 1.25%.

In a multiply connected space, the two twins of the special relativity twin paradox move with constant relative speed and meet a second time without acceleration. The twins' situations appear to be symmetrical despite the need for one to be younger due to time dilation. Here, the suggestion that the apparent symmetry is broken by homotopy classes of the twins' worldlines is reexamined using space-time diagrams. It is found that each twin finds her own spatial path to have zero winding index and that of the other twin to have unity winding index, i.e. the twins' worldlines' relative homotopy classes are symmetrical. Although the twins' apparent symmetry is in fact broken by the need for the non-favoured twin to non-simultaneously identify spatial domain boundaries, the non-favoured twin {\em cannot} detect her disfavoured state if she only measures the homotopy classes of the two twins' projected worldlines, contrary to what was previously suggested. We also note that for the non-favoured twin, the fundamental domain can be chosen by identifying time boundaries (with a spatial offset) instead of space boundaries (with a temporal offset).

Several studies have proposed that the shape of the Universe may be a Poincare dodecahedral space (PDS) rather than an infinite, simply connected, flat space. Both models assume a close to flat FLRW metric of about 30% matter density. We study two predictions of the PDS model. (i) For the correct model, the spatial two-point cross-correlation function, $\ximc$, of temperature fluctuations in the covering space, where the two points in any pair are on different copies of the surface of last scattering (SLS), should be of a similar order of magnitude to the auto-correlation function, $\xisc$, on a single copy of the SLS. (ii) The optimal orientation and identified circle radius for a "generalised" PDS model of arbitrary twist $\phi$, found by maximising $\ximc$ relative to $\xisc$ in the WMAP maps, should yield $\phi \in \{\pm 36\deg\}$. We optimise the cross-correlation at scales < 4.0 h^-1 Gpc using a Markov chain Monte Carlo (MCMC) method over orientation, circle size and $\phi$. Both predictions were satisfied: (i) an optimal "generalised" PDS solution was found, with a strong cross-correlation between points which would be distant and only weakly correlated according to the simply connected hypothesis, for two different foreground-reduced versions of the WMAP 3-year all-sky map, both with and without the kp2 Galaxy mask: the face centres are $(l,b)_{i=1,6}\approx (184d, 62d), (305d, 44d), (46d, 49d), (117d, 20d), (176d, -4d), (240d, 13d) to within ~2d, and their antipodes; (ii) this solution has twist \phi= (+39 \pm 2.5)d, in agreement with the PDS model. The chance of this occurring in the simply connected model, assuming a uniform distribution $\phi \in [0,2\pi]$, is about 6-9%.

It has been suggested by Roukema and coworkers (hereafter R04) that the topology of the Universe as probed by the ``matched circles'' method using the first year release of the WMAP CMB data, might be that of the Poincar\'e dodecahedral space (PDS) model. An excess in the correlation of the ``identified circles'' was reported by R04, for circles of angular radius of ~11 deg for a relative phase twist -36deg, hinting that this could be due to a Clifford translation, if the hypothesized model were true. R04 did not however specify the statistical significance of the correlation signal. We investigate the statistical significance of the signal using Monte Carlo CMB simulations in a simply connected Universe, and present an updated analysis using the three-year WMAP data. We find that our analyses of the first and three year WMAP data provide results that are consistent with the simply connected space at a confidence level as low as 68%.

Could cosmic topology imply dark energy? We use a weak field (Newtonian) approximation of gravity and consider the gravitational effect from distant, multiple copies of a large, collapsed (virialised) object today (i.e. a massive galaxy cluster), taking into account the finite propagation speed of gravity, in a flat, multiply connected universe, and assume that due to a prior epoch of fast expansion (e.g. inflation), the gravitational effect of the distant copies is felt locally, from beyond the naively calculated horizon. We find that for a universe with a $T^1xR^2$ spatial section, the residual Newtonian gravitational force (to first order) provides an anisotropic effect that repels test particles from the cluster in the compact direction, in a way algebraically similar to that of dark energy. For a typical test object at comoving distance $\chi$ from the nearest dense nodes of the cosmic web of density perturbations, the pressure-to-density ratio $w$ of the equation of state in an FLRW universe, is w \sim - (\chi/L)^3, where $L$ is the size of the fundamental domain, i.e. of the universe. Clearly, |w|<<1. For a T^3 spatial section of exactly equal fundamental lengths, the effect cancels to zero. For a T^3 spatial section of unequal fundamental lengths, the acceleration effect is anisotropic in the sense that it will *tend to equalise the three fundamental lengths*. Provided that at least a modest amount of inflation occurred in the early Universe, and given some other conditions, multiple connectedness does generate an effect similar to that of dark energy, but the amplitude of the effect at the present epoch is too small to explain the observed dark energy density and its anisotropy makes it an unrealistic candidate for the observed dark energy.

Several analyses of the cosmic microwave background map made by the satellite WMAP suggest that the global shape of a spatial section of the Universe is that of a Poincare dodecahedral space. A summary of some of these analyses and a description of independent tests which should be able to either increase confidence in the hypothesis or else refute it to extremely high significance will be presented.

Withdrawn due to a mistake in jets' Doppler factor modeling. An updated version will appear after revisiting this work.

Three-dimensional catalogues of objects at cosmological distances can potentially yield candidate topologically lensed pairs of sets of objects, which would be a sign of the global topology of the Universe. In the spherical case, a necessary condition, which does not exist for either null or negative curvature, can be used to falsify such hypotheses, without needing to loop through a list of individual spherical 3-manifolds. This condition is that the isometry between the two sets of objects must be a root of the identity isometry in the covering space S^3. This enables numerical falsification of topological lensing hypotheses without needing to assume any particular spherical 3-manifold. By embedding S^3 in R^4, this condition can be expressed as the requirement that M^n = I for an integer n, where M is the matrix representation of the hypothesised lensing isometry and I is the identity. Moreover, this test becomes even simpler with the requirement that the two rotation angles, theta, phi, corresponding to the given isometry, satisfy 2\pi / \theta, 2\pi / \phi \in Z. The calculation of this test involves finding the two eigenplanes of the matrix M. A GNU General Public Licence numerical package, called eigenplane, is made available at http://cosmo.torun.pl/GPLdownload/eigen/ for finding the rotation angles and eigenplanes of an arbitrary isometry M of S^3.

Cosmic crystallography is based on the principle that peaks in the pair separation histogram (PSH) of objects in a catalogue should be induced by the high number of topologically lensed pairs that are separated by Clifford translations, in excess to ``random'' pairs of objects. Here we present modifications of this method that successively improve the signal-to-noise ratio by removing a large part of the noise and then false signals induced by selection effects. Given the transient nature of the most readily available tracer objects, active galactic nuclei (AGNs), the former is possible because a natural filter for removing many of the noise pairs is available: when counting pairs of objects in order to create PSHs, only those with nearly identical redshifts need to be counted. This redshift filter (a maximum value of $\Delta z/z = 0.005$) was applied to a compilation of AGN catalogues. Further noise was removed by applying a second filter, a maximum angle $\Delta \theta =0.075$ rad, and a minimum number of pairs $\protect\npairs=3$ to find each ``bunch of pairs'' (BoP) where the {\em vectors} (in Euclidean comoving space) defined by pairs are required to be nearly equal, whereas in the PSH only the {\em lengths} must be nearly equal. These filters reveal significant signals, which, however, are due to selection effects. A third filter, a minimum length $\Lselec=150$ {\hMpc} between the (parallel) vectors in a BoP, is found to effectively remove these selection effect pairs. After application of these successive filters, no significant topological signal was found.

Gorski et al (1999b) have earlier presented the outline of a pixelisation-to-spherical-coordinate transformation scheme which simultaneously satisfies three properties which are especially useful for rapid analyses of maps on a sphere: (i) equal spacing of pixels along lines of constant latitude, (ii) equal pixel `areas' (solid angles) and (iii) hierarchical scaling with increasing numbers of pixels. Their outline is based on the division of the sphere into twelve regions covering equal solid angles, which are hierarchically subdivided in a way compatible with these three criteria. In this paper, a complete derivation of this scheme is presented, including, in particular, (1) the angle theta^* defining the limit between polar and equatorial regions, and (2) the transformations from the unit interval [0,1] \wedge [0,1] to spherical coordinates in a polar region.

Luminet et al. (2003) suggested that WMAP data are better matched by a Poincar\'e dodecahedral FLRW model of global geometry, rather than by an infinite flat model. The analysis by Cornish et al. (2003) for angular radii 25-90 degrees failed to support this. Here, a matched circles analysis specifically designed to detect dodecahedral patterns of matched circles is performed over angular radii in the range 1-40 degrees on the one-year WMAP ILC map, using a correlation statistic and an rms difference statistic. Extreme value distributions of these statistics are calculated for left-handed and right-handed 36 degree `screw motions' (Clifford translations) when matching circles and for a zero (unphysical) rotation. The most correlated circles appear for circle radii of 11\pm1 degrees, for the left-handed screw motion, but not for the right-handed one, nor for the zero rotation. The favoured six dodecahedral face centres in galactic coordinates are (l,b)= (252, +65), (51, +51), (144,+38), (207,+10), (271,+3), (332,+25) and their opposites. The six pairs of circles_independently_ each favour a circle angular radius of 11\pm1 degrees. Whether or not these six circle pairs centred on dodecahedral faces match via a 36 degree rotation only due to unexpected statistical properties of the WMAP ILC map, or whether they match due to global geometry, it is clear that the WMAP ILC map has some unusual statistical properties which mimic a potentially interesting cosmological signal. The software for reading the WMAP data and for carrying out this analysis are released under the GNU General Public License.

Structures on very large scales (> 100 Mpc) have negligible peculiar motions, and are thus roughly fixed in comoving space. We looked for significant peaks at very large separation in the two-point correlation function -- corrected for redshift selection effects -- of a well convered subsample of 2378 quasars of the recently released 10k sample of the 2dF quasar survey. Dividing our sample in three redshift intervals, we find a peak at 244/h Mpc, which is perfectly comoving for a restricted set of cosmological parameters, namely Omega_m = 0.25+/-0.15 and Omega_Lambda=0.65+/-0.35 (both at 95% confidence). Assuming a flat Universe, we constrain the quintessence parameter w_Q < -0.35 (95% confidence). We discuss the compatibility of our analysis with possible peaks in the power spectrum.

The Universe is a physical object. Physical objects have shapes and sizes. General relativity is insufficient to describe the global shape and size of the Universe: the Hilbert-Einstein equations only treat limiting quantities towards an arbitrary point. Empirical work on measuring the shape and size of the Universe (formally: the ``3-manifold of the spatial hypersurface at constant cosmological time'', and, e.g. the ``injectivity diameter'' respectively) has progressed significantly in the late 1980's and the 1990's, using observational catalogues of galaxy clusters, of quasars and of the microwave background, though the analyses are still hindered by simplifying (and often observationally unsupported) assumptions. A review of the different observational strategies and claimed constraints was presented at the meeting.

Where's the matter? The answer seems to be: ``distributed according to a power spectrum with at least one feature (local maximum) near 120-130h^-1Mpc''. Analyses of the Iovino, Clowes & Shaver quasar candidate catalogue at z~2 and the 2dF Quasar Survey 10K Release (2QZ-10K) support this claim, which has previously been made both for low redshift survey analyses and for high z~3) redshift surveys will be presented. This feature (i) offers a comoving standard ruler which can lift the matter density-cosmological constant (Omega_m-Omega_Lambda) degeneracy and (ii) might be due either to baryonic acoustic oscillations or to Planck epoch physics which survives through inflation. (i) The 95% confidence constraint from the 2QZ-10K is Omega_m= 0.25\pm0.15, Omega_Lambda=0.60{\pm0.35}. This constraint is independent of cosmic microwave background constraints and type Ia supernovae constraints. The only assumptions required are (a) that the Universe satisfies a perturbed Friedmann-Lema\^{\i}tre-Robertson-Walker model with a possibly non-zero cosmological constant, (b) that the density perturbations in this model on large scales (>> 10h^-1Mpc) remain small (``linear'') and approximately spatially fixed in comoving coordinates, (c) that the statistics (power spectrum or correlation function) of the perturbations are redshift independent, and (d) that quasar redshifts are cosmological.

Detections of local maxima in the density perturbation spectrum at characteristic comoving scales L~100-200h^{-1}Mpc have previously been claimed. Here, this cosmic standard ruler is sought in the ``10K'' release of the 2dF QSO Redshift Survey (2QZ-10K), by estimating the comoving, spatial correlation functions \xi(r) of the three-dimensional distribution of the N=2378 quasars in the most completely observed and ``covered'' sky regions, over the three redshift ranges 0.6 < z < 1.1, 1.1 < z < 1.6 and 1.6 < z < 2.2. Because of the selection method of the survey and sparsity of the data, the analysis was done conservatively to avoid non-cosmological artefacts. (i) Avoiding 'a priori' estimates of the length scales of features, local maxima in \xi(r) are found in all three redshift ranges. The requirement that a local maximum be present in all three redshift ranges at a fixed comoving length scale implies strong, purely geometric constraints on the local cosmological parameters. The length scale of the local maximum common to the three redshift ranges is 2L= (244\pm17)h^{-1}Mpc. (ii) For a standard FLRW model, the matter density \Omm and cosmological constant \Omega_\Lambda are constrained to \Omm= 0.25\pm0.10, \Omega_\Lambda=0.65\pm0.25 (68% confidence), \Omm= 0.25\pm0.15, \Omega_\Lambda=0.60\pm0.35 (95%), respectively, *from the 2QZ-10K alone*. Independently of the SNe Ia data, the zero cosmological constant model (\Omega_\Lambda=0) is rejected at the 99.7% confidence level. (iii) For an effective quintessence (w_Q) model and zero curvature, w_Q<-0.5 (68%), w_Q<-0.35 (95%) are found, again *from the 2QZ-10K alone*.

If the Universe satisfies a perturbed Friedmann-Lema\^{\i}tre model, then the bright matter (e.g. radio-loud active galactic nuclei, RLAGNs) may be topologically lensed by global geometry. The generation of candidate topological lensing pairs of RLAGNs, which may be double images of single objects seen at very different celestial positions, provides very easily falsifiable (given moderate telescope time) cosmological geometry hypotheses.

The GAIA satellite will observe the Galaxy and its closest satellites in great detail. This should allow (i) dating past events of dwarf galaxies merging or interacting with the Galaxy, and much improved orbital parameters of the nine dwarf spheroidals and the Magellanic Clouds and (ii) dating of active galactic nucleus (AGN) phases of the Galaxy several Gigayears in the past, by detecting coeval star formation (e.g. open star clusters) that occurred along kpc scale Galactic jets. Both predictions (i) and (ii) of past history will be highly valuable for searching for or excluding topologically lensed images of the Milky Way at high redshift.

As the largest gravitationally collapsed objects, and as objects with a relatively low space density, clusters of galaxies offer one of the best sets of standard candles for trying to measure basic cosmological parameters such as the injectivity diameter 2r_{inj} (the shortest distance between two topologically lensed images of any object, e.g. cluster) and the out-diameter 2r_{+} (the maximum `size' of the Universe). Present constraints indicate that either of these may be smaller or larger than the horizon diameter.

Bottom-up hierarchical formation of dark matter haloes is not as monotonic as implicitly assumed in the Press-Schechter formalism: matter can be ejected into tidal tails, shells or low density `atmospheres'. The implications that the possible truncation of star formation in these tidal `debris' may have for observational galaxy statistics are examined here using the ArFus N-body plus semi-analytical galaxy modelling software. Upper and lower bounds on stellar losses implied by a given set of N-body simulation output data can be investigated by choice of the merging/identity criterion of haloes between successive N-body simulation output times. A median merging/identity criterion is defined and used to deduce an upper estimate of possible star formation and stellar population losses. A largest successor merging/identity criterion is defined to deduce an estimate which minimises stellar losses. In the N-body simulations studied, the debris losses are short range in length and temporary; maximum loss is around 16%. The induced losses for star formation and luminosity functions are strongest (losses of 10%-30%) for low luminosity galaxies and at intermediate redshifts (1 < z < 3). This upper bound on likely losses is smaller than present observational uncertainties. Hence, Press-Schechter based galaxy formation models are approximately valid despite ignoring loss of debris, provided that dwarf galaxy statistics are not under study.

Several recent observations using standard rulers and standard candles now suggest, either individually or in combination, that the Universe is close to flat, i.e. that the curvature radius is about as large as the horizon radius (\sim 10h^{-1}Gpc) or larger. Here, a method of distinguishing an almost flat universe from a precisely flat universe using a single observational data set, without using any microwave background information, is presented. The method (i) assumes that a standard ruler should have no preferred orientation (radial versus tangential) to the observer, and (ii) requires that the (comoving) length of the standard ruler be known independently (e.g. from low redshift estimates). The claimed feature at fixed {\em comoving} length in the power spectrum of density perturbations, detected among quasars, Lyman break galaxies or other high redshift objects, would provide an adequate standard candle to prove that the Universe is curved, if indeed it is curved. For example, a combined intrinsic and measurement uncertainty of 1% in the length of the standard ruler $\llss$ applied at a redshift of $z=3$ would distinguish an hyperbolic $(\Omega_m=0.2,\Omega_\Lambda=0.7)$ or a spherical $(\Omega_m=0.4,\Omega_\Lambda=0.7)$ universe from a flat one to 1-P > 95% confidence.

The inner product provides a conceptually and algorithmically simple method for calculating the comoving distance between two cosmological objects given their redshifts, right ascension and declination, and arbitrary constant curvature. The key to this is that just as a distance between two points `on' the surface of the ordinary 2-sphere ${\cal S}^2$ is simply an arc-length (angle multiplied by radius) in ordinary Euclidean 3-space ${\cal E}^3$, the distance between two points `on' a 3-sphere ${\cal S}^3$ (a 3-hyperboloid ${\cal H}^3$) is simply an `arc-length' in Euclidean 4-space ${\cal E}^4$ (Minkowski 4-space ${\cal M}^4$), i.e. an `hyper-angle' multiplied by the curvature radius of the 3-sphere (3-hyperboloid).

It is often stated that homogeneity and isotropy of the Universe are assumptions of the almost Friedmann-Lema^itre (FL) model (the hot big bang model), inspired from the Copernican Principle. However, only local homogeneity and isotropy are required by the model: multiply connected almost FL models are locally homogeneous and isotropic, but they can be globally anisotropic and/or globally inhomogeneous. Toy models are used here to show how global anisotropy and/or global inhomogeneity of an almost FL model could be shown directly in observations. This approach may avoid having to make any assumptions regarding global anisotropy and inhomogeneity.

In the almost Friedmann-Lema^itre model of the Universe, the density parameter, Omega_matter, and the cosmological constant, Omega_Lambda, measure curvature. Several linearly degenerate relations between these two parameters have recently been measured. Here, large scale structure correlations at \~100-150 h^-1 Mpc are found in the comoving three-dimensional separations of redshift z \approx 2 quasars. These function as a comoving standard rod of length L_LSS \approx 130+-10 h^-1 Mpc. A local maximum in the correlation function at \approx L_LSS/2 also appears to be significant. By combining separate radial and tangential standard ruler analyses, the lifting of the Omega_matter - Omega_Lambda linear degeneracy within a single data set is demonstrated for the first time.

Many different and complementary strategies for translating the basic principle of multiple topological imaging into observational analysis are now available, both for three-dimensional and two-dimensional catalogues.

The Hilbert-Einstein equations are insufficient to describe the geometry of the Universe, as they only constrain a local geometrical property: curvature. A global knowledge of the geometry of space, if possible, would require measurement of the topology of the Universe. Since the subject was discussed in 1900 by Schwarzschild, observational attempts to measure global topology have been rare for most of this century, but have accelerated in the 1990's due to the rapidly increasing amount of observations of non-negligible fractions of the observational sphere. A brief review of basic concepts of cosmic topology and of the rapidly growing gamut of diverse and complementary observational strategies for measuring the topology of the Universe is provided here.

Quel est, ou pourrait \^etre, la topologie globale de la partie spatiale de l'Univers ? L'Univers entier (pr\'ecis\'ement, l'hypersurface spatiale de celui-ci) est-il observable ? Les math\'ematiciens, les physiciens et les cosmologistes observationnels ont des approches diff\'erentes pour aborder ces questions qui restent ouvertes. What is, or could be, the global topology of spatial sections of the Universe? Is the entire Universe (spatial hypersurface thereof) observable ? Mathematicians, physicists and observational cosmologists have different strategies to approaching these questions which are not yet fully answered.

It has been suggested that if the Universe satisfies a flat, multiply connected, perturbed Friedmann-Lema^itre model, then cosmic microwave background data from the COBE satellite implies that the minimum size of the injectivity diameter (shortest closed spatial geodesic) must be larger than about two fifths of the horizon diameter. To show that this claim is misleading, a simple $T^2 \times R$ universe model of injectivity diameter a quarter of this size, i.e. a tenth of the horizon diameter, is shown to be consistent with COBE four year observational maps of the cosmic microwave background. This is done using the identified circles principle.

Several observational analyses suggest that matter is spatially structured at $\approx 130h^{-1}Mpc$ at low redshifts. This peak in the power spectrum provides a standard ruler in comoving space which can be used to compare the local geometry at high and low redshifts, thereby constraining the curvature parameters. It is shown here that this power spectrum peak is present in the observed quasar distribution at $z\sim 2$: qualitatively, via wedge diagrams which clearly show a void-like structure, and quantitatively, via one-dimensional Fourier analysis of the quasars' tangential distribution. The sample studied here contains 812 quasars. The method produces strong constraints (68% confidence limits) on the density parameter $\Omega_0$ and weaker constraints on the cosmological constant $\lambda_0$, which can be expressed by the relation $ \Omega_0 = (0.24 \pm0.15) + (0.10\pm0.08) \lambda_0 $. Independently of $\lambda_0$ (in the range $\lambda_0 \in [0,1]$), the constraint is $0.1 < \Omega_0 < 0.45$. Combination of the present results with SN Type Ia results yields $\Omega_0 = (0.30\pm0.11) + (0.57\pm0.11) (\lambda_0-0.7),$ $0.55 < \lambda_0 < 0.95,$ (68% confidence limits). This strongly supports the possibility that the observable universe satisfies a nearly flat, perturbed Friedmann-Lema\^{\i}tre-Robertson-Walker model, independently of any cosmic microwave background observations.

An initially highly biased value of the spatial correlation function, which decreases with time up to a transition redshift z_t ~ 1-2, was predicted theoretically at least as early as 1993, and shown to be consistent with the HDF-N angular correlation function estimates in 1997. Observational analyses are presented here which show (i) an HDF-N estimate of the correlation function amplitude of galaxies selected at z ~ 2 (by the UV drop-{\em in} technique), which supports the estimate z_t ~ 2 and (ii) an HDF-N estimate of the correlation function amplitude of galaxies selected at z ~ 3.7 by using photometric redshifts, which suggests that the correlation function amplitude evolves as (1+z)^{2\pm3.5} during epochs earlier than z_t.

The significance to which the cosmic microwave background observations by the satellite COBE can be used to refute a specific observationally based hypothesis for the global topology (3-manifold) of the Universe is investigated, by a new method of applying the principle of matched circle pairs. Moreover, it is shown that this can be done without assuming Gaussian distributions for the density perturbation spectrum. The Universe is assumed to correspond to a flat Friedmann-Lema\^{\i}tre model with a zero value of the cosmological constant. The 3-manifold is hypothesised to be a 2-torus in two directions, with a third axis larger than the horizon diameter. The positions and lengths of the axes are determined by the relative positions of the galaxy clusters Coma, RX~J1347.5-1145 and CL 09104+4109, assumed to be multiple topological images of a single, physical cluster. If the following two assumptions are valid: (i) the error estimates in the COBE DMR data are accurate estimates of the total random plus systematic error; and (ii) the temperature fluctuations are dominated by the na\"{\i}ve Sachs-Wolfe effect; then the distribution of the temperature differences between multiply imaged pixels is significantly wider than the uncertainty in the differences, and the candidate is rejected at the 94% level. This result is valid for either the `subtracted' or `combined' Analysed Science Data Sets, for either $10\deg$ or $20\deg$ smoothing, and is strengthened if suspected contaminated regions from the galactic centre and the Ophiuchus and Orion complexes are removed.

If the assumption that physical space has a trivial topology is dropped, then the Universe may be described by a multiply connected Friedmann-Lema\^{\i}tre model on a sub-horizon scale. Specific candidates for the multiply connected space manifold have already been suggested. How precisely would a significant detection of multiple topological images of a single object, or a region on the cosmic microwave background, (due to photons arriving at the observer by multiple paths which have crossed the Universe in different directions), constrain the values of the curvature parameters $\Omega_0$ and $\lambda_0$? The way that the constraints on $\Omega_0$ and $\lambda_0$ depend on the redshifts of multiple topological images and on their radial and tangential separations is presented and calculated. The tangential separations give the tighter constraints: multiple topological images of known types of astrophysical objects at redshifts $z \ltapprox 3$ would imply values of $\Omega_0$ and $\lambda_0$ preciser than $\sim 1%$ and $\sim 10%$ respectively. Cosmic microwave background `spots' identified with lower redshift objects by the Planck or MAP satellites would provide similar precision. This method is purely geometrical: no dynamical assumptions (such as the virial theorem) are required and the constraints are independent of the Hubble constant, $H_0.$

``One could imagine that as a result of enormously extended astronomical experience, the entire universe consists of countless identical copies of our Milky Way, that the infinite space can be partitioned into cubes each containing an exactly identical copy of our Milky Way. Would we really cling on to the assumption of infinitely many identical repetitions of the same world? >... We would be much happier with the view that these repetitions are illusory, that in reality space has peculiar connection properties so that if we leave any one cube through a side, then we immediately reenter it through the opposite side.'' (Schwarzschild 1900, translation 1998) Developments in the theoretical and observational sides of cosmic topology were slow for most of the century, but are now progressing rapidly, at the scale of most interest which is 1-10 h^-1 Gpc rather than 10kpc. The historical, mathematical and observational sides of this subject are briefly reviewed in this course.

Is the Universe (a spatial section thereof) finite or infinite? Knowing the global geometry of a Friedmann-Lema\^{\i}tre (FL) universe requires knowing both its curvature and its topology. A flat or hyperbolic (``open'') FL universe is {\em not} necessarily infinite in volume. Multiply connected flat and hyperbolic models are, in general, as consistent with present observations on scales of 1-20{\hGpc} as are the corresponding simply connected flat and hyperbolic models. The methods of detecting multiply connected models (MCM's) are presently in their pioneering phase of development and the optimal observationally realistic strategy is probably yet to be calculated. Constraints against MCM's on ~1-4 h^{-1} Gpc scales have been claimed, but relate more to inconsistent assumptions on perturbation statistics rather than just to topology. Candidate 3-manifolds based on hypothesised multiply imaged objects are being offered for observational refutation. The theoretical and observational sides of this rapidly developing subject have yet to make any serious contact, but the prospects of a significant detection in the coming decade may well propel the two together.

The study of the kinematics of galaxies within clusters or groups has the limitation that only one of the three velocity components and only two of the three spatial components of a galaxy position in six-dimensional phase space can normally be measured. However, if multiple topological images of a cluster exist, then the radial positions and sky plane mean velocities of galaxies in the cluster may also be measurable from photometry of the two cluster images. The vector arithmetic and principles of the analysis are presented. These are demonstrated by assuming the suggested topological identification of the clusters RX J1347.5-1145 and CL 09104+4109 to be correct and deducing the sky-plane relative velocity component along the axis common to both images of this would-be single cluster. Three out of four of the inferred transverse velocities are consistent with those expected in a rich cluster. A control sample of random `common' sky-plane axes, independent of the topological hypothesis, implies that this is not surprising. This shows that while galaxy kinematics are deducible from knowledge of cosmological topology, it is not easy to use them to refute a specific candidate manifold.

The amplitude of the angular two-point galaxy correlation function w(\theta) for galaxies at z~2 is estimated for galaxies in the Hubble Deep Field by using a U < 27 complete sub-sample. (i) It is confirmed that the amplitude of the correlation can be corrected for the integral constraint without having to make assumptions about the shape of the correlation function and by avoiding the introduction of linear error terms. The estimate using this technique is w(5'') = 0.10 \pm 0.09. (ii) If the biases introduced in faint galaxy selection due to obscuration by large objects are not corrected for by masking areas around them, then the estimate would be w(5'') =0.16\pm 0.07. (iii) The effective (3-D) galaxy pair separation at 5'' and this redshift range is ~ 25-250 /h kpc, so the correction to the spatial correlation function \xi(r) due to exclusion of overlapping galaxy dark matter haloes should be considered. For clustering stable in proper units in an \Omega=1,\lambda=0 universe, our w(5\arcs) estimate (a) implies a present-day correlation length of r_0 ~ 2.6^{+1.1}_{-1.7}/h Mpc if halo overlapping is ignored, but (b) for a present-day correlation length of r_0=5.5/h Mpc implies that a typical halo exclusion radius is r_halo=70^{+420}_{-30}/h kpc. (iv) The decreasing correlation period (DCP) of a high initial bias in the spatial correlation function is not detected at this redshift. For an \Omega=1,\lambda=0 universe and (proper) stable clustering, possible detections of the DCP in other work would imply that \xi at redshifts greater than z_t = 1.7\pm0.9 would be [(1+z)/(1+z_t)]^{2.1\pm3.6} times higher than at z_t, which is consistent with our lack of a detection at z ~ 2.

A local property of standard cosmology - growth of primordial perturbations - implies at least six different formation mechanisms of dwarf galaxies. A global property - the topology of the Universe - may enable direct observational study of the aging of individual galaxies. Dwarfs (by mass) can form by star formation in low mass dark matter haloes, which generally collapse early in hierarchical galaxy formation scenarios. However, several other formation paths exist, including that of low redshift, recently collapsed, low mass dark haloes containing low metallicity stellar populations. The latter possibility is purely a property of gravitational collapse --- no special mechanisms to retard star formation are required. Dwarfs by ``nurture'', as opposed to ``nature'', are also possible.

The space-like hypersurface of the Universe at the present cosmological time is a three-dimensional manifold. A non-trivial global topology of this space-like hypersurface would imply that the apparently observable universe (the sphere of particle horizon radius) could contain several images of the single, physical Universe. Recent three-dimensional techniques for constraining and/or detecting this topology are reviewed. Initial applications of these techniques using X-ray bright clusters of galaxies and quasars imply (weak) candidates for a non-trivial topology.

Several recent large length scale methods of detecting the topology of the Universe are briefly reviewed. Three-dimensional methods include those of Lehoucq, Luminet & Lachi\`eze-Rey (1996), Roukema (1996) and Roukema & Edge (1997). The two-dimensional ``circles'' method of Cornish, Spergel & Starkman (1996b) is also mentioned.

A method of deriving and using merging history trees of dark matter galaxy haloes directly from pure gravity N-body simulations is presented. This combines the full non-linearity of N-body simulations with the flexibility of the semi-analytical approach. Merging history trees derived from power-law initial perturbation spectrum simulations (for indices $n=-2$ and $n=0$) by Warren et al. (1992) are shown. As an example of a galaxy formation model, these are combined with evolutionary stellar population synthesis, via simple scaling laws for star formation rates, showing that if most star formation occurs during merger-induced bursts, then a nearly flat faint-end slope of the galaxy luminosity function may be obtained in certain cases. Interesting properties of hierarchical halo formation are noted: (1) In a given model, merger rates may vary widely between individual haloes, and typically 20%-30% of a halo's mass may be due to infall of uncollapsed material. (2) Small mass haloes continue to form at recent times: as expected, the existence of young, low redshift, low metallicity galaxies (e.g., Izotov et al. 1997) is consistent with hierarchical galaxy formation models. (3) For $n=-2,$ the halo spatial correlation function can have a very high initial bias due to the high power on large scales.

The topology of the observable Universe is not yet known. The most significant observational sign of a non-trivial topology would be multiple images (``ghosts'') of a single object at (in general) different sky positions and redshifts. It is pointed out that the previous search by Gott (1980) for ghost images of the Coma cluster can be extended by using highly X-ray luminous clusters of galaxies. This is likely to be more efficient than with other astrophysical objects viewable on these scales since (1) X-ray clusters would be at least as easy to identify if viewed from other angles as any other objects and (2) the X-ray emitting thermally heated gas is likely to be simpler than for other objects. Possibilities that the highly luminous cluster RXJ~1347.5-1145 ($z=0\.45$) has a ``ghost image'' at lower redshift are analysed. It is noted that RXJ~1347.5-1145, the Coma cluster and the cluster CL~09104+4109 form nearly a right angle ($\approx 88\deg$) with arms of nearly identical length ($970h^{-1}$ and $960h^{-1}$ Mpc respectively) for $\Omega_0=1, \lambda_0=0$ curvature ($h\equiv H_0/100 kms^{-1} Mpc^{-1}$). This is a clue that the three clusters could be ghost images of one and the same cluster, for a hypertoroidal topology. However, several arguments are presented that this relation is not physical.

Several recent measurements have been made of the angular correlation function $w(\theta,m)$ of faint galaxies in deep surveys (e.g., Hubble Deep Field, HDF). Are the measured correlations indicative of gravitational growth of primordial perturbations or of the relationship between galaxies and (dark matter dominated) galaxy haloes? A first step in answering this question is to determine the typical spatial separations of galaxies whose spatial correlations, $\xi(r,z),$ contribute most to w. The median spatial separation of galaxy pairs contributing to a fraction $p$ of $w$ for a galaxy survey is denoted by r_{eff} and compared with the perpendicular distance, $r_\perp,$ at the median redshift, z_{median}, of the galaxies. Over a wide range in $\xi$ growth rates $\epsilon$ and in $z_{median}, \r_eff(50%) \ltsim r_\perp$ and $\r_eff(90%) \sim 4 r_\perp.$ Values of r_eff for typical redshift distributions indicate that many $w$ measurements correspond to spatial correlations at comoving length scales well below 1 h^-1 Mpc. For $\Omega_0=1$ and $\lambda_0=0$, the correlation signal at 4\arcsec predominant in the Villumsen et al. (1996) estimates of $w$ for faint HDF galaxies corresponds to r_eff(50%)(HDF) \approx 40 h^-1 kpc; other cosmologies and angles up to 10\arcsec can increase this to r_eff(50%)(HDF) $\ltsim 200 h^-1 kpc. The proper separations are (1+z) times smaller. These scales are small: the faint galaxy $w$ measurements are at scales where halo and/or galaxy existence, let alone interactions, may modify the spatial correlation function. These measurements could maybe be used to probe the radial extent of haloes at high redshift.

During the epoch of galaxy formation, the formation of the first galaxies in regions of high overdensity may lead to an initial amplitude of the spatial correlation function of galaxies, $\xi_0,$ much higher than that expected in linear perturbation theory. This initially high ``bias'' would consequently decrease to the near-unityvalues expected from local observations. Such a ``Decreasing Correlation Period'' (DCP), detected in N-body simulations under certain conditions by several authors, is parametrised in a simple way, ending at a ``transition'' redshift $z_t.$ Consequences on $w_0,$ the amplitude of the angular correlation function of faint galaxies (for a fixed apparent magnitude range), have been estimated. The sensitivity of $w_0$ to $z_t,$ and how this is affected by the redshift distribution of the galaxies in the sample, by low redshift spatial correlation function behaviour, by cosmology, and by the rate at which $\xi_0$ decreases during the DCP, are shown and discussed. This simple model is intended as a first exploration of the effects of the DCP for a likely range of parameter space. Comparison with the Hubble Deep Field (HDF) estimate of $w_0$ (Villumsen et al. 1996) indicates that the DCP is compatible with observation for values of $z_t$ in the theoretically expected range $z_t \gtapprox 1.$ However, tighter observational constraints on the HDF redshift distribution and on low redshift correlation function growth will be needed either to detect the DCP unambiguously or to constrain $z_t$ to redshifts greater than those of galaxies in existing surveys.

Hot big bang cosmology says nothing about the topology of the Universe. A topology-independent algorithm is presented which is complementary to that of Lehoucq et al. 1996 and which searches for evidence of multi-connectedness using catalogues of astrophysically observed objects. The basis of this algorithm is simply to search for a quintuplet of quasars (over a region of a few hundred comoving Mpc) which can be seen in two different parts of our past time cone, allowing for a translation, an arbitrary rotation and possibly a reflection. This algorithm is demonstrated by application to the distribution of quasars between redshifts of $z=1$ and $z\approx4,$ i.e., at a comoving distance from the observer $1700 h^{-1}\mbox{\rm Mpc} \ltapprox d \ltapprox 3300 h^{-1}\mbox{\rm Mpc}.$ Two pairs of isometric quintuplets separated by more than {$300$\hMpc} are found. This is consistent with the number expected from Monte Carlo simulations in a simply connected Universe if the detailed anisotropy of sky coverage by the individual quasar surveys is taken into account. The linear transformation in (flat) comoving space from one quintuplet to another requires translations of {$353$\hMpc} and {4922\hMpc} respectively, plus a reflection in the former case, and plus both a rotation and a reflection in the latter. Since reflections are required, if these two matches were due to multi-connectedness, then the Universe would be non-orientable.