We consider Higgs inflation with an $\alpha R^2$ term. It adds a new scalar degree of freedom, which leads to a two-field model of inflation. We do a complete slow-roll analysis of the three-dimensional parameter space of the $R^2$ coefficient $\alpha$, the non-minimal coupling $\xi$ and the Higgs self-coupling $\lambda$. We find three classes of inflationary solutions, but only pure $R^2$ and attractor solutions fit observations. We find that pure Higgs inflation is impossible when the $R^2$ term is present regardless of how small $\alpha$ is. However, we can have Higgs-like inflation, where the amplitude of the perturbations does not depend on $\alpha$ and the predictions as a function of e-folds are the same as in Higgs inflation, although the inflationary trajectory is curved in field space. The spectral index is $0.939 < n_R < 0.967$, and constraining it to the observed range, the tensor-to-scalar ratio varies from $3.8\times10^{-3}$ to the maximum allowed by observations, $0.079$. Observational constraints on isocurvature perturbations contribute to these limits, whereas non-Gaussianity is automatically in the range allowed by observations.

We revisit a cosmological constraint on dark matter decaying into dark radiation at late times. In Enqvist et al. (2015), we mainly focused on the effects of decaying dark matter (DDM) on the cosmic microwave background (CMB) and nonlinear matter power spectrum. Extending our previous analysis, here we use N-body simulation to investigate how DDM affects the halo mass function. This allows us to incorporate the cluster counts observed by the Sunyaev-Zel'dovich effect to study a bound on the lifetime of DDM. We also update the data of CMB and cosmic shear power spectrum with the Planck 2015 results and KiDS450 observations, respectively. From these cosmological observations, we obtain an lower bound on the lifetime $\Gamma^{-1}\ge 175\,$Gyr from the Planck2015 results (CMB+SZ cluster count) combined with the KiDS450 and the recent measurements of the baryon acoustic scale.

We discuss structure formation in models with a spectator field in small-field inflation which accommodate a secondary period of inflation. In such models, subgalactic scale primordial fluctuations can be much suppressed in comparison to the usual power-law $\Lambda$CDM model while the large scale fluctuations remain consistent with current observations. We discuss how a secondary inflationary epoch may give rise to observable features in the small scale power spectrum and hence be tested by the structures in the Local Universe.

We study scalar field inflation in $F(R)$ gravity in the Palatini formulation of general relativity. Unlike in the metric formulation, in the Palatini formulation $F(R)$ gravity does not introduce new degrees of freedom. However, it changes the relations between existing degrees of freedom, including the inflaton and spacetime curvature. Considering the case $F(R)=R+\alpha R^2$, we find that the $R^2$ term decreases the height of the effective inflaton potential. By adjusting the value of $\alpha$, this mechanism can be used to suppress the tensor-to-scalar ratio $r$ without limit in any scalar field model of inflation without affecting the spectrum of scalar perturbations.

We study inflation with the non-minimally coupled Standard Model Higgs in the case when quantum corrections generate a hilltop in the potential. We consider both the metric and the Palatini formulation of general relativity. We investigate hilltop inflation in different parts of the Higgs potential and calculate predictions for CMB observables. We run the renormalization group equations up from the electroweak scale and down from the hilltop, adding a jump in-between to account for unknown corrections in the intermediate regime. Within our approximation, no viable hilltop inflation is possible for small field values, where the non-minimal coupling has no role, nor for intermediate field values. For large field values, hilltop inflation works. We find the spectral index to be $n_s\leq0.96$ in both the metric and the Palatini formulation, the upper bound coinciding with the tree-level result. The tensor-to-scalar ratio is $r\leq1.2\times10^{-3}$ in the metric case and $r\leq2.2\times10^{-9}$ in the Palatini case. Successful inflation is possible even when the renormalization group running is continuous with no jumps. In the metric formulation, $r$ is smaller than in Higgs inflation on the tree-level plateau or at the critical point, making it possible to distinguish hilltop inflation from these scenarios with next-generation CMB experiments.

We show that in the Feebly Interacting Massive Particle (FIMP) model of Dark Matter (DM), one may express the inflationary energy scale $H_*$ as a function of three otherwise unrelated quantities, the DM isocurvature perturbation amplitude, its mass and its self-coupling constant, independently of the tensor-to-scalar ratio. The FIMP model assumes that there exists a real scalar particle that alone constitutes the DM content of the Universe and couples to the Standard Model via a Higgs portal. We consider carefully the various astrophysical, cosmological and model constraints, accounting also for variations in inflationary dynamics and the reheating history, to derive a robust estimate for $H_*$ that is confined to a relatively narrow range. We point out that, within the context of the FIMP DM model, one may thus determine $H_*$ reliably even in the absence of observable tensor perturbations.

The Higgs-inflaton coupling plays an important role in the Higgs field dynamics in the early Universe. Even a tiny coupling generated at loop level can have a dramatic effect on the fate of the electroweak vacuum. Such Higgs-inflaton interaction is present both at the trilinear and quartic levels in realistic reheating models. In this work, we examine the Higgs dynamics during the preheating epoch, focusing on the effects of the parametric and tachyonic resonances. We use lattice simulations and other numerical tools in our studies. We find that the resonances can induce large fluctuations of the Higgs field which destabilize the electroweak vacuum. Our considerations thus provide an upper bound on quartic and trilinear interactions between the Higgs and the inflaton. We conclude that there exists a favourable range of the couplings within which the Higgs field is stabilized during both inflation and preheating epochs.

We investigate the dependency of Higgs inflation on the non-renormalisable matching between the low energy Standard Model limit and the inflationary regime at high energies. We show that for the top mass range $m_t \gtrsim 171.8$ GeV the scenario robustly predicts the spectral index $n_s \simeq 0.97$ and the tensor-to-scalar ratio $r\simeq 0.003$. The matching is however non-trivial, even the best-fit values $m_h=125.09$ GeV and $m_t=173.21$ GeV require a jump $\delta \lambda \sim 0.01$ in the Higgs coupling below the inflationary scale. For $m_t\lesssim 171.8$ GeV, the matching may generate a feature in the inflationary potential. In this case the predicted values of $n_s$ and $r$ vary but the model is still falsifiable. For example, a detection of negative running of spectral index at level $\alpha_s \lesssim -0.01$ would rule out Higgs inflation.

We study the resonant decay of the primordial Standard Model Higgs condensate after inflation into $SU(2)$ gauge bosons on the lattice. We find that the non-Abelian interactions between the gauge bosons quickly extend the momentum distribution towards high values, efficiently destroying the condensate after the onset of backreaction. For the inflationary scale $H = 10^8$ GeV, we find that 90% of the Higgs condensate has decayed after $n \sim 10$ oscillation cycles. This differs significantly from the Abelian case where, given the same coupling strengths, most of the condensate would persist after the resonance.

We discuss mixed inflaton and spectator field models where both the fields are responsible for the observed density fluctuations. We use the current CMB data to constrain both the general mixed model as well as some specific representative scenarios, and collate the results with the model predictions for the CMB spectral $\mu$ distortion. We find the posterior distribution of $\mu$ using MCMC chains and demonstrate that the standard single-field inflaton model typically predicts $\mu \sim 10^{-8}$ with a relatively narrow distribution, whereas for the mixed models, the distribution turns out to be much broader, and $\mu$ could be larger by almost an order of magnitude. Hence future experiments of $\mu$ distortion could provide a tool for the critical testing of the mixed source models of the primordial perturbation.

We consider decaying dark matter (DDM) as a resolution to the possible tension between cosmic microwave background (CMB) and weak lensing (WL) based determinations of the amplitude of matter fluctuations, $\sigma_8$. We perform N-body simulations in a model where dark matter decays into dark radiation and develop an accurate fitting formula for the non-linear matter power spectrum, which enables us to test the DDM model by the combined measurements of CMB, WL and the baryon acoustic oscillation (BAO). We employ a Markov chain Monte Carlo analysis to examine the overlap of posterior distributions of the cosmological parameters, comparing CMB alone with WL+BAO. We find an overlap that is significantly larger in the DDM model than in the standard CDM model. This may be hinting at DDM, although current data is not constraining enough to unambiguously favour a non-zero dark matter decay rate $\Gamma$. From the combined CMB+WL data, we obtain a lower bound $\Gamma^{-1}\ge 97$ Gyr at 95 % C.L, which is less tight than the constraint from CMB alone.

