We present a description of the Australian Dark Energy Survey (OzDES) and summarise the results from its six years of operations. Using the 2dF fibre positioner and AAOmega spectrograph on the 3.9-metre Anglo-Australian Telescope, OzDES has monitored 771 AGN, classified hundreds of supernovae, and obtained redshifts for thousands of galaxies that hosted a transient within the 10 deep fields of the Dark Energy Survey. We also present the second OzDES data release, containing the redshifts of almost 30,000 sources, some as faint as $r_{\mathrm AB}=24$ mag, and 375,000 individual spectra. These data, in combination with the time-series photometry from the Dark Energy Survey, will be used to measure the expansion history of the Universe out to $z\sim1.2$ and the masses of hundreds of black holes out to $z\sim4$. OzDES is a template for future surveys that combine simultaneous monitoring of targets with wide-field imaging cameras and wide-field multi-object spectrographs.

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.

The standard model of cosmology is based on two unknown dark components that are uncoupled from each other. In this paper we investigate whether there is evidence for an interaction between these components of cold dark matter (CDM) and dark energy (DE). In particular, we focus on a minimal extension and reconstruct the interaction history at low-redshifts non-parametrically using a variation of the commonly used principal component analysis. Although we focus on the interaction in the dark sector, any significant deviation from the standard model that changes the expansion history of the Universe, should leave imprints detectable by our analysis. Thus, detecting signatures of interaction could also be indicative of other non-standard phenomena even if they are not the results of the interaction. It is thus interesting to note that the results presented in this paper do not provide support for the interaction in the dark sector, although the uncertainty is still quite large. In so far as interaction is present but undetectable using current data, we show from a Fisher forecast that forthcoming LSST and DESI surveys will be able to constrain a DM-DE coupling at $20\%$ precision --- enough to falsify the non-interacting scenario, assuming the presence of a modest amount of interaction.

Observations of the large-scale structure (LSS) implicitly assume an ideal FLRW observer with the ambient structure having no influence on the observer. However, due to correlations in the LSS, cosmological observables are dependent on the position of an observer. We investigate this influence in full generality for a weakly non-Gaussian random field, for which we derive expressions for angular spectra of large-scale structure observables conditional on a property of the large-scale structure that is typical for the observer's location. As an application, we then apply to the formalism to angular spectra of the weak gravitational lensing effect and provide numerical estimates for the resulting change on the spectra using linear structure formation. For angular weak lensing spectra we find the effect to be of order of a few percent, for instance we estimate for an overdensity of $\delta =0.5$ and multipoles up to $\ell=100$ the change in the weak lensing spectra to be approximately 4 percent. We show that without accounting for correlation between the density at observer's location and the weak gravitational lensing spectra, the values of the parameters $\Omega_m$ and $\sigma_8$ are underestimated by a few percent. Thus, this effect will be important when analysing data from future surveys such as Euclid, which aim at the percent-level precision. The effect is difficult to capture in simulations, as estimates of the number of numerical simulations necessary to quantify the effect are high.

The redshift drift, the change of cosmological redshift with time, is a direct consequence of the expansion of the Universe. Thus the measurement of the cosmological redshift drift will offer a direct test of our models of cosmology. The magnitude of the effect is very small, i.e. the spectral shift is of order of $10^{-10} - 10^{-9}$ over the period of a decade, but the next generation facilities such as ELT and SKA will be able to directly detect the expansion of our Universe by the year 2040. In this paper we focus on detectebility of this effect, including strategies of overcoming the kinematic contamination of the cosmological signal. We also show that the redshift drift directly impacts the change of flux. Thus apart from the redshift drift, measurements of the flux drift will provide an additional tool of detecting the expansion of the universe, including its acceleration. We discuss the strategies of detecting the flux drift and show that by including the flux drift signal to the redshift drift signal we boost the chances of a direct detection of the expansion of the Universe. We show that if only the stability of flux is at the level of $\Delta F/F \sim 10^{-6}$ then the SKA1-mid Array should be able to detect these effects, before the ETL and the full SKA. Thus, by including the flux drift into the SKA1-mid Array's analysis pipeline, we could be able to provide by mid-2030s a direct evidence of the expansion of the universe including its accelerating phase.

This paper investigates the growth of the most massive cosmological objects. We utilize the Simsilun simulation, which is based on the approximation of the silent universe. In the limit of spatial homogeneity and isotropy the silent universes reduce to the standard FLRW models. We show that within the approximation of the silent universe the formation of the most massive cosmological objects differs from the standard background-dependent approaches. For objects with masses above $10^{15} M_\odot$ the effect of spatial curvature (overdense regions are characterized by positive spatial curvature) leads to measurable effects. The effect is analogous to the effect that the background cosmological model has on the formation of these objects (i.e. the higher matter density and spatial curvature the faster the growth of cosmic structures). We measure this by the means of the mass function and show that the mass function obtained from the Simsilun simulation has a higher amplitude at the high-mass end compared to standard mass function such as the Press-Schechter or the Tinker mass function. For comparison, we find that the expected mass of most massive objects using the Tinker mass function is $4.4^{+0.8}_{-0.6} \, \times 10^{15} M_\odot$, whereas for the Simsilun simulation is $6.3^{+1.0}_{-0.8} \, \times 10^{15} M_\odot$. We conclude that the nonlinear relativistic effects could affect the formation of the most massive cosmological objects, leading to a relativistic environment-dependence of the growth rate of the most massive clusters.

