Minkowski functionals and Shapefinders shed light on the connectedness of large-scale structure by determining its topology and morphology. We use a sophisticated code, SURFGEN2, to measure the Minkowski functionals and Shapefinders of individual clusters by modelling cluster surfaces using the 'Marching Cube 33' triangulation algorithm. In this paper, we study the morphology of simulated neutral hydrogen (HI) density fields using Shapefinders at various stages of reionization from the excursion set approach. Accompanying the Shapefinders, we also employ the 'largest cluster statistic' (LCS) to understand the percolation process. Percolation curves demonstrate that the non-Gaussianity in the HI field increases as reionization progresses. The large clusters in both the HI overdense and underdense excursion sets possess similar values of "thickness" ($T$), as well as "breadth" ($B$), but their third Shapefinder - "length" ($L$) - becomes almost proportional to their volume. The large clusters in both HI overdense and underdense segments are overwhelmingly filamentary. The 'cross-section' of a filamentary cluster can be estimated using the product of the first two Shapefinders, $T \times B$. Hence the cross sections of the large clusters at the onset of percolation do not vary much with volume and their sizes only differ in terms of their lengths. This feature appears more vividly in HI overdense regions than in underdense regions and is more pronounced at lower redshifts which correspond to an advanced stage of reionization.

Scalar perturbations during inflation can be substantially amplified by tiny features in the inflaton potential. A bump-like feature behaves like a local speed-breaker and lowers the speed of the scalar field, thereby locally enhancing the scalar power spectrum. A bump-like feature emerges naturally if the base inflaton potential $V_b(\phi)$ contains a local correction term such as $V_b(\phi)\left[1+\varepsilon(\phi)\right]$ at $\phi=\phi_0$. The presence of such a localised correction term at $\phi_0$ leads to a large peak in the curvature power spectrum and to an enhanced probability of black hole formation. Remarkably this does not significantly affect the scalar spectral index $n_{_S}$ and tensor to scalar ratio $r$ on CMB scales. Consequently such models can produce higher mass primordial black holes ($M_{\rm PBH}\geq 1 M_{\odot}$) in contrast to models with `near inflection-point potentials' in which generating higher mass black holes severely affects $n_{_S}$ and $r$. With a suitable choice of the base potential - such as the string theory based (KKLT) inflation or the $\alpha$-attractor models - the amplification of primordial scalar power spectrum can be as large as $10^7$ which leads to a significant contribution of primordial black holes (PBHs) to the dark matter density today, $f_{\rm PBH} = \Omega_{0,\rm PBH}/\Omega_{0,\rm DM} \sim O(1)$. Interestingly, our results remain valid if the bump is replaced by a dip. In this case the base inflaton potential $V_b(\phi)$ contains a negative local correction term such as $V_b(\phi)\left[1-\varepsilon(\phi)\right]$ at $\phi=\phi_0$ which leads to an enhanced probability of PBH formation. We conclude that primordial black holes in the mass range $10^{-17} M_{\odot} \leq M_{\rm PBH} \leq 100\, M_{\odot}$ can easily form in single field inflation in the presence of small bump-like and dip-like features in the inflaton potential.

We investigate constraints on some key cosmological parameters by confronting metastable dark energy models with different combinations of the most recent cosmological observations. Along with the standard $\Lambda$CDM model, two phenomenological metastable dark energy models are considered: (\romannumeral1) DE decays exponentially, (\romannumeral2) DE decays into dark matter. We find that: (1) when considering the most recent supernovae and BAO data, and assuming a fiducial $\Lambda$CDM model, the inconsistency in the estimated value of the $\Omega_{\rm{m,0}}h^2$ parameter obtained by either including or excluding Planck CMB data becomes very much substantial and points to a clear tension~\citep{sahni2014model,zhao2017dynamical}; (2) although the two metastable dark energy models that we study provide greater flexibility in fitting the data, and they indeed fit the SNe Ia+BAO data substantially better than $\Lambda$CDM, they are not able to alleviate this tension significantly when CMB data are included; (3) while local measurements of the Hubble constant are significantly higher relative to the estimated value of $H_0$ in our models (obtained by fitting to SNe Ia and BAO data), the situation seems to be rather complicated with hints of inconsistency among different observational data sets (CMB, SNe Ia+BAO and local $H_0$ measurements). Our results indicate that we might not be able to remove the current tensions among different cosmological observations by considering simple modifications of the standard model or by introducing minimal dark energy models. A complicated form of expansion history, different systematics in different data and/or a non-conventional model of the early Universe might be responsible for these tensions.

We examine the class of initial conditions which give rise to inflation. Our analysis is carried out for several popular models including: Higgs inflation, Starobinsky inflation, chaotic inflation, axion monodromy inflation and non-canonical inflation. In each case we determine the set of initial conditions which give rise to sufficient inflation, with at least $60$ e-foldings. A phase-space analysis has been performed for each of these models and the effect of the initial inflationary energy scale on inflation has been studied numerically. This paper discusses two scenarios of Higgs inflation: (i) the Higgs is coupled to the scalar curvature, (ii) the Higgs Lagrangian contains a non-canonical kinetic term. In both cases we find Higgs inflation to be very robust since it can arise for a large class of initial conditions. One of the central results of our analysis is that, for plateau-like potentials associated with the Higgs and Starobinsky models, inflation can be realised even for initial scalar field values which lie close to the minimum of the potential. This dispels a misconception relating to plateau potentials prevailing in the literature. We also find that inflation in all models is more robust for larger values of the initial energy scale.

Using Shapefinders, which are ratios of Minkowski functionals, we study the morphology of neutral hydrogen (HI) density fields, simulated using semi-numerical technique (inside-out), at various stages of reionization. Accompanying the Shapefinders, we also employ the 'largest cluster statistic' (LCS), originally proposed in Klypin and Shandarin (1993), to study the percolation in both neutral and ionized hydrogen. We find that the largest ionized region is percolating below the neutral fraction $x_{HI} \lesssim 0.728$ (or equivalently $z \lesssim 9$). The study of Shapefinders reveals that the largest ionized region starts to become highly filamentary with non-trivial topology near the percolation transition. During the percolation transition, the first two Shapefinders - 'thickness' ($T$) and 'breadth' ($B$) - of the largest ionized region do not vary much, while the third Shapefinder - 'length' ($L$) - abruptly increases. Consequently, the largest ionized region tends to be highly filamentary and topologically quite complex. The product of the first two Shapefinders, $T\times B$, provides a measure of the 'cross-section' of a filament-like ionized region. We find that, near percolation, the value of $T\times B$ for the largest ionized region remains stable at $\sim 7$ Mpc$^2$ (in comoving scale) while its length increases with time. Interestingly all large ionized regions have similar cross-sections. However their length shows a power-law dependence on their volume, $L\propto V^{0.72}$, at the onset of percolation.

We derive an analytical expression for the growth rate of matter density perturbations on the phantom brane (which is the normal branch of the Dvali-Gabadadze-Porrati model). This model is characterized by a phantomlike effective equation of state for dark energy at the present epoch. It agrees very well with observations. We demonstrate that the traditional parametrization $f=\Omega_m^\gamma$ with a quasiconstant growth index $\gamma$ is not successful in this case. Based on a power series expansion at large redshifts, we propose a different parametrization for this model: $f=\Omega_m^\gamma\left(1+\frac{b}{\ell H}\right)^\beta$, where $\beta$ and $b$ are constants. Our numerical simulations demonstrate that this new parametrization describes the growth rate with great accuracy - the maximum error being $\leq 0.1\%$ for parameter values consistent with observations.

We describe a new class of dark energy (DE) models which behave like cosmological trackers at early times. These models are based on the $\alpha$-attractor set of potentials, originally discussed in the context of inflation. The new models allow the current acceleration of the universe to be reached from a wide class of initial conditions. Prominent examples of this class of models are the potentials $\coth\varphi$ and $\cosh\varphi$. A remarkable feature of this new class of models is that they lead to large enough negative values of the equation of state at the present epoch, consistent with the observations of accelerated expansion of the universe, from a very large initial basin of attraction. They therefore avoid the fine tuning problem which afflicts many models of DE.

In [Schmidt PRD 80 123003 (2009)], the author suggested that dynamical dark energy (DDE) propagating on the phantom brane could mimick $\Lambda$CDM. Schmidt went on to derive a phenomenological expression for $\rho_{\rm DE}$ which could achieve this. We demonstrate that while Schmidt's central premise is correct, the expression for $\rho_{\rm DE}$ derived in Schmidt (2009) is flawed. We derive the correct expression for $\rho_{\rm DE}$ which leads to $\Lambda$CDM-like expansion on the phantom brane. We also show that DDE on the brane can be associated with a Quintessence field and derive a closed form expression for its potential $V(\phi)$. Interestingly the $\alpha$-attractor based potential $V(\phi) \propto \coth^2{\lambda\phi}$ makes braneworld expansion resemble $\Lambda$CDM. However the two models can easily be distinguished on the basis of density perturbations which grow at different rates on the braneworld and in $\Lambda$CDM.

Non-canonical scalar fields with the Lagrangian ${\cal L} = X^\alpha - V(\phi)$, possess the attractive property that the speed of sound, $c_s^{2} = (2\,\alpha - 1)^{-1}$, can be exceedingly small for large values of $\alpha$. This allows a non-canonical field to cluster and behave like warm/cold dark matter on small scales. We demonstrate that simple potentials including $V = V_0\coth^2{\phi}$ and a Starobinsky-type potential can unify dark matter and dark energy. Cascading dark energy, in which the potential cascades to lower values in a series of discrete steps, can also work as a unified model. In all of these models the kinetic term $X^\alpha$ plays the role of dark matter, while the potential term $V(\phi)$ plays the role of dark energy.

We propose a new class of metastable dark energy (DE) phenomenological models in which the DE decay rate does not depend on external parameters such as the scale factor or the curvature of the Universe. Instead, the DE decay rate is assumed to be a constant depending only on intrinsic properties of DE and the type of a decay channel, similar to case of the radioactive decay of unstable particles and nuclei. As a consequence, the DE energy density becomes a function of the proper time elapsed since its formation, presumably in the very early Universe. Such a natural type of DE decay can profoundly affect the expansion history of the Universe and its age. Metastable DE can decay in three distinct ways: (i) exponentially, (ii) into dark matter, (iii) into dark radiation. Testing metastable DE models with observational data we find that the decay half-life must be many times larger than the age of the Universe. Models in which dark energy decays into dark matter lead to lower values of the Hubble parameter at large redshifts relative to $\Lambda$CDM. Consequently these models provide a better fit to cosmological BAO data (especially data from from high redshift quasars) than concordance ($\Lambda$CDM) cosmology.

Recently, Kallosh and Linde have drawn attention to a new family of superconformal inflationary potentials, subsequently called $\alpha$-attractors. The $\alpha$-attractor family can interpolate between a large class of inflationary models. It also has an important theoretical underpinning within the framework of supergravity. We demonstrate that the $\alpha$-attractors have an even wider appeal since they may describe dark matter and perhaps even dark energy. The dark matter associated with the $\alpha$-attractors, which we call $\alpha$-dark matter ($\alpha$DM), shares many of the attractive features of fuzzy dark matter, with $V(\varphi) = \frac{1}{2}m^2\varphi^2$, while having none of its drawbacks. Like fuzzy dark matter, $\alpha$DM can have a large Jeans length which could resolve the cusp-core and substructure problems faced by standard cold dark matter. $\alpha$DM also has an appealing tracker property which enables it to converge to the late-time dark matter asymptote, $\langle w\rangle \simeq 0$, from a wide range of initial conditions. It thus avoids the enormous fine-tuning problems faced by the $m^2\varphi^2$ potential in describing dark matter.

The phantom brane has several important distinctive features: (i) Its equation of state is phantom-like, but there is no future `big rip' singularity, (ii) the effective cosmological constant on the brane is dynamically screened, because of which the expansion rate is {\em smaller} than that in $\Lambda$CDM at high redshifts. In this paper, we constrain the Phantom braneworld using distance measures such as Type Ia supernovae (SNeIa), Baryon Acoustic Oscillations (BAO), and the compressed Cosmic Microwave Background (CMB) data. We find that the simplest braneworld models provide a good fit to the data. For instance, BAO +SNeIa data can be accommodated by the braneworld for a large region in parameter space $0 < \Omega_l < 0.3$ at $1\sigma$. The Hubble parameter can be as high as $H_0 < 78$ km/s/Mpc, and the effective equation of state at present can show phantom-like behaviour with $w_0 < -1.2$ at $1\sigma$. We note a correlation between $H_0$ and $w_0$, with higher values of $H_0$ leading to a lower, and more phantom-like, value of $w_0$. Inclusion of CMB data provides tighter constraints $\Omega_l < 0.1$. (Here $\Omega_l$ encodes the ratio of the five and four dimensional Planck mass.) The Hubble parameter in this case is more tightly constrained to $H_0 < 71$ km/s/Mpc, and the effective equation of state to $w_0 < -1.1$. Interestingly, we find that the universe is allowed be closed or open, with $-0.5 < \Omega_{\kappa} < 0.5$, even on including the compressed CMB data. There appears to be some tension in the low and high $z$ BAO data which may either be resolved by future data, or act as a pointer to interesting new cosmology.

