This article investigates the stability of a generic Kasner spacetime to linear perturbations, both at late and early times. It demonstrates that the perturbation of the Weyl tensor diverges at late time in all cases but in the particular one in which the Kasner spacetime is the product of a two-dimensional Milne spacetime and a two-dimensional Euclidean space. At early times, the perturbation of the Weyl tensor also diverges unless one imposes a condition on the perturbations so as to avoid the most divergent modes to be excited.

Bulk dark energy properties are determined by the redshift evolution of its pressure-to-density ratio, $w_{de}(z)$. An experimental goal is to decide if the dark energy is dynamical, as in the quintessence (and phantom) models treated here. We show that a three-parameter approximation $w_{de}(z; \epsilon_s, \epsilon_{\phi\infty}, \zeta_s)$ fits well the ensemble of trajectories for a wide class of late-inflaton potentials $V(\phi)$. Markov Chain Monte Carlo probability calculations are used to confront our $w_{de}(z)$ trajectories with current observational information on Type Ia supernova, Cosmic Microwave Background, galaxy power spectra, weak lensing and the Lyman-${\alpha}$ forest. We find the best constrained parameter is a low redshift slope parameter, $\epsilon_s \propto (\partial \ln V / \partial \phi)^2$ when the dark energy and matter have equal energy densities. A tracking parameter $\epsilon_{\phi\infty}$ defining the high-redshift attractor of $1+w_{de}$ is marginally constrained. Poorly determined is $\zeta_s$, characterizing the evolution of $\epsilon_s$, and a measure of $\partial^2 \ln V / \partial \phi^2$ . The constraints we find already rule out some popular quintessence and phantom models, or restrict their potential parameters. We also forecast how the next generation of cosmological observations improve the constraints: by a factor of about five on $\epsilon_s$ and $\epsilon_{\phi\infty}$, but with $\zeta_s$ remaining unconstrained (unless the true model significantly deviates from $\Lambda$CDM). Thus potential reconstruction beyond an overall height and a gradient is not feasible for the large space of late-inflaton models considered here.

We study in detail (p)reheating after multi-field inflation models with a particular focus on N-flation. We consider a variety of different couplings between the inflatons and the matter sector, including both quartic and trilinear interactions with a light scalar field. We show that the presence of multiple oscillating inflatons makes parametric resonance inefficient in the case of the quartic interactions. Moreover, perturbative processes do not permit a complete decay of the inflaton for this coupling. In order to recover the hot big bang, we must instead consider trilinear couplings. In this case we show that strong nonperturbative preheating is possible via multi-field tachyonic resonance. In addition, late-time perturbative effects do permit a complete decay of the condensate. We also study the production of gauge fields for several prototype couplings, finding similar results to the trilinear scalar coupling. During the course of our analysis we develop the mathematical theory of the quasi-periodic Mathieu equation, the multi-field generalization of the Floquet theory familiar from preheating after single field inflation. We also elaborate on the theory of perturbative decays of a classical inflaton condensate, which is applicable in single-field models also.

We study (p)reheating in modular (closed string) inflationary scenarios, with a special emphasis on Kahler moduli/Roulette models. It is usually assumed that reheating in such models occurs through perturbative decays. However, we find that there are very strong non-perturbative preheating decay channels related to the particular shape of the inflaton potential (which is highly nonlinear and has a very steep minimum). Preheating after modular inflation, proceeding through a combination of tachyonic instability and broad-band parametric resonance, is perhaps the most violent example of preheating after inflation known in the literature. Further, we consider the subsequent transfer of energy to the standard model sector in scenarios where the standard model particles are confined to a D7-brane wrapping the inflationary blow-up cycle of the compactification manifold or, more interestingly, a non-inflationary blow up cycle. We explicitly identify the decay channels of the inflaton in these two scenarios. We also consider the case where the inflationary cycle shrinks to the string scale at the end of inflation; here a field theoretical treatment of reheating is insufficient and one must turn instead to a stringy description. We estimate the decay rate of the inflaton and the reheat temperature for various scenarios.

We propose a qualitatively new mechanism for generating cosmological fluctuations from inflation. The non-equilibrium excitation of interacting scalar fields often evolves into infra-red (IR) and ultra-violet (UV) cascading, resulting in an intermediate scaling regime. We observe elements of this phenomenon in a simple model with inflaton \phi and iso-inflaton \chi fields interacting during inflation via the coupling g^2 (\phi-\phi_0)^2 \chi^2. Iso-inflaton particles are created during inflation when they become instantaneously massless at \phi=\phi_0, with occupation numbers not exceeding unity. We point out that very quickly the produced \chi particles become heavy and their multiple re-scatterings off the homogeneous condensate \phi(t) generates bremschtrahlung radiation of light inflaton IR fluctuations with high occupation numbers. The subsequent evolution of these IR fluctuations is qualitatively similar to that of the usual inflationary fluctuations, but their initial amplitude is different. The IR cascading generates a bump-shaped contribution to the cosmological curvature fluctuations, which can even dominate over the usual fluctuations for g^2>0.06. The IR cascading curvature fluctuations are significantly non-gaussian and the strength and location of the bump are model-dependent, through g^2 and \phi_0. The effect from IR cascading fluctuations is significantly larger than that from the momentary slowing-down of \phi(t). With a sequence of such bursts of particle production, the superposition of the bumps can lead to a new broad band non-gaussian component of cosmological fluctuations added to the usual fluctuations. Such a sequence of particle creation events can, but need not, lead to trapped inflation.

A new and powerful probe of the origin and evolution of structures in the Universe has emerged and been actively developed over the last decade. In the coming decade, non-Gaussianity, i.e., the study of non-Gaussian contributions to the correlations of cosmological fluctuations, will become an important probe of both the early and the late Universe. Specifically, it will play a leading role in furthering our understanding of two fundamental aspects of cosmology and astrophysics: (i) the physics of the very early universe that created the primordial seeds for large-scale structures, and (ii) the subsequent growth of structures via gravitational instability and gas physics at later times. To date, observations of fluctuations in the Cosmic Microwave Background (CMB) and the Large-Scale Structure of the Universe (LSS) have focused largely on the Gaussian contribution as measured by the two-point correlations (or the power spectrum) of density fluctuations. However, an even greater amount of information is contained in non-Gaussianity and a large discovery space therefore still remains to be explored. Many observational probes can be used to measure non-Gaussianity, including CMB, LSS, gravitational lensing, Lyman-alpha forest, 21-cm fluctuations, and the abundance of rare objects such as clusters of galaxies and high-redshift galaxies. Not only does the study of non-Gaussianity maximize the science return from a plethora of present and future cosmological experiments and observations, but it also carries great potential for important discoveries in the coming decade.

A new class of non-Gaussian curvature fluctuations \zeta_{pr} (\bx) \equiv \delta N(\chi_i) arises from the post-inflation preheating behaviour of a non-inflaton field \chi_i. Its billiard-like chaotic dynamics imprints regular log-spaced narrow spikes in the number of preheating e-folds N(\chi_i). We perform highly accurate lattice simulations of SUSY-inspired quartic inflaton and coupling potentials in a separate-universe approximation to compute N(\chi_i) as a function of the (nearly homogeneous) initial condition \chi_i. The super-horizon modes of \chi_i(\bx) result in positive spiky excursions in \zeta_{pr} and hence negative gravitational potential fluctuations added to the usual sign-independent inflaton-induced perturbations, observably manifested in large cosmic structures and as (polarized) temperature CMB cold spots.

We study the stochastic background of gravitational waves produced from preheating in hybrid inflation models. We investigate different dynamical regimes of preheating in these models and we compute the resulting gravity wave spectra using analytical estimates and numerical simulations. We discuss the dependence of the gravity wave frequencies and amplitudes on the various potential parameters. We find that large regions of the parameter space leads to gravity waves that may be observable in upcoming interferometric experiments, including Advanced LIGO, but this generally requires very small coupling constants.

We show that expanding or contracting Kasner universes are unstable due to the amplification of gravitational waves (GW). As an application of this general relativity effect, we consider a pre-inflationary anisotropic geometry characterized by a Kasner-like expansion, which is driven dynamically towards inflation by a scalar field. We investigate the evolution of linear metric fluctuations around this background, and calculate the amplification of the long-wavelength GW of a certain polarization during the anisotropic expansion (this effect is absent for another GW polarization, and for scalar fluctuations). These GW are superimposed to the usual tensor modes of quantum origin from inflation, and are potentially observable if the total number of inflationary e-folds exceeds the minimum required to homogenize the observable universe only by a small margin. Their contribution to the temperature anisotropy angular power spectrum decreases with the multipole l as l^(-p), where p depends on the slope of the initial GW power-spectrum. Constraints on the long-wavelength GW can be translated into limits on the total duration of inflation and the initial GW amplitude. The instability of classical GW (and zero-vacuum fluctuations of gravitons) during Kasner-like expansion (or contraction) may have other interesting applications. In particular, if GW become non-linear, they can significantly alter the geometry before the onset of inflation.