The measurements of the Hubble constant reveal a tension between high-redshift (CMB) and low-redshift (distance ladder) constraints. So far neither observational systematics nor new physics has been successfully implemented to explain this tension away. This paper present a new solution to the Hubble constant problem. The solution is based on the Simsilun simulation (relativistic simulation of the large scale structure of the Universe) with the ray-tracing algorithm implemented. The initial conditions for the Simsilun simulation were set up as perturbations around the $\Lambda$CDM model. However, unlike in the Standard Cosmological Model (i.e. $\Lambda$CDM model + perturbations), within the Simsilun simulation relativistic and nonlinear evolution of cosmic structures leads to the phenomenon of emerging spatial curvature, where the mean spatial curvature evolves from spatial flatness of the early universe towards slightly curved present-day universe. Consqeuently, the present-day expansion rate is slightly faster compared to the spatially flat $\Lambda$CDM model. The results of the ray-tracing analysis show that the universe which starts with initial conditions consistent with the Planck constraints should have the Hubble constant $H_0 = 72.5 \pm 2.1$ km s$^{-1}$ Mpc$^{-1}$. When the Simsilun simulation was re-run with no inhomogeneities imposed, the Hubble constant inferred within such a homogeneous simulation was $H_0 = 68.1 \pm 2.0$ km s$^{-1}$ Mpc$^{-1}$. Thus, the inclusion of nonlinear relativistic evolution that leads to the emergence of the spatial curvature can explain why the low-redshift measurements favour higher values compared to high-redshift constraints and alleviate the tension between the CMB and distance ladder measurements of the Hubble constant.

The gravitational entropy and no-hair conjectures seems to predict contradictory future states of our Universe. The growth of the gravitational entropy is associated with the growth of inhomogeneity, while the no-hair conjecture argues that a universe dominated by dark energy should asymptotically approach a homogeneous and isotropic de Sitter state. The aim of this paper is to study these two conjectures. The investigation is based on the Simsilun simulation, which simulates the universe using the approximation of the Silent Universe. The Silent Universe is a solution to the Einstein equations that assumes irrotational, non-viscous, and insulated dust, with vanishing magnetic part of the Weyl curvature. The initial conditions for the Simsilun simulation are sourced from the Millennium simulation, which results with a realistically appearing but relativistic at origin simulation of a universe. The Simsilun simulation is evolved from the early universe (t = 25 Myr) till far future (t = 1000 Gyr). The results of this investigation show that both conjectures are correct. On global scales, a universe with a positive cosmological constant and non-positive spatial curvature does indeed approach the de Sitter state. At the same time it keeps generating the gravitational entropy.

Taipan is a multi-object spectroscopic galaxy survey starting in 2017 that will cover 2pi steradians over the southern sky, and obtain optical spectra for about two million galaxies out to z<0.4. Taipan will use the newly-refurbished 1.2m UK Schmidt Telescope at Siding Spring Observatory with the new TAIPAN instrument, which includes an innovative 'Starbugs' positioning system capable of rapidly and simultaneously deploying up to 150 spectroscopic fibres (and up to 300 with a proposed upgrade) over the 6-deg diameter focal plane, and a purpose-built spectrograph operating from 370 to 870nm with resolving power R>2000. The main scientific goals of Taipan are: (i) to measure the distance scale of the Universe (primarily governed by the local expansion rate, H_0) to 1% precision, and the structure growth rate of structure to 5%; (ii) to make the most extensive map yet constructed of the mass distribution and motions in the local Universe, using peculiar velocities based on improved Fundamental Plane distances, which will enable sensitive tests of gravitational physics; and (iii) to deliver a legacy sample of low-redshift galaxies as a unique laboratory for studying galaxy evolution as a function of mass and environment. The final survey, which will be completed within 5 years, will consist of a complete magnitude-limited sample (i<17) of about 1.2x10^6 galaxies, supplemented by an extension to higher redshifts and fainter magnitudes (i<18.1) of a luminous red galaxy sample of about 0.8x10^6 galaxies. Observations and data processing will be carried out remotely and in a fully-automated way, using a purpose-built automated 'virtual observer' software and an automated data reduction pipeline. The Taipan survey is deliberately designed to maximise its legacy value, by complementing and enhancing current and planned surveys of the southern sky at wavelengths from the optical to the radio.

Relativistic numerical cosmology is most often based either on the exact solutions of the Einstein equations, or perturbation theory, or weak-field limit, or the BSSN formalism. The Silent Universe provides an alternative approach to investigate relativistic evolution of cosmological systems. The silent universe is based on the solution of the Einstein equations in 1+3 comoving coordinates with additional constraints imposed. These constraints include: the gravitational field is sourced by dust and cosmological constant only, both rotation and magnetic part of the Weyl tensor vanish, and the shear is diagnosable. This paper describes the code simsilun (free software distributed under the terms of the reposi General Public License), which implements the equations of the Silent Universe. The paper also discusses applications of the Silent Universe and it uses the Millennium simulation to set up the initial conditions for the code simsilun. The simulation obtained this way consists of 16,777,216 worldlines, which are evolved from z=80 to z=0. Initially, the mean evolution (averaged over the whole domain) follows the evolution of the background $\Lambda$CDM model. However, once the evolution of cosmic structures becomes nonlinear, the spatial curvature evolves from $\Omega_K =0$ to $\Omega_K \approx 0.1$ at the present day. The emergence of the spatial curvature is associated with $\Omega_M$ and $\Omega_\Lambda$ being smaller by approximately 0.05 compared to the $\Lambda$CDM.

