We provide a detailed study of reheating in the kination regime where particle content is created by gravitational production of massive scalars mutually interacting with a massless scalar field. The produced particles subsequently decay into massless particles eventually reheating the Universe. We aim for a more precise picture using Boltzmann equations and decay rates obtained by methods of quantum field theory in curved spacetime. By numerical calculations it is found that after inflation the Universe ends up being dominated by ordinary matter for a while before the radiation dominated era. The reheating temperature itself is found to be in the $10^6-10^{12}\ \text{GeV}$ regime.

We derive a curved space generalization of a scalar to fermion decay rate with a Yukawa coupling in expanding Friedmann-Robertson-Walker universes. This is done using the full theory of quantum fields in curved spacetime and the added-up transition probability method. It is found that in an expanding universe the usual Minkowskian decay rates are considerably modified for early times. For conformally coupled scalars the decay rate is modified by a positive additive term proportional to the inverse of mass and related to the expansion rate of the Universe. We compare and contrast our results with previous studies on scalar to scalar decay and find that in general the decay channel into fermions is the dominant channel of decay in the very early Universe.

The lack of energy conservation introduces new particle processes in curved spacetime that are forbidden in flat space. Therefore one has to be very cautious about using the results calculated in Minkowskian space in early universe applications. This is true for particle decay rates in particular, which need to be calculated using quantum field theory in curved spacetime. Previous studies are usually restricted to using minimal or conformal coupling for the decaying particle, while using a more general coupling would give deeper insight into particle decay. This paper presents the results we obtained for a massive particle decaying in a general power-law universe with arbitrary coupling to gravity. We find that depending on the value of the gravitational coupling, the effect of gravitation may either strengthen or weaken the decay. The analysis further reveals that, apart from radiation dominated universe, there are values of the coupling constant for which the decay rate is exactly Minkowskian for all universe types. Because the decay rate may be considerably modified in curved space, these issues need to be considered when doing precise cosmological calculations.

According to the standard model of cosmology, the early universe has been dominated by radiation or non-relativistic matter in several eras of its history. However, many cosmological calculations involving particle processes are commonly done using Minkowskian results for them, although, for more precise treatment, quantum field theory in curved spacetime is needed. This paper aims to fill this gap by presenting decay rates for matter and radiation dominated universes in this more precise treatment. We provide a study of the average decay rates for a process where a conformally coupled massive scalar field decays into massless scalar particles. It is found that the presence of a curved spacetime modifies the Minkowskian result considerably for early times but asymptotically only by an additive term proportional to the inverse of mass and interaction time. Thus, the correction is small for large time scales, but on the time scales of the order of $m\sim t$, the relative correction term may be of importance.

In the presence of a gravitational field decay rates may significantly differ from flat space equivalent. By studying mutually interacting quantum fields the decay rates can be calculated on a given spacetime. This paper presents the calculation of the transition probability for the decay of a massive scalar particle in a stiff matter dominated universe. We find that due to the precence of a gravitational field a finite correction to the transition probability is added which depends inversely on the mass. Moreover the decay rate is smaller and lifetime of the particles is longer compared to flat space. The mass dependence is such that the lifetime of lighter particles is prolonged more compared to heavier particles. This result may be of significance when studying cosmological situations involving stiff matter.

A scenario for gravitational particle creation in a stiff matter dominated flat Friedmann-Robertson-Walker universe is presented. The primary creation of scalar particles is calculated using quantum field theory in curved spacetime and it is found to be strongly dependent on the scalar mass and the expansion parameter of the universe. The particle creation is most effective for a very massive scalar field and large expansion parameter. We apply the results to cosmology and calculate an upper bound for the equilibrium temperature of the secondary particles created by the scalar field decay.

We examine the effects of $f(R)$ gravity on Jeans analysis of collapsing dust clouds. We provide a method for testing modified gravity models by their effects on star formation as the presence of $f(R)$ gravity is found to modify the limit for collapse. In this analysis we add perturbations to a de Sitter background. Depending on the characteristics of a chosen $f(R)$ model, the appearance of new limits is possible. The physicality of these limits is further examined. We find the asymptotic Jeans masses for $f(R)$ theories compared to standard Jeans mass. Through this ratio, the effects of the $f(R)$ modified Jeans mass for viable theories are examined in molecular clouds. Bok globules have a mass range comparable to Jeans masses in question and are therefore used for comparing different $f(R)$ models. Viable theories are found to assist in star formation.