We obtain a closed system of equations for scalar perturbations in a multi-component braneworld. Our braneworld possesses a phantom-like equation of state at late times, $w_{\rm eff} < -1$, but no big-rip future singularity. In addition to matter and radiation, the braneworld possesses a new effective degree of freedom - the 'Weyl fluid' or 'dark radiation'. Setting initial conditions on super-Hubble spatial scales at the epoch of radiation domination, we evolve perturbations of radiation, pressureless matter and the Weyl fluid until the present epoch. We observe a gradual decrease in the amplitude of the Weyl-fluid perturbations after Hubble-radius crossing, which results in a negligible effect of the Weyl fluid on the evolution of matter perturbations on spatial scales relevant for structure formation. Consequently, the quasi-static approximation of Koyama and Maartens provides a good fit to the exact results during the matter-dominated epoch. We find that the late-time growth of density perturbations on the brane proceeds at a faster rate than in $\Lambda$CDM. Additionally, the gravitational potentials $\Phi$ and $\Psi$ evolve differently on the brane than in $\Lambda$CDM, for which $\Phi = \Psi$. On the brane, by contrast, the ratio $\Phi/\Psi$ exceeds unity during the late matter-dominated epoch ($z \lesssim 50$). These features emerge as "smoking gun" tests of phantom brane cosmology and allow predictions of this scenario to be tested against observations of galaxy clustering and large-scale structure.

It is well known that a canonical scalar field is able to describe either dark matter or dark energy but not both. We demonstrate that a non-canonical scalar field can describe both dark matter and dark energy within a unified setting. We consider the simplest extension of the canonical Lagrangian ${\cal L} \propto X^\alpha - V(\phi)$ where $\alpha \geq 1$ and $V$ is a sufficiently flat potential. In this case the kinetic term in the Lagrangian behaves just like a perfect fluid, whereas the potential term mimicks the cosmological constant. For very large values, $\alpha \gg 1$, the equation of state of the kinetic term drops to zero and the expansion rate of the universe mimicks $\Lambda$CDM. The velocity of sound in this model, and the associated gravitational clustering, is sensitive to the value of $\alpha$. For very large values of $\alpha$ the clustering properties of our model resemble those of cold dark matter (CDM). But for smaller values of $\alpha$, gravitational clustering on small scales is suppressed, and our model has properties resembling those of warm dark matter (WDM). Therefore our non-canonical model has an interesting new property: while the background universe expands like $\Lambda$CDM, its clustering properties can resemble those of either cold or warm dark matter.

It is generally believed that a cosmological arrow of time must be associated with entropy production. Indeed, in his seminal work on cyclic cosmology, Tolman introduced a viscous fluid in order to make successive expansion/contraction cycles larger than previous ones, thereby generating an arrow of time. However, as we demonstrate in this letter, the production of entropy is not the only means by which a cosmological arrow of time may emerge. Remarkably, systems which are dissipationless may nevertheless demonstrate a preferred direction of time provided they possess attractors. An example of a system with well defined attractors is scalar-field driven cosmology. In this case, for a wide class of potentials (especially those responsible for inflation), the attractor equation of state during expansion can have the form $p \simeq -\rho$, and during contraction $p \simeq \rho$. If the resulting cosmology is cyclic, then the presence of cosmological hysteresis, $\oint p~dV \neq 0$ during successive cycles, causes an arrow of time to emerge in a system which is formally dissipationless. An important analogy is drawn between the arrow of time in cyclic cosmology and an arrow of time in an $N$-body system of gravitationally interacting particles. We find that, like the $N$-body system, a cyclic universe can evolve from a single past into two futures with oppositely directed arrows of time.

Solar models suggest that four billion years ago the young Sun was about 25% fainter than it is today, rendering Earth's oceans frozen and lifeless. However, there is ample geophysical evidence that Earth had a liquid ocean teeming with life 4 Gyr ago. Since ${\cal L_\odot} \propto G^7M_\odot^5$, the Sun's luminosity ${\cal L_\odot}$ is exceedingly sensitive to small changes in the gravitational constant $G$. We show that a percent-level increase in $G$ in the past would have prevented Earth's oceans from freezing, resolving the faint young Sun paradox. Such small changes in $G$ are consistent with observational bounds on ${\Delta G}/G$. Since ${\cal L}_{\rm SNIa} \propto G^{-3/2}$, an increase in $G$ leads to fainter supernovae, creating tension between standard candle and standard ruler probes of dark energy. Precisely such a tension has recently been reported by the Planck team.

Baryon Acoustic Oscillations (BAO) allow us to determine the expansion history of the Universe, thereby shedding light on the nature of dark energy. Recent observations of BAO's in the SDSS DR9 and DR11 have provided us with statistically independent measurements of $H(z)$ at redshifts of 0.57 and 2.34, respectively. We show that these measurements can be used to test the cosmological constant hypothesis in a model independent manner by means of an improved version of the $Om$ diagnostic. Our results indicate that the SDSS DR11 measurement of $H(z) = 222 \pm 7$ km/sec/Mpc at $z = 2.34$, when taken in tandem with measurements of $H(z)$ at lower redshifts, imply considerable tension with the standard $\Lambda$CDM model. Our estimation of the new diagnostic $Omh^2$ from SDSS DR9 and DR11 data, namely $Omh^2 \approx 0.122 \pm 0.01$, which is equivalent to $\Omega_{0m}h^2$ for the spatially flat $\Lambda$CDM model, is in tension with the value $\Omega_{0m}h^2 = 0.1426 \pm 0.0025$ determined for $\Lambda$CDM from Planck+WP. This tension is alleviated in models in which the cosmological constant was dynamically screened (compensated) in the past. Such evolving dark energy models display a pole in the effective equation of state of dark energy at high redshifts, which emerges as a smoking gun test for these theories.

We explore the possibility of emergent cosmology using the effective potential formalism. We discover new models of emergent cosmology which satisfy the constraints posed by the cosmic microwave background (CMB). We demonstrate that, within the framework of modified gravity, the emergent scenario can arise in a universe which is spatially open/closed. By contrast, in general relativity (GR) emergent cosmology arises from a spatially closed past-eternal Einstein Static Universe (ESU). In GR the ESU is unstable, which creates fine tuning problems for emergent cosmology. However, modified gravity models including Braneworld models, Loop Quantum Cosmology (LQC) and Asymptotically Free Gravity result in a stable ESU. Consequently, in these models emergent cosmology arises from a larger class of initial conditions including those in which the universe eternally oscillates about the ESU fixed point. We demonstrate that such an oscillating universe is necessarily accompanied by graviton production. For a large region in parameter space graviton production is enhanced through a parametric resonance, casting serious doubts as to whether this emergent scenario can be past-eternal.

The Javalambre-Physics of the Accelerated Universe Astrophysical Survey (J-PAS) is a narrow band, very wide field Cosmological Survey to be carried out from the Javalambre Observatory in Spain with a purpose-built, dedicated 2.5m telescope and a 4.7 sq.deg. camera with 1.2Gpix. Starting in late 2015, J-PAS will observe 8500sq.deg. of Northern Sky and measure $0.003(1+z)$ photo-z for $9\times10^7$ LRG and ELG galaxies plus several million QSOs, sampling an effective volume of $\sim 14$ Gpc$^3$ up to $z=1.3$ and becoming the first radial BAO experiment to reach Stage IV. J-PAS will detect $7\times 10^5$ galaxy clusters and groups, setting constrains on Dark Energy which rival those obtained from its BAO measurements. Thanks to the superb characteristics of the site (seeing ~0.7 arcsec), J-PAS is expected to obtain a deep, sub-arcsec image of the Northern sky, which combined with its unique photo-z precision will produce one of the most powerful cosmological lensing surveys before the arrival of Euclid. J-PAS unprecedented spectral time domain information will enable a self-contained SN survey that, without the need for external spectroscopic follow-up, will detect, classify and measure $\sigma_z\sim 0.5\%$ redshifts for $\sim 4000$ SNeIa and $\sim 900$ core-collapse SNe. The key to the J-PAS potential is its innovative approach: a contiguous system of 54 filters with $145\AA$ width, placed $100\AA$ apart over a multi-degree FoV is a powerful "redshift machine", with the survey speed of a 4000 multiplexing low resolution spectrograph, but many times cheaper and much faster to build. The J-PAS camera is equivalent to a 4.7 sq.deg. "IFU" and it will produce a time-resolved, 3D image of the Northern Sky with a very wide range of Astrophysical applications in Galaxy Evolution, the nearby Universe and the study of resolved stellar populations.

