The effects implied for the structure of compact objects by the modification of General Relativity produced by the generalization of the Lagrangian density to the form f (R) = R + αR 2 , where R is the Ricci curvature scalar, have been recently explored. It seems likely that this squared-gravity may allow heavier Neutron Stars (NSs) than GR. In addition, these objects can be useful to constrain free parameters of modifiedgravity theories. The differences between alternative gravity theories is enhanced in the strong gravitational regime. In this regime, because of the complexity of the field equations, perturbative methods become a good choice to treat the problem. Following previous works in the field, we performed a numerical integration of the structure equations that describe NSs in f (R)-gravity, recovering their mass-radius relations, but focusing on particular features that arise from this approach in the profiles of the NS interior.We show that these profiles run in correlation with the second-order derivative of the analytic approximation to the Equation of State (EoS), which leads to regions where the enclosed mass decreases with the radius in a counter-intuitive way. We reproduce all computations with a simple polytropic EoS to separate zerothorder modified gravity effects.

We consider a cosmological model where dark matter and dark energy feature a coupling that only affects their momentum transfer in the corresponding Euler equations. We perform a fit to cosmological observables and confirm previous findings within these scenarios that favour the presence of a coupling at more than 3σ. This improvement is mainly driven by cluster counts from Planck Sunyaev-Zeldovich data that we include as a certain prior. We subsequently perform a forecast for future J-PAS data and find that clustering measurements will permit to clearly discern the presence of an interaction within a few percent level with the uncoupled case at more than 10σ when the complete survey, covering 8500 sq. deg., is considered. We found that the inclusion of weak lensing measurements will not help to further constrain the coupling parameter. For completeness, we compare to forecasts for DESI and Euclid, which provide similar discriminating power.

We consider a scenario where dark energy and baryons are dynamically coupled without any energy transfer. In this scenario, the background cosmology is unaffected and, at the perturbations level, the coupling only appears through the corresponding Euler equations of dark energy and baryons. We then explore some phenomenological consequences of this scenario and their signatures in several cosmological observables. In particular, we show its ability to suppress the growth of cosmic structures. We also constrain the parameters of the model with cosmological data and show that an interaction of dark energy with baryons on cosmological scales is mildly favoured.

We analyse configurations of neutron stars in the so-called R-squared gravity in the Palatini formalism. Using a realistic equation of state we show that the mass-radius configurations are lighter than their counterparts in General Relativity. We also obtain the internal profiles, which run in strong correlation with the derivatives of the equation of state, leading to regions where the mass parameter decreases with the radial coordinate in a counter-intuitive way. In order to analyse such correlation, we introduce a parametrisation of the equation of state given by multiple polytropes, which allows us to explicitly control its derivatives. We show that, even in a limiting case where hard phase transitions in matter are allowed, the internal profile of the mass parameter still presents strange features and the calculated M −R configurations also yield neutron stars lighter than those obtained in General Relativity.2 Teppa Pannia et al.

We consider a cosmological scenario where the dark sector is described by two perfect fluids that interact through a velocity-dependent coupling. This coupling gives rise to an interaction in the dark sector driven by the relative velocity of the components, thus making the background evolution oblivious to the interaction and only the perturbed Euler equations are affected at first order. We obtain the equations governing this system with the Schutz-Sorkin Lagrangian formulation for perfect fluids and derive the corresponding stability conditions to avoid ghosts and Laplacian instabilities. As a particular example, we study a model where dark energy behaves as a radiation fluid at high redshift while it effectively becomes a cosmological constant in the late Universe. Within this scenario, we show that the interaction of both dark components leads to a suppression of the dark matter clustering at late times. We also argue the possibility that this suppression of clustering together with the additional dark radiation at early times can simultaneously alleviate the σ8 and H0 tensions.