There are two strong clues about the quantum structure of spacetime and the gravitational dynamics, which are almost universally ignored in the conventional approaches to quantize gravity. The first clue is that null surfaces exhibit (observer dependent) thermal properties and possess a heat density. This suggests that spacetime, like matter, has microscopic degrees of freedom and its long wavelength limit should be described in thermodynamic language and not in a geometric language. Second clue is related to the existence of the cosmological constant. Its understanding from first principles will require the dynamical principles of the theory to be invariant under the shift $T^a_b \to T^a_b + (constant) \delta^a_b$. This puts strong constraints on the nature of gravitational dynamics and excludes metric tensor as a fundamental dynamical variable. In fact, these two clues are closely related to each other. When the dynamical principles are recast, respecting the symmetry $T^a_b \to T^a_b + (constant) \delta^a_b$, they automatically acquire a thermodynamic interpretation related to the first clue. The first part of this review provides a pedagogical introduction to thermal properties of the horizons, including some novel derivations. The second part describes some aspects of cosmological constant problem and the last part provides a perspective on gravity which takes into account these principles.

At the mesoscopic scales --- which interpolate between the macroscopic, classical, geometry and the microscopic, quantum, structure of spacetime --- one can identify the density of states of the geometry which arises from the existence of a zero-point length in the spacetime. This spacetime discreteness also associates an internal degree of freedom with each event, in the form of a fluctuating vector of constant norm. The equilibrium state, corresponding to the extremum of the total density of states of geometry plus matter, leads precisely to Einstein's equations. In fact, the field equation can now be reinterpreted as a zero-heat dissipation principle. The analysis of fluctuations around the equilibrium state (described by Einstein's equations), will provide new insights about quantum gravity.

A unique feature of gravity is its ability to control the information accessible to any specific observer. We quantify the notion of cosmic information ('CosmIn') for an eternal observer in the universe. Demanding the finiteness of CosmIn requires the universe to have a late-time accelerated expansion. Combining the introduction of CosmIn with generic features of the quantum structure of spacetime (e.g., the holographic principle), we present a holistic model for cosmology. We show that (i) the numerical value of the cosmological constant, as well as (ii) the amplitude of the primordial, scale invariant, perturbation spectrum can be determined in terms of a single free parameter, which specifies the energy scale at which the universe makes a transition from a pre-geometric phase to the classical phase. For a specific value of the parameter, we obtain the correct results for both (i) and (ii). This formalism also shows that the quantum gravitational information content of spacetime can be tested using precision cosmology.

Our knowledge about the universe has increased tremendously in the last three decades or so --- thanks to the progress in observations --- but our understanding has improved very little. There are several fundamental questions about our universe for which we have no answers within the current, operationally very successful, approach to cosmology. Worse still, we do not even know how to address some of these issues within the conventional approach to cosmology. This fact suggests that we are missing some important theoretical ingredients in the overall description of the cosmos. I will argue that these issues --- some of which are not fully appreciated or emphasized in the literature --- demand a paradigm shift: We should not think of the universe as described by a specific solution to the gravitational field equations; instead, it should be treated as a special physical system governed by a different mathematical description, rooted in the quantum description of spacetime. I will outline how this can possibly be done.

I describe an approach which connects classical gravity with the quantum microstructure of spacetime. The field equations arise from maximizing the density of states of matter plus geometry. The former is identified using the thermodynamics of null surfaces while the latter arises due to the existence of a zero-point length in the spacetime. The resulting field equations remain invariant when a constant is added to the matter Lagrangian, which is a symmetry of the matter sector. Therefore, the cosmological constant arises as an integration constant. A non-zero value $(\Lambda)$ of the cosmological constant renders the amount of cosmic information $(I_c)$ accessible to an eternal observer finite and hence is directly related to it. This relation allows us to determine the numerical value of $(\Lambda)$ from the quantum structure of spacetime.

