Using a joint statistical analysis, we test a four-dimensional FLWR model embedded in a five-dimensional bulk based on the Nash-Greene embedding theorem. Performing a Markov Chain Monte Carlo (MCMC) modelling, we combine observational data sets as those of the recent ``Gold 2018'' growth data, the best-fit Planck2018/$\Lambda$CDM parameters on the Cosmic Microwave Background (CMB), the Baryon Acoustic Oscillations (BAO) measurements, the Pantheon Supernovae type Ia and the Hubble parameter data with redshift ranging from $0.01 < z < 2.3$ to impose restrictions on the model. From linear Nash-Green fluctuations of the metric, we obtain the related perturbed equations in longitudinal Newtonian gauge to obtain the evolution of growth matter. We show that a mild alleviation is obtained from the $\sigma$ tension between the growth amplitude factor and the matter content $(\sigma_8$-$\Omega_m)$ of the observations from CMB and Large Scale Structure (LSS) probes with degeneracies on the parameters. We also apply the Om(z) diagnosis to distinguish the model from $\Lambda$CDM. The Akaike Information Criterion (AIC) is used and we find a relative statistical consistence of the present model with both $\Lambda$CDM and $w$CDM models up to 1$\%$ of percentage difference at early times.

The acceleration of the universe is described as a consequence of the extrinsic curvature of the space-time embedded in a bulk space, defined by the Einstein-Hilbert. Using the linear approximation of Nash-Green theorem, we obtain the related perturbed equations in which just the gravitational-tensor field equations contribute to propagation of cosmological perturbations. In accordance with Big Bang Nucleosynthesis and solar constraints, we calculate numerically the effective Newtonian function Geff to constrain the related parameters of the model. We numerically solve the growth density equation for two possible family of solutions leading to an interesting overdensity and, in some cases, a mild damping of the growth profiles, with a top amplification of the growth perturbations around 14% in comparison with LCDM model and quintessence. The effective gravitational Phi and Newtonian curvature Psi are also analysed showing mild perturbations in early times induced only by the extrinsic curvature differently from the LambdaCDM standards.

The Hubble parameter H(z), as a function of redshift, is modified by the presence of a new term originated from the extrinsic curvature in an embedded space-time. Assuming an asymptotic expansion factor $a\sim 0$ or, equivalently $z\rightarrow \infty$, it is obtained a nearly resemblance of the present model ($\beta$-model) with quintessence $w$CDM model with $H(z)\sim H_{wCDM}(z)$ at background level. In conjunction with $\Lambda$CDM, we test the models using a pack of recent datasets like that of the ``Gold 2018'' growth data, the best-fit Planck2018/$\Lambda$CDM parameters on the Cosmic Microwave Background (CMB), the Baryon Acoustic Oscillations (BAO) measurements, the Pantheon Supernovae type Ia and the Hubble parameter data with redshift ranging from $0.01 < z < 2.3$. Performing the Akaike Information Criterion (AIC) to ascertain the statistical viability of the model from Jeffreys' scale, we apply a joint likelihood analysis to the data with the Markov Chain Monte Carlo (MCMC) method. We find that the present model is in very good agreement with observations with a close statistical equivalence with the $\Lambda$CDM and $w$CDM cosmologies at 1-$\sigma$ level. We also show that a mild alleviation of the $\sigma$ tension between the growth amplitude factor and the matter content $(\sigma_8$-$\Omega_m)$ of the observations from CMB and Large Scale Structure (LSS) probes. A comparison of the aforementioned full pack of data in MCMC chains is made with the resulting MCMC from Pantheon SNIa+$H(z)$ in order to analyse the sensitivity of the models and how they respond to a cosmography analysis on the evolution of $H(z)$ and the deceleration parameter $q(z)$. In this sense, we find that only the $\beta$-model can be unaffected to the variations of the previous datasets due to its several minima in the likelihood (degeneracies) of the related parameters with an overall percentage relative difference only up to 4$\%$.

By using geodesic equations to obtain a gravitational potential generated from a point-like source, we end up in the concept of a nearly Newtonian gravity to analyse effective potentials of quasi-circular orbits. By means of an approximate solution from an axially static and symmetric Weyl metric, we study an effective gravitational potential to obtain its related rotation curves, orbital planes and orbits. Moreover, using as initial condition a Plummer sphere, some prospects on star cluster disruption are also discussed in this framework.