This paper investigates the phenomenon of emergence of spatial curvature. This phenomenon is absent in the Standard Cosmological Model, which has a flat and fixed spatial curvature (small perturbations are considered in the Standard Cosmological Model but their global average vanishes, leading to spatial flatness at all times). This paper shows that with the nonlinear growth of cosmic structures the global average deviates from zero. The analysis is based on the {\em silent universes} (a wide class of inhomogeneous cosmological solutions of the Einstein equations). The initial conditions are set in the early universe as perturbations around the $\Lambda$CDM model with $\Omega_m = 0.31$, $\Omega_\Lambda = 0.69$, and $H_0 = 67.8$ km s$^{-1}$ Mpc$^{-1}$. As the growth of structures becomes nonlinear, the model deviates from the $\Lambda$CDM model, and at the present instant if averaged over a domain ${\cal D}$ with volume $V = (2150\,{\rm Mpc})^3$ (at these scales the cosmic variance is negligibly small) gives: $\Omega_m^{\cal D} = 0.22$, $\Omega_\Lambda^{\cal D} = 0.61$, $\Omega_{\cal R}^{\cal D} = 0.15$ (in the FLRW limit $\Omega_{\cal R}^{\cal D} \to \Omega_k$), and $\langle H \rangle_{\cal D} = 72.2$ km s$^{-1}$ Mpc$^{-1}$. Given the fact that low-redshift observations favor higher values of the Hubble constant and lower values of matter density, compared to the CMB constraints, the emergence of the spatial curvature in the low-redshift universe could be a possible solution to these discrepancies.

This paper discusses the phenomenon of backreaction within the Szekeres model. Cosmological backreaction describes how the mean global evolution of the Universe deviates from the Friedmannian evolution. The analysis is based on models of a single cosmological environment and the global ensemble of the Szekeres models (of the Swiss-Cheese-type and Styrofoam-type). The obtained results show that non-linear growth of cosmic structures is associated with the growth of the spatial curvature $\Omega_{\cal R}$ (in the FLRW limit $\Omega_{\cal R} \to \Omega_k$). If averaged over global scales the result depends on the assumed global model of the Universe. Within the Swiss-Cheese model, which does have a fixed background, the volume average follows the evolution of the background, and the global spatial curvature averages out to zero (the background model is the $\Lambda$CDM model, which is spatially flat). In the Styrofoam-type model, which does not have a fixed background, the mean evolution deviates from the spatially flat $\Lambda$CDM model, and the mean spatial curvature evolves to from $\Omega_{\cal R} =0 $ at the CMB to $\Omega_{\cal R} \sim 0.1$ at $z =0$. If the Styrofoam-type model correctly captures evolutionary features of the real Universe then one should expect that in our Universe, the spatial curvature should build up (local growth of cosmic structures) and its mean global average should deviate from zero (backreaction). As a result, this paper predicts that the low-redshift Universe should not be spatially flat (i.e. $\Omega_k \ne 0$, even if in the early Universe $\Omega_k = 0$) and therefore when analysing low-$z$ cosmological data one should keep $\Omega_k$ as a free parameter and independent from the CMB constraints.

Astronomical observations reveal hierarchical structures in the Universe, from galaxies, groups of galaxies, clusters and superclusters, to filaments and voids. On the largest scales it seems that some kind of statistical homogeneity can be observed. As a result, modern cosmological models are based on spatially homogeneous and isotropic solutions of the Einstein equations, and the evolution of the universe is approximated by the Friedmann equations. In parallel to standard homogeneous cosmology, the field of inhomogeneous cosmology and backreaction is being developed. This field investigates whether small scale inhomogeneities via non-linear effects can backreact and alter the properties of the Universe on its largest scales, leading to a non-Friedmannian evolution. This paper presents the current status of inhomogeneous cosmology and backreaction. It also discusses future prospects of the field of inhomogeneous cosmology, which is based on a survey of 50 academics working in the field of inhomogeneous cosmology.

The Universe on scales $10-100~h^{-1}$ Mpc is dominated by a cosmic web of voids, filaments, sheets and knots of galaxy clusters. These structures participate differently in the global expansion of the Universe: from non-expanding clusters to the above average expansion rate of voids. In this paper we characterize Hubble expansion anisotropies in the COMPOSITE sample of 4534 galaxies and clusters. We concentrate on the dipole and quadrupole in the rest frame of the Local Group. These both have statistically significant amplitudes. These anisotropies, and their redshift dependence, cannot be explained solely by a boost of the Local Group in the Friedmann-Lema\^itre-Robertson-Walker (FLRW) model which expands isotropically in the rest frame of the cosmic microwave background (CMB) radiation. We simulate the local expansion of the Universe with inhomogeneous Szekeres models, which match the standard FLRW model on $> 100~h^{-1}$ Mpc scales but exhibit nonkinematic relativistic differential expansion on small scales. We restrict models to be consistent with observed CMB temperature anisotropies, while simultaneously fitting the redshift variation of the Hubble expansion dipole. We include features to account for both the Local Void and the "Great Attractor". While this naturally accounts for the Hubble expansion and CMB dipoles, the simulated quadrupoles are smaller than observed. Further refinement to incorporate additional structures may improve this. This would enable a test of the hypothesis that some large angle CMB anomalies result from failing to treat the relativistic differential expansion of the background geometry, a natural feature of solutions to Einstein's equations not included in the current standard model of cosmology.

We present a new test of the validity of the Friedmann-Lemaitre-Robertson-Walker (FLRW) metric, based on comparing the distance from redshift 0 to $z_1$ and from $z_1$ to $z_2$ to the distance from $0$ to $z_2$. If the universe is described by the FLRW metric, the comparison provides a model-independent measurement of spatial curvature. The test relies on geometrical optics, it is independent of the matter content of the universe and the applicability of the Einstein equation on cosmological scales. We apply the test to observations, using the Union2.1 compilation of supernova distances and Sloan Lens ACS Survey galaxy strong lensing data. The FLRW metric is consistent with the data, and the spatial curvature parameter is constrained to be $-1.22<\Omega_{K0}<0.63$, or $-0.08<\Omega_{K0}<0.97$ with a prior from the cosmic microwave background and the local Hubble constant, though modelling of the lenses is a source of significant systematic uncertainty.