We discuss correlations among spectral observables as a new tool for differentiating between models for the primordial perturbation. We show that if generated in the isocurvature sector, a running of the scalar spectral index is correlated with the statistical properties of non-Gaussianities. In particular, we find a large running will inevitably be accompanied by a large running of $f_{\rm NL}$ and enhanced $g_{\rm NL}$, with $g_{\rm NL}\gg f_{\rm NL}^2$. If the tensor to scalar ratio is large, a large negative running must turn positive on smaller scales. Interestingly, the characteristic scale of the transition could potentially distinguish between the inflaton and isocurvature fields.

We study the resonant decay of the Higgs condensate into weak gauge bosons after inflation and estimate the corrections arising from the non-Abelian self-interactions of the gauge fields. We find that non-Abelian interaction terms induce an effective mass which tends to shut down the resonance. For the broad resonance relevant for the Standard Model Higgs the produced gauge particles backreact on the dynamics of the Higgs condensate before the non-Abelian terms grow large. The non-Abelian terms can however significantly affect the final stages of the resonance after the backreaction. In the narrow resonance regime, which may be important for extensions of the Standard Model, the non-Abelian terms affect already the linear stage and terminate the resonance before the Higgs condensate is affected by the backreaction of decay products.

We study the post-inflationary dynamics of the Standard Model Higgs and a real singlet scalar $s$, coupled together through a renormalizable coupling $\lambda_{sh}h^2s^2$, in a $Z_2$ symmetric model that may explain the observed dark matter abundance and/or the origin of baryon asymmetry. The initial values for the Higgs and $s$ condensates are given by inflationary fluctuations, and we follow their dissipation and relaxation to the low energy vacua. We find that both the lowest order perturbative and the non-perturbative decays are blocked by thermal effects and large background fields and that the condensates decay by two-loop thermal effects. Assuming instant reheating at $T=10^{16}$ GeV, the characteristic temperature for the Higgs condensate thermalization is found to be $T_h \sim 10^{14}$ GeV, whereas $s$ thermalizes typically around $T_s \sim 10^{6}$ GeV. By that time, the amplitude of the singlet is driven very close to the vacuum value by the expansion of the universe, unless the portal coupling takes a value $\lambda_{sh}\lesssim 10^{-7}$ and the singlet $s$ never thermalizes. With these values of the coupling, it is possible to slowly produce a sizeable fraction of the observed dark matter abundance via singlet condensate fragmentation and thermal Higgs scattering. Physics also below the electroweak scale can therefore be affected by the non-vacuum initial conditions generated by inflation.

We investigate inflationary Higgs dynamics and constraints on the Standard Model parameters assuming the Higgs potential, computed to next-to-next leading order precision, is not significantly affected by new physics. For a high inflationary scale $H\sim 10^{14}$ GeV suggested by BICEP2, we show that the Higgs is a light field subject to fluctuations which affect its dynamics in a stochastic way. Starting from its inflationary value the Higgs must be able to relax to the Standard Model vacuum well before the electroweak scale. We find that this is consistent with the high inflationary scale only if the top mass $m_t$ is significantly below the best fit value. The region within $2\sigma$ errors of the measured $m_t$, the Higgs mass $m_h$ and the strong coupling $\alpha_s$ and consistent with inflation covers approximately the interval $m_t \lesssim 171.8\,{\rm GeV} + 0.538(m_h-125.5\,{\rm GeV})$ with $125.4\,{\rm GeV}\lesssim m_h\lesssim 126.3\,{\rm GeV}$. If the low top mass region could be definitively ruled out, the observed high inflationary scale alone, if confirmed, would seem to imply new physics necessarily modifying the Standard Model Higgs potential below the inflationary scale.

We present the first full study of the minimal curvaton-higgs (MCH) model, which is a minimal interpretation of the curvaton scenario with one real scalar coupled to the standard model Higgs boson. The standard model coupling allows the dynamics of the model to be determined in detail, including effects from the thermal background and from radiative corrections to the potential. The relevant mechanisms for curvaton decay are incomplete non-perturbative decay (delayed by thermal blocking), followed by decay via a dimension-5 non-renormalisable operator. To avoid spoiling the predictions of big bang nucleosynthesis, we find the "bare" curvaton mass to be m_\sigma > 8 x 10^4 GeV. To match observational data from Planck there is an upper limit on the curvaton-higgs coupling g, between 10^-3 and 10^-2, depending on the mass. This is due to interactions with the thermal background. We find that typically non-Gaussianities are small but that if fnl is observed in the near future then m_\sigma < 5 x 10^9 GeV, depending on Hubble scale during inflation. In a thermal dark matter model, the lower bound on m_\sigma\ can increase substantially. The parameter space may also be affected once the baryogenesis mechanism is specified.

We consider the modifications of monomial chaotic inflation models due to radiative corrections induced by inflaton couplings to bosons and/or fermions necessary for reheating. To the lowest order, ignoring gravitational corrections and treating the inflaton as a classical background field, they are of the Coleman-Weinberg type and parametrized by the renormalization scale $\mu$. In cosmology, there are not enough measurements to fix $\mu$ so that we end up with a family of models, each having a slightly different slope of the potential. We demonstrate by explicit calculation that within the family of chaotic $\phi^2$ models, some may be ruled out by Planck whereas some remain perfectly viable. In contrast, radiative corrections do not seem to help chaotic $\phi^{4}$ models to meet the Planck constraints.

The behaviour of oscillating scalar spectator fields after inflation depends on the thermal background produced by inflaton decay. Resonant decay of the spectator is often blocked by large induced thermal masses. We account for the finite decay width of the inflaton and the protracted build-up of the thermal bath to determine the early evolution of a homogeneous spectator field, \sigma, coupled to the Higgs Boson, \Phi, through the term g^2 \sigma^2 \Phi^2, the only renormalisable coupling of a new scalar to the Standard Model. We find that for very large higgs-spectator coupling g > 10^{-3}, the resonance is not always blocked as was previously suggested. As a consequence, the oscillating spectator can decay quickly. For other parameter values, we find that although qualitative features of the thermal blocking still hold, the dynamics are altered compared to the instant decay case. These findings are important for curvaton models, where the oscillating field must be relatively long lived in order to produce the curvature perturbation. They are also relevant for other spectator fields, which must decay sufficiently early to avoid spoiling the predictions of baryogenesis and nucleosynthesis.

We investigate the dynamics of the Standard Model higgs with a minimal coupling to gravity during and after inflation. In the regime where the Standard Model vacuum is stable, we find that the higgs becomes a light spectator field after about 30 efolds of inflation, irrespectively of its initial value. Once the higgs has become light, its root-mean-square value h_* relaxes to equilibrium in about 85 efolds for the inflationary scale of H_*=10^4 GeV and in 20 efolds for H_*=10^10 GeV. The equilibrium value is given by h_*~0.36 lambda^{-1/4}H_*, where lambda=0.09 ... 0.0005 is the higgs self coupling at the scales H_*=10^4 ... 10^10 GeV. We show that the main decay channel of the higgs condensate after inflation is the resonant production of Standard Model gauge bosons. For a set of parameters we find that a significant part of the condensate has decayed in between 340 and 630 Hubble times after the onset of higgs oscillations, depending on H_* in a non-trivial way. The higgs perturbations correspond to isocurvature modes during inflation but they could generate significant adiabatic perturbations at a later stage for example through a modulation of the reheating stage. However, this requires that the inflaton(s) decay no later than a few hundred Hubble times after the onset of higgs oscillations.

