A massive, nonrelativistic scalar field in an expanding spacetime is usually approximated by a pressureless perfect fluid, which leads to the standard conclusion that such a field can play the role of cold dark matter. In this paper, we systematically study these approximations, incorporating subleading corrections. We provide two equivalent effective descriptions of the system, each of which offers its own advantages and insights: (i) A nonrelativistic effective field theory (EFT) with which we show that the relativistic corrections induce an effective self-interaction for the nonrelativistic field. As a byproduct, our EFT also allows one to construct the exact solution, including oscillatory behavior, which is often difficult to achieve from the exact equations. (ii) An effective (imperfect) fluid description, with which we demonstrate that, for a perturbed Friedmann-Lema\^{i}tre-Robertson-Walker (FLRW) universe: (a) The pressure is small but nonzero (and positive), even for a free theory with no tree-level self-interactions. (b) The sound speed of small fluctuations is also nonzero (and positive), reproducing already known leading-order results, correcting a subdominant term, and identifying a new contribution that had been omitted in previous analyses. (c) The fluctuations experience a negative effective bulk viscosity. The positive sound speed and the negative bulk viscosity act in favor of and against the growth of overdensities, respectively. The net effect may be considered a smoking gun for ultra-light dark matter.

In this work we study vacuum decay and bubble nucleation in models of $f(R)$ higher curvature gravity. Building upon the analysis of Coleman-De Luccia (CDL), we present the formalism to calculate the Euclidean action and the bounce solution for a general $f(R)$ gravity in the thin wall approximation. We calculate the size of the nucleated bubble and the decay exponent for the Starobinsky model and its higher power extensions. We have shown that in the Starobinsky model with a typical potential the nucleated bubble has a larger size in comparison to the CDL bubble and with a lower tunneling rate. However, for higher power extension of the Starobinsky model the size of the bubble and the tunneling exponent can be larger or smaller than the CDL bubble depending on the model parameters. As a counterintuitive example, we have shown that a bubble with a larger size than the CDL bubble but with a higher nucleation rate can be formed in $f(R)$ gravity.