Do current observational data confirm the assumptions of the cosmological principle, or is there statistical evidence for deviations from spatial homogeneity on large scales? To address these questions, we developed a flexible framework based on spherically symmetric, but radially inhomogeneous Lemaitre-Tolman-Bondi (LTB) models with synchronous Big Bang. We expanded the (local) matter density profile in terms of flexible interpolation schemes and orthonormal polynomials. A Monte Carlo technique in combination with recent observational data was used to systematically vary the shape of these profiles. In the first part of this article, we reconsider giant LTB voids without dark energy to investigate whether extremely fine-tuned mass profiles can reconcile these models with current data. While the local Hubble rate and supernovae can easily be fitted without dark energy, however, model-independent constraints from the Planck 2013 data require an unrealistically low local Hubble rate, which is strongly inconsistent with the observed value; this result agrees well with previous studies. In the second part, we explain why it seems natural to extend our framework by a non-zero cosmological constant, which then allows us to perform general tests of the cosmological principle. Moreover, these extended models facilitate explorating whether fluctuations in the local matter density profile might potentially alleviate the tension between local and global measurements of the Hubble rate, as derived from Cepheid-calibrated type Ia supernovae and CMB experiments, respectively. We show that current data provide no evidence for deviations from spatial homogeneity on large scales. More accurate constraints are required to ultimately confirm the validity of the cosmological principle, however.

Doppler lensing is the apparent change in object size and magnitude due to peculiar velocities. Objects falling into an overdensity appear larger on its near side, and smaller on its far side, than typical objects at the same redshifts. This effect dominates over the usual gravitational lensing magnification at low redshift. Doppler lensing is a promising new probe of cosmology, and we explore in detail how to utilize the effect with forthcoming surveys. We present cosmological simulations of the Doppler and gravitational lensing effects based on the Millennium simulation. We show that Doppler lensing can be detected around stacked voids or unvirialised over-densities. New power spectra and correlation functions are proposed which are designed to be sensitive to Doppler lensing. We consider the impact of gravitational lensing and intrinsic size correlations on these quantities. We compute the correlation functions and forecast the errors for realistic forthcoming surveys, providing predictions for constraints on cosmological parameters. Finally, we demonstrate how we can make 3-D potential maps of large volumes of the Universe using Doppler lensing.

We study the nonlinear regime of Unified Dark Energy models, using Generalized Chaplygin Gas cosmologies as a representative example, and introduce a new parameter characterizing the level of small scale clustering in these scenarios. We show that viable Generalized Chaplygin Gas cosmologies, consistent with the most recent observational constraints, may be constructed for any value of the Generalized Chaplygin Gas parameter by considering models with a sufficiently high level of nonlinear clustering.

This paper studies intermediate homogenization of inhomogeneous cosmological models. It shows that spherically symmetric models, regardless of the equation of state, can undergo intermediate homogenization, i.e. a model can approach a homogeneous and isotropic state (which acts as a saddle point) from a relatively wide range of initial inhomogeneous conditions. The homogenization is not permanent - just temporary. Eventually the model evolves toward a future inhomogeneous state. We also looked at the problem of the gravitational entropy. All definitions of entropy that we checked give decreasing gravitational entropy during the homogenization process. Thus, we should either accept that gravitational entropy can decrease or try to define it in other ways than just via density gradients, as these decrease during homogenization.

More than half of the volume of our Universe is occupied by cosmic voids. The lensing magnification effect from those under-dense regions is generally thought to give a small dimming contribution: objects on the far side of a void are supposed to be observed as slightly smaller than if the void were not there, which together with conservation of surface brightness implies net reduction in photons received. This is predicted by the usual weak lensing integral of the density contrast along the line of sight. We show that this standard effect is swamped at low redshifts by a relativistic Doppler term that is typically neglected. Contrary to the usual expectation, objects on the far side of a void are brighter than they would be otherwise. Thus the local dynamics of matter in and near the void is crucial and is only captured by the full relativistic lensing convergence. There are also significant nonlinear corrections to the relativistic linear theory, which we show actually under-predicts the effect. We use exact solutions to estimate that these can be more than 20% for deep voids. This remains an important source of systematic errors for weak lensing density reconstruction in galaxy surveys and for supernovae observations, and may be the cause of the reported extra scatter of field supernovae located on the edge of voids compared to those in clusters.

The notion of an apparent horizon (AH) in a collapsing object can be carried over from the Lema\^{\i}tre -- Tolman (L--T) to the quasispherical Szekeres models in three ways: 1. Literally by the definition -- the AH is the boundary of the region, in which every bundle of null geodesics has negative expansion scalar. 2. As the locus, at which null lines that are as nearly radial as possible are turned toward decreasing areal radius $R$. These lines are in general nongeodesic. The name "absolute apparent horizon" (AAH) is proposed for this locus. 3. As the boundary of a region, where null \textit{geodesics} are turned toward decreasing $R$. The name "light collapse region" (LCR) is proposed for this region (which is 3-dimensional in every space of constant $t$); its boundary coincides with the AAH. The AH and AAH coincide in the L--T models. In the quasispherical Szekeres models, the AH is different from (but not disjoint with) the AAH. Properties of the AAH and LCR are investigated, and the relations between the AAH and the AH are illustrated with diagrams using an explicit example of a Szekeres metric. It turns out that an observer who is already within the AH is, for some time, not yet within the AAH. Nevertheless, no light signal can be sent through the AH from the inside. The analogue of the AAH for massive particles is also considered.