We investigate the possibility that the primordial perturbation has two sources: the inflaton and a spectator field, which is not dynamically important during inflation but which after inflation can contribute to the curvature perturbation. The recent Planck results on the power spectrum and non-Gaussianity allow us to put constraints on such mixed models. In the generic case, where no specific model for the inflaton or the spectator is assumed, one finds that in the mixed scenario it is possible to have a large trispectrum with tau_NL >> (f_NL)^2. The constraints on inflation models in the plane of the spectral index and tensor-to-scalar ratio are modified by the presence of a spectator and depend also on the ratio of the spectator-to-inflaton power R. If one chooses the spectator to be the curvaton with a quadratic potential, non-Gaussianities can be computed and imply restrictions on possible values of R. We also consider a mixed curvaton and chaotic inflation model and show that even quartic chaotic inflation is still feasible in the context of mixed models.

We consider a model of preheating where the coupling of the inflaton to the preheat field is modulated by an additional scalar field which is light during inflation. We establish that such a model produces the observed curvature perturbation analogously to the modulated reheating scenario. The contribution of modulated preheating to the power spectrum and to non-Gaussianity can however be significantly larger compared to modulated perturbative reheating. We also consider the implications of the current constraints on isocurvature perturbations in case where the modulating field is responsible for cold dark matter. We find that existing bounds on CDM isocurvature perturbations imply that modulated preheating is unlikely to give a dominant contribution to the curvature perturbation and that the same bounds suggest important constraints on non-Gaussianity and the amount of primordial gravitational waves.

We investigate the stochastic dynamics of the long wavelength modes of a generic light scalar field that during inflation is coupled to another scalar field. The coupling plays an important role for the fluctuation of the field amplitude and may block its initial growth. We find that such a blocking is avoided, albeit only temporarily, if the light scalar has an initial non-vanishing expectation value <\phi> larger than a certain critical value, for which we provide an estimate. We also show that the field fluctuations will eventually reach an equilibrium amplitude provided inflation is sufficiently long-lasting. We present a novel, general expression for the variance <\phi^2> that takes into account the coupling of the massless field and describes the growth during the epoch of quasi-free fluctuations as well as the late-time approach to the equilibrium.

We investigate in detail a model where the curvaton is coupled to the Standard Model higgs. Parametric resonance might be expected to cause a fast decay of the curvaton, so that it would not have time to build up the curvature perturbation. However, we show that this is not the case, and that the resonant decay of the curvaton may be delayed even down to electroweak symmetry breaking. This delay is due to the coupling of the higgs to the thermal background, which is formed by the Standard Model degrees of freedom created from the inflaton decay. We establish the occurrence of the delay by considering the curvaton evolution and the structure of the higgs resonances. We then provide analytical expressions for the delay time, and for the subsequent resonant production of the higgs, which ultimately leads to the curvaton effective decay width. Contrary to expectations, it is possible to obtain the observed curvature perturbation for values of the curvaton-higgs coupling as large as 0.1. Our calculations also apply in the general case of curvaton decay into any non Standard Model species coupled to the thermal background.

We investigate the stochastic behaviour of long wavelength modes of light spectator scalar fields during inflation. When starting from a classical field value, the probability distribution for the spectator both spreads out and moves towards an equilibrium distribution. We study the timescales for a mixed quadratic and quartic potential. The timescale of equilibration depends on the parameters of the model, and can be surprisingly large, even much more than thousands of e-folds. These results imply that the initial conditions for spectator fields are not automatically erased during inflation. Applying the results to the curvaton model, we calculate the probability distribution of the curvature perturbation and discuss 'typical' Universes.

Preheating and other particle production phenomena in the early Universe can give rise to high- energy out-of-equilibrium fermions with an anisotropic stress. We develop a formalism to calculate the spectrum of gravitational waves due to fermions, and apply it to a variety of scenarios after inflation. We pay particular attention to regularization issues. We show that fermion production sources a stochastic background of gravitational waves with a significant amplitude, but we find that typical frequencies of this new background are not within the presently accessible direct detec- tion range. However, small-coupling scenarios might still produce a signal observable by planned detectors, and thus open a new window into the physics of the very early Universe.

We address the issue of cosmological backreaction from non-linear structure formation by constructing an approximation for the time evolved metric of a dust dominated universe based on a gradient expansion. Our metric begins as a perturbation of a flat Friedmann-Robertson-Walker state described by a nearly scale invariant, Gaussian, power-law distribution, and evolves in time until non-linear structures have formed. After describing and attempting to control for certain complications in the implementation of this approach, this metric then forms a working model of the universe. We numerically calculate the evolution of the average scale factor in this model and hence the backreaction. We argue that, despite its limitations, this model is more realistic than previous models that have confronted the issue of backreaction. We find that the \emph{instantaneous} effects of backreaction in this model could be as large as $\sim10%$ of the background. This suggests that a proper understanding of the \emph{cumulative} effects of backreaction could be crucial for precision cosmology and any future exploration of the dark sector.

We study the scale dependence of the non-linearity parameters f_NL and g_NL in curvaton models with self-interactions. We show that the spectral indices n_fNL=d ln|f_NL|/(d ln k) and n_gNL=d ln |g_NL|/(d ln k) can take values much greater than the slow--roll parameters and the spectral index of the power spectrum. This means that the scale--dependence of the bi and trispectrum could be easily observable in this scenario with Planck, which would lead to tight additional constraints on the model. Inspite of the highly non-trivial behaviour of f_NL and g_NL in the curvaton models with self-interactions, we find that the model can be falsified if g_NL(k) is also observed.

In an inflationary cosmology, the observed primoridal density perturbation could come from the quantum fluctuations of another light 'curvaton' field, rather than the inflaton. In this case, it is essential that the curvaton decays, converting its perturbation to an adiabatic perturbation. For the first time, we consistently account for this decay in the simplest curvaton model V(\sigma) = (m^2\sigma^2)/2 and point out that it gives rise to an important logarithmic correction to the potential. Moreover, the potential will also receive a correction from the thermal bath. As a consequence, the dynamics of the curvaton are substantially altered compared to the original model in the majority of the parameter space. It will therefore be necessary to re-calculate all the predictions of the original model.

There are observations of 15 high-redshift massive galaxy clusters, which have an extremely small probability with a purely Gaussian initial curvature perturbation. Here we revisit the estimation of the contribution of non-Gaussianities to the cluster mass function and point out serious problems that have resulted in the application of the mass function out of the range of its validity. We remedy the situation and show that the values of f_NL previously claimed to completely reconcile (i.e. at ~100% confidence) the existence of the clusters with LambdaCDM are unphysically small. However, for WMAP cosmology and at 95% confidence, we arrive at the limit f_NL>411, which is similar to previous estimates. We also explore the possibility of a large g_NL as the reason for the observed excess of the massive galaxy clusters. This scenario, g_NL>2*10^6, appears to be in more agreement with CMB and LSS limits for the non-Gaussianity parameters and could also provide an explanation for the overabundance of large voids in the early universe.

We consider a set of scalar fields, consisting of a single flat direction and one or several non-flat directions. We take our cue from the MSSM, considering separately D-flat and F-flat directions, but our results apply to any supersymmetric scenario containing flat directions. We study the field fluctuations during pure de Sitter Inflation, following the evolution of the infrared modes by numerically solving the appropriate Langevin equations. We demonstrate that for the Standard Model U(1), SU(2) or SU(3) gauge couplings, as well as for large enough Yukawa couplings, the fluctuations along the non-flat directions effectively block the fluctuations along the flat directions. The usual expected behaviour <{\phi}^2> \propto N, with N the number of efolds, may be strongly violated, depending on the coupling strengths. As a consequence, those cosmological considerations, which are derived assuming that during inflation flat directions fluctuate freely, should be revised.

We consider inflation within the context of what is arguably the simplest non-metric extension of Einstein gravity. There non-metricity is described by a single graviscalar field with a non-minimal kinetic coupling to the inflaton field $\Psi$, parameterized by a single parameter $\gamma$. We discuss the implications of non-metricity for chaotic inflation and find that it significantly alters the inflaton dynamics for field values $\Psi \gtrsim M_P/\gamma$, dramatically changing the qualitative behaviour in this regime. For potentials with a positive slope non-metricity imposes an upper bound on the possible number of e-folds. For chaotic inflation with a monomial potential, the spectral index and the tensor-to-scalar ratio receive small corrections dependent on the non-metricity parameter. We also argue that significant post-inflationary non-metricity may be generated.