Ya.B. Zeldovich was a pre-eminent Soviet physicist whose seminal contributions spanned many fields ranging from physical chemistry to nuclear and particle physics, and finally astrophysics and cosmology. March 8, 2014 marks Zeldovich's birth centenary, and this article attempts to convey the zest with which Zeldovich did science, and the important role he played in fostering and mentoring a whole generation of talented Scientists.

It is well known that a canonical scalar field with an exponential potential can drive power law inflation (PLI). However, the tensor-to-scalar ratio in such models turns out to be larger than the stringent limit set by recent Planck results. We propose a new model of power law inflation for which the scalar spectra index, the tensor-to-scalar ratio and the non-gaussianity parameter $f_{_{\mathbf{NL}}}^{\mathrm{equil}}$ are in excellent agreement with Planck results. Inflation, in this model, is driven by a non-canonical scalar field with an inverse power law potential. The Lagrangian for our model is structurally similar to that of a canonical scalar field and has a power law form for the kinetic term. A simple extension of our model resolves the graceful exit problem which usually afflicts models of power law inflation.

We propose a novel approach to the problem of cosmological perturbations in a braneworld model with induced gravity, which leads to a closed system of equations on the brane. We focus on a spatially closed brane that bounds the interior four-ball of the bulk space. The background cosmological evolution on the brane is now described by the normal branch, and the boundary conditions in the bulk become the regularity conditions for the metric everywhere inside the four-ball. In this approach, there is no spatial infinity or any other boundary in the bulk space since the spatial section is compact, hence, we term this setup as a no-boundary proposal. Assuming that the bulk cosmological constant is absent and employing the Mukohyama master variable, we argue that the effects of nonlocality on brane perturbations may be ignored if the brane is marginally closed. In this case, there arises a relation that closes the system of equations for perturbations on the brane. Perturbations of pressureless matter and dark radiation can now be described by a system of coupled second-order differential equations. Remarkably, this system can be exactly solved in the matter-dominated and de Sitter regimes. In this case, apart from the usual growing and decaying modes, we find two additional modes that behave monotonically on super-Hubble spatial scales and exhibit rapid oscillations with decaying amplitude on sub-Hubble spatial scales.

We compare two different probes of the expansion history of the universe, namely, luminosity distances from type Ia supernovae and angular diameter distances from galaxy clusters, using the Bayesian interpretation of Crossing statistic [1, 2] in conjunction with the assumption of cosmic duality relation. Our analysis is conducted independently of any a-priori assumptions about the nature of dark energy. The model independent method which we invoke searches for inconsistencies between SNIa and galaxy cluster data sets. If detected such an inconsistency would imply the presence of systematics in either of the two data sets. Simulating observations based on expected WFIRST supernovae data and X-ray eROSITA + SZ Planck cluster data, we show that our method allows one to detect systematics with high precision and without advancing any hypothesis about the nature of dark energy.

Future space-based gravity wave experiments such as the Big Bang Observatory (BBO), with their excellent projected, one sigma angular resolution, will measure the luminosity distance to a large number of gravity wave (GW) sources to high precision, and the redshift of the single galaxies in the narrow solid angles towards the sources will provide the redshifts of the gravity wave sources. One sigma BBO beams contain the actual source only in 68 per cent cases; the beams that do not contain the source may contain a spurious single galaxy, leading to misidentification. To increase the probability of the source falling within the beam, larger beams have to be considered, decreasing the chances of finding single galaxies in the beams. Saini, Sethi and Sahni (2010) argued, largely analytically, that identifying even a small number of GW source galaxies furnishes a rough distance-redshift relation, which could be used to further resolve sources that have multiple objects in the angular beam. In this work we further develop this idea by introducing a self-calibrating iterative scheme which works in conjunction with Monte-Carlo simulations to determine the luminosity distance to GW sources with progressively greater accuracy. This iterative scheme allows one to determine the equation of state of dark energy to within an accuracy of a few percent for a gravity wave experiment possessing a beam width an order of magnitude larger than BBO (and therefore having a far poorer angular resolution). This is achieved with no prior information about the nature of dark energy from other data sets such as SN Ia, BAO, CMB etc.

In a series of recent papers, including arXiv:1210.1183, it was claimed that large-scale magnetic fields generated during inflation in a spatially open universe could remain astrophysically significant at the present time since they experienced superadiabatic amplification specific to an open universe. We reexamine this assertion and show that, on the contrary, large-scale magnetic fields in a realistic open universe decay in much the same manner as they would in a spatially flat universe. Consequently, their amplitude today is extremely small (B0 < 10^{-59} G) and is unlikely to be of astrophysical significance.

Baryon Acoustic Oscillations (BAO) provide an important standard ruler which can be used to probe the recent expansion history of our universe. We show how a simple extension of the Om diagnostic, which we call Om3, can combine standard ruler information from BAO with standard candle information from type Ia supernovae (SNIa) to yield a powerful novel null diagnostic of the cosmological constant hypothesis. A unique feature of Om3 is that it requires minimal cosmological assumptions since its determination does not rely upon prior knowledge of either the current value of the matter density and the Hubble constant, or the distance to the last scattering surface. Observational uncertainties in these quantities therefore do not affect the reconstruction of Om3. We reconstruct Om3 using the Union 2.1 SNIa data set and BAO data from SDSS, WiggleZ and 6dFGS. Our results are consistent with dark energy being the cosmological constant. We show how Om and Om3 can be used to obtain accurate model independent constraints on the properties of dark energy from future data sets such as BigBOSS.

This paper revisits the Inflationary scenario within the framework of scalar field models possessing a non-canonical kinetic term. We obtain closed form solutions for all essential quantities associated with chaotic inflation including slow roll parameters, scalar and tensor power spectra, spectral indices, the tensor-to-scalar ratio, etc. We also examine the Hamilton-Jacobi equation and demonstrate the existence of an inflationary attractor. Our results highlight the fact that non-canonical scalars can significantly improve the viability of inflationary models. They accomplish this by decreasing the tensor-to-scalar ratio while simultaneously increasing the value of the scalar spectral index, thereby redeeming models which are incompatible with the cosmic microwave background (CMB) in their canonical version. For instance, the non-canonical version of the chaotic inflationary potential, $V(\phi) \sim \lambda\phi^4$, is found to agree with observations for values of $\lambda$ as large as unity ! The exponential potential can also provide a reasonable fit to CMB observations. A central result of this paper is that {\em steep potentials} (such as $V \propto \phi^{-n}$) usually associated with dark energy, can drive inflation in the non-canonical setting. Interestingly, non-canonical scalars violate the consistency relation $r = -8n_T$, which emerges as a {\em smoking gun} test for this class of models.

A Universe filled with a homogeneous scalar field exhibits `Cosmological hysteresis'. Cosmological hysteresis is caused by the asymmetry in the equation of state during expansion and contraction. This asymmetry results in the formation of a hysteresis loop: $\oint pdV$, whose value can be non-vanishing during each oscillatory cycle. For flat potentials, a negative value of the hysteresis loop leads to the increase in amplitude of consecutive cycles and to a universe with older and larger successive cycles. Such a universe appears to possess an arrow of time even though entropy production is absent and all of the equations respect time-reversal symmetry ! Cosmological hysteresis appears to be widespread and exists for a large class of scalar field potentials and mechanisms for making the universe bounce. For steep potentials, the value of the hysteresis loop can be positive as well as negative. The expansion factor in this case displays quasi-periodic behaviour in which successive cycles can be both larger as well as smaller than previous ones. This quasi-regular pattern resembles the phenomenon of BEATS displayed by acoustic systems. Remarkably, the expression relating the increase/decrease in oscillatory cycles to the quantum of hysteresis appears to be model independent. The cyclic scenario is extended to spatially anisotropic models and it is shown that the anisotropy density decreases during successive cycles if the hysteresis loop is negative.

We show how higher derivatives of the expansion factor can be developed into a null diagnostic for concordance cosmology (LCDM). It is well known that the Statefinder -- the third derivative of the expansion factor written in dimensionless form, a^{(3)}/aH^3, equals unity for LCDM. We generalize this result to higher derivatives of the expansion factor and demonstrate that the hierarchy, a^{(n)}/aH^n, can be converted to a form that stays pegged at unity in concordance cosmology. This remarkable property of the Statefinder hierarchy enables it to be used as an extended null diagnostic for the cosmological constant. The Statefinder hierarchy combined with the growth rate of matter perturbations defines a composite null diagnostic which can distinguish evolving dark energy from LCDM.

We generalize the cosmic energy equation to the case when massive particles interact via a modified gravitational potential of the form phi(a, r), which is allowed to explicitly depend upon the cosmological time through the expansion factor a(t). Using the nonrelativistic approximation for particle dynamics, we derive the equation for the cosmological expansion which has the form of the Friedmann equation with a renormalized gravitational constant. The generalized Layzer-Irvine cosmic energy equation and the associated cosmic virial theorem are applied to some recently proposed modifications of the Newtonian gravitational interaction between dark-matter particles. We also draw attention to the possibility that the cosmic energy equation may be used to probe the expansion history of the universe thereby throwing light on the nature of dark matter and dark energy.

By observing mergers of compact objects, future gravity wave experiments would measure the luminosity distance to a large number of sources to a high precision but not their redshifts. Given the directional sensitivity of an experiment, a fraction of such sources (gold plated -- GP) can be identified optically as single objects in the direction of the source. We show that if an approximate distance-redshift relation is known then it is possible to statistically resolve those sources that have multiple galaxies in the beam. We study the feasibility of using gold plated sources to iteratively resolve the unresolved sources, obtain the self-calibrated best possible distance-redshift relation and provide an analytical expression for the accuracy achievable. We derive lower limit on the total number of sources that is needed to achieve this accuracy through self-calibration. We show that this limit depends exponentially on the beam width and give estimates for various experimental parameters representative of future gravitational wave experiments DECIGO and BBO.

We discuss the cosmological evolution of a braneworld in five dimensional Gauss-Bonnet gravity. Our discussion allows the fifth (bulk) dimension to be space-like as well as time-like. The resulting equations of motion have the form of a cubic equation in the (H^2,(\rho+\sigma)^2) plane, where \sigma is the brane tension and \rho is the matter density. This allows us to conduct a comprehensive pictorial analysis of cosmological evolution for the Gauss-Bonnet brane. The many interesting properties of this braneworld include the possibility of accelerated expansion at late times. For a finite region in parameter space the accelerated expansion can be phantom-like so that w < -1. At late times, this branch approaches de Sitter space (w = -1) and avoids the big-rip singularities usually present in phantom models. For a time-like extra dimension the Gauss-Bonnet brane can bounce and avoid the initial singularity.

We investigate the course of cosmic expansion in its `recent past' using the Constitution SN Ia sample (which includes CfA data at low redshifts), jointly with signatures of baryon acoustic oscillations (BAO) in the galaxy distribution and fluctuations in the cosmic microwave background (CMB). Earlier SN Ia data sets could not address this issue because of a paucity of data at low redshifts. Allowing the equation of state of dark energy (DE) to vary, we find that a coasting model of the universe (q_0=0) fits the data about as well as LCDM. This effect, which is most clearly seen using the recently introduced `Om' diagnostic, corresponds to an increase of Om(z) and q(z) at redshifts z \lleq 0.3. In geometrical terms, this suggests that cosmic acceleration may have already peaked and that we are currently witnessing its slowing down. The case for evolving DE strengthens if a subsample of the Constitution set consisting of SNLS+ESSENCE+CfA SN Ia data is analysed in combination with BAO+CMB using the same statistical methods. The effect we observe could correspond to DE decaying into dark matter (or something else). A toy model which mimics this process agrees well with the combined SN Ia+BAO+CMB data.