The Universe is inhomogeneous, and yet it seems to be incredibly well-characterised by a homogeneous relativistic model. One of the current challenges is to accurately characterise the properties of such a model. In this paper we explore how inhomogeneities may affect the overall optical properties of the Universe by quantifying how they can bias the redshift-distance relation in a number of toy models that mimic the real Universe. The models that we explore are statistically homogeneous on large scales. We find that the effect of inhomogeneities is of order of a few percent, which can be quite important in precise estimation of cosmological parameters. We discuss what lessons can be learned to help us tackle a more realistic inhomogeneous universe.

Recently, inhomogeneous generalisations of the Friedmann-Lemaitre-Robertson-Walker cosmological models have gained interest in the astrophysical community and are more often employed to study cosmological phenomena. However, in many papers the inhomogeneous cosmological models are treated as an alternative to the FLRW models. In fact, they are not an alternative, but an exact perturbation of the latter, and are gradually becoming a necessity in modern cosmology. The assumption of homogeneity is just a first approximation introduced to simplify equations. So far this assumption is commonly believed to have worked well, but future and more precise observations will not be properly analysed unless inhomogeneities are taken into account. This paper reviews recent developments in the field and shows the importance of an inhomogeneous framework in the analysis of cosmological observations.

This paper presents the algorithm for determining the Lemaitre-Tolman (LT) model that best fits given datasets for maximum stellar ages, and SNIa luminosities, both as functions of redshift. It then applies it to current cosmological data. Special attention must be given to the handling of the origin, and the region of the maximum diameter distances. As with a previous combination of datasets (galaxy number counts and luminosity distances versus redshift), there are relationships that must hold at the region of the maximum diameter distance, which are unlikely to be obeyed exactly by real data. We show how to make corrections that enable a self-consistent solution to be found. We address the questions of the best way to approximate discrete data with smooth functions, and how to estimate the uncertainties of the output - the 3 free functions that determine a specific LT metric. While current data does not permit any confidence in our results, we show that the method works well, and reasonable LT models do fit with or without a cosmological constant.

We use different particular classes of axially symmetric Szekeres Swiss-cheese models for the study of the apparent dimming of the supernovae of type Ia. We compare the results with those obtained in the corresponding Lemaitre-Tolman Swiss-cheese models. Although the quantitative picture is different the qualitative results are comparable, i.e, one cannot fully explain the dimming of the supernovae using small scale ~50 Mpc inhomogeneities. To fit successfully the data we need structures of order of 500 Mpc size or larger. However, this result might be an artifact due to the use of axial light rays in axially symmetric models. Anyhow, this work is a first step in trying to use Szekeres Swiss-cheese models in cosmology and it will be followed by the study of more physical models with still less symmetry.

Inhomogeneous cosmological models are able to fit cosmological observations without dark energy under the assumption that we live close to the "center" of a very large-scale under-dense region. Most studies fitting observations by means of inhomogeneities also assume spherical symmetry, and thus being at (or very near) the center may imply being located at a very special and unlikely observation point. We argue that such spherical voids should be treated only as a gross first approximation to configurations that follow from a suitable smoothing out of the non-spherical part of the inhomogeneities on angular scales. In this Letter we present a toy construction that supports the above statement. The construction uses parts of the Szekeres model, which is inhomogeneous and anisotropic thus it also addresses the limitations of spherical inhomogeneities. By using the thin-shell approximation (which means that the Israel-Darmois continuity conditions are not fulfilled between the shells) we construct a model of evolving cosmic structures, containing several elongated supercluster-like structures with underdense regions between them, which altogether provides a reasonable coarse-grained description of cosmic structures. While this configuration is not spherically symmetric, its proper volume average yields a spherical void profile of 250 Mpc that roughly agrees with observations. Also, by considering a non-spherical inhomogeneity, the definition of a "center" location becomes more nuanced, and thus the constraints placed by fitting observations on our position with respect to this location become less restrictive.

The standard analysis of the CMB data assumes that the distance to the last scattering surface can be calculated using the distance-redshift relation as in the Friedmann model. However, in the inhomogeneous universe, even if <\delta\rho> =0, the distance relation is not the same as in the unperturbed universe. This can be of serious consequences as a change of distance affects the mapping of CMB temperature fluctuations into the angular power spectrum. In addition, if the change of distance is relatively uniform no new temperature fluctuations are generated. It is therefore a different effect than the lensing or ISW effects which introduce additional CMB anisotropies. This paper shows that the accuracy of the CMB analysis can be impaired by the accuracy of calculation of the distance within the cosmological models. Since this effect has not been fully explored before, to test how the inhomogeneities affect the distance-redshift relation, several methods are examined: the Dyer-Roeder relation, lensing approximation, and non-linear Swiss-Cheese model. In all cases, the distance to the last scattering surface is different than when homogeneity is assumed. The difference can be as low as 1% and as high as 80%. Excluding extreme cases, the distance changes by about 20-30%. Since the distance to the last scattering surface is set by the position of the CMB peaks, in order to have a good fit, the distance needs to be adjusted. After correcting the distance, the cosmological parameters change. Therefore, a not properly estimated distance to the last scattering surface can be a major source of systematics. This paper shows that if inhomogeneities are taken into account when calculating the distance then models with positive spatial curvature and with \Omega_\Lambda ~ 0.8-0.9 are preferred. The \Lambda CDM model in most cases, is at odds with the current data.

