Minkowski functionals and Shapefinders shed light on the connectedness of large-scale structure by determining its topology and morphology. We use a sophisticated code, SURFGEN2, to measure the Minkowski functionals and Shapefinders of individual clusters by modelling cluster surfaces using the 'Marching Cube 33' triangulation algorithm. In this paper, we study the morphology of simulated neutral hydrogen (HI) density fields using Shapefinders at various stages of reionization from the excursion set approach. Accompanying the Shapefinders, we also employ the 'largest cluster statistic' (LCS) to understand the percolation process. Percolation curves demonstrate that the non-Gaussianity in the HI field increases as reionization progresses. The large clusters in both the HI overdense and underdense excursion sets possess similar values of "thickness" ($T$), as well as "breadth" ($B$), but their third Shapefinder - "length" ($L$) - becomes almost proportional to their volume. The large clusters in both HI overdense and underdense segments are overwhelmingly filamentary. The 'cross-section' of a filamentary cluster can be estimated using the product of the first two Shapefinders, $T \times B$. Hence the cross sections of the large clusters at the onset of percolation do not vary much with volume and their sizes only differ in terms of their lengths. This feature appears more vividly in HI overdense regions than in underdense regions and is more pronounced at lower redshifts which correspond to an advanced stage of reionization.

Using Shapefinders, which are ratios of Minkowski functionals, we study the morphology of neutral hydrogen (HI) density fields, simulated using semi-numerical technique (inside-out), at various stages of reionization. Accompanying the Shapefinders, we also employ the 'largest cluster statistic' (LCS), originally proposed in Klypin and Shandarin (1993), to study the percolation in both neutral and ionized hydrogen. We find that the largest ionized region is percolating below the neutral fraction $x_{HI} \lesssim 0.728$ (or equivalently $z \lesssim 9$). The study of Shapefinders reveals that the largest ionized region starts to become highly filamentary with non-trivial topology near the percolation transition. During the percolation transition, the first two Shapefinders - 'thickness' ($T$) and 'breadth' ($B$) - of the largest ionized region do not vary much, while the third Shapefinder - 'length' ($L$) - abruptly increases. Consequently, the largest ionized region tends to be highly filamentary and topologically quite complex. The product of the first two Shapefinders, $T\times B$, provides a measure of the 'cross-section' of a filament-like ionized region. We find that, near percolation, the value of $T\times B$ for the largest ionized region remains stable at $\sim 7$ Mpc$^2$ (in comoving scale) while its length increases with time. Interestingly all large ionized regions have similar cross-sections. However their length shows a power-law dependence on their volume, $L\propto V^{0.72}$, at the onset of percolation.

We derive an analytical expression for the growth rate of matter density perturbations on the phantom brane (which is the normal branch of the Dvali-Gabadadze-Porrati model). This model is characterized by a phantomlike effective equation of state for dark energy at the present epoch. It agrees very well with observations. We demonstrate that the traditional parametrization $f=\Omega_m^\gamma$ with a quasiconstant growth index $\gamma$ is not successful in this case. Based on a power series expansion at large redshifts, we propose a different parametrization for this model: $f=\Omega_m^\gamma\left(1+\frac{b}{\ell H}\right)^\beta$, where $\beta$ and $b$ are constants. Our numerical simulations demonstrate that this new parametrization describes the growth rate with great accuracy - the maximum error being $\leq 0.1\%$ for parameter values consistent with observations.

We describe a new class of dark energy (DE) models which behave like cosmological trackers at early times. These models are based on the $\alpha$-attractor set of potentials, originally discussed in the context of inflation. The new models allow the current acceleration of the universe to be reached from a wide class of initial conditions. Prominent examples of this class of models are the potentials $\coth\varphi$ and $\cosh\varphi$. A remarkable feature of this new class of models is that they lead to large enough negative values of the equation of state at the present epoch, consistent with the observations of accelerated expansion of the universe, from a very large initial basin of attraction. They therefore avoid the fine tuning problem which afflicts many models of DE.