The presence of Kaluza-Klein particles in the universe is a potential manifestation of string theory cosmology. In general, they can be present in the high temperature bath of the early universe. In particular examples, string theory inflation often ends with brane-antibrane annihilation followed by the energy cascading through massive closed string loops to KK modes which then decay into lighter standard model particles. However, massive KK modes in the early universe may become dangerous cosmological relics if the inner manifold contains warped throat(s) with approximate isometries. In the complimentary picture, in the AdS/CFT dual gauge theory with extra symmetries, massive glueballs of various spins become the dangerous cosmological relics. The decay of these angular KK modes/glueballs, located around the tip of the throat, is caused by isometry breaking which results from gluing the throat to the compact CY manifold. We address the problem of these angular KK particles/glueballs, studying their interactions and decay channels, from the theory side, and the resulting cosmological constraints on the warped compactification parameters, from the phenomenology side. The abundance and decay time of the long-lived non-relativistic angular KK modes depend strongly on the parameters of the warped geometry, so that observational constraints rule out a significant fraction of the parameter space. In particular, the coupling of the angular KK particles can be weaker than gravitational.

String theory compactification involves manifolds with multiple warp factors. For cosmological applications, one often introduces a short, high-energy inflationary throat, and a long,low-energy Standard Model throat. It is assumed that at the end of inflation,the excited Kaluza-Klein modes from the Inflationary throat tunnel to the SM throat and reheat Standard Model degrees of freedom, which are attached to probe brane(s). However, the huge hierarchy of energy scales can result in a highly dynamic transition of the throat geometry. We point out that in such a cosmological scenario the Standard Model throat (together with SM brane) will be cloaked by a Schwarzschild horizon, produced by the Kaluza-Klein modes tunneling from the short throat. The Black Brane formation is dual to the first order chiral phase transition of the cascading gauge theory. We calculate the critical energy density corresponding the formation of the the BH horizon in the long throat. We discuss the duality between "Black Universe" cosmology and an expanding universe driven by the hot gauge theory radiation. We address the new problem of the hierarchical multiple-throat scenarios: SM brane disappearance after the decay of the BH horizon.

Usual inflation is realized with a slow rolling scalar field minimally coupled to gravity. In contrast, we consider dynamics of a scalar with a flat effective potential, conformally coupled to gravity. Surprisingly, it contains an attractor inflationary solution with the rapidly rolling inflaton field. We discuss models with the conformal inflaton with a flat potential (including hybrid inflation). There is no generation of cosmological fluctuations from the conformally coupled inflaton. We consider realizations of modulated (inhomogeneous reheating) or curvaton cosmological fluctuations in these models. We also implement these unusual features for the popular string-theoretic warped inflationary scenario, based on the interacting D3-anti D3 branes. The original warped brane inflation suffers a large inflaton mass due to conformal coupling to 4-dimensional gravity. Instead of considering this as a problem and trying to cure it with extra engineering, we show that warped inflation with the conformally coupled, rapidly rolling inflaton is yet possible with N=37 efoldings, which requires low energy scales 1-100 TeV of inflation. Coincidentally, the same warping numerology can be responsible for the hierarchy. It is shown that the scalars associated with angular isometries of the warped geometry of compact manifold (e.g. S^3 of KS geometry) have solutions identical to conformally coupled modes and also cannot be responsible for cosmological fluctuations. We discuss other possibilities.

Preheating after inflation involves large, time-dependent field inhomogeneities, which act as a classical source of gravitational radiation. The resulting spectrum might be probed by direct detection experiments if inflation occurs at a low enough energy scale. In this paper, we develop a theory and algorithm to calculate, analytically and numerically, the spectrum of energy density in gravitational waves produced from an inhomogeneous background of stochastic scalar fields in an expanding universe. We derive some generic analytical results for the emission of gravity waves by stochastic media of random fields, which can test the validity/accuracy of numerical calculations. We contrast our method with other numerical methods in the literature, and then we apply it to preheating after chaotic inflation. In this case, we are able to check analytically our numerical results, which differ significantly from previous works. We discuss how the gravity wave spectrum builds up with time and find that the amplitude and the frequency of its peak depend in a relatively simple way on the characteristic spatial scale amplified during preheating. We then estimate the peak frequency and amplitude of the spectrum produced in two models of preheating after hybrid inflation, which for some parameters may be relevant for gravity wave interferometric experiments.

We study 2-field inflation models based on the ``large-volume'' flux compactification of type IIB string theory. The role of the inflaton is played by a K\"ahler modulus \tau corresponding to a 4-cycle volume and its axionic partner \theta. The freedom associated with the choice of Calabi Yau manifold and the non-perturbative effects defining the potential V(\tau, \theta) and kinetic parameters of the moduli bring an unavoidable statistical element to theory prior probabilities within the low energy landscape. The further randomness of (\tau, \theta) initial conditions allows for a large ensemble of trajectories. Features in the ensemble of histories include ``roulette tractories'', with long-lasting inflations in the direction of the rolling axion, enhanced in number of e-foldings over those restricted to lie in the \tau-trough. Asymptotic flatness of the potential makes possible an eternal stochastic self-reproducing inflation. A wide variety of potentials and inflaton trajectories agree with the cosmic microwave background and large scale structure data. In particular, the observed scalar tilt with weak or no running can be achieved in spite of a nearly critical de Sitter deceleration parameter and consequently a low gravity wave power relative to the scalar curvature power.

Recently a non-inflationary mechanism of generation of scale-free cosmological perturbations of metric was proposed by Brandenberger, Nayeri, and Vafa in the context of the string gas cosmology. We discuss various problems of their model and argue that the cosmological perturbations of metric produced in this model have blue spectrum with a spectral index n = 5, which strongly disagrees with observations. We conclude that this model in its present form is not a viable alternative to inflationary cosmology.

We consider the nonlinear dynamics of inflaton fragmentation during and after preheating in the simplest model of chaotic inflation. While the earlier regime of parametric resonant particle production and the later turbulent regime of interacting fields evolving towards equilibrium are well identified and understood, the short intermediate stage of violent nonlinear dynamics remains less explored. Lattice simulations of fully nonlinear preheating dynamics show specific features of this intermediate stage: occupation numbers of the scalar particles are peaked, scalar fields become significantly non-gaussian and the field dynamics become chaotic and irreversible. Visualization of the field dynamics in configuration space reveals that nonlinear interactions generate non-gaussian inflaton inhomogeneities with very fast growing amplitudes. The peaks of the inflaton inhomogeneities coincide with the peaks of the scalar field(s) produced by parametric resonance. When the inflaton peaks reach their maxima, they stop growing and begin to expand. The subsequent dynamics is determined by expansion and superposition of the scalar waves originating from the peaks. Multiple wave superposition results in phase mixing and turbulent wave dynamics. Thus, the short intermediate stage is defined by the formation, expansion and collision of bubble-like field inhomogeneities associated with the peaks of the original gaussian field. This process is qualitatively similar to the bubble-like inflaton fragmentation that occurs during tachyonic preheating after hybrid or new inflation.

We investigate the effects of bosonic trilinear interactions in preheating after chaotic inflation. A trilinear interaction term allows for the complete decay of the massive inflaton particles, which is necessary for the transition to radiation domination. We found that typically the trilinear term is subdominant during early stages of preheating, but it actually amplifies parametric resonance driven by the four-legs interaction. In cases where the trilinear term does dominate during preheating, the process occurs through periodic tachyonic amplifications with resonance effects, which is so effective that preheating completes within a few inflaton oscillations. We develop an analytic theory of this process, which we call tachyonic resonance. We also study numerically the influence of trilinear interactions on the dynamics after preheating. The trilinear term eventually comes to dominate after preheating, leading to faster rescattering and thermalization than could occur without it. Finally, we investigate the role of non-renormalizable interaction terms during preheating. We find that if they are present they generally dominate (while still in a controllable regime) in chaotic inflation models. Preheating due to these terms proceeds through a modified form of tachyonic resonance.

In String theory realizations of inflation, the end point of inflation is often brane-anti brane annihilation. We consider the processes of reheating of the Standard Model universe after brane inflation. We identify the channels of inflaton energy decay, cascading from tachyon annihilation through massive closed string loops, KK modes, and brane displacement moduli to the lighter standard model particles. Cosmological data constrains scenarios by putting stringent limits on the fraction of reheating energy deposited in gravitons and nonstandard sector massive relics. We estimate the energy deposited into various light degrees of freedom in the open and closed string sectors, the timing of reheating, and the reheating temperature. Production of gravitons is significantly suppressed in warped inflation. However, we predict a residual gravitational radiation background at the level $\Omega_{GW} \sim 10^{-8}$ of the present cosmological energy density. We also extend our analysis to multiple throat scenarios. A viable reheating would be possible in a single throat or in a certain subclass of multiple throat scenarios of the KKLMMT type inflation model, but overproduction of massive KK modes poses a serious problem. The problem is quite severe if some inner manifold comes with approximate isometries (angular KK modes) or if there exists a throat of modest length other than the standard model throat, possibly associated with some hidden sector (low-lying KK modes).