The distance-redshift relation plays an important role in cosmology. In the standard approach to cosmology it is assumed that this relation is the same as in the homogeneous universe. As the real universe is not homogeneous there are several methods to calculate the correction. The weak lensing approximation and the Dyer-Roeder relation are one of them. This paper establishes a link between these two approximations. It is shown that if the universe is homogeneous with only small, vanishing after averaging, density fluctuations along the line of sight, then the distance correction is negligible. It is also shown that a vanishing 3D average of density fluctuations does not imply that the mean of density fluctuations along the line of sight is zero. In this case, even within the linear approximation, the distance correction is not negligible. The modified version of the Dyer-Roeder relation is presented and it is shown that this modified relation is consistent with the correction obtained within the weak lensing approximation. The correction to the distance for a source at z ~ 2 is of order of a few percent. Thus, with an increasing precision of cosmological observations an accurate estimation of the distance is essential. Otherwise errors due to miscalculation the distance can become a major source of systematics.

Dark energy is now one of the most important and topical problems in cosmology. The first step to reveal its nature is to detect the evolution of dark energy or to prove beyond doubt that the cosmological constant is indeed constant. However, in the standard approach to cosmology, the Universe is described by the homogeneous and isotropic Friedmann models. Here we show that in the perturbed universe (even if perturbations vanish if averaged over sufficiently large scales) the distance-redshift relation is not the same as in the unperturbed universe. This has a serious consequence when studying the nature of dark energy and, as shown here, can impair the analysis and studies of dark energy. The analysis is based on two methods: the linear lensing approximation and the non-linear Szekeres Swiss-Cheese model. The inhomogeneity scale is ~50 Mpc, and both models have the same density fluctuations along the line of sight. The comparison between linear and non-linear methods shows that non-linear corrections are not negligible. When inhomogeneities are present the distance changes by several percent. To show how this change influences the measurements of dark energy, ten future observations with 2% uncertainties are generated. It is shown the using the standard methods (i.e. under the assumption of homogeneity) the systematics due to inhomogeneities can distort our analysis, and may lead to a conclusion that dark energy evolves when in fact it is constant (or vice versa). Therefore, if future observations are analysed only within the homogeneous framework then the impact of inhomogeneities (such as voids and superclusters) can be mistaken for evolving dark energy. Since the robust distinction between the evolution and non-evolution of dark energy is the first step to understanding the nature of dark energy a proper handling of inhomogeneities is essential.

Lema\^itre - Tolman (L--T) toy models with a central observer have been used to study the effect of large scale inhomogeneities on the SN Ia dimming. Claims that a giant void is mandatory to explain away dark energy in this framework are currently dominating. Our aim is to show that L-T models exist that reproduce a few features of the $\Lambda$CDM model, but do not contain the giant cosmic void. We propose to use two sets of data - the angular diameter distance together with the redshift-space mass-density and the angular diameter distance together with the expansion rate - both defined on the past null cone as functions of the redshift. We assume that these functions are of the same form as in the $\Lambda$CDM model. Using the Mustapha-Hellaby-Ellis algorithm, we numerically transform these initial data into the usual two L-T arbitrary functions and solve the evolution equation to calculate the mass distribution in spacetime. For both models, we find that the current density profile does not exhibit a giant void, but rather a giant hump. However, this hump is not directly observable, since it is in a spacelike relation to a present observer. The alleged existence of the giant void was a consequence of the L-T models used earlier because their generality was limited a priori by needless simplifying assumptions, like, for example, the bang-time function being constant. Instead, one can feed any mass distribution or expansion rate history on the past light cone as initial data to the L-T evolution equation. When a fully general L-T metric is used, the giant void is not implied.

The present-day Universe appears to be homogeneous on very large scales. Yet when the casual structure of the early Universe is considered, it becomes apparent that the early Universe must have been highly inhomogeneous. The current paradigm attempts to answer this problem by postulating the inflation mechanism However, inflation in order to start requires a homogeneous patch of at least the horizon size. This paper examines if dynamical processes of the early Universe could lead to homogenization. In the past similar studies seem to imply that the set of initial conditions that leads to homogenization is of measure zero. This essay proves contrary: a set of initial conditions for spontaneous homogenization of cosmological models can form a set of non-zero measure.