For a large class of dark energy (DE) models, for which the effective gravitational constant is a constant and there is no direct exchange of energy between DE and dark matter (DM), knowledge of the expansion history suffices to reconstruct the growth factor of linearized density perturbations in the non-relativistic matter component on scales much smaller than the Hubble distance. In this paper we develop a non-parametric method for extracting information about the perturbative growth factor from data pertaining to the luminosity or angular size distances. A comparison of the reconstructed density contrast with observations of large scale structure and gravitational lensing can help distinguish DE models such as the cosmological constant and quintessence from models based on modified gravity theories as well as models in which DE and DM are either unified, or interact directly. We show that for current SNe data, the linear growth factor at z = 0.3 can be constrained to 5%, and the linear growth rate to 6%. With future SNe data, such as expected from the JDEM mission, we may be able to constrain the growth factor to 2-3% and the growth rate to 3-4% at z = 0.3 with this unbiased, model-independent reconstruction method. For future BAO data which would deliver measurements of both the angular diameter distance and Hubble parameter, it should be possible to constrain the growth factor at z = 2.5 to 9%. These constraints grow tighter with the errors on the datasets. With a large quantity of data expected in the next few years, this method can emerge as a competitive tool for distinguishing between different models of dark energy.

A sudden small change in the second derivative of the inflaton potential can result in a universal local feature in the spectrum of primordial perturbations generated during inflation. The exact solution describing this feature \cite{minu} is characterized by a step in the spectral index modulated by characteristic oscillations and results in a large running of the spectral index localized over a few e-folds of scale. In this paper we confront this step-like feature with the 5 year WMAP results and demonstrate that it provides a better fit to this data than a featureless initial spectrum. If such a feature exists at all, then it should lie at sufficiently large scales $k_0 \lesssim 0.003 {\rm Mpc}^{-1}$ corresponding to $l \lesssim 40$. The sign of the effect is shown to correspond to the negative running of $n_s$ localized near this scale. This feature could arise as a result of a `mini-waterfall'-type fast second order phase transition experienced by an auxiliary heavy field during inflation, in a model similar to hybrid inflation (though for a different choice of parameters). If this is the case, then the auxiliary field should be positively coupled to the inflaton.

We investigate the cosmological properties of an "induced gravity" brane scenario in the absence of mirror symmetry with respect to the brane. We find that brane evolution can proceed along one of four distinct branches. By contrast, when mirror symmetry is imposed, only two branches exist, one of which represents the self-accelerating brane, while the other is the so-called normal branch. This model incorporates many of the well-known possibilities of brane cosmology including phantom acceleration (w < -1), self-acceleration, transient acceleration, quiescent singularities, and cosmic mimicry. Significantly, the absence of mirror symmetry also provides an interesting way of inducing a sufficiently small cosmological constant on the brane. A small (positive) Lambda-term in this case is induced by a small asymmetry in the values of bulk fundamental constants on the two sides of the brane.

Brane cosmology presents many interesting possibilities including: phantom acceleration (w<-1), self-acceleration, unification of dark energy with inflation, transient acceleration, loitering cosmology, new singularities at which the Hubble parameter remains finite, cosmic mimicry, etc. The existence of a time-like extra dimension can result in a singularity-free cyclic cosmology.

We introduce two new diagnostics of dark energy (DE). The first, Om, is a combination of the Hubble parameter and the cosmological redshift and provides a "null test" of dark energy being a cosmological constant. Namely, if the value of Om(z) is the same at different redshifts, then DE is exactly cosmological constant. The slope of Om(z) can differentiate between different models of dark energy even if the value of the matter density is not accurately known. For DE with an unevolving equation of state, a positive slope of Om(z) is suggestive of Phantom (w < -1) while a negative slope indicates Quintessence (w > -1). The second diagnostic, "acceleration probe"(q-probe), is the mean value of the deceleration parameter over a small redshift range. It can be used to determine the cosmological redshift at which the universe began to accelerate, again without reference to the current value of the matter density. We apply the "Om" and "q-probe" diagnostics to the Union data set of type Ia supernovae combined with recent data from the cosmic microwave background (WMAP5) and baryon acoustic oscillations.

A model is developed in which the inflaton potential experiences a sudden small change in its second derivative (the effective mass of the inflaton). An exact treatment demonstrates that the resulting density perturbation has a quasi-flat power spectrum with a break in its slope (a step in n_s). The step in the spectral index is modulated by characteristic oscillations and results in large running of the spectral index localised over a few e-folds of scales. A field-theoretic model giving rise to such behaviour of the inflationary potential is based on a fast phase transition experienced by a second scalar field weakly coupled to the inflaton. Such a transition is similar to that which terminates inflation in the hybrid inflationary scenario. This scenario suggests that the observed running of the spectral index in the WMAP data may be caused by a fast second order phase transition which occurred during inflation.

A proposal is made to test Newton's inverse-square law using the perihelion shift of test masses (planets) in free fall within a spacecraft located at the Earth-Sun L2 point. Such an Artificial Planetary System In Space (APSIS) will operate in a drag-free environment with controlled experimental conditions and minimal interference from terrestrial sources of contamination. We demonstrate that such a space experiment can probe the presence of a "hidden" fifth dimension on the scale of a micron, if the perihelion shift of a "planet" can be measured to sub-arc-second accuracy. Some suggestions for spacecraft design are made.

An outstanding issue in braneworld theory concerns the setting up of proper boundary conditions for the brane-bulk system. Boundary conditions (BC's) employing regulatory branes or demanding that the bulk metric be nonsingular have yet to be implemented in full generality. In this paper, we take a different route and specify boundary conditions directly on the brane thereby arriving at a local and closed system of equations (on the brane). We consider a one-parameter family of boundary conditions involving the anisotropic stress of the projection of the bulk Weyl tensor on the brane and derive an exact system of equations describing scalar cosmological perturbations on a generic braneworld with induced gravity. Depending upon our choice of boundary conditions, perturbations on the brane either grow moderately (region of stability) or rapidly (instability). In the instability region, the evolution of perturbations usually depends upon the scale: small scale perturbations grow much more rapidly than those on larger scales. This instability is caused by a peculiar gravitational interaction between dark radiation and matter on the brane. Generalizing the boundary conditions obtained by Koyama and Maartens, we find for the Dvali-Gabadadze-Porrati model an instability, which leads to a dramatic scale-dependence of the evolution of density perturbations in matter and dark radiation. A different set of BC's, however, leads to a more moderate and scale-independent growth of perturbations. For the mimicry braneworld, which expands like LCDM, this class of BC's can lead to an earlier epoch of structure formation.

We propose a non-parametric method of smoothing supernova data over redshift using a Gaussian kernel in order to reconstruct important cosmological quantities including H(z) and w(z) in a model independent manner. This method is shown to be successful in discriminating between different models of dark energy when the quality of data is commensurate with that expected from the future SuperNova Acceleration Probe (SNAP). We find that the Hubble parameter is especially well-determined and useful for this purpose. The look back time of the universe may also be determined to a very high degree of accuracy (\lleq 0.2 %) in this method. By refining the method, it is also possible to obtain reasonable bounds on the equation of state of dark energy. We explore a new diagnostic of dark energy-- the `w-probe'-- which can be calculated from the first derivative of the data. We find that this diagnostic is reconstructed extremely accurately for different reconstruction methods even if \Omega_m is marginalized over. The w-probe can be used to successfully distinguish between $\Lambda$CDM and other models of dark energy to a high degree of accuracy.

We reconstruct dark energy properties from two complementary supernova datasets -- the newly released Gold+HST sample and SNLS. The results obtained are consistent with standard $\Lambda$CDM model within $2\sigma$ error bars although the Gold+HST data favour evolving dark energy slightly more than SNLS. Using complementary data from baryon acoustic oscillations and the cosmic microwave background to constrain dark energy, we find that our results in this case are strongly dependent on the present value of the matter density $\Omega_m$. Consequently, no firm conclusions regarding constancy or variability of dark energy density can be drawn from these data alone unless the value of $\Omega_m$ is known to an accuracy of a few percent. However, possible variability is significantly restricted if this data is used in conjunction with supernova data.

This review summarizes recent attempts to reconstruct the expansion history of the Universe and to probe the nature of dark energy. Reconstruction methods can be broadly classified into parametric and non-parametric approaches. It is encouraging that, even with the limited observational data currently available, different approaches give consistent results for the reconstruction of the Hubble parameter $H(z)$ and the effective equation of state $w(z)$ of dark energy. Model independent reconstruction using current data allows for modest evolution of dark energy density with redshift. However, a cosmological constant (= dark energy with a constant energy density) remains an excellent fit to the data. Some pitfalls to be guarded against during cosmological reconstruction are summarized and future directions for the model independent reconstruction of dark energy are explored.