In [Schmidt PRD 80 123003 (2009)], the author suggested that dynamical dark energy (DDE) propagating on the phantom brane could mimick $\Lambda$CDM. Schmidt went on to derive a phenomenological expression for $\rho_{\rm DE}$ which could achieve this. We demonstrate that while Schmidt's central premise is correct, the expression for $\rho_{\rm DE}$ derived in Schmidt (2009) is flawed. We derive the correct expression for $\rho_{\rm DE}$ which leads to $\Lambda$CDM-like expansion on the phantom brane. We also show that DDE on the brane can be associated with a Quintessence field and derive a closed form expression for its potential $V(\phi)$. Interestingly the $\alpha$-attractor based potential $V(\phi) \propto \coth^2{\lambda\phi}$ makes braneworld expansion resemble $\Lambda$CDM. However the two models can easily be distinguished on the basis of density perturbations which grow at different rates on the braneworld and in $\Lambda$CDM.

The phantom brane has several important distinctive features: (i) Its equation of state is phantom-like, but there is no future `big rip' singularity, (ii) the effective cosmological constant on the brane is dynamically screened, because of which the expansion rate is {\em smaller} than that in $\Lambda$CDM at high redshifts. In this paper, we constrain the Phantom braneworld using distance measures such as Type Ia supernovae (SNeIa), Baryon Acoustic Oscillations (BAO), and the compressed Cosmic Microwave Background (CMB) data. We find that the simplest braneworld models provide a good fit to the data. For instance, BAO +SNeIa data can be accommodated by the braneworld for a large region in parameter space $0 < \Omega_l < 0.3$ at $1\sigma$. The Hubble parameter can be as high as $H_0 < 78$ km/s/Mpc, and the effective equation of state at present can show phantom-like behaviour with $w_0 < -1.2$ at $1\sigma$. We note a correlation between $H_0$ and $w_0$, with higher values of $H_0$ leading to a lower, and more phantom-like, value of $w_0$. Inclusion of CMB data provides tighter constraints $\Omega_l < 0.1$. (Here $\Omega_l$ encodes the ratio of the five and four dimensional Planck mass.) The Hubble parameter in this case is more tightly constrained to $H_0 < 71$ km/s/Mpc, and the effective equation of state to $w_0 < -1.1$. Interestingly, we find that the universe is allowed be closed or open, with $-0.5 < \Omega_{\kappa} < 0.5$, even on including the compressed CMB data. There appears to be some tension in the low and high $z$ BAO data which may either be resolved by future data, or act as a pointer to interesting new cosmology.

We obtain a closed system of equations for scalar perturbations in a multi-component braneworld. Our braneworld possesses a phantom-like equation of state at late times, $w_{\rm eff} < -1$, but no big-rip future singularity. In addition to matter and radiation, the braneworld possesses a new effective degree of freedom - the 'Weyl fluid' or 'dark radiation'. Setting initial conditions on super-Hubble spatial scales at the epoch of radiation domination, we evolve perturbations of radiation, pressureless matter and the Weyl fluid until the present epoch. We observe a gradual decrease in the amplitude of the Weyl-fluid perturbations after Hubble-radius crossing, which results in a negligible effect of the Weyl fluid on the evolution of matter perturbations on spatial scales relevant for structure formation. Consequently, the quasi-static approximation of Koyama and Maartens provides a good fit to the exact results during the matter-dominated epoch. We find that the late-time growth of density perturbations on the brane proceeds at a faster rate than in $\Lambda$CDM. Additionally, the gravitational potentials $\Phi$ and $\Psi$ evolve differently on the brane than in $\Lambda$CDM, for which $\Phi = \Psi$. On the brane, by contrast, the ratio $\Phi/\Psi$ exceeds unity during the late matter-dominated epoch ($z \lesssim 50$). These features emerge as "smoking gun" tests of phantom brane cosmology and allow predictions of this scenario to be tested against observations of galaxy clustering and large-scale structure.

We explore the possibility of emergent cosmology using the effective potential formalism. We discover new models of emergent cosmology which satisfy the constraints posed by the cosmic microwave background (CMB). We demonstrate that, within the framework of modified gravity, the emergent scenario can arise in a universe which is spatially open/closed. By contrast, in general relativity (GR) emergent cosmology arises from a spatially closed past-eternal Einstein Static Universe (ESU). In GR the ESU is unstable, which creates fine tuning problems for emergent cosmology. However, modified gravity models including Braneworld models, Loop Quantum Cosmology (LQC) and Asymptotically Free Gravity result in a stable ESU. Consequently, in these models emergent cosmology arises from a larger class of initial conditions including those in which the universe eternally oscillates about the ESU fixed point. We demonstrate that such an oscillating universe is necessarily accompanied by graviton production. For a large region in parameter space graviton production is enhanced through a parametric resonance, casting serious doubts as to whether this emergent scenario can be past-eternal.