We study the out-of-equilibrium nonlinear dynamics of fields after post-inflationary preheating. During preheating, the energy in the homogeneous inflaton is exponentially rapidly transfered into highly occupied out-of-equilibrium inhomogeneous modes, which subsequently evolve towards equilibrium. The infrared modes excited during preheating evolve towards a saturated distribution long before thermalization completes. We compute the equation of state during and immediately after preheating. It rapidly evolves towards radiation domination long before the actual thermal equilibrium is established. The exact time of this transition is a non-monotonic function of the coupling between the inflaton and the decay products, and it varies only very weakly (around 10^(-35) s) as this coupling changes over several orders of magnitude. This result is applied to refine the relation between the number of efoldings N and the physical wavelength of perturbations generated during inflation. We also discuss the implications for the theory of modulated perturbations from preheating. We finally argue that many questions of the thermal history of the universe should be addressed in terms of pre-thermalization, illustrating this point with a calculation of perturbative production of gravitinos immediately after chaotic inflation. We also highlight the effects of three-legs inflaton interactions on the dynamics of preheating and thermalization in an expanding universe.

We study classical dynamics of a probe Dp-brane moving in a background sourced by a stack of Dp-branes. In this context the physics is similar to that of the effective action for open-string tachyon condensation, but with a power-law runaway potential. We show that small inhomogeneous ripples of the probe brane embedding grow with time, leading to folding of the brane as it moves. We give a full nonlinear analytical treatment of inhomogeneous brane dynamics, suitable for the Dirac-Born-Infeld + Wess-Zumino theory with arbitrary runaway potential, in the case where the source branes are BPS. In the near-horizon geometry, the inhomogeneous brane motion has a dual description in terms of free streaming of massive relativistic test particles originating from the initial hypersurface of the probe brane. We discuss limitations of the effective action description around loci of self-crossing of the probe brane (caustics). We also discuss the effect of brane folding in application to the theory of cosmological fluctuations in string theory inflation.

We introduce the general N=1 gauge theory superconformally coupled to supergravity. The theory has local SU(2,2|1) symmetry and no dimensional parameters. The superconformal origin of the Fayet-Iliopoulos terms is clarified. The phase of this theory with spontaneously broken conformal symmetry gives various formulations of N=1 supergravity interacting with matter, depending on the choice of the R-symmetry fixing. We have found that the locally superconformal theory is useful for describing the physics of the early universe with a conformally flat FRW metric. Few applications of superconformal theory to cosmology include the study of i) particle production after inflation, particularly the non-conformal helicity 1/2 states of gravitino, ii) the super-Higgs effect in cosmology and the derivation of the equations for the gravitino interacting with any number of chiral and vector multiplets in the gravitational background with varying scalar fields, iii) the weak coupling limit of supergravity and gravitino-goldstino equivalence. This explains why gravitino production in the early universe is not suppressed in the limit of weak gravitational coupling. We discuss the possible existence of an unbroken phase of the superconformal theories, interpreted as a strong coupling limit of supergravity.

We study quantum effects on moduli dynamics arising from the production of particles which are light at special points in moduli space. The resulting forces trap the moduli at these points, which often exhibit enhanced symmetry. Moduli trapping occurs in time-dependent quantum field theory, as well as in systems of moving D-branes, where it leads the branes to combine into stacks. Trapping also occurs in an expanding universe, though the range over which the moduli can roll is limited by Hubble friction. We observe that a scalar field trapped on a steep potential can induce a stage of acceleration of the universe, which we call trapped inflation. Moduli trapping ameliorates the cosmological moduli problem and may affect vacuum selection. In particular, rolling moduli are most powerfully attracted to the points with the largest number of light particles, which are often the points of greatest symmetry. Given suitable assumptions about the dynamics of the very early universe, this effect might help to explain why among the plethora of possible vacuum states of string theory, we appear to live in one with a large number of light particles and (spontaneously broken) symmetries. In other words, some of the surprising properties of our world might arise not through pure chance or miraculous cancellations, but through a natural selection mechanism during dynamical evolution.

We describe an algorithm and a C++ implementation that we have written and made available for calculating the fully nonlinear evolution of 5D braneworld models with scalar fields. Bulk fields allow for the stabilization of the extra space. However, they complicate the dynamics of the system, so that analytic calculations (performed within an effective 4D theory) are typically only reliable close to stabilized configurations or when the evolution of the extra space is negligible. In the general case, a numerical study of the 5D equations is necessary, and the algorithm and code we describe are the first ones designed for this task. The program and its full documentation are available on the Web at http://www.cita.utoronto.ca/~jmartin/BRANECODE/. In this paper we provide a brief overview of what the program does and how to use it.

We consider warped compactifications in (4+d)-dimensional theories, with four dimensional de Sitter dS_4 vacua (with Hubble parameter H) and with a compact internal space. After introducing a gauge-invariant formalism for the generic metric perturbations of these backgrounds, we focus on modes which are scalar with respect to dS_4. The physical eigenmasses of these modes acquire a large universal tachyonic contribution -12d/(d+2) H^2, independently of the stabilization mechanism for the compact space, in addition to the usual KK masses, which instead encode the effects of the stabilization. General arguments, as well as specific examples, lead us to conjecture that, for sufficiently large dS curvature, the compactified geometry becomes gravitationally unstable due to the tachyonic growth of the scalar perturbations. This mean that for any stabilization mechanism the curvature of the dS geometry cannot exceed some critical value. We relate this effect to the anisotropy of the bulk geometry and suggest the end points of the instability. Of relevance for inflationary cosmology, the perturbations of the bulk metric inevitably induce a new modulus field, which describes the conformal fluctuations of the 4 dimensional metric. If this mode is light during inflation, the induced conformal fluctuations will be amplified with a scale free spectrum and with an amplitude which is disentangled from the standard result of slow-roll inflation. The conformal 4d metric fluctuations give rise to a very generic realization of the mechanism of modulated cosmological fluctuations, related to spatial variation of couplings during (p)reheating after inflation.

Inflation universally produces classical almost scale free Gaussian inhomogeneities of any light scalars. Assuming the coupling constants at the time of inflation depend on some light moduli fields, we encounter the generation of modulated cosmological fluctuations from (p)reheating. This is an alternative mechanism to generate observable (almost) scale free adiabatic metric perturbations. We extend this idea to the class of hybrid inflation, where the bifurcation value of the inflaton is modulated by the spatial inhomogeneities of the couplings. As a result, the symmetry breaking after inflation occurs not simultaneously in space but with the time laps in different Hubble patches inherited from the long-wavelength moduli inhomogeneities. To calculate modulated fluctuations we introduce techniques of general relativistic matching conditions for metric perturbations at the time hypersurface where the equation of state after inflation undergoes a jump, without evoking the detailed microscopic physics, as far as it justifies the jump. We apply this theory to the modulated fluctuations from the hybrid and chaotic inflations. We discuss what distinguish the modulated from the inflation-driven fluctuations, in particular, their spectral index, modification of the consistency relation and the issue of weak non-Gaussianity.

Dirac-Born-Infeld type effective actions reproduce many aspects of string theory classical tachyon dynamics of unstable Dp-branes. The inhomogeneous tachyon field rolling from the top of its potential forms topological defects of lower codimensions. In between them, as we show, the tachyon energy density fragments into a p-dimensional web-like high density network evolving with time. We present an analytic asymptotic series solution of the non-linear equations for the inhomogeneous tachyon and its stress energy. The generic solution for a tachyon field with a runaway potential in arbitrary dimensions is described by the free streaming of noninteracting massive particles whose initial velocities are defined by the gradients of the initial tachyon profile. Thus, relativistic particle mechanics is a dual picture of the tachyon field effective action. Implications of this picture for inflationary models with a decaying tachyon field are discussed.

Braneworld models in extra dimensions can be tested in laboratory by the coupling of the radion to the Standard Model fields. The identification of the radion as a canonically normalized field involves a careful General Relativity treatment: if a bulk scalar is responsible for the stabilization of the system, its fluctuations are entangled with the perturbations of the metric and they also have to be taken into account (similarly to the well-developed theory of scalar metric perturbations in 4D cosmology with a scalar field). Extracting a proper dynamical variable in a warped geometry/scalar setting is a nontrivial task, performed so far only in the limit of negligible backreaction of the scalar field on the background geometry. We perform the general calculation, diagonalizing the action up to second order in the perturbations and identifying the physical eigenmodes of the system for any amplitude of the bulk scalar. This computation allows us to derive a very simple expression for the exact coupling of the eigenmodes to the Standard Model fields on the brane, valid for an arbitrary background configuration. As an application, we discuss the Goldberger-Wise mechanism for the stabilization of the radion in the Randall-Sundrum type models. The existing studies, limited to small amplitude of the bulk scalar field, are characterized by a radion mass which is significantly below the physical scale at the observable brane. We extend them beyond the small backreaction regime. For intermediate amplitudes, the radion mass approaches the electroweak scale, while its coupling to the observable brane remains nearly constant. At very high amplitudes, the radion mass instead decreases, while the coupling sharply increases. Severe experimental constraints are expected in this regime.

We give a full nonlinear numerical treatment of time-dependent 5d braneworld geometry, which is determined self-consistently by potentials for the scalar field in the bulk and at two orbifold branes, supplemented by boundary conditions at the branes. We describe the BraneCode, an algorithm which we designed to solve the dynamical equations numerically. We applied the BraneCode to braneworld models and found several novel phenomena of the brane dynamics. Starting with static warped geometry with de Sitter branes, we found numerically that this configuration is often unstable due to a tachyonic mass of the radion during inflation. If the model admits other static configurations with lower values of de Sitter curvature, this effect causes a violent re-structuring towards them, flattening the branes, which appears as a lowering of the 4d effective cosmological constant. Braneworld dynamics can often lead to brane collisions. We found that in the presence of the bulk scalar field, the 5d geometry between colliding branes approaches a universal, homogeneous, anisotropic strong gravity Kasner-like asymptotic, irrespective of the bulk/brane potentials. The Kasner indices of the brane directions are equal to each other but different from that of the extra dimension.

