In the present work we present the full treatment of scalar and vector cosmological perturbations in a non-singular bouncing universe in the context of metric $f(R)$ cosmology. Scalar metric perturbations in $f(R)$ cosmology were previously calculated in the Jordan frame, in the present paper we successfully use the Einstein frame to calculate the scalar metric perturbations where the cosmological bounce takes place in the Jordan frame. The Einstein frame picture presented corrects and completes a previous calculation of scalar perturbations and adds new information. Behavior of fluid velocity potential and the pure vector fluid velocity terms are elaborately calculated for the first time in $f(R)$ cosmology for a bouncing universe in presence of exponential gravity. It is shown that the vector perturbations can remain bounded and almost constant during the non-singular bounce in $f(R)$ gravity unlike general relativistic models where we expect the vector perturbations to be growing during the contracting phase and decaying during the expanding phase. The paper shows that the Einstein frame can be used for calculation of scalar and vector metric perturbations in a bouncing universe for most of the cases except the case of asymmetric non-singular bounces.

In this work we present some cosmologically relevant solutions using the spatially flat Friedmann-Lemaitre-Robertson-Walker (FLRW) spacetime in metric $f(R)$ gravity where the form of the gravitational Lagrangian is given by $\frac{1}{\alpha}e^{\alpha R}$. In the low curvature limit this theory reduces to ordinary Einstein-Hilbert Lagrangian together with a cosmological constant term. Precisely because of this cosmological constant term this theory of gravity is able to support nonsingular bouncing solutions in both matter and vacuum background. Since for this theory of gravity $f^{\prime}$ and $f^{\prime\prime}$ is always positive, this is free of both ghost instability and tachyonic instability. Moreover, because of the existence of the cosmological constant term, this gravity theory also admits a de-Sitter solution. Lastly we hint towards the possibility of a new type of cosmological solution that is possible only in higher derivative theories of gravity like this one.

As our understanding of the past in a bouncing universe is limited, it becomes difficult to propose a cosmological model which can give some understanding of the causal structure of the bouncing universe. In this article we address the issue related to the particle horizon problem in the bouncing universe models. It is shown that in many models the particle horizon does not exist, and consequently the horizon problem is trivially solved. In some cases a bouncing universe can have a particle horizon and we specify the conditions for its existence. In the absence of a particle horizon the Hubble surface specifies the causal structure of a bouncing universe. We specify the complex relationship between the Hubble surface and the particle horizon when the particle horizon exists. The article also address the issue related to the event horizon in a bouncing universe. A toy example of a bouncing universe is first presented where we specify the conditions which dictate the presence of a particle horizon. Next we specify the causal structures of three widely used bouncing models. The first case is related to quintom matter bounce model, the second one is loop quantum cosmology based bounce model and lastly $f(R)$ gravity induced bounce model. We present a brief discussion on the horizon problem in bouncing cosmologies. We point out that the causal structure of the various bounce models fit our general theoretical predictions.