We investigate effective equations governing the volume expansion of spatially averaged portions of inhomogeneous cosmologies in spacetimes filled with an arbitrary fluid. This work is a follow-up to previous studies focused on irrotational dust models (Paper I) and irrotational perfect fluids (Paper II) in flow-orthogonal foliations of spacetime. It complements them by considering arbitrary foliations, arbitrary lapse and shift, and by allowing for a tilted fluid flow with vorticity. As for the first studies, the propagation of the spatial averaging domain is chosen to follow the congruence of the fluid, which avoids unphysical dependencies in the averaged system that is obtained. We present two different averaging schemes and corresponding systems of averaged evolution equations providing generalizations of Papers I and II. The first one retains the averaging operator used in several other generalizations found in the literature. We extensively discuss relations to these formalisms and pinpoint limitations, in particular regarding rest mass conservation on the averaging domain. The alternative averaging scheme that we subsequently introduce follows the spirit of Papers I and II and focuses on the fluid flow and the associated 1+3 threading congruence, used jointly with the 3+1 foliation that builds the surfaces of averaging. This results in compact averaged equations with a minimal number of cosmological backreaction terms. We highlight that this system becomes especially transparent when applied to a natural class of foliations which have constant fluid proper time slices.

Context. Atmospheric superrotating flows at the equator are an almost ubiquitous result of simulations of hot Jupiters, and a theory explaining how this zonally coherent flow reaches an equilibrium has been developed in the literature. However, this understanding relies on the existence of either an initial superrotating or a sheared flow, coupled with a slow evolution such that a linear steady state can be reached. Aims. A consistent physical understanding of superrotation is needed for arbitrary drag and radiative timescales, and the relevance of considering linear steady states needs to be assessed. Methods. We obtain an analytical expression for the structure, frequency and decay rate of propagating waves in hot Jupiter atmospheres around a state at rest in the 2D shallow water beta plane limit. We solve this expression numerically and confirm the robustness of our results with a 3D linear wave algorithm. We then compare with 3D simulations of hot Jupiter atmospheres and study the non linear momentum fluxes. Results. We show that under strong day night heating the dynamics does not transit through a linear steady state when starting from an initial atmosphere in solid body rotation. We further show that non linear effects favour the initial spin up of superrotation and that the acceleration due to the vertical component of the eddy momentum flux is critical to the initial development of superrotation. Conclusions. Overall, we describe the initial phases of the acceleration of superrotation, including consideration of differing radiative and drag timescales, and conclude that eddy-momentum driven superrotating equatorial jets are robust, physical phenomena in simulations of hot Jupiter atmospheres.

We introduce a generalization of the 4-dimensional averaging window function of Gasperini, Marozzi and Veneziano (2010) that may prove useful for a number of applications. The covariant nature of spatial scalar averaging schemes to address the averaging problem in relativistic cosmology is an important property that is implied by construction, but usually remains implicit. We employ here the approach of Gasperini et al. for two reasons. First, the formalism and its generalization presented here are manifestly covariant. Second, the formalism is convenient for disentangling the dependencies on foliation, volume measure, and boundaries in the averaged expressions entering in scalar averaging schemes. These properties will prove handy for simplifying expressions, but also for investigating extremal foliations and for comparing averaged properties of different foliations directly. The proposed generalization of the window function allows for choosing the most appropriate averaging scheme for the physical problem at hand, and for distinguishing between the role of the foliation itself and the role of the volume measure in averaged dynamic equations. We also show that one particular window function obtained from this generalized class results in an averaging scheme corresponding to that of a recent investigation by Buchert, Mourier and Roy (2018) and, as a byproduct, we explicitly show that the general equations for backreaction derived therein are covariant.

The subject of cosmological backreaction in General Relativity is often approached by coordinate-dependent and metric-based analyses. We present in this letter an averaging formalism for the scalar parts of Einstein's equations that is coordinate-independent and only functionally depends on a metric. This formalism is applicable to general 3+1 foliations of spacetime for an arbitrary fluid with tilted flow. We clarify the dependence on spacetime foliation and argue that this dependence is weak in cosmological settings. We also introduce a new set of averaged equations that feature a global cosmological time despite the generality of the setting, and we put the statistical nature of effective cosmologies into perspective.

We extend the general relativistic Lagrangian perturbation theory, recently developed for the formation of cosmic structures in a dust continuum, to the case of model universes containing a single fluid with a single-valued analytic equation of state. Using a coframe-based perturbation approach, we investigate evolution equations for structure formation in pressure-supported irrotational fluids that generate their rest-frame spacetime foliation. We provide master equations to first order for the evolution of the trace and traceless parts of barotropic perturbations that evolve in the perturbed space, where the latter describes the propagation of gravitational waves in the fluid. We illustrate the trace evolution for a linear equation of state and for a model equation of state describing isotropic velocity dispersion, and we discuss differences to the dust matter model, to the Newtonian case, and to standard perturbation approaches.

