Halo model is a physically intuitive method for modelling the non-linear power spectrum, especially for the alternatives to the standard $\Lambda$CDM models. In this paper, we exam the Sheth-Tormen barrier formula adopted in the previous \texttt{CHAM} method \citep{2018MNRAS.476L..65H}. As an example, we model the ellipsoidal collapse of top-hat dark matter haloes in $f(R)$ gravity. A good agreement between Sheth-Tormen formula and our result is achieved. The relative difference in the ellipsoidal collapse barrier is less than or equal to $1.6\%$. Furthermore, we verify that, for F4 and F5 cases of Hu-Sawicki $f(R)$ gravity, the screening mechanism do not play a crucial role in the non-linear power spectrum modelling up to $k\sim1[h/{\rm Mpc}]$. We compare two versions of modified gravity modelling, namely with/without screening. We find that by treating the effective Newton constant as constant number ($G_{\rm eff}=4/3G_N$) is acceptable. The scale dependence of the gravitational coupling is sub-relevant. The resulting spectra in F4 and F5, are in $0.1\%$ agreement with the previous \texttt{CHAM} results. The published code is accelerated significantly. Finally, we compare our halo model prediction with N-body simulation. We find that the general spectrum profile agree, qualitatively. However, via the halo model approach, there exists a systematic under-estimation of the matter power spectrum in the co-moving wavenumber range between $0.3 h/{\rm Mpc}$ and $3 h/{\rm Mpc}$. These scales are overlapping with the transition scales from two halo term dominated regimes to those of one halo term dominated.

We present a model-independent measurement of spatial curvature $\Omega_{k}$ in the Friedmann-Lema\^itre-Robertson-Walker (FLRW) universe, based on observations of the Hubble parameter $H(z)$ using cosmic chronometers, and a Gaussian Process (GP) reconstruction of the HII galaxy Hubble diagram. We show that the imposition of spatial flatness (i.e., $\Omega_k=0$) easily distinguishes between the Hubble constant measured with {\it Planck} and that based on the local distance ladder. We find an optimized curvature parameter $\Omega_{k} = -0.120^{+0.168}_{-0.147}$ when using the former (i.e., $H_0=67.66\pm0.42 \, \mathrm{km}\,\mathrm{s}^{-1} \,\mathrm{Mpc}^{-1}$), and $\Omega_{k} = -0.298^{+0.122}_{-0.088}$ for the latter ($H_0=73.24\pm 1.74 \,\mathrm{km}\,\mathrm{s}^{-1} \,\mathrm{Mpc}^{-1}$). The quoted uncertainties are extracted by Monte Carlo sampling, taking into consideration the covariances between the function and its derivative reconstructed by GP. These data therefore reveal that the condition of spatial flatness favours the {\it Planck} measurement, while ruling out the locally inferred Hubble constant as a true measure of the large-scale cosmic expansion rate at a confidence level of $\sim 3\sigma$.

A validation of the cosmic distance duality (CDD) relation, eta(z)=(1+z)^2 d_A(z)/d_L(z)=1, coupling the luminosity (d_L) and angular-diameter (d_A) distances, is crucial because its violation would require exotic new physics. We present a model-independent test of the CDD, based on strong lensing and a reconstruction of the HII galaxy Hubble diagram using Gaussian Processes, to confirm the validity of the CDD at a very high level of confidence. Using parameterizations eta(z) = 1 + eta_0 z and eta(z) = 1 + eta_1 z + eta_2 z^2, our best-fit results are eta_0 = 0.0147^{+0.056}_{-0.066}, and eta_1 = 0.1091^{+0.1680}_{-0.1568} and eta_2 = -0.0603^{+0.0999}_{-0.0988}, respectively. In spite of these strong constraints, however, we also point out that the analysis of strong lensing using a simplified single isothermal sphere (SIS) model for the lens produces some irreducible scatter in the inferred CDD data. The use of an extended SIS approximation, with a power-law density structure, yields very similar results, but does not lessen the scatter due to its larger number of free parameters, which weakens the best-fit constraints. Future work with these strong lenses should therefore be based on more detailed ray-tracing calculations to determine the mass distribution more precisely.