We examine a class of braneworld models in which the expanding universe encounters a "quiescent" future singularity. At a quiescent singularity, the energy density and pressure of the cosmic fluid as well as the Hubble parameter remain finite while all derivatives of the Hubble parameter diverge (i.e., ${\dot H}$, ${\ddot H}$, etc. $\to \infty$). Since the Kretschmann invariant diverges ($R_{iklm}R^{iklm} \to \infty$) at the singularity, one expects quantum effects to play an important role as the quiescent singularity is approached. We explore the effects of vacuum polarization due to massless conformally coupled fields near the singularity and show that these can either cause the universe to recollapse or, else, lead to a softer singularity at which $H$, ${\dot H}$, and ${\ddot H}$ remain finite while ${\dddot H}$ and higher derivatives of the Hubble parameter diverge. An important aspect of the quiescent singularity is that it is encountered in regions of low density, which has obvious implications for a universe consisting of a cosmic web of high and low density regions -- superclusters and voids. In addition to vacuum polarization, the effects of quantum particle production of non-conformal fields are also likely to be important. A preliminary examination shows that intense particle production can lead to an accelerating universe whose Hubble parameter shows oscillations about a constant value.

Braneworld cosmology has several attractive and distinctive features. For instance the effective equation of state in braneworld models can be both quintessence-like (w_0 \geq -1) as well as phantom-like (w_0 \leq -1). Models with w_0 \geq -1 (w_0 \leq -1) are referred to as Brane 2 (Brane 1) and correspond to complementary embeddings of the brane in the bulk. (The equation of state in Brane 1 can successfully cross the `phantom divide' at w = -1.) In this paper we compare the predictions of braneworld models to two recently released supernova data sets: the `Gold' data (Riess et al, 2004) and the data from the Supernova Legacy Survey (SNLS) (Astier et al, 2005). We also incorporate the recent discovery of the baryon acoustic peak in the Sloan Digital Sky Survey (Eisenstein et al, 2005) into our analysis. Our main results are that braneworld models satisfy both sets of SNe data. Brane 1 (with w_0 \leq -1) shows very good agreement with data for values of the matter density bounded {\em from below}: \Omega_m \geq 0.25 (Gold) and $\Omega_m \geq 0.2 (SNLS). On the other hand Brane 2 (with w_0 \geq -1) shows excellent agreement with data for values of the matter density which are bounded {\em from above}: \Omega_m \leq 0.45 (Gold) and \Omega_m \leq 0.35 (SNLS). The DGP model is excluded at 3\sigma by SNLS and at 1\sigma by the Gold dataset. Braneworld models with future `quiescent' singularities (at which the Hubble parameter and the matter density remain finite but higher derivatives of the expansion factor diverge) are excluded by both datasets.

For a broad range of parameter values, braneworld models display a remarkable property which we call cosmic mimicry. Cosmic mimicry is characterized by the fact that, at low redshifts, the Hubble parameter in the braneworld model is virtually indistinguishable from that in the LCDM cosmology. An important point to note is that the \Omega_m parameters in the braneworld model and in the LCDM cosmology can nevertheless be quite different. Thus, at high redshifts (early times), the braneworld asymptotically expands like a matter-dominated universe with the value of \Omega_m inferred from the observations of the local matter density. At low redshifts (late times), the braneworld model behaves almost exactly like the LCDM model but with a renormalized value of the cosmological density parameter \Omega_m^{LCDM}. The redshift which characterizes cosmic mimicry is related to the parameters in the higher-dimensional braneworld Lagrangian. Cosmic mimicry is a natural consequence of the scale-dependence of gravity in braneworld models. The change in the value of the cosmological density parameter is shown to be related to the spatial dependence of the effective gravitational constant in braneworld theory. A subclass of mimicry models lead to an older age of the universe and also predict a redshift of reionization which is lower than z_{reion} \simeq 17 in the LCDM cosmology. These models might therefore provide a background cosmology which is in better agreement both with the observed quasar abundance at z \gsim 4 and with the large optical depth to reionization measured by the Wilkinson Microwave Anisotropy Probe.

Properties of Braneworld models of dark energy are reviewed. Braneworld models admit the following interesting possibilities: (i) The effective equation of state can be w < -1 as well as w > -1. In the former case the expansion of the universe is well behaved at all times and the universe does not run into a future `Big Rip' singularity which is usually encountered by Phantom models. (ii) For a class of Braneworld models the acceleration of the universe can be a transient phenomenon. In this case the current acceleration of the universe is sandwiched between two matter dominated epochs. Such a braneworld does not have a horizon in contrast to LCDM and most Quintessence models. (iii) For a specific set of parameter values the universe can either originate from, or end its existence in a Quiescent singularity, at which the density, pressure and Hubble parameter remain finite, while the deceleration parameter and all invariants of the Riemann tensor diverge to infinity within a finite interval of cosmic time. (iv) Braneworld models of dark energy can loiter at high redshifts: $6 \lleq z \lleq 40$. The Hubble parameter decreasesduring the loitering epoch relative to its value in LCDM. As a result the age of the universe at loitering dramatically increases and this is expected to boost the formation of high redshift gravitationally bound systems such as $10^9 M_\odot$ black holes at $z \sim 6$ and lower-mass black holes and/or Population III stars at $z > 10$, whose existence could be problematic within the LCDM scenario. (v) Braneworld models with a time-like extra dimension bounce at early times thereby avoiding the initial `Big Bang singularity'. (vi) Both Inflation and Dark Energy can be successfully unified within a single scheme (Quintessential Inflation).

We show that loitering at high redshifts ($z \gsim 6$) can easily arise in braneworld models of dark energy which, in addition to being spatially flat, also accelerate at late times. Loitering is characterized by the fact that the Hubble parameter dips in value over a narrow redshift range which we shall refer to as the `loitering epoch'. During loitering, density perturbations are expected to grow rapidly. In addition, since the expansion of the universe slows down, its age near loitering dramatically increases. An early epoch of loitering is expected to boost the formation of high redshift gravitationally bound systems such as $10^9 M_\odot$ black holes at $z \sim 6$ and lower-mass black holes and/or Population III stars at $z > 10$, whose existence could be problematic within the LCDM scenario. Loitering models also help to reduce the redshift of reionization from its currently (high) value of $z_{\rm reion} \simeq 17$ in LCDM cosmology, thus alleviating a significant source of tension between observations of the high-redshift universe and theoretical model building. Currently a loitering universe accelerates with an effective equation of state $w < -1$ thus mimicking phantom dark energy. Unlike phantom, however, the late-time expansion of the universe in our model is singularity free, and a universe that loitered in the past will approach a LCDM model asymptotically in the distant future.

Modern redshift surveys such as the 2 degree field Galaxy Redshift Survey (2dFGRS) and the Sloan Digital Sky Survey (SDSS) reveal the fully 3 dimensional distribution of a million or so galaxies over a large cosmological volume. Visually galaxies appear to be distributed along sheet-like and/or filamentary superclusters. The CfA Great Wall, Southern Great Wall and the recently discovered SDSS Great Wall are very spectacular superclusters. Clearly theoretical predictions for galaxy clustering must be tested against these rich datasets. This can be achieved by means of the Minkowski Functionals (MFs). A MF-based approach provides an excellent description of superclusters and voids and allows one to quantify the properties of the cosmic web. In this review we give a summary of the progress made in this direction. After reviewing the status of observations and of numerical simulations, we comment upon the nature of bias which serves as a link between theoretical predictions and observations. We also summarise the methods developed for efficient numerical estimation of MFs for cosmological datasets and list several important results obtained using these methods. Specifically, we stress the discriminatory power of MFs and of the derived morphological statistics, the Shapefinders. Shapefinders are an excellent tool with which to study the shapes and sizes of superclusters and voids. We also discuss some of the important effects of scale-dependent bias which are brought out by a MF-based study of the mock catalogues of galaxies. Such effects, we note, should be accounted for before comparing theoretical models with observations.

Quintessential inflation describes a scenario in which both inflation and dark energy (quintessence) are described by the same scalar field. In conventional braneworld models of quintessential inflation gravitational particle production is used to reheat the universe. This reheating mechanism is very inefficient and results in an excessive production of gravity waves which violate nucleosynthesis constraints and invalidate the model. We describe a new method of realizing quintessential inflation on the brane in which inflation is followed by `instant preheating' (Felder, Kofman & Linde 1999). The larger reheating temperature in this model results in a smaller amplitude of relic gravity waves which is consistent with nucleosynthesis bounds. The relic gravity wave background has a `blue' spectrum at high frequencies and is a generic byproduct of successful quintessential inflation on the brane.

We reconstruct the equation of state $w(z)$ of dark energy (DE) using a recently released data set containing 172 type Ia supernovae without assuming the prior $w(z) \geq -1$ (in contrast to previous studies). We find that dark energy evolves rapidly and metamorphoses from dust-like behaviour at high $z$ ($w \simeq 0$ at $z \sim 1$) to a strongly negative equation of state at present ($w \lleq -1$ at $z \simeq 0$). Dark energy metamorphosis appears to be a robust phenomenon which manifests for a large variety of SNe data samples provided one does not invoke the weak energy prior $\rho + p \geq 0$. Invoking this prior considerably weakens the rate of growth of $w(z)$. These results demonstrate that dark energy with an evolving equation of state provides a compelling alternative to a cosmological constant if data are analysed in a prior-free manner and the weak energy condition is not imposed by hand.

In a recent paper (astro-ph/0311364) Alam et al argued that the SNe data of Tonry et al 2003 and Barris et al 2003 appear to favour DE which evolves in time, provided no other priors are invoked. (The effect of invoking priors such as the age of the Universe, the values of $H_0$ and $\Omega_m$ and CMB/LSS observations could modify this conclusion, as demonstrated in astro-ph/0403687 and other recent papers.) The approach adopted by Alam et al to reconstruct the properties of DE was severely (and, as we shall show below -- unfairly) criticized by J\"onsson et al in astro-ph/0404468. In this paper we re-examine the parametrisation used in astro-ph/0311364 and show that, contrary to the claims of J\"onsson et al, the results obtained from this reconstruction are robust and therefore representative of the true nature of dark energy.

We investigate the behaviour of dark energy using the recently released supernova data of Riess et al ~(2004) and a model independent parameterization for dark energy (DE). We find that, if no priors are imposed on $\Omega_{0m}$ and $h$, DE which evolves with time provides a better fit to the SNe data than $\Lambda$CDM. This is also true if we include results from the WMAP CMB data. From a joint analysis of SNe+CMB, the best-fit DE model has $w_0 < -1$ at the present epoch and the transition from deceleration to acceleration occurs at $z_T = 0.39 \pm 0.03$. However, DE evolution becomes weaker if the $\Lambda$CDM based CMB results $\Omega_{0m} = 0.27 \pm 0.04$, $h = 0.71 \pm 0.06$ are incorporated in the analysis. In this case, $z_T = 0.57 \pm 0.07$. Our results also show that the extent of DE evolution is sensitive to the manner in which the supernova data is sampled.

I briefly review our current understanding of dark matter and dark energy. The first part of this paper focusses on issues pertaining to dark matter including observational evidence for its existence, current constraints and the `abundance of substructure' and `cuspy core' issues which arise in CDM. I also briefly describe MOND. The second part of this review focusses on dark energy. In this part I discuss the significance of the cosmological constant problem which leads to a predicted value of the cosmological constant which is almost $10^{123}$ times larger than the observed value $\la/8\pi G \simeq 10^{-47}$GeV$^4$. Setting $\la$ to this small value ensures that the acceleration of the universe is a fairly recent phenomenon giving rise to the `cosmic coincidence' conundrum according to which we live during a special epoch when the density in matter and $\la$ are almost equal. Anthropic arguments are briefly discussed but more emphasis is placed upon dynamical dark energy models in which the equation of state is time dependent. These include Quintessence, Braneworld models, Chaplygin gas and Phantom energy. Model independent methods to determine the cosmic equation of state and the Statefinder diagnostic are also discussed. The Statefinder has the attractive property $\atridot/a H^3 = 1 $ for LCDM, which is helpful for differentiating between LCDM and rival dark energy models. The review ends with a brief discussion of the fate of the universe in dark energy models.

We explore a new class of braneworld models in which the scalar curvature of the (induced) brane metric contributes to the brane action. The scalar curvature term arises generically on account of one-loop effects induced by matter fields residing on the brane. Spatially flat braneworld models can enter into a regime of accelerated expansion at late times. This is true even if the brane tension and the bulk cosmological constant are tuned to satisfy the Randall--Sundrum constraint on the brane. Braneworld models admit a wider range of possibilities for dark energy than standard LCDM. In these models the luminosity distance can be both smaller and larger than the luminosity distance in LCDM. Whereas models with $d_L \leq d_L(\rm LCDM)$ imply $w = p/\rho \geq -1$ and have frequently been discussed in the literature, models with $d_L > d_L(\rm LCDM)$ have traditionally been ignored, perhaps because within the general-relativistic framework, the luminosity distance has this property {\em only if} the equation of state of matter is strongly negative ($w < -1$). Matter with $w < -1$ is beset with a host of undesirable properties, which makes this model of dark energy unattractive within the conventional framework. Braneworld models, on the other hand, have the capacity to endow dark energy with exciting new possibilities without suffering from the problems faced by models with $w < -1$. For a subclass of parameter values, braneworld dark energy and the acceleration of the universe are {\em transient} phenomena. In these models, the universe, after the current period of acceleration, re-enters the matter dominated regime so that the deceleration parameter $q(t) \to 0.5$ when $t \gg t_0$, where $t_0$ is the present epoch. Such models could help reconcile an accelerating universe with the requirements of string/M-theory.

We report here the first systematic study of the supercluster-void network in the LCDM concordance cosmology in which voids and superclusters are treated on an equal footing. Superclusters are defined as individual members of an over-dense excursion set and voids are defined as individual members of a complementary under-dense excursion set at the same density threshold. We determine the geometric, topological and morphological properties of the cosmic web at a large set of density levels by computing Minkowski functionals for every supercluster and void using SURFGEN. The properties of the largest (percolating) supercluster and the complementary void are found to be very different from properties of individual superclusters and voids. Individual superclusters totally occupy no more than about 5% of the total volume and contain no more than 20% of mass if the largest supercluster is excluded. Likewise, individual voids totally occupy no more than 14% of volume and contain no more than 4% of mass if the largest void is excluded. Although superclusters are more massive and voids are more voluminous the difference in maximum volumes is not greater than by an order of magnitude. The genus value of individual superclusters can be 5 while the genus of individual voids can reach 40, implying significant amount of substructure in superclusters and especially in voids. One of our main results is that large voids, as defined through the density field (read dark matter distribution) can be distinctly non-spherical.

The coming few years are likely to witness a dramatic increase in high quality Sn data as current surveys add more high redshift supernovae to their inventory and as newer and deeper supernova experiments become operational. Given the current variety in dark energy models and the expected improvement in observational data, an accurate and versatile diagnostic of dark energy is the need of the hour. This paper examines the Statefinder diagnostic in the light of the proposed SNAP satellite which is expected to observe about 2000 supernovae per year. We show that the Statefinder is versatile enough to differentiate between dark energy models as varied as the cosmological constant on the one hand, and quintessence, the Chaplygin gas and braneworld models, on the other. Using SNAP data, the Statefinder can distinguish a cosmological constant ($w=-1$) from quintessence models with $w \geq -0.9$ and Chaplygin gas models with $\kappa \leq 15$ at the $3\sigma$ level if the value of $\om$ is known exactly. The Statefinder gives reasonable results even when the value of $\om$ is known to only $\sim 20%$ accuracy. In this case, marginalizing over $\om$ and assuming a fiducial LCDM model allows us to rule out quintessence with $w \geq -0.85$ and the Chaplygin gas with $\kappa \leq 7$ (both at $3\sigma$). These constraints can be made even tighter if we use the Statefinders in conjunction with the deceleration parameter. The Statefinder is very sensitive to the total pressure exerted by all forms of matter and radiation in the universe. It can therefore differentiate between dark energy models at moderately high redshifts of $z \lleq 10$.

We explore the fate of the universe given the possibility that the density associated with `dark energy' may decay slowly with time. Decaying dark energy is modeled by a homogeneous scalar field which couples minimally to gravity and whose potential has {\em at least one} local quadratic maximum. Dark energy decays as the scalar field rolls down its potential, consequently the current acceleration epoch is a transient. We examine two models of decaying dark energy. In the first, the dark energy potential is modeled by an analytical form which is generic close to the potential maximum. The second potential is the cosine, which can become negative as the field evolves, ensuring that a spatially flat universe collapses in the future. We examine the feasibility of both models using observations of high redshift type Ia supernovae. A maximum likelihood analysis is used to find allowed regions in the $\lbrace m, \phi_0\rbrace$ plane ($m$ is the tachyon mass modulus and $\phi_0$ the initial scalar field value; $m\sim H_0$ and $\phi_0\sim M_P$ by order of magnitude). For the first model, the time for the potential to drop to half its maximum value is larger than $\sim 8$ Gyrs. In the case of the cosine potential, the time left until the universe collapses is always greater than $\sim 18$ Gyrs (both estimates are presented for $\om = 0.3$, $m/H_0 \sim 1$, $H_0 \simeq 70$ km/sec/Mpc, and at the 95.4% confidence level).

We present a new ansatz which can successfully be used to determine the morphological properties of the supercluster-void network. The ansatz is based on a surface modelling scheme SURFGEN, which generates a triangulated surface from a discrete data set representing (say) the distribution of galaxies in real (or redshift) space. Four Minkowski functionals -- surface area, volume, extrinsic curvature and genus -- describe the geometry and topology of the supercluster-void network. Ratio's of Minkowski functionals -- Shapefinders -- provide us with an excellent diagnostic of three dimensional shapes of clusters, superclusters and voids. Minkowski functionals and Shapefinders are determined for a triangulated iso-density surface using SURFGEN. SURFGEN is tested against both simply and multiply connected eikonal surfaces such as triaxial ellipsoids and tori. Remarkably, the first three Minkowski functionals are computed to better than 1% accuracy while the fourth (genus) is known exactly. SURFGEN also gives excellent results when applied to Gaussian random fields. Our results indicate that the surface modelling scheme SURFGEN is accurate and robust and can successfully be used to quantify the topology and morphology of the supercluster-void network in the universe. We apply SURFGEN to three cosmological models, $\L$CDM, $\T$CDM and SCDM and obtain interesting new results pertaining to the geometry, morphology and topology of large scale structure.

We introduce a new cosmological diagnostic pair $\lbrace r,s\rbrace$ called Statefinder. The Statefinder is dimensionless and, like the Hubble and deceleration parameters $H(z)$ and $q(z)$, is constructed from the scale factor of the Universe and its derivatives only. The parameter $r(z)$ forms the next step in the hierarchy of geometrical cosmological parameters used to study the Universe after $H$ and $q$, while the parameter $s(z)$ is a linear combination of $q$ and $r$ chosen in such a way that it does not depend upon the dark energy density $\Omega_X(z)$. The Statefinder pair $\lbrace r,s\rbrace$ is algebraically related to the the dark energy pressure-to-energy ratio $w=p/\epsilon$ and its time derivative, and sheds light on the nature of dark energy/quintessence. Its properties allow to usefully differentiate between different forms of dark energy with constant and variable $w$, including a cosmological constant ($w = -1$). The Statefinder pair can be determined to very good accuracy from a SNAP type experiment.

We study cosmological braneworld models with a single timelike extra dimension. Such models admit the intriguing possibility that a contracting braneworld experiences a natural bounce without ever reaching a singular state. This feature persists in the case of anisotropic braneworlds under some additional and not very restrictive assumptions. Generalizing our study to braneworld models containing an induced brane curvature term, we find that a FRW-type singularity is once again absent if the bulk extra dimension is timelike. In this case, the universe either has a non-singular origin or commences its expansion from a quasi-singular state during which both the Hubble parameter and the energy density and pressure remain finite while the curvature tensor diverges. The non-singular and quasi-singular behaviour which we have discovered differs both qualitatively and quantitatively from what is usually observed in braneworld models with spacelike extra dimensions and could have interesting cosmological implications.

Observations of high redshift supernovae indicate that the universe is accelerating. The hypothesis of `Dark energy' (cosmological constant, scalar field tracker potentials, braneworld models, etc.) has been advanced to explain this phenomenon. Sensitive tests of dark energy which can differentiate between rival models are clearly the need of the hour. The statefinder pair $\lbrace r,s\rbrace$ is a geometrical diagnostic which can play this role. $r$ & $s$ depend upon the third time derivative of the scale factor $\stackrel{...}{a}$ and provide the next logical step in the hierarchy of the cosmological parameter set after $H$ and $q$. The statefinder pair $\lbrace r,s\rbrace$ can be determined to high accuracy from a SNAP type experiment and allows us to successfully differentiate between dark energy models having constant as well as time-varying equations of state.

Braneworld models of dark energy are examined in the light of observations of high redshift type Ia supernovae. Braneworld models admit several novel and even exotic possibilities which include: (i) The effective equation of state of dark energy can be more negative than in LCDM: $w \leq -1$; (ii) A class of braneworld models can encounter a `quiescent' future singularity at which the energy density and the Hubble parameter remain well behaved, but higher derivatives of the expansion factor ($\stackrel{..}{a}$, $\stackrel{...}{a}$ etc.) diverge when the future singularity is reached; (iii) The current acceleration of the universe is a {\em transient feature} in a class of models in which dark energy `disappears' after a certain time, giving rise to a matter dominated universe in the future. Since horizons are absent in such a space-time, a braneworld model with {\em transient acceleration} might help reconcile current supernova-based observations of an accelerating universe with the demands of string/M-theory. A maximum likelihood analysis reveals that braneworld models satisfy the stringent demands imposed by high redshift supernovae and a large region in parameter space agrees marginally better with current observations than LCDM. For instance, models with $w < -1 (> -1)$ provide better agreement with data than LCDM for $\Omega_m \ggeq 0.3 (\lleq 0.25)$.

Traditionally, higher-dimensional cosmological models have sought to provide a description of the fundamental forces in terms of a unifying geometrical construction. In this essay we discuss how, in their present incarnation, higher-dimensional `braneworld' models might provide answers to a number of cosmological puzzles including the issue of dark energy and the nature of the big-bang singularity.

Higher-dimensional braneworld models which contain both bulk and brane curvature terms in the action admit cosmological singularities of rather unusual form and nature. These `quiescent' singularities, which can occur both during the contracting as well as the expanding phase, are characterised by the fact that while the matter density and Hubble parameter remain finite, all higher derivatives of the scale factor ($\stackrel{..}{a}$, $\stackrel{...}{a}$ etc.) diverge as the cosmological singularity is approached. The singularities are the result of the embedding of the (3+1)-dimensional brane in the bulk and can exist even in an empty homogeneous and isotropic (FRW) universe. The possibility that the present universe may expand into a singular state is discussed.

I briefly review the cosmological constant problem and the issue of dark energy (or quintessence). Within the framework of quantum field theory, the vacuum expectation value of the energy momentum tensor formally diverges as $k^4$. A cutoff at the Planck or electroweak scale leads to a cosmological constant which is, respectively, $10^{123}$ or $10^{55}$ times larger than the observed value, $\l/8\pi G \simeq 10^{-47}$ GeV$^4$. The absence of a fundamental symmetry which could set the value of $\l$ to either zero or a very small value leads to {\em the cosmological constant problem}. Most cosmological scenario's favour a large time-dependent $\l$-term in the past (in order to generate inflation at $z \gg 10^{10}$), and a small $\l$-term today, to account for the current acceleration of the universe at $z \lleq 1$. Constraints arising from cosmological nucleosynthesis, CMB and structure formation constrain $\l$ to be sub-dominant during most of the intermediate epoch $10^{10} < z < 1$. This leads to the {\em cosmic coincidence} conundrum which suggests that the acceleration of the universe is a recent phenomenon and that we live during a special epoch when the density in $\l$ and in matter are almost equal. Time varying models of dark energy can, to a certain extent, ameliorate the fine tuning problem (faced by $\l$), but do not resolve the puzzle of cosmic coincidence. I briefly review tracker models of dark energy, as well as more recent brane inspired ideas and the issue of horizons in an accelerating universe. Model independent methods which reconstruct the cosmic equation of state from supernova observations are also assessed. Finally, a new diagnostic of dark energy -- `Statefinder', is discussed.

We discuss a scenario in which extra dimensional effects allow a scalar field with a steep potential to play the dual role of the inflaton as well as dark energy (quintessence). The post-inflationary evolution of the universe in this scenario is generically characterised by a `kinetic regime' during which the kinetic energy of the scalar field greatly exceeds its potential energy resulting in a `stiff' equation of state for scalar field matter $P_\phi \simeq \rho_\phi$. The kinetic regime precedes the radiation dominated epoch and introduces an important new feature into the spectrum of relic gravity waves created quantum mechanically during inflation. The gravity wave spectrum increases with wavenumber for wavelengths shorter than the comoving horizon scale at the commencement of the radiative regime. This `blue tilt' is a generic feature of models with steep potentials and imposes strong constraints on a class of inflationary braneworld models. Prospects for detection of the gravity wave background by terrestrial and space-borne gravity wave observatories such as LIGO II and LISA are discussed.

We examine the behavior of an anisotropic brane-world in the presence of inflationary scalar fields. We show that, contrary to naive expectations, a large anisotropy does not adversely affect inflation. On the contrary, a large initial anisotropy introduces more damping into the scalar field equation of motion, resulting in greater inflation. The rapid decay of anisotropy in the brane-world significantly increases the class of initial conditions from which the observed universe could have originated. This generalizes a similar result in general relativity. A unique feature of Bianchi I brane-world cosmology appears to be that for scalar fields with a large kinetic term the initial expansion of the universe is quasi-isotropic. The universe grows more anisotropic during an intermediate transient regime until anisotropy finally disappears during inflationary expansion.

We examine the behaviour of a closed oscillating universe filled with a homogeneous scalar field and find that, contrary to naive expectations, such a universe expands to larger volumes during successive expansion epochs. This intriguing behaviour introduces an arrow of time in a system which is time-reversible. The increase in the maximum size of the universe is closely related to the work done on/by the scalar field during one complete oscillatory cycle which, in turn, is related to the asymmetry in the scalar field equation of state during expansion and collapse. Our analysis shows that scalar fields with polynomial potentials $V(\phi) = \lambda \phi^q$, $q > 1$ lead to a growing oscillation amplitude for the universe: the increase in amplitude between successive oscillations is more significant for smaller values of $q$. Such behaviour allows for the effective recycling of the universe. A recycled universe can be quite old and can resolve the flatness problem. These results have strong bearing on cosmological models in which the role of dark matter is played by a scalar field. They are also relevant for chaotic inflationary models of the early universe since they demonstrate that, even if the universe fails to inflate the first time around, it will eventually do so during future oscillatory cycles. Thus, the space of initial conditions favourable for chaotic inflation increases significantly.

I describe a new class of quintessence+CDM models in which late time scalar field oscillations can give rise to both quintessence and cold dark matter. Additionally, a versatile ansatz for the luminosity distance is used to reconstruct the quintessence equation of state in a model independent manner from observations of high redshift supernovae.

Observations of high-redshift supernovae indicate that the universe is accelerating. Here we present a {\em model-independent} method for estimating the form of the potential $V(\phi)$ of the scalar field driving this acceleration, and the associated equation of state $w_\phi$. Our method is based on a versatile analytical form for the luminosity distance $D_L$, optimized to fit observed distances to distant supernovae and differentiated to yield $V(\phi)$ and $w_\phi$. Our results favor $w_\phi\simeq -1$ at the present epoch, steadily increasing with redshift. A cosmological constant is consistent with our results.

Recent observations of Type 1a supernovae indicating an accelerating universe have once more drawn attention to the possible existence, at the present epoch, of a small positive Lambda-term (cosmological constant). In this paper we review both observational and theoretical aspects of a small cosmological Lambda-term. We discuss the current observational situation focusing on cosmological tests of Lambda including the age of the universe, high redshift supernovae, gravitational lensing, galaxy clustering and the cosmic microwave background. We also review the theoretical debate surrounding Lambda: the generation of Lambda in models with spontaneous symmetry breaking and through quantum vacuum polarization effects -- mechanisms which are known to give rise to alarge value of Lambda hence leading to the `cosmological constant problem'. More recent attempts to generate a small cosmological constant at the present epoch using either field theoretic techniques, or by modeling a dynamical Lambda-term by scalar fields are also extensively discussed. Anthropic arguments favouring a small cosmological constant are briefly reviewed. A comprehensive bibliography of recent work on Lambda is provided.

We propose a new class of quintessence models in which late times oscillations of a scalar field give rise to an effective equation of state which can be negative and hence drive the observed acceleration of the universe. Our ansatz provides a unified picture of quintessence and a new form of dark matter we call "Frustrated Cold Dark Matter" (FCDM). FCDM inhibits gravitational clustering on small scales and could provide a natural resolution to the core density problem for disc galaxy halos. Since the quintessence field rolls towards a small value, constraints on slow-roll quintessence models are safely circumvented in our model.

We apply Shapefinders, statistical measures of `shape' constructed from two dimensional partial Minkowski functionals, to study the degree of filamentarity in the Las Campanas Redshift Survey (LCRS). In two dimensions, three Minkowski functionals characterise the morphology of an object, they are: its perimeter (L), area (S), and genus. Out of L and S a single dimensionless Shapefinder Statistic, F can be constructed (0 <=F <=1). F acquires extreme values on a circle (F = 0) and a filament (F = 1). Using F, we quantify the extent of filamentarity in the LCRS by comparing our results with a Poisson distribution with similar geometrical properties and having the same selection function as the survey. Our results unambiguously demonstrate that the LCRS displays a high degree of filamentarity both in the Northern and Southern galactic sections a result that is in general agreement with the visual appearance of the catalogue. It is well known that gravitational clustering from Gaussian initial conditions gives rise to the development of non-Gaussianity reflected in the formation of a network-like filamentary structure on supercluster scales. Consequently the fact that the smoothed LCRS catalogue shows properties consistent with those of a Gaussian random field (Colley 1997) whereas the unsmoothed catalogue demonstrates the presence of filamentarity lends strong support to the conjecture that the large scale clustering of galaxies is driven by gravitational instability.

We apply Minkowski functionals and various derived measures to decipher the morphological properties of large-scale structure seen in simulations of gravitational evolution. Minkowski functionals of isodensity contours serve as tools to test global properties of the density field. Furthermore, we identify coherent objects at various threshold levels and calculate their partial Minkowski functionals. We propose a set of two derived dimensionless quantities, planarity and filamentarity, which reduce the morphological information in a simple and intuitive way. Several simulations of the gravitational evolution of initial power-law spectra provide a framework for systematic tests of our method.

We study a class of FRW spacetimes with a nonminimally coupled light massive scalar field. Values of the coupling parameter xi<0 enhance long range power in the vacuum expectation value of the energy-momentum tensor and fundamentally alter the nature of inflationary particle production: the energy density of created particles behaves like an effective cosmological constant, leading generically to Omega_m < 1 in clustered matter and providing a possible resolution of the ``Omega problem'' for low density cosmological models.

We analyse shapes of overdense regions (clusters and superclusters) in controlled N-body simulations of gravitational clustering with power law initial spectra P(k) \propto k^n, n = -3, -2, -1, 0. At values of the density just above the percolation transition the number of distinct (isolated) clusters peaks and we use this `natural threshold' to study the shapes and multiplicity function of clusters and superclusters. We find that the extent of both filamentarity and pancakeness increases as the simulation evolves, the former being appreciably larger than the latter at virtually all epochs and for all spectra considered by us. Our results also show that high density regions within very massive clusters/superclusters are likely to be noticeably filamentary or pancake/ribbon-like when compared to the less dense regions within these objects. We make a detailed study of two, moment-based `shape statistics' proposed, respectively, by Babul & Starkman (BS) and Luo & Vishniac (LV) and find that both LV and BS have certain built-in limitations: LV is biased towards oblate structures and tends to overemphasise this property; neither BS nor LV correctly describe the shape of strongly curved or topologically nontrivial objects. For instance, a thin filamentary torus and a ribbon are both described by BS and LV as being pancakes ! By contrast Shapefinders, a new shape diagnostic not constructed from density moments but from Minskowski functionals, does not suffer from these limitations and appears to faithfully reproduce the shapes of both simple and topologically complex objects.

We explore shapes of clusters and superclusters in the IRAS 1.2 Jy redshift survey with three reconstructions spanning the range \beta = 0.1, 0.5, 1.0, where \beta = \Omega^{0.6}/b, b is the bias factor and \Omega the present value of the dimensionless matter density. Comparing our results to Gaussian randomized reconstructions of the IRAS catalogue, we find structures having both planar and filamentary properties. For \beta = 0.5, 1.0 the largest structures in the survey have a distinct tendency to be filament-like} in general agreement with the results of N-body simulations.