We investigate the scale-dependence of f_NL in the self-interacting curvaton model. We show that the scale-dependence, encoded in the spectral index n_{f_NL}, can be observable by future cosmic microwave background observations, such as CMBpol, in a significant part of the parameter space of the model. We point out that together with information about the trispectrum g_NL, the self-interacting curvaton model parameters could be completely fixed by observations. We also discuss the scale-dependence of g_NL and its implications for the curvaton model, arguing that it could provide a complementary probe in cases where the theoretical value of n_{f_NL} is below observational sensitivity.

The evolution of the curvature perturbation is highly non-trivial for curvaton models with self-interactions and is very sensitive to the parameter values. The final perturbation depends also on the curvaton decay rate $\Gamma$. As a consequence, non-gaussianities can be greatly different from the purely quadratic case, even if the deviation is very small. Here we consider a class of polynomial curvaton potentials and discuss the dynamical behavior of the curvature perturbation. We point out that, for example, it is possible that the non-gaussianity parameter $\fnl\simeq 0$ while $\gnl$ is non-zero. In the case of a curvaton with mass $m\sim {\cal O}(1)$ TeV we show that one cannot ignore non-quadratic terms in the potential, and that only a self-interaction of the type $V_{\rm int}=\sigma^8/M^4$ is consistent with various theoretical and observational constraints. Moreover, the curvaton decay rate should then be in the range $\Gamma=10^{-15}- 10^{-17}$ GeV.

We consider a novel class of $f(\R)$ gravity theories where the connection is related to the conformally scaled metric $\hat g_{\mu\nu}=C(\R)g_{\mu\nu}$ with a scaling that depends on the scalar curvature $\R$ only. We call them C-theories and show that the Einstein and Palatini gravities can be obtained as special limits. In addition, C-theories include completely new physically distinct gravity theories even when $f(\R)=\R$. With nonlinear $f(\R)$, C-theories interpolate and extrapolate the Einstein and Palatini cases and may avoid some of their conceptual and observational problems. We further show that C-theories have a scalar-tensor formulation, which in some special cases reduces to simple Brans-Dicke-type gravity. If matter fields couple to the connection, the conservation laws in C-theories are modified. The stability of perturbations about flat space is determined by a simple condition on the lagrangian.

We construct a model for delayed electroweak symmetry breaking that takes place in a cold Universe with T<<100 GeV and which proceeds by a fast quench rather than by a conventional, slow, phase transition. This is achieved by coupling the Standard Model Higgs to an additional scalar field. We show that the quench transition can be made fast enough for successful Cold Electroweak Baryogenesis, while leaving known electroweak physics unchanged.

Cosmic structures at small non-linear scales $k>L\sim 0.2 h $ Mpc$^{-1}$ have an impact on the longer (quasi-)linear wavelengths with $k<L$ via non-linear UV-IR mode coupling. We evaluate this effect for a $\Lambda$CDM universe applying the effective fluid method of Baumann, Nicolis, Senatore and Zaldarriaga. For $k<L$ the $\Lambda$CDM growth function for the density contrast is found to receive a scale dependent correction and an effective anisotropic stress sources a shift between the two gravitational potentials, setting $\phi$ - $\psi \neq 0$. Since such a situation is generically considered as a signature of modified gravity and/or dark energy, these effects should be taken into account before any conclusions on the dark sector are drawn from the interpretation of future observations.

We consider the constraints for a curvaton with mass m ~ 1 TeV and show that they are not consistent with a purely quadratic potential. Even if the curvaton self-interactions were very weak, they must be accounted for as they affect the dynamical evolution of the curvature perturbation. We show that the only TeV-mass curvaton interaction potential that yields the correct perturbation amplitude, decays before the dark matter freeze-out, and does not give rise to non-Gaussian perturbations that are in conflict with the present limits, is given by V_int= sigma^8/M^4. The decay width of the curvaton should be in the range Gamma= 10^-15...10^-17 GeV. The model typically predicts large non-linearity parameters f_NL and g_NL that should be observable by the Planck satellite. We also discuss various physical possibilities to obtain the required small curvaton decay rate.

We present a systematic study of the amplitude of the primordial perturbation in curvaton models with self-interactions, treating both renormalizable and non-renormalizable interactions. In particular, we consider the possibility that the curvaton energy density is subdominant at the time of the curvaton decay. We find that large regions in the parameter space give rise to the observed amplitude of primordial perturbation even for non-renormalizable curvaton potentials, for which the curvaton energy density dilutes fast. At the time of its decay, the curvaton energy density may typically be subdominant by a relative factor of 10^-3 and still produce the observed perturbation. Field dynamics turns out to be highly non-trivial, and for non-renormalizable potentials and certain regions of the parameter space we observe a non-monotonous relation between the final curvature perturbation and the initial curvaton value. In those cases, the time evolution of the primordial perturbation also displays an oscillatory behaviour before the curvaton decay.

We investigate non-Gaussianities in self-interacting curvaton models treating both renormalizable and non-renormalizable polynomial interactions. We scan the parameter space systematically and compute numerically the non-linearity parameters f_NL and g_NL. We find that even in the interaction dominated regime there are large regions consistent with current observable bounds. Whenever the interactions dominate, we discover significant deviations from the relations f_NL ~ 1/r_decay and g_NL ~ 1/r_decay valid for quadratic curvaton potentials, where r_decay measures the curvaton contribution to the total energy density at the time of its decay. Even if r_decay << 1, there always exists regions with f_NL ~ 0 since the sign of f_NL oscillates as a function of the parameters. While g_NL can also change sign, typically g_NL is non-zero in the low-f_NL regions. Hence, for some parameters the non-Gaussian statistics is dominated by g_NL rather than by f_NL. Due to self-interactions, both the relative signs of f_NL and g_NL and the functional relation between them is typically modified from the quadratic case, offering a possible experimental test of the curvaton interactions.

We propose to address the fine tuning problem of inflection point inflation by the addition of extra vacuum energy that is present during inflation but disappears afterwards. We show that in such a case, the required amount of fine tuning is greatly reduced. We suggest that the extra vacuum energy can be associated with an earlier phase transition and provide a simple model, based on extending the SM gauge group to SU(3)_C \times SU(2)_L\times U(1)_Y\times U(1)_{B-L}, where the Higgs field of U(1)_{B-L} is in a false vacuum during inflation. In this case, there is virtually no fine tuning of the soft SUSY breaking parameters of the flat direction which serves as the inflaton. However, the absence of radiative corrections which would spoil the flatness of the inflaton potential requires that the U(1)_{B-L} gauge coupling should be small with g_{B-L}\leq 10^{-4}.

We show that a spherically symmetric perturbation of a dust dominated $\Omega=1$ FRW universe in the Newtonian gauge can lead to an apparent acceleration of standard candles and provide a fit to the magnitude-redshift relation inferred from the supernovae data, while the perturbation in the gravitational potential remains small at all scales. We also demonstrate that the supernovae data does not necessarily imply the presence of some additional non-perturbative contribution by showing that any Lemaitre-Tolman-Bondi model fitting the supernovae data (with appropriate initial conditions) will be equivalent to a perturbed FRW spacetime along the past light cone.

We investigate how the background evolution affects the curvature perturbations generated by the curvaton, assuming a curvaton potential that may deviate slightly from the quadratic one, and parameterizing the background fluid density as \rho\propto a^{-\alpha}, where a is the scale factor, and \alpha depends on the background fluid. It turns out that the more there is deviation from the quadratic case, the more pronounced is the dependence of the curvature perturbation on \alpha. We also show that the background can have a significant effect on the nonlinearity parameters f_NL and g_NL. As an example, if at the onset of the curvaton oscillation there is a dimension 6 contribution to the potential at 5 % level and the energy fraction of the curvaton to the total one at the time of its decay is at 1 %, we find variations \Delta f_NL \sim \mathcal{O}(10) and \Delta g_NL \sim \mathcal{O}(10^4) between matter and radiation dominated backgrounds. Moreover, we demonstrate that there is a relation between f_NL and g_NL that can be used to probe the form of the curvaton potential and the equation of state of the background fluid.