The recent LIGO observation sparked interest in the field of gravity wave signals. Besides the gravity wave observation the LIGO collaboration used the inspiraling black hole pair to constrain the graviton mass. Unlike general relativity, $f(R)$ theories have a characteristic non-zero mass graviton. We apply the constraint on the graviton mass to viable $f(R)$ models to find the effects on model parameters. We find it possible to constrain the parameter space with the gravity wave based observations. We make a case study for the popular Hu-Sawicki model and find a parameter bracket. The result generalizes to other $f(R)$ theories and can be used to contain the parameter space.

The Lema\'itre-Toman-Bondi (LTB) models have reported to suffer from incompatibility with cosmological observations and fine-tuning of the observer's location. Further analysis of these issues indicates that they could be resolved by models that are compatible with the supernova Ia data, but less inhomogeneous than those that have been presented in the literature so far. We study if such models exist by employing the degrees of freedom of the LTB models in a novel manner. We discovered two scenarios which may meet the expectations, but extensive numerical and analytical investigation showed them inviable. We extended our studies to the $\Lambda$LTB models, which generalizes the LTB models by including a non-zero cosmological constant $\Lambda$ in Einsteins equations. This adds an additional degree of freedom for the earlier scenarios and introduces a new scenario capable of meeting the expectations. However, extensive numerical and analytical investigation reveals that inclusion of $\Lambda$ does not enhance the viability of the models. We identify the lack of degrees of freedom to be the reason for the unviability. However, the method presented here can be generalized to models including more degrees of freedom, like the Szekeres models, which have more promise to overcome the issues in the LTB models.

We investigate Lema\^itre-Tolman-Bondi (LTB) models, whose early time evolution and bang time are homogeneous and the distance - redshift relation and local Hubble parameter are inherited from the $\Lambda$CDM model. We show that the obtained LTB models and the $\Lambda$CDM model predict different relative local expansion rates and that the Hubble functions of the models diverge increasingly with redshift. The LTB models show tension between low redshift baryon acoustic oscillation and supernova observations and including Lyman-$\alpha$ forest or cosmic microwave background observations only accentuates the better fit of the $\Lambda$CDM model compared to the LTB model. The result indicates that additional degrees of freedom are needed to explain the observations, for example by renouncing spherical symmetry, homogeneous bang time, negligible effects of pressure, or the early time homogeneity assumption.

We study the scalar field potential $V(\phi)$ in the scalar-tensor gravity with self-consistent polytropic stellar configurations. Without choosing a particular potential, we numerically derive the potential inside various stellar objects. We restrict the potential to conform to general relativity or to $f(R)$ gravity inside and require the solution to arrive at SdS vacuum at the surface. The studied objects are required to obtain observationally valid masses and radii corresponding to solar type stars, white dwarfs and neutron stars. We find that the resulting scalar-tensor potential $V(\phi)$ for the numerically derived polytrope that conforms to general relativity, in each object class, is highly dependent on the matter configuration as well as on the vacuum requirement at the boundary. As a result, every stellar configuration arrives at a potential $V(\phi)$ that is not consistent with the other stellar class potentials. Therefore, a general potential that conforms to all these polytropic stellar classes could not be found.

Inhomogeneous cosmological models are often reported to suffer from a fine-tuning problem because of the observer's location. We study if this is a generic feature in the Lema\^{i}tre-Tolman (LT) models, by investigating if there are models with freedom in the initial state. In these cases, the present fine-tuned location would be evolved from a non-fine-tuned initial state and thus vanishing the problem. In this paper, we show that this is not a generic problem and we give the condition when the LT models do not have fine-tuned initial state. The physical meaning of this condition, however, requires more investigation. We investigate if this condition can be found from a special case: homogeneous models with matter, dark, and curvature density as parameters. We found that with any reasonable density values, these models do not satisfy this condition and thus do not have freedom in the initial state. We interpret this to be linked with the fine-tuning problem of the initial state of the homogeneous models, when the early time inflation is not included to them. We discuss of the condition in the context of non-homogeneous models.

It is stated that a class of $f(R)$ gravity models seem to obtain $\Lambda$CDM behaviour for high redshifts and general relativistic behaviour locally at high curvatures. In the present paper, we numerically study polytropic configurations that resemble stars like young sun with Hu and Sawicki $f(R)$ gravity field equations and compare the spacetime at the boundary to the general relativistic counterpart. These polytropes are stationary spherically symmetric configurations and have regular metrics at the origin. Since Birkhoff's theorem does not apply for modified gravity, the solution outside may deviate from Schwarzschild-de Sitter spacetime. At the boundary, Post-Newtonian parametrization was used to determine how much the studied model deviates from the general relativistic $\Lambda$CDM model.

