Aims. To date, more than 600 multi-planetary systems have been discovered. Due to the limitations of the detection methods, our knowledge of the systems is usually far from complete. In particular, for planetary systems discovered with the radial velocity (RV) technique, the inclinations of the orbital planes, and thus the mutual inclinations and planetary masses, are unknown. Our work aims to constrain the spatial configuration of several RV-detected extrasolar systems that are not in a mean-motion resonance. Methods. Through an analytical study based on a first-order secular Hamiltonian expansion and numerical explorations performed with a chaos detector, we identified ranges of values for the orbital inclinations and the mutual inclinations, which ensure the long-term stability of the system. Our results were validated by comparison with n-body simulations, showing the accuracy of our analytical approach up to high mutual inclinations (approx. 70{\deg}-80{\deg}). Results. We find that, given the current estimations for the parameters of the selected systems, long-term regular evolution of the spatial configurations is observed, for all the systems, i) at low mutual inclinations (typically less than 35{\deg}) and ii) at higher mutual inclinations, preferentially if the system is in a Lidov-Kozai resonance. Indeed, a rapid destabilisation of highly mutually inclined orbits is commonly observed, due to the significant chaos that develops around the stability islands of the Lidov-Kozai resonance. The extent of the Lidov-Kozai resonant region is discussed for ten planetary systems (HD 11506, HD 12661, HD 134987, HD 142, HD 154857, HD 164922, HD 169830, HD 207832, HD 4732, and HD 74156).

Spatial averaging and time evolving are non-commutative operations in General Relativity, which questions the reliability of the FLRW model. The long standing issue of the importance of backreactions induced by cosmic inhomogeneities is addressed for a toy model assuming a peak in the primordial spectrum of density perturbations and a simple CDM cosmology. The backreactions of initial Hubble-size inhomogeneities are determined in a fully relativistic framework, from a series of simulations using the BSSN formalism of numerical relativity. In the FLRW picture, these backreactions can be effectively described by two so-called morphon scalar fields, one of them acting at late time like a tiny cosmological constant. Initial density contrasts ranging from $10^{-2}$ down to $10^{-4}$, on scales crossing the Hubble radius between $z\sim 45 $ and $z\sim 1000$ respectively, i.e. comoving gigaparsec scales, mimic a Dark Energy (DE) component that can reach $\Omega_{\mathrm{DE}} \approx 0.7$ when extrapolated until today. A similar effect is not excluded for lower density contrasts but our results are then strongly contaminated by numerical noise and thus hardly reliable. A potentially detectable signature of this scenario is a phantom-like equation of state $w< -1$, at redshifts $z\gtrsim 4$ for a density contrast of $10^{-2}$ initially, relaxing slowly to $w \approx -1$ today. This new class of scenarios would send the fine-tuning and coincidence issues of Dark energy back to the mechanism at the origin of the primordial power spectrum enhancement, possibly in the context of inflation.