We argue that the curvaton decay takes place most naturally by way of a broad parametric resonance. The mechanism is analogous to resonant inflaton decay but does not require any tuning of the curvaton coupling strength to other scalar fields. For low scale inflation and a correspondingly low mass scale for the curvaton, we speculate on observable consequences including the possibility of stochastic gravitational waves

We discuss the signatures of non-Gaussianity in the curvaton model where the potential includes also a non-quadratic term. In such a case the non-linearity parameter f_NL can become very small, and we show that non-Gaussianity is then encoded in the non-reducible non-linearity parameter g_NL of the trispectrum, which can be very large. Thus the place to look for the non-Gaussianity in the curvaton model may be the trispectrum rather than the bispectrum. We also show that g_NL measures directly the deviation of the curvaton potential from the purely quadratic form. While g_NL depends on the strength of the non-quadratic terms relative to the quadratic one, we find that for reasonable cases roughly g_NL\sim O(-10^4)-O(-10^5), which are values that may well be accessible by future observations.

We discuss the infrared divergences that appear to plague cosmological perturbation theory. We show that within the stochastic framework they are regulated by eternal inflation so that the theory predicts finite fluctuations. Using the $\Delta N$ formalism to one loop, we demonstrate that the infrared modes can be absorbed into additive constants and the coefficients of the diagrammatic expansion for the connected parts of two and three-point functions of the curvature perturbation. As a result, the use of any infrared cutoff below the scale of eternal inflation is permitted, provided that the background fields are appropriately redefined. The natural choice for the infrared cutoff would of course be the present horizon; other choices manifest themselves in the running of the correlators. We also demonstrate that it is possible to define observables that are renormalization group invariant. As an example, we derive a non-perturbative, infrared finite and renormalization point independent relation between the two-point correlators of the curvature perturbation for the case of the free single field.

We consider the supergravity origin of the recently proposed MSSM inflationary model, which relies on the existence of a saddle point along a dimension six flat direction. We derive the conditions that the Kahler potential has to satisfy for the saddle point to exist irrespective of the hidden sector vevs. We show that these conditions are satisfied by a simple class of Kahler potentials, which we find to have a similar form as in various string theory compactifications. For these potentials, slow roll MSSM inflation requires no fine tuning of the soft supersymmetry breaking parameters.

We model the essential features of eternal inflation on the landscape of a dense discretuum of vacua by the potential $V(\phi)=V_{0}+\delta V(\phi)$, where $|\delta V(\phi)|\ll V_{0}$ is random. We find that the diffusion of the distribution function $\rho(\phi,t)$ of the inflaton expectation value in different Hubble patches may be suppressed due to the effect analogous to the Anderson localization in disordered quantum systems. At $t \to \infty$ only the localized part of the distribution function $\rho (\phi, t)$ survives which leads to dynamical selection principle on the landscape. The probability to measure any but a small value of the cosmological constant in a given Hubble patch on the landscape is exponentially suppressed at $t\to \infty$.

We investigate decaying D-branes as the origin of the thermal string gas of string gas cosmology. We consider initial configurations of low-dimensional branes and argue that they can time evolve to thermal string gas. We find that there is a range in the weak string coupling and fast brane decay time regimes, where the initial configuration could drive the evolution of the dilaton to values, where exactly three spacelike directions grow large.

I discuss the spherically symmetric but inhomogeneous Lemaitre-Tolman- Bondi (LTB) metric, which provides an exact toy model for an inhomogeneous universe. Since we observe light rays from the past light cone, not the expansion of the universe, spatial variation in matter density and Hubble rate can have the same effect on redshift as acceleration in a perfectly homogeneous universe. As a consequence, a simple spatial variation in the Hubble rate can account for the distant supernova data in a dust universe without any dark energy. I also review various attempts towards a semirealistic description of the universe based on the LTB model.

We propose a novel mechanism for dark energy, based on an extended seesaw for scalar fields, which does not require any new physics at energies below the TeV scale. A very light quintessence mass is usually considered to be technically unnatural, unless it is protected by some symmetry broken at the new very light scale. We propose that one can use an extended seesaw mechanism to construct technically natural models for very light fields, protected by SUSY softly broken above a TeV.

We consider an inhomogeneous but spherically symmetric Lemaitre-Tolman-Bondi model to demonstrate that spatial variations of the expansion rate can have a significant effect on the cosmological supernova observations. A model with no dark energy but a local Hubble parameter about 15% larger than its global value fits the supernova data better than the homogeneous model with the cosmological constant. The goodness of the fit is not sensitive to inhomogeneities in the present-day matter density, and our best fit model has Omega_M ~ 0.3, in agreement with galaxy surveys. We also compute the averaged expansion rate, defined by the Buchert equations, of the best fit model and show explicitly that there is no average acceleration.

We present an approach to cosmological perturbations based on a covariant perturbative expansion between two worldlines in the real inhomogeneous universe. As an application, at an arbitrary order we define an exact scalar quantity which describes the inhomogeneities in the number of e-folds on uniform density hypersurfaces and which is conserved on all scales for a barotropic ideal fluid. We derive a compact form for its conservation equation at all orders and assign it a simple physical interpretation. To make a comparison with the standard perturbation theory, we develop a method to construct gauge-invariant quantities in a coordinate system at arbitrary order, which we apply to derive the form of the n-th order perturbation in the number of e-folds on uniform density hypersurfaces and its exact evolution equation. On large scales, this provides the gauge-invariant expression for the curvature perturbation on uniform density hypersurfaces and its evolution equation at any order.

We consider inflaton couplings to MSSM flat directions and the thermalization of the inflaton decay products, taking into account gauge symmetry breaking due to flat direction condensates. We then search for a suitable curvaton candidate among the flat directions, requiring an early thermally induced start for the flat direction oscillations to facilitate the necessary curvaton energy density dominance. We demonstrate that the supersymmetry breaking $A$-term is crucial for achieving a successful curvaton scenario. Among the many possible candidates, we identify the ${\bf u_1dd}$ flat direction as a viable MSSM curvaton.

We argue that all the necessary ingredients for successful inflation are present in the flat directions of the Minimally Supersymmetric Standard Model. We show that out of many gauge invariant combinations of squarks, sleptons and Higgses, there are two directions, ${\bf LLe}$, and ${\bf udd}$, which are promising candidates for the inflaton. The model predicts more than $10^3$ e-foldings with an inflationary scale of $H_{\rm inf}\sim {\cal O}(1-10)$ GeV, provides a tilted spectrum with an amplitude of $\delta_H\sim 10^{-5}$ and a negligible tensor perturbation. The temperature of the thermalized plasma could be as low as $T_{rh}\sim {\cal O}(1-10)$~TeV. Parts of the inflaton potential can be determined independently of cosmology by future particle physics experiments.

We consider curvaton models with potentials that depart slightly from the quadratic form. We show that although such a small departure does not modify significantly the Gaussian part of the curvature perturbation, it can have a pronounced effect on the level of non-Gaussianity. We find that unlike in the quadratic case, the limit of small non-Gaussianity, $|f_{NL}|\ll1$, is quite possible even with small curvaton energy density $r\ll1$ . Furthermore, non-Gaussianity does not imply any strict bounds on $r$ but the bounds depend on the assumptions about the higher order terms in the curvaton potential.

We consider one-loop thermal effects on the decay of a scalar field zero mode initially dominating the energy density of the universe. We assume fermionic decay channels and take into account the effects due to both particle and hole excitations, and present approximate expressions for the absorption and decay rates. We apply the results to the inflaton and solve the Boltzmann equations to find the temperature evolution of the fermionic plasma. We show that the reheat temperature can be greater than the inflaton mass if the inflaton decays into more than one fermionic species.

We consider a two-field model for inflation where the second order metric perturbations can be amplified by a parametric resonance during preheating. We demonstrate that there can arise a considerable enhancement of non-Gaussianity sourced by the local terms generated through the coupled perturbations. We argue that the non-Gaussianity parameter could be as large as f_{NL} ~ 50. Our results may provide a useful test of preheating in future CMB experiments.