I clarify the differences between various approaches in the literature which attempt to link gravity and thermodynamics. I then describe a new perspective based on the following features: (1) As in the case of any other matter field, the gravitational field equations should also remain unchanged if a constant is added to the Lagrangian; in other words, the field equations of gravity should remain invariant under the transformation $T^a_b \to T^a_b + \delta^a_b $(constant). (2) Each event of spacetime has a certain number ($f$) of microscopic degrees of freedom (`atoms of spacetime'). This quantity $f$ is proportional to the area measure of an equi-geodesic surface, centered at that event, when the geodesic distance tends to zero. The spacetime should have a zero-point length in order for $f$ to remain finite. (3) The dynamics is determined by extremizing the heat density at all events of the spacetime. The heat density is the sum of a part contributed by matter and a part contributed by the atoms of spacetime, with the latter being $L_P^{-4} f$. The implications of this approach are discussed.

General Relativity (GR) revolutionized the way we thought about gravity. After briefly describing the key successes of GR and its impact, I will discuss the major conceptual challenges it faces today. I conclude by outlining the prospective future directions of development, which hold the promise of deepening our understanding of the nature of gravity.

I introduce a covariant four-vector $\mathcal{G}^a[v]$, which can be interpreted as the momentum density attributed to the spacetime geometry by an observer with velocity $v^a$, and describe its properties: (a) Demanding that the total momentum of matter plus geometry is conserved for all observers, leads to the gravitational field equations. Thus, how matter curves spacetime is entirely determined by this principle of momentum conservation. (b) The $\mathcal{G}^a[v]$ can be related to the gravitational Lagrangian in a manner similar to the usual definition of Hamiltonian in, say, classical mechanics. (c) Geodesic observers in a spacetime will find that the conserved total momentum vanishes on-shell. (d) The on-shell, conserved, total energy in a region of space, as measured by the comoving observers, will be equal to the total heat energy of the boundary surface. (e) The off-shell gravitational energy in a region will be the sum of the ADM energy in the bulk plus the thermal energy of the boundary. These results suggest that $\mathcal{G}^a[v]$ can be a useful physical quantity to probe the gravitational theories.

I show that the gravitational dynamics in a bulk region of space can be connected to a thermodynamic description in the boundary of that region, thereby providing clear physical interpretations of several mathematical features of classical general relativity: (1) The Noether charge contained in a bulk region, associated with a specific time evolution vector field, has a direct thermodynamic interpretation as the gravitational heat content of the boundary surface. (2) This result, in turn, shows that all static spacetimes maintain holographic equipartition; in these spacetimes, the number of degrees of freedom in the boundary is equal to the number of degrees of freedom in the bulk. (3) In a general, evolving spacetime, the rate of change of gravitational momentum is related to the difference between the number of bulk and boundary degrees of freedom. It is this departure from the holographic equipartition which drives the time evolution of the spacetime. (4) When the equations of motion hold, the (naturally defined) total energy of the gravity plus matter within a bulk region, will be equal to the boundary heat content. (5) After motivating the need for an alternate description of gravity (if we have to solve the cosmological constant problem), I describe a thermodynamic variational principle based on null surfaces to achieve this goal. The concept of gravitational heat density of the null surfaces arises naturally from the Noether charge associated with the null congruence. The null surface variational principle, in fact, extremises the total heat content of the matter plus gravity system. Several variations on this theme and implications are described. [Abridged]

I show that in a general, evolving spacetime, the rate of change of gravitational momentum is related to the difference between the number of degrees of freedom in the bulk and the boundary of a region. This expresses the gravitational field equation in the thermodynamic language which is the natural description, if gravity is an emergent phenomenon. In all static spacetimes, the number of degrees of freedom in the boundary is equal to the number of degrees of freedom in the bulk; i.e., these spacetimes maintain holographic equipartition. It is the departure from holographic equipartition that drives the time evolution of the spacetime. This result, which is equivalent to Einstein's equations, provides an elegant, holographic, description of spacetime dynamics.

Observations indicate that our universe is characterized by a late-time accelerating phase, possibly driven by a cosmological constant $\Lambda$, with the dimensionless parameter $\Lambda L_P^2 \simeq 10^{-122}$, where $L_P = (G \hbar /c^3)^{1/2}$ is the Planck length. In this review, we describe how the emergent gravity paradigm provides a new insight and a possible solution to the cosmological constant problem. After reviewing the necessary background material, we identify the necessary and sufficient conditions for solving the cosmological constant problem. We show that these conditions are naturally satisfied in the emergent gravity paradigm in which (i) the field equations of gravity are invariant under the addition of a constant to the matter Lagrangian and (ii) the cosmological constant appears as an integration constant in the solution. The numerical value of this integration constant can be related to another dimensionless number (called CosMIn) that counts the number of modes inside a Hubble volume that cross the Hubble radius during the radiation and the matter dominated epochs of the universe. The emergent gravity paradigm suggests that CosMIn has the numerical value $4 \pi$, which, in turn, leads to the correct, observed value of the cosmological constant. Further, the emergent gravity paradigm provides an alternative perspective on cosmology and interprets the expansion of the universe itself as a quest towards holographic equipartition. We discuss the implications of this novel and alternate description of cosmology.

Lanczos-Lovelock models of gravity represent a natural and elegant generalization of Einstein's theory of gravity to higher dimensions. They are characterized by the fact that the field equations only contain up to second derivatives of the metric even though the action functional can be a quadratic or higher degree polynomial in the curvature tensor. Because these models share several key properties of Einstein's theory they serve as a useful set of candidate models for testing the emergent paradigm for gravity. This review highlights several geometrical and thermodynamical aspects of Lanczos-Lovelock models which have attracted recent attention.

We study the evolution of a quantum scalar field in a toy universe which has three stages of evolution, viz., (i) an early (inflationary) de Sitter phase (ii) radiation dominated phase and (iii) late-time (cosmological constant dominated) de Sitter phase. Using Schr\"odinger picture, the scalar field equations are solved separately for the three stages and matched at the transition points. The boundary conditions are chosen so that field modes in the early de Sitter evolves from Bunch-Davies vacuum state. We determine the (time-dependent) particle content of this quantum state for the entire evolution of the universe and describe the various features both numerically and analytically. We also describe the quantum to classical transition in terms of a {\it classicality parameter} which tracks the particle creation and its effect on phase space correlation of the quantum field.

If gravity is an emergent phenomenon, as suggested by several recent results, then the structure of the action principle for gravity should encode this fact. With this motivation we study several features of the Einstein-Hilbert action and establish direct connections with horizon thermodynamics. We begin by introducing the concept of holographically conjugate variables (HCVs) in terms of which the surface term in the action has a specific relationship with the bulk term. In addition to g_{ab} and its conjugate momentum \sqrt{-g} M^{cab}, this procedure allows us to (re)discover and motivate strongly the use of f^{ab}=\sqrt{-g}g^{ab} and its conjugate momentum N^c_{ab}. The gravitational action can then be interpreted as a momentum space action for these variables. We also show that many expressions in classical gravity simplify considerably in this approach. For example, the field equations can be written in a form analogous to Hamilton's equations for a suitable Hamiltonian if we use these variables. More importantly, the variation of the surface term, evaluated on any null surface which acts a local Rindler horizon can be given a direct thermodynamic interpretation. The term involving the variation of the dynamical variable leads to T\delta S while the term involving the variation of the conjugate momentum leads to S\delta T. We have found this correspondence only for the choice of variables (g_{ab}, \sqrt{-g} M^{cab}) or (f^{ab}, N^c_{ab}). We use this result to provide a direct thermodynamical interpretation of the boundary condition in the action principle, when it is formulated in a spacetime region bounded by the null surfaces. We analyse these features from several different perspectives and provide a detailed description, which offers insights about the nature of classical gravity and emergent paradigm.

Thermal properties of a static horizon, (like the entropy S, heat content TS etc.) can be obtained either from the surface term of the Einstein-Hilbert action or by evaluating the Noether charge, corresponding to the diffeomorphisms generated by the timelike Killing vector field. We show that, for a wide class of geometries, the same results can be obtained using the vector field which produces an infinitesimal coordinate transformation between two physically relevant reference frames, viz. the freely falling frame near the horizon and the static, accelerated, frame. In particular, the infinitesimal coordinate transformation from inertial coordinates to uniformly accelerated frame can be used to obtain the heat content and entropy of the Rindler horizon. This result offers insight into understanding the observer dependent degrees of freedom which contribute to the entropy of null surfaces.

The current acceleration of the universe can be modeled in terms of a cosmological constant. We show that the extremely small value of \Lambda L_P^2 ~ 3.4 x 10^{-122}, the holy grail of theoretical physics, can be understood in terms of a new, dimensionless, conserved number CosMIn (N), which counts the number of modes crossing the Hubble radius during the three phases of evolution of the universe. Theoretical considerations suggest that N ~ 4\pi. This single postulate leads us to the correct, observed numerical value of the cosmological constant! This approach also provides a unified picture of cosmic evolution relating the early inflationary phase to the late-time accelerating phase.

We study quasi-stationary physical process for black holes within the context of Lanczos-Lovelock gravity. We show that the Wald entropy of stationary black holes in Lanczos-Lovelock gravity monotonically increases for quasi-stationary physical processes in which the horizon is perturbed by the accretion of positive energy matter and the black hole ultimately settles down to a stationary state. This result reinforces the physical interpretation of Wald entropy for Lanczos-Lovelock models and takes a step towards proving the analogue of the black hole area increase-theorem in a wider class of gravitational theories.

An interesting feature of the Davies-Unruh effect is that an uniformly accelerated observer sees an isotropic thermal spectrum of particles even though there is a preferred direction in this context, determined by the direction of the acceleration g. We investigate the thermal fluctuations in the Unruh bath by studying the Brownian motion of particles in the bath, especially as regards to isotropy. We find that the thermal fluctuations are anisotropic and induce different frictional drag forces on the Brownian particle depending on whether it has a drift velocity along the direction of acceleration g or in a direction transverse to it. Using the fluctuation-dissipation theorem, we argue that this anisotropy arises due to quantum correlations in the fluctuations at large correlation time scales.

Observations indicate that the evolution of our universe can be divided into three epochs consisting of early time inflation, radiation (and matter) domination and the late time acceleration. One can associate with each of these epochs a number N which is the phase space volume of the modes which cross the Hubble radius during the corresponding epoch. This number turns out to be (approximately) the same for the cosmologically relevant ranges of the three epochs. When the initial de Sitter space is characterized by the Planck length, the natural value for N is 4\pi. This allows us to determine the cosmological constant which drives the late time acceleration, to be \Lambda L_P^2 = 3 \exp(-24\pi^2 \mu) where \mu\ is a number of order unity. This expression leads to the observed value of cosmological constant for \mu ~ 1.19. The implications are discussed.

I describe several conceptual aspects of a particular paradigm which treats the field equations of gravity as emergent. These aspects are related to the features of classical gravitational theories which defy explanation within the conventional perspective. The alternative interpretation throws light on these features and could provide better insights into possible description of quantum structure of spacetime. This review complements the discussion in arXiv:1207.0505, which describes space itself as emergent in the cosmological context.

There is sufficient amount of internal evidence in the nature of gravitational theories to indicate that gravity is an emergent phenomenon like, e.g, elasticity. Such an emergent nature is most apparent in the structure of gravitational dynamics. It is, however, possible to go beyond the field equations and study the space itself as emergent in a well-defined manner in (and possibly only in) the context of cosmology. In the first part of this review, I describe various pieces of evidence which show that gravitational field equations are emergent. In the second part, I describe a novel way of studying cosmology in which I interpret the expansion of the universe as equivalent to the emergence of space itself. In such an approach, the dynamics evolves towards a state of holographic equipartition, characterized by the equality of number of bulk and surface degrees of freedom in a region bounded by the Hubble radius. This principle correctly reproduces the standard evolution of a Friedmann universe. Further, (a) it demands the existence of an early inflationary phase as well as late time acceleration for its successful implementation and (b) allows us to link the value of late time cosmological constant to the e-folding factor during inflation.

One possible interpretation of the holographic principle is the equality of the number of degrees of freedom in a bulk region of space and the number of degrees of freedom on the boundary surface. It is known that such an equality is maintained on equipotential surfaces in any static spacetime in the form of an equipartition law N_{bulk}= N_{sur}. In the cosmological context, the de Sitter universe obeys the same holographic equipartition. I argue that the difference between the surface degrees of freedom and the bulk degrees of freedom in a region of space (which has already emerged) drives the accelerated expansion of the universe through a simple equation dV/dt = (N_{sur} - N_{bulk}) where V is the Hubble volume in Planck units and t is the cosmic time in Planck units. This equation reproduces the standard evolution of the universe. This approach provides a novel paradigm to study the emergence of space and cosmology and has far reaching implications.

A class of theories of gravity based on a Lagrangian which depends on the curvature and metric - but not on the derivatives of the curvature tensor - is of interest in several contexts including in the development of the paradigm that treats gravity as an emergent phenomenon. This class of models contains, as an important subset, all Lanczos-Lovelock models of gravity. I derive several identities and properties which are useful in the study of these models and clarify some of the issues that seem to have received insufficient attention in the past literature.

I present the theoretical evidence which suggests that gravity is an emergent phenomenon like gas dynamics or elasticity with the gravitational field equations having the same status as, say, the equations of fluid dynamics/elasticity. This paradigm views a wide class of gravitational theories - including Einstein's theory - as describing the thermodynamic limit of the statistical mechanics of "atoms of spacetime". The evidence for this paradigm is hidden in several classical features of the gravitational theories and depends on just one quantum mechanical input, viz. the existence of Davies-Unruh temperature of horizons. I discuss several conceptual ingredients of this approach.

I show that the principle of equipartition, applied to area elements of a surface which are in equilibrium at the local Davies-Unruh temperature, allows one to determine the surface number density of the microscopic spacetime degrees of freedom in any diffeomorphism invariant theory of gravity. The entropy associated with these degrees of freedom matches with the Wald entropy for the theory. This result also allows one to attribute an entropy density to the spacetime in a natural manner. The field equations of the theory can then be obtained by extremising this entropy. Moreover, when the microscopic degrees of freedom are in local thermal equilibrium, the spacetime entropy of a bulk region resides on its boundary.

We study the observational constraints of CMB temperature and polarization anisotropies on models of dark energy, with special focus on models with variation in properties of dark energy with time. We demonstrate that the key constraint from CMB observations arises from the location of acoustic peaks. An additional constraint arises from the limits on Omega_nr from the relative amplitudes of acoustic peaks. Further, we show that the distance to the last scattering surface is not how the CMB observations constrain the combination of parameters for models of dark energy. We also use constraints from Supernova observations and show that unlike the Gold and Silver samples, the SNLS sample prefers a region of parameter space that has a significant overlap with the region preferred by the CMB observations. This is a verification of a conjecture made by us in an earlier work. We discuss combined constraints from WMAP5 and SNLS observations. We find that models with w ~ -1 are preferred for models with a constant equation of state parameters. In case of models with a time varying dark energy, we show that constraints on evolution of dark energy density are almost independent of the type of variation assumed for the equation of state parameter. This makes it easy to get approximate constraints from CMB observations on arbitrary models of dark energy. Constraints on models with a time varying dark energy are predominantly due to CMB observations, with Supernova constraints playing only a marginal role.

The purpose of the paper is five-fold: (a) Argue that the question in the title can be presented in a meaningful manner and that it requires an answer. (b) Discuss the conventional answers and explain why they are unsatisfactory. (c) Suggest that a key ingredient in the answer could be the instability arising due to the `wrong' sign in the Hilbert action for the kinetic energy term corresponding to expansion factor. (d) Describe how this idea connects up with another peculiar feature of our universe, viz. it spontaneously became more and more classical in the course of evolution. (e) Provide a speculative but plausible scenario, based on the thermodynamic perspective of gravity, in which one has the hope for relating the thermodynamic and cosmological arrows of time.

The fact that one can associate thermodynamic properties with horizons brings together principles of quantum theory, gravitation and thermodynamics and possibly offers a window to the nature of quantum geometry. This review discusses certain aspects of this topic concentrating on new insights gained from some recent work. After a brief introduction of the overall perspective, Sections 2 and 3 provide the pedagogical background on the geometrical features of bifurcation horizons, path integral derivation of horizon temperature, black hole evaporation, structure of Lanczos-Lovelock models, the concept of Noether charge and its relation to horizon entropy. Section 4 discusses several conceptual issues introduced by the existence of temperature and entropy of the horizons. In Section 5 we take up the connection between horizon thermodynamics and gravitational dynamics and describe several peculiar features which have no simple interpretation in the conventional approach. The next two sections describe the recent progress achieved in an alternative perspective of gravity. In Section 6 we provide a thermodynamic interpretation of the field equations of gravity in any diffeomorphism invariant theory and in Section 7 we obtain the field equations of gravity from an entropy maximization principle. The last section provides a summary.

I describe the conceptual and mathematical basis of an approach which describes gravity as an emergent phenomenon. Combining principle of equivalence and principle of general covariance with known properties of local Rindler horizons, perceived by observers accelerated with respect to local inertial frames, one can prove that the field equations describing gravity in any diffeomorphism invariant theory can be given a thermodynamic re-interpretation. This fact, in turn, leads us to the possibility of deriving the field equations of gravity by maximising a suitably defined entropy functional, without using the metric tensor as a dynamical variable. The approach synthesises concepts from quantum theory, thermodynamics and gravity leading to a fresh perspective on the nature of gravity. The description is presented here in the form of a dialogue, thereby addressing several frequently-asked-questions.

It is possible to provide a thermodynamic interpretation for the field equations in any diffeomorphism invariant theory of gravity. This insight, in turn, leads us to the possibility of deriving the gravitational field equations from another variational principle without using the metric as a dynamical variable. I review this approach and discuss its implications.

We provide a detailed analysis of quantum field theory around a collapsing shell and discuss several conceptual issues related to the emission of radiation flux and formation of black holes. Explicit calculations are performed using a model for a collapsing shell which turns out to be analytically solvable. We use the insights gained in this model to draw reliable conclusions regarding more realistic models. We first show that any shell of mass $M$ which collapses to a radius close to $r=2M$ will emit approximately thermal radiation for a period of time. In particular, a shell which collapses from some initial radius to a final radius $2M(1-\epsilon^2)^{-1}$ (where $\epsilon \ll 1$) without forming a black hole, will emit thermal radiation during the period $M\lesssim t \lesssim M\ln (1/\epsilon^2)$. Later on ($t\gg M \ln(1/\epsilon^2)$), the flux from such a shell will decay to zero exponentially. We next study the effect of backreaction computed using the vacuum expectation value of the stress tensor on the collapse. We find that, in any realistic collapse scenario, the backreaction effects do \emph{not} prevent the formation of the event horizon. The time at which the event horizon is formed is, of course, delayed due to the radiated flux -- which decreases the mass of the shell -- but this effect is not sufficient to prevent horizon formation. We also clarify several conceptual issues and provide pedagogical details of the calculations in the Appendices to the paper.

The epoch of reionization and formation of first stars are inter-linked topics that are of considerable interest. We use a simplified approach for studying formation of stars in collapsed haloes, and the resulting ionization of the inter-galactic medium (IGM). We consider a set of LCDM models allowed by observations of CMB temperature and polarization anisotropies for this study. We constrain parameters related to star formation with the help of observations. We constrain subsets of these parameters independently by using the observed metallicity of the inter-galactic medium at z ~ 5 and the requirement that the Thomson scattering optical depth due to an ionized IGM as determined for the model from CMB observations be reproduced. We consider a range of initial metalicities for star forming gas, and some variations of the initial mass function of stars. We find that a "normal" initial mass function (IMF) may satisfy these two constraints with a raised efficiency of star formation as compared to that seen in the local universe. Observations require a significant fraction of metals to escape from haloes to the IGM. We can also place constraints on the ratio of escape fraction for metals and ionizing photons, we find that this ratio is of order unity for most models. This highlights the importance of using the constraints arising from enrichment of the inter-galactic medium. (Abridged)

I argue that the field equations of any theory of gravity which is diffeomorphism invariant must be expressible as a thermodynamic identity, TdS=dE around any event in the spacetime. This fact can be demonstrated explicitly (and rather easily) if: (a) one accepts the Noether current of the theory as providing the definition for local entropy density and (b) one is allowed to introduce the local notions of a Rindler frame, acceleration horizon and a Killing vector (related to translation in Rindler time) around any event. It is conceptually incorrect - in general - to invert this argument and obtain the field equations of the theory from the thermodynamic identity. I discuss under what conditions this may be possible. Several subtleties related to these arguments are described.

I review several issues related to statistical description of gravitating systems in both static and expanding backgrounds. After briefly reviewing the results for the static background, I concentrate on gravitational clustering of collisionless particles in an expanding universe. In particular, I describe (a) how the non linear mode-mode coupling transfers power from one scale to another in the Fourier space if the initial power spectrum is sharply peaked at a given scale and (b) the asymptotic characteristics of gravitational clustering which are independent of the initial conditions. Numerical simulations as well as analytic work shows that power transfer leads to a universal power spectrum at late times, somewhat reminiscent of the existence of Kolmogorov spectrum in fluid turbulence.

The cosmological constant is the most economical candidate for dark energy. No other approach really alleviates the difficulties faced by the cosmological constant because, in all other attempts to model the dark energy, one still has to explain why the bulk cosmological constant (treated as a low-energy parameter in the action principle) is zero. I argue that until the theory is made invariant under the shifting of the Lagrangian by a constant, one cannot obtain a satisfactory solution to the cosmological constant problem. This is impossible in any generally covariant theory with the conventional low-energy matter action, if the metric is varied in the action to obtain the field equations. I review an alternative perspective in which gravity arises as an emergent, long wavelength phenomenon and can be described in terms of an effective theory using an action associated with null vectors in the spacetime. This action is explicitly invariant under the shift of the energy momentum tensor $T_{ab}\to T_{ab}+\Lambda g_{ab}$ and any bulk cosmological constant can be gauged away. Such an approach seems to be necessary for addressing the cosmological constant problem and can easily explain why its bulk value is zero. I describe some possibilities for obtaining its observed value from quantum gravitational fluctuations.

In Einstein's gravity, the entropy of horizons is proportional to their area. Several arguments given in the literature suggest that, in this context, both area and entropy should be quantized with an equally spaced spectrum for large quantum numbers. But in more general theories (like, for e.g, in the black hole solutions of Gauss-Bonnet or Lanczos-Lovelock gravity) the horizon entropy is \emph{not} proportional to area and the question arises as to which of the two (if at all) will have this property. We give a general argument that in all Lanczos-Lovelock theories of gravity, it is the \emph{entropy} that has equally spaced spectrum. In the case of Gauss-Bonnet gravity, we use the asymptotic form of quasi normal mode frequencies to explicitly demonstrate this result. Hence, the concept of a quantum of area in Einstein Hilbert (EH) gravity needs to be replaced by a concept of \emph{quantum of entropy} in a more general context.

We adopt the general formalism, which was developed in Paper I (arXiv:0708.1233) to analyze the evolution of a quantized time-dependent oscillator, to address several questions in the context of quantum field theory in time dependent external backgrounds. In particular, we study the question of emergence of classicality in terms of the phase space evolution and its relation to particle production, and clarify some conceptual issues. We consider a quantized scalar field evolving in a constant electric field and in FRW spacetimes which illustrate the two extreme cases of late time adiabatic and highly non-adiabatic evolution. Using the time-dependent generalizations of various quantities like particle number density, effective Lagrangian etc. introduced in Paper I, we contrast the evolution in these two limits bringing out key differences between the Schwinger effect and evolution in the de Sitter background. Further, our examples suggest that the notion of classicality is multifaceted and any one single criterion may not have universal applicability. For example, the peaking of the phase space Wigner distribution on the classical trajectory \emph{alone} does not imply transition to classical behavior. An analysis of the behavior of the \emph{classicality parameter}, which was introduced in Paper I, leads to the conclusion that strong particle production is necessary for the quantum state to become highly correlated in phase space at late times.

The quantum theory of a harmonic oscillator with a time dependent frequency arises in several important physical problems, especially in the study of quantum field theory in an external background. While the mathematics of this system is straightforward, several conceptual issues arise in such a study. We present a general formalism to address some of the conceptual issues like the emergence of classicality, definition of particle content, back reaction etc. In particular, we parametrize the wave function in terms of a complex number (which we call excitation parameter) and express all physically relevant quantities in terms it. Many of the notions -- like those of particle number density, effective Lagrangian etc., which are usually defined using asymptotic in-out states -- are generalized as time-dependent concepts and we show that these generalized definitions lead to useful and reasonable results. Having developed the general formalism we apply it to several examples. Exact analytic expressions are found for a particular toy model and approximate analytic solutions are obtained in the extreme cases of adiabatic and highly non-adiabatic evolution. We then work out the exact results numerically for a variety of models and compare them with the analytic results and approximations. The formalism is useful in addressing the question of emergence of classicality of the quantum state, its relation to particle production and to clarify several conceptual issues related to this. In Paper II (arXiv:0708.1237), which is a sequel to this, the formalism will be applied to analyze the corresponding issues in the context of quantum field theory in background cosmological models and electric fields.

The motion of a black hole about the centre of gravity of its host galaxy induces a strong response from the surrounding stellar population. We treat the case of a harmonic potential analytically and show that half of the stars on circular orbits in that potential shift to an orbit of lower energy, while the other half receive a positive boost and recede to a larger radius. The black hole itself remains on an orbit of fixed amplitude and merely acts as a catalyst for the evolution of the stellar energy distribution function f(E). We show that this effect is operative out to a radius of approx 3 to 4 times the hole's influence radius, R_bh. We use numerical integration to explore more fully the response of a stellar distribution to black hole motion. We consider orbits in a logarithmic potential and compare the response of stars on circular orbits, to the situation of a `warm' and `hot' (isotropic) stellar velocity field. While features seen in density maps are now wiped out, the kinematic signature of black hole motion still imprints the stellar line-of-sight mean velocity to a magnitude ~18% the local root mean-square velocity dispersion sigma.

I describe several broad features of a programme to understand gravity as an emergent, long wavelength, phenomenon (like elasticity) and discuss one concrete framework for realizing this paradigm in the backdrop of several recent results.

I review the problem of dark energy focusing on the cosmological constant as the candidate and discuss its implications for the nature of gravity. Part 1 briefly overviews the currently popular `concordance cosmology' and summarises the evidence for dark energy. It also provides the observational and theoretical arguments in favour of the cosmological constant as the candidate and emphasises why no other approach really solves the conceptual problems usually attributed to the cosmological constant. Part 2 describes some of the approaches to understand the nature of the cosmological constant and attempts to extract the key ingredients which must be present in any viable solution. I argue that (i)the cosmological constant problem cannot be satisfactorily solved until gravitational action is made invariant under the shift of the matter lagrangian by a constant and (ii) this cannot happen if the metric is the dynamical variable. Hence the cosmological constant problem essentially has to do with our (mis)understanding of the nature of gravity. Part 3 discusses an alternative perspective on gravity in which the action is explicitly invariant under the above transformation. Extremizing this action leads to an equation determining the background geometry which gives Einstein's theory at the lowest order with Lanczos-Lovelock type corrections. (Condensed abstract).

The null surfaces of a spacetime act as one-way membranes and can block information for a corresponding family of observers (time-like curves). Since lack of information can be related to entropy, this suggests the possibility of assigning an entropy to the null surfaces of a spacetime. We motivate and introduce such an entropy functional for any vector field in terms of a fourth-rank divergence free tensor P_{ab}^{cd} with the symmetries of the curvature tensor. Extremising this entropy for all the null surfaces then leads to equations for the background metric of the spacetime. When P_{ab}^{cd} is constructed from the metric alone, these equations are identical to Einstein's equations with an undetermined cosmological constant (which arises as an integration constant). More generally, if P_{ab}^{cd} is allowed to depend on both metric and curvature in a polynomial form, one recovers the Lanczos-Lovelock gravity. In all these cases: (a) We only need to extremise the entropy associated with the null surfaces; the metric is not a dynamical variable in this approach. (b) The extremal value of the entropy agrees with standard results, when evaluated on-shell for a solution admitting a horizon. The role of full quantum theory of gravity will be to provide the specific form of P_{ab}^{cd} which should be used in the entropy functional. With such an interpretation, it seems reasonable to interpret the Lanczos-Lovelock type terms as quantum corrections to classical gravity.

There is an intriguing analogy between the gravitational dynamics of the horizons and thermodynamics. In case of general relativity, as well as for a wider class of Lanczos-Lovelock theories of gravity, it is possible to interpret the field equations near any spherically symmetric horizon as a thermodynamic identity TdS = dE + PdV. We study this approach further and generalize the results to two more generic cases within the context of general relativity: (i) stationary axis-symmetric horizons and (ii) time dependent evolving horizons. In both the cases, the near horizon structure of Einstein equations can be expressed as a thermodynamic identity under the virtual displacement of the horizon. This result demonstrates the fact that the thermodynamic interpretation of gravitational dynamics is not restricted to spherically symmetric or static horizons but is quite generic in nature and indicates a deeper connection between gravity and thermodynamics.

The equations of motion for matter fields are invariant under the shift of the matter lagrangian by a constant. Such a shift changes the energy momentum tensor of matter by T^a_b --> T^a_b +\rho \delta^a_b. In the conventional approach, gravity breaks this symmetry and the gravitational field equations are not invariant under such a shift of the energy momentum tensor. I argue that until this symmetry is restored, one cannot obtain a satisfactory solution to the cosmological constant problem. I describe an alternative perspective to gravity in which the gravitational field equations are [G_{ab} -\kappa T_{ab}] n^an^b =0 for all null vectors n^a. This is obviously invariant under the change T^a_b --> T^a_b +\rho \delta^a_b and restores the symmetry under shifting the matter lagrangian by a constant. These equations are equivalent to G_{ab} = \kappa T_{ab} + Cg_{ab} where C is now an integration constant so that the role of the cosmological constant is very different in this approach. The cosmological constant now arises as an integration constant, somewhat like the mass M in the Schwarzschild metric, the value of which can be chosen depending on the physical context. These equations can be obtained from a variational principle which uses the null surfaces of spacetime as local Rindler horizons and can be given a thermodynamic interpretation. This approach turns out to be quite general and can encompass even the higher order corrections to Einstein's gravity and suggests a principle to determine the form of these corrections in a systematic manner.

Einstein-Hilbert (EH) action can be separated into a bulk and a surface term, with a specific ("holographic") relationship between the two, so that either can be used to extract information about the other. The surface term can also be interpreted as the entropy of the horizon in a wide class of spacetimes. Since EH action is likely to just the first term in the derivative expansion of an effective theory, it is interesting to ask whether these features continue to hold for more general gravitational actions. We provide a comprehensive analysis of lagrangians of the form L=Q_a^{bcd}R^a_{bcd}, in which Q_a^{bcd} is a tensor with the symmetries of the curvature tensor, made from metric and curvature tensor and satisfies the condition \nabla_cQ^{abcd}=0, and show that they share these features. The Lanczos-Lovelock lagrangians are a subset of these in which Q^{abcd} is a homogeneous function of the curvature tensor. They are all holographic, in a specific sense of the term, and -- in all these cases -- the surface term can be interpreted as the horizon entropy. The thermodynamics route to gravity, in which the field equations are interpreted as TdS=dE+pdV, seems to have greater degree of validity than the field equations of Einstein gravity itself. The results suggest that the holographic feature of EH action could also serve as a new symmetry principle in constraining the semiclassical corrections to Einstein gravity. The implications are discussed.

Spacetimes with horizons show a resemblance to thermodynamic systems and one can associate the notions of temperature and entropy with them. In the case of Einstein-Hilbert gravity, it is possible to interpret Einstein's equations as the thermodynamic identity TdS = dE + PdV for a spherically symmetric spacetime and thus provide a thermodynamic route to understand the dynamics of gravity. We study this approach further and show that the field equations for Lanczos-Lovelock action in a spherically symmetric spacetime can also be expressed as TdS = dE + PdV with S and E being given by expressions previously derived in the literature by other approaches. The Lanczos-Lovelock Lagrangians are of the form L=Q_a^{bcd}R^a_{bcd} with \nabla_b Q^{abcd}=0. In such models, the expansion of Q^{abcd} in terms of the derivatives of the metric tensor determines the structure of the theory and higher order terms can be interpreted quantum corrections to Einstein gravity. Our result indicates a deep connection between the thermodynamics of horizons and the allowed quantum corrections to standard Einstein gravity, and shows that the relation TdS = dE + PdV has a greater domain of validity that Einstein's field equations.

Nearly seventy per cent of the energy density in the universe is unclustered and exerts negative pressure. This conclusion -- now supported by numerous observations -- poses the greatest challenge for theoretical physics today. I discuss this issue with special emphasis on the cosmological constant as the possible choice for the dark energy. Several curious features of a universe with a cosmological constant are described and some possible approaches to understand the nature of the cosmological constant are reviewed. In particular, I show how some of the recent ideas, related to a thermodynamic route to gravity, allow us to: (i) create a paradigm in which the bulk value of cosmological constant is irrelevant and (ii) obtain the correct, observed, value for the cosmological constant from vacuum fluctuations in a region confined by the deSitter horizon.

The Einstein-Hilbert action has a bulk term and a surface term (which arises from integrating a four divergence). I show that one can obtain Einstein's equations from the surface term alone. This leads to: (i) a novel, completely self contained, perspective on gravity and (ii) a concrete mathematical framework in which the description of spacetime dynamics by Einstein's equations is similar to the description of a continuum solid in the thermodynamic limit.

A general paradigm for describing classical (and semiclassical) gravity is presented. This approach brings to the centre-stage a holographic relationship between the bulk and surface terms in a general class of action functionals and provides a deeper insight into several aspects of classical gravity which have no explanation in the conventional approach. After highlighting a series of unresolved issues in the conventional approach to gravity, I show that (i) principle of equivalence, (ii) general covariance and (iii)a reasonable condition on the variation of the action functional, suggest a generic Lagrangian for semiclassical gravity of the form $L=Q_a^{bcd}R^a_{bcd}$ with $\nabla_b Q_a^{bcd}=0$. The expansion of $Q_a^{bcd}$ in terms of the derivatives of the metric tensor determines the structure of the theory uniquely. The zeroth order term gives the Einstein-Hilbert action and the first order correction is given by the Gauss-Bonnet action. Any such Lagrangian can be decomposed into a surface and bulk terms which are related holographically. The equations of motion can be obtained purely from a surface term in the gravity sector. Hence the field equations are invariant under the transformation $T_{ab} \to T_{ab} + \lambda g_{ab}$ and gravity does not respond to the changes in the bulk vacuum energy density. The cosmological constant arises as an integration constant in this approach. The implications are discussed.

It has been speculated that the zero-point energy of the vacuum, regularized due to the existence of a suitable ultraviolet cut-off scale, could be the source of the non-vanishing cosmological constant that is driving the present acceleration of the universe. We show that the presence of such a cut-off can significantly alter the results for the Casimir force between parallel conducting plates and even lead to repulsive Casimir force when the plate separation is smaller than the cut-off scale length. Using the current experimental data we rule out the possibility that the observed cosmological constant arises from the zero-point energy which is made finite by a suitable cut-off. Any such cut-off which is consistent with the observed Casimir effect will lead to an energy density which is about 10^{12} times larger than the observed one, if gravity couples to these modes. The implications are discussed.

In the local universe, the masses of Super-Massive Black-Holes (SMBH) appear to correlate with the physical properties of their hosts, including the mass of the dark-matter halo. Using these clues as a starting point many studies have produced models that can explain phenomena like the quasar luminosity function. The shortcoming of this approach is that working models are not unique, and as a result it is not always clear what input physics is being constrained. Here we take a different approach. We identify critical parameters that describe the evolution of SMBHs at high redshift, and constrain their parameter space based on observations of high redshift quasars from the Sloan Digital Sky Survey. We find that the luminosity function taken in isolation is somewhat limited in its ability to constrain SMBH evolution due to some strong degeneracies. This explains the presence in the literature of a range of equally successful models based on different physical hypotheses. Including the constraint of the local SMBH to halo mass ratio breaks some of the degeneracies, and our results suggest halo masses at z~4.8 of 10^{12.5+/-0.3}M_solar (with 90% confidence), with a SMBH to halo mass ratio that decreases with time (>99%). We also find a quasar luminosity to halo mass ratio that increases with halo mass (>99%). These features need to be incorporated in all successful models of SMBH evolution. On the other hand current observations do not permit any conclusions regarding the evolution of quasar lifetime, or the SMBH occupation fraction in dark matter halos.

These lecture notes provide a concise, rapid and pedagogical introduction to several advanced topics in contemporary cosmology. The discussion of thermal history of the universe, linear perturbation theory, theory of CMBR temperature anisotropies and the inflationary generation of perturbation are presented in a manner accessible to someone who has done a first course in cosmology. The discussion of dark energy is more research oriented and reflects the personal bias of the author. Contents: (I) The cosmological paradigm and Friedmann model (II) Thermal history of the universe (III) Structure formation and linear perturbation theories (IV) Perturbations in dark matter and radiation (V) Transfer function for matter perturbations (VI) Temperature anisotropies of CMBR (VII) Generation of initial perturbations from inflation (VIII) The dark energy.

We study the non linear gravitational clustering of collisionless particles in an expanding background using an integro-differential equation for the gravitational potential. In particular, we address the question of how the nonlinear mode-mode coupling transfers power from one scale to another in the Fourier space if the initial power spectrum is sharply peaked at a given scale. We show that the dynamical equation allows self similar evolution for the gravitational potential $\phi_{\bf k}(t)$ in Fourier space of the form $\phi_{\bf k}(t) = F(t)D({\bf k}) $ where the function $F(t)$ satisfies a second order non-linear differential equation. We provide a complete analysis of the relevant solutions of this equation thereby determining the asymptotic time evolution of the gravitational potential and density contrast. The analysis shows that both $F(t) $ and $D({\bf k})$ have well defined asymptotic forms indicating that the power transfer leads to a universal power spectrum at late times. The analytic results are compared with numerical simulations showing good agreement.

This article addresses some issues related to statistical description of gravitating systems in an expanding backgrounds. In particular, I describe (a) how the non linear mode-mode coupling transfers power from one scale to another in the Fourier space if the initial power spectrum is sharply peaked at a given scale and (b) what are the asymptotic characteristics of gravitational clustering that are independent of the initial conditions. The analysis uses a new approach based on an integro-differential equation for the evolution of the gravitational potential in the Fourier space. I show how this equation allows one to understand several aspects of nonlinear gravitational clustering and provides insight in to the transfer of power from one scale to another through nonlinear mode coupling. Numerical simulations as well as analytic work shows that power transfer leads to a universal power spectrum at late times, somewhat reminiscent of the existence of Kolmogorov spectrum in fluid turbulence.

A new approach to study the nonlinear phase of gravitational clustering in an expanding universe is explored. This approach is based on an integro-differential equation, for the evolution of the gravitational potential in the Fourier space, which is obtained by using a physically motivated closure condition. I show how this equation allows one to understand several aspects of nonlinear gravitational clustering and provides insight in to the transfer of power from one scale to another through nonlinear mode coupling.

Observations suggest that nearly seventy per cent of the energy density in the universe is unclustered and exerts negative pressure. Theoretical understanding of this component (`dark energy'), which is driving an accelerated expansion of the universe, is {\it the} problem in cosmology today. I discuss this issue with special emphasis on the cosmological constant as the possible choice for the dark energy. Several curious features of a universe with a cosmological constant are described and some possible approaches to understand the nature of the cosmological constant are discussed.

The dark energy component of the universe is often interpreted either in terms of a cosmological constant or as a scalar field. A generic feature of the scalar field models is that the equation of state parameter w= P/rho for the dark energy need not satisfy w=-1 and, in general, it can be a function of time. Using the Markov chain Monte Carlo method we perform a critical analysis of the cosmological parameter space, allowing for a varying w. We use constraints on w(z) from the observations of high redshift supernovae (SN), the WMAP observations of CMB anisotropies and abundance of rich clusters of galaxies. For models with a constant w, the LCDM model is allowed with a probability of about 6% by the SN observations while it is allowed with a probability of 98.9% by WMAP observations. The LCDM model is allowed even within the context of models with variable w: WMAP observations allow it with a probability of 99.1% whereas SN data allows it with 23% probability. The SN data, on its own, favors phantom like equation of state (w<-1) and high values for Omega_NR. It does not distinguish between constant w (with w<-1) models and those with varying w(z) in a statistically significant manner. The SN data allows a very wide range for variation of dark energy density, e.g., a variation by factor ten in the dark energy density between z=0 and z=1 is allowed at 95% confidence level. WMAP observations provide a better constraint and the corresponding allowed variation is less than a factor of three. Allowing for variation in w has an impact on the values for other cosmological parameters in that the allowed range often becomes larger. (Abridged)

We investigate non-linear scaling relations for two-dimensional gravitational collapse in an expanding background using a 2D TreePM code and study the strongly non-linear regime ($\bar\xi \leq 200$) for power law models. Evolution of these models is found to be scale-invariant in all our simulations. We find that the stable clustering limit is not reached, but there is a model independent non-linear scaling relation in the asymptotic regime. This confirms results from an earlier study which only probed the mildly non-linear regime ($\bar\xi \leq 40$). The correlation function in the extremely non-linear regime is a less steep function of scale than reported in earlier studies. We show that this is due to coherent transverse motions in massive haloes. We also study density profiles and find that the scatter in the inner and outer slopes is large and that there is no single universal profile that fits all cases. We find that the difference in typical density profiles for different models is smaller than expected from similarity solutions for halo profiles and transverse motions induced by substructure are a likely reason for this difference being small.

Last couple of decades have been the golden age for cosmology. High quality data confirmed the broad paradigm of standard cosmology but have thrusted upon us a preposterous composition for the universe which defies any simple explanation, thereby posing probably the greatest challenge theoretical physics has ever faced. Several aspects of these developments are critically reviewed, concentrating on conceptual issues and open questions. [Topics discussed include: Cosmological Paradigm, Growth of structures in the universe, Inflation and generation of initial perturbations, Temperature anisotropies of the CMBR, Dark energy, Cosmological Constant, Deeper issues in cosmology.]

Certain peculiar features of Einstein-Hilbert (EH) action provide clues towards a holographic approach to gravity which is independent of the detailed microstructure of spacetime. These features of the EH action include: (a) the existence of second derivatives of dynamical variables; (b) a non trivial relation between the surface term and the bulk term; (c) the fact that surface term is non analytic in the coupling constant, when gravity is treated as a spin-2 perturbation around flat spacetime and (d) the form of the variation of the surface term under infinitesimal coordinate transformations. The surface term can be derived directly from very general considerations and using (d) one can obtain Einstein's equations {\it just from the surface term of the action}. Further one can relate the bulk term to the surface term and derive the full EH action based on purely thermodynamic considerations. The features (a), (b) and (c) above emerge in a natural fashion in this approach. It is shown that action $A_{grav}$ splits into two terms $-S+\beta E$ in a natural manner \textit{in any stationary spacetime with horizon}, where $E$ is essentially an integral over ADM energy density and $S$ arises from the integral of the surface gravity over the horizon. This analysis shows that the true degrees of freedom of gravity reside in the surface term of the action, making gravity intrinsically holographic. It also provides a close connection between gravity and gauge theories, and highlights the subtle role of the singular coordinate transformations.

Spacetimes with horizons show a resemblance to thermodynamic systems and it is possible to associate the notions of temperature and entropy with them. Several aspects of this connection are reviewed in a manner appropriate for broad readership. The approach uses two essential principles: (a) the physical theories must be formulated for each observer entirely in terms of variables any given observer can access and (b) consistent formulation of quantum field theory requires analytic continuation to the complex plane. These two principles, when used together in spacetimes with horizons, are powerful enough to provide several results in a unified manner. Since spacetimes with horizons have a generic behaviour under analytic continuation, standard results of quantum field theory in curved spacetimes with horizons can be obtained directly (Sections III to VII). The requirements (a) and (b) also put strong constraints on the action principle describing the gravity and, in fact, one can obtain the Einstein-Hilbert action from the thermodynamic considerations. The latter part of the review (Sections VIII to X) investigates this deeper connection between gravity, spacetime microstructure and thermodynamics of horizons. This approach leads to several interesting results in the semiclassical limit of quantum gravity, which are described.

In any spacetime, it is possible to have a family of observers following a congruence of timelike curves such that they do not have access to part of the spacetime. This lack of information suggests associating a (congruence dependent) notion of entropy with the horizon that blocks the information from these observers. While the blockage of information is absolute in classical physics, quantum mechanics will allow tunneling across the horizon. This process can be analysed in a simple, yet general, manner and we show that the probability for a system with energy $E$ to tunnel across the horizon is $P(E)\propto\exp[-(2\pi/\kappa)E)$ where $\kappa$ is the surface gravity of the horizon. If the surface gravity changes due to the leakage of energy through the horizon, then one can associate an entropy $S(M)$ with the horizon where $dS = [ 2\pi / \kappa (M) ] dM$ and $M$ is the active gravitational mass of the system. Using this result, we discuss the conditions under which, a small patch of area $\Delta A$ of the horizon contributes an entropy $(\Delta A/4L_P^2)$, where $L_P^2$ is the Planck area.

In the local universe, the masses of Super-Massive Black-Holes (SMBH) appear to correlate with physical properties of their hosts, including the mass of the dark-matter halos. At higher redshifts, we observe the growth of SMBHs indirectly through identification of high redshift quasars. However information on their hosts is difficult to obtain. In this paper we determine the masses of halos that host high redshift quasars (at z>4) by comparing the rate of growth of quasar density with that predicted by the Press-Schechter mass function. The host mass determined depends on how the ratio between SMBH and host halo mass evolves with redshift. Assuming the ratio between SMBH and halo mass does not evolve with redshift, we find a host halo mass of log(M)=11.7+/-0.3 solar masses. Even if the quasars shine at their Eddington limit, this mass is significantly smaller than that seen in the local universe. Indeed the null-hypothesis, of a constant SMBH to halo mass ratio at all redshifts, can be ruled out at greater than a 5-sigma level. SMBHs must therefore have contributed a larger fraction to the host mass in the past. Including redshift evolution of the SMBH to halo mass ratio, we find larger halo masses of log(M)=12.4+/-0.3 solar masses, and a ratio between SMBH and halo mass that increases with redshift as ~(1+z)^1.5 are required to be consistent with both local and high redshift observations. We investigate restrictions placed on the critical linear overdensity of quasar hosts at their epoch of virialisation and find that it cannot exceed the traditional value of delta_c=1.69 by more than a factor of two. Finally, we find that the high redshift quasars are hosted by fluctuations on scales that have a variance of (delta M/M) = 2-3, corresponding to (3-4.5)-sigma fluctuations in the density field.

In this review we discuss Cosmological N-Body codes with a special emphasis on Particle Mesh codes. We present the mathematical model for each component of N-Body codes. We compare alternative methods for computing each quantity by calculating errors for each of the components. We suggest an optimum set of components that can be combined reduce overall errors in N-Body codes.

The conceptual difficulties associated with a cosmological constant have led to the investigation of alternative models in which the equation of state parameter, $w=p/\rho$, of the dark energy evolves with time. We show that combining the supernova type Ia observations {\it with the constraints from WMAP observations} restricts large variation of $\rho(z)$ at low redshifts. The combination of these two observational constraints is stronger than either one. The results are completely consistent with the cosmological constant as the source of dark energy.

Recent cosmological observations suggest that nearly seventy per cent of the energy density in the universe is unclustered and has negative pressure. Several conceptual issues related to the modeling of this component (`dark energy'), which is driving an accelerated expansion of the universe, are discussed with special emphasis on the cosmological constant as the possible choice for the dark energy. Some curious geometrical features of a universe with a cosmological constant are described and a few attempts to understand the nature of the cosmological constant are reviewed.

There is a general belief, reinforced by statements in standard textbooks, that: (i) one can obtain the full non-linear Einstein's theory of gravity by coupling a massless, spin-2 field $h_{ab}$ self-consistently to the total energy momentum tensor, including its own; (ii) this procedure is unique and leads to Einstein-Hilbert action and (iii) it only uses standard concepts in Lorentz invariant field theory and does not involve any geometrical assumptions. After providing several reasons why such beliefs are suspect -- and critically re-examining several previous attempts -- we provide a detailed analysis aimed at clarifying the situation. First, we prove that it is \textit{impossible} to obtain the Einstein-Hilbert (EH) action, starting from the standard action for gravitons in linear theory and iterating repeatedly. Second, we use the Taylor series expansion of the action for Einstein's theory, to identify the tensor $\mathcal{S}^{ab}$, to which the graviton field $h_{ab}$ couples to the lowest order. We show that the second rank tensor $\mathcal{S}^{ab}$ is {\it not} the conventional energy momentum tensor $T^{ab}$ of the graviton and provide an explanation for this feature. Third, we construct the full nonlinear Einstein's theory with the source being spin-0 field, spin-1 field or relativistic particles by explicitly coupling the spin-2 field to this second rank tensor $\mathcal{S}^{ab}$ order by order and summing up the infinite series. Finally, we construct the theory obtained by self consistently coupling $h_{ab}$ to the conventional energy momentum tensor $T^{ab}$ order by order and show that this does {\it not} lead to Einstein's theory. (condensed).

We extend our previous analysis of cosmological supernova data (Padmanabhan & Choudhury 2003) to include three recent compilation of data sets. Our analysis ignores systematic effects in the data and concentrates on some key theoretical issues. The first data set consists of 194 points while the second discards some of the points from the first one because of large uncertainties and thus consists of 142 points. The third data set is obtained from the second by adding the latest 14 points observed through HST. A careful comparison of these different data sets help us to draw the following conclusions: (i) All the data sets strongly rule out non-accelerating models. Interestingly, the first and the second sets favour a closed universe; i.e., the probability P(Omega_{tot} > 1) \gtrsim 0.97. Hence these sets are in mild disagreement with the ``concordance'' flat model. However, this disagreement is reduced [P(Omega_{tot} > 1) \approx 0.9] for the third data set, which includes the points observed by HST around 1 < z < 1.6. (ii) When the first data set is divided into two separate subsets consisting of low (z < 0.34) and high (z > 0.34) redshift supernova, it turns out that these two subsets, individually, admit non-accelerating models. However, these models seem to be ruled out using only the low redshift data for the other two data sets, which have less uncertainties. (iii) It is quite difficult to measure the evolution of the dark energy equation of state w_X though its present value can be constrained quite well. The best-fit value seems to mildly favour a dark energy component with current equation of state w_X < -1, thus opening the possibility of existence of more exotic forms of matter. However, the data is still consistent with the the standard cosmological constant at 99% confidence level (abridged).

It is very likely that the quantum description of spacetime is quite different from what we perceive at large scales, $l\gg (G\hbar/c^3)^{1/2}$. The long wave length description of spacetime, based on Einstein's equations, is similar to the description of a continuum solid made of a large number of microscopic degrees of freedom. This paradigm provides a novel interpretation of coordinate transformations as deformations of "spacetime solid" and allows one to obtain Einstein's equations as a consistency condition in the long wavelength limit. The entropy contributed by the microscopic degrees of freedom reduces to a pure surface contribution when Einstein's equations are satisfied. The horizons arises as "defects" in the "spacetime solid" (in the sense of well defined singular points) and contributes an entropy which is one quarter of the horizon area. Finally, the response of the microstructure to vacuum energy leads to a near cancellation of the cosmological constant, leaving behind a tiny fluctuation which matches with the observed value.

The standard, scale-invariant, inflationary perturbation spectrum will be modified if the effects of trans-Planckian physics are incorporated into the dynamics of the matter field in a phenomenological manner, say, by the modification of the dispersion relation. The spectrum also changes if we retain the standard dynamics but modify the initial quantum state of the matter field. We show that, given {\it any} spectrum of perturbations, it is possible to choose a class of initial quantum states which can exactly reproduce this spectrum with the standard dynamics. We provide an explicit construction of the quantum state which will produce the given spectrum. We find that the various modified spectra that have been recently obtained from `trans-Planckian considerations' can be constructed from suitable squeezed states above the Bunch-Davies vacuum in the standard theory. Hence, the CMB observations can, at most, be useful in determining the initial state of the matter field in the standard theory, but it can {\it not} help us to discriminate between the various Planck scale models of matter fields. We study the problem in the Schrodinger picture, clarify various conceptual issues and determine the criterion for negligible back reaction due to modified initial conditions.

The energy density associated with Planck length is $\rho_{uv}\propto L_P^{-4}$ while the energy density associated with the Hubble length is $\rho_{ir}\propto L_H^{-4}$ where $L_H=1/H$. The observed value of the dark energy density is quite different from {\it either} of these and is close to the geometric mean of the two: $\rho_{vac}\simeq \sqrt{\rho_{uv} \rho_{ir}}$. It is argued that classical gravity is actually a probe of the vacuum {\it fluctuations} of energy density, rather than the energy density itself. While the globally defined ground state, being an eigenstate of Hamiltonian, will not have any fluctuations, the ground state energy in the finite region of space bounded by the cosmic horizon will exhibit fluctuations $\Delta\rho_{\rm vac}(L_P, L_H)$. When used as a source of gravity, this $\Delta \rho$ should lead to a spacetime with a horizon size $L_H$. This bootstrapping condition leads naturally to an effective dark energy density $\Delta\rho\propto (L_{uv}L_H)^{-2}\propto H^2/G$ which is precisely the observed value. The model requires, either (i) a stochastic fluctuations of vacuum energy which is correlated over about a Hubble time or (ii) a semi- anthropic interpretation. The implications are discussed.

In case of spacetimes with single horizon, there exist several well-established procedures for relating the surface gravity of the horizon to a thermodynamic temperature. Such procedures, however, cannot be extended in a straightforward manner when a spacetime has multiple horizons. In particular, it is not clear whether there exists a notion of global temperature characterizing the multi-horizon spacetimes. We examine the conditions under which a global temperature can exist for a spacetime with two horizons using the example of Schwarzschild-De Sitter (SDS) spacetime. We systematically extend different procedures (like the expectation value of stress tensor, response of particle detectors, periodicity in the Euclidean time etc.) for identifying a temperature in the case of spacetimes with single horizon to the SDS spacetime. This analysis is facilitated by using a global coordinate chart which covers the entire SDS manifold. We find that all the procedures lead to a consistent picture characterized by the following features: (a) In general, SDS spacetime behaves like a non-equilibrium system characterized by two temperatures. (b) It is not possible to associate a global temperature with SDS spacetime except when the ratio of the two surface gravities is rational (c) Even when the ratio of the two surface gravities is rational, the thermal nature depends on the coordinate chart used. There exists a global coordinate chart in which there is global equilibrium temperature while there exist other charts in which SDS behaves as though it has two different temperatures. The coordinate dependence of the thermal nature is reminiscent of the flat spacetime in Minkowski and Rindler coordinate charts. The implications are discussed.

It is known that the imaginary parts of the quasi normal mode (QNM) frequencies for the Schwarzschild black hole are evenly spaced with a spacing that depends only on the surface gravity. On the other hand, for massless minimally coupled scalar fields, there exist no QNMs in the pure DeSitter spacetime. It is not clear what the structure of the QNMs would be for the Schwarzschild-DeSitter (SDS) spacetime, which is characterized by two different surface gravities. We provide a simple derivation of the imaginary parts of the QNM frequencies for the SDS spacetime by calculating the scattering amplitude in the first Born approximation and determining its poles. We find that, for the usual set of boundary conditions in which the incident wave is scattered off the black hole horizon, the imaginary parts of the QNM frequencies have a equally spaced structure with the level spacing depending on the surface gravity of the black hole. Several conceptual issues related to the QNM are discussed in the light of this result and comparison with previous work is presented.

A class of metrics $g_{ab}(x^i)$ describing spacetimes with horizons (and associated thermodynamics) can be thought of as a limiting case of a family of metrics $g_{ab}(x^i;\lambda)$ {\it without horizons} when $\lambda\to 0$. I construct specific examples in which the curvature corresponding $g_{ab}(x^i;\lambda)$ becomes a Dirac delta function and gets concentrated on the horizon when the limit $\lambda\to 0$ is taken, but the action remains finite. When the horizon is interpreted in this manner, one needs to remove the corresponding surface from the Euclidean sector, leading to winding numbers and thermal behaviour. In particular, the Rindler spacetime can be thought of as the limiting case of (horizon-free) metrics of the form [$g_{00}=\epsilon^2+a^2x^2; g_{\mu\nu}=-\delta_{\mu\nu}$] or [$g_{00} = - g^{xx} = (\epsilon^2 +4 a^2 x^2)^{1/2}, g_{yy}=g_{zz}=-1]$ when $\epsilon\to 0$. In the Euclidean sector, the curvature gets concentrated on the origin of $t_E-x$ plane in a manner analogous to Aharanov-Bohm effect (in which the the vector potential is a pure gauge everywhere except at the origin) and the curvature at the origin leads to nontrivial topological features and winding number.

It is known that the imaginary parts of the frequencies of the quasi normal modes of the Schwarzschild black hole are equally spaced, with the level spacing dependent only on the surface gravity. We generalize this result to a wider class of spacetimes and provide a simple derivation of the imaginary parts of the frequencies. The analysis shows that the result is closely linked to the thermal nature of horizons and arises from the exponential redshift of the wave modes close to the horizon.

In any space-time, it is possible to have a family of observers who have access to only part of the space-time manifold, because of the existence of a horizon. We demand that \emph{physical theories in a given coordinate system must be formulated entirely in terms of variables that an observer using that coordinate system can access}. In the coordinate frame in which these observers are at rest, the horizon manifests itself as a (coordinate) singularity in the metric tensor. Regularization of this singularity removes the inaccessible region, and leads to the following consequences: (a) The non-trivial topological structure for the effective manifold allows one to obtain the standard results of quantum field theory in curved space-time. (b) In case of gravity, this principle requires that the effect of the unobserved degrees of freedom should reduce to a boundary contribution $A_{\rm boundary}$ to the gravitational action. When the boundary is a horizon, $A_{\rm boundary}$ reduces to a single, well-defined term proportional to the area of the horizon. Using the form of this boundary term, it is possible to obtain the full gravitational action in the semiclassical limit. (c) This boundary term must have a quantized spectrum with uniform spacing, $\Delta A_{boundary}=2\pi\hbar$, in the semiclassical limit. This, in turn, yields the following results for semiclassical gravity: (i) The area of any one-way membrane is quantized. (ii) The information hidden by a one-way membrane amounts to an entropy, which is always one-fourth of the area of the membrane in the leading order. (iii) In static space-times, the action for gravity can be given a purely thermodynamic interpretation and the Einstein equations have a formal similarity to laws of thermodynamics.

A brief summary of several topics in the study of gravitational many body problem is given. The discussion covers both static backgrounds (applicable to astrophysical systems) as well as clustering in an expanding background (relevant for cosmology)

Current cosmological observations show a strong signature of the existence of a dark energy component with negative pressure. The most obvious candidate for this dark energy is the cosmological constant (with the equation of state w_X=p/\rho=-1), which, however, raises several theoretical difficulties. This has led to models for dark energy component which evolves with time. We discuss certain questions related to the determination of the nature of dark energy component from observations of high redshift supernova. The main results of our analysis are: (i) Even if the precise value of w_X is known from observations, it is not possible to determine the nature of the unknown dark energy source using only kinematical and geometrical measurements. We have given explicit examples to show that different types of sources can give rise to a given w_X. (ii) Although the full data set of supernova observations (which are currently available) strongly rule out models without dark energy, the high (z>0.25) and low (z<0.25) redshift data sets, individually, admit decelerating models with zero dark energy. Any possible evolution in the absolute magnitude of the supernovae, if detected, might allow the decelerating models to be consistent with the data. (iii) We have introduced two parameters, which can be obtained entirely from theory, to study the sensitivity of the luminosity distance on w_X. Using these two parameters, we have argued that although one can determine the present value of w_X accurately from the data, one cannot constrain the evolution of w_X.

We show that the distance measures (such as the luminosity and angular diameter distances) are linear functionals of the equation-of-state function w(z) of the dark energy to a fair degree of accuracy in the regimes of interest. That is, the distance measures can be expressed as a sum of (i) a constant and (ii) an integral of a weighting function multiplied by the equation of state parameter w(z). The existence of such an accurate linear response approximation has several important implications: (a) Fitting a constant w model to the data drawn from an evolving model has a simple interpretation as a weighted average of w(z). (b) Any polynomial (or other expansion coefficients can also be expressed as weighted sums of the true w(z). (c) A replacement for the commonly used heuristic equation for the effective w(z), as determined by the CMB, can be derived and the result is found to be quite close to the heuristic expression commonly used. (d) The reconstruction of w(z) by Huterer et al. (2002) can be expressed as a matrix inversion. In each case the limitations of the linear response approximation are explored and found to be surprisingly small.

The principle of equivalence provides a description of gravity in terms of the metric tensor and determines how gravity affects the light cone structure of the space-time. This, in turn, leads to the existence of observers (in any space-time) who do not have access to regions of space-time bounded by horizons. To take into account this generic possibility, it is necessary to demand that \emph{physical theories in a given coordinate system must be formulated entirely in terms of variables that an observer using that coordinate system can access}. This principle is powerful enough to obtain the following results: (a) The action principle of gravity must be of such a structure that, in the semiclassical limit, the action of the unobserved degrees of freedom reduces to a boundary contribution $A_{\rm boundary}$ obtained by integrating a four divergence. (b) When the boundary is a horizon, $A_{\rm boundary}$ essentially reduces to a single, well-defined, term. (c) This boundary term must have a quantized spectrum with uniform spacing, $\Delta A_{boundary}=2\pi\hbar$, in the semiclassical limit. Using this principle in conjunction with the usual action principle in gravity, we show that: (i) The area of any one-way membrane is quantized. (ii) The information hidden by a one-way membrane leads to an entropy which is always one-fourth of the area of the membrane, in the leading order. (iii) In static space-times, the action for gravity can be given a purely thermodynamic interpretation and the Einstein equations have a formal similarity to laws of thermodynamics.

We present a detailed study of cosmological effects of homogeneous tachyon matter coexisting with non-relativistic matter and radiation, concentrating on the inverse square potential and the exponential potential for the tachyonic scalar field. A distinguishing feature of these models (compared to other cosmological models) is that the matter density parameter and the density parameter for tachyons remain comparable even in the matter dominated phase. For the exponential potential, the solutions have an accelerating phase, followed by a phase with a(t) ~ t^{2/3} as t goes to infinity. This eliminates the future event horizon present in LCDM models and is an attractive feature from the string theory perspective. A comparison with supernova Ia data shows that for both the potentials there exists a range of models in which the universe undergoes an accelerated expansion at low redshifts and are also consistent with requirements of structure formation. They do require fine tuning of parameters but not any more than in the case of LCDM or quintessence models.

Recent cosmological observations suggest the existence of a positive cosmological constant $\Lambda$ with the magnitude $\Lambda(G\hbar/c^3) \approx 10^{-123}$. This review discusses several aspects of the cosmological constant both from the cosmological (sections 1-6) and field theoretical (sections 7-11) perspectives. The first section introduces the key issues related to cosmological constant and provides a brief historical overview. This is followed by a summary of the kinematics and dynamics of the standard Friedmann model of the universe paying special attention to features involving the cosmological constant. Section 3 reviews the observational evidence for cosmological constant, especially the supernova results, constraints from the age of the universe and a few others. Theoretical models (quintessence, tachyonic scalar field, ...) with evolving cosmological `constant' are described from different perspectives in the next section. Constraints on dark energy from structure formation and from CMBR anisotropies are discussed in the next two sections. The latter part of the review (sections 7-11) concentrates on more conceptual and fundamental aspects of the cosmological constant. Section 7 provides some alternative interpretations of the cosmological constant which could have a bearing on the possible solution to the problem. Several relaxation mechanisms have been suggested in the literature to reduce the cosmological constant to the currently observed value and some of these attempts are described in section 8. Next section gives a brief description of the geometrical structure of the de Sitter spacetime and the thermodynamics of the de Sitter universe is taken up in section 10. The last section deals with the role of string theory in the cosmological constant problem.

We study the cosmological evolution of a scalar field that couples to the trace $T=T^{a}_a$ of energy momentum tensor of all the fields (including itself). In the case of a shallow exponential potential, the presence of coupling to the trace $T$ in the field equation makes the energy density of the scalar field decrease faster thereby hastening the commencement of radiation domination. This effect gradually diminishes at later epochs allowing the scalar field to dominate the energy density again. We interpret this phase as the current epoch of cosmic acceleration with $\Omega_{\phi}=0.7$. A variant of this model can lead to accelerated expansion at the present epoch followed by a $a(t)\propto t^{2/3}$ behaviour as $t\to \infty$, making the model free from future event horizon. The main features of the model are independent of initial conditions. However, fine tuning of parameters is necessary for viable evolution.

The Einstein-Hilbert action (and thus the dynamics of gravity) can be obtained by combining the principle of equivalence, special relativity and quantum theory in the Rindler frame and postulating that the horizon area must be proportional to the entropy. This approach uses the local Rindler frame as a natural extension of the local inertial frame, and leads to the interpretation that the gravitational action represents the free energy of the spacetime geometry. As an aside, one obtains an insight into the peculiar structure of Einstein-Hilbert action and a natural explanation to the questions:(i) Why does the covariant action for gravity contain second derivatives of the metric tensor? (ii) Why is the gravitational coupling constant is positive ? Some geometrical features of gravitational action are clarified.

Cosmological observations suggest the existence of two different kinds of energy densities dominating at small ($ \lesssim 500$ Mpc) and large ($\gtrsim 1000 $ Mpc) scales. The dark matter component, which dominates at small scales, contributes $\Omega_m \approx 0.35$ and has an equation of state $p=0$ while the dark energy component, which dominates at large scales, contributes $\Omega_V \approx 0.65$ and has an equation of state $p\simeq -\rho$. It is usual to postulate wimps for the first component and some form of scalar field or cosmological constant for the second component. We explore the possibility of a scalar field with a Lagrangian $L =- V(\phi) \sqrt{1 - \del^i \phi \del_i \phi}$ acting as {\it both} clustered dark matter and smoother dark energy and having a scale dependent equation of state. This model predicts a relation between the ratio $ r = \rho_V/\rho_{\rm DM}$ of the energy densities of the two dark components and expansion rate $n$ of the universe (with $a(t) \propto t^n$) in the form $n = (2/3) (1+r) $. For $r \approx 2$, we get $n \approx 2$ which is consistent with observations.

A general formalism for understanding the thermodynamics of horizons in spherically symmetric spacetimes is developed. The formalism reproduces known results in the case of black hole spacetimes. But its power lies in being able to handle more general situations like: (i) spacetimes which are not asymptotically flat (like the de Sitter spacetime) and (ii) spacetimes with multiple horizons having different temperatures (like the Schwarzschild-de Sitter spacetime) and provide a consistent interpretation for temperature, entropy and energy. I show that it is possible to write Einstein's equations for a spherically symmetric spacetime in the form $TdS-dE=PdV$ near {\it any} horizon of radius $a$ with $S=(1/4)(4\pi a^2), |E| = (a/2)$ and the temperature $T$ determined from the surface gravity at the horizon. The pressure $P$ is provided by the source of the Einstein's equations and $dV$ is the change in the volume when the horizon is displaced infinitesimally. The same results can be obtained by evaluating the quantum mechanical partition function {\it without using Einstein's equations or WKB approximation for the action}. Both the classical and quantum analysis provide a simple and consistent interpretation of entropy and energy for de Sitter spacetime as well as for $(1+2)$ dimensional gravity. For the Rindler spacetime the entropy per unit transverse area turns out to be $(1/4)$ while the energy is zero. The approach also shows that the de Sitter horizon -- like the Schwarzschild horizon -- is effectively one dimensional as far as the flow of information is concerned, while the Schwarzschild-de Sitter, Reissner-Nordstrom horizons are not. The implications for spacetimes with multiple horizons are discussed.

The current observations seem to suggest that the universe has a positive cosmological constant of the order of $H_0^2$ while the most natural value for the cosmological constant will be $L_P^{-2}$ where $L_P = (G\hbar/c^3)^{1/2}$ is the Planck length. This reduction of the cosmological constant from $L_P^{-2} $ to $L_P^{-2}(L_PH_0)^2$ may be interpreted as due to the ability of quantum micro structure of spacetime to readjust itself and absorb bulk vacuum energy densities. Being a quantum mechanical process, such a cancellation cannot be exact and the residual quantum fluctuations appear as the ``small'' cosmological constant. I describe the features of a toy model for the spacetime micro structure which could allow for the bulk vacuum energy densities to be canceled leaving behind a small residual value of the the correct magnitude. Some other models (like the ones based on canonical ensemble for the four volume or quantum fluctuations of the horizon size) lead to an insignificantly small value of $H_0^2(L_PH_0)^n$ with $n=0.5-1$ showing that obtaining the correct order of magnitude for the residual fluctuations in the cosmological constant is a nontrivial task, becaue of the existence of the small dimensionless number $H_0L_P$ .

A large class of evolutionary processes can be modeled by a rule which involves self-replication of some physical quantity with a non local rescaling. I show that a class of such models are exactly solvable -- in the discrete as well as continuum limit -- and can represent several physical situations as varied from the formation of galaxies in some cosmological models to growth of bacterial cultures. This class of models, in general, has no steady state solution and evolve unstably as $t \to \infty$ for generic initial conditions. They can however exhibit (unstable) power law correlation function in the continuum limit, for an intermediate range of times and length scales.

This pedagogical review addresses several issues related to statistical description of gravitating systems in both static and expanding backgrounds, focusing on the latter. After briefly reviewing the results for the static background, I describe the key issues and recent progress in the context of gravitational clustering of collision-less particles in an expanding universe. The questions addressed include: (a) How does the power injected into the system at a given wave number spread to smaller and larger scales? (b) How does the power spectrum of density fluctuations behave asymptotically at late times? (c) What are the universal characteristics of gravitational clustering that are independent of the initial conditions and manifest at the late time evolution of the system? The review is intended for non cosmologists and will be of interest to people working in fluid mechanics, non linear dynamics and condensed matter physics.

The two surprising features of gravity are (a) the principle of equivalence and (b) the connection between gravity and thermodynamics. Using principle of equivalence and special relativity in the {\it local inertial frame}, one could obtain the insight that gravity must possess a geometrical description. I show that, using the same principle of equivalence, special relativity and quantum theory in the {\it local Rindler frame} one can obtain the Einstein-Hilbert action functional for gravity and thus the dynamics of the spacetime. This approach, which essentially involves postulating that the horizon area must be proportional to the entropy, uses the local Rindler frame as a natural extension of the local inertial frame and leads to the interpretation that the gravitational action represents the free energy of the spacetime geometry. As an aside, one also obtains a natural explanation as to: (i) why the covariant action for gravity contains second derivatives of the metric tensor and (ii) why the gravitational coupling constant is positive. The analysis suggests that gravity is intrinsically holographic and even intrinsically quantum mechanical.

I provide a general proof of the conjecture that one can attribute an entropy to the area of {\it any} horizon. This is done by constructing a canonical ensemble of a subclass of spacetimes with a fixed value for the temperature $T=\beta^{-1}$ and evaluating the {\it exact} partition function $Z(\beta)$. For spherically symmetric spacetimes with a horizon at $r=a$, the partition function has the generic form $Z\propto \exp[S-\beta E]$, where $S= (1/4) 4\pi a^2$ and $|E|=(a/2)$. Both $S$ and $E$ are determined entirely by the properties of the metric near the horizon. This analysis reproduces the conventional result for the black-hole spacetimes and provides a simple and consistent interpretation of entropy and energy for De Sitter spacetime. For the Rindler spacetime the entropy per unit transverse area turns out to be $(1/4)$ while the energy is zero. Further, I show that the relationship between entropy and area allows one to construct the action for the gravitational field on the bulk and thus the full theory. In this sense, gravity is intrinsically holographic.

It is an accepted practice in cosmology to invoke a scalar field with potential $V(\phi)$ when observed evolution of the universe cannot be reconciled with theoretical prejudices. Since one function-degree-of-freedom in the expansion factor $a(t)$ can be traded off for the function $V(\phi)$, it is {\it always} possible to find a scalar field potential which will reproduce a given evolution. I provide a recipe for determining $V(\phi)$ from $a(t)$ in two cases:(i) Normal scalar field with Lagrangian ${\cal L} = (1/2)\partial_a\phi \partial^a\phi - V(\phi)$ used in quintessence/dark energy models. (ii) A tachyonic field with Lagrangian ${\cal L} = -V(\phi) [ 1- \partial_a\phi \partial^a\phi]^{1/2} $, motivated by recent string theoretic results. In the latter case, it is possible to have accelerated expansion of the universe during the late phase in certain cases. This suggests a string theory based interpretation of the current phase of the universe with tachyonic condensate acting as effective cosmological constant.

A simple analytical model for estimating the fraction (\Omega_{gas}) of matter in gaseous form within the collapsed dark matter (DM) haloes is presented. The model is developed using (i) the Press-Schechter formalism to estimate the fraction of baryons in DM haloes, and (ii) the observational estimates of the star formation rate at different redshifts. The prediction for \Omega_{gas} from the model is in broad agreement with the observed abundance of the damped Ly-\alpha systems. Furthermore, it can be used for estimating the circular velocities of the collapsed haloes at different redshifts, which could be compared with future observations.

The dynamics of our universe is characterised by the density parameters for cosmological constant ($\Omega_V$), nonbaryonic darkmatter($\Omega_{\rm wimp}$), radiation ($\Omega_R$) and baryons ($\Omega_B$). To these parameters -- which describe the smooth background universe -- one needs to add at least another dimensionless number ($\sim 10^{-5}$) characterising the strength of primordial fluctuations in the gravitational potential, in order to ensure formation of structures by gravitational instability. I discuss several issues related to the description of the universe in terms of these numbers and argue that we do not yet have a fundamental understanding of these issues.

We point out that a pseudo-Newtonian interpretations of the $(1/3)$ factor in the Sachs-Wolfe effect, which relates the fluctuations in temperature and potential, $(\delta T/ T) = (1/ 3) \delta \Phi$, is not supported by the General Relativistic analysis. Dividing the full gravitational effect into separate parts depends on the choice of time slicing (gauge) and there exist infinitely many different choices. More importantly, interpreting the parts as being due to the gravitational redshift and the time dilation is not justified in the rigorous relativistic perturbation theory. We suggest to regard the $(1/3)$ factor as the General Relativistic result which applies in a restricted situation of adiabatic perturbation in the $K = 0 = \Lambda$ model with the last scattering occuring in the matter dominated era. For an isocurvature initial condition the corresponding result, $(\delta T/ T) = 2 \delta \Phi$, has a different numerical coefficient.

The issues related to bringing together the principles of general relativity and quantum theory are discussed. After briefly summarising the points of conflict between the two formalisms I focus on four specific themes in which some contact has been established in the past between GR and quantum field theory: (i) The role of planck length in the microstructure of spacetime (ii) The role of quantum effects in cosmology and origin of the universe (iii) The thermodynamics of spacetimes with horizons and especially the concept of entropy related to spacetime geometry (iv) The problem of the cosmological constant.

Following Bi & Davidsen (1997), we perform one dimensional semi analytic simulations along the lines of sight to model the intergalactic medium (IGM). Since this procedure is computationally efficient in probing the parameter space -- and reasonably accurate -- we use it to recover the values of various parameters related to the IGM (for a fixed background cosmology) by comparing the model predictions with different observations. For the currently favoured LCDM model (\Omega_m=0.4, \Omega_{\Lambda}=0.6 and h=0.65), we obtain, using statistics obtained from the transmitted flux, constraints on (i) the combination f=(\Omega_B h^2)^2/J_{-12}, where \Omega_B is the baryonic density parameter and J_{-12} is the total photoionisation rate in units of 10^{-12} s^{-1}, (ii) temperature T_0 corresponding to the mean density and (iii) the slope \gamma of the effective equation of state of the IGM at a mean redshift z \simeq 2.5. We find that 0.8 <(T_0/10^4 K)< 2.5 and 1.3<\gamma<2.3. while the constraint obtained on f is 0.020^2<f<0.032^2. A reliable lower bound on J_{-12} can be used to put a lower bound on \Omega_B h^2, which can be compared with similar constraints obtained from Big Bang Nucleosynthesis (BBN) and CMBR studies. We find that if J_{-12}>1.2, the lower bound on \Omega_B h^2 is in violation of the BBN value.

We use non-linear scaling relations (NSRs) to investigate the effects arising from the existence of negative correlations on the evolution of gravitational clustering in an expanding universe. It turns out that such anti-correlated regions have important dynamical effects on {\it all} scales. In particular, the mere existence of negative values for the linear two-point correlation function $\xib_L$ over some range of scales starting from $l = L_o$, implies that the non-linear correlation function is bounded from above at {\it all} scales $x < L_o$. This also results in the relation $\xib \propto x^{-3}$, at these scales, at late times, independent of the original form of the correlation function. Current observations do not rule out the existence of negative $\xib$ for $ 200 h^{-1}$ Mpc $\la \xib \la 1000 h^{-1}$ Mpc; the present work may thus have relevance for the real Universe. The only assumption made in the analysis is the {\it existence} of the NSR; the results are independent of the form of the NSR as well as of the stable clustering hypothesis.

Analytic derivations of the correlation function and the column density distribution for neutral hydrogen in the IGM are presented, assuming that the non-linear baryonic mass density distribution in the IGM is lognormal. This ansatz was used earlier by Bi & Davidsen (1997) to perform 1D simulations of lines-of-sight and analyse the properties of absorption systems. Our approach is completely analytic, which allows us to explore a wide region of the parameter space for our model. The analytic results have been compared with observations, whenever possible. Two kinds of correlation functions are defined: along the line-of-sight (LOS) and across the transverse direction. We find that the effects on the LOS correlation due to change in cosmology and the slope of the equation of state of the IGM, \gamma are of the same order, which means that we cannot constrain both the parameters simultaneously. However, it is possible to constrain \gamma and its evolution using the observed LOS correlation function at different epochs, provided one knows the background cosmology. We suggest that the constraints on the evolution of \gamma obtained using the LOS correlation can be used as an independent tool to probe the reionisation history of the universe. From the transverse correlation function, we find that the excess probability, over random, of finding two neutral hydrogen overdense regions separated by an angle \theta, is always less than 1 per cent for redshifts greater than 2. Our models also reproduce the observed column density distribution for neutral hydrogen and the shape of the distribution depends on \gamma. Our calculations suggest that one can rule out \gamma > 1.6 for z \simeq 2.31 using the column density distribution. However, one cannot rule higher values of \gamma at higher redshifts.

We analyse the BOOMERANG and MAXIMA results in the context of models with $\Omega_{\rm Total} = 1$ and $n_s = 1$. We attempt to constrain three other parameters---$h$, $\Omega_B$, and $\Omega_m$---from these observations. We show that: (a) the value of $\Omega_B h^2$ is too high to be compatible with primordial nucleosynthesis observations at 95% confidence level (b) universe with age greater than $ 12 \rm Gyr$ is ruled out at 95% confidence level and (c) the value of $\Omega_m h$ is too high to be compatible with the shape of the power spectrum of gravitational clustering. In effect, our analysis shows that models with $\Omega_{\rm Total} = 1$ and $n_s = 1$ are ruled out by BOOMERANG and MAXIMA observations.

We investigate the evolution of non-linear density perturbations by taking into account the effects of deviations from spherical symmetry of a system. Starting from the standard spherical top hat model in which these effects are ignored, we introduce a physically motivated closure condition which specifies the dependence of the additional terms on the density contrast, $\delta$. The modified equation can be used to model the behaviour of an overdense region over a sufficiently large range of $\delta$. The key new idea is a Taylor series expansion in ($1/\delta$) to model the non-linear epoch. We show that the modified equations quite generically lead to the formation of stable structures in which the gravitational collapse is halted at around the virial radius. The analysis also allows us to connect up the behaviour of individual overdense regions with the non-linear scaling relations satisfied by the two point correlation function.

Several issues related to the gravitational clustering of collisionless dark matter in an expanding universe is discussed. The discussion is pedagogical but the emphasis is on semianalytic methods and open questions - rather than on well established results.

We consider the problem of gravitational clustering in a D-dimensional expanding Universe and derive scaling relations connecting the exact mean two-point correlation function with the linear mean correlation function, in the quasi-linear and non-linear regimes, using the standard paradigms of scale-invariant radial collapse and stable clustering. We show that the existence of scaling laws is a generic feature of gravitational clustering in an expanding background, in all dimensions except D=2 and comment on the special nature of the 2-dimensional case. The D-dimensional scaling laws derived here reduce, in the 3-dimensional case, to scaling relations obtained earlier from N-body simulations. Finally, we consider the case of clustering of 2-dimensional particles in a 2-D expanding background, governed by a force -GM/R, and show that the correlation function does not grow (to first order) until much after the recollapse of any shell.

Points of conflict between the principles of general relativity and quantum theory are highlighted. I argue that the current language of QFT is inadequete to deal with gravity and review attempts to identify some of the features which are likely to present in the correct theory of quantum gravity.

Cosmological perturbations with wavelengths smaller than Hubble radius can be handled in the context of Newtonian theory with very high accuracy. The application of this Newtonian approximation, however, is restricted to nonrelativistic matter and cannot be used for relativistic matter. Recently, by modifying the continuity equation, Lima, et. al., extended the domain of applicability of Newtonian cosmology to radiation dominated phase. We adopted this continuity equation to re-examine linear cosmological perturbation theory for a two fluid universe with uniform pressure. We study the evolution equations for density contrasts and their validity in different epochs and on scales larger than Hubble radius and compare the results with the full relativistic approach. The comparison shows the high accuracy of this approximation.

Several investigations in the study of cosmological structure formation use numerical simulations in both two and three dimensions. In this paper we address the subtle question of ambiguities in the nature of two dimensional gravity in an expanding background. We take a detailed and formal approach by deriving the equations describing gravity in (D+1) dimensions using the action principle of Einstein. We then consider the Newtonian limit of these equations and finally obtain the necessary fluid equations required to describe structure formation. These equations are solved for the density perturbation in both the linearised form and in the spherical top hat model of nonlinear growth. We find that, when the special case of D=2 is considered, no structures can grow. We therefore conclude that, within the frame work of Einstein's theory of gravity in (2+1) dimensions, formation of structures cannot take place. Finally, we indicate the different possible ways of getting around this difficulty so that growing structures can be obtained in two dimensional cosmological gravitational simulations and discuss their implications.

It is known that radial collapse around density peaks can explain the key features of evolution of correlation function in gravitational clustering in three dimensions. The same model also makes specific predictions for two dimensions. In this paper we test these predictions in two dimensions with the help of N-Body simulations. We find that there is no stable clustering in the extremely non-linear regime, but a nonlinear scaling relation does exist and can be used to relate the linear and the non-linear correlation function. In the intermediate regime, the simulations agree with the model.

Padmanabhan (1996) has suggested a model to relate the nonlinear two - point correlation function to the linear two - point correlation function. In this paper, we extend this model in two directions: (1) By averaging over the initial Gaussian distribution of density contrasts, we estimate the spectral dependence of the scaling between nonlinear and linear correlation functions. (2) By using a physically motivated ansatz, we generalise the model to N-point correlation functions and relate the nonlinear, volume averaged, N-point correlation function $\bar\xi_N(x,a)$ with linearly extrapolated volume averaged 2-point correlation function $\bar \xi_2(l,a)$ evaluated at a different scale. We compare the point of transition between different regimes obtained from our model with numerical simulations and show that the spectral dependence of the scaling relations seen in the simulations can be easily understood. Comparison of the calculated form of $\bar\xi_N$ with the simulations show reasonable agreement. We discuss several implications of the results.

We discuss the late time evolution of the gravitational clustering in an expanding universe, based on the nonlinear scaling relations (NSR) which connect the nonlinear and linear two point correlation functions. The existence of critical indices for the NSR suggests that the evolution may proceed towards a universal profile which does not change its shape at late times. We begin by clarifying the relation between the density profiles of the individual halo and the slope of the correlation function and discuss the conditions under which the slopes of the correlation function at the extreme nonlinear end can be independent of the initial power spectrum. If the evolution should lead to a profile which preserves the shape at late times, then the correlation function should grow as $a^2$ [in a $\Omega=1$ universe] een at nonlinear scales. We prove that such exact solutions do not exist; however, ther e exists a class of solutions (``psuedo-linear profiles'', PLP's for short) which evolve as $a^2$ to a good approximation. It turns out that the PLP's are the correlation functions which arise if the individual halos are assumed to be isothermal spheres. They are also configurations of mass in which the nonlinear effects of gravitational clustering is a minimum and hence can act as building blocks of the nonlinear universe. We discuss the implicatios of this result.

Large inhomogeneities in neutral hydrogen in the universe can be detected at redshifts $z \leq 10$ using the redshifted 21 cm line emission. We use cosmological N-Body simulations for dark matter and a simple model for baryonic collapse to estimate the signal expected from structures like proto-clusters of galaxies at high redshifts. We study: (i) the standard CDM model, (ii) a modified CDM model with less power at small scales, and (iii) a $\Lambda+$CDM model in a universe with $\Omega_0 + \Omega_\Lambda = 1$. We show that it should be possible for the next generation radio telescopes to detect such structures at the redshift 3.34 with an integration of about 100 hours. We also discuss possible schemes for enhancing signal to noise ratio to detect proto-condensates at high redshifts.

We investigate the transfer of power between different scales and coupling of modes during non-linear evolution of gravitational clustering in an expanding universe. We start with a power spectrum of density fluctuations that is exponentially damped outside a narrow range of scales and use numerical simulations to study evolution of this power spectrum. Non-Linear effects generate power at other scales with most power flowing from larger to smaller scales. The ``cascade'' of power leads to equipartition of energy at smaller scales, implying a power spectrum with index $n\approx -1$. We find that such a spectrum is produced in the range $1 < \delta < 200$ for density contrast $\delta$. This result continues to hold even when small scale power is added to the initial power spectrum. Semi-analytic models for gravitational clustering suggest a tendency for the effective index to move towards a critical index $n_c\approx -1$ in this range. For n<n_c, power in this range grows faster than linear rate, while if n>n_c, it grows at a slower rate - thereby changing the index closer to n_c. At scales larger than the narrow range of scales with initial power, a k^4 tail is produced. We demonstrate that non-linear small scales do not effect the growth of perturbations at larger scales.

The gravitational clustering of collisionless particles in an expanding universe is modelled using some simple physical ideas. I show that it is possible to understand the nonlinear clustering in terms of three well defined regimes: (1) linear regime; (2) quasilinear regime which is dominated by scale-invariant radial infall and (3) nonlinear regime dominated by nonradial motions and mergers. Modelling each of these regimes separately I show how the nonlinear two point correlation function can be related to the linear correlation function in hierarchical models. This analysis leads to results which are in good agreement with numerical simulations thereby providing an explanation for numerical results. Using this model and some simple extensions, it is possible to understand the transfer of power from large to small scales and the behaviour of higher order correlation functions. The ideas presented here will also serve as a powerful analytical tool to investigate nonlinear clustering in different models.

The nonlinear clustering of dark matter particles in an expanding universe is usually studied by N-body simulations. One can gain some insight into this complex problem if simple relations between physical quantities in the linear and nonlinear regimes can be extracted from the results of N-body simulations. Hamilton et al. (1991) and Nityananda and Padmanabhan (1994) have made an attempt in this direction by relating the mean relative pair velocities to the mean correlation function in a useful manner. We investigate this relation and other closely related issues in detail for the case of six different power spectra: power laws with spectral indexes $n=-2,-1$, cold dark matter (CDM), and hot dark matter models with density parameter $\Omega=1$; CDM including a cosmological constant ($\Lambda$) with $\Omega_{CDM}=0.4$, $\Omega_{\Lambda}=0.6$; and $n=-1$ model with $\Omega=0.1$. We find that: (i) Power law spectra lead to self-similar evolution in an $\Omega=1$ universe. (ii) Stable clustering does not hold in an $\Omega=1$ universe to the extent our simulations can ascertain. (iii) Stable clustering is a better approximation in the case of $\Omega<1$ universe in which structure formation freezes out at some low redshift. (iv) The relation between dimensionless pair velocity and the mean correlation function, $\bar\xi$, is only approximately independent of the shape of the power spectrum. At the nonlinear end, the asymptotic value of the dimensionless pair velocity decreases with increasing small scale power, because the stable clustering assumption is not universally true. (v) The relation between the evolved $\bar\xi$ and the linear regime $\bar\xi$ is also not universal but shows a weak spectrum dependence. We present simple theoretical arguments for these conclusions.

Thanks to new technology of observations and fresh inputs from particle physics, cosmology has advanced on both observational and theoretical fronts. It is therefore opportune that we take stock of the cosmological situation today and examine the observational and theoretical constraints as they are now. The bottom line in this review is that despite the availability of the cosmological constant as an extra parameter for flat Friedmann models, the allowed parameter space for such models is very small. The observations that we consider here include the ages of globular clusters, measurement of Hubble's constant, abundance of rich clusters of galaxies, fraction of mass contributed by baryons in rich clusters and abundance of high redshift objects.

The dynamical evolution of collisionless particles in an expanding background is described. After discussing qualitatively the key features, the gravitational clustering of collisionless particles in an expanding universe is modelled using some simple physical ideas. I show that it is indeed possible to understand the nonlinear clustering in terms of three well defined regimes: (1) linear regime (2) quasilinear regime which is dominated by scale-invariant radial infall and (3) nonlinear regime dominated by nonradial motions and mergers. Modelling each of these regimes separately I show how the nonlinear two point correlation function can be related to the linear correlation function in hierarchical models. This analysis leads to results which are in good agreement with numerical simulations thereby providing an explanation for numerical results. The ideas presented here will also serve as a powerful analytical tool to investigate nonlinear clustering in different models. Several implications of the result are discussed.

In an expanding universe, velocity field and gravitational force field are proportional to each other in the linear regime. Neither of these quantities evolve in time and these can be scaled suitably so that the constant of proportionality is unity and velocity and force field are equal. The Zeldovich approximation extends this feature beyond the linear regime, until formation of pancakes. Nonlinear clustering which takes place {\it after} the breakdown of Zeldovich approximation, breaks this relation and the mismatch between these two vectors increases as the evolution proceeds. We suggest that the difference of these two vectors could form the basis for a powerful, new, statistical indicator of nonlinear clustering. We define an indicator called velocity contrast, study its behaviour using N-Body simulations and show that it can be used effectively to delineate the regions where nonlinear clustering has taken place. We discuss several features of this statistical indicator and provide simple analytic models to understand its behaviour. Particles with velocity contrast higher than a threshold have a correlation function which is biased with respect to the original sample. This bias factor is scale dependent and tends to unity at large scales.

The gravitational clustering of collisionless particles in an expanding universe is modelled using some simple physical ideas. I show that it is indeed possible to understand the nonlinear clustering in terms of three well defined regimes: (1) linear regime (2) quasilinear regime which is dominated by scale-invariant radial infall and (3) nonlinear regime dominated by nonradial motions and mergers. Modelling each of these regimes separately I show how the nonlinear two point correlation function can be related to the linear correlation function in heirarchical models. This analysis leads to results which are in good agreement with numerical simulations thereby providing an explanation for numerical results. The ideas presented here will also serve as a powerful anlytical tool to investigate nonlinear clustering in different models. Several implications of the result are discussed.

The constraints on the models for the structure formation arising from various cosmological observations at different length scales are reviewed. The status of different models for structure formation is examined critically in the light of these observations.

We study the evolution of the power spectrum of gravitational potential during the nonlinear clustering in an $\Omega=1$ matter dominated phase. N-body simulations suggest that the potential does not evolve in time even in the quasilinear phase for an $n=-1$ power spectrum. For $n=-2$, the potential evolves in the quasilinear phase but not in the extreme nonlinear phase. Becauase of these facts, the evolution of the gravitational potential in spectra like CDM is less than what would have been expected naively. We discuss a class of CDM like models in which such an interesting conspiracy between the amplitude and local index occurs.

Models for structure formation attempt to predict the power spectrum of density perturbations in the present universe from the initial power spectrum and the nature of dark matter. Observational constraints on the power spectrum at different scales in the present epoch can, therefore, be used to eliminate (or choose between) different theoretical models. Such a comparison is fairly easy at large scales (at which linear theory is valid), and one can use observations like the MBR anisotropy, large scale steaming motions etc to constrain the models. But to discriminate between the models effectivley, it is necessay to constrain the power spectrum at small scales. The most reliable constraints on the power spectra at small scales come from the predicted abundance of bound systems which can be estimated reasonably accurately using Press-Schecter (or similar) methods$^1$. In the past, this method has been used in conjunction with the quasar abundance$^{2-4}$ and cluster abundance$^{5-7}$. We show here that the abundance of damped lyman alpha systems (DLAS, hereafter), provides a far stronger constraint on the models for structure formation. Models with a mixture of hot and cold dark matter $^{8-11}$ (which are consistent with large scale observations) are strongly ruled out by the DLAS constraints while models with cosmological constant $^{12}$ are marginally inconsistent. It is also possible to combine the constraints from the abundance of clusters, DLAS and QSO's to obtain model-independent bounds on the power spectrum at the nonlinear scales. These bounds are to be respected by any viable model for structure formation.

The recent detection of microlensing of stars of LMC by compact objects in the halo of our galaxy suggests that our galaxy is surrounded by a non-luminous halo made of compact objects with mass of about $(0.03-0.5) \msun$. The rate of detection could be consistent with the assumption that these halo objects are distributed with a softened isothermal profile with a core radius of $(2-8)$Kpc and asymptotic circular velocity of $220\kms$. Taken in isolation, this observation is consistent with a universe having only baryonic dark matter (BDM, hereafter) contributing $\Omega_b=\Omega_{total}\simeq0.06$ which will, however, violently contradict several other large scale observations, notably the COBE-DMR results. Alternatively one can assume that galaxies like ours are surrounded by both BDM and non-baryonic dark matter (NBDM, hereafter). A model with a single component for NBDM with, say, $\Omega_b\simeq0.06,\Omega_{cdm}\simeq0.94$, is also ruled out if we demand that: (i) at least thirty percent of the dark matter density within $100$ kpc is baryonic and (ii) galactic structures should have collapsed by redshift of $z=1$. Further if half or more of the dark matter within $100 kpc$ is baryonic then the maximum value of $\Omega_{dm}$, contributed by NBDM clustered at galactic scales, can only be about $\Omega_{max}\simeq0.29$. Models with C+HDM cannot satisfy this constraint but $\Lambda$+ CDM models are still

: From the epoch of recombination $(z\approx 10^3)$ till today, the typical density contrasts have grown by a factor of about $10^6$ in a Friedmann universe with $\Omega=1$. However, during the same epoch the typical gravitational potential has grown only by a factor of order unity. We present theoretical arguments explaining the origin of this approximate constancy of gravitational potential. This fact can be exploited to provide a new, powerful, approximation scheme to study the formation of nonlinear structures in the universe. The essential idea of this method is to evolve the initial distribution of particles using a gravitational potential frozen in time (Frozen Potential Approximation). This approximation provides valuable insight into understanding various features of nonlinear evolution; for example, it provides a simple explanation as to why pancakes remain thin during the evolution even in the absence of any artificial, adhesion-like, damping terms. We compare the trajectories of particles in various approximations. We also discuss a few applications of the frozen potential approximation.