We probe gravitational clustering in N-body simulations using geometrical descriptors sensitive to `connectedness': the genus curve, percolation and shape statistics. We find that both genus and percolation curves provide complementary probes of large scale structure topology and could be used to discriminate between models of structure formation and the analysis of observational data such as galaxy catalogs and MBR maps. An analysis of `shapes' in N-body simulations has shown that filaments are more pronounced than pancakes. To probe shapes of clusters and superclusters more rigorously we propose a new shape statistic which does not fit isodensity surfaces by ellipsoids (as done earlier). Our shape statistic is derived from fundamental properties of a compact body: its Minkowski functionals. The new shape statistic gives sensible results for topologically simple surfaces such as the ellipsoid, and for more complicated surfaces such as the torus. (Invited talk, to appear in: Proceedings of the IAU Symposium No. 183, Kyoto, Japan Aug. 1997, ed. K. Sato, Kluwer Academic Publ.)

We construct a set of shape-finders which determine shapes of compact surfaces (iso-density surfaces in galaxy surveys or N-body simulations) without fitting them to ellipsoidal configurations as done earlier. The new indicators arise from simple, geometrical considerations and are derived from fundamental properties of a surface such as its volume, surface area, integrated mean curvature and connectivity characterized by the Genus. These `Shapefinders' could be used to diagnose the presence of filaments, pancakes and ribbons in large scale structure. Their lower-dimensional generalization may be useful for the study of two-dimensional distributions such as temperature maps of the Cosmic Microwave Background.

We study topological properties of large scale structure in a set of scale free N-body simulations using the genus and percolation curves as topological characteristics. Our results show that as gravitational clustering advances, the density field shows an increasingly pronounced departure from Gaussianity reflected in the changing shape of the percolation curve and the changing amplitude and shape of the genus curve. Both genus and percolation curves differentiate between the connectedness of overdense and underdense regions if plotted against the density. When plotted against the filling factor the percolation curve alone retains this property. The genus curve shows a pronounced decrease in amplitude caused by phase correlations in the non-linear regime. Both genus and percolation curves provide complementary probes of large scale structure topology and can be used to discriminate between models of structure formation and the analysis of observational data such as galaxy catalogs and MBR maps.

The morphological nature of structures that form under gravitational instability has been of central interest to cosmology for over two decades. A remarkable feature of large scale structures in the Universe is that they occupy a relatively small fraction of the volume and yet show coherence on scales comparable to the survey size. With the aid of a useful synthesis of percolation analysis and shape statistics we explore the evolution of morphology of isolated density clumps in {\it real space} and that of the cluster distribution as a whole in scale-invariant cosmological models of gravitational instability. Our results, based on an exhaustive statistical analysis, indicate that at finite density thresholds one-dimensional filaments are more abundant than two-dimensional sheets (pancakes) at most epochs and for all spectra although the first singularities could be pancake-like. Both filamentarity and pancakeness of structures grow with time (in scale-free models this is equivalent to an increase in resolution) leading to the development of a long coherence length scale in simulations.

We study the behaviour of spherical Voids in Lagrangian perturbation theories L(n), of which the Zel'dovich approximation is the lowest order solution L(1). We find that at early times higher order L(n) give an increasingly accurate picture of Void expansion. However at late times particle trajectories in L(2) begin to turnaround and converge leading to the {\em contraction} of a Void, a sign of pathological behaviour. By contrast particle trajectories in L(3) are well behaved and this approximation gives results in excellent agreement with the exact top-hat solution as long as the Void is not too underdense. For very underdense Voids, L(3) evacuates the Void much too rapidly leading us to conclude that the Zel'dovich approximation L(1), remains the best approximation to apply to the late time study of Voids. The behavior of high order approximations in spherical voids is typical for asymptotic series and may be generic for Lagrangian perturbation theory.

We discuss various analytical approximation methods for following the evolution of cosmological density perturbations into the strong (i.e. nonlinear) clustering regime. These methods can be classified into five types: (i) simple extrapolations from linear theory, such as the high--peak model and the lognormal model; (ii) {\em dynamical} approximations, including the Zel'dovich approximation and its extensions; (iii) non--linear models based on purely geometric considerations, of which the main example is the Voronoi model; (iv) statistical solutions involving scaling arguments, such as the hierarchical closure {\em ansatz} for BBGKY, fractal models and the thermodynamic model of Saslaw; (v) numerical techniques based on particles and/or hydrodynamics. We compare the results of full dynamical evolution using particle codes and the various other approximation schemes. To put the models we discuss into perspective, we give a brief review of the observed properties of galaxy clustering and the statistical methods used to quantify it, such as correlation functions, power spectra, topology and spanning trees.

The potential is a constant to linear order in cosmological gravitational clustering. In this Letter we present results of testing the conjecture, proposed by Pauls and Melott (1995), that the effect of nonlinear evolution on the potential can be better described by smoothing it on the scale of nonlinearity. We show two-dimensional simulations consisting of an N-body code in which particles are accelerated not by their mutual attraction, but by the gradient of the initial potential smoothed on the current scale of nonlinearity. This approximation produces results considerably improved over using a constant potential to move particles, and it is generally better than any approximation we have tested, lending further support to the smoothing description of the evolved potential.

We study the development of gravitational instability in the strongly non-linear regime. For this purpose we use a number of statistical indicators such as filamentary statistics, spectrum of overdense/underdense regions and the void probability function, each of which probes a particular aspect of gravitational clustering. We use these statistical indicators to discriminate between different approximations to gravitational instability which we test against N-body simulations. The approximations which we test are, the truncated Zel'dovich approximation (TZ), the adhesion model (AM), and the frozen flow (FF) and linear potential (LP) approximations. Of these we find that FF and LP break down relatively early, soon after the non-linear length scale exceeds $R_*$ -- the mean distance between peaks of the gravitational potential. The reason for this break down is easy to understand, particles in FF are constrained to follow the streamlines of the initial velocity field. Shell crossing is absent in this case and structure gradually freezes as particles begin to collect near minima of the gravitational potential. In LP particles follow the lines of force of the primordial potential, oscillating about its minima at late times when the non-linear length scale $k_{\rm NL}^{-1}\simeq R_*$. Unlike FF and LP the adhesion model (and to some extent TZ) continues to give accurate results even at late times when $k_{\rm NL}^{-1} \ge R_*$. This is because both AM and TZ use

We apply the adhesion approximation to study the formation and evolution of voids in the Universe. Our simulations -- carried out using 128$^3$ particles in a cubical box with side 128 Mpc -- indicate that the void spectrum evolves with time and that the mean void size in the standard COBE-normalised Cold Dark Matter (hereafter CDM) model with $h_{50} = 1,$ scales approximately as $\bar D(z) = {\bar D_0\over \sqrt {1+z}},$ where $\bar D_0 \simeq 10.5$ Mpc. Interestingly, we find a strong correlation between the sizes of voids and the value of the primordial gravitational potential at void centers. This observation could in principle, pave the way towards reconstructing the form of the primordial potential from a knowledge of the observed void spectrum. Studying the void spectrum at different cosmological epochs, for spectra with a built in $k$-space cutoff we find that, the number of voids in a representative volume evolves with time. The mean number of voids first increases until a maximum value is reached (indicating that the formation of cellular structure is complete), and then begins to decrease as clumps and filaments merge leading to hierarchical clustering and the subsequent elimination of small voids. The cosmological epoch characterizing the completion of cellular structure occurs when the length scale going nonlinear approaches the mean distance between peaks of the gravitational potential. A central result of this paper is that

We compare different nonlinear approximations to gravitational clustering in the weakly nonlinear regime, using as a comparative statistic the evolution of non-Gaussianity which can be characterised by a set of numbers $S_p$ describing connected moments of the density field at the lowest order in $<\delta^2>$: $<\delta^n>_c \simeq S_n<\delta^2>^{n-1}$. Generalizing earlier work by Bernardeau (1992) we develop an ansatz to evaluate all $S_p$ in a given approximation by means of a generating function which can be shown to satisfy the equations of motion of a homogeneous spherical density enhancement in that approximation. On the basis of the values of we show that approximations formulated in Lagrangian space (such as the Zeldovich approximation and its extensions) are considerably more accurate than those formulated in Eulerian space such as the Frozen Flow and Linear Potential approximations. In particular we find that the $n$th order Lagrangian perturbation approximation correctly reproduces the first $n+1$ parameters $S_n$. We also evaluate the density probability distribution function for the different approximations in the quasi-linear regime and compare our results with an exact analytic treatment in the case of the Zeldovich approximation.

We asses the contribution to the observed large scale anisotropy of the cosmic microwave background radiation, arising from both gravity waves as well as adiabatic density perturbations, generated by a common inflationary mechanism in the early Universe. We find that for inflationary models predicting power law primordial spectra %$a(t) \propto t^p, p > 1$, $\dk^2 \propto k^n$, the relative contribution to the quadrupole anisotropy from gravity waves and scalar density perturbations, depends crucially upon $n$. For $n < 0.87$, gravity waves perturb the CMBR by a larger amount than density perturbations, whereas for $n > 0.87$ the reverse is true. %For nearly exponential inflation, on the other hand, the contribution to the %microwave anisotropy from scalar density perturbations always exceeds %that due to gravity waves. Normalising the amplitude of the density perturbation spectrum at large scales, using the observed value of the COBE quadrupole, we determine $(\delta M /M)_{16}$ -- the rms density contrast on scales $\sim 16 h_{50}^{-1}$ Mpc, for cosmological models with cold dark matter. We find that for $n < 0.75$, a large amount of biasing is required in order to reconcile theory with observations.