We investigate scalar perturbations from inflation in braneworld cosmologies with extra dimensions. For this we calculate scalar metric fluctuations around five dimensional warped geometry with four dimensional de Sitter slices. The background metric is determined self-consistently by the (arbitrary) bulk scalar field potential, supplemented by the boundary conditions at both orbifold branes. Assuming that the inflating branes are stabilized (by the brane scalar field potentials), we estimate the lowest eigenvalue of the scalar fluctuations - the radion mass. In the limit of flat branes, we reproduce well known estimates of the positive radion mass for stabilized branes. Surprisingly, however, we found that for de Sitter (inflating) branes the square of the radion mass is typically negative, which leads to a strong tachyonic instability. Thus, parameters of stabilized inflating braneworlds must be constrained to avoid this tachyonic instability. Instability of "stabilized" de Sitter branes is confirmed by the BraneCode numerical calculations in the accompanying paper hep-th/0309001. If the model's parameters are such that the radion mass is smaller than the Hubble parameter, we encounter a new mechanism of generation of primordial scalar fluctuations, which have a scale free spectrum and acceptable amplitude.

The Cosmic Microwave Background (CMB) anisotropy power on the largest angular scales observed both by WMAP and COBE DMR appears to be lower than the one predicted by the standard model of cosmology with almost scale free primordial perturbations arising from a period of inflation \cite{cobe,Bennett:2003bz,Spergel,Peiris}. One can either interpret this as a manifestation of cosmic variance or as a physical effect that requires an explanation. We discuss various mechanisms that could be responsible for the suppression of such low $\ell$ multipoles. Features in the late time evolution of metric fluctuations may do this via the integral Sachs-Wolfe effect. Another possibility is a suppression of power at large scales in the primordial spectrum induced by a fast rolling stage in the evolution of the inflaton field at the beginning of the last 65 e-folds of inflation. We illustrate this effect in a simple model of inflation and fit the resulting CMB spectrum to the observed temperature-temperature (TT) power spectrum. We find that the WMAP observations suggest a cutoff at $k_c=4.9^{+1.3}_{-1.6}\times 10^{-4}$Mpc$^{-1}$ at 68% confidence, while only an upper limit of $k_c < 7.4\times 10^{-4}$Mpc$^{-1}$ at 95%. Thus, although it improves the fit of the data, the presence of a cutoff in power spectrum is only required at a level close to $2\sigma$. This is obtained with a prior which corresponds to equal distribution wrt $k_c$. We discuss how other choices (such as an equal distribution wrt $\ln k_c$ which is natural in the context of inflation) can affect the statistical interpretation.

Superstring theory, models with extra dimensions and other SUSY models generically predict that the coupling constants are in fact vacuum expectation values of fields like the dilaton, moduli etc. Assuming some of these fields are light during inflation, we get generation of small classical inhomogeneities in these fields from inflation. Consequently, coupling constants inherit small inhomogeneities at scales much larger than the causal horizon in the early universe. After the moduli get pinned down to their minima, the spatial variations of coupling constants in the late time universe will be erased. However, inhomogeneities in coupling constants in the very early universe would generate modulated large scale fluctuations in all relic species that are produced due to interactions and freezing out. Moreover (p)reheating of the inflaton field results in modulated curvature fluctuations. Even if the standard inflaton fluctuations are suppressed, in this picture we may have pure curvature cosmological fluctuations entirely generated by the modulated spatial variations of the coupling constants during preheating.

We consider the quasi-de Sitter geometry of the inflationary universe. We calculate the energy flux of the slowly rolling background scalar field through the quasi-de Sitter apparent horizon and set it equal to the change of the entropy (1/4 of the area) multiplied by the temperature, dE=TdS. Remarkably, this thermodynamic law reproduces the Friedmann equation for the rolling scalar field. The flux of the slowly rolling field through the horizon of the quasi-de Sitter geometry is similar to the accretion of a rolling scalar field onto a black hole, which we also analyze. Next we add inflaton fluctuations which generate scalar metric perturbations. Metric perturbations result in a variation of the area entropy. Again, the equation dE=TdS with fluctuations reproduces the linearized Einstein equations. In this picture as long as the Einstein equations hold, holography does not put limits on the quantum field theory during inflation. Due to the accumulating metric perturbations, the horizon area during inflation randomly wiggles with dispersion increasing with time. We discuss this in connection with the stochastic decsription of inflation. We also address the issue of the instability of inflaton fluctuations in the ``hot tin can'' picture of de Sitter horizon.

We investigate the generation of primordial gravitational waves from inflation in braneworld cosmologies with extra dimensions. Advantage of using primordial gravitational waves to probe extra dimensions is that their theory depends only on the geometry, not on the microscopic models of inflation and stabilization. D(D-3)/2 degrees of freedom of the free bulk gravitons are projected onto the 3d brane as tensor, vector and scalar modes. We found the following no-go results for a generic geometry of a five (or D) dimensional warped metric with four dimensional de Sitter (inflationary) slices and two (or one) edge of the world branes: Massive KK graviton modes are not generated from inflation (with the Hubble parameter H) due to the gap in the KK spectrum; the universal lower bound on the gap is sqrt{3/2} H. Massless scalar and vector projections of the bulk gravitons are absent, unlike in geometries with KK compactification. A massless 4d tensor mode is generated from inflation with the amplitude H/M_P, where M_P is the effective Planck mass during inflation, derived from the D dimensional fundamental mass M_S and the volume of the inner dimensions. However, M_P for a curved dS braneworld may differ from that of the flat brane at low energies, either due to the H-dependence of the inner space volume or variations in the brane separation before stabilization. Thus the amplitude of gravitational waves from inflation in braneworld cosmology may be different from that predicted by inflation in 4d theory.

Recently Hollands and Wald argued that inflation does not solve any of the major cosmological problems. We explain why we disagree with their arguments. They also proposed a new speculative mechanism of generation of density perturbations. We show that in their scenario the inhomogeneities responsible for the large scale structure observed today were generated at an epoch when the energy density of the hot universe was 10^{95} times greater than the Planck density. The only way to avoid this problem is to assume that there was a stage of inflation in the early universe.

We consider scalar Born-Infeld type theories with arbitrary potentials V(T) of a scalar field T. We find that for models with runaway potentials V(T) the generic inhomogeneous solutions after a short transient stage can be very well approximated by the solutions of a Hamilton-Jacobi equation that describes free streaming wave front propagation. The analytic solution for this wave propagation shows the formation of caustics with multi-valued regions beyond them. We verified that these caustics appear in numerical solutions of the original scalar BI non-linear equations. Our results include the scalar BI model with an exponential potential, which was recently proposed as an effective action for the string theory tachyon in the approximation where high-order spacetime derivatives of T are truncated. Since the actual string tachyon dynamics contain derivatives of all orders, the tachyon BI model with an exponential potential becomes inadequate when the caustics develop because high order spatial derivatives of T become divergent. BI type tachyon theory with a potential decreasing at large T could have interesting cosmological applications because the tachyon field rolling towards its ground state at infinity acts as pressureless dark matter. We find that inhomogeneous cosmological tachyon fluctuations rapidly grow and develop multiple caustics. Any considerations of the role of the tachyon field in cosmology will have to involve finding a way to predict the behavior of the field at and beyond these caustics.

We consider cosmological consequences of string theory tachyon condensation. We show that it is very difficult to obtain inflation in the simplest versions of this theory. Typically, inflation in these theories could occur only at super-Planckian densities, where the effective 4D field theory is inapplicable. Reheating and creation of matter in models where the tachyon potential V(T) has a minimum at infinitely large T is problematic because the tachyon field in such theories does not oscillate. If the universe after inflation is dominated by the energy density of the tachyon condensate, it will always remain dominated by the tachyons. It might happen that string condensation is responsible for a short stage of inflation at a nearly Planckian density, but one would need to have a second stage of inflation after that. This would imply that the tachyon played no role in the post-inflationary universe until the very late stages of its evolution. These problems do not appear in the recently proposed models of hybrid inflation where the complex tachyon field has a minimum at T << M_p.

We study the dynamical equations for extra-dimensional dependence of a warp factor and a bulk scalar in 5d brane world scenarios with induced brane metric of constant curvature. These equations are similar to those for the time dependence of the scale factor and a scalar field in 4d cosmology, but with the sign of the scalar field potential reversed. Based on this analogy, we introduce novel methods for studying the warped geometry. We construct the full phase portraits of the warp factor/scalar system for several examples of the bulk potential. This allows us to view the global properties of the warped geometry. For flat branes, the phase portrait is two dimensional. Moving along typical phase trajectories, the warp factor is initially increasing and finally decreasing. All trajectories have timelike gradient-dominated singularities at one or both of their ends, which are reachable in a finite distance and must be screened by the branes. For curved branes, the phase portrait is three dimensional. However, as the warp factor increases the phase trajectories tend towards the two dimensional surface corresponding to flat branes. We discuss this property as a mechanism that may stretch the curved brane to be almost flat, with a small cosmological constant. Finally, we describe the embedding of branes in the 5d bulk using the phase space geometric methods developed here. In this language the boundary conditions at the branes can be described as a 1d curve in the phase space. We discuss the naturalness of tuning the brane potential to stabilize the brane world system.