We study non-Gaussianity in two distinct models of preheating: instant and tachyonic. In instant preheating non-Gaussianity is sourced by the local terms generated through the coupled perturbations of the two scalar fields. We find that the non-Gaussianity parameter is given by $f_{NL}^{\phi}\sim 2g < O(1)$, where $g$ is a coupling constant, so that instant preheating is unlikely to be constrained by WMAP or Planck. In the case of tachyonic preheating non-Gaussianity arises solely from the instability of the tachyon matter and is found to be large. We find that for single field inflation the present WMAP data implies a bound $V_{0}^{1/4}/M_{P}\leq 10^{-4}$ on the scale of tachyonic instability. We argue that the tachyonic preheating limits are useful also for string-motivated inflationary models.

We demonstrate that string motivated inflation ending via tachyonic instability leaves a detectable imprint on the cosmic microwave background (CMB) radiation by virtue of the excitation of non-Gaussian gravitational fluctuations. The present WMAP bound on non-Gaussianity is shown to constrain the string scale by $M_S/M_P\leq 10^{-4}$ for string coupling $g_{s}<0.1$, hence improving the existing bounds. If tachyon fluctuations during inflation are not negligible, we find the stringent constraint $g_s\sim 10^{-9}$ for $M_{S}/M_{P}<10^{-3}$. This case may soon be ruled out by the forthcoming CMB non-Gaussinianity bounds.

We consider the angular power spectrum in a finite universe with different boundary conditions and perform a fit to the CMB, LSS and supernova data. A finite universe could be the consequence of a holographic constraint, giving rise to an effective IR cutoff at the future event horizon. In such a model there is a cosmic duality relating the dark energy equation of state and the power spectrum, which shows a suppression and oscillatory behaviour that is found to describe the low l features extremely well. However, much of the discussion here will also apply if we actually live inside an expanding bubble that describes our universe. The best fit to the CMB and LSS data turns out to be better than in the standard Lambda-CDM model, but when combined with the supernova data, the holographic model becomes disfavored. We speculate on the possible implications.

We investigate a possible connection between the suppression of the power at low multipoles in the CMB spectrum and the late time acceleration. We show that, assuming a cosmic IR/UV duality between the UV cutoff and a global infrared cutoff given by the size of the future event horizon, the equation of state of the dark energy can be related to the apparent cutoff in the CMB spectrum. The present limits on the equation of state of dark energy are shown to imply an IR cutoff in the CMB multipole interval of 9>l>8.5.

We consider second order inflationary perturbations in the case of two scalar fields, $\sigma$ and the inflaton $\phi$. We derive an expression for the non-Gaussianity of perturbations and apply the results to hybrid inflation. We isolate the contributions due to $\sigma$ and evaluate the resulting terms to show that $\sigma$-induced non-Gaussianities dominate over inflaton-induced non-Gaussianities when $m^2_{\sigma} \gtrsim \eta H^2$, where $\eta$ is the slow-roll parameter. This may provide a useful test of hybrid inflation in forthcoming CMB experiments.

We demonstrate that the MSSM flat directions can naturally account for the seed magnetic fields in the early Universe. The non-zero vacuum expectation value of an MSSM flat direction condensate provides masses to the gauge fields and thereby breaks conformal invariance. During inflation the condensate receives spatial perturbations and $SU(2) x U(1)_Y$ gauge currents are generated together with (hyper)magnetic fields. When these long wavelength vector perturbations reenter our horizon they give rise to $U(1)_{em}$ magnetic fields with an amplitude of $10^{-30}$ Gauss, as required by the dynamo mechanism.

We argue that a brane world with a warped, infinite extra dimension allows for the inflaton to decay into the bulk so that after inflation, the effective dark energy disappears from our brane. This is achieved by the redshifting of the decay products into infinity of the 5th dimension. As a consequence, all matter and CMB density perturbations could have their origin in the decay of a MSSM flat direction rather than the inflaton. We also discuss a string theoretical model where reheating after inflation may not affect the observable brane.

Although often chosen because of simplicity, a single scalar field does not provide a general parametrization of an MSSM flat direction. We derive a formalism for a class of gauge invariant polynomials which result in a multifield description of the flat directions. In contrast to the single field case, the vanishing of the gauge currents yields an important dynamical constraint in the multifield framework. We consider in detail the example of the H_uL flat direction and study the dynamical evolution during and after inflation. We highlight the differences between the single and the multifield flat directions. We show that in the multifield case the field space has an intrinsic curvature and hence unsuppressed non-minimal kinetic terms for the flat direction scalars arise. Also the phases of the individual components non-trivially right after inflation, charging the components of the condensate and producing an enhanced entropy after the decay of the condensate, which is due to cross-coupling of different lepton flavours in the F term. However, the qualitative features of the single field Affleck-Dine baryogenesis, such as the produced total charge, remain largely unchanged.

We consider the possibility that the dark energy responsible for inflation is deposited into extra dimensions outside of our observable universe. Reheating and all matter can then be obtained from the MSSM flat direction condensate involving the Higgses $H_u$ and $H_d$, which acquires large amplitude by virtue of quantum fluctuations during inflation. The reheat temperature is $T_{RH} \lesssim 10^9$ GeV so that there is no gravitino problem. We find a spectral index $n_s\approx 1$ with a very weak dependence on the Higgs potential.

Relativistic scalar field theories with a conserved global charge Q possess often (meta)stable spherically symmetric soliton solutions, called Q-balls. We elaborate on the perfect formal analogy which exists between Q-balls, and spherically symmetric solitons in certain non-relativistic atomic Bose-Einstein condensates, for which the dominant interatomic interaction can be tuned attractive. In a harmonic trap, present in existing experiments, the Q-ball solution is modified in an essential way. If the trap is significantly prolongated in one direction, however, then genuine solitons do appear, and actual experimental data can be obtained for some of the Q-ball properties studied numerically in the relativistic cosmological context, such as their formation and collisions. We also suggest conditions under which the same cosmologically relevant analogies could be extended to the fully three-dimensional case.

We discuss the possibility of generating adiabatic density perturbations from spatial fluctuations in the inflaton decay rate which are due to quantum fluctuations of light moduli fields coupling to the inflaton. We point out that non-renormalizable operators, which lift the flatness of the moduli potential, play an important role for the density perturbations. In particular, the non-renormalizable terms give rise to a considerable damping of the fluctuations and thereby pose an obstruction to the construction of possible models.

We study in detail the possibility that the flat directions of the Minimal Supersymmetric Standard Model (MSSM) could act as a curvaton and generate the observed adiabatic density perturbations. For that the flat direction energy density has to dominate the Universe at the time when it decays. We point out that this is not possible if the inflaton decays into MSSM degrees of freedom. If the inflaton is completely in the hidden sector, its decay products do not couple to the flat direction, and the flat direction curvaton can dominate the energy density. This requires the absence of a Hubble-induced mass for the curvaton, e.g. by virtue of the Heisenberg symmetry. In the case of hidden radiation, $n=9$ is the only admissible direction; for other hidden equations of state, directions with lower $n$ may also dominate. We show that the MSSM curvaton is further constrained severely by the damping of the fluctuations, and as an example, demonstrate that in no-scale supergravity it would fragment into $Q$ balls rather than decay. Damping of fluctuations can be avoided by an initial condition, which for the $n=9$ direction would require an initial curvaton amplitude of $\sim 10^{-2}M_p$, thereby providing a working example of the MSSM flat direction curvaton.

We review the cosmological implications of the flat directions of the Minimally Supersymmetric Standard Model (MSSM). We describe how field condensates are created along the flat directions because of inflationary fluctuations. The post-inflationary dynamical evolution of the field condensate can charge up the condensate with B or L in a process known as Affleck-Dine baryogenesis. Condensate fluctuations can give rise to both adiabatic and isocurvature density perturbations and could be observable in future cosmic microwave experiments. In many cases the condensate is however not the state of lowest energy but fragments, with many interesting cosmological consequences. Fragmentation is triggered by inflation-induced perturbations and the condensate lumps will eventually form non-topological solitons, known as Q-balls. Their properties depend on how supersymmetry breaking is transmitted to the MSSM; if by gravity, then the Q-balls are semi-stable but long-lived and can be the source of all the baryons and LSP dark matter; if by gauge interactions, the Q-balls can be absolutely stable and form dark matter that can be searched for directly. We also discuss some cosmological applications of generic flat directions and Q-balls in the context of self-interacting dark matter, inflatonic solitons and extra dimensions.