The evolution of the scale parameter in the Hu-Sawicki model is examined. We search the parameter area for instabilities. It turns out the parameter area of physically meaningful evolution is non-existent. For greater stability a stabilising term such as R^2 should be added. With the addition, the area mainly coincides with phantom models. We also find that the evolution of the scale parameter can be examined through a linear expansion and not only in the Hu-Sawicki case. The oscillatory behaviour found in the evolution raises questions of the physicality of the model. The behaviour can also cause current observations of the Hubble parameter to be misleading. This paper has been withdrawn by the authords due to a critical error in equation (10)

We apply phase space analysis to inhomogeneous cosmological model given by Lema\^itre-Tolman model. We describe some general conditions required to interpret the model stable enough and, in the present paper, apply them to two special cases: dust filled homogeneous model with and without cosmological constant. We find that such stability explaining all present astrophysical observations can not be achieved due to instabilities in phase space. This hints that non-homogeneous models are not likely to be physically viable, although any conclusive analysis requires more realistic modeling of non-homogeneous universe.

The framework for soft leptogenesis minimally extended with a DM sector is studied. A heavy singlet neutrino superfield acts as the source for (s)lepton asymmetry and by coupling to the singlet DM superfield it produces a DM particle density through decays. The nature of DM generated is twofold depending on whether the Yukawa and DM couplings are either small or large. With sufficiently small Yukawa and DM couplings DM annihilations into MSSM particles are slow and as a consequence all DM particles form the DM component. The solutions to Boltzmann equations are given and the dependence between the DM masses and coupling are presented in this weak coupling regime. Also, the behavior of the efficiency of producing asymmetric DM is determined with weak couplings. We note that a different outcome arises if the couplings are larger because then the ADM component is dominant due to the effectiveness of DM decays into the MSSM sector.

We consider an inflationary curvaton scenario, where the curvaton decays into two non-interacting relativistic fluids and later during the cosmological evolution one of them becomes non-relativistic, forming dark matter component of the universe. We study the thermic properties and the generation of non-gaussianity in this three fluid curvaton model. By solving the evolution of the system and using several cosmological conditions we find that the allowed parameter space is strongly constrained. The naturalness of this curvaton scenario is also discussed.

Hamiltonian perturbation theory is used to analyse the stability of f(R) models. The Hamiltonian equations for the metric and its momentum conjugate are written for f(R) Lagrangian in the presence of perfect fluid matter. The perturbations examined are perpendicular to R. As perturbations are added to the metric and momentum conjugate to the induced metric instabilities are found, depending on the form of f(R). Thus the examination of these instabilities is a way to rule out certain f(R) models.

The affine connection in a space-time with a maximally symmetric spatial subspace is derived using the properties of maximally symmetric tensors. The number of degrees of freedom in metric-affine gravity is thereby considerably reduced while the theory allows spatio-temporal torsion and remains non-metric. The Ricci tensor and scalar are calculated in terms of the connection and the field equations derived for the Einstein-Hilbert as wells as for f(R) Lagrangians. By considering specific forms of f(R), we demonstrate that the resulting Friedmann equations in Palatini formalism without torsion and metric-affine formalism with maximal symmetry are in general different in the presence of matter.

The generation of non-gaussianity is studied in a three fluid curvaton model. By utilizing second order perturbation theory we derive general formulae for the large scale temperature fluctuation and non-gaussianity parameter, $f_{NL}$, that includes the possibility of a non-adiabatic final state. In the adiabatic limit we recover previously known results. The results are applied to a three fluid curvaton model where the curvaton decays into radiation and matter. We find that the amount of non-gaussianity decreases as the final state of the system becomes more adiabatic and that the generated non-gaussianity in the scenario is small, $|f_{NL}| \sim \mathcal{O}(1)$.

Despite the fact that modified theories of gravity, in particular the f(R) gravity models have attracted much attention in the last years, the problem of the energy localization in the framework of these models has not been addressed. In the present work the concept of energy-momentum complexes is presented in this context. We generalize the Landau-Lifshitz prescription of calculating the energy-momentum complex to the framework of f(R) gravity. As an important special case, we explicitly calculate the energy-momentum complex for the Schwarzschild-de Sitter metric for a general f(R) theory as well as for a number of specific, popular choices of f(R).

We consider general metric $f(R)$ theories of gravity by solving the field equations in the presence of a spherical static mass distribution by analytical perturbative means. Expanding the field equations systematically in $\cO(G)$, we solve the resulting set of equations and show that $f(R)$ theories which attempt to solve the dark energy problem very generally lead to $\gamma_{PPN}=1/2$ in the solar system. This excludes a large class of theories as possible explanations of dark energy. We also present the first order correction to $\gamma_{PPN}$ and show that it cannot have a significant effect.