We consider the evolution of FRW cosmological models and linear perturbations of tachyon matter rolling towards a minimum of its potential. The tachyon coupled to gravity is described by an effective 4d field theory of string theory tachyon. In the model where a tachyon potential $V(T)$ has a quadratic minimum at finite value of the tachyon field $T_0$ and $V(T_0)=0$, the tachyon condensate oscillates around its minimum with a decreasing amplitude. It is shown that its effective equation of state is $p=-\epsilon/3$. However, linear inhomogeneous tachyon fluctuations coupled to the oscillating background condensate are exponentially unstable due to the effect of parametric resonance. In another interesting model, where tachyon potential exponentially approaches zero at infinity of $T$, rolling tachyon condensate in an expanding universe behaves as pressureless fluid. Its linear fluctuations coupled with small metric perturbations evolve similar to these in the pressureless fluid. However, this linear stage changes to a strongly non-linear one very early, so that the usual quasi-linear stage observed at sufficiently large scales in the present Universe may not be realized in the absence of the usual particle-like cold dark matter.

Four dimensional gravity in the low energy limit of a higher dimensional theory has been expected to be a (generalized) Brans-Dicke theory. A subtle point in brane world scenarios is that the system of four dimensional effective gravitational equations is not closed due to bulk gravitational waves and bulk scalars. Nonetheless, weak gravity on the brane can be analyzed completely. We revisit the theory of weak brane gravity using gauge-invariant gravitational and scalar perturbations around a background warped geometry with a bulk scalar between two flat branes. We obtain a simple condition for the radion stabilization in terms of the scalar field potentials. We show that for general potentials of the scalar field which provides radion stabilization and a general conformal transformation to a frame in which matter on the branes are minimally coupled to the metric, 4-dimensional Einstein gravity, not BD gravity, is restored at low energies on either brane. In contrast, in RS brane world scenario without a bulk scalar, low energy gravity is BD one. We conjecture that in general brane world scenarios with more than one scalar field, one will again encounter the situation that low energy gravity is not described by the Einstein theory. Equipped with the weak gravity results, we discuss the properties of 4d brane gravitational equations, in particular, the value and sign of 4d Newton's gravitational coupling.

We investigate cosmological evolution in models where the effective potential V(\phi) may become negative for some values of the field \phi. Phase portraits of such theories in space of variables (\phi,\dot\phi,H) have several qualitatively new features as compared with phase portraits in the theories with V(\phi) > 0. Cosmological evolution in models with potentials with a "stable" minimum at V(\phi)<0 is similar in some respects to the evolution in models with potentials unbounded from below. Instead of reaching an AdS regime dominated by the negative vacuum energy, the universe reaches a turning point where its energy density vanishes, and then it contracts to a singularity with properties that are practically independent of V(\phi). We apply our methods to investigation of the recently proposed cyclic universe scenario. We show that in addition to the singularity problem there are other problems that need to be resolved in order to realize a cyclic regime in this scenario. We propose several modifications of this scenario and conclude that the best way to improve it is to add a usual stage of inflation after the singularity and use that inflationary stage to generate perturbations in the standard way.

I review the theory of preheating after inflation, focusing on the recently found tachyonic preheating in the theories with spontaneous symmetry breaking. This occurs due to the tachyonic instability of the scalar field near the top of its effective potential. Contrary to the common expectation, tachyonic instability converts most of the energy into that of colliding classical waves very rapidly, within a single oscillation. Efficient tachyonic preheating is typical for the hybrid inflationary scenario, including SUSY motivated and brane inflation models.

One of the central points of the ekpyrotic cosmological scenario based on Horava-Witten theory is that we live on a negative tension brane. However, the tension of the visible brane is positive in the usual HW phenomenology with stronger coupling on the hidden brane, both for standard and non-standard embedding. To make ekpyrotic scenario realistic one must solve the problem of the negative cosmological constant on the visible brane and fine-tune the bulk brane potential with an accuracy of $10^{-50}$. In terms of a canonically normalized scalar field $\phi$ describing the position of the brane, this potential must take a very unusual form $V(\phi)\sim e^{-{5000 \phi\over M_p}}$. We describe the problems which appear when one attempts to obtain this potential in string theory. The mechanism for the generation of density perturbations in this scenario is not brane-specific; it is a particular limiting case of the mechanism of tachyonic preheating. Unlike inflation, this mechanism exponentially amplifies not only quantum fluctuations, but also initial inhomogeneities. As a result, to solve the homogeneity problem in this scenario, one would need the branes to be parallel to each other with an accuracy better than $10^{-60}$ on a scale $10^{30}$ times greater than the distance between the branes. Thus, at present, inflation remains the only robust mechanism that produces density perturbations with a flat spectrum and simultaneously solves all major cosmological problems.

The possibility to study the M/string theory cosmology via 5d bulk & brane action is investigated. The role of the 4-form field in the theory of BPS branes in 5d is clarified. We describe arguments suggesting that the effective 4d description of the universe in the ekpyrotic scenario (hep-th/0103239) should lead to contraction rather than expansion of the universe. To verify these arguments, we study the full 5d action prior to its integration over the 5th dimension. We show that if one adds the potential V(Y) to the action of the bulk brane, then the metric ansatz used in the ekpyrotic scenario does not solve the dilaton and gravitational equations. To find a consistent cosmological solution one must use a more general metric ansatz and a complete 5d description of the brane interaction instead of simply adding an effective 4d bulk brane potential V(Y).

Spontaneous symmetry breaking usually occurs due to the tachyonic (spinodal) instability of a scalar field near the top of its effective potential at $\phi = 0$. Naively, one might expect the field $\phi$ to fall from the top of the effective potential and then experience a long stage of oscillations with amplitude O(v) near the minimum of the effective potential at $\phi = v$ until it gives its energy to particles produced during these oscillations. However, it was recently found that the tachyonic instability rapidly converts most of the potential energy V(0) into the energy of colliding classical waves of the scalar field. This conversion, which was called "tachyonic preheating," is so efficient that symmetry breaking typically completes within a single oscillation of the field distribution as it rolls towards the minimum of its effective potential. In this paper we give a detailed description of tachyonic preheating and show that the dynamics of this process crucially depend on the shape of the effective potential near its maximum. In the simplest models where $V(\phi) \sim -m^2\phi^2$ near the maximum, the process occurs solely due to the tachyonic instability, whereas in the theories $-\lambda\phi^n$ with n > 2 one encounters a combination of the effects of tunneling, tachyonic instability and bubble wall collisions.

We present a fully nonlinear study of the development of equilibrium after preheating. Preheating is the exponentially rapid transfer of energy from the nearly homogeneous inflaton field to fluctuations of other fields and/or the inflaton itself. This rapid transfer leaves these fields in a highly nonthermal state with energy concentrated in infrared modes. We have performed lattice simulations of the evolution of interacting scalar fields during and after preheating for a variety of inflationary models. We have formulated a set of generic rules that govern the thermalization process in all of these models. Notably, we see that once one of the fields is amplified through parametric resonance or other mechanisms it rapidly excites other coupled fields to exponentially large occupation numbers. These fields quickly acquire nearly thermal spectra in the infrared, which gradually propagates into higher momenta. Prior to the formation of total equilibrium, the excited fields group into subsets with almost identical characteristics (e.g. group effective temperature). The way fields form into these groups and the properties of the groups depend on the couplings between them. We also studied the onset of chaos after preheating by calculating the Lyapunov exponent of the scalar fields.

This contribution to the proceedings collects new recent results on preheating after inflation. We discuss tachyonic preheating in the SUSY motivated hybrid inflation; development of equilibrium after preheating; theory of fermionic preheating and the problem of gravitino overproduction from preheating.

We reconsider the old problem of the dynamics of spontaneous symmetry breaking using 3d lattice simulations, and develop a theory of tachyonic preheating, which occurs due to the spinodal instability of the scalar field. Tachyonic preheating is so efficient that symmetry breaking typically completes within a single oscillation of the field distribution as it rolls towards the minimum of its effective potential. As an application of this theory we consider preheating in the hybrid inflation scenario, including SUSY-motivated F-term and D-term inflationary models. We show that preheating in hybrid inflation is typically tachyonic and the stage of oscillations of a homogeneous component of the scalar fields driving inflation ends after a single oscillation. Our results may also be relevant for the theory of the formation of disoriented chiral condensates in heavy ion collisions.

We investigate the production of gravitinos in a cosmological background. Gravitinos can be produced during preheating after inflation due to a combined effect of interactions with an oscillating inflaton field and absence of conformal invariance. In order to get insight on conformal properties of gravitino we reformulate phenomenological supergravity in SU(2,2|1)-symmetric way. The Planck mass and F- and D-terms appear via the gauge-fixed value of a superfield that we call conformon. We find that in general the probability of gravitino production is not suppressed by the small gravitational coupling. This may lead to a copious production of gravitinos after inflation. Efficiency of the new non-thermal mechanism of gravitino production is very sensitive to the choice of the underlying theory. This may put strong constraints on certain classes of inflationary models.