We propose that the inflaton is coupled to ordinary matter only gravitationally and that it decays into a completely hidden sector. In this scenario both baryonic and dark matter originate from the decay of a flat direction of the minimal supersymmetric standard model (MSSM), which is shown to generate the desired adiabatic perturbation spectrum via the curvaton mechanism. The requirement that the energy density along the flat direction dominates over the inflaton decay products fixes the flat direction almost uniquely. The present residual energy density in the hidden sector is typically shown to be small.

We describe a new mechanism for reheating the Universe through evaporation of a surface charge of a fragmented inflaton condensate. We show that for a range of Yukawa coupling of the inflaton to the matter sector evaporation gives rise to a much smaller reheat temperature compared to the standard perturbative decay. As a consequence, reheating through a surface effect could solve the gravitino and moduli over production problem in inflationary models without fine tuning the Yukawa sector.

We consider the pre-big bang scenario with a massive axion field which starts to dominate energy density when oscillating in an instanton-induced potential and subsequently reheats the universe as it decays into photons, thus creating adiabatic CMB perturbations. We find that the fluctuations in the axion field can give rise to a nearly flat spectrum of adiabatic perturbations with a spectral tilt $\Delta n$ in the range $-0.1 \lesssim \Delta n \lesssim 1$.

We consider pure isocurvature cold dark matter models in the case of open and closed universe. We allow for a large spectral tilt and scan the 6-dimensional parameter space for the best fit to the COBE, Boomerang, and Maxima-1 data. Taking into account constraints from large-scale structure and big bang nucleosynthesis, we find a best fit with $\chi^2 = 121$, which is to be compared to $\chi^2 = 44$ of a flat adiabatic reference model. Hence the current data strongly disfavour pure isocurvature perturbations.

We consider different types of Q-balls as self-interacting dark matter. For the Q-balls to act as the dark matter of the universe they should not evaporate, which requires them to carry very large charges; depending on the type, the minimum charge could be as high as Q \sim 10^{33} or the Q-ball coupling to ordinary matter as small as \sim 10^{-35}. The cross-section-to-mass ratio needed for self-interacting dark matter implies a mass scale of m \sim O(1) MeV for the quanta that the Q-balls consist of, which is very difficult to achieve in the MSSM.

We consider heavy stable neutral particles in the context of supergravity and show that a gravitationally suppressed inflaton decay can produce such particles in cosmologically interesting abundances within a wide mass range $10^3 {\rm GeV} \leq m_X \leq 10^{11} {\rm GeV}$. In gravity-mediated supersymmetry breaking models, a heavy particle can decay into its superpartner and a photon-photino pair or a gravitino. Such decays only change the identity of a possible dark matter candidate. However, for $10^3 {\rm GeV} \leq m_X \leq 10^7 {\rm GeV}$, astrophysical bounds from gamma-ray background and photodissociation of light elements can be more stringent than the overclosure bound, thus ruling out the particle as a dark matter candidate.

We describe a simple and predictive scenario for baryogenesis in theories with large extra dimensions which resembles Affleck-Dine baryogenesis. The Affleck-Dine field is a complex scalar field carrying a $U(1)_{\chi}$ charge which is dynamically broken after the end of inflation. This generates an excess of $\chi$ over $\bar \chi$, which then decays into Standard Model fermions to produce an excess of baryons over anti-baryons. Our model is very constrained because the Affleck-Dine field has to be sufficiently flat during inflation. It is also a source for density fluctuations which can be tested in the coming satellite and balloon experiments.

We consider active-sterile neutrino oscillations in the early universe in an inhomogeneous isocurvature background. We show that very small initial baryonic seed-inhomogeneities can trigger a growth of very large amplitude spatial fluctuations in lepton asymmetry. Domains of varying asymmetry are observed to persist for a long time despite dissipative effects. The space dependent asymmetry profiles give rise to MSW-resonances within the domain boundaries, enhancing dramatically the equilibration of the sterile neutrino species. According to our one-dimensional toy-model, the effect is so strong that almost the entire parameter space where exponential growth of asymmetry occurs would be ruled out by nucleosynthesis.

We present numerical simulations of fragmentation of the Affleck-Dine condensate in two spatial dimensions. We argue analytically that the final state should consist of both Q-balls and anti-Q-balls in a state of maximum entropy, with most of the balls small and relativistic. Such a behaviour is found in simulations on a 100x100 lattice with cosmologically realistic parameter values. During fragmentation process, we observe filament-like texture in the spatial distribution of charge. The total charge in Q-balls is found to be almost equal to the charge in anti-Q-balls and typically orders of magnitude larger than charge asymmetry. Analytical considerations indicate that, apart from geometrical factors, the results of the simulated two dimensional case should apply also to the fully realistic three dimensional case.

We present the constraints on isocurvature fluctuations for a flat universe implied by the Boomerang and Maxima-1 data on the anisotropy of the cosmic microwave background. Because the new data defines the shape of the angular power spectrum in the region of the first acoustic peaks much more clearly than earlier data, even a tilted pure isocurvature model is now ruled out. However, a mixed model with a sizable isocurvature contribution remains allowed. We consider primordial fluctuations with different spectral indices for the adiabatic and isocurvature perturbations, and find that the 95% C.L. upper limit to the isocurvature contribution to the low multipoles is $\alpha \leq 0.63$. The upper limit to the contribution in the ${\l}\sim 200$ region is $\alpha_{200} \leq 0.13$.

We consider the evolution of the Affleck-Dine scalar during D-term and F-term inflation and solve the combined slow-roll equations of motion. We show that for a typical case, where both the Affleck-Dine scalar and inflaton initially have large values, in D-term inflation the Affleck-Dine scalar is driven to a fixed value, with only a very slight dependence on the number of e-foldings. As a result, there is a definite prediction for the ratio of the baryonic isocurvature perturbation to the adiabatic perturbation. In minimal (d=4) Affleck-Dine baryogenesis the relative isocurvature contribution to the CMB angular power spectrum amplitude is predicted to be in the range $0.01-0.1$, which can account for present large-scale structure observations and should be observable by PLANCK. In a very general case, scale-invariance of the adiabatic perturbations from the Affleck-Dine scalar imposes a lower bound of about 0.01 for d=4. For d=6 the isocurvature perturbation may just be observable, although this is less certain. We also consider F-term inflation and show that the magnitude of the baryonic isocurvature perturbation is fixed by the value of $H$ during inflation. For typical values of $H$ the isocurvature perturbation could be close to present observational limits.

We consider the detection possibilities of isocurvature fluctuations in the future CMB satellite experiments MAP and Planck for different cosmological reference models. We present a simultaneous 10 parameter fit (8 for the case of open model) to determine the correlations between the cosmological parameters, including isocurvature cold dark matter contribution to the anisotropy. Assuming that polarization information can be fully exploited, we find that an isocurvature perturbation can be detected by the Planck Surveyor if the ratio of the initial isocurvature and adiabatic perturbation amplitude is larger than 0.07. In the absence of polarization data, the signal from isocurvature perturbations can be confused with tensor perturbations or early reionization effects, and the limit is larger by almost an order of magnitude.

We consider numerically the growth of neutrino asymmetry in active-sterile neutrino oscillations in the early universe. It is shown that the final sign of the asymmetry can be highly sensitive to small variations of the oscillation parameters. We find regions which are completely or partially chaotic, but also regions where the sign remains very robust. The consequences for atmospheric neutrino oscillations and primordial nucleosynthesis are then discussed. In the completely chaotic region the predicted He(4)-abundance has an inherent arbitrariness \Delta Y \simeq 0.01.

We consider the B-ball cosmology of the MSSM in the context of D-term inflation models where the reheating temperature is determined by the Affleck-Dine mechanism to be of the order of 1 GeV. We show that such a low reheating temperature can arise quite naturally as the result of a symmetry which is required to maintain the flatness of the inflaton potential. In this case the B-balls will decay well after the freeze-out temperature of the LSP, allowing baryons and cold dark matter to originate primarily from B-ball decays.