Self-consistent treatment of cosmological structure formation and expansion within the context of classical general relativity may lead to "extra" expansion above that expected in a structureless universe. We argue that in comparison to an early-epoch, extrapolated Einstein-de Sitter model, about 10-15% "extra" expansion is sufficient at the present to render superfluous the "dark energy" 68% contribution to the energy density budget, and that this is observationally realistic.

In relativistic inhomogeneous cosmology, structure formation couples to average cosmological expansion. A conservative approach to modelling this assumes an Einstein--de Sitter model (EdS) at early times and extrapolates this forward in cosmological time as a "background model" against which average properties of today's Universe can be measured. This requires adopting an early-epoch--normalised background Hubble constant $H_1^{bg}$. Here, we show that the $\Lambda$CDM model can be used as an observational proxy to estimate $H_1^{bg}$ rather than choose it arbitrarily. We assume (i) an EdS model at early times; (ii) a zero dark energy parameter; (iii) bi-domain scalar averaging---division of the spatial sections into over- and underdense regions; and (iv) virialisation (stable clustering) of collapsed regions. We find $H_1^{bg}= 37.7 \pm 0.4$ km/s/Mpc (random error only) based on a Planck $\Lambda$CDM observational proxy. Moreover, since the scalar-averaged expansion rate is expected to exceed the (extrapolated) background expansion rate, the expected age of the Universe should be much less than $2/(3 H_1^{bg}) = 17.3$ Gyr. The maximum stellar age of Galactic Bulge microlensed low-mass stars (most likely: 14.7 Gyr; 68\% confidence: 14.0--15.0 Gyr) suggests an age about a Gyr older than the (no-backreaction) $\Lambda$CDM estimate.

The admitted, conventional scenario to explain the complex spectral evolution of brown dwarfs (BD) since their first detections twenty years ago, has always been the key role played by micron-size condensates, called "dust" or "clouds", in their atmosphere. This scenario, however, faces major problems, in particular the J-band brightening and the resurgence of FeH absorption at the L to T transition, and a physical first-principle understanding of this transition is lacking. In this paper, we propose a new, completely different explanation for BD and extrasolar giant planet (EGP) spectral evolution, without the need to invoke clouds. We show that, due to the slowness of the CO/CH4 and N2/NH3 chemical reactions, brown dwarf (L and T, respectively) and EGP atmospheres are subject to a thermo-chemical instability similar in nature to the fingering or chemical convective instability present in Earth oceans and at the Earth core/mantle boundary. The induced small-scale turbulent energy transport reduces the temperature gradient in the atmosphere, explaining the observed increase in near infrared J-H and J-K colors of L dwarfs and hot EGPs, while a warming up of the deep atmosphere along the L to T transition, as the CO/CH4 instability vanishes, naturally solves the two aforementioned puzzles, and provides a physical explanation of the L to T transition. This new picture leads to a drastic revision of our understanding of BD and EGP atmospheres and their evolution.

We present an analysis of Spitzer/IRAC primary transit and secondary eclipse lightcurves measured for HD209458b, using Gaussian process models to marginalise over the intrapixel sensitivity variations in the 3.6 micron and 4.5 micron channels and the ramp effect in the 5.8 micron and 8.0 micron channels. The main advantage of this approach is that we can account for a broad range of degeneracies between the planet signal and systematics without actually having to specify a deterministic functional form for the latter. Our results do not confirm a previous claim of water absorption in transmission. Instead, our results are more consistent with a featureless transmission spectrum, possibly due to a cloud deck obscuring molecular absorption bands. For the emission data, our values are not consistent with the thermal inversion in the dayside atmosphere that was originally inferred from these data. Instead, we agree with another re-analysis of these same data, which concluded a non-inverted atmosphere provides a better fit. We find that a solar-abundance clear-atmosphere model without a thermal inversion underpredicts the measured emission in the 4.5 micron channel, which may suggest the atmosphere is depleted in carbon monoxide. An acceptable fit to the emission data can be achieved by assuming that the planet radiates as an isothermal blackbody with a temperature of $1484\pm 18$ K.

This work aims to improve the current understanding of the atmospheres of brown dwarfs, especially cold ones with spectral type T and Y, whose modeling is a current challenge. Silicate and iron clouds are believed to disappear at the photosphere at the L/T transition, but cloudless models fail to reproduce correctly the spectra of T dwarfs, advocating for the addition of more physics, e.g. other types of clouds or internal energy transport mechanisms. We use a one-dimensional (1D) radiative/convective equilibrium code ATMO to investigate this issue. This code includes both equilibrium and out-of-equilibrium chemistry and solves consistently the PT structure. Included opacity sources are H2-H2, H2-He, H2O, CO, CO2, CH4, NH3, K, Na, and TiO, VO if they are present in the atmosphere. We show that the spectra of Y dwarfs can be accurately reproduced with a cloudless model if vertical mixing and NH3 quenching are taken into account. T dwarf spectra still have some reddening in e.g. J - H compared to cloudless models. This reddening can be reproduced by slightly reducing the temperature gradient in the atmosphere. We propose that this reduction of the stabilizing temperature gradient in these layers, leading to cooler structures, is due to the onset of fingering convection, triggered by the destabilizing impact of condensation of very thin dust.