Preheating after inflation may lead to nonthermal phase transitions with symmetry restoration. These phase transitions may occur even if the total energy density of fluctuations produced during reheating is relatively small as compared with the vacuum energy in the state with restored symmetry. As a result, in some inflationary models one encounters a secondary, nonthermal stage of inflation due to symmetry restoration after preheating. We review the theory of nonthermal phase transitions and make a prediction about the expansion factor during the secondary inflationary stage. We then present the results of lattice simulations which verify these predictions, and discuss possible implications of our results for the theory of formation of topological defects during nonthermal phase transitions.

A theory of gravitational production of light scalar particles during and after inflation is investigated. We show that in the most interesting cases where long-wavelength fluctuations of light scalar fields can be generated during inflation, these fluctuations rather than quantum fluctuations produced after inflation give the dominant contribution to particle production. In such cases a simple analytical theory of particle production can be developed. Application of our results to the theory of quantum creation of moduli fields demonstrates that if the moduli mass is smaller than the Hubble constant then these fields are copiously produced during inflation. This gives rise to the cosmological moduli problem even if there is no homogeneous component of the classical moduli field in the universe. To avoid this version of the moduli problem it is necessary for the Hubble constant H during the last stages of inflation and/or the reheating temperature T_R after inflation to be extremely small.

In this talk recent developments of the theory of preheating after inflation are briefly reviewed. In inflationary cosmology, the particles constituting the Universe are created after inflation due to their interaction with moving inflaton field(s) in the process of reheating. In inflationary models motivated by supergravity, both bosons and fermions are created. In the bosonic sector, the leading channel of particle production is the non-perturbative regime of parametric resonance dominated by those bosons which are created exponentially fast with the largest characteristic exponent. In the fermionic sector, the leading channel corresponds to the regime of parametric excitation of fermions, which respects Pauli blocking but differs significantly from the perturbative expectation. In supergravity we also have to consider production of gravitinos and moduli fields, which are cosmologically dangerous relics. We discuss the derivation of the gravitino equations in curved space-time with moving background scalars. We describe recent results on the production of gravitinos from preheating, which may put strong constraints on the inflationary models.

We study inflationary models in which the effective potential of the inflaton field does not have a minimum, but rather gradually decreases at large $\phi$. In such models the inflaton field does not oscillate after inflation, and its effective mass becomes vanishingly small, so the standard theory of reheating based on the decay of the oscillating inflaton field does not apply. For a long time the only mechanism of reheating in such non-oscillatory (NO) models was based on gravitational particle production in an expanding universe. This mechanism is very inefficient. We will show that it may lead to cosmological problems associated with large isocurvature fluctuations and overproduction of dangerous relics such as gravitinos and moduli fields. We also note that the setting of initial conditions for the stage of reheating in these models should be reconsidered. All of these problems can be resolved in the context of the recently proposed scenario of instant preheating if there exists an interaction ${g^2} \phi^2\chi^2$ of the inflaton field $\phi$ with another scalar field $\chi$. We show that the mechanism of instant preheating in NO models is much more efficient than the usual mechanism of gravitational particle production even if the coupling constant $g^2$ is extremely small, $10^{-14} \ll g^2 \ll 1$.

We examine several aspects of redshift distortions by expressing the redshift-space density in terms of the eigenvalues and orientation of the local Lagrangian deformation tensor. We explore the importance of multiple-streaming using the Zel'dovich approximation (ZA), and compute the average number of streams in real and redshift-space. It is found that multiple-streaming can be significant in redshift-space but negligible in real-space, even at moderate values of the linear fluctuation amplitude ($\sigma < 1$). Moreover, unlike their real-space counter-parts, redshift-space multiple-streams can flow past each other with minimal interactions. Such nonlinear redshift-space effects, which operate even when the real-space density field is quite linear, could suppress the classic compression of redshift-structures predicted by linear theory (Kaiser 1987). We also compute using the ZA the probability distribution function (PDF) of density, as well as $S_3$, in real and redshift-space, and compare it with the PDF measured from N-body simulations. The role of caustics in defining the character of the high density tail is examined. It is found that (non-Lagrangian) smoothing, due to both finite resolution or discreteness and small-scale velocity dispersions, is very effective in erasing caustic structures, unless the initial power spectrum is sufficiently truncated.

We describe a new efficient mechanism of reheating. Immediately after rolling down the rapidly moving inflaton field $\phi$ produces particles $\chi$, which may be either bosons or fermions. This is a nonperturbative process which occurs almost instantly; no oscillations or parametric resonance is required. The effective masses of the $\chi$ particles may be very small at the moment when they are produced, but they ``fatten'' when the field $\phi$ increases. When the particles $\chi$ become sufficiently heavy, they rapidly decay to other, lighter particles. This leads to an almost instantaneous reheating accompanied by the production of particles with masses which may be as large as $10^{17}-10^{18}$ GeV. This mechanism works in the usual inflationary models where $V(\phi)$ has a minimum, where it takes only a half of a single oscillation of the inflaton field $\phi$, but it is especially efficient in models with effective potentials slowly decreasing at large $\phi$ as in the theory of quintessence.

We review the main theoretical aspects of the structure formation paradigm which impinge upon wide angle surveys: the early universe generation of gravitational metric fluctuations from quantum noise in scalar inflaton fields; the well understood and computed linear regime of CMB anisotropy and large scale structure (LSS) generation; the weakly nonlinear regime, where higher order perturbation theory works well, and where the cosmic web picture operates, describing an interconnected LSS of clusters bridged by filaments, with membranes as the intrafilament webbing. Current CMB+LSS data favour the simplest inflation-based $\Lambda$CDM models, with a primordial spectral index within about 5% of scale invariant and $\Omega_\Lambda \approx 2/3$, similar to that inferred from SNIa observations, and with open CDM models strongly disfavoured. The attack on the nonlinear regime with a variety of N-body and gas codes is described, as are the excursion set and peak-patch semianalytic approaches to object collapse. The ingredients are mixed together in an illustrative gasdynamical simulation of dense supercluster formation.

Simple analytic arguments are used to understand the predominantly filamentary web in the large-scale distribution of galaxies. Numerical simulations of superclusters are performed to assess the feasibility of directly mapping the intracluster webbing, in particular with weak gravitational lensing using Wide Field Cameras.