In D-term inflation models, Affleck-Dine baryogenesis produces isocurvature density fluctuations. These can be perturbations in the baryon number, or, in the case where the present neutralino density comes directly from B-ball decay, perturbations in the number of dark matter neutralinos. The latter case results in a large enhancement of the isocurvature perturbation. The requirement that the deviation of the adiabatic perturbations from scale invariance due to the Affleck-Dine field is not too large then imposes a lower bound on the magnitude of the isocurvature fluctuation of about $10^{-2}$ times the adiabatic perturbation. This should be observable by MAP and PLANCK.

We demonstrate that B-ball decay in the MSSM can naturally solve the puzzle of why the densities of baryons and dark matter in the Universe are similar. This requires that the B-balls survive thermalization and decay below the freeze-out temperature of the neutralino LSP, typically 1-10GeV. It is shown that this can happen if the baryon asymmetry originates from a squark condensate along the d=6 u^{c}d^{c}d^{c} D-flat direction of the MSSM scalar potential. For this to work the reheating temperature after inflation must be no greater than 10^3-10^5 GeV.

We consider the effect of the presence of a hypermagnetic field at the electroweak phase transition. Screening of the Z-component inside a bubble of the broken phase delays the phase transition and makes it stronger first order. We show that the sphaleron constraint can be evaded for m_H up to 100 GeV if a B_Y\gsim 0.3 T^2 exists at the time of the EW phase transition, thus resurrecting the possibility for baryogenesis within the minimal standard model (provided enough CP violation can be obtained). We estimate that for m_H\gsim 100 GeV the Higgs condendsate behaves like a type II superconductor with Z-vortices penetrating the bubble. Also, for such high Higgs masses the minimum B_Y field required for a strong first order phase transition is large enough to render the W-field unstable towards forming a condensate which changes the simple picture of the symmetry breaking.

We propose to test the magnetic transition moments of Majorana neutrinos by comparing the fluxes of different flavours of neutrinos coming from active galactic nuclei (AGN). We show that, with reasonable assumptions about the magnetic field of the AGN, it is possible to obtain limits on $\nu_{\tau}\nu_{e}$ and $\nu_{\tau}\nu_{\mu}$ transition moments which are three to five orders of magnitude better than the laboratory limits. We also point out that with certain parameter values the ratio $\nu_{\tau}/\nu_{e,\mu}$, when measured from different sources, is expected to vary from zero to values somewhat higher than one, providing an unambigious signal of a magnetic transition within the AGN which cannot be explained by neutrino oscillations.

The explanation of the observed galactic magnetic fields may require the existence of a primordial magnetic field. Such a field may arise during the early cosmological phase transitions, or because of other particle physics related phenomena in the very early universe reviewed here. The turbulent evolution of the initial, randomly fluctuating microscopic field to a large-scale macroscopic field can be described in terms of a shell model, which provides an approximation to the complete magnetohydrodynamics. The results indicate that there is an inverse cascade of magnetic energy whereby the coherence of the magnetic field is increased by many orders of magnitude. Cosmological seed fields roughly of the order of $10^{-20}$ G at the scale of protogalaxy, as required by the dynamo explanation of galactic magnetic fields, thus seem plausible.

We take the first steps towards the complete computation of the thermalization rate of the supersymmetric particles involved in electroweak baryogenesis by computing the thermalization rate of the right-handed stop from the imainary part of the two-point Green function. We use improved propagators including resummation of hard thermal loops. The thermalization rate is computed at the one-loop level in the high temperature approximation as a function of $M_{\tilde t_R}(T)$. We also give an estimate for the magnitude of the two-loop contributions which dominate the rate for small $M_{\tilde t_R}(T)$. If the stop is non-relativistic with $M_{\tilde t_R}(T)\gg T$, thermalization takes place by decay and is very fast.

The observed galactic magnetic fields may have a primordial origin. I briefly review the observations, their interpretation in terms of the dynamo theory, and the current limits on cosmological magnetic fields. Several possible mechanisms for generating a primordial magnetic field are then discussed. Turbulence and the evolution of the microscopic fields to macroscopic fields is described in terms of a shell model, which provides an approximation to the full magnetohydrodynamics and indicates the existence of an inverse cascade of magnetic energy. Cosmological seed fields roughly of the order of $10^{-20}$ G at the scale of protogalaxy, as required by the dynamo explanation of galactic magnetic fields, seem rather plausible.

Neutrinos with a magnetic moment $\mu$ change their helicity when interacting with an electromagnetic field. Various aspects of this effect have been described as spin precession, spin-flip scattering, and magnetic Cherenkov radiation. These perspectives are unified in an expression for the $\nu_L\to\nu_R$ transition rate which involves the correlators of the electromagnetic field distribution. Our general formula corrects a previous result and generalizes it to the case where the fields cannot be viewed as classical and where the momentum transfers need not be small. We evaluate our result explicitly for a relativistic QED plasma and determine the depolarization rate to leading order in the fine structure constant. Assuming that big-bang nucleosynthesis constraints do not allow a right-handed neutrino in equilibrium we derive the limit $\mu<6.2\EE{-11}\mu_B$ on the neutrino magnetic moment. Bounds on $\mu$ from a possible large scale magnetic fields are found to be more stringent even for very weak fields

We consider the formation of a ring-like magnetic field in collisions of bubbles of broken phase in an abelian Higgs model. Particular attention is paid on multiple collisions. The small collision velocity limit, appropriate to the electroweak phase transition, is discussed. We argue that after the completion of the electroweak phase transition, when averaged over nucleation center distances, there exists a mean magnetic field $B\simeq 2.0\times 10^{20}$ G with a coherence length $9.1 \times 10^3 GeV^{-1}$ (for m_H=68 GeV). Because of the ring-like nature of B, the volume average behaves as $B \sim 1/L$. Taking into account the turbulent enhancement of the field by inverse cascade, we estimate that colliding electroweak bubbles would give rise to a mean field $B_{rms}\simeq 10^{-21}$ G at 10 Mpc comoving scale today.

We study the effects of plasma viscosity on the dynamics of primordial magnetic fields by simulating magnetohydrodynamics in the early universe by appropriate non-linear cascade models. We find numerically that even in the presence of large kinetic viscosity, magnetic energy is transferred to large length scales. There are indications, however, that the inverse cascade stops at a given time which depends on the magnitude of viscosity. For realistic viscosities we do not find equipartition between magnetic and kinetic energies.

We study the equilibration of the right-helicity states of light Dirac neutrinos in the early universe by solving the momentum dependent Boltzmann equations numerically. We show that the main effect is due to electroweak gauge boson poles, which enhance thermalization rates by some three orders of magnitude. The right-helicity states of tau neutrinos will be brought in equilibrium independently of their initial distribution at a temperature above the poles if the tau neutrino mass is larger than about 10 keV.

We investigate hydromagnetic turbulence of primordial magnetic fields using magnetohydrodynamics (MHD) in an expanding universe. We present the basic, covariant MHD equations, find solutions for MHD waves in the early universe, and investigate the equations numerically for random magnetic fields in two spatial dimensions. We find the formation of magnetic structures at larger and larger scales as time goes on. In three dimensions we use a cascade (shell) model, that has been rather successful in the study of certain aspects of hydrodynamic turbulence. Using such a model we find that after ${\cal O}(10^9)$ times the initial time the scale of the magnetic field fluctuation (in the comoving frame) has increased by 4-5 orders of magnitude as a consequence of an inverse cascade effect (i.e. transfer of energy from smaller to larger scales). Thus {\it at large scales} primordial magnetic fields are considerably stronger than expected from considerations which do not take into account the effects of MHD turbulence.

In the presence of primordial magnetic fields the oscillating cosmic axion field drives an oscillating electric field. The ensuing dissipation of axions is found to be inversely proportional to the conductivity of the primordial plasma. This counterintuitive result is essentially equivalent to ``Zeno's paradox'' or the ``watched-pot effect'' of quantum mechanics. It implies that the standard predictions of the cosmic axion density remain unaltered even if primordial magnetic fields are strong.