A number of misunderstandings about modeling the apparent accelerated expansion of the Universe, and about the `weak singularity' are clarified: 1. Of the five definitions of the deceleration parameter given by Hirata and Seljak (HS), only $q_1$ is a correct invariant measure of acceleration/deceleration of expansion. The $q_3$ and $q_4$ are unrelated to acceleration in an inhomogeneous model. 2. The averaging over directions involved in the definition of $q_4$ does not correspond to what is done in observational astronomy. 3. HS's equation (38) connecting $q_4$ to the flow invariants gives self-contradictory results when applied at the centre of symmetry of the Lema\^{\i}tre-Tolman (L-T) model. The intermediate equation (31) that determines $q_{3'}$ is correct, but approximate, so it cannot be used for determining the sign of the deceleration parameter. Even so, at the centre of symmetry of the L-T model, it puts no limitation on the sign of $q_{3'}(0)$. 4. The `weak singularity' of Vanderveld {\it et al.} is a conical profile of mass density at the centre - a perfectly acceptable configuration. 5. The so-called `critical point' in the equations of the `inverse problem' for a central observer in an L-T model is a manifestation of the apparent horizon - a common property of the past light cones in zero-lambda L-T models, perfectly manageable if the equations are correctly integrated.

We present the Bayesian analysis of four different types of backreation models, which are based on the Buchert equations. In this approach, one considers a solution to the Einstein equations for a general matter distribution and then an average of various observable quantities is taken. Such an approach became of considerable interest when it was shown that it could lead to agreement with observations without resorting to dark energy. In this paper we compare the LambdaCDM model and the backreation models with SNIa, BAO, and CMB data, and find that the former is favoured. However, the tested models were based on some particular assumptions about the relation between the average spatial curvature and the backreaction, as well as the relation between the curvature and curvature index. In this paper we modified the latter assumption, leaving the former unchanged. We find that, by varying the relation between the curvature and curvature index, we can obtain a better fit. Therefore, some further work is still needed -- in particular the relation between the backreaction and the curvature should be revisited in order to fully determine the feasibility of the backreaction models to mimic dark energy.

Inhomogeneous cosmological models have had significant success in explaining cosmological observations without the need for dark energy. Generally, these models imply inhomogeneous matter distributions alter the observable relations that are taken for granted when assuming the Universe evolves according to the standard Friedmann equations. Moreover, it has recently been shown that both inhomogeneous matter and pressure distributions are required in both early and late stages of cosmological evolution. These associated pressure gradients are required in the early Universe to sufficiently describe void formation, whilst late-stage pressure gradients stop the appearance of anomalous singularities. In this paper we investigate the effect of pressure gradients on cosmological observations by deriving the luminosity distance - redshift relations in spherically symmetric, inhomogeneous spacetimes endowed with a perfect fluid. By applying this to a specific example for the energy density distribution and using various equations of state, we are able to explicitly show that pressure gradients may have a non-negligble effect on cosmological observations. In particular, we show that a non-zero pressure gradient can imply significantly different residual Hubble diagrams for $z\lesssim1$ compared to when the pressure is ignored. This paper therefore highlights the need to properly consider pressure gradients when interpreting cosmological observations.

Observations of distances to Type-Ia supernovae can be explained by cosmological models that include either a gigaparsec-scale void, or a cosmic flow, without the need for Dark Energy. Instead of invoking dark energy, these inhomogeneous models instead violate the Copernican Principle. we show that current cosmological observations (Supernovae, Baryon Acoustic Oscillations and estimates of the Hubble parameters based on the age of the oldest stars) are not able to rule out inhomogeneous anti-Copernican models. The next generation of surveys for baryonic acoustic oscillations will be sufficiently precise to either validate the Copernican Principle or determine the existence of a local Gpc scale inhomogeneity.

This paper presents the application of the Szekeres Swiss Cheese model to the analysis of observations of the cosmic microwave background (CMB) radiation. The impact of inhomogeneous matter distribution on the CMB observations is in most cases studied within the linear perturbations of the Friedmann model. However, since the density contrast and the Weyl curvature within the cosmic structures are large, this issue is worth studying using another approach. The Szekeres model is an inhomogeneous, non-symmetrical and exact solution of the Einstein equations. In this model, light propagation and matter evolution can be exactly calculated, without such approximations as small amplitude of the density contrast. This allows to examine in more realistic manner the contribution of the light propagation effect to the measured CMB temperature fluctuations. The results of such analysis show that small-scale, non-linear inhomogeneities induce, via Rees-Sciama effect, temperature fluctuations of amplitude $10^{-7} - 10^{-5}$ on angular scale $\theta< 0.24^{\circ}$ ($\ell > 750$). This is still much smaller than the measured temperature fluctuations on this angular scale.However, local and uncompensated inhomogeneities can induce temperature fluctuations of amplitude as large as $10^{-3}$, and thus can be responsible the low multipoles anomalies observed in the angular CMB power spectrum.

This paper considers the volume averaging in the quasispherical Szekeres model. The volume averaging became of considerable interest after it was shown that the volume acceleration calculated within the averaging framework can be positive even though the local expansion rate is always decelerating. This issue was intensively studied within spherically symmetric models. However, since our Universe is not spherically symmetric similar analysis is needed in non symmetrical models. This papers presents the averaging analysis within the quasispherical Szekeres model which is a non-symmetrical generalisation of the spherically symmetric Lema\^itre--Tolman family of models. Density distribution in the quasispherical Szekeres has a structure of a time-dependent mass dipole superposed on a monopole. This paper shows that when calculating the volume acceleration, $\ddot{a}$, within the Szekeres model, the dipole does not contribute to the final result, hence $\ddot{a}$ only depends on a monopole configuration. Thus, the volume averaging within the Szekeres model leads to literally the same solutions as obtained within the Lema\^itre--Tolman model.

In this paper we consider the relation between the volume deceleration parameter obtained within the Buchert averaging scheme and the deceleration parameter derived from the supernova observation. This work was motivated by recent findings that showed that there are models which despite $\Lambda=0$ have volume deceleration parameter $q^{vol} < 0$. This opens the possibility that backreaction and averaging effects may be used as an interesting alternative explanation to the dark energy phenomenon. We have calculated $q^{vol}$ in some Lema\^itre--Tolman models. For those models which are chosen to be realistic and which fit the supernova data, we find that $q^{vol} > 0$, while those models which we have been able to find which exhibit $q^{vol} < 0$ turn out to be unrealistic. This indicates that care must be exercised in relating the deceleration parameter to observations.

Inhomogeneous cosmological models have recently become a very interesting alternative to standard cosmology. This is because these models are able to fit cosmological observations without the need for dark energy. However, due to inhomogeneity and pressure-less matter content, these models can suffer from shell crossing singularities. These singularities occur when two shell of dust collide with each other leading to infinite values of the density. In this Letter we show that if inhomogeneous pressure is included then these singularities can be prevented from occurring over the period of structure formation. Thus, a simple incorporation of a gradient of pressure allows for more comprehensive studies of inhomogeneous cosmological models and their application to cosmology.

This paper analyses cosmological observations within inhomogeneous and exact solutions of the Einstein equations. In some way the analyses presented here can be freed from assumptions such as small amplitude of the density contrast. The supernova observations are analysed using the Lema\itre-Tolman model and the CMB observations are analysed using the quasispherical Szekeres model. The results show that it is possible to fit the supernova data without the cosmological constant. However if inhomogeneities of sizes and amplitudes as observed in the local Universe are considered, their impact on cosmological observations is small.

The aim of this paper is to check if the models with realistic inhomogeneous matter distribution and without cosmological constant can explain the dimming of the supernovae in such a way that it can be interpreted as an acceleration of the Universe. Employing the simplest inhomogeneous model, i.e. Lemaitre-Tolman model, this paper examines the impact of inhomogeneous matter distribution on light propagation. These analyses show that realistic matter fluctuations on small scales induce brightness fluctuations in the residual Hubble diagram of amplitude around 0.15 mag, and thus can mimic acceleration. However, it is different on large scales. All these brightness fluctuations decrease with distance and hence cannot explain the dimmining of supernovae for high redshift without without invoking the cosmological constant.This paper concludes that models with realistic matter distribution (i.e. where variation of the density contrast is similar to what is observed in the local Universe) cannot explain the observed dimming of supernovae without the cosmological constant.

This paper investigates the evolution of a void and an adjacent galaxy supercluster. For this purpose the quasispherical Szekeres model is employed. The Szekeres model is an exact solution of the Einstein field equations. In this way investigations of the evolution of the cosmic structures presented here can be freed from such assumptions as small amplitude of the density contrast. Studying the evolution of a void and an adjacent supercluster results not only in better understanding the evolution of cosmic structures but also in acquainting us with the Szekeres model. The main results include the conclusion that small voids surrounded by large overdensities evolve slower than large, isolated voids do. On the other hand, large voids enhance the evolution of adjacent superclusters which evolve much faster than isolated galaxy superclusters.

This paper investigates evolution of cosmic structures in different environments. For this purpose the quasispherical Szekeres model is employed. The Szekeres model is an exact solution of the Einstein field equations within which it is possible to describe more than one structure. In this way investigations of the evolution of the cosmic structures presented here can be freed from such assumptions as small value of the density contrast. Also, studying the evolution of two or three structures within one framework enables us to follow the interaction between these structures and their impact on the evolution. Main findings include a conclusion that small voids surrounded by large overdensities evolve slower than large, isolated voids do. On the other hand, large voids enhance the evolution of adjacent galaxy superclusters which evolve much faster than isolated superclusters

This paper presents the cosmological applications of the quasispherical Szekeres model. The quasispherical Szekeres model is an exact solution of the Einstein field equations, which represents a time-dependent mass dipole superposed on a monopole and therefore is suitable for modelling double structures such as voids and adjourning galaxy superclusters. Moreover, as the Szekeres model is an exact solution of the Einstein equations it enables tracing light and estimation of the impact of cosmic structures on light propagation. This paper presents the evolution of a void and adjourning supercluster and also reports on how the Szekeres model might be employed either for the estimation of mass of galaxies clusters or for the estimation of the luminosity distance.

Structure formation in the Szekeres model is investigated. Since the Szekeres model is an inhomogeneous model with no symmetries, it is possible to examine the interaction of neighboring structures and its impact on the growth of a density contrast. It has been found that the mass flow from voids to clusters enhances the growth of the density contrast. In the model presented here, the growth of the density contrast is almost 8 times faster than in the linear approach.

This paper aims to investigate the influence of inhomogeneous radiation in void formation. Since the process of void formation is non-linear, a fully relativistic model, which simulates the evolution of voids from the moment of the last scattering until the present instant, is presented. It is found that in order to obtain a model of a void which evolves from $\delta_0 = - 10^{-5}$ at the last scattering moment to the present day $\delta = - 0.94$, the existence of radiation must be taken into account. The ratio of radiation energy density to matter energy density in CDM models at the moment of last scattering is 1/5. This paper proves that such a value of radiation energy density cannot be neglected and influences the first stages of void evolution. Namely, it is important to the process of structure formation and hence significantly influences the dynamics of the Universe in the first millions of years after the last scattering. From the fact that the evolution of voids proceeds differently in various cosmological background models we use the process of void formation to put some limits on values of cosmological parameters. We find that the model with $\Omega_{mat} \approx 0.3$ fits the observational data best.

We develop a non-linear relativistic model of the Shapley Concentration (SC) and its environs, including the Great Attractor (GA) and the Local Group (LG). We take the Shapley concentration as a major attractive centre, and we use the Lemaitre-Tolman model. We constrain our model with present day observations, plus the requirement that it have a physically reasonable evolution from small perturbations at last scattering. We investigate possible mass and velocity distributions, and we find that the peculiar velocity maximum near the SC is \~800 km/s inwards, the density between GA and SC must be about ~0.9 times background, the mass of the GA is probably 4 - 6 x 10^15 M_sun, the SC's contribution to the LG motion is negligible, and the value of the cosmological constant is not significant on this scale.

We develop models of void formation starting from a small initial fluctuation at recombination and growing to a realistic present day density profile in agreement with observations of voids. The model construction is an extension of previously developed algorithms for finding a Lemaitre-Tolman metric that evolves between two profiles of either density or velocity specified at two times. Of the 4 profiles of concern -- those of density and velocity at recombination and at the present day -- two can be specified and the other two follow from the derived model. We find that, in order to reproduce the present-day void density profiles, the initial velocity profile is more important than the initial density profile. Extrapolation of current CMB observations to the scales relevant to proto-voids is very uncertain. Even so, we find that it is very difficult to make both the initial density and velocity fluctuation amplitudes small enough, and still obtain a realistic void by today.