We compute the inertial modes of a freely rotating two-layer planetary model with an ellipsoidal inviscid fluid core and a perfectly rigid mantle. We present a method to derive analytical formulae for the frequencies of the Free Core Nutation (FCN) and Chandler Wobble (CW) which are valid to all orders of the dynamical flattening of the core and mantle, and we show how the FCN and CW are the direct generalisation of the purely fluid Spin-Over mode (SO) and of the Eulerian Wobble (EW) to the case where the mantle can oscillate freely around a state of steady rotation. Through a numerical computation for an axisymmetric (oblate spheroidal) planet, we demonstrate that all other inertial modes of the steadily rotating fluid core are also free modes of the freely rotating two-layer planet.

The Earth is a rapidly rotating body. The centrifugal pull makes its shape resemble a flattened ellipsoid and Coriolis forces support waves in its fluid core, known as inertial waves. These waves can lead to global oscillations, or modes, of the fluid. Periodic variations of the Earth's rotation axis (nutations) can lead to an exchange of angular momentum between the mantle and the fluid core and excite these inertial modes. In addition to viscous torques that exist regardless of the shape of the boundaries, the small flattening of the core-mantle boundary (CMB) allows inertial modes to exert pressure torques on the mantle. These torques effectively couple the rigid-body dynamics of the Earth with the fluid dynamics of the fluid core. Here we present the first high resolution numerical model that solves simultaneously the rigid body dynamics of the mantle and the Navier-Stokes equation for the liquid core. This method takes naturally into account dissipative processes in the fluid that are ignored in current nutation models. We find that the Free Core Nutation (FCN) mode, mostly a toroidal fluid flow if the mantle has a large moment of inertia, enters into resonance with nearby modes if the mantle's moment of inertia is reduced. These mode interactions seem to be completely analogous to the ones discovered by Schmitt (2006) in a uniformly rotating ellipsoid with varying flattening.

We propose a numerical method to compute the inertial modes of a container with near-spherical geometry based on the fully spectral discretisation of the angular and radial directions using spherical harmonics and Gegenbauer polynomial expansion respectively. This allows to solve simultaneously the Poincare equation and the no penetration condition as an algebraic polynomial eigenvalue problem. The inertial modes of an exact oblate spheroid are recovered to machine precision using an appropriate set of spheroidal coordinates. We show how other boundaries that deviate slightly from a sphere can be accommodated for with the technique of equivalent spherical boundary and we demonstrate the convergence properties of this approach for the triaxial ellipsoid.

Enceladus is characterised by a south polar hot spot associated with a large outflow of heat, the source of which remains unclear. We compute the viscous dissipation resulting from tidal and libration forcing in the moon's subsurface ocean using the linearised Navier-Stokes equation in a 3-dimensional spherical model. We conclude that libration is the dominant cause of dissipation at the linear order, providing up to about 0.001 GW of heat to the ocean, which remains insufficient to explain the (about) 10 GW observed by Cassini. We also illustrate how resonances with inertial modes can significantly augment the dissipation. Our work is an extension to Rovira-Navarro et al. [2019] to include the effects of libration. The model developed here is readily applicable to the study of other moons and planets.

We present a study of the fully relativistic spherical collapse in presence of quintessence using on Numerical Relativity, following the method proposed by the authors in a previous article [arXiv:1409.3476]. We ascertain the validity of the method by studying the evolution of a spherically symmetric quintessence inhomogeneity on a de Sitter background and we find that it has an impact on the local expansion around the centre of coordinates. We then proceed to compare the results of our method to those of the more largely adopted top-hat model. We find that quintessence inhomogeneities do build up under the effect that matter inhomogeneities have on the local space-time yet remain very small due to the presence of momentum transfer from the over-dense to the background regions. We expect that these might have an even more important role in modified theories of gravitation.

We present a fully relativistic numerical method for the study of cosmological problems using the Baumgarte-Shapiro-Shibata-Nakamura formalism on a dynamical Friedmann-Lema\^itre-Robertson-Walker background. This has many potential applications including the study of the growth of structures beyond the linear regime. We present one such application by reproducing the Lema\^itre-Tolman-Bondi solution for the collapse of pressureless matter with arbitrary lapse function. The regular and smooth numerical solution at the center of coordinates proceeds in a natural way by relying on the Partially Implicit Runge-Kutta algorithm described in Montero and Cordero-Carri\'on [arXiv:1211.5930]. We generalize the usual radiative outer boundary condition to the case of a dynamical background and show the stability and convergence properties of the method in the study of pure gauge dynamics on a de Sitter background.

We revisit the status of hybrid inflation in the light of Planck and recent BICEP2 results, taking care of possible transient violations of the slow-roll conditions as the field passes from the large field to the vacuum dominated phase. The usual regime where observable scales exit the Hubble radius in the vacuum dominated phase predicts a blue scalar spectrum, which is ruled out. But whereas assuming slow-roll one expects this regime to be generic, by solving the exact dynamics we identify the parameter space for which the small field phase is naturally avoided due to slow-roll violations at the end of the large field phase. When the number of e-folds generated at small field is negligible, the model predictions are degenerated with those of a quadratic potential. There exists also a transitory case for which the small field phase is sufficiently long to affect importantly the observable predictions. Interestingly, in this case the spectral index and the tensor to scalar ratio agree respectively with the best fit of Planck and BICEP2. This results in a \Delta \chi^2 \simeq 5.0 in favor of hybrid inflation for Planck+BICEP2 (\Delta \chi^2 \simeq 0.9 for Planck only). The last considered regime is when the critical point at which inflation ends is located in the large field phase. It is constrained to be lower than about ten times the reduced Planck mass. The analysis has been conducted with the use of Markov-Chain-Monte-Carlo bayesian method, in a reheating consistent way, and we present the posterior probability distributions for all the model parameters.

The dark energy scalar field is here presented as a mean-field effect arising from the collective motion of interacting structures on an expanding lattice. This cosmological analogue to solid-state soft phonons in an unstable crystal network is shown to produce cosmic acceleration while mimicking phantom equation of state. From an analysis of the Hubble diagram of type Ia supernovae, we present constraints and discussion on the parameters of the cosmic Lagrange chain, as well as on time-variation of the soft phonon equation of state.