A three fluid system describing the decay of the curvaton is studied by numerical and analytical means. We place constraints on the allowed interaction strengths between the fluids and initial curvaton density by requiring that the curvaton decays before nucleosynthesis while nucleosynthesis, radiation-matter equality and decoupling occur at correct temperatures. We find that with a continuous, time-independent interaction, a small initial curvaton density is naturally preferred along with a low reheating temperature. Allowing for a time-dependent interaction, this constraint can be relaxed. In both cases, a purely adiabatic final state can be generated, but not without fine-tuning. Unlike in the two fluid system, the time-dependent interactions are found to have a small effect on the curvature perturbation itself due to the different nature of the system. The presence of non-gaussianity in the model is discussed.

We study stellar configurations and the space-time around them in metric $f(R)$ theories of gravity. In particular, we focus on the polytropic model of the Sun in the $f(R)=R-\mu^4/R$ model. We show how the stellar configuration in the $f(R)$ theory can, by appropriate initial conditions, be selected to be equal to that described by the Lane-Emden -equation and how a simple scaling relation exists between the solutions. We also derive the correct solution analytically near the center of the star in $f(R)$ theory. Previous analytical and numerical results are confirmed, indicating that the space-time around the Sun is incompatible with Solar System constraints on the properties of gravity. Numerical work shows that stellar configurations, with a regular metric at the center, lead to $\gamma_{PPN}\simeq1/2$ outside the star ie. the Schwarzschild-de Sitter -space-time is not the correct vacuum solution for such configurations. Conversely, by selecting the Schwarzschild-de Sitter -metric as the outside solution, we find that the stellar configuration is unchanged but the metric is irregular at the center. The possibility of constructing a $f(R)$ theory compatible with the Solar System experiments and possible new constraints arising from the radius-mass -relation of stellar objects is discussed.

The large-scale dynamics of a two-fluid system with a time dependent interaction is studied analytically and numerically. We show how a rapid transition can significantly suppress the large-scale curvature perturbation and present approximative formulae for estimating the effect. By comparing to numerical results, we study the applicability of the approximation and find good agreement with exact calculations.

Static spherically symmetric perfect fluid solutions are studied in metric $f(R)$ theories of gravity. We show that pressure and density do not uniquely determine $f(R)$ ie. given a matter distribution and an equation state, one cannot determine the functional form of $f(R)$. However, we also show that matching the outside Schwarzschild-de Sitter-metric to the metric inside the mass distribution leads to additional constraints that severely limit the allowed fluid configurations.

Spherically symmetric static empty space solutions are studied in f(R) theories of gravity. We reduce the set of modified Einstein's equations to a single equation and show how one can construct exact solutions in different f(R) models. In particular, we show that for a large class models, including e.g. the f(R)=R-\mu^4/R model, the Schwarzschild-de Sitter metric is an exact solution of the field equations. The significance of these solutions is discussed in light of solar system constraints on $f(R)$ theories of gravity.

A class of supergravity inspired quintessence models is studied by comparing to cosmological data. The set of considered models includes several previously studied quintessential potentials, as well as the $\Lambda$CDM model. We find that even though the commonly studied supergravity inspired quintessence models fit the data better than the $\Lambda$CDM model, they are a relatively poor fit when compared to the best fit model in the studied class. Our results suggest a low energy scale, less than ${\cal{M}}\sim 1$ TeV, for the effective supergravity potential.

Modified theories of gravity have recently been studied by several authors as possibly viable alternatives to the cosmological concordance model. Such theories attempt to explain the accelerating expansion of the universe by changing the theory of gravity, instead of introducing dark energy. In particular, a class of models based on higher order curvature invariants, so-called $f(R)$ gravity models, has drawn attention. In this letter we show that within this framework, the expansion history of the universe does not uniquely determine the form of the gravitational action and it can be radically different from the standard Einstein-Hilbert action. We demonstrate that for any barotropic fluid, there always exists a class of $f(R)$ models that will have exactly the same expansion history as that arising from the Einstein-Hilbert action. We explicitly show how one can extend the Einstein-Hilbert action by constructing a $f(R)$ theory that is equivalent on the classical level. Due to the classical equivalence between $f(R)$ theories and Einstein-Hilbert gravity with an extra scalar field, one can also hence construct equivalent scalar-tensor theories with standard expansion.

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 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.

The evaporation of the black holes during the very early universe is studied. Starting from black hole filled universe, the distiributions of particle species are calculated and showed, that they differ remarkably from the corresponding equilibrium distributions. This may have great impact to the physics of the very early universe. Also the evolution of the universe during the evaporation has been studied.

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.