We perform a weak lensing and photometric study of the z=0.42 supercluster MS0302+17 using deep I and V band images taken with the UH8K CCD mosaic camera at the CFHT. We use archival ROSAT HRI data to estimate fluxes, gas masses and, in one case, the binding mass of the three major clusters. We then use our CCD data to determine the optical richness and luminosities of the clusters and to map out the spatial distribution of the early type galaxies in the supercluster and in other foreground and background structures. We measure the gravitational shear from a sample of \simeq 30,000 faint background galaxies in the range 22 < m_I < 26 and find this correlates strongly with that predicted from the early type galaxies if they trace the mass with M/L_B \simeq 250 h. We make 2-dimensional reconstructions of the mass surface density. These recover all of the major concentrations of galaxies and indicate that most of the supercluster mass, like the early type galaxies, is concentrated in the three X-ray clusters, and we obtain mean mass-to-light ratios for the clusters of M/L_B \simeq 260 h. Cross-correlation of the measured mass surface density with that predicted from the early type galaxy distribution shows a strong peak at zero lag (significant at the \simeq 9-sigma level), and that at separations \gsim 200 h^{-1}$kpc the early galaxies trace the mass very accurately. This conclusion is supported by cross-correlation in Fourier space; we see little evidence for any variation of M/L or `bias' with scale, and from the longest wavelength modes with \lambda = 1.5-6 h^{-1}Mpc we find M/L \simeq (280 \pm 40)h, quite similar to that obtained for the cluster centers. We discuss the implication of these results for the cosmological density parameter.

We consider the instability of fluctuations in an oscillating scalar field which obeys the Sine-Gordon equation. We present simple closed-form analytic solutions describing the parametric resonance in the Sine-Gordon model. The structure of the resonance differs from that obtained with the Mathieu equation which is usually derived with the small angle approximation to the equation for fluctuations. The results are applied to axion cosmology, where the oscillations of the classical axion field, with a Sine-Gordon self-interaction potential, constitute the cold dark matter of the universe. When the axion misalignment angle at the QCD epoch, $\theta_0$, is small, the parametric resonance of the axion fluctuations is not significant. However, in regions of larger $\theta_0$ where axion miniclusters would form, the resonance may be important. As a result, axion miniclusters may disintegrate into finer, denser clumps. We also apply the theory of Sine-Gordon parametric resonance to reheating in the Natural Inflation scenario. The decay of the inflaton field due to the self-interaction alone is ineffective, but a coupling to other bosons can lead to preheating in the broad resonance regime. Together with the preheating of fermions, this can alter the reheating scenario for Natural Inflation.

In inflationary cosmology, the particles constituting the Universe are created after inflation in the process of reheating due to their interaction with the oscillating inflaton field. In the bosonic sector, the leading channel of particle production is the non-perturbative regime of parametric resonance, preheating, during which bosons are created exponentially fast. Pauli blocking prohibits the unbounded creation of fermions. For this reason, it has been silently assumed that the creation of fermions can be treated with perturbation theory for the decay of individual inflatons. We consider the production of fermions interacting with the coherently oscillating inflatons. We find that the actual particle production occurs in a regime of the parametric excitation of fermions, leading to preheating of fermions. Fermion preheating differs significantly from the perturbative expectation. It turns out that the number density of fermions varies periodically with time. The total number of fermions quickly saturates to an average value within a broad range of momenta $\propto q^{1/4}$, where $q$ is the usual resonance parameter. The resonant excitation of fermions may affect the transfer inflaton energy, estimations of the reheating temperature, and the abundance of superheavy fermions and gravitinos. Back in the bosonic sector, outside of the parametric resonance bands there is an additional effect of parametric excitation of bosons with bounded occupation number in the momentum range $\propto q^{1/4}$.

We investigate nonthermal phase transitions that may occur after post-inflationary preheating in a simple model of a two-component scalar field with the effective potential $\lambda (\phi_i^2 - {\rm v}^2)^2/4$, where $\phi_1$ is identified with the inflaton field. We use three-dimensional lattice simulations to investigate the full nonlinear dynamics of the model. Fluctuations of the fields generated during and after preheating temporarily make the effective potential convex in the $\phi_1$ direction. The subsequent nonthermal phase transition with symmetry breaking leads to formation of cosmic strings even for ${\rm v} \gg 10^{16}$ GeV. This mechanism of string formation, in a modulated (by the oscillating field $\phi_{1}$) phase transition, is different from the usual Kibble mechanism.

During preheating after inflation, parametric resonance rapidly generates very large fluctuations of scalar fields. In models where the inflaton field $\phi$ oscillates in a double-well potential and interacts with another scalar field $X$, fluctuations of X can keep the \phi to -\phi symmetry temporarily restored. If the coupling of \phi to X is much stronger than the inflaton self-coupling, the subsequent symmetry breaking is a first-order phase transition. We demonstrate the existence of this nonthermal phase transition with lattice simulations of the full nonlinear dynamics of the interacting fields. In particular, we observe nucleation of an expanding bubble.

In inflationary cosmology, the particles constituting the Universe are created after inflation in the process of reheating due to the interaction with the oscillating inflaton field. We briefly review the basics of the slow reheating, and the stage of fast preheating, when the particles are created explosively in the regime of parametric resonance. The non-perturbative, out-of-equilibrium character of the parametric resonance changes many features of reheating. For these proceedings, we will highlight a few aspects of preheating: the structural dependence of the parametric resonance on the inflationary model, including $V(\phi)={m^2 \over 2}\phi^2$, ${\lambda\over 4}\phi^4$, $1-\cos {\phi \over f}$; ``rescattering'' of created particles; and phase transitions after inflation.

It is assumed that during inflation, all energy was contained in a slow-rolling inflaton field $\phi$. The particles constituting the Universe are created due to interactions with the field $\phi$ coherently oscillating after inflation. The leading channel of the particle production would be the non-perturbative regime of parametric resonance, preheating. Some of the recent developments in the theory of preheating are briefly reviewed. As a prototype we use the model ${\lambda\over 4}\phi^4 + {m^2 \over 2}\phi^2+ {g^2\over 2}\phi^2\chi^2$ in an expanding universe. In different domains of parameters the character of preheating $(\phi \to \chi)$ or $(\phi \to \delta \phi)$ is different, ranging from regular Lam\'{e} to stochastic parametric resonances. Fortunately, in many important cases simple analytic resonant solutions are found. Eventually, a picture which unifies the different regimes emerges.

Reheating after inflation occurs due to particle production by the oscillating inflaton field. In this paper we describe the perturbative approach to reheating, and then concentrate on effects beyond the perturbation theory. They are related to the stage of parametric resonance called preheating. It may occur in an expanding universe if the initial amplitude of oscillations of the inflaton field is large enough. We investigate a simple model of a massive inflaton field coupled to another scalar field X. Parametric resonance in this model is very broad. It occurs in a very unusual stochastic manner, which is different from the parametric resonance in the case when the expansion of the universe is neglected. Quantum fields interacting with the oscillating inflaton field experience a series of kicks which occur with phases uncorrelated to each other. We call this process stochastic resonance. We develop the theory of preheating taking into account the expansion of the universe and backreaction of produced particles, including the effects of rescattering. The process of preheating can be divided into several distinct stages. At the first stage the backreaction of created particles is not important. At the second stage backreaction increases the frequency of oscillations of the inflaton field, which makes the process even more efficient than before. Then the effects related to scattering of X-particles terminate the resonance. We calculate the density of X-particles and their quantum fluctuations with all backreaction effects taken into account. This allows us to find the range of masses and coupling constants for which one has efficient preheating. In particular, under certain conditions this process may produce particles with a mass much greater than the mass of the inflaton field.

We consider preheating in the theory $1/4 \lambda \phi^4 + 1/2 g^2\phi^2\chi^2 $, where the classical oscillating inflaton field $\phi$ decays into $\chi$-particles and $\phi$-particles. The parametric resonance which leads to particle production in this conformally invariant theory is described by the Lame equation. It significantly differs from the resonance in the theory with a quadratic potential. The structure of the resonance depends in a rather nontrivial way on the parameter $g^2/\lambda$. We construct the stability/instability chart in this theory for arbitrary $g^2/\lambda$. We give simple analytic solutions describing the resonance in the limiting cases $g^2/\lambda\ll 1$ and $g^2/\lambda \gg 1$, and in the theory with $g^2=3\lambda$, and with $g^2 =\lambda$. From the point of view of parametric resonance for $\chi$, the theories with $g^2=3\lambda$ and with $g^2 =\lambda$ have the same structure, respectively, as the theory $1/4 \lambda \phi^4$, and the theory $\lambda /(4 N) (\phi^2_i)^2$ of an N-component scalar field $\phi_i$ in the limit $N \to \infty$. We show that in some of the conformally invariant theories such as the simplest model $1/4 \lambda\phi^4$, the resonance can be terminated by the backreaction of produced particles long before $<\chi^2>$ or $<\phi^2 >$ become of the order $\phi^2$. We analyze the changes in the theory of reheating in this model which appear if the inflaton field has a small mass.

In Paper I we studied the theory of gravitational microlensing for a planar distribution of point masses. In this second paper, we extend the analysis to a three-dimensional lens distribution. First we study the lensing properties of three-dimensional lens distributions by considering in detail the critical curves, the caustics, the illumination patterns, and the magnification cross-sections sigma(A) of multiplane configurations with two, three, and four point masses. For N_\ast point masses that are widely separated in Lagrangian space, we find that there are ~ 2^{N_\ast} - 1 critical curves in total, but that only ~ N_\ast of these produce prominent caustic-induced features at the high magnification end of sigma(A). In the case of a low optical depth random distribution of point masses, we show that the multiplane lens equation near a point mass can be reduced to the single plane equation of a point mass perturbed by weak shear. This allows us to calculate the caustic-induced feature in the macroimage magnification distribution P(A) as a weighted sum of the semi-analytic feature derived in Paper I for a planar lens distribution. The resulting semi-analytic caustic-induced feature is similar to the feature in the planar case, but it does not have any simple scaling properties, and it is shifted to higher magnification. The semi-analytic distribution is compared to the results of previous numerical simulations for optical depth tau ~ 0.1, and they are in better agreement than a similar comparison in the planar case. We explain this by estimating the fraction of caustics of individual lenses that merge with those of their neighbors. For tau=0.1, the fraction is approx 20%, much less than the approx 55% for the planar case.

(Abridged) In this paper we refine the theory of microlensing for a planar distribution of point masses. We derive the macroimage magnification distribution P(A) at high magnification (A-1 >> tau^2) for a low optical depth (tau << 1) lens distribution by modeling the illumination pattern as a superposition of the patterns due to individual ``point mass plus weak shear'' lenses. We show that a point mass plus weak shear lens produces an astroid- shaped caustic and that the magnification cross-section obeys a simple scaling property. By convolving this cross-section with the shear distribution, we obtain a caustic-induced feature in P(A) which also exhibits a simple scaling property. This feature results in a 20% enhancement in P(A) at A approx 2/tau. In the low magnification (A-1 << 1) limit, the macroimage consists of a bright primary image and a large number of faint secondary images formed close to each of the point masses. Taking into account the correlations between the primary and secondary images, we derive P(A) for low A. The low-A distribution has a peak of amplitude ~ 1/tau^2 at A-1 ~ tau^2 and matches smoothly to the high-A distribution. We combine the high- and low-A results and obtain a practical semi-analytic expression for P(A). This semi-analytic distribution is in qualitative agreement with previous numerical results, but the latter show stronger caustic-induced features at moderate A for tau as small as 0.1. We resolve this discrepancy by re-examining the criterion for low optical depth. A simple argument shows that the fraction of caustics of individual lenses that merge with those of their neighbors is approx 1-exp(-8 tau). For tau=0.1, the fraction is surprisingly high: approx 55%. For the purpose of computing P(A) in the manner we did, low optical depth corresponds to tau << 1/8.

In the inflationary scenario all the matter constituting the universe was created from the process of reheating after inflation. Recent development of the theory of reheating is briefly reviewed. The list of topics includes elementary (perturbative) theory of reheating; quantum field theory in a time-varying background; parametric resonance and explosive particle creation; non-thermal phase transitions from reheating; baryogenesis from reheating; residual oscillations of the scalar field, and other cosmological applications.

At the first stage of reheating after inflation, parametric resonance may rapidly transfer most of the energy of an inflaton field $\phi$ to the energy of other bosons. We show that quantum fluctuations of scalar and vector fields produced at this stage are much greater than they would be in a state of thermal equilibrium. This leads to cosmological phase transitions of a new type, which may result in a copious production of topological defects and in a secondary stage of inflation after reheating.

Observations indicate galaxies are distributed in a filament-dominated web-like structure. Numerical experiments at high and low redshift of viable structure formation theories also show filament-dominance. We present a simple quantitative explanation of why this is so, showing that the final-state web is actually present in embryonic form in the overdensity pattern of the initial fluctuations, with nonlinear dynamics just sharpening the image. The web is largely defined by the position and primordial tidal fields of rare events in the medium, with the strongest filaments between nearby clusters whose tidal tensors are nearly aligned. Applications of the cosmic web theory to observations include probing cluster-cluster bridges by weak gravitational lensing, X-rays, and the Sunyaev-Zeldovich effect and probing high redshift galaxy-galaxy bridges by low column density Lyman alpha absorption lines.

We discuss the structure of clusters in a class of flat cosmological models with the fraction of mass \Omega_{CDM} ~0.8 in cold dark matter, and the rest in hot dark matter in the form of massive neutrinos. We consider such Cold+Hot Dark Matter (CHDM) models with one, two or three massive neutrinos, with total mass ~4.6eV. Neutrinos of such low mass cannot constitute halos of galaxies and groups, but only of clusters of galaxies. The limit on the density of neutrinos in the central parts of galaxy clusters is estimated from the phase space den- sity constraints. The ratio of the neutrino density to that of CDM through the cluster is found analytically. It appears that the density of neutrinos is suppressed within the Abell radius. However, neutrinos contribute ~20% of the mass density to the cluster halo. Our numerical simulations match analytical results. The simulations indicate that the cluster halo dark matter density profile has the power-law slope ~-2.5 which is close to that in the model with cosmological constant. We also found that in the CHDM models the velocity dispersion is almost constant across the cluster. This is quite different from the model with cosmological constant or the open model where the velocity dispersion falls in the cluster outskirts. We discuss X-ray emission and weak gravitational lensing by clusters in the model. We input the found spherically symmetrical fit to the CHDM mass density profile and the X-ray surface brightness for the cluster A2256 into simple equation of hydrostatic equilibrium of the hot gas. X-ray temperature derived this way departures from both the data and actual prediction of the model, which give almost constant temperature. We found also that the problem of high baryonic fraction in clusters is not resolved in the CHDM models.

This contribution to the Proceedings is based on the talk given at the Conference on Birth of the Universe and Fundamental Physics, Rome, May 18-21, 1994. Some selected topics of the subject are reviewed: Models of Primordial Fluctuations; Reconstruction of the Cosmological Density Probability Distribution Function (PDF) from Cumulants; PDFs from the Zel'dovich Approximation and from Summarizing Perturbation Series; Fitting by the Log-normal Distribution.

The theory of reheating of the Universe after inflation is developed. The transition from inflation to the hot Universe turns out to be strongly model-dependent and typically consists of several stages. Immediately after inflation the field $\phi$ begins rapidly rolling towards the minimum of its effective potential. Contrary to some earlier expectations, particle production during this stage does not lead to the appearance of an extra friction term $\Gamma\dot\phi$ in the equation of motion of the field $\phi$. Reheating becomes efficient only at the next stage, when the field $\phi$ rapidly oscillates near the minimum of its effective potential. We have found that typically in the beginning of this stage the classical inflaton field $\phi$ very rapidly (explosively) decays into $\phi$-particles or into other bosons due to broad parametric resonance. This stage cannot be described by the standard elementary approach to reheating based on perturbation theory. The bosons produced at this stage, as well as some part of the classical field $\phi$ which survives the stage of explosive reheating, should further decay into other particles, which eventually become thermalized. The last stages of decay can be described in terms of perturbation theory. Complete reheating is possible only in those theories where a single massive $\phi$-particle can decay into other particles. This imposes strong constraints on the structure of inflationary models. On the other hand, this means that a scalar field can be a cold dark matter candidate even if it is strongly coupled to other fields.

Few recent generations of cosmologists have solved non-local newtonian equations of the gravitational instability in an expanding universe. In this approach pancaking is the predominant form of first collapsing objects. Relativistic counterparts of these equations contain the electric and magnetic parts of the Weyl tensor. In the linear theory the magnetic part is associated with gravitational waves. If the magnetic part is ignored, then the newtonian limit of the relativistic equations is reduced to the closed set of the local Lagrangian equations. Recently this fact drew much attention since the gravitational instability in that form would greatly simplify the study of cosmic structure formation. In particular, the filamentary structure of collapsing is predicted. In this paper we resolve the contradiction between the newtonian theory and relativistic version adopted in some recent papers. We show that dropping the magnetic part from the basic relativistic equations is {\it incorrect}. The correct newtonian limit is derived by the $1/c$-expansion of the GR equations and the Bianchi identities for the Weyl tensor. The last ones begin with $\sim 1/c^3$ order, therefore one {\it must} take into account the magnetic part in the post newtonian order $\sim 1/c^3$, which contains non-local terms, related to the non-local gravitational interaction. For the first time we rigorously show that the GR equations with the magnetic part are reduced precisely to the canonic newtonian non-local equations. Thus, the correct treatment of the relativistic version of the gravitational instability resurrects the canonic picture of the structure formation.

The properties of the probability distribution function of the cosmological continuous density field are studied. We present further developments and compare dynamically motivated methods to derive the PDF. One of them is based on the Zel'dovich approximation (ZA). We extend this method for arbitrary initial conditions, regardless of whether they are Gaussian or not. The other approach is based on perturbation theory with Gaussian initial fluctuations. We include the smoothing effects in the PDFs. We examine the relationships between the shapes of the PDFs and the moments. It is found that formally there are no moments in the ZA, but a way to resolve this issue is proposed, based on the regularization of integrals. A closed form for the generating function of the moments in the ZA is also presented, including the smoothing effects. We suggest the methods to build PDFs out of the whole series of the moments, or out of a limited number of moments -- the Edgeworth expansion. The last approach gives us an alternative method to evaluate the skewness and kurtosis by measuring the PDF around its peak. We note a general connection between the generating function of moments for small r.m.s $\sigma$ and the non-linear evolution of the overdense spherical fluctuation in the dynamical models. All these approaches have been applied in 1D case where the ZA is exact, and simple analytical results are obtained. The 3D case is analyzed in the same manner and we found a mutual agreement in the PDFs derived by different methods in the the quasi-linear regime. Numerical CDM simulation was used to validate the accuracy of considered approximations. We explain the successful log-normal fit of the PDF from that simulation at moderate $\sigma$ as mere fortune, but not as a universal form of density PDF in general.

We study the quasilinear evolution of the one-point probability density functions (PDFs) of the smoothed density and velocity fields in a cosmological gravitating system beginning with Gaussian initial fluctuations. Our analytic results are based on the Zel'dovich approximation and laminar flow. A numerical analysis extends the results into the multistreaming regime using the smoothed fields of a CDM N-body simulation. We find that the PDF of velocity, both Lagrangian and Eulerian, remains Gaussian under the laminar Zel'dovich approximation, and it is almost indistinguishable from Gaussian in the simulations. The PDF of mass density deviates from a normal distribution early in the quasilinear regime and it develops a shape remarkably similar to a lognormal distribution with one parameter, the \rms density fluctuation $\sigma$. Applying these results to currently available data we find that the PDFs of the velocity and density fields, as recovered by the \pot\ procedure from observed velocities assuming $\Omega=1$, or as deduced from a redshift survey of \iras\ galaxies assuming that galaxies trace mass, are consistent with Gaussian initial fluctuations.

The evolution of probability distribution functions (PDFs) of continuous density, velocity and velocity derivatives ( deformation tensor) fields in the theory of cosmological gravitational instability are considered. We show that in the Newtonian theory the dynamical equations cannot be reduced to the closed set of Lagrangian equations. Since continuous fields from galaxy surveys need sufficiently large smoothing which exceeds the scale of nonlinearity, one can use the Zel'dovich approximation to describe the mildly non-linear matter evolution, which allows the closed set of Lagrangian equations. The closed kinetic equation for the joint PDF of cosmological continuous fields is derived in this approximation. The analytical theory of the cosmological PDFs with arbitrary (including Gaussian) initial statistics is developed, based on the solution on the kinetic equaiton. For Gaussian initial fluctuations, the PDFs are parametrized by only linear {\it rms} fluctuations $\sigma$ on given filtering scale. Density PDF $P(\rho, t)$ and PDF $M(\lambda_1, \lambda_2, \lambda_3; t)$ of eigenvalues of the deformation tensor field in the Eulerian space evolve very rapidly in non-linear regime. On the contrary, velocity PDF $Q(\vec v, t)$ remains invariant under non-linear evolution. For small $\sigma$ the Edgworth series is suggested to reconstruct