We investigate cosmological models described by a scalar field with an exponential potential, and apply the stochastic formalism, which allows us to study how quantum field fluctuations give rise to stochastic noise. This modifies the classical dynamics of the scalar field at large scales, above a coarse-graining scale. In particular we explore how quantum field fluctuations perturb the equation of state on large scales which can lead to a quantum instability of the classical collapse solution below the Planck scale in the case of a pressureless collapse.

We address the issue of the existence of inequivalent definitions of gravitational mass in $R^{2}$-gravity. We present several definitions of gravitational mass, and discuss the formal relations between them. We then consider the concrete case of a static and spherically symmetric neutron star, and solve numerically the equations of motion for several values of the free parameter of the model. We compare the features of the mass-radius relations obtained for each definition of gravitational mass, and we comment on their dependence on the free parameter. We then argue that $R^{2}$-gravity is a valuable proxy to discuss the existence of inequivalent definitions of gravitational mass in a generic modified gravity theory, and present some comments on the general case.

In this paper we establish \textit{formulae} for the inflationary slow-roll parameters $\epsilon$, $\eta$ and $\zeta$ as functions of the Ricci scalar $R$ for $f(R)$ theories of gravity. As examples, we present the analytic and numerical solutions of $\epsilon$, $\eta$ and $\zeta$ as functions of the number of e-folds $N$ in two important instances: for the Starobinsky model and for a $f(R)$ reconstruction of the $\alpha$-Attractors. The highlight of our proposal is to rewrite the slow-roll parameters in terms of $f(R)$, which allows to find directly $n_{\rm s}$, $r$, $\alpha_{\rm s}$ and $f^{\text{equil}}_{\text{NL}}$ as functions of $R$ itself. We obtain that both models indicate a small contribution to the non-Gaussianity parameters, which are in good agreement with current observational constraints.

We show an analogy at high curvature between a $f(R) = R + aR^{n - 1} + bR^2$ theory and the $\alpha$-Attractors. We calculate the expressions of the parameters $a$, $b$ and $n$ as functions of $\alpha$ and the predictions of the model $f(R) = R + aR^{n - 1} + bR^2$ on the scalar spectral index $n_{\rm s}$ and the tensor-to-scalar ratio $r$. We find that the power law correction $R^{n - 1}$ allows for a production of gravitational waves enhanced with respect to the one in the Starobinsky model, while maintaining a viable prediction on $n_{\rm s}$. We numerically reconstruct the full $\alpha$-Attractors class of models testing the goodness of our high-energy approximation $f(R) = R + aR^{n - 1} + bR^2$. Moreover, we also investigate the case of a single power law $f(R) = \gamma R^{2 - \delta}$ theory, with $\gamma$ and $\delta$ free parameters. We calculate analytically the predictions of this model on the scalar spectral index $n_{\rm s}$ and the tensor-to-scalar ratio $r$ and the values of $\delta$ which are allowed from the current observational results. We find that $-0.015 < \delta < 0.016$, confirming once again the excellent agreement between the Starobinsky model and observation.

We analyse $f(R)$ theories of gravity from a dynamical system perspective, showing how the $R^2$ correction in Starobinsky's model plays a crucial role from the viewpoint of the inflationary paradigm. Then, we propose a modification of Starobinsky's model by adding an exponential term in the $f(R)$ Lagrangian. We show how this modification could allow to test the robustness of the model by means of the predictions on the scalar spectral index $n_s$.

The theoretical determination of braking indices of pulsars is still an open problem. In this paper we report results of a study concerning such determination based on a modification of the canonical model, which admits that pulsars are rotating magnetic dipoles, and on data from the seven pulsars with known braking indices. In order to test the modified model we predict ranges for the